MeteoInfoLab脚本示例：Himawari 8 真彩色图合成

MeteoInfo

这里示例用Himawari 8卫星的第1（蓝色）、2（绿色）和3（红色）波段的Albedo数据合成真彩色图，Albedo数据的范围应该是0-1，这里都出来的数据有极少数格点的值大于1，对此进行处理将其改为1。imshow和imshowm函数支持RBG(A)三维数组，或者R, G, B数组列表，将其合成为真彩色图。

1. #Add data file
   
2. fn = r'C:\Temp\himawari8\NC_H08_20170508_0040_r14_FLDK.02701_02601.nc'
   
3. f = addfile(fn)
   
4. #Read data
   
5. bdata = f['albedo_01'][:,:]
   
6. gdata = f['albedo_02'][:,:]
   
7. rdata = f['albedo_03'][:,:]
   
8. bdata[bdata>1] = 1
   
9. gdata[gdata>1] = 1
   
10. rdata[rdata>1] = 1
    
11. #Plot
    
12. axesm()
    
13. mlayer = shaperead('D:/Temp/map/country1.shp')
    
14. geoshow(mlayer, edgecolor='g')
    
15. layer = imshowm([rdata,gdata,bdata])
    
16. #Adjust image
    
17. imagelib.hsb_adjust(layer, h=0, s=0.1, b=0.2)
    
18. title('Himawari 8 true color image example')

#Add data file<br /> fn = r'C:\Temp\himawari8\NC_H08_20170508_0040_r14_FLDK.02701_02601.nc'<br /> f = addfile(fn)<br /> #Read data<br /> bdata = f['albedo_01'][:,:]<br /> gdata = f['albedo_02'][:,:]<br /> rdata = f['albedo_03'][:,:]<br /> bdata[bdata>1] = 1<br /> gdata[gdata>1] = 1<br /> rdata[rdata>1] = 1<br /> #Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer, edgecolor='g')<br /> layer = imshowm([rdata,gdata,bdata])<br /> #Adjust image<br /> imagelib.hsb_adjust(layer, h=0, s=0.1, b=0.2)<br /> title('Himawari 8 true color image example')

zxholystar

150度E图片和边框之间的白色间隙好解决吗？

MeteoInfo

zxholystar 发表于 2017-5-27 12:58
150度E图片和边框之间的白色间隙好解决吗？

xlim, ylim函数可以自定义显示范围

duduer

  王老师，请问出现这个错误是怎么回事呢？
File "<iostream>", line 15, in <module>
  File "D:\meteoinfo\MeteoInfo_Java_1.4.3R3_Files\MeteoInfo\pylib\mipylib\plotlib\miplot.py", line 4326, in imshowm
    ls = LegendManage.createImageLegend(gdata, cmap)
    at org.meteoinfo.legend.LegendManage.createImageLegend(LegendManage.java:1483)
    at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
    at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
    at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
    at java.lang.reflect.Method.invoke(Unknown Source)

java.lang.NullPointerException: java.lang.NullPointerException

出现在layer=imshowm这行

MeteoInfo

duduer 发表于 2017-11-24 23:27
王老师，请问出现这个错误是怎么回事呢？
File "", line 15, in
  File "D:\meteoinfo\MeteoInfo_Jav ...

建议你使用最新版本的MeteoInfo，如果还有问题再说。

duduer

MeteoInfo 发表于 2017-11-25 10:17
建议你使用最新版本的MeteoInfo，如果还有问题再说。

果然问题解决了，谢谢王老师！

duduer

MeteoInfo 发表于 2017-11-25 10:17
建议你使用最新版本的MeteoInfo，如果还有问题再说。

NameError: name 'imfilter' is not defined
王老师您好，请问这个错误怎么解决呢？

MeteoInfo

duduer 发表于 2017-11-25 17:42
NameError: name 'imfilter' is not defined
王老师您好，请问这个错误怎么解决呢？

imfilter改为imagelib

wuwei2163

王老师，一直是您忠实的粉丝，这两天再跟您学习葵花8号数据，在官网上看到同一时刻有三种数据，但后面的字母以及数字不明白是什么意思？？？
请王老师指导？？？

MeteoInfo

你下载数据的地方应该有文件名的说明文件，你仔细看看