MeteoInfoLab脚本示例：Hamawari-8 netCDF data

示例数据：ftp://ftp.bom.gov.au/anon/sample/catalogue/Satellite/IDE00220.201507140300.nc

 

该数据的分辨率很高（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)