如何将绘图坐标转换为地理坐标？附Basemap经纬度绘图代码

嘿，你已经搞定了Basemap的地图可视化，现在要解决绘图坐标转地理坐标的问题对吧？其实Basemap本身就提供了非常直接的方法来实现这个转换，下面给你详细讲两种常用方式：

1. 直接使用Basemap的逆转换功能 #

Basemap的核心__call__方法支持通过inverse=True参数，把投影后的绘图坐标直接转换成经纬度（地理坐标）。

单个点转换 #

假设你有一个绘图坐标点（比如从交互点击获取的x_plot, y_plot），可以这样写：

# 假设x_plot和y_plot是你的绘图坐标值
lon, lat = m(x_plot, y_plot, inverse=True)
print(f"对应的地理坐标：经度{lon:.4f}，纬度{lat:.4f}")

批量点转换 #

如果是一组绘图坐标数组，同样可以批量处理：

import numpy as np
# 示例批量绘图坐标
x_plots = np.array([120000, 240000, 360000])
y_plots = np.array([50000, 150000, 250000])
# 批量转换为地理坐标
lons, lats = m(x_plots, y_plots, inverse=True)

2. 使用invert()方法 #

Basemap还提供了专门的invert()方法，功能和上面的方式完全一致，写法更直观：

# 单个点转换
lon, lat = m.invert(x_plot, y_plot)

# 批量点转换
lons, lats = m.invert(x_plots, y_plots)

注意：处理像素/归一化坐标 #

如果你获取的是像素坐标（比如从鼠标点击事件中得到的屏幕像素位置）或者归一化轴坐标（范围0-1），需要先把它们转换成Basemap的投影坐标，再进行地理坐标转换：

import matplotlib.pyplot as plt

# 假设你拿到了屏幕像素坐标(px, py)
fig, ax = plt.gcf(), plt.gca()
# 将像素坐标转换为Basemap的投影坐标
x_plot, y_plot = ax.transData.inverted().transform((px, py))
# 再转成地理坐标
lon, lat = m(x_plot, y_plot, inverse=True)

结合你现有的代码，只需要在绘图完成后加入上述转换逻辑就能实现需求啦~

内容的提问来源于stack exchange，提问作者GeoMonkey