智能书签

那里示例用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')