MeteoInfoLab脚本示例：AOD

cat1106123

谢谢老师，MOD08_D3数据文件太大了，有90多MB（所以我怀疑是不是下载错了），不让上传...MOD04_L2的可以上传，麻烦您了

cat1106123

下载金钱-5不知道怎么弄掉，不好意思哈王老师

MeteoInfo

cat1106123 发表于 2016-2-8 18:06
下载金钱-5不知道怎么弄掉，不好意思哈王老师

没事，论坛的money我是不缺的。大文件你可以通过百度云之类的网络空间共享。

cat1106123

MeteoInfo 发表于 2016-2-8 18:17
没事，论坛的money我是不缺的。大文件你可以通过百度云之类的网络空间共享。

王老师，MOD08_D3发到百度云了http://pan.baidu.com/s/1o7b721G

cat1106123

MeteoInfo 发表于 2016-2-8 18:17
没事，论坛的money我是不缺的。大文件你可以通过百度云之类的网络空间共享。

王老师，MOD08_D3发到百度云了http://pan.baidu.com/s/1o7b721G

MeteoInfo

cat1106123 发表于 2016-2-8 19:28
王老师，MOD08_D3发到百度云了http://pan.baidu.com/s/1o7b721G

需要下载最新版本的MeteoInfo。
MOD08：

1. f = addfile('D:/Temp/hdf/MOD08_D3.A2015307.051.2015308233133.hdf')
   
2. v_aod = f['Optical_Depth_Land_And_Ocean_Mean']
   
3. v_aod.setdimrev(0, True)
   
4. aod = v_aod[[15,55],[70,140]]
   
5. #Plot
   
6. axesm()
   
7. china = shaperead('D:/Temp/map/bou2_4p.shp')
   
8. geoshow(china, linecolor='gray')
   
9. world = shaperead('D:/Temp/map/country1.shp')
   
10. geoshow(world)
    
11. levels = arange(0, 1.0, 0.05)
    
12. layer = contourfm(aod, levels, interpolate=True)
    
13. #layer = imshowm(aod, levels)
    
14. title('Aerosol Optical Depth at 550 nm')
    
15. colorbar(layer)

f = addfile('D:/Temp/hdf/MOD08_D3.A2015307.051.2015308233133.hdf')<br /> v_aod = f['Optical_Depth_Land_And_Ocean_Mean']<br /> v_aod.setdimrev(0, True)<br /> aod = v_aod[[15,55],[70,140]]<br /> #Plot<br /> axesm()<br /> china = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(china, linecolor='gray')<br /> world = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(world)<br /> levels = arange(0, 1.0, 0.05)<br /> layer = contourfm(aod, levels, interpolate=True)<br /> #layer = imshowm(aod, levels)<br /> title('Aerosol Optical Depth at 550 nm')<br /> colorbar(layer)

MOD04:

1. f = addfile('D:/Temp/hdf/MOD04_L2.A2015307.0225.006.2015309143224.hdf')
   
2. lon = f['Longitude'][:,:]
   
3. lat = f['Latitude'][:,:]
   
4. data = f['Optical_Depth_Land_And_Ocean'][:,:]
   
5. axesm()
   
6. mlayer = shaperead('D:/Temp/map/country1.shp')
   
7. geoshow(mlayer, edgecolor='k')
   
8. layer = surfacem(lon, lat, data)
   
9. colorbar(layer)

f = addfile('D:/Temp/hdf/MOD04_L2.A2015307.0225.006.2015309143224.hdf')<br /> lon = f['Longitude'][:,:]<br /> lat = f['Latitude'][:,:]<br /> data = f['Optical_Depth_Land_And_Ocean'][:,:]<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer, edgecolor='k')<br /> layer = surfacem(lon, lat, data)<br /> colorbar(layer)

cat1106123

老师，一再麻烦您真是不好意思哈...我下载的版本是下面图上红圈里的MeteoInfo，然后又用您刚刚发的代码试了一下，是图2这个样子......，属性错误。已经一头雾水了，盼您再次指点......非常感谢！

cat1106123

而且根据图上来看，不知是因为缺测还是别的原因，空白的地方好多啊。请教一下老师这是什么原因呢？另外请问老师有没有什么相关教材之类的，我好买来系统学习一下？再次感谢您~

MeteoInfo

cat1106123 发表于 2016-2-9 17:48
而且根据图上来看，不知是因为缺测还是别的原因，空白的地方好多啊。请教一下老师这是什么原因呢？另外请问 ...

我测试过，应该不会出现你说的问题。卫星数据会有很多缺测值，这很正常（有云等情况）。所有的教程都在本论坛和MeteoInfo的网站上了。

cat1106123

好的O(∩_∩)O好的~再次感谢老师！！打扰这么久，您辛苦啦！！不过还是想问，
Traceback (most recent call last):
  File "<iostream>", line 3, in <module>
AttributeError: DimVariable instance has no attribute 'setdimrev'
这是为什么呢? setdimrev 出了什么问题呢？是因为我下载的软件版本不对嘛？