python - 如何使用 python 在地理地图上将边缘绘制为 3D 弧？

Question

我想创建一个类似于

在蟒蛇。到目前为止，我能够在 2D 上绘制它。边的密度使图形变得不整齐。在 3-D 中绘制边缘将增强美感并使其更具吸引力。任何帮助将不胜感激。提前致谢

标签： pythonmatplotlib-basemapnetwork-analysis

Answer 1

原则上，你想做的事情是可能的Basemap。您可以按照本教程Basemap在 3D 轴上绘制 a ，然后在其上绘制弧线。但是，如果您想按照示例图所暗示的方式为不同的国家着色，则必须开始处理（例如，参见此处或此处）。无论如何，下面是一个在 3D 图中绘制底图和弧线的示例代码。在 3D 轴上绘制底图的部分直接取自我上面链接的教程。shapefiles

import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from mpl_toolkits.basemap import Basemap
from matplotlib.collections import PolyCollection
import numpy as np

m = Basemap()

fig = plt.figure()
ax = Axes3D(fig)

ax.azim = 270
ax.elev = 50
ax.dist = 8

ax.add_collection3d(m.drawcoastlines(linewidth=0.25))
ax.add_collection3d(m.drawcountries(linewidth=0.35))

polys = []
for polygon in m.landpolygons:
    polys.append(polygon.get_coords())


lc = PolyCollection(polys, edgecolor='black',facecolor='#DDDDDD', closed=False)
ax.add_collection3d(lc)

def plot_arc(ax,p1,p2,height,N=100):
    x = np.linspace(p1[0],p2[0],N)
    y = np.linspace(p1[1],p2[1],N)
    z = (1-np.linspace(-1,1,N)**2)*height
    ax.plot(x,y,z)

plot_arc(ax,(-70,0),(80,50),1)
plot_arc(ax,(80,50),(125,-25),1)
plot_arc(ax,(125,-25),(-70,0),1)

plt.show()

这里是示例代码的结果：

希望这可以帮助。