通过http://epsg.io轻松访问epsg数据库/
pyepsg的Python项目详细描述
的简单接口
例如,我们可以请求投影坐标系的详细信息 由EPSG代码21781识别,又名。“瑞士CH1903/LV03”:
>>> projection = pyepsg.get(21781) >>> projection <ProjectedCRS: 21781, CH1903 / LV03> >>> projection.domain_of_validity() [5.97, 10.49, 45.83, 47.81] >>> projection.cartesian_cs.axes [<Axis: east / metre>, <Axis: north / metre>]
投影坐标系和大地坐标系也可以转换为 其他形式:
>>> print(projection.as_proj4()[:70] + '...') +proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333 +k_0=1 ... >>> print(projection.as_wkt()[:70] + '...') PROJCS["CH1903 / LV03",GEOGCS["CH1903",DATUM["CH1903",SPHEROID["Bessel...