GDAL / Python: Bagaimana cara saya mendapatkan nama sistem koordinat dari SpatialReference?

Dalam Python, menggunakan GDAL, saya telah mengekstrak proyeksi raster sebagai string WKT sebagai berikut:

wkt = dataset.GetProjection()
# wkt is 'PROJCS["GDA_1994_Transverse_Mercator",GEOGCS["GDA_1994",DATUM["GDA_1994",SPHEROID["GRS_1980",6378137,298.2572221010002],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",117],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["Meter",1]]'

Menggunakan string WKT, saya bisa membuat contoh SpatialReference sebagai berikut:

src = osr.SpatialReference()
src.ImportFromWkt(wkt)

Ini mudah sekali. Saya dapat mengekstrak beberapa parameter proyeksi dari srcseperti zona UTM, dll., Cukup mudah. Tapi saya tidak tahu bagaimana cara mengekstrak sesuatu seperti nama proyeksi, yaitu "GDA_1994_Transverse_Mercator". Ini pasti harus dimungkinkan, tetapi dokumentasi API Python mungkin juga tidak ada untuk semua penggunaannya.

Bagaimana cara mengekstrak nama proyeksi dan sistem koordinat geografis?

Lihat tutorial Proyeksi OGR dan kelas OGRSpatialReference . Secara khusus, metode GetAttrValue .

Inilah contoh yang berhasil.

from osgeo import gdal,osr
ds=gdal.Open(r'SOMERASTER.TIF')
prj=ds.GetProjection()
print prj

srs=osr.SpatialReference(wkt=prj)
if srs.IsProjected:
    print srs.GetAttrValue('projcs')
print srs.GetAttrValue('geogcs')

Untuk raster saya, ini dicetak:

PROJCS["WGS 84 / UTM zone 55N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",147],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","32655"]]
'WGS 84 / UTM zone 55N'
'WGS 84'