Modis FMF資料在MeteoInfo應用

MeteoInfo

KIMO23 发表于 2016-7-20 17:42
老师您好~~~又有问题要请教了，我打算要求取3笔Modis卫星数据的平均，所以撰写了下述的脚本，执行结果如 ...

v1 = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean']
改为
v1 = f1['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean']
v2也要改

这样平均是可以的，不过卫星数据有很多缺测值，这样平均的结果格点数据中的缺测值格点是所有参与平均的格点数据中的缺测值格点并集。也就是说参与平均的任一格点数据的缺测值格点在结果数据中也是缺测值格点。

还有一种方式能够很多程度上降低缺测值格点数，结果中某个格点只有所有参与平均的数据中该格点都是缺测值时才是缺测值格点。
dlist = [data, data1, data2]
dataavg = mean(dlist)

KIMO23

MeteoInfo 发表于 2016-7-20 18:07
v1 = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean']
改为
v1 = f1['Optical_Depth_Ratio_Small_ ...

感谢老师~~~刚刚修改了一下语法以及运行结果如下~~~~老师于脚本中有提到dlist = [data,data1, data2]  dataavg =mean(dlist)此种指令，能否告知我其意义及原理为和呢???谢谢您

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 = f1['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 = f2['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. dlist = [data, data1, data2]
    
24. dataavg = mean(dlist)
    
25. #dataavg=(data+data1+data2)/3
    
26. #Plot
    
27. #subplot(2,1,1)
    
28. axesm()
    
29. china = shaperead('C:/MeteoInfo_Java_1.3.4_Files/MeteoInfo/map/bou2_4p.shp')
    
30. geoshow(china)
    
31. world = shaperead('C:/MeteoInfo_Java_1.3.4_Files/MeteoInfo/map/country1.shp')
    
32. geoshow(world, edgecolor='k')
    
33. levels = arange(0, 1., 0.05)
    
34. layer = imshowm(dataavg, levels)
    
35. #geoshow(lat, lon, size=2, color='k')
    
36. #title('Ratio of small mode aerosol optical depth at 0.55 micron')
    
37. title('Ratio of small mode aerosol optical depth at 0.55 micron')
    
38. colorbar(layer)

复制代码

MeteoInfo

KIMO23 发表于 2016-7-20 20:34
感谢老师~~~刚刚修改了一下语法以及运行结果如下~~~~老师于脚本中有提到dlist = [data,data1, data2]  da ...

参考此贴：MeteoInfo脚本示例：格点数据平均
http://bbs.06climate.com/forum.p ... 604&fromuid=106
(出处: 气象家园)