क्षमा याचना यदि निम्न प्रश्न कुछ हद तक मूर्खतापूर्ण है, लेकिन मैं इस पूरी जीआईएस बात के लिए केवल बहुत ही नया हूं।
मैं अजगर में गदल का उपयोग करके कुछ अनुमानित भू-आकृति छवियों को WGS84 में बदलने की कोशिश कर रहा हूं। मैंने एक पोस्ट पाया है जो अनुमानित GeoTiffs के भीतर निम्नलिखित के समान कुछ का उपयोग करके बिंदुओं को बदलने की प्रक्रिया को रेखांकित करता है:
from osgeo import osr, gdal
# get the existing coordinate system
ds = gdal.Open('path/to/file')
old_cs= osr.SpatialReference()
old_cs.ImportFromWkt(ds.GetProjectionRef())
# create the new coordinate system
wgs84_wkt = """
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]]"""
new_cs = osr.SpatialReference()
new_cs .ImportFromWkt(wgs84_wkt)
# create a transform object to convert between coordinate systems
transform = osr.CoordinateTransformation(old_cs,new_cs)
#get the point to transform, pixel (0,0) in this case
width = ds.RasterXSize
height = ds.RasterYSize
gt = ds.GetGeoTransform()
minx = gt[0]
miny = gt[3] + width*gt[4] + height*gt[5]
#get the coordinates in lat long
latlong = transform.TransformPoint(x,y) मेरा सवाल है, अगर मैं इन बिंदुओं को परिवर्तित करना चाहता हूं और एक नई WGS84 जियोफाई फाइल बनाना चाहता हूं, तो क्या इसके बारे में जाने का सबसे अच्छा तरीका है? क्या कोई फ़ंक्शन मौजूद है जो 1 चरण में कार्य जैसे करेगा?
धन्यवाद!