地理坐标系转换 python

在 Python 中，可以使用多种库来进行地理坐标系的转换，其中最常用的是 `pyproj` 和 `geopy`。下面是使用这两个库进行地理坐标系转换的示例代码：

使用 `pyproj` 库进行坐标系转换的示例代码：

from pyproj import Transformer

# 定义输入和输出坐标系
input_crs = 'EPSG:4326'  # WGS84 坐标系
output_crs = 'EPSG:3857'  # Web Mercator 坐标系

# 创建转换器
transformer = Transformer.from_crs(input_crs, output_crs, always_xy=True)

# 输入经纬度坐标
lon, lat = 116.397458, 39.909716

# 坐标系转换
x, y = transformer.transform(lon, lat)

# 输出结果
print(f'转换后的坐标：{x}, {y}')

使用 `geopy` 库进行坐标系转换的示例代码：

from geopy.point import Point

# 定义输入和输出坐标系
input_crs = 'EPSG:4326'  # WGS84 坐标系
output_crs = 'EPSG:3857'  # Web Mercator 坐标系

# 创建输入点对象
input_point = Point(39.909716, 116.397458, 0)  # 输入经纬度坐标（纬度在前，经度在后）

# 坐标系转换
output_point = input_point.to(output_crs)

# 输出结果
print(f'转换后的坐标：{output_point.latitude}, {output_point.longitude}')

以上代码中，`EPSG:4326` 表示 WGS84 坐标系，也就是经纬度坐标系，`EPSG:3857` 表示 Web Mercator 坐标系，也就是墨卡托投影坐标系。你可以根据你的实际需求选择不同的输入和输出坐标系。