MeteoInfoLab脚本示例：FY-3C全球火点HDF数据

MeteoInfo

FY-3C全球火点HDF数据包含一个FIRES二维变量，第一维是火点数，第二维是一些属性，其中第3、4列分别是火点的纬度和经度。下面的脚本示例读出所有火点经纬度并绘图。

脚本程序：

1. #Add data file
   
2. fn = 'D:/Temp/hdf/FY3C_VIRRX_GBAL_L2_GFR_MLT_GLL_20150811_POAD_1000M_MS.HDF'
   
3. f = addfile(fn)
   
4. #Get data variable
   
5. v = f['FIRES']
   
6. #Get data array
   
7. data = v[:,5]
   
8. lat = v[:,3]
   
9. lon = v[:,4]
   
10. #Plot
    
11. axesm()
    
12. world = shaperead('D:/Temp/map/country1.shp')
    
13. china = shaperead('D:/Temp/map/bou2_4p.shp')
    
14. geoshow(china, edgecolor='gray')
    
15. geoshow(world, edgecolor=(100,100,100))
    
16. layer = scatterm(lon, lat, data, s=3, colors=['r'], edge=False, marker='+')
    
17. title('FY-3C GFR')
    
18. axism()

#Add data file <br /> fn = 'D:/Temp/hdf/FY3C_VIRRX_GBAL_L2_GFR_MLT_GLL_20150811_POAD_1000M_MS.HDF' <br /> f = addfile(fn) <br /> #Get data variable <br /> v = f['FIRES'] <br /> #Get data array <br /> data = v[:,5] <br /> lat = v[:,3] <br /> lon = v[:,4] <br /> #Plot <br /> axesm() <br /> world = shaperead('D:/Temp/map/country1.shp') <br /> china = shaperead('D:/Temp/map/bou2_4p.shp') <br /> geoshow(china, edgecolor='gray') <br /> geoshow(world, edgecolor=(100,100,100)) <br /> layer = scatterm(lon, lat, data, s=3, colors=['r'], edge=False, marker='+') <br /> title('FY-3C GFR') <br /> axism()

午后阳光

学无止境！谢谢楼主分享！

DDQZXB

这几天深受王老师教诲，受益匪浅！
能否向王老师请教，如何使用Meteoinfo读取FY4A沙尘监测的Netcode格式数据？
学生在代码方面没有经验，一无所知！
拜谢！