<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=utf-8">
<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;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
.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">
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US">Hi Pierre,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US">Your horizontal transformation is in geographic (angular) coordinates. So you have to do that before the projection to coordinates in metres.
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US">Regard, Jochem<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span lang="EN-US">From:</span></b><span lang="EN-US"> PROJ <proj-bounces@lists.osgeo.org>
<b>On Behalf Of </b>Pierre Chaponniere<br>
<b>Sent:</b> maandag 10 februari 2020 10:11<br>
<b>To:</b> proj@lists.osgeo.org<br>
<b>Subject:</b> [PROJ] Transformation for local grid<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">Dear all, <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I am struggling to elaborate a pipeline with proj v6.3 that would transform lat long coordinates to a local grid system.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">From what I see, I have all the required parameters, so I did a 3 step pipeline :<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">1/ transform from lat long wgs84 to projected using tmerc.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">2/ run horizontal transformation which is a rotation around an origin - I was thinking of using the the Molobadekas transform for this purpose<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">3/ vertical shift for which I have an origin as well, east and north slope entries in ppm - i still need to figure how I'll pass this one...<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">So my proj line so far looks like this :<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">proj +proj=pipeline +step +proj=tmerc +lat_0=-1.90 +lon_0=115.90 +y_0=538.374 +x_0=5200.27 +ellps=GRS80 +step +proj=molobadekas +convention=coordinate_frame +x=-9.85881089465147 +y=-0.172038726321603 +s=17.60710142 +rz=1485.3205006576768028
+px=5225.30897756133 +py=-830.070127940344<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">But I keep on having errors as 'malformed pipeline' or 'non-angular input coordinates'...<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Even trying to pass the Molobadekas transform on its own, I struggle...<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Do you see anything that I am doing totally wrong here ?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Would be really appreciated to have your insight on this !<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks heaps,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Best regards, <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Pierre C<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</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></body>
</html>