MeteoInfo脚本示例：读取FY3A AOD HDF文件

MeteoInfo

FY3A卫星有AOD产品数据，HDF格式，这里示例用MeteoInfo脚本程序读取和显示该类数据。


脚本程序如下：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2015-3-18
   
4. # Purpose: Read FY3A AOD hdf5 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 LegendScheme
    
12. from org.meteoinfo.shape import ShapeTypes
    
13. from ucar.nc2 import NetcdfFile
    
14. import os.path
    
15. import jarray
    
16. 
    
17. #Set data directory
    
18. dataDir = 'D:/Temp/Hdf'
    
19. 
    
20. #Create MeteoDataInfo object
    
21. mdi = MeteoDataInfo()
    
22. 
    
23. #Read hdf data file
    
24. fn = os.path.join(dataDir, 'FY3A_MERSI_GBAL_L3_ASL_MLT_GLL_20140331_AOAM_5000M_MS.HDF')
    
25. if os.path.isfile(fn):
    
26.     print fn
    
27.     mdi.openNetCDFData(fn)
    
28.     dataInfo = mdi.getDataInfo()
    
29.     #print mdi.getInfoText()
    
30.     xmin = -30.0
    
31.     ymin = -90.0
    
32.     xnum = 7200
    
33.     ynum = 3600
    
34.     xdelt = 0.05
    
35.     ydelt = 0.05
    
36.     xlist = []
    
37.     ylist = []
    
38.     for i in range(0,xnum):
    
39.         xlist.append(xmin + xdelt * i)
    
40.     for i in range(0,ynum):
    
41.         ylist.append(ymin + ydelt * i)
    
42. 
    
43.     X = jarray.array(xlist, 'd')
    
44.     Y = jarray.array(ylist, 'd')
    
45.     xDim = Dimension(DimensionType.X)
    
46.     xDim.setValues(X)
    
47.     dataInfo.setXDimension(xDim)
    
48.     yDim = Dimension(DimensionType.Y)
    
49.     yDim.setValues(Y)
    
50.     dataInfo.setYDimension(yDim)
    
51.     var = dataInfo.getVariable('Aerosol_Optical_Thickness_of_MERSI_550nm')
    
52.     print var.getName()
    
53.     dimList = [yDim, xDim]
    
54.     var.setDimensions(dimList)
    
55.     gData = mdi.getGridData(var.getName())      
    
56.     gData.yReverse()
    
57.     gData = gData.div(10000.0)
    
58.     gData.missingValue = -3.2767
    
59.     #aLS = LegendManage.createLegendSchemeFromGridData(gData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
60.     aLS = LegendScheme(ShapeTypes.Polygon)
    
61.     aLS.importFromXMLFile(os.path.join(dataDir, 'AOD.lgs'))
    
62.     aLayer = DrawMeteoData.createRasterLayer(gData, "Test_HDF", aLS)
    
63.     mf = miapp.getMapDocument().getActiveMapFrame()
    
64.     mf.addLayer(aLayer)
    
65.     mf.moveLayer(aLayer, 0)
    
66. 
    
67. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2015-3-18<br /> # Purpose: Read FY3A AOD hdf5 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 LegendScheme<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, 'FY3A_MERSI_GBAL_L3_ASL_MLT_GLL_20140331_AOAM_5000M_MS.HDF')<br /> if os.path.isfile(fn):<br />     print fn<br />     mdi.openNetCDFData(fn)<br />     dataInfo = mdi.getDataInfo()<br />     #print mdi.getInfoText()<br />     xmin = -30.0<br />     ymin = -90.0<br />     xnum = 7200<br />     ynum = 3600<br />     xdelt = 0.05<br />     ydelt = 0.05<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('Aerosol_Optical_Thickness_of_MERSI_550nm')<br />     print var.getName()<br />     dimList = [yDim, xDim]<br />     var.setDimensions(dimList)<br />     gData = mdi.getGridData(var.getName())      <br />     gData.yReverse()<br />     gData = gData.div(10000.0)<br />     gData.missingValue = -3.2767<br />     #aLS = LegendManage.createLegendSchemeFromGridData(gData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br />     aLS = LegendScheme(ShapeTypes.Polygon)<br />     aLS.importFromXMLFile(os.path.join(dataDir, 'AOD.lgs'))<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!'

兰溪之水

王老师真及时，专门间接帮助本科生做论文

星期八

王老师，我依葫芦画瓢学着读取FY3A AOD HDF文件，我将老师的代码复制进去，修改了路径和名字，系统提示说找不到指定文件

D:\Temp\Hdf\FY3A_MERSI_GBAL_L3_ASL_MLT_GLL_20140430_AOAM_5000M_MS.HDF
Aerosol_Optical_Thickness_of_MERSI_550nm
Traceback (most recent call last):
  File "<iostream>", line 57, in <module>
at java.io.FileInputStream.open(Native Method)
at java.io.FileInputStream.<init>(Unknown Source)
at java.io.FileInputStream.<init>(Unknown Source)
at sun.net.www.protocol.file.FileURLConnection.connect(Unknown Source)
at sun.net.www.protocol.file.FileURLConnection.getInputStream(Unknown Source)
at com.sun.org.apache.xerces.internal.impl.XMLEntityManager.setupCurrentEntity(Unknown Source)
at com.sun.org.apache.xerces.internal.impl.XMLVersionDetector.determineDocVersion(Unknown Source)
at com.sun.org.apache.xerces.internal.parsers.XML11Configuration.parse(Unknown Source)
at com.sun.org.apache.xerces.internal.parsers.XML11Configuration.parse(Unknown Source)
at com.sun.org.apache.xerces.internal.parsers.XMLParser.parse(Unknown Source)
at com.sun.org.apache.xerces.internal.parsers.DOMParser.parse(Unknown Source)
at com.sun.org.apache.xerces.internal.jaxp.DocumentBuilderImpl.parse(Unknown Source)
at javax.xml.parsers.DocumentBuilder.parse(Unknown Source)
at org.meteoinfo.legend.LegendScheme.importFromXMLFile(LegendScheme.java:887)
at org.meteoinfo.legend.LegendScheme.importFromXMLFile(LegendScheme.java:872)
at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
at java.lang.reflect.Method.invoke(Unknown Source)
java.io.FileNotFoundException: java.io.FileNotFoundException: D:\Temp\Hdf\AOD.lgs (系统找不到指定的文件。)

想请教老师，是哪里出错了。

MeteoInfo

应该是没有AOD.lgs文件

星期八

谢谢老师，我从您其他的帖子里面下载到了

mzy1027

大家来顶帖啊。。。

zhangyanlin

请问哪里看一下到AOD.lgs文件

MeteoInfo

zhangyanlin 发表于 2015-4-6 15:21
请问哪里看一下到AOD.lgs文件

见此贴：http://bbs.06climate.com/forum.p ... 6&highlight=AOD