<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Frank, <br>
<br>
I am using GDAL<br>
Version : 1.9.2<br>
<br>
so maybe that is why I do not see the method SetGCPProjection().<br>
<br>
I think you have convinced me, though, to calculate the 6
coefficients directly,<br>
and then call SetGeoTransform(). Will let you guys know how this
works out<br>
for me.<br>
<br>
Thank you,<br>
Norm<br>
<br>
<br>
On 11/06/2013 01:51 PM, Frank Warmerdam wrote:<br>
</div>
<blockquote
cite="mid:CA+YzLBdy+UVXK=fjZt6TNtYLtA2izcL2-hP20gecqv77De9jZA@mail.gmail.com"
type="cite">
<div dir="ltr">Norman,
<div><br>
</div>
<div>I believe you want to call SetGCPProjection() instead of
SetProjection() when using GCPs instead of an affine
transform. </div>
<div><br>
</div>
<div>Best regards,</div>
<div>Frank</div>
<div><br>
</div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Wed, Nov 6, 2013 at 1:23 PM, Norman
Goldstein <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:normvcr@telus.net" target="_blank">normvcr@telus.net</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">I have
created a geotiff file using the "GTiff" driver, and have
successfully set options for it to be a strips file with no
compression. Also, successfully called the functions<br>
<br>
dataset->SetMetadataItem( "AREA_OR_POINT",<br>
"Point",<br>
nullptr ) )<br>
<br>
and<br>
<br>
dataset->SetGCPs( 3,<br>
gcps,<br>
"Richmond" )<br>
<br>
Here is the listgeo dump:<br>
<br>
############# listgeo dump ##############<br>
Geotiff_Information:<br>
Version: 1<br>
Key_Revision: 1.0<br>
Tagged_Information:<br>
ModelTiepointTag (6,3):<br>
0 249 0<br>
0 0 0<br>
399 249 0<br>
1000 0 0<br>
0 0 0<br>
0 2000 0<br>
End_Of_Tags.<br>
Keyed_Information:<br>
GTRasterTypeGeoKey (Short,1): RasterPixelIsPoint<br>
End_Of_Keys.<br>
End_Of_Geotiff.<br>
<br>
<br>
Corner Coordinates:<br>
... unable to transform points between pixel/line and PCS
space<br>
#####################################<br>
<br>
I also set the reference system using the following code:<br>
<br>
/////////////////////// c++ code
////////////////////////////<br>
OGRSpatialReference oSRS;<br>
oSRS.SetProjCS( "NoWhere" );<br>
oSRS.SetWellKnownGeogCS( "WGS84" );<br>
oSRS.SetEquirectangular( 0.0, // Centre lat<br>
0.0, // Centre lon<br>
0.0, // False Easting<br>
0.0 ); // False Northing<br>
<br>
char* wkt = nullptr;<br>
<br>
if( OGRERR_NONE != oSRS.exportToPrettyWkt( &wkt ) )<br>
{<br>
error...<br>
}<br>
<br>
if( CE_Failure == dataset->SetProjection( wkt ) )<br>
{<br>
error...<br>
}<br>
///////////////////////////////////////////////////////////////<br>
<br>
So, why is listgeo not able to transform points between
pixel/line and PCS space?<br>
<br>
I am happy to upload a full working example if needed.<br>
<br>
Thank you.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a moz-do-not-send="true"
href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a moz-do-not-send="true"
href="http://lists.osgeo.org/mailman/listinfo/gdal-dev"
target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up | Frank Warmerdam, <a
moz-do-not-send="true" href="mailto:warmerdam@pobox.com"
target="_blank">warmerdam@pobox.com</a><br>
light and sound - activate the windows | <a
moz-do-not-send="true" href="http://pobox.com/%7Ewarmerdam"
target="_blank">http://pobox.com/~warmerdam</a><br>
and watch the world go round - Rush | Geospatial Software
Developer<br>
</div>
</blockquote>
<br>
</body>
</html>