GDAL / Python: Як отримати ім'я системи координат від SpatialReference?
У Python, використовуючи GDAL, я витягнув проекцію растру як рядок WKT таким чином: 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]]' Використовуючи рядок WKT, я можу створити екземпляр SpatialReference наступним чином: src = osr.SpatialReference() src.ImportFromWkt(wkt) Це легко-горох. Я можу досить легко витягти кілька параметрів проекції з srcтаких, як зона UTM тощо. …