在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区域等中提取投影的几个参数。但是我不知道如何提取类似投影名称的内容,即"GDA_1994_Transverse_Mercator"
。这肯定是有可能的,但是就其所有用途而言,Python API文档也可能不存在。
如何提取投影和地理坐标系的名称?