MeteoInfoLab脚本示例：AOD

MeteoInfo

cat1106123 发表于 2016-2-9 21:16
好的O(∩_∩)O好的~再次感谢老师！！打扰这么久，您辛苦啦！！不过还是想问，
Traceback (most r ...

事实上MeteoInfo Java 1.3R7里刚给DimVariable类增加了setdimrev()函数，你可以尝试重新下载1.3R7试试。你还可以在MeteoInfo目录下的pylib\mipylib子目录中打开DimVariable.py文件（用任何一款文本编辑软件即可），看看里面是否有下面的语句:
    def setdimrev(self, idx, reverse):
        self.dims[idx].setReverse(reverse)

实际上是因为数据在Y轴方向（纬度方向）数据是反向的造成之前的错误，要求数据纬度方向是从南向北排列，实际数据是从北向南，因此必须将数据Y轴反向才可以。所以增加了setdimrev()这个函数来处理这个问题。

cat1106123

终于顺利出图啦！！！简直感动哭！！谢谢老师的耐心指导~前几天下载的版本里确实没有您说的setdimrev()函数，导致程序识别不出总报错，今天又重新下载了一次，dimvariable.py中有这个函数啦，大成功！真是个方便的好软件呢！以后也要多多学习！再次感谢老师的详细讲解，新年最棒的礼物！

wanghp

很好很强大

931404656

我画的怎么没有空白

chenfei0430

正在学习下meteoinfo 出来AOD

chenfei0430

王老师，调试了好几次都调不出来，请您帮我看看问题所在。

MeteoInfo

chenfei0430 发表于 2016-7-1 13:39
王老师，调试了好几次都调不出来，请您帮我看看问题所在。

1、此数据是Swath，不是格点数据，需要不同的处理方式。
2、文件中没有Optical_Depth_Land_And_Ocean_Mean变量名，有Optical_Depth_Land_And_Ocean，但是变量中所有值都是缺测值（没有有意义的数据）。
3、此类文件处理的参考脚本：

1. fn = 'E:/Temp/MYD04_L2.A2016121.1755.006.2016122152651.hdf'
   
2. f = addfile(fn)
   
3. lon = f['Longitude'][:,:]
   
4. lat = f['Latitude'][:,:]
   
5. data = f['Optical_Depth_Land_And_Ocean'][:,:]
   
6. axesm()
   
7. mlayer = shaperead('D:/Temp/map/country1.shp')
   
8. geoshow(mlayer, edgecolor='k')
   
9. layer = surfacem(lon, lat, data)
   
10. colorbar(layer)
    
11. xlim(-180, 180)
    
12. ylim(-90, 90)

fn = 'E:/Temp/MYD04_L2.A2016121.1755.006.2016122152651.hdf'<br /> f = addfile(fn)<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)<br /> xlim(-180, 180)<br /> ylim(-90, 90)

chenfei0430

MeteoInfo 发表于 2016-7-1 14:25
1、此数据是Swath，不是格点数据，需要不同的处理方式。
2、文件中没有Optical_Depth_Land_And_Ocean_Me ...

王老师，我这个数据是对的吗？能运行起来不

MeteoInfo

chenfei0430 发表于 2016-7-1 17:00
王老师，我这个数据是对的吗？能运行起来不

这个有数据。

1. fn = 'C:/Temp/MYD04_L2.A2016153.0420.006.2016158224104.hdf'
   
2. f = addfile(fn)
   
3. lon = f['Longitude'][:,:]
   
4. lat = f['Latitude'][:,:]
   
5. data = f['Optical_Depth_Land_And_Ocean'][:,:]
   
6. data[data<0] = nan
   
7. axesm()
   
8. mlayer = shaperead('D:/Temp/map/country1.shp')
   
9. geoshow(mlayer, edgecolor='k')
   
10. layer = surfacem(lon, lat, data)
    
11. colorbar(layer)
    
12. xlim(110, 150)
    
13. ylim(15, 40)

fn = 'C:/Temp/MYD04_L2.A2016153.0420.006.2016158224104.hdf'<br /> f = addfile(fn)<br /> lon = f['Longitude'][:,:]<br /> lat = f['Latitude'][:,:]<br /> data = f['Optical_Depth_Land_And_Ocean'][:,:]<br /> data[data<0] = nan<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer, edgecolor='k')<br /> layer = surfacem(lon, lat, data)<br /> colorbar(layer)<br /> xlim(110, 150)<br /> ylim(15, 40)

chenfei0430

MeteoInfo 发表于 2016-7-1 19:41
这个有数据。

多谢，王老师，已经算出来了，请教meteoinfo有没有批处理功能？