meteoinfo可以读取trmm的HDF降水资料么？

123木头兔

如题，我发现meteoinfo有读取HDF文件的选项，但是却无法读取trmm的HDF降水资料，不知为何，求解。谢谢

mofangbao

这种情况最好给个非保密数据的样例来测试

MeteoInfo

1、用MeteoInfo Java版
2、写脚本程序，示例如下：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-7-2
   
4. # Purpose: Read TRMM precipitation hdf4 data
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. from org.meteoinfo.data.meteodata import MeteoDataInfo
   
8. from org.meteoinfo.data.meteodata import Dimension
   
9. from org.meteoinfo.data.meteodata import DimensionType
   
10. from org.meteoinfo.data.meteodata import DrawMeteoData
    
11. from org.meteoinfo.legend import LegendType
    
12. from org.meteoinfo.legend import LegendManage
    
13. from org.meteoinfo.shape import ShapeTypes
    
14. from ucar.nc2 import NetcdfFile
    
15. import os.path
    
16. import jarray
    
17. 
    
18. #Set data directory
    
19. dataDir = 'D:/Temp/hdf'
    
20. 
    
21. #Create MeteoDataInfo object
    
22. mdi = MeteoDataInfo()
    
23. 
    
24. #Read hdf data file
    
25. fn = os.path.join(dataDir, '3B43.19980101.7.HDF')
    
26. if os.path.isfile(fn):
    
27.         print fn
    
28.         mdi.openNetCDFData(fn)
    
29.         dataInfo = mdi.getDataInfo()
    
30.         #print mdi.getInfoText()
    
31.         xmin = -180.0
    
32.         ymin = -50.0
    
33.         xnum = 1440
    
34.         ynum = 400
    
35.         xdelt = 0.25
    
36.         ydelt = 0.25
    
37.         xlist = []
    
38.         ylist = []
    
39.         for i in range(0,xnum):
    
40.                 xlist.append(xmin + xdelt * i)
    
41.         for i in range(0,ynum):
    
42.                 ylist.append(ymin + ydelt * i)
    
43. 
    
44.         X = jarray.array(xlist, 'd')
    
45.         Y = jarray.array(ylist, 'd')
    
46.         xDim = Dimension(DimensionType.X)
    
47.         xDim.setValues(X)
    
48.         dataInfo.setXDimension(xDim)
    
49.         yDim = Dimension(DimensionType.Y)
    
50.         yDim.setValues(Y)
    
51.         dataInfo.setYDimension(yDim)
    
52.         var = dataInfo.getVariable('precipitation')
    
53.         print var.getName()
    
54.         dimList = [xDim, yDim]
    
55.         var.setDimensions(dimList)
    
56.         gData = mdi.getGridData(var.getName())       
    
57.         aLS = LegendManage.createLegendSchemeFromGridData(gData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
58.         aLayer = DrawMeteoData.createRasterLayer(gData, "Test_HDF", aLS)
    
59.         mf = miapp.getMapDocument().getActiveMapFrame()
    
60.         mf.addLayer(aLayer)
    
61.         mf.moveLayer(aLayer, 0)       
    
62. 
    
63. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-7-2<br /> # Purpose: Read TRMM precipitation hdf4 data<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.data.meteodata import Dimension<br /> from org.meteoinfo.data.meteodata import DimensionType<br /> from org.meteoinfo.data.meteodata import DrawMeteoData<br /> from org.meteoinfo.legend import LegendType<br /> from org.meteoinfo.legend import LegendManage<br /> from org.meteoinfo.shape import ShapeTypes<br /> from ucar.nc2 import NetcdfFile<br /> import os.path<br /> import jarray<br /> <br /> #Set data directory<br /> dataDir = 'D:/Temp/hdf'<br /> <br /> #Create MeteoDataInfo object<br /> mdi = MeteoDataInfo()<br /> <br /> #Read hdf data file<br /> fn = os.path.join(dataDir, '3B43.19980101.7.HDF')<br /> if os.path.isfile(fn):<br />         print fn<br />         mdi.openNetCDFData(fn)<br />         dataInfo = mdi.getDataInfo()<br />         #print mdi.getInfoText()<br />         xmin = -180.0<br />         ymin = -50.0<br />         xnum = 1440<br />         ynum = 400<br />         xdelt = 0.25<br />         ydelt = 0.25<br />         xlist = []<br />         ylist = []<br />         for i in range(0,xnum):<br />                 xlist.append(xmin + xdelt * i)<br />         for i in range(0,ynum):<br />                 ylist.append(ymin + ydelt * i)<br /> <br />         X = jarray.array(xlist, 'd')<br />         Y = jarray.array(ylist, 'd')<br />         xDim = Dimension(DimensionType.X)<br />         xDim.setValues(X)<br />         dataInfo.setXDimension(xDim)<br />         yDim = Dimension(DimensionType.Y)<br />         yDim.setValues(Y)<br />         dataInfo.setYDimension(yDim)<br />         var = dataInfo.getVariable('precipitation')<br />         print var.getName()<br />         dimList = [xDim, yDim]<br />         var.setDimensions(dimList)<br />         gData = mdi.getGridData(var.getName())        <br />         aLS = LegendManage.createLegendSchemeFromGridData(gData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br />         aLayer = DrawMeteoData.createRasterLayer(gData, "Test_HDF", aLS)<br />         mf = miapp.getMapDocument().getActiveMapFrame()<br />         mf.addLayer(aLayer)<br />         mf.moveLayer(aLayer, 0)        <br /> <br /> print 'Finished!'

3、结果：

123木头兔

MeteoInfo 发表于 2014-7-2 23:46
1、用MeteoInfo Java版
2、写脚本程序，示例如下：

太好了，谢谢王老师。

123木头兔

mofangbao 发表于 2014-7-2 23:22
这种情况最好给个非保密数据的样例来测试

谢谢关注和提醒

letme

MeteoInfo 发表于 2014-7-2 23:46
1、用MeteoInfo Java版
2、写脚本程序，示例如下：

王老师，使用您编的脚本,运行显示
run script...
E:/SYY/TRMM\3B43.20010801.7A.HDF
precipitation
Traceback (most recent call last):
  File "<iostream>", line 47, in <module>
NameError: name 'miapp' is not defined

请问怎么解决呢？

MeteoInfo

letme 发表于 2016-8-7 22:10
王老师，使用您编的脚本,运行显示
run script...
E:/SYY/TRMM\3B43.20010801.7A.HDF

参考此贴：MeteoInfoLab脚本示例：TRMM 3B43 HDF数据
http://bbs.06climate.com/forum.p ... 579&fromuid=106
(出处: 气象家园)

呼小喵喵

王老师太棒了