MeteoInfoLab脚本示例：CALIPSO三维曲面图

MeteoInfo

用surf函数，facecolor='texturemap' ，结合 cdata 数据参数可以快速绘制3维坐标系中的2维曲面图。

1. f = addfile('D:/Temp/satellite/calipso/CAL_LID_L2_05kmAPro-Standard-V4-20.2006-09-01T08-55-57ZD_Subset.hdf')
   
2. EC = f['Extinction_Coefficient_532'][:]
   
3. EC = EC.T
   
4. EC[EC==nan] = -1
   
5. Lat = f['Latitude'][:,0]
   
6. Lon = f['Longitude'][:,0]
   
7. 
   
8. alt1 = arange(-0.5,20.08,0.06).tolist()
   
9. alt2 = arange1(20.2,56,0.18).tolist()
   
10. alt = array(alt1 + alt2)
    
11. 
    
12. Lon = Lon[np.newaxis,:]
    
13. Lon = Lon.repeat(len(alt),0)
    
14. Lat = Lat[np.newaxis,:]
    
15. Lat = Lat.repeat(len(alt),0)
    
16. alt = alt[:,np.newaxis]
    
17. alt = alt.repeat(Lon.shape[1],1)
    
18. 
    
19. fn = os.path.join(migl.get_map_folder(), 'world.topo.bathy.jpg')
    
20. land = georead(fn)
    
21. 
    
22. #plot
    
23. levs = arange(0,0.41,0.01)
    
24. cols = makecolors(len(levs)+1, cmap='rainbow')
    
25. cols[0] = miutil.getcolor(cols[0], alpha=0.5)
    
26. ax = axes3d()
    
27. grid(False)
    
28. #lighting()
    
29. geoshow(land)
    
30. surf(Lon[::10,::10], Lat[::10,::10], alt[::10,::10], levs, colors=cols,
    
31.     facecolor='texturemap', cdata=EC, edgecolor=None)
    
32. colorbar(aspect=30, ticks=[0,0.1,0.2,0.3,0.4,0.5])
    
33. zlim(0,8)
    
34. xlim(-100,140)
    
35. ylim(-90,90)
    
36. xlabel('Longitue')
    
37. ylabel('Latitude')
    
38. zlabel('Height (km)')

f = addfile('D:/Temp/satellite/calipso/CAL_LID_L2_05kmAPro-Standard-V4-20.2006-09-01T08-55-57ZD_Subset.hdf')<br /> EC = f['Extinction_Coefficient_532'][:]<br /> EC = EC.T<br /> EC[EC==nan] = -1<br /> Lat = f['Latitude'][:,0]<br /> Lon = f['Longitude'][:,0]<br /> <br /> alt1 = arange(-0.5,20.08,0.06).tolist()<br /> alt2 = arange1(20.2,56,0.18).tolist()<br /> alt = array(alt1 + alt2)<br /> <br /> Lon = Lon[np.newaxis,:]<br /> Lon = Lon.repeat(len(alt),0)<br /> Lat = Lat[np.newaxis,:]<br /> Lat = Lat.repeat(len(alt),0)<br /> alt = alt[:,np.newaxis]<br /> alt = alt.repeat(Lon.shape[1],1)<br /> <br /> fn = os.path.join(migl.get_map_folder(), 'world.topo.bathy.jpg')<br /> land = georead(fn)<br /> <br /> #plot<br /> levs = arange(0,0.41,0.01)<br /> cols = makecolors(len(levs)+1, cmap='rainbow')<br /> cols[0] = miutil.getcolor(cols[0], alpha=0.5)<br /> ax = axes3d()<br /> grid(False)<br /> #lighting()<br /> geoshow(land)<br /> surf(Lon[::10,::10], Lat[::10,::10], alt[::10,::10], levs, colors=cols,<br />     facecolor='texturemap', cdata=EC, edgecolor=None)<br /> colorbar(aspect=30, ticks=[0,0.1,0.2,0.3,0.4,0.5])<br /> zlim(0,8)<br /> xlim(-100,140)<br /> ylim(-90,90)<br /> xlabel('Longitue')<br /> ylabel('Latitude')<br /> zlabel('Height (km)')

                               
登录/注册后可看大图

MeteoInfo

Lizhengpeng 发表于 2022-9-5 09:54
感谢楼主的分享，请问如果要在图里加上后向轨迹该怎么做呢？

参考这里：http://www.meteothink.org/exampl ... tory/traj_plot.html

Lizhengpeng

感谢楼主的分享，请问如果要在图里加上后向轨迹该怎么做呢？

Lizhengpeng

MeteoInfo 发表于 2022-9-5 12:00
参考这里：http://www.meteothink.org/examples/meteoinfolab/trajectory/traj_plot.html

谢谢老师！刚刚接触MeteoInfo,前两行读取的轨迹文件也是MetoInfo计算生成的吗？像这样吗？
http://www.meteothink.org/exampl ... ctory/traj_cal.html