জিডিএল / পাইথন: আমি কীভাবে স্থানীয় স্থানীয় রেফারেন্স থেকে সিস্টেমের নাম সমন্বয় করব?

Question

পাইথনে, জিডিএল ব্যবহার করে, আমি নীচে নীচে 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]]'

ডাব্লুকেটি স্ট্রিং ব্যবহার করে, আমি নিম্নলিখিত হিসাবে একটি স্পেশিয়াল রেফারেন্স উদাহরণ তৈরি করতে পারি:

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

এটি সহজ-উদ্বেগজনক। আমি প্রজেকশনটির বেশ কয়েকটি পরামিতি srcযেমন ইউটিএম জোন ইত্যাদি থেকে সহজেই বের করতে পারি। তবে আমি বুঝতে পারি না যে কীভাবে প্রজেকশনটির নামের মতো কিছু বের করা যায় "GDA_1994_Transverse_Mercator"। এটি অবশ্যই সম্ভব হয়ে উঠেছে, তবে পাইথন এপিআই ডকুমেন্টেশন পাশাপাশি এটির সমস্ত ব্যবহারের অস্তিত্ব থাকতে পারে।

আমি কীভাবে প্রজেকশন এবং ভৌগলিক সমন্বয় সিস্টেমের নামগুলি বের করব?

Answer 1

দেখুন OGR অভিক্ষেপ টিউটোরিয়াল এবং OGRSpatialReference বর্গ। বিশেষত, GetAttrValue পদ্ধতি।

এখানে একটি কাজের উদাহরণ।

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')

আমার রাস্টার জন্য এই মুদ্রণ:

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'