Python中使用墨卡托投影方法把经纬度转换成x和y的坐标

可以使用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