qiu_xingye的博客

代码示例import pyprojgeod = pyproj.Geod(ellps="WGS84")#由经纬度，方位，距离计算另一个点的经纬度lon, lat, back_az = geod.fwd(120, 0, 180, 100000)#输出：(120.0, -0.9043687229127633, 0.0)#由两点经纬度计算距离和方位forward_az, back_az, dist = geod.inv(120, 0, lon, lat)#输出：(180.0, 0.0, 99999.

pyproj.Proj也支持直接通过epsg代码设置投影坐标系。EPSP的英文全称是European Petroleum Survey Group，中文名称为欧洲石油调查组织。它负责维护并发布坐标参照系统的数据集参数，以及坐标转换描述，该数据集被广泛接受并使用，通过一个Web发布平台进行分发，同时提供了微软Acess数据库的存储文件，通过SQL脚本文件，mySQL， Oracle和PostgreSQL等数据库也可使用。目前已有的参考椭球体，大地坐标系，投影坐标系等不同组合都对应着不同的epsg ID号，

代码示例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)基本知识投影地