MeteoInfoLab脚本示例：3维流线

MeteoInfo

streamplot3函数可以追踪并绘制三维流线，也可以自定义流线追踪的起始点。

1. # Make the grid
   
2. x, y, z = meshgrid(arange(-1.5, 1.5, 0.1),
   
3.                    arange(-1.5, 1.5, 0.1),
   
4.                    arange(-1.5, 1.5, 0.1))
   
5. 
   
6. # Make the direction data for the arrows
   
7. u = x + cos(4*x) + 3         # x-component of vector field
   
8. v = sin(4*x) - sin(2*y)      # y-component of vector field
   
9. w = -z                       # z-component of vector field
   
10. speed = sqrt(u*u + v*v + w*w)
    
11. 
    
12. sx, sy, sz = meshgrid([-1.5], [-1,0,1], [-1,0,1])
    
13. 
    
14. qq = streamplot3(x, y, z, u, v, w, speed, linewidth=2,
    
15.     density=4, interval=10, start_x=sx, start_y=sy, start_z=sz)
    
16. scatter3(sx, sy, sz, c='k')
    
17. colorbar(qq)
    
18. xlim(-1.5, 1.5)
    
19. ylim(-1.5, 1.5)
    
20. zlim(-1.5, 1.5)
    
21. antialias(True)

# Make the grid<br /> x, y, z = meshgrid(arange(-1.5, 1.5, 0.1),<br />                    arange(-1.5, 1.5, 0.1),<br />                    arange(-1.5, 1.5, 0.1))<br /> <br /> # Make the direction data for the arrows<br /> u = x + cos(4*x) + 3         # x-component of vector field<br /> v = sin(4*x) - sin(2*y)      # y-component of vector field<br /> w = -z                       # z-component of vector field<br /> speed = sqrt(u*u + v*v + w*w)<br /> <br /> sx, sy, sz = meshgrid([-1.5], [-1,0,1], [-1,0,1])<br /> <br /> qq = streamplot3(x, y, z, u, v, w, speed, linewidth=2,<br />     density=4, interval=10, start_x=sx, start_y=sy, start_z=sz)<br /> scatter3(sx, sy, sz, c='k')<br /> colorbar(qq)<br /> xlim(-1.5, 1.5)<br /> ylim(-1.5, 1.5)<br /> zlim(-1.5, 1.5)<br /> antialias(True)

streamslice函数可以在三维数据的某个轴切片上绘制流线图。

1. # Make the grid
   
2. x, y, z = meshgrid(arange(-1.5, 1.6, 0.1),
   
3.                    arange(-1.5, 1.6, 0.1),
   
4.                    arange(-1.5, 1.6, 0.1))
   
5. 
   
6. # Make the direction data for the arrows
   
7. u = x + cos(4*x) + 3         # x-component of vector field
   
8. v = sin(4*x) - sin(2*y)      # y-component of vector field
   
9. w = -z                       # z-component of vector field
   
10. speed = sqrt(u*u + v*v + w*w)
    
11. 
    
12. streamslice(x, y, z, u, v, w, xslice=1.5, yslice=1.5, zslice=-1.5,
    
13.     color='b', linewidth=1, density=4, interval=5)
    
14. xlim(-1.5, 1.5)
    
15. ylim(-1.5, 1.5)
    
16. zlim(-1.5, 1.5)
    
17. 
    
18. plt.antialias(True)

# Make the grid<br /> x, y, z = meshgrid(arange(-1.5, 1.6, 0.1),<br />                    arange(-1.5, 1.6, 0.1),<br />                    arange(-1.5, 1.6, 0.1))<br /> <br /> # Make the direction data for the arrows<br /> u = x + cos(4*x) + 3         # x-component of vector field<br /> v = sin(4*x) - sin(2*y)      # y-component of vector field<br /> w = -z                       # z-component of vector field<br /> speed = sqrt(u*u + v*v + w*w)<br /> <br /> streamslice(x, y, z, u, v, w, xslice=1.5, yslice=1.5, zslice=-1.5,<br />     color='b', linewidth=1, density=4, interval=5)<br /> xlim(-1.5, 1.5)<br /> ylim(-1.5, 1.5)<br /> zlim(-1.5, 1.5)<br /> <br /> plt.antialias(True)

freekiller

{:eb502:}{:eb502:}

晴天

尊敬的王老师：请问calipso卫星VFM数据的的'Feature_Classification_Flags'参数，可以求平均吗？我求了多年的平均，然后用您的程序作图，感觉有问题，特地请教您这个参数是不是不可以求平均值。它的值是分类的特征代码，不像光学厚度是十进制的数值。