আমি একজন রাস্টার ব্যবহার করে python
এবং জিওরিফারেন্স করতে চাই GDAL
। আমার বর্তমান পদ্ধতির কল gdal_translate
এবং gdalwarp
ব্যবহার os.system
এবং স্থল নিয়ন্ত্রণ পয়েন্টগুলির একটি কুরুচিপূর্ণ তালিকা। আমি সত্যিই অভ্যন্তরীণভাবে এটি করার একটি উপায় চাই python
।
এটি আমি ব্যবহার করছি বর্তমান প্রক্রিয়া:
import os
os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 -7.9164e+06 5.21094e+06 -gcp 865.827 107.699 -7.91651e+06 5.21104e+06 "inraster.tif" "outraster1.tif"')
os.system('gdalwarp -r bilinear -tps -co COMPRESS=NONE "outraster2.tif" "outraster3.tif"')
পূর্ববর্তী নেই প্রশ্ন এবং উত্তর 2012 থেকে যা যে gdal_translate
আমদানির পর অ্যাক্সেস করা যেতে পারে gdal
। অপ্রচলিত কিনা তা আমি নিশ্চিত নই, বা এটি ভুল কিনা তবে আমি যখন দৌড়ে from osgeo import gdal
যাই তখন আমি gdal.gdal_translate
কোনও বিকল্প হিসাবে দেখি না ।
আমি জানি না এটি বিদ্যমান কিনা তবে আমি পছন্দ করতাম যদি আমি অজগর উপায়ে রাস্টারদের অনুবাদ ও পুনরায় প্রচার করতে পারি। উদাহরণ স্বরূপ:
# translate
gcp_points = [(1251.92, 414.538), (-7.9164e+06, 5.21094e+06)]
gdal.gdal_translate(in_raster, gcp_points, out_raster1)
# warp
gdal.gdalwarp(out_raster1, out_raster2, 'bilinear', args*)
এ জাতীয় পন্থা কি সম্ভব?