ESRI中WKT格式的空间投影坐标转换为GDAL中的Proj4格式

首先是自定义的空间投影坐标，可以在在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/，可以查看更多投影坐标格式的转换。