matplotlib 底图绘制经度/经度坐标不正确



我正在尝试学习matplotlib和映射函数底图,并尝试在地图上绘制一个简单的经度/经度坐标列表。但是,底图坐标转换为简单的 x/y 坐标正在绘制地图之外的点,我无法弄清楚原因。我的代码如下:

locationData = parseFile(locationFile)
fig = plt.figure(figsize=(8,8))
m = Basemap(projection='merc', resolution='h',
llcrnrlat=49, llcrnrlon=-13.5,
urcrnrlat=59.5, urcrnrlon=4)
m.drawcoastlines()
m.drawmapboundary(fill_color='aqua')
m.fillcontinents(color='green', lake_color='aqua')
index=0
for entry in locationData:
lat = entry["latitudeE7"]/1E7
long = entry["longitudeE7"]/1E7
x,y = m(lat, long)
print("Plotting {} {} to x{} y{}".format(lat,long,x,y))
plt.plot(x,y, 'ok', markersize=20)
# break at 30 points for testing
if index>30: break
index+=1
plt.show()

这是我从 print 语句中获得的输出:

Plotting 52.------- 1.------- to x79364--.------- y-31476--.-------
Plotting 52.------- 1.------- to x79368--.------- y-31475--.-------
Plotting 52.------- 1.------- to x79362--.------- y-31471--.-------
Plotting 52.------- 1.------- to x79361--.------- y-31472--.-------
Plotting 52.------- 1.------- to x79360--.------- y-31475--.-------
Plotting 52.------- 1.------- to x79365--.------- y-31476--.-------
Plotting 52.------- 1.------- to x79361--.------- y-31476--.-------
...

出于显而易见的原因,我已经审查了确切的值,但您可以看到指向英国的 52°N 1°E 的预期值被大大偏离了图表,这是英国的图表。将地图扩展到整个地球,可以发现它正在绘制马达加斯加北海岸的点。

我从Google位置历史记录下载中获取坐标,然后将其除以10 ^7,因为它们存储为整数。print 语句显示正在正确解析这些内容。

我是 matplotlib 和底图的新手,在 Windows 10 上运行 Python 3.6。有什么帮助吗?

编辑 - 发布我的旧代码我的错误

来自底图文档

使用参数 lon 调用底图类实例时,lat 会将 lon/lat(以度为单位)转换为 x/y 地图投影坐标(以米为单位)。如果可选关键字 inverse 为 True(默认值为 False),则执行从 x/y 到 lon/lat 的逆变换。

而不是x,y = m(lat, long)因此你需要

x,y = m(long, lat)

最新更新