用MeteoInfo Java版尝试绘制FY3B火点数据

MeteoInfo

FY3B火点数据是HDF5格式，用了一个2维数据变量来表示火点，用HDFView打开数据如下：


由于数据本身的自描述信息比较缺乏，目前也没有该数据的详细说明，只能从数据本身猜测第3列和第4列分别为火点的纬度和经度。基于上述猜测写了一个脚本来读取火点数据并绘制分布图。


脚本代码：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2013-10-30
   
4. # Purpose: Read FY3B fire point 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.data import StationData
    
12. from org.meteoinfo.projection import ProjectionInfo
    
13. from org.meteoinfo.legend import LegendManage
    
14. from org.meteoinfo.legend import LegendType
    
15. from org.meteoinfo.shape import ShapeTypes
    
16. import os.path
    
17. import jarray
    
18. 
    
19. #Set data directory
    
20. dataDir = 'D:/Temp/hdf/'
    
21. 
    
22. #Create MeteoDataInfo object
    
23. mdi = MeteoDataInfo()
    
24. 
    
25. #Read hdf data file
    
26. fn = dataDir + 'FY3B_VIRRX_GBAL_L2_GFR_MLT_GLL_20130529_POAD_1000M_MS.HDF'
    
27. if os.path.isfile(fn):
    
28.         print fn
    
29.         mdi.openNetCDFData(fn)
    
30.         dataInfo = mdi.getDataInfo()
    
31.         print mdi.getInfoText()
    
32.         xnum = 9
    
33.         ynum = 8338
    
34.         xlist = []
    
35.         ylist = []
    
36.         for i in range(0,xnum):
    
37.                 xlist.append(i)
    
38.         for i in range(0,ynum):
    
39.                 ylist.append(i)
    
40. 
    
41.         X = jarray.array(xlist, 'd')
    
42.         Y = jarray.array(ylist, 'd')
    
43.         xDim = Dimension(DimensionType.X)
    
44.         xDim.setValues(X)
    
45.         dataInfo.setXDimension(xDim)
    
46.         yDim = Dimension(DimensionType.Y)
    
47.         yDim.setValues(Y)
    
48.         dataInfo.setYDimension(yDim)
    
49.         var = dataInfo.getVariable('FIRES')
    
50.         dimList = [yDim, xDim]
    
51.         var.setDimensions(dimList)
    
52.         gData = mdi.getGridData(var.getName())        
    
53.         sData = StationData()
    
54.         for d in range(0,ynum):
    
55.                 stid = str(d + 1)
    
56.                 lat = gData.data[d][3]
    
57.                 lon = gData.data[d][4]
    
58.                 value = gData.data[d][5]
    
59.                 sData.addData(stid, lon, lat, value)
    
60.                
    
61.         aLS = LegendManage.createLegendSchemeFromStationData(sData, LegendType.GraduatedColor, ShapeTypes.Point)
    
62.         aLayer = DrawMeteoData.createSTPointLayer(sData, aLS, 'FY3B_Fire', 'Fire')
    
63.         #aLayer.setProjInfo(dataInfo.getProjectionInfo())
    
64.         mf = miapp.getMapDocument().getActiveMapFrame()
    
65.         mf.addLayer(aLayer)
    
66.         #mf.moveLayer(aLayer, 0)
    
67. 
    
68. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2013-10-30<br /> # Purpose: Read FY3B fire point 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.data import StationData<br /> from org.meteoinfo.projection import ProjectionInfo<br /> from org.meteoinfo.legend import LegendManage<br /> from org.meteoinfo.legend import LegendType<br /> from org.meteoinfo.shape import ShapeTypes<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 = dataDir + 'FY3B_VIRRX_GBAL_L2_GFR_MLT_GLL_20130529_POAD_1000M_MS.HDF'<br /> if os.path.isfile(fn):<br />         print fn<br />         mdi.openNetCDFData(fn)<br />         dataInfo = mdi.getDataInfo()<br />         print mdi.getInfoText()<br />         xnum = 9<br />         ynum = 8338<br />         xlist = []<br />         ylist = []<br />         for i in range(0,xnum):<br />                 xlist.append(i)<br />         for i in range(0,ynum):<br />                 ylist.append(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('FIRES')<br />         dimList = [yDim, xDim]<br />         var.setDimensions(dimList)<br />         gData = mdi.getGridData(var.getName())        <br />         sData = StationData()<br />         for d in range(0,ynum):<br />                 stid = str(d + 1)<br />                 lat = gData.data[d][3]<br />                 lon = gData.data[d][4]<br />                 value = gData.data[d][5]<br />                 sData.addData(stid, lon, lat, value)<br />                 <br />         aLS = LegendManage.createLegendSchemeFromStationData(sData, LegendType.GraduatedColor, ShapeTypes.Point)<br />         aLayer = DrawMeteoData.createSTPointLayer(sData, aLS, 'FY3B_Fire', 'Fire')<br />         #aLayer.setProjInfo(dataInfo.getProjectionInfo())<br />         mf = miapp.getMapDocument().getActiveMapFrame()<br />         mf.addLayer(aLayer)<br />         #mf.moveLayer(aLayer, 0)<br /> <br /> print 'Finished!'

Mc.Fish

非常感谢王老师.
数据是在卫星中心数据网上下的.
业务急需呀.解了燃眉之急.

topmad

串了一个字  看成了点火数据...突然感觉高大上