MeteoInfoLab脚本示例：TRMM 3B43 HDF数据

MeteoInfo

TRMM 3B43是卫星观测月平均降水量产品，是HDF的格点数据。需要注意的是数据中降水变量维的顺序里经度维在前纬度维在后，这与通常的设置（纬度维在前经度维在后）相反，需要对获取的二维数组进行转置，使用transpose函数。数据的经纬度数组也需要定义（文件中没有）。数据中的0值为缺测值（fill_value），并不是一个好的缺测值设置，因此改为-9999.0。

脚本程序：

1. #Add data file
   
2. folder = 'D:/Temp/hdf/'
   
3. fns = '3B43.100301.6A.HDF'
   
4. fn = folder + fns
   
5. f = addfile(fn)
   
6. vname = 'precipitation'
   
7. t = f[vname]
   
8. rain = t[0,:,:]
   
9. rain = transpose(rain)
   
10. rain[rain==0] = -9999.0
    
11. rain.fill_value = -9999.0
    
12. lat = arange(-49.875, 49.875, 0.249375)
    
13. lon = arange(-179.875, 179.876, 0.25)
    
14. #Plot
    
15. axesm()
    
16. mlayer = shaperead('D:/Temp/map/country1.shp')
    
17. geoshow(mlayer)
    
18. layer = imshowm(lon, lat, rain, 20)
    
19. colorbar(layer, orientation='horizontal')
    
20. title(vname + 'at scan=0')
    
21. axism()

#Add data file<br /> folder = 'D:/Temp/hdf/'<br /> fns = '3B43.100301.6A.HDF'<br /> fn = folder + fns<br /> f = addfile(fn)<br /> vname = 'precipitation'<br /> t = f[vname]<br /> rain = t[0,:,:]<br /> rain = transpose(rain)<br /> rain[rain==0] = -9999.0<br /> rain.fill_value = -9999.0<br /> lat = arange(-49.875, 49.875, 0.249375)<br /> lon = arange(-179.875, 179.876, 0.25)<br /> #Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> layer = imshowm(lon, lat, rain, 20)<br /> colorbar(layer, orientation='horizontal')<br /> title(vname + 'at scan=0')<br /> axism()

a1s1

赞，谢谢分享  

喃音

{:5_235:}围观

先驱vx

你好，数据中的0值为缺测值（fill_value），这个有官方的说明文档吗？
我在数据hdf的元数据内没有看见有描述，找了下官网，http://disc.sci.gsfc.nasa.gov/ad ... n_faq.shtml#v7_chgs
有如下，
What are the values for latitude and longitude in TRMM 3B42 and 3B43 HDF data files?

The data in the HDF file are written in the following order, (-49.875, -179.875), (-49.625, -179.875), (-49.375, -179.875), (-49.125, -179.875)......(49.875, -179.875), (-49.875, -179.625), (-49.625, -179.625)...... The dimension for longitude is 1440 and the dimension for latitude is 400. The value used to indicate missing data is "-9999.9".

缺省为0的话，应该不明智，这样就区分不了无降水了。

MeteoInfo

先驱vx 发表于 2015-9-10 18:45
你好，数据中的0值为缺测值（fill_value），这个有官方的说明文档吗？
我在数据hdf的元数据内没有看见有描 ...

这个数据文件是从HDFEOS网站下载的，参考网站上MatLab代码中将0作为缺测值。也许数据官网有改进吧，没研究过。

先驱vx

MeteoInfo 发表于 2015-9-10 20:25
这个数据文件是从HDFEOS网站下载的，参考网站上MatLab代码中将0作为缺测值。也许数据官网有改进吧，没研 ...

好的，谢谢！

梁的丰

这是不是说数据不能在meteoinfo直接打开和绘图，也不能直接插值到站点了？

MeteoInfo

梁的丰 发表于 2015-9-28 17:05
这是不是说数据不能在meteoinfo直接打开和绘图，也不能直接插值到站点了？

可以用MeteoInfoLab脚本处理。

梁的丰

MeteoInfo 发表于 2015-9-28 19:12
可以用MeteoInfoLab脚本处理。

因为之前看过老师的帖子（http://bbs.06climate.com/forum.php?mod=viewthread&tid=1791），发现很多格式的文件在meteoinfo里都可以很方便的显示，并插值到站点，所以来问问是不是最新版的meteoinfo支持TRMM 3B43  HDF的直接显示，谢谢老师！！！

MeteoInfo

梁的丰 发表于 2015-9-28 19:22
因为之前看过老师的帖子（http://bbs.06climate.com/forum.php?mod=viewthread&tid=1791），发现很多格式 ...

象netCDF和HDF之类的数据格式非常灵活，在MeteoInfo桌面软件中直接打开需要遵循某种约定，比如CF Conventions。否则就只能通过程序来处理。