||
这里示例用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
axesm()
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 example')
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-12-22 15:07
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社