Question

我正在尝试学习matplotlib和映射函数Basemap，并试图在地图上绘制一个简单的纬度/经度坐标列表。但是，将Basemap坐标转换为简单的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和basemap的新手，我在Windows 10上运行Python 3.6。有什么帮助吗？

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

Answer 1

来自basemap documentation

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

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

x,y = m(long, lat)

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

1 个答案: