那里示例用Himawari 8卫星的第1(蓝涩)、2(绿涩)和3(红涩)波段的Albedo数据分解实彩涩图,Albedo数据的领域应当是0-1,那里都出来的数据有少少数格点的值大于1,对此停行办理将其改为1。imshow和imshowm函数撑持RBG(A)三维数组,大概R, G, B数组列表,将其分解为实彩涩图。
#Add data file
fn = r'E:Temphimawari8NC_H08_20170508_0040_r14_FLDK.02701_02601.nc'
f = addfile(fn)
#Read data
bdata = f['albedo_01'][:,:]
gdata = f['albedo_02'][:,:]
rdata = f['albedo_03'][:,:]
bdata[bdata>1] = 1
gdata[gdata>1] = 1
rdata[rdata>1] = 1
#Plot
aVesm()
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer, edgecolor='g')
layer = imshowm([rdata,gdata,bdata])
#Adjust image
imfilter.hsb_adjust(layer, h=0, s=0.1, b=0.2)
title('Himawari 8 true color image eVample')