我正在尝试使用 gdal 版本重新投影 shapefile：GDAL 3.0.2，发布于 2019/10/28 和 python 版本：Python 3.7.11。
下面的代码：
# set the desired projection (proj4 and WKT)
proj4 = '+proj=utm +zone=%s +ellps=WGS84 +datum=WGS84 +units=m +no_defs' %(zone)
print(proj4)
srs = osr.SpatialReference()
print(srs)
srs.ImportFromProj4(proj4)
projwkt = srs.ExportToWkt()
print(projwkt)
projectedAOI = aoishp[0:-4] + '_projected.shp'
print(projectedAOI)
# Reproject to desired
os.system('/home/test/anaconda3/envs/forest/bin/ogr2ogr -t_srs %s -overwrite %s %s' %(projwkt, projectedAOI, aoishp))
但我得到这个错误：
PROJCS[unknown,GEOGCS[unknown,DATUM[WGS_1984,SPHEROID[WGS 84,6378137,298.257223563,AUTHORITY[EPSG,7030]],AUTHORITY[EPSG,6326]],PRIMEM[Greenwich,0,AUTHORITY[EPSG,8901]],UNIT[degree,0.0174532925199433,AUTHORITY[EPSG,9122]]],PROJECTION[Transverse_Mercator],PARAMETER[latitude_of_origin,0],PARAMETER[central_meridian,-69],PARAMETER[scale_factor,0.9996],PARAMETER[false_easting,500000],PARAMETER[false_northing,0],UNIT[metre,1,AUTHORITY[EPSG,9001]],AXIS[Easting,EAST],AXIS[Northing,NORTH]]