Modis FMF資料在MeteoInfo應用

KIMO23

MeteoInfo 发表于 2016-7-15 15:55
有点复杂，仔细看看下面的代码。

感谢老师，MeteoInfo软件实在强大，不知能否绘制一条以上的剖面曲线在同一张图上呢???此外MeteoInfoLab所撰写的语法是与Matlab相似吗???

MeteoInfo

KIMO23 发表于 2016-7-16 10:47
感谢老师，MeteoInfo软件实在强大，不知能否绘制一条以上的剖面曲线在同一张图上呢???此外MeteoInfoLab所 ...

MeteoInfo网站上和论坛里有一些脚本示例程序，里面有相关代码。
http://www.meteothinker.com/examples/index.html
http://bbs.06climate.com/forum.p ... &extra=page%3D1

很多函数和MatLab类似，语法是Python。

KIMO23

MeteoInfo 发表于 2016-7-16 20:14
MeteoInfo网站上和论坛里有一些脚本示例程序，里面有相关代码。
http://www.meteothinker.com/examples/ ...

感谢老师的分享与指导

KIMO23

MeteoInfo 发表于 2016-7-15 15:55
有点复杂，仔细看看下面的代码。

老师您好
      又再打扰了，我引用您的叙述脚本(版本为1.3.4)结果执行时发生错误讯息如图1，于是我将您的叙述脚本中的第30-31行移除后就能执行如图2，不适是哪边出问题呢????

MeteoInfo

KIMO23 发表于 2016-7-18 15:45
老师您好
      又再打扰了，我引用您的叙述脚本(版本为1.3.4)结果执行时发生错误讯息如图1，于是我将您 ...

KIMO23

MeteoInfo 发表于 2016-7-18 16:24

老师  感谢您  OK了，我想再请教一下为何我的2个图形显示大小不一呢????  {:eb502:}

MeteoInfo

KIMO23 发表于 2016-7-18 16:36
老师  感谢您  OK了，我想再请教一下为何我的2个图形显示大小不一呢????

对于地图缩放时需要保持它的纵横比例一致，其它图形可以随着缩放变形。所以会不一样，你可以把图形窗口点出来成为单独一个窗口（图形窗口右上角中间的按钮），然后通过缩放调整。

KIMO23

MeteoInfo 发表于 2016-7-18 16:48
对于地图缩放时需要保持它的纵横比例一致，其它图形可以随着缩放变形。所以会不一样，你可以把图形窗口点 ...

感谢老师研发与无私分享的这套软件~~~

MeteoInfo

KIMO23 发表于 2016-7-18 16:58
感谢老师研发与无私分享的这套软件~~~

另外还可以用figure函数固定图形的大小，参考此贴：MeteoInfoLab脚本示例：图形版面、点标注
http://bbs.06climate.com/forum.p ... 624&fromuid=106
(出处: 气象家园)

KIMO23

MeteoInfo 发表于 2016-7-18 18:32
另外还可以用figure函数固定图形的大小，参考此贴：MeteoInfoLab脚本示例：图形版面、点标注
http://bbs ...

老师您好~~~又有问题要请教了，我打算要求取3笔Modis卫星数据的平均，所以撰写了下述的脚本，执行结果如图，但不知道这样是否正确敬请老师指导一下，谢谢您

1. #Read data1
   
2. f = addfile('C:/MeteoInfo_Java_1.3.4_Files/MeteoInfo/sample/HDF/Modis/MOD08_D3.A2014017.051.2014018224401.hdf')
   
3. v = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean']
   
4. data = v[::-1,:]
   
5. data = data[[15,35],[105,135]]
   
6. data[data==-9.999] = nan
   
7. #Read data-2
   
8. f1 = addfile('C:/MeteoInfo_Java_1.3.4_Files/MeteoInfo/sample/HDF/Modis/MOD08_D3.A2014018.051.2014019101002.hdf')
   
9. v1 = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean']
   
10. data1 = v1[::-1,:]
    
11. data1 = data1[[15,35],[105,135]]
    
12. data1[data1==-9.999] = nan
    
13. #Read data-3
    
14. f2 = addfile('C:/MeteoInfo_Java_1.3.4_Files/MeteoInfo/sample/HDF/Modis/MOD08_D3.A2014019.051.2014020090624.hdf')
    
15. v2 = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean']
    
16. data2 = v2[::-1,:]
    
17. data2 = data2[[15,35],[105,135]]
    
18. data2[data2==-9.999] = nan
    
19. 
    
20. lon = [100,120]
    
21. lat = [20,30]
    
22. #Calculate
    
23. dataavg=(data+data1+data2)/3
    
24. #Plot
    
25. #subplot(2,1,1)
    
26. axesm()
    
27. china = shaperead('C:/MeteoInfo_Java_1.3.4_Files/MeteoInfo/map/bou2_4p.shp')
    
28. geoshow(china)
    
29. world = shaperead('C:/MeteoInfo_Java_1.3.4_Files/MeteoInfo/map/country1.shp')
    
30. geoshow(world, edgecolor='k')
    
31. levels = arange(0, 1., 0.05)
    
32. layer = imshowm(dataavg, levels)
    
33. #geoshow(lat, lon, size=2, color='k')
    
34. #title('Ratio of small mode aerosol optical depth at 0.55 micron')
    
35. title('Ratio of small mode aerosol optical depth at 0.55 micron(19)')
    
36. colorbar(layer)

复制代码