W Pythonie, używając GDAL, wyodrębniłem projekcję rastra jako ciąg WKT w następujący sposób:
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]]'
Za pomocą ciągu WKT mogę utworzyć instancję SpatialReference w następujący sposób:
src = osr.SpatialReference()
src.ImportFromWkt(wkt)
To jest łatwe. Mogę src
dość łatwo wyodrębnić kilka parametrów projekcji, takich jak strefa UTM itp. Ale nie mogę wymyślić, jak wyodrębnić coś takiego jak nazwa projekcji, tj "GDA_1994_Transverse_Mercator"
. Z pewnością jest to możliwe, ale dokumentacja API języka Python może być nieistniejąca dla wszystkich jej zastosowań.
Jak wyodrębnić nazwy rzutu i geograficzny układ współrzędnych?
python
gdal
coordinate-system
python-2.7
hendra
źródło
źródło