MeteoInfoLab脚本示例：3维图形

MeteoInfo

先创建一个3维坐标系： ax = axes3d()，然后用3维坐标系对象（ax）里的绘图函数绘图。绘图函数有plot, scatter, plot_surface, plot_wireframe, plot_layer, contour, contourf, imshow等。

线条图：

1. z = linspace(0, 1, 100)
   
2. x = z * np.sin(20 * z)
   
3. y = z * np.cos(20 * z)
   
4. 
   
5. #Plot
   
6. ax = axes3d()
   
7. ax.plot(x, y, z, '-b')
   
8. title('3D plot example')

z = linspace(0, 1, 100)<br /> x = z * np.sin(20 * z)<br /> y = z * np.cos(20 * z)<br /> <br /> #Plot<br /> ax = axes3d()<br /> ax.plot(x, y, z, '-b')<br /> title('3D plot example')

散点图：

1. z = linspace(0, 1, 100)
   
2. x = z * np.sin(20 * z)
   
3. y = z * np.cos(20 * z)
   
4. c = x + y
   
5. 
   
6. #Plot
   
7. ax = axes3d()
   
8. points = ax.scatter(x, y, z, c=c)
   
9. colorbar(points,shrink=0.8)
   
10. title('Point 3D plot example')

z = linspace(0, 1, 100)<br /> x = z * np.sin(20 * z)<br /> y = z * np.cos(20 * z)<br /> c = x + y<br /> <br /> #Plot<br /> ax = axes3d()<br /> points = ax.scatter(x, y, z, c=c)<br /> colorbar(points,shrink=0.8)<br /> title('Point 3D plot example')

Wireframe和等值线图：

1. alpha = 0.7
   
2. phi_ext = 2 * pi * 0.5
   
3. N = 25
   
4. x1 = linspace(0, 2*pi, N)
   
5. y1 = linspace(0, 2*pi, N)
   
6. x,y = meshgrid(x1, y1)
   
7. z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)
   
8. z = z.T
   
9. 
   
10. #Plot
    
11. ax = axes3d()
    
12. lines = ax.contour(x1, y1, z, 10, offset=-2)
    
13. ax.plot_wireframe(x, y, z, color='b')
    
14. colorbar(lines)
    
15. title('Wireframe and contour 3D plot example')

alpha = 0.7<br /> phi_ext = 2 * pi * 0.5<br /> N = 25<br /> x1 = linspace(0, 2*pi, N)<br /> y1 = linspace(0, 2*pi, N)<br /> x,y = meshgrid(x1, y1)<br /> z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)<br /> z = z.T<br /> <br /> #Plot<br /> ax = axes3d()<br /> lines = ax.contour(x1, y1, z, 10, offset=-2)<br /> ax.plot_wireframe(x, y, z, color='b')<br /> colorbar(lines)<br /> title('Wireframe and contour 3D plot example')

Wireframe和等值线填色图：

1. alpha = 0.7
   
2. phi_ext = 2 * pi * 0.5
   
3. N = 25
   
4. x1 = linspace(0, 2*pi, N)
   
5. y1 = linspace(0, 2*pi, N)
   
6. x,y = meshgrid(x1, y1)
   
7. z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)
   
8. z = z.T
   
9. 
   
10. #Plot
    
11. ax = axes3d()
    
12. lines = ax.contourf(x1, y1, z, 10, offset=-2)
    
13. ax.plot_wireframe(x, y, z, color='b')
    
14. colorbar(lines)
    
15. title('Wireframe and contourf 3D plot example')

alpha = 0.7<br /> phi_ext = 2 * pi * 0.5<br /> N = 25<br /> x1 = linspace(0, 2*pi, N)<br /> y1 = linspace(0, 2*pi, N)<br /> x,y = meshgrid(x1, y1)<br /> z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)<br /> z = z.T<br /> <br /> #Plot<br /> ax = axes3d()<br /> lines = ax.contourf(x1, y1, z, 10, offset=-2)<br /> ax.plot_wireframe(x, y, z, color='b')<br /> colorbar(lines)<br /> title('Wireframe and contourf 3D plot example')

Surface图：

1. alpha = 0.7
   
2. phi_ext = 2 * pi * 0.5
   
3. x = linspace(0, 2*pi, 100)
   
4. y = linspace(0, 2*pi, 100)
   
5. x,y = meshgrid(x, y)
   
6. z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)
   
7. z = z.T
   
8. 
   
9. #Plot
   
10. ax = axes3d()
    
11. ls = ax.plot_surface(x, y, z, 20, edge=False)
    
12. colorbar(ls,shrink=0.8)
    
13. title('Surface 3D plot example')

alpha = 0.7<br /> phi_ext = 2 * pi * 0.5<br /> x = linspace(0, 2*pi, 100)<br /> y = linspace(0, 2*pi, 100)<br /> x,y = meshgrid(x, y)<br /> z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)<br /> z = z.T<br /> <br /> #Plot<br /> ax = axes3d()<br /> ls = ax.plot_surface(x, y, z, 20, edge=False)<br /> colorbar(ls,shrink=0.8)<br /> title('Surface 3D plot example')

图像：

1. fn = 'D:/Temp/nc/air_clm.nc'
   
2. f = addfile(fn)
   
3. ps = f['aveair'][0,:,:,'120']
   
4. yy = linspace(0, 1., ps.shape[0])
   
5. ps.setdimvalue(0, yy)
   
6. 
   
7. #Map layer
   
8. layer = shaperead('D:/Temp/map/110m-land.shp')
   
9. 
   
10. #Plot
    
11. ax = axes3d(bbox=True)
    
12. ax.plot_layer(layer, color='c', edgecolor='b')
    
13. ls = ax.imshow(ps, 10, offset=120, zdir='x', alpha=0.8)
    
14. colorbar(ls)
    
15. zlim(0, 1)
    
16. xlim(0, 180)
    
17. title('3D imshow x direction example')

fn = 'D:/Temp/nc/air_clm.nc'<br /> f = addfile(fn)<br /> ps = f['aveair'][0,:,:,'120']<br /> yy = linspace(0, 1., ps.shape[0])<br /> ps.setdimvalue(0, yy)<br /> <br /> #Map layer<br /> layer = shaperead('D:/Temp/map/110m-land.shp')<br /> <br /> #Plot<br /> ax = axes3d(bbox=True)<br /> ax.plot_layer(layer, color='c', edgecolor='b')<br /> ls = ax.imshow(ps, 10, offset=120, zdir='x', alpha=0.8)<br /> colorbar(ls)<br /> zlim(0, 1)<br /> xlim(0, 180)<br /> title('3D imshow x direction example')

轨迹图：

1. #Open trajectory data and get trajectory layer
   
2. fn = 'D:/Temp/HYSPLIT/traj_20131211_00'
   
3. f = addfile_hytraj(fn)
   
4. tlayer = f.trajlayer()
   
5. stlayer = f.trajsplayer()
   
6. 
   
7. #Map layer
   
8. layer = shaperead('D:/Temp/map/110m-land.shp')
   
9. 
   
10. #Relief data
    
11. fn = 'D:/Temp/nc/elev.0.25-deg.nc'
    
12. f = addfile(fn)
    
13. elev = f['data'][0,::8,::8]
    
14. elev = elev[:,'0:180']
    
15. elev[elev<0] = 0
    
16. 
    
17. #Plot
    
18. ax = axes3d()
    
19. ls = ax.plot_surface(elev, 20, cmap='MPL_terrain', edge=False)
    
20. ax.plot_layer(layer, edgecolor='g')
    
21. ax.plot_layer(tlayer)
    
22. ax.plot_layer(stlayer, fill=False)
    
23. zlim(0, 10000)
    
24. xlim(0, 180)
    
25. colorbar(ls)
    
26. title('3D trajectory example')

#Open trajectory data and get trajectory layer<br /> fn = 'D:/Temp/HYSPLIT/traj_20131211_00'<br /> f = addfile_hytraj(fn)<br /> tlayer = f.trajlayer()<br /> stlayer = f.trajsplayer()<br /> <br /> #Map layer<br /> layer = shaperead('D:/Temp/map/110m-land.shp')<br /> <br /> #Relief data<br /> fn = 'D:/Temp/nc/elev.0.25-deg.nc'<br /> f = addfile(fn)<br /> elev = f['data'][0,::8,::8]<br /> elev = elev[:,'0:180']<br /> elev[elev<0] = 0<br /> <br /> #Plot<br /> ax = axes3d()<br /> ls = ax.plot_surface(elev, 20, cmap='MPL_terrain', edge=False)<br /> ax.plot_layer(layer, edgecolor='g')<br /> ax.plot_layer(tlayer)<br /> ax.plot_layer(stlayer, fill=False)<br /> zlim(0, 10000)<br /> xlim(0, 180)<br /> colorbar(ls)<br /> title('3D trajectory example')

931404656

好想学习哦

MeteoInfo

rcdxph 发表于 2017-8-2 13:10
老师，实例的 数据没有啊，能否提供下载？

有些数据太大，无法在这里提供下载。你可以加入MeteoInfo QQ群，里面有。

migrant

太酷了。真心不错啊。
厉害

♂雨已~○

绘制3D图需要更到新文件吗？在哪儿下载最新文件？

♂雨已~○

王老师3D图能不能做雷达的呢？要是王老师能将以前那个grib格式雷达反射率演示一个3D效果图那肯定用处越来越大啊！

MeteoInfo

♂雨已~○ 发表于 2017-8-14 19:22
王老师3D图能不能做雷达的呢？要是王老师能将以前那个grib格式雷达反射率演示一个3D效果图那肯定用处越来越 ...

你如果有雷达数据可以试试

♂雨已~○

MeteoInfo 发表于 2017-8-14 23:00
你如果有雷达数据可以试试

我有一个雷达数据，怎么传给你呢，是nc格式的雷达拼图，它里面有高度层次的雷达反射率数据，希望王老师能做一个例子…

MeteoInfo

♂雨已~○ 发表于 2017-8-14 23:25
我有一个雷达数据，怎么传给你呢，是nc格式的雷达拼图，它里面有高度层次的雷达反射率数据，希望王老师能 ...

你可以加入MeteoInfo QQ群： 361327727

♂雨已~○

好的…我马上加进去