<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Just out of curiosity, I ran your code on my machine and I got exactly the same results as Even.<br>I do not think that Endianity will cause such wrong results. Are you sure you are not changing the values<br>of new_x and new_y before printing them out?<br><br>~Belaid.<br><br>&gt; From: even.rouault@mines-paris.org<br>&gt; To: gdal-dev@lists.osgeo.org<br>&gt; Subject: Re: [gdal-dev] Help understanding GDALApplyGeoTransform and        GDALInvGeoTransform<br>&gt; Date: Fri, 31 Jul 2009 20:20:46 +0200<br>&gt; CC: BCassanova@weather.com<br>&gt; <br>&gt; Bill,<br>&gt; <br>&gt; that's weird. Your code looks correct and I've just compiled it and I get the <br>&gt; following results :<br>&gt; <br>&gt; lat-lon -127, 51<br>&gt; 7.24754e-13, 4.12115e-13<br>&gt; <br>&gt; The new_x, new_y values are almost 0 as expected (the small difference is due <br>&gt; to numerical imprecision when computing the inverse geotransform).<br>&gt; <br>&gt; I can't understand why you get such wrong values.<br>&gt; <br>&gt; Le Friday 31 July 2009 19:56:39 Cassanova, Bill, vous avez écrit&nbsp;:<br>&gt; &gt; double adfGeoTransform[6] = { -127, .011571482, 0, 51, 0, -.008993 };<br>&gt; &gt;<br>&gt; &gt; double adfReverseGeoTransform[6] = { 0, 0, 0, 0, 0, 0 };<br>&gt; &gt;<br>&gt; &gt; GDALInvGeoTransform( adfGeoTransform, adfReverseGeoTransform );<br>&gt; &gt;<br>&gt; &gt; double dfPixel = 0;<br>&gt; &gt; double dfLine = 0;<br>&gt; &gt; double x;<br>&gt; &gt; double y;<br>&gt; &gt;<br>&gt; &gt; GDALApplyGeoTransform( &nbsp; adfGeoTransform, dfPixel, dfLine, &amp;x, &amp;y );<br>&gt; &gt;<br>&gt; &gt; std::cout &lt;&lt; "lat-lon " &lt;&lt; x &lt;&lt; ", " &lt;&lt; y &lt;&lt; std::endl;<br>&gt; &gt;<br>&gt; &gt; double new_x;<br>&gt; &gt; double new_y;<br>&gt; &gt;<br>&gt; &gt; GDALApplyGeoTransform( &nbsp; adfReverseGeoTransform, x, y, &amp;new_x, &amp;new_y );<br>&gt; &gt;<br>&gt; &gt; std::cout &lt;&lt; new_x &lt;&lt; ", " &lt;&lt; new_y &lt;&lt; std::endl;<br>&gt; <br>&gt; <br>&gt; _______________________________________________<br>&gt; gdal-dev mailing list<br>&gt; gdal-dev@lists.osgeo.org<br>&gt; http://lists.osgeo.org/mailman/listinfo/gdal-dev<br><br /><hr />More storage. Better anti-spam and antivirus protection. <a href='http://go.microsoft.com/?linkid=9671352' target='_new'>Hotmail makes it simple.</a></body>
</html>