<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:12.0pt;
font-family:"Calibri",sans-serif;
mso-ligatures:standardcontextual;
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;
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="EN-GB" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">I have coordinates of a triangulation pilar in EPSG:10645 (compound projected CRS: DPnet Saba + Saba height) in an input file:<o:p></o:p></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">echo 4740.3800 1890.1200 870.3900 > input.txt<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The transformation to EPSG:10638 (3D geographic CRS: BES2020 Saba) with PROJ 9.6.0 gives correct results:<o:p></o:p></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">cs2cs epsg:10645 epsg:10638 --only-best -f %.9f input.txt<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">17.635528781 -63.237251291 827.6282</span>[...]<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I saved the output coordinates to file:<o:p></o:p></p>
<p class="MsoNormal"><span style="font-family:Consolas"> echo 17.635528781 -63.237251291 827.6282 > output.txt<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The used transformation is:<o:p></o:p></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">projinfo -s epsg:10645 -t epsg:10638 -o proj<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">Candidate operations found: 4<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">-------------------------------------<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">Operation No. 1:<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas"><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">unknown id, Inverse of Saba Transverse Mercator 2020 + Inverse of Saba to Saba height (1) + Saba to BES2020 Saba (1),
<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt">[…]<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">However, the location of input.txt and output.txt is visually 0.1 m different in QGIS (version 3.42.2):<o:p></o:p></p>
<p class="MsoNormal">- open input.txt in QGIS (field_1, field_2, field_3 as EPSG:10645)<o:p></o:p></p>
<p class="MsoNormal">- open output.txt in QGIS (field_2, field_1, field_3 as EPSG:10638), where I selected this transformation (the inverse of the transformation used by cs2cs):<o:p></o:p></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">Inverse of Saba to BES2020 Saba (1) + Saba to Saba height (1) + Saba Transverse Mercator 2020<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The wrong visualisation of the location of output.txt can be exported in QGIS as CSV geometry AS_XYZ:<o:p></o:p></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">4740.3647</span>[...]<span style="font-family:Consolas">,1890.2050</span>[...]<span style="font-family:Consolas">,827.6282</span>[...]<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">This is exactly the xy that I get when I transform output.txt back to EPSG:10645 with zero height:<o:p></o:p></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">echo 17.635528781 -63.237251291 0 | cs2cs epsg:10638 epsg:10645 --only-best -f %.4f
<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas">4740.3647 1890.2052 42.4286<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-family:Consolas"><o:p> </o:p></span></p>
<p class="MsoNormal">Is this the intended behaviour for QGIS? I don’t like it…<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Jochem<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<br>
<br>
<font size="2">Disclaimer:<br>
De inhoud van deze e-mail is vertrouwelijk en uitsluitend bestemd voor de geadresseerde(n).<br>
Gebruik, openbaarmaking, vermenigvuldiging, verspreiding en/of verstrekking van deze informatie aan derden is niet toegestaan.<br>
Op al onze producten en diensten zijn onze algemene leveringsvoorwaarden van toepassing<br>
[https://www.kadaster.nl/algemene-leveringsvoorwaarden].<br>
<br>
Disclaimer:<br>
This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed.<br>
If you are not the intended recipient, you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited.<br>
Our general terms and conditions of delivery apply to all our products and services<br>
[https://www.kadaster.com/general-terms-and-conditions]. <br>
</font>
</body>
</html>