python epsg 互转

Python中可以使用pyproj库来进行EPSG（European Petroleum Survey Group）投影坐标系的互转操作。pyproj是一个基于Proj4库的Python封装库，它用于进行地理空间坐标转换和投影操作。

首先，我们需要安装pyproj库。可以通过在终端或命令提示符中运行以下命令来安装：

pip install pyproj

接下来，我们可以使用pyproj库的功能来进行EPSG互转。假设我们要将经纬度坐标（EPSG 4326）转换为UTM投影坐标（EPSG 32613），可以使用以下代码：

import pyproj

# 定义输入坐标系和输出坐标系的EPSG编码
src_crs = pyproj.CRS.from_epsg(4326)
dst_crs = pyproj.CRS.from_epsg(32613)

# 创建坐标转换器
transformer = pyproj.Transformer.from_crs(src_crs, dst_crs)

# 输入经纬度坐标
lon = 120.123456
lat = 30.654321

# 进行转换
x, y = transformer.transform(lon, lat)

# 打印输出UTM投影坐标
print(x, y)

如果我们要将UTM投影坐标（EPSG 32613）转换为经纬度坐标（EPSG 4326），可以使用以下代码：

import pyproj

# 定义输入坐标系和输出坐标系的EPSG编码
src_crs = pyproj.CRS.from_epsg(32613)
dst_crs = pyproj.CRS.from_epsg(4326)

# 创建坐标转换器
transformer = pyproj.Transformer.from_crs(src_crs, dst_crs)

# 输入UTM投影坐标
x = 500000
y = 4000000

# 进行转换
lon, lat = transformer.transform(x, y)

# 打印输出经纬度坐标
print(lon, lat)

通过使用pyproj库的相关功能，我们可以方便地进行EPSG投影坐标系的互转操作。