मैं एक रेखापुंज का उपयोग कर रेखांकन करना चाहता हूं 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*)
क्या ऐसा दृष्टिकोण संभव है?