Python坐标系转换

1. 相关名词

• 地心地固坐标系（Earth-Centered, Earth-Fixed，ECEF），简称地心坐标系。
• 地理坐标系统（Geographic Coordinate System，GCS），坐标系是地心坐标系，用经纬度表示球面上的点。
• 世界大地测量系统（World Geodetic System, WGS），比如WGS84，是一种地理坐标系统，用于全球定位系统（GPS）。
• 投影坐标系统（Projection Coordinate System，PCS），在二维平面上用米表示位置。
• 通用横轴墨卡托投影（Universal Transverse Mercator，UTM），是一种投影方法。

坐标系及其对应 WKID（可以理解为唯一标识符）：

• WGS84：4326
• Web_Mercator_Auxiliary_Sphere（Web墨卡托投影坐标系）：3857

关于Web墨卡托投影：是Web地图应用的事实标准。

1. Web墨卡托投影维基百科
2. ArcMap关于墨卡托投影介绍

Python 提供了第三方库 pyproj 可以方便地进行坐标转换。下面主要介绍其中的一个函数：

pip3 install pyproj

transform(p1, p2, x, y, z=None, radians=False)
# 示例
x2, y2, z2 = transform(p1, p2, x1, y1, z1, radians=False)

函数标识在 p1 坐标系和 p2 坐标系之间进行坐标转换，x1, y1, z1 是由 p1 坐标系定义的坐标，z 为高度，单位是米。x2, y2, z2是由 p2 坐标系定义的坐标，它们是经过转化后返回的，默认 z1 = none。Randians 参数标识是否用弧度返回值。

参考文章：

1. Python各坐标系之间的转换
2. python–坐标系转换算法
3. Python地理坐标系和投影坐标系转换
4. 你必须知道的地理坐标系和投影坐标系
5. 利用Pyproj进行地理投影坐标系转换