ESRI中WKT格式的空间投影坐标转换为GDAL中的Proj4格式
原创
©著作权归作者所有:来自51CTO博客作者echohye的原创作品,请联系作者获取转载授权,否则将追究法律责任
首先是自定义的空间投影坐标,可以在在ArcGIS中自定义空间投影坐标,导出得到ESRI标准的.prj后缀的文件,该格式为 Well-Known Text (WKT) 文本。
而在Python的gdal包中,所采用的空间投影坐标的格式为EPSG code或者是proj4格式的字符串,标准格式的EPSG代码可以从https://epsg.io/这个网址查找,而自定义的WKT格式则需要单独转换。
https://spatialreference.org/
解决方法:采用Python包pycrs进行转换,具体代码如下
import pycrs
# wkt text可以直接从ArcGIS导出
wkt_text = """
PROJCS["WGS 84 / World Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Mercator"],PARAMETER["central_meridian",0],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]
"""
srs_proj4 = pycrs.parse.from_esri_wkt(wkt_text).to_proj4()
print(srs_proj4)
pycrs的官方文档为:https://karimbahgat.github.io/PyCRS/,可以查看更多投影坐标格式的转换。
箴言:因为这些东西是非常简单的。不要抱怨自己学不会,那是因为你没有足够用心。