三维飞行路线动图

MeteoInfo

三维模型加载后可以用model函数形成一个Model对象，绘图缩放时大小固定，可以设置位置、旋转角度、缩放比例等。下面的例子加载了一个飞机三维模型以及包括经纬度和高度数组的飞行路线，在循环中设置飞机的位置，根据当前点和下一个的坐标确定旋转角度（只旋转Z轴，飞机保持水平）。运行脚本可以看到飞机沿设定的路线飞行的三维动图。可以用鼠标旋转、缩放三维图形进行更多角度的观察。

数据文件： guiji.csv (428.5 KB, 下载次数: 3)

2025-1-14 08:22 上传

点击文件名下载附件     11805_airplane_v2_L2.obj (6.96 MB, 下载次数: 2)

2025-1-14 08:25 上传

点击文件名下载附件

1. import time
   
2. 
   
3. # Read airplane trajectory data
   
4. fn = 'D:/Temp/ascii/guiji.csv'
   
5. df = DataFrame.read_table(fn, delimiter=',', format='%i%{yyyy/M/dd H:mm}D%5f')
   
6. lon = df['Longitude'].values
   
7. lat = df['Latitude'].values
   
8. alt = df['Altitude'].values
   
9. n = len(lon)
   
10. 
    
11. # Read airplane model data
    
12. mfn = 'D:/Temp/3d/11805_airplane_v2_L2.obj'
    
13. T, x, y, z, normal = plt.load_obj_model(mfn)
    
14. 
    
15. # Plot
    
16. ax = axes3d(axes_zoom=True)
    
17. lighting(mat_specular=1)
    
18. #geoshow('cn_province', edgecolor='b')
    
19. ax.plot(lon, lat, alt, linewidth=2, color='r')
    
20. xlim(111, 113)
    
21. ylim(40, 41)
    
22. zlim(0, 5000)
    
23. xlabel('Longitude')
    
24. ylabel('Latitude')
    
25. zlabel('Altitude (m)')
    
26. antialias()
    
27. ap = model(T, x, y, z, facecolor='c', scale=0.3, location=[lon[0],lat[0],alt[0]])
    
28. n = len(lon)
    
29. for j in range(0, n-1, 20):
    
30.     print(j)
    
31.     ap.setLocation([lon[j], lat[j], alt[j]])
    
32.     a = np.rad2deg(np.atan2(lat[j+1]-lat[j], lon[j+1]-lon[j]))
    
33.     ap.setRotation([0,0,a-90])
    
34.     plt.draw()
    
35.     time.sleep(0.05)

import time<br /> <br /> # Read airplane trajectory data<br /> fn = 'D:/Temp/ascii/guiji.csv'<br /> df = DataFrame.read_table(fn, delimiter=',', format='%i%{yyyy/M/dd H:mm}D%5f')<br /> lon = df['Longitude'].values<br /> lat = df['Latitude'].values<br /> alt = df['Altitude'].values<br /> n = len(lon)<br /> <br /> # Read airplane model data<br /> mfn = 'D:/Temp/3d/11805_airplane_v2_L2.obj'<br /> T, x, y, z, normal = plt.load_obj_model(mfn)<br /> <br /> # Plot<br /> ax = axes3d(axes_zoom=True)<br /> lighting(mat_specular=1)<br /> #geoshow('cn_province', edgecolor='b')<br /> ax.plot(lon, lat, alt, linewidth=2, color='r')<br /> xlim(111, 113)<br /> ylim(40, 41)<br /> zlim(0, 5000)<br /> xlabel('Longitude')<br /> ylabel('Latitude')<br /> zlabel('Altitude (m)')<br /> antialias()<br /> ap = model(T, x, y, z, facecolor='c', scale=0.3, location=[lon[0],lat[0],alt[0]])<br /> n = len(lon)<br /> for j in range(0, n-1, 20):<br />     print(j)<br />     ap.setLocation([lon[j], lat[j], alt[j]])<br />     a = np.rad2deg(np.atan2(lat[j+1]-lat[j], lon[j+1]-lon[j]))<br />     ap.setRotation([0,0,a-90])<br />     plt.draw()<br />     time.sleep(0.05)

20230406_142141.mp4 (8.17 MB, 下载次数: 16)

2023-4-6 16:14 上传

点击文件名下载附件

jimeking

{:eb502:}{:eb502:}{:eb502:}{:eb502:}

yangmt

老师您好，请问能补充下测试数据吗？

MeteoInfo

yangmt 发表于 2025-1-13 23:12
老师您好，请问能补充下测试数据吗？

已经放在帖子一楼了。