利用pyproj将经纬度投影为平面坐标以及地理坐标系背景知识

代码示例

import pyproj
#定义投影坐标系
proj = pyproj.Proj(proj='tmerc', lon_0=120, lat_0=15, preserve_units=False)
#将经纬度投影为平面坐标
print(proj(120, 15))#结果为(0.0, 0.0)
#将平面坐标反投影为经纬度坐标
print(proj(0, 0, inverse=True))#结果为(119.99999999999999, 14.999999999999998)

背景知识

地球的三级逼近

地球的自然表面有高山也有洼地，是崎岖不平的，我们要使用数学法则来描述他，就必须找到一个相对规则的数学面。为此，通过大地水准面 (geoid)，参考椭球体（Reference ellipsoid），大地基准面（Geodetic datum）实现对地球表面的三级逼近。

一级逼近：大地水准面

大地水准面是地球表面的第一级逼近。假设当海水处于完全静止的平衡状态时，从海平面延伸到所有大陆下部，而与地球重力方向处处正交的一个连续、闭合的曲面，这就是大地水准面。

二级逼近：参考椭球体

大地水准面表面仍然是“高低起伏”的，但可以近似成一个规则椭球体，其形状接近一个扁率极小的椭圆绕短轴旋转所形成的规则椭球体，这个椭球体