标签：
FY-3C全球火点HDF数据包含一个FIRES二维变量，第一维是火点数，第二维是一些属性，其中第3、4列分别是火点的纬度和经度。下面的脚本示例读出所有火点经纬度并绘图。
脚本程序：
#Add data file fn = ‘D:/Temp/hdf/FY3C_VIRRX_GBAL_L2_GFR_MLT_GLL_20150811_POAD_1000M_MS.HDF‘ f = addfile(fn) #Get data variable v = f[‘FIRES‘] #Get data array data = v[:,5] lat = v[:,3] lon = v[:,4] #Plot axesm() world = shaperead(‘D:/Temp/map/country1.shp‘) china = shaperead(‘D:/Temp/map/bou2_4p.shp‘) geoshow(china, edgecolor=‘gray‘) geoshow(world, edgecolor=(100,100,100)) layer = scatterm(lon, lat, data, s=3, colors=[‘r‘], edge=False, marker=‘+‘) title(‘FY-3C GFR‘) axism()
MeteoInfoLab脚本示例：FY-3C全球火点HDF数据
标签：
原文地址：http://www.cnblogs.com/yaqiang/p/4728015.html
