首先是自定義的空間投影坐標贴谎,可以在在ArcGIS中自定義空間投影坐標汞扎,導(dǎo)出得到ESRI標準的.prj后綴的文件,該格式為 Well-Known Text (WKT) 文本擅这。
而在Python的gdal包中澈魄,所采用的空間投影坐標的格式為EPSG code或者是proj4格式的字符串,標準格式的EPSG代碼可以從https://epsg.io/這個網(wǎng)址查找仲翎,而自定義的WKT格式則需要單獨轉(zhuǎn)換痹扇。
解決方法:采用Python包pycrs進行轉(zhuǎn)換,具體代碼如下
import pycrs
# wkt text可以直接從ArcGIS導(dǎo)出
wkt_text = 'PROJCS["China_Lambert",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,' '298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION[' '"Lambert_Conformal_Conic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],' 'PARAMETER["Central_Meridian",105.0],PARAMETER["Standard_Parallel_1",25.0],PARAMETER[' '"Standard_Parallel_2",47.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]'
srs_proj4 = pycrs.parse.from_esri_wkt(wkt_text).to_proj4()
此后溯香,srs_proj4 就可以用在gdal.Warp中的 dstSRS 參數(shù)中了鲫构。
pycrs的官方文檔為:https://karimbahgat.github.io/PyCRS/,可以查看更多投影坐標格式的轉(zhuǎn)換玫坛。