<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:10pt;color:#003366;font-family:Tahoma,Geneva,sans-serif;" dir="ltr">

<p></p>

<div>Dear gdal community,<br>

<br>

I'm currently working with a set of satellite images that are geo-registered with a set of GCPs. To compute a geo-transform from the GCPs, I've been using gdalwarp with:<br>

<br>

gdalwarp -of GTiff path/to/original/file path/to/warped/file<br>

<br>

While "gdalwarp" correctly computes a transformation from the GCPs, at the same it also performs some kind of rectification / north-up image transformation. I'm not sure about the applied operation.<br>

<br>

Presumably, the reason for this is that the original GCPs reflect a transformation with skew factors.<br>

See the following transformations of the original and the warped image. <br>

<br>

>>> import gdal<br>

>>> dataset = gdal.Open("path/to/original/file")<br>

>>> print(gdal.GCPsToGeoTransform(dataset.GetGCPs()))<br>

(-58.57294039342205, -3.991598267026281e-06, 0.0, -34.451175177186684, -3.6492767328575985e-07, 3.1987810722132996e-06)<br>

<br>

>>> import gdal<br>

>>> dataset = gdal.Open("path/to/warped/file")<br>

>>> print(dataset.GetGeoTransform())<br>

(-58.577917916461026, 3.435846284510554e-06, 0.0, -34.44693359348493, 0.0, -3.435846284510554e-06)<br>

<br>

The "rectification" result is very convenient, since it allows to use the images with an existing satellite image segmentation pipeline (which tiles the satellite images using geo-tiles)<br>

However, in order to overlay the segmentations (performed on the "rectified" images) with the original images, I need some information about the transformation applied during "gdalwarp".<br>

<br>

I've searched the documentation for more information - unfortunately without success.

<br>

For example, parameters like "-novshiftgrid" do not affect this transformation.<br>

<br>

<span></span>It would be very helpful, if someone could shed some light on the used "rectification".<br>

What kind of operation is applied? What would be the correct term to search for?<br>

Is it "only" a decomposition of the transform? Or does <span>"gdalwarp</span>" something else too?<br>

Is there a possibility to access this information, while using<span> "gdalwarp</span>"? (Or can this operation also be done by hand?)<br>

Does GDAL contain any functions to invert this "rectification" step - given the corresponding transformation used by "gdalwarp"?<br>

<br>

Best regards,<br>

Sebastian</div>

<br>

<p></p>

<p><br>

</p>

<div id="Signature">

<div style="font-family:Tahoma; font-size:13px">

<div style="font-family:Tahoma; font-size:13px">

<div style="font-family:Tahoma; font-size:13px">

<div style="font-family:Tahoma; font-size:13px"><font face="Tahoma" color="003366"><span lang="de">

<div style="margin:0"><font face="Tahoma" color="003366"><span lang="de"><font size="2"><span style="font-size:11pt"><span lang="en-US"><span lang="de"><font size="2" face="Calibri,sans-serif"><span style="font-size:11pt"><span lang="en-US">--</span></span></font></span></span></span></font></span></font><br>

<font size="2"><span style="font-size:11pt"><span lang="en-US">Dr. Sebastian Bullinger<br>

</span></span></font></div>

<div style="margin:0"><font size="2"><span style="font-size:11pt"><span lang="en-US">Department Object Recognition</span></span></font></div>

<div style="margin:0"><font size="2"><span style="font-size:11pt"><span lang="en-US">Fraunhofer Institute of

</span></span></font></div>

<div style="margin:0"><font size="2"><span style="font-size:11pt"><span lang="en-US">Optronics, Sytem Technologies and Image Exploitation IOSB</span></span></font></div>

<div style="margin:0"><font size="2"><span style="font-size:11pt">Gutleuthausstr. 1, 76275 Ettlingen, Germany</span></font></div>

<div style="margin:0"><font size="2"><span style="font-size:11pt">Phone +49 7243 992-197</span></font></div>

<div style="margin:0"><font size="2"><span style="font-size:11pt"><a href="mailto:sebastian.bullinger@iosb.fraunhofer.de" id="LPNoLP">sebastian.bullinger@iosb.fraunhofer.de</a><br>

</span></font></div>

<div style="margin:0"><font size="2"><span style="font-size:11pt"><a href="https://webmail.iosb.fraunhofer.de/owa/redir.aspx?C=GM65JhP8Gk6mqlokIktWtqcWz_hm-dIIiaYxp7LDE5w39h4r54rwzTEQR1dXSKQtdNkO601Flpk.&URL=http%3a%2f%2fwww.iosb.fraunhofer.de%2f" target="_blank" id="LPNoLP">www.iosb.fraunhofer.de</a></span></font></div>

</span></font></div>

</div>

</div>

</div>

</div>

</div>

</body>

</html>