Python 经纬度坐标转换墨卡托投影 - 使用 pyproj 库

可以使用 pyproj 库将经纬度坐标转换为墨卡托投影坐标。

首先安装 pyproj 库：

pip install pyproj

然后可以使用以下代码进行转换：

import pyproj

# 定义原始坐标系和目标坐标系
wgs84 = pyproj.Proj('+proj=longlat +datum=WGS84')
mercator = pyproj.Proj('+proj=merc')

# 输入经纬度坐标，输出墨卡托投影坐标
x, y = pyproj.transform(wgs84, mercator, longitude, latitude)

其中 longitude 和 latitude 是经纬度值，x 和 y 是墨卡托投影坐标值。

需要注意的是，墨卡托投影坐标系的单位是米，因此得到的 x 和 y 值需要进行单位转换，如：

x /= 1000  # 将单位转换为千米
y /= 1000

完整代码示例：

import pyproj

# 定义原始坐标系和目标坐标系
wgs84 = pyproj.Proj('+proj=longlat +datum=WGS84')
mercator = pyproj.Proj('+proj=merc')

# 输入经纬度坐标，输出墨卡托投影坐标
longitude = 116.3975
latitude = 39.9086
x, y = pyproj.transform(wgs84, mercator, longitude, latitude)

# 将单位转换为千米
x /= 1000
y /= 1000

print(x, y)

输出结果：

12954.312922686071 4861.075151702097