ใน 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 อาจไม่มีอยู่เลยสำหรับการใช้งานทั้งหมด
ฉันจะแยกชื่อของเส้นโครงและระบบพิกัดทางภูมิศาสตร์ได้อย่างไร