Bug reading a byn file for Canadian geoid model?

--===============4215115532229865296==
Content-Type: multipart/alternative; boundary="0000000000007059aa06562f2c68"

--0000000000007059aa06562f2c68
Content-Type: text/plain; charset="UTF-8"

Hi

I am trying to read this file, in theory based in ITRF
https://webapp.csrs-scrs.nrcan-rncan.gc.ca/geod/process/download-helper.php?file_id=GSD95
(there is another one based in NAD83xxx)

However in with gdalinfo I see this:
"Warning 1: CRS EPSG:4140 is deprecated. Its non-deprecated replacement
EPSG:4617 will be used instead. To use the original CRS, set the
OSR_USE_NON_DEPRECATED configuration option to NO."

EPSG:4140 https://spatialreference.org/ref/epsg/4140/ is NAD83(CSRS98),
that is not ITRF, as far as I know.

However in byndataset.h there is this line:

constexpr int BYN_DATUM_0 = 4140;  /* ITRF2008 (GRS80 based WGS84) */

https://github.com/OSGeo/gdal/blob/master/frmts/raw/byndataset.h#L168

Are we parsing the header wrongly?

There is some documentation in
https://webapp.csrs-scrs.nrcan-rncan.gc.ca/geod/tools-outils/sample_doc_files/NRCan_Geoid_GridFormat_EN_2023.pdf
I think what I say above applies to the field 16, Datum.

Thanks
Javier

--0000000000007059aa06562f2c68
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr"><div>Hi</div><div><br></div><div>I am trying to read this =
file, in theory based in ITRF</div><div><a href=3D"https://webapp.csrs-scrs=
.nrcan-rncan.gc.ca/geod/process/download-helper.php?file_id=3DGSD95">https:=
//webapp.csrs-scrs.nrcan-rncan.gc.ca/geod/process/download-helper.php?file_=
id=3DGSD95</a></div><div>(there is another one based in NAD83xxx)<br><br></=
div><div>However in with gdalinfo I see this:</div><div>&quot;Warning 1: CR=
S EPSG:4140 is deprecated. Its non-deprecated replacement EPSG:4617 will be=
 used instead. To use the original CRS, set the OSR_USE_NON_DEPRECATED conf=
iguration option to NO.&quot;</div><div><br></div><div>EPSG:4140 <a href=3D=
"https://spatialreference.org/ref/epsg/4140/">https://spatialreference.org/=
ref/epsg/4140/</a> is=C2=A0NAD83(CSRS98), that is not ITRF, as far as I kno=
w.</div><div><br></div><div>However in byndataset.h there is this line:</di=
v><div><br></div><div></div><div>constexpr int BYN_DATUM_0 =3D 4140; =C2=A0=
/* ITRF2008 (GRS80 based WGS84) */</div><div><br></div><div><a href=3D"http=
s://github.com/OSGeo/gdal/blob/master/frmts/raw/byndataset.h#L168">https://=
github.com/OSGeo/gdal/blob/master/frmts/raw/byndataset.h#L168</a></div><div=
><br></div><div>Are we parsing the header wrongly?</div><div><br></div><div=
>There is some documentation in=C2=A0<a href=3D"https://webapp.csrs-scrs.nr=
can-rncan.gc.ca/geod/tools-outils/sample_doc_files/NRCan_Geoid_GridFormat_E=
N_2023.pdf">https://webapp.csrs-scrs.nrcan-rncan.gc.ca/geod/tools-outils/sa=
mple_doc_files/NRCan_Geoid_GridFormat_EN_2023.pdf</a></div><div>I think wha=
t I say above applies to the field 16, Datum.</div><div><br></div><div>Than=
ks</div><div>Javier</div></div>

--0000000000007059aa06562f2c68--

--===============4215115532229865296==
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev

--===============4215115532229865296==--