নিম্নলিখিত প্রশ্নটি কিছুটা বোকা হলে দুঃখিত, তবে আমি এই পুরো জিআইএস জিনিসটিতে কেবলমাত্র খুব নতুন।
আমি পাইথনের জিডিএল ব্যবহার করে কিছু প্রত্যাশিত জিওটিফ চিত্রগুলি WGS84 এ রূপান্তর করার চেষ্টা করছি। আমি একটি পোস্ট পেয়েছি যা নীচের মতো কিছু ব্যবহার করে অনুমিত জিওটিফগুলির মধ্যে পয়েন্টগুলি রূপান্তর করার প্রক্রিয়াটির রূপরেখা দেয়:
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) আমার প্রশ্ন হ'ল আমি যদি এই পয়েন্টগুলি রূপান্তর করতে এবং একটি নতুন ডাব্লু জিএস 84 জিওটিফ ফাইল তৈরি করতে চাই তবে এটি কি সবচেয়ে ভাল উপায়? এমন কোনও ফাংশন আছে যা 1 টি পদক্ষেপে যেমন কাজ করবে?
ধন্যবাদ!