MeteoInfoLab脚本示例：TRMM 2A12 HDF数据

MeteoInfo

TRMM 2A12 HDF数据是卫星观测的SWATH数据（轨道数据），比格点数据处理起来要麻烦一些。数据的经纬度保存在geolocation变量中，需要先将经纬度数据读出来（均为2维数组），然后读取云水含量数据（cldWater）。虽然都是2维数组，但并不是格点数据，相当于2维的散点数据，点数相当多（3019*208）如果用散点图来绘制的话会非常的慢，需要将其插值为格点数据。插值方法建议选择nearest，该方法速度最快。插值后的格点数据用imshowm函数显示为图像就很快了。

脚本程序：

1. #Add data file
   
2. folder = 'D:/Temp/hdf/'
   
3. fns = '2A12_090711_66394_6.HDF'
   
4. fn = folder + fns
   
5. f = addfile(fn)
   
6. #Get data variable
   
7. geo = f['geolocation']
   
8. lat = geo[:,:,0]
   
9. lon = geo[:,:,1]
   
10. water = f['cldWater'][:,:,9]
    
11. water = water / 1000
    
12. #Interpolate to grid data
    
13. gwater, x_g, y_g = griddata((lon, lat), water, method='neareast')
    
14. #Plot
    
15. axesm()
    
16. mlayer = shaperead('D:/Temp/map/country1.shp')
    
17. geoshow(mlayer)
    
18. levs = arange(0, 10, 0.5)
    
19. layer = imshowm(x_g, y_g, gwater, levs)
    
20. colorbar(layer, orientation='horizontal')
    
21. title(fns + ' (cldWater level=9)')
    
22. axism([-180, 180, -90, 90])

#Add data file<br /> folder = 'D:/Temp/hdf/'<br /> fns = '2A12_090711_66394_6.HDF'<br /> fn = folder + fns<br /> f = addfile(fn)<br /> #Get data variable<br /> geo = f['geolocation']<br /> lat = geo[:,:,0]<br /> lon = geo[:,:,1]<br /> water = f['cldWater'][:,:,9]<br /> water = water / 1000<br /> #Interpolate to grid data<br /> gwater, x_g, y_g = griddata((lon, lat), water, method='neareast')<br /> #Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> levs = arange(0, 10, 0.5)<br /> layer = imshowm(x_g, y_g, gwater, levs)<br /> colorbar(layer, orientation='horizontal')<br /> title(fns + ' (cldWater level=9)')<br /> axism([-180, 180, -90, 90])

付亚男

王老师辛苦，看您写的东西总是获益良多。

kongfeng0824

太强大了               

kongfeng0824

王老师有出书的计划吗？感觉这么多教程，绝对可以出书了。

MeteoInfo

kongfeng0824 发表于 2015-7-14 14:57
王老师有出书的计划吗？感觉这么多教程，绝对可以出书了。

软件还处在很初级的阶段，谈不上出书

dannybear

请问一下，王老师，您使用的是资料是2A12的hdf.v7吗、

MeteoInfo

dannybear 发表于 2016-4-20 20:12
请问一下，王老师，您使用的是资料是2A12的hdf.v7吗、

这个记不清了

dannybear

MeteoInfo 发表于 2016-4-20 20:44
这个记不清了

哦哦，好的，谢谢王老师，是这样的，现在这个资料貌似变成hdf.V7，里面并没有geolocation这个变量

MeteoInfo

dannybear 发表于 2016-4-20 21:11
哦哦，好的，谢谢王老师，是这样的，现在这个资料貌似变成hdf.V7，里面并没有geolocation这个变量

你在Console里输入

>>> f

可以看到文件里的信息，包括变量名。

dannybear

MeteoInfo 发表于 2016-4-20 21:16
你在Console里输入

>>> f

恩恩，谢谢王老师