<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@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;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></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]-->
</head>
<body lang="NL" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-GB">Hi list members,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">I have a transformation pipeline with a horizontal and a vertical grid shift file. For a point outside both grids, I want the horizontal component to apply a zero shift and the vertical component to return a no data value.
I had this working in PROJ 6, but I don’t get it working in PROJ 7 or 8. For points outside the vertical grid area, I get: “TRANSFORMATION ERROR: 52 9 0 (point not within available datum shift grids)”. Is there a way to get horizontal coordinates outside the
vertical grid area without getting a height value? I was thinking that I could create a nodata.tif covering the entire world, but I hope there is a more elegant way to make the PROJ pipeline (see below) for the national CRS of the Netherlands give the results
recommended by the national authorities.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Regards, Jochem<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">cct +proj=pipeline
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=unitconvert +xy_in=deg +z_in=m +xy_out=rad +z_out=m
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=axisswap +order=2,1,3,4
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=vgridshift +grids=nlgeo2018.gtx
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=push +v_3
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=set +v_3=43 +omit_inv
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=cart +ellps=GRS80
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=helmert +x=-565.7346 +y=-50.4058 +z=-465.2895 +rx=-0.395023 +ry=0.330776 +rz=-1.876073 +s=-4.07242 +convention=coordinate_frame +exact
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=cart +inv +ellps=bessel
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=hgridshift +inv +grids=rdcorr2018.gsb,null
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=sterea +lat_0=52.156160556 +lon_0=5.387638889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:Consolas">+step +proj=set +v_3=0 +omit_fwd
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:Consolas">+step +proj=pop +v_3<o:p></o:p></span></p>
</div>
<br>
<br>
<font size="2">Disclaimer:<br>
De inhoud van dit bericht is uitsluitend bestemd voor geadresseerde.<br>
Gebruik van de inhoud van dit bericht door anderen zonder toestemming van het Kadaster<br>
is onrechtmatig. Mocht dit bericht ten onrechte bij u terecht komen, dan verzoeken wij u<br>
dit direct te melden aan de verzender en het bericht te vernietigen.<br>
Aan de inhoud van dit bericht kunnen geen rechten worden ontleend.<br>
<br>
Disclaimer:<br>
The content of this message is meant to be received by the addressee only.<br>
Use of the content of this message by anyone other than the addressee without the consent<br>
of the Kadaster is unlawful. If you have received this message, but are not the addressee,<br>
please contact the sender immediately and destroy the message.<br>
No rights can be derived from the content of this message.<br>
</font>
</body>
</html>