Modis FMF資料在MeteoInfo應用

KIMO23

想请教老师几个Modis在MeteoInfo应用的问题，

(1)您有用过modis的（气溶胶光学厚度）细粒子比FMF资料？这个是我在论坛看到网友下载的地址ftp://ladsftp.nascom.nasa.gov/allData/51/MOD08_M3/
网友也提到文献中这个资料都是指modis的三级产品中的月平均数据中就有，也就是在MOD08_M3里面可以找到(optical_depth_ratio_small_land_and_ocean_mean_mean这个参数似乎就是要的FMF（fine modefraction）)。不知道您能否针对这个变量的分析撰写MeteoInfoLab相关分析脚本吗???

(2)该数据分析完毕后MeteoInfo软件能绘制该数据剖面折线图吗??例如说我固定经度(110 E)后显示各纬度(Modis卫星网格分辨率)上的变量值???

                                                                  

                                                                                 感谢您

MeteoInfo

参考下面的代码：

1. f = addfile('D:/Temp/hdf/MOD08_M3.A2010032.005.2010062234836.hdf')
   
2. v = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean_Mean']
   
3. data = v[::-1,:]
   
4. data[data==-9.999] = nan
   
5. dd = data[:,[110]]
   
6. #Plot
   
7. subplot(2,1,1)
   
8. axesm()
   
9. china = shaperead('D:/Temp/map/bou2_4p.shp')
   
10. geoshow(china, linecolor='gray')
    
11. world = shaperead('D:/Temp/map/country1.shp')
    
12. geoshow(world)
    
13. levels = arange(0, 1.0, 0.05)
    
14. layer = imshowm(data, levels)
    
15. title('Ratio of small mode aerosol optical depth at 0.55 micron')
    
16. colorbar(layer)
    
17. 
    
18. subplot(2,1,2)
    
19. plot(dd)
    
20. xlabel('Latitude')
    
21. title('Value along longitude of 110E')

f = addfile('D:/Temp/hdf/MOD08_M3.A2010032.005.2010062234836.hdf')<br /> v = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean_Mean']<br /> data = v[::-1,:]<br /> data[data==-9.999] = nan<br /> dd = data[:,[110]]<br /> #Plot<br /> subplot(2,1,1)<br /> axesm()<br /> china = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(china, linecolor='gray')<br /> world = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(world)<br /> levels = arange(0, 1.0, 0.05)<br /> layer = imshowm(data, levels)<br /> title('Ratio of small mode aerosol optical depth at 0.55 micron')<br /> colorbar(layer)<br /> <br /> subplot(2,1,2)<br /> plot(dd)<br /> xlabel('Latitude')<br /> title('Value along longitude of 110E')

KIMO23

MeteoInfo 发表于 2016-7-14 13:56
参考下面的代码：

谢谢老师的分享~~~学习了~~~~

兰溪之水

MeteoInfo 发表于 2016-7-14 13:56
参考下面的代码：

王老师又热情帮助台湾同胞解决问题

MeteoInfo

兰溪之水 发表于 2016-7-15 00:29
王老师又热情帮助台湾同胞解决问题

用几个繁体字就是台湾同胞了？

KIMO23

兰溪之水 发表于 2016-7-15 00:29
王老师又热情帮助台湾同胞解决问题

兰溪版主实在眼利，知道我来自台湾，所以在用语问话方面上不到位的地方请各位多海涵，尔后请大家多多指教，一同交流学习。

KIMO23

MeteoInfo 发表于 2016-7-14 13:56
参考下面的代码：

老师，不好意思再请教一些问题。
(1)若今天是想要绘制一个跨经纬度的曲线剖面图(例如:20N,100E~30N,120E),是要如何dd=data[:,[110]]这个部份上撰写???
(2)若曲线图上我想要绘制1条以上的曲线，脚本将要如何宣告??
                                                                                                                             感谢您

兰溪之水

MeteoInfo 发表于 2016-7-15 08:33
用几个繁体字就是台湾同胞了？

哈哈，还有IP佐证

931404656

都是中华儿女！！！！

MeteoInfo

KIMO23 发表于 2016-7-15 11:37
老师，不好意思再请教一些问题。
(1)若今天是想要绘制一个跨经纬度的曲线剖面图(例如:20N,100E~30N,120E ...

有点复杂，仔细看看下面的代码。

1. f = addfile('D:/Temp/hdf/MOD08_M3.A2010032.005.2010062234836.hdf')
   
2. v = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean_Mean']
   
3. data = v[::-1,:]
   
4. data = data[[10,40],[90,130]]
   
5. data[data==-9.999] = nan
   
6. 
   
7. lon = [100,120]
   
8. lat = [20,30]
   
9. 
   
10. #Plot
    
11. subplot(2,1,1)
    
12. axesm()
    
13. china = shaperead('D:/Temp/map/bou2_4p.shp')
    
14. geoshow(china)
    
15. world = shaperead('D:/Temp/map/country1.shp')
    
16. geoshow(world, edgecolor='k')
    
17. levels = arange(0, 1.0, 0.05)
    
18. layer = imshowm(data, levels)
    
19. geoshow(lat, lon, size=2, color='k')
    
20. title('Ratio of small mode aerosol optical depth at 0.55 micron')
    
21. colorbar(layer)
    
22. 
    
23. subplot(2,1,2)
    
24. slope, intercept, r = linregress(lon, lat)[:3]
    
25. x = arange(100, 120, 1)
    
26. y = polyval([slope, intercept], x)
    
27. dd = data.tostation(x,y)
    
28. plot(dd)
    
29. xts = []
    
30. for i in range(0, len(x)):
    
31.     xts.append('(%.1fE,%.1fN)' % (x, y))
    
32. xticks(arange(len(x)), xts)
    
33. title('Value between (100E, 20N) and (120E, 30N)')

f = addfile('D:/Temp/hdf/MOD08_M3.A2010032.005.2010062234836.hdf')<br /> v = f['Optical_Depth_Ratio_Small_Land_And_Ocean_Mean_Mean']<br /> data = v[::-1,:]<br /> data = data[[10,40],[90,130]]<br /> data[data==-9.999] = nan<br /> <br /> lon = [100,120]<br /> lat = [20,30]<br /> <br /> #Plot<br /> subplot(2,1,1)<br /> axesm()<br /> china = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(china)<br /> world = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(world, edgecolor='k')<br /> levels = arange(0, 1.0, 0.05)<br /> layer = imshowm(data, levels)<br /> geoshow(lat, lon, size=2, color='k')<br /> title('Ratio of small mode aerosol optical depth at 0.55 micron')<br /> colorbar(layer)<br /> <br /> subplot(2,1,2)<br /> slope, intercept, r = linregress(lon, lat)[:3]<br /> x = arange(100, 120, 1)<br /> y = polyval([slope, intercept], x)<br /> dd = data.tostation(x,y)<br /> plot(dd)<br /> xts = []<br /> for i in range(0, len(x)):<br />     xts.append('(%.1fE,%.1fN)' % (x<i>, y<i>))<br /> xticks(arange(len(x)), xts)<br /> title('Value between (100E, 20N) and (120E, 30N)')