pyproj是一个非常好用的地图投影和坐标转换库。
pip install pyproj
构造Transformer
from pyproj import Transformer
transformer = Transformer.from_crs(4326, 32651)
transformer = Transformer.from_crs("EPSG:4326", "EPSG:32651")
CRS(Coordinate Reference System)用于指定转换的源坐标系和目标坐标系。4326、32650是坐标系的唯一编号(WKID,Well Known ID),4326对应的就是WGS84地理坐标系统,32650对应WGS_1984_UTM_Zone_51N。
如果对自己需要进行转换的坐标系的WKID不了解,可以从以下两个网站进行查询:
实现坐标变换
utm_x, utm_y = transformer.transform(30.90786255,121.84786891)
print("utm_x:{}, utm_y:{}".format(utm_x, utm_y))
完整的转换代码如下:
from pyproj import Transformer
transformer = Transformer.from_crs(4326, 32651)
transformer = Transformer.from_crs("EPSG:4326", "EPSG:32651")
utm_x, utm_y = transformer.transform(30.90786255,121.84786891)
p