在Python中使用fiona
和osgeo
处理坐标系时,似乎有很多方法可以通过导入/导出不同的crs格式来定义坐标系,例如:
菲奥娜:
from fiona.crs import from_epsg,from_string,to_string
# Import crs from different formats:
wgs = from_epsg(4326)
wgs = from_string("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs ")
# Export crs as proj4 string
wgs_proj4_string = to_string(wgs)
OSGEO:
^{pr2}$然而,我注意到,这两个库都允许通过EPSG代码轻松定义crs,但是它们都缺少一个反函数(将crs导出为ESPG代码)。在
我最接近的EPSG代码是:
srs.AutoIdentifyEPSG()
epsg = srs.GetAuthorityCode(None)
但是它似乎不那么可靠,而且其他提议的solutions似乎包含了大量的调整或者至少是{a2}依赖关系。在
问题:
有人能给我演示一个简单、直接的方法,用python将CRS导出为EPSG代码吗?像Fiona
中的to_epsg()
或{ExportToEPSG()
?
有人能解释一下在整个互联网上EPSG出口可能性如此短缺的理论背景,特别是与EPSG代码导入的容易程度相比。EPSG代码的全部意义不在于使坐标系易于识别和使用,而不适合没有高级地理空间专业知识的人使用?它不应该像CRS的身份证一样容易检索吗?
目前没有回答
相关问题 更多 >
编程相关推荐