MeteoInfoLab脚本示例：Hamawari-8 netCDF data

示例数据：

该数据的分辨率很高（22000*22000），如果全部读入一个数组需要很大的内存，可以在读入时设置step（比如4，以间隔4读取数据），或者限定更小的空间范围。

脚本程序：

#Add data file

fn = 'D:/Temp/nc/IDE00220.201507140300.nc'

f = addfile(fn)

#Get data variable

v = f['channel_0003_brf']

data = v[0,::4,::4]

data = data[::-1,:]

#Plot

ax, proj = axesm(proj='geos', lon_0=104.7, h=35785863, gridlabel=True, gridline=True, frameon=False)

mlayer = shaperead('D:/Temp/map/country1.shp')

geoshow(mlayer)

levs = arange(0, 1, 0.1)

layer = imshowm(data, levs, proj=proj)

colorbar(layer)