嗨,我正在尝试使用 Basemap(matplot lib 的一部分)中的大圆函数将 ip 地址图映射到世界 map 上,但是每次我连接跨越太平洋的两个点时(即美国西海岸的某个地方和澳大利亚)有一条水平线出现在我的情节中。
我知道这是因为这个问题:
Note Cannot handle situations in which the great circle intersects the edge of the map projection domain, and then re-enters the domain. (http://matplotlib.org/basemap/api/basemap_api.html#mpl_toolkits.basemap.Basemap.drawgreatcircle)
只是想知道是否有人知道如何修复它或知道另一个没有这个问题的软件包。谢谢!
最佳答案
要在特定情况下在 Basemap 中解决此问题,您可以检查结果路径中 x 顶点的差异。从那里,我们可以将坏路径分成两部分。我已经将我自己的示例与 basemap 放在一起:
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
m = Basemap(projection='cyl', lon_0=0, resolution='c')
m.fillcontinents(color='coral',lake_color='aqua')
m.drawmapboundary(fill_color='aqua')
places = {'Mexico City': (19.05, -99.366667),
'London': (51.507222, -0.1275),
'Sydney': (-33.859972, 151.211111),
'Cape Town': (-33.925278, 18.423889),
'Delhi': (28.61, 77.23),
}
network = [
('London', 'Delhi'),
('Mexico City', 'Sydney'),
]
for source, target in network:
lat1, lon1 = places[source]
lat2, lon2 = places[target]
line, = m.drawgreatcircle(lon1, lat1, lon2, lat2, lw=3)
p = line.get_path()
# find the index which crosses the dateline (the delta is large)
cut_point = np.where(np.abs(np.diff(p.vertices[:, 0])) > 200)[0]
if cut_point:
cut_point = cut_point[0]
# create new vertices with a nan inbetween and set those as the path's vertices
new_verts = np.concatenate(
[p.vertices[:cut_point, :],
[[np.nan, np.nan]],
p.vertices[cut_point+1:, :]]
)
p.codes = None
p.vertices = new_verts
plt.show()
结果:
这不是一个通用的解决方案,只有当你有纬度和经度时才有效。解决一般问题是 很多, 更难,是 cartopy (http://scitools.org.uk/cartopy/docs/latest/) 的重点。在接下来的几天里,我将发布一个在 cartopy 中执行完全相同绘图的示例。
HTH
关于Python basemap drawgreatcircle 函数,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/13888566/