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

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()