<p dir="ltr">On Nov 10, 2013 1:44 AM, "Even Rouault" <<a href="mailto:even.rouault@mines-paris.org">even.rouault@mines-paris.org</a>> wrote:<br>
><br>
> Le dimanche 10 novembre 2013 00:14:42, Trent Piepho a écrit :<br>
> > I've found that unless you call SetGeoTransform() and give an affine<br>
> > transform, most apps, including listgeo and gdalsrsinfo, aren't<br>
> > entirely happy with the georeferencing.  GDAL has a function that<br>
> > computes a transform from GCPs, but it needs to be part of the GDAL<br>
> > code for the dataset driver.  GDAL doesn't automatically do it when a<br>
> > user of GDAL wants an affine transform from a dataset.<br>
><br>
> True. This is (I think) a design decision, as generally GDAL tries to expose<br>
> information as close as possible to what is in the dataset. GCPs alone don't<br>
> provide an unambiguous way of transforming image pixel coordinates into real-<br>
> world coordinates. You have to specify which extrapolation method you want to<br>
> use : polynomial at which order, thin plate spline, ..</p>
<p dir="ltr">From what I've seen many apps just want to translate between pixel,line and geographic coordinates.  The method to do this appears to be to get the SRS from the dataset, use one of it's methods to transform between geographic coordinates and projected coordinates, then get the geotransform coefficients and preform the affine transform with them.</p>

<p dir="ltr">With no geotransform, this doesn't work.  Instead one is supposed to....I don't know.  I don't see any OGR or dataset function that will do the coordinate transform.  I've looked at other apps, and none of them can handle this.   So you end up with the only way to georeference a dataset, and have anything work, is to use a affine transform.<br>

</p>