Oh so sorry Jochem<br><br>Thank you for the explanations. I use my own software and I had never realized proj. I am sorry to hear a so usefull program like proj has these bugs. Thank you for informing me.&nbsp;<div id="yMail_cursorElementTracker_1490256558665">Regards Stefanos<br><div id="ymail_android_signature"><a href="https://overview.mail.yahoo.com/mobile/?.src=Android">Στάλθηκε από το Ταχυδρομείο Yahoo σε Android</a></div> <br> <blockquote style="margin: 0 0 20px 0;"> <header style="font-family:Roboto, sans-serif; color:#6D00F6;"> <div>Στις Πέμ, 23 Μαρ, 2017 στις 2:48, ο χρήστηςStefanos Beligiannis</div><div>&lt;smpel@yahoo.com&gt; έγραψε:</div> </header> <div style="padding: 10px 0 0 20px; margin: 10px 0 0 0; border-left: 1px solid #6D00F6;"> <div id="yiv9262072406"><div><div style="color:#000;background-color:#fff;font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px;"><div id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674">Hi Jochem,</div><div id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674"><br clear="none"></div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674">if I understood well you need a method to transform your coordinates from one system to another.</div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674">7 parameters transformation (or Bursa-Wolf) is a geodetic coordinates transformation method, so the first question you must answer is what type of coordinates do you have on both systems and you must describe these systems.</div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674">A second issue that arises is the accuracy. Martin already explained you that 7 param. method gives transformations with errors of even &gt;10 m big. So big errors are expected when we deal with datum (ellipsoids) that have parameters with big differences between them.&nbsp;</div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674">There is a confuse also, if I got the point, on the proj.4 calculation rounding errors that you say are 0.1 mm and the standard deviation of your transformed coordinates which you say are 4.5 cm.</div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674">In turn we can't help you unless you explain clearly : "what you have" and "what you need".</div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674"><br clear="none"></div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674">regards Stefanos Beligiannis &nbsp;</div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674">&nbsp;</div><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329674"><br clear="none"></div><div class="yiv9262072406qtdSeparateBR" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329709"><br clear="none"><br clear="none"></div><div class="yiv9262072406yqt5273599583 yQTDBase" id="yiv9262072406yqt47268"><div class="yiv9262072406yahoo_quoted" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329693" style="display:block;">  <div id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329692" style="font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px;"> <div id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329691" style="font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px;"> <div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329708"> <font id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329710" size="2" face="Arial"> </font><hr id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329711" size="1"> <b id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_330219"><span id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_330218" style="font-weight:bold;">From:</span></b> Martin Desruisseaux &lt;martin.desruisseaux@geomatys.com&gt;<br clear="none"> <b id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_330221"><span id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_330220" style="font-weight:bold;">To:</span></b> proj@lists.maptools.org <br clear="none"> <b><span style="font-weight:bold;">Sent:</span></b> Thursday, March 23, 2017 12:33 AM<br clear="none"> <b id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_330184"><span id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_330183" style="font-weight:bold;">Subject:</span></b> Re: [Proj] +towgs84 approximation error<br clear="none">  </div> <div class="yiv9262072406y_msg_container" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329690"><br clear="none"><div dir="ltr" id="yiv9262072406yui_3_16_0_ym19_1_1490185827564_329689">Hello Jochem<br clear="none"><br clear="none">The 7 parameter coordinate transformations are for datum shifts, which<br clear="none">are themselves approximative operations by nature. They have stochastic<br clear="none">error unrelated to the Proj.4 accuracy, much greater than the sine<br clear="none">approximation. For example the "NAD27 to WGS 84 (4)" transformation over<br clear="none">USA (EPSG:1173) has an accuracy of 10 metres. So in this context I<br clear="none">wonder what would be the purpose of aiming for an "accuracy" of 0.1 mm.<br clear="none">It can not be a positional accuracy (in the geodetic sense) at least,<br clear="none">but maybe the intend was something else?<br clear="none"><br clear="none">&nbsp; &nbsp; Martin<br clear="none"><br clear="none"><div class="yiv9262072406yqt7429681066" id="yiv9262072406yqtfd98030"><br clear="none">Le 22/03/2017 à 22:41, Jochem a écrit :<br clear="none">&gt; Hello, <br clear="none">&gt;<br clear="none">&gt; I joined the mailing list because I noticed that the PROJ.4 library uses<br clear="none">&gt; approximate formulas for the 7 parameter coordinate transformation defined<br clear="none">&gt; by +towgs84 (technical details: The elements of the rotation matrix are the<br clear="none">&gt; rotation angles instead of the product of sines and cosines of the angles).<br clear="none">&gt; For really small rotation angles this gives the same result, but not for<br clear="none">&gt; larger rotations. <br clear="none">&gt;<br clear="none">&gt; For my current application I want to transform a map in QGIS with 1 mm<br clear="none">&gt; accuracy. However, the approximation causes an error of 4 cm. Since PROJ.4<br clear="none">&gt; normally gives results within 0.1 mm accuracy, I did not expect such<br clear="none">&gt; approximation.<br clear="none">&gt;<br clear="none">&gt; What can I do to help to fix this?<br clear="none">&gt;<br clear="none">&gt; Regards, Jochem<br clear="none"><br clear="none">_______________________________________________<br clear="none">Proj mailing list<br clear="none"><a rel="nofollow" shape="rect" ymailto="mailto:Proj@lists.maptools.org" target="_blank" href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br clear="none"><a rel="nofollow" shape="rect" target="_blank" href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</a></div></div><br clear="none"><br clear="none"></div> </div> </div>  </div></div></div></div></div> </div> </blockquote></div>