<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p><br>
</p>
<div class="moz-cite-prefix">Le 26/01/2022 à 20:54, Lesparre, Jochem
a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:AM0P192MB0482F4D50704BFD65CF704EDEF209@AM0P192MB0482.EURP192.PROD.OUTLOOK.COM">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<meta name="Generator" content="Microsoft Word 15 (filtered
medium)">
<style>@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
font-size:10.0pt;
font-family:"Courier New";}p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:36.0pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;}span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:"Calibri",sans-serif;}span.pl-s
{mso-style-name:pl-s;}span.EmailStyle25
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}div.WordSection1
{page:WordSection1;}ol
{margin-bottom:0cm;}ul
{margin-bottom:0cm;}</style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB">Hi Even,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB">Thanks again for your quick reply, but I think
you misunderstood the problem I tried to describe in my last
email. I’ll try to explain our problem a bit more
comprehensive:
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB">We distribute two grid files for the
quasi-geoid of the Netherlands: NLGEO2018 and NAPTRANS2018.
Both expect ETRS89 (EPSG:4937) ellipsoidal height as input
and give the same NAP national levelling system height
(EPSG:5709) as output. However, the latlon coordinates of
NLGEO2018 are in ETRS89 and the latlon coordinates of
NAPTRANS2018 are in the national triangulation system
“Amersfoort” (EPSG:4289). We prefer the use of the
NLGEO2018. We supply NAPTRANS2018 only because this is the
type of VDatum grid file that was needed for 3D
transformation with PROJ4 (see [1] in my previous email).</span></p>
</div>
</blockquote>
ok got it<br>
<blockquote type="cite"
cite="mid:AM0P192MB0482F4D50704BFD65CF704EDEF209@AM0P192MB0482.EURP192.PROD.OUTLOOK.COM">
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB">Creating a Geodetic TIFF Grid file for
NLGEO2018 is straight-forward [A]. However, I can’t create a
Geodetic TIFF Grid file for NAPTRANS2018.</span></p>
</div>
</blockquote>
<br>
<blockquote type="cite"
cite="mid:AM0P192MB0482F4D50704BFD65CF704EDEF209@AM0P192MB0482.EURP192.PROD.OUTLOOK.COM">
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB"> If I use
</span><span
style="font-size:10.0pt;font-family:Consolas;mso-fareast-language:EN-US"
lang="EN-GB">GEOGRAPHIC_TO_VERTICAL</span><span
style="mso-fareast-language:EN-US" lang="EN-GB"> the
</span><span
style="font-size:10.0pt;font-family:Consolas;mso-fareast-language:EN-US"
lang="EN-GB">--interpolation-crs</span><span
style="mso-fareast-language:EN-US" lang="EN-GB"> is ignored
[B].</span></p>
</div>
</blockquote>
<p>Yes, this is indeed ignored. That type is modeled from the EPSG
corresponding transformation method "Geographic3D to
GravityRelatedHeight", and it doesn't accept an interpolation CRS.
Your best option is to use gdal_translate -mo
interpolation_crs_wkt=.... to add it<br>
</p>
<br>
<blockquote type="cite"
cite="mid:AM0P192MB0482F4D50704BFD65CF704EDEF209@AM0P192MB0482.EURP192.PROD.OUTLOOK.COM">
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-GB"> If I use
</span><span
style="font-size:10.0pt;font-family:Consolas;mso-fareast-language:EN-US"
lang="EN-GB">VERTICAL_TO_VERTICAL</span><span
style="mso-fareast-language:EN-US" lang="EN-GB">, I get
</span><span
style="font-size:10.0pt;font-family:Consolas;mso-fareast-language:EN-US"
lang="EN-GB">AssertionError</span><span
style="mso-fareast-language:EN-US" lang="EN-GB"> due to the
use of a 3D instead of a vertical
</span><span
style="font-size:10.0pt;font-family:Consolas;mso-fareast-language:EN-US"
lang="EN-GB">--source-crs</span><span
style="font-size:10.0pt;mso-fareast-language:EN-US"
lang="EN-GB">
</span><span style="mso-fareast-language:EN-US" lang="EN-GB">[C].
</span></p>
</div>
</blockquote>
<p>Using VERTICAL_TO_VERTICAL is definitely not appropriate for your
use case.</p>
<br>
<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>