<div dir="ltr">After also finding this thread...<div><br></div><div>  <a href="https://github.com/OSGeo/gdal/issues/2321">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">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 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></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 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>