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

转载于:https://www.cnblogs.com/yaqiang/p/4728015.html