<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<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 class="moz-txt-link-freetext" href="https://github.com/OSGeo/gdal/pull/4975">https://github.com/OSGeo/gdal/pull/4975</a><br>
</p>
<div class="moz-cite-prefix">Le 11/12/2021 à 00:05, Simon Eves a
écrit :<br>
</div>
<blockquote type="cite"
cite="mid:CAJf0KTSAZZBQc1afLghw0FzcsRA4KjNrE=W7PtiCjTqKhXLbLw@mail.gmail.com">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
<div dir="ltr">After also finding this thread...
<div><br>
</div>
<div> <a href="https://github.com/OSGeo/gdal/issues/2321"
moz-do-not-send="true">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"
moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true"><img
src="http://www2.omnisci.com/l/298412/2018-09-18/3sqpg/298412/61753/OmniSci_Email_Header2.png"
moz-do-not-send="true"></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"
moz-do-not-send="true">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" class="gmail_signature">
<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" moz-do-not-send="true"><img
src="http://www2.omnisci.com/l/298412/2018-09-18/3sqpg/298412/61753/OmniSci_Email_Header2.png"
moz-do-not-send="true"></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"
moz-do-not-send="true">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 class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>