<div dir="ltr">Dear Even,<div><br></div><div>Thank you, yet again, for a quick resolution. I cherry-picked that change over to my 3.4.0 build and it certainly gets rid of the Proj errors.</div><div><br></div><div>I also found the mistake in the rest of my code that was creating the large values (I had bDstToSrc set to true on the GDALGCPTransform call for some reason).</div><div><br></div><div>All working now! :)</div><div><br></div><div>SE</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Dec 10, 2021 at 3:40 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<p>Those files use a invalid GeogEllipsoidGeoKey=4326 for the WGS 84
ellipsoid. The correct code should be 7030. Sometimes should
report that to the data producer.<br>
</p>
<p>Warning silenced in <a href="https://github.com/OSGeo/gdal/pull/4975" target="_blank">https://github.com/OSGeo/gdal/pull/4975</a><br>
</p>
<div>Le 11/12/2021 à 00:05, Simon Eves a
écrit :<br>
</div>
<blockquote type="cite">
<div dir="ltr">After also finding this thread...
<div><br>
</div>
<div> <a href="https://github.com/OSGeo/gdal/issues/2321" target="_blank">https://github.com/OSGeo/gdal/issues/2321</a><br>
</div>
<div><br>
</div>
<div>...I installed geotiff-bin and ran listgeo on the file.</div>
<div><br>
</div>
<div>It also spits the same message (three times).</div>
<div><br>
</div>
<div>The output is attached.</div>
<div><br>
</div>
<div>S</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Fri, Dec 10, 2021 at 2:56
PM Simon Eves <<a href="mailto:simon.eves@omnisci.com" target="_blank">simon.eves@omnisci.com</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">Further to my earlier post about GCP
Transforms, I now have that code written, and it seems to be
functioning fine, although whether it's generating the right
values is as yet unclear.
<div><br>
</div>
<div>The problem is that I am having trouble creating the
secondary transformation from whatever the output space of
the GCP transform is, to my desired space (4326 in this
case).<br clear="all">
<div><br>
</div>
<div>Whenever I call GetGCPSpatialRef() on the datasource,
I get the following error (twice)</div>
<div><br>
</div>
<div> PROJ: proj_create_from_database: ellipsoid not
found (1)<br>
</div>
<div><br>
</div>
<div>I get the same error (also twice) when I run gdalinfo
on the file, but then with a plausible looking SR dump,
which may of course just be a default:</div>
<div>______________________________________________________</div>
<div><br>
</div>
<div>$ gdalinfo
s1b-iw-grd-vv-20211110t033513-20211110t033538-029520-0385e8-001.tiff<br>
ERROR 1: PROJ: proj_create_from_database: ellipsoid not
found<br>
ERROR 1: PROJ: proj_create_from_database: ellipsoid not
found<br>
Driver: GTiff/GeoTIFF<br>
Files:
s1b-iw-grd-vv-20211110t033513-20211110t033538-029520-0385e8-001.tiff<br>
Size is 25541, 16650<br>
GCP Projection = <br>
GEOGCRS["WGS 84",<br>
DATUM["World Geodetic System 1984",<br>
ELLIPSOID["unnamed",6378137,298.25722356049,<br>
LENGTHUNIT["metre",1]]],<br>
PRIMEM["Greenwich",0,<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
CS[ellipsoidal,2],<br>
AXIS["geodetic latitude (Lat)",north,<br>
ORDER[1],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
AXIS["geodetic longitude (Lon)",east,<br>
ORDER[2],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
ID["EPSG",4326]]<br>
Data axis to CRS axis mapping: 2,1<br>
GCP[ 0]: Id=1, Info=<br>
(0,0) ->
(45.0323375081341,62.6475859336067,149.992539359257)<br>
</div>
<div>[etc.]</div>
<div>
<div>______________________________________________________</div>
<div><br>
</div>
</div>
<div>My constructed OGRCoordinateTransformation for the
second transform is therefore a no-op, but the final
values (which are therefore just the values out of the
GCP transform) are definitely not 4326.</div>
<div>
<div>______________________________________________________</div>
<div><br>
</div>
</div>
<div>omnisql> select * from s1 where rowid<10;<br>
</div>
<div>raster_lon|raster_lat<br>
476788.1263799404|747980.9704571336<br>
460813.8583986053|747120.1413700345<br>
444852.7791957706|746184.4282221086<br>
428904.8887714364|745173.8310133558<br>
412970.1871256025|744088.349743776<br>
397048.6742582691|742927.9844133693<br>
381140.3501694361|741692.7350221358<br>
365245.2148591036|740382.6015700753<br>
349363.2683272713|738997.5840571877<br>
333494.5105739397|737537.6824834733<br>
</div>
<div>
<div>______________________________________________________</div>
<div><br>
</div>
</div>
<div>Those are the first ten pixels of the first scanline
of the (~25K x ~16K) file, btw, so intuitively there is
something worrying me about those values anyway (they
are moving very quickly and non-linearly considering
that's just ten pixels).</div>
<div><br>
</div>
<div>Anyway, some Googling led me to...</div>
<div><br>
</div>
<div> <a href="https://gitlab.orfeo-toolbox.org/s1-tiling/s1tiling/-/issues/46" target="_blank">https://gitlab.orfeo-toolbox.org/s1-tiling/s1tiling/-/issues/46</a><br>
</div>
<div><br>
</div>
<div>...which seems very related, and even mentions a
patch to libGeoTIFF by Even. I checked, and that patch
IS in the (internal, as configured in the build) GeoTIFF
in our GDAL.</div>
<div><br>
</div>
<div>This is with GDAL 3.4.0 and Proj 7.2.1.</div>
<div><br>
</div>
<div>Thanks in advance,</div>
<div><br>
</div>
<div>Simon</div>
<div><br>
</div>
-- <br>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div style="margin:0px;padding:0px 0px 20px;width:2544px;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:medium">
<div>
<div style="font-size:12.8px;margin:8px 0px 0px;padding:0px">
<div>
<div dir="ltr"><span><font color="#888888">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<table style="font-family:Times;width:2544px" cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="vertical-align:top;font-size:0px" align="left">
<table cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="padding:0px 15px 0px 0px;vertical-align:middle" align="left"><font size="2" face="arial,
helvetica,
sans-serif"><a href="http://www.omnisci.com/" target="_blank"><img src="http://www2.omnisci.com/l/298412/2018-09-18/3sqpg/298412/61753/OmniSci_Email_Header2.png"></a><br>
</font></td>
<td style="padding:0px 0px 0px 15px;vertical-align:top" align="left">
<table style="width:215px" cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="vertical-align:top" align="left"><span style="white-space:nowrap;color:rgb(0,0,1)"><span style="color:rgb(14,76,144);font-weight:700"><font size="2" face="arial,
helvetica,
sans-serif">Simon
Eves</font></span></span></td>
</tr>
<tr>
<td style="vertical-align:top" align="left">
<table cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="vertical-align:top" align="left"><span style="white-space:nowrap;color:rgb(0,0,1)"><font size="2" face="arial,
helvetica,
sans-serif">Senior
Graphics
Engineer,
Rendering
Group<br>
100 Montgomery
St (5th
Floor), San
Francisco, CA
94104, USA<br>
</font></span></td>
</tr>
</tbody>
</table>
</td>
</tr>
<tr>
<td style="vertical-align:top" align="left">
<table cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="padding:0px;vertical-align:top" align="left"><br>
</td>
<td style="padding:0px;vertical-align:top" align="left"><br>
</td>
</tr>
<tr>
<td style="padding:0px;vertical-align:top" align="left"><span style="white-space:nowrap;color:rgb(0,0,1)"><font size="2" face="arial,
helvetica,
sans-serif">Email: <a href="mailto:simon.eves@omnisci.com" target="_blank">simon.eves@omnisci.com</a> |
Cell: </font></span></td>
<td style="padding:0px;vertical-align:top" align="left"><span style="white-space:nowrap;color:rgb(0,0,1)"><font size="2" face="arial,
helvetica,
sans-serif">+1
(415) 902-1996</font></span></td>
</tr>
</tbody>
</table>
<br>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</font></span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
-- <br>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div style="margin:0px;padding:0px 0px 20px;width:2544px;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:medium">
<div>
<div style="font-size:12.8px;margin:8px 0px 0px;padding:0px">
<div>
<div dir="ltr"><span><font color="#888888">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<table style="font-family:Times;width:2544px" cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="vertical-align:top;font-size:0px" align="left">
<table cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="padding:0px 15px 0px 0px;vertical-align:middle" align="left"><font size="2" face="arial,
helvetica,
sans-serif"><a href="http://www.omnisci.com/" target="_blank"><img src="http://www2.omnisci.com/l/298412/2018-09-18/3sqpg/298412/61753/OmniSci_Email_Header2.png"></a><br>
</font></td>
<td style="padding:0px 0px 0px 15px;vertical-align:top" align="left">
<table style="width:215px" cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="vertical-align:top" align="left"><span style="white-space:nowrap;color:rgb(0,0,1)"><span style="color:rgb(14,76,144);font-weight:700"><font size="2" face="arial,
helvetica,
sans-serif">Simon
Eves</font></span></span></td>
</tr>
<tr>
<td style="vertical-align:top" align="left">
<table cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="vertical-align:top" align="left"><span style="white-space:nowrap;color:rgb(0,0,1)"><font size="2" face="arial,
helvetica,
sans-serif">Senior
Graphics
Engineer,
Rendering
Group<br>
100 Montgomery
St (5th
Floor), San
Francisco, CA
94104, USA<br>
</font></span></td>
</tr>
</tbody>
</table>
</td>
</tr>
<tr>
<td style="vertical-align:top" align="left">
<table cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td style="padding:0px;vertical-align:top" align="left"><br>
</td>
<td style="padding:0px;vertical-align:top" align="left"><br>
</td>
</tr>
<tr>
<td style="padding:0px;vertical-align:top" align="left"><span style="white-space:nowrap;color:rgb(0,0,1)"><font size="2" face="arial,
helvetica,
sans-serif">Email: <a href="mailto:simon.eves@omnisci.com" target="_blank">simon.eves@omnisci.com</a> |
Cell: </font></span></td>
<td style="padding:0px;vertical-align:top" align="left"><span style="white-space:nowrap;color:rgb(0,0,1)"><font size="2" face="arial,
helvetica,
sans-serif">+1
(415) 902-1996</font></span></td>
</tr>
</tbody>
</table>
<br>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</font></span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset></fieldset>
<pre>_______________________________________________
gdal-dev mailing list
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
</blockquote>
<pre cols="72">--
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div style="margin:0px;padding:0px 0px 20px;width:2544px;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:medium"><div><div style="font-size:12.8px;margin:8px 0px 0px;padding:0px"><div><div dir="ltr"><span><font color="#888888"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><table cellpadding="0" cellspacing="0" border="0" style="font-family:Times;width:2544px"><tbody><tr><td align="left" style="vertical-align:top;font-size:0px"><table cellpadding="0" cellspacing="0" border="0"><tbody><tr><td align="left" style="padding:0px 15px 0px 0px;vertical-align:middle"><font face="arial, helvetica, sans-serif" size="2"><a href="http://www.omnisci.com/" target="_blank"><img src="http://www2.omnisci.com/l/298412/2018-09-18/3sqpg/298412/61753/OmniSci_Email_Header2.png"></a><br></font></td><td align="left" style="padding:0px 0px 0px 15px;vertical-align:top"><table cellpadding="0" cellspacing="0" border="0" style="width:215px"><tbody><tr><td align="left" style="vertical-align:top"><span style="white-space:nowrap;color:rgb(0,0,1)"><span style="color:rgb(14,76,144);font-weight:700"><font face="arial, helvetica, sans-serif" size="2">Simon Eves</font></span></span></td></tr><tr><td align="left" style="vertical-align:top"><table cellpadding="0" cellspacing="0" border="0"><tbody><tr><td align="left" style="vertical-align:top"><span style="white-space:nowrap;color:rgb(0,0,1)"><font face="arial, helvetica, sans-serif" size="2">Senior Graphics Engineer, Rendering Group<br>100 Montgomery St (5th Floor), San Francisco, CA 94104, USA<br></font></span></td></tr></tbody></table></td></tr><tr><td align="left" style="vertical-align:top"><table cellpadding="0" cellspacing="0" border="0"><tbody><tr><td align="left" style="padding:0px;vertical-align:top"><br></td><td align="left" style="padding:0px;vertical-align:top"><br></td></tr><tr><td align="left" style="padding:0px;vertical-align:top"><span style="white-space:nowrap;color:rgb(0,0,1)"><font face="arial, helvetica, sans-serif" size="2">Email: <a href="mailto:simon.eves@omnisci.com" target="_blank">simon.eves@omnisci.com</a> | Cell: </font></span></td><td align="left" style="padding:0px;vertical-align:top"><span style="white-space:nowrap;color:rgb(0,0,1)"><font face="arial, helvetica, sans-serif" size="2">+1 (415) 902-1996</font></span></td></tr></tbody></table><br></td></tr></tbody></table></td></tr></tbody></table></td></tr></tbody></table><br></div></div></div></div></div></div></div></div></div></div></font></span></div><div></div></div></div><div></div></div></div></div></div></div></div></div></div></div></div>