<div dir="ltr"><div>Hi Jochem<br></div><div><br></div><div>After looking at the code, this error happens when the coordinates are out of range. The default value for that range is 500000</div><div>You can change it with <a href="https://proj.org/en/stable/operations/transformations/horner.html#cmdoption-arg-range">https://proj.org/en/stable/operations/transformations/horner.html#cmdoption-arg-range</a></div><div><br></div><div>I hope it helps<br></div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Wed, 11 Dec 2024 at 11:54, Lesparre, Jochem via PROJ <<a href="mailto:proj@lists.osgeo.org">proj@lists.osgeo.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-1915382970844817410">
<div lang="EN-GB" style="overflow-wrap: break-word;">
<div class="m_-1915382970844817410WordSection1">
<p class="MsoNormal">Hi list,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I tried a polynomial transformation published in 1925 and it works fine for the forward transformation:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Consolas">(echo 30000.00 100000.00 & echo 30000.00 120000.00) | cct -d 2 -z 0 +proj=horner +deg=3 +fwd_origin=0,0 +fwd_u=+120025.8705,+0.999986846,-0.0250515e-8,+0.1471965e-14,+0.006918172,-0.0860518e-8,-0.0398590e-14,+0.0011249e-8,+0.4556545e-14,+0.1991105e-14
+fwd_v=+390181.7849,+0.999993024,-0.0338879e-8,+0.1184845e-14,-0.006907254,-0.0492600e-8,-0.1945916e-14,+0.0222878e-8,+0.9244571e-14,-0.1091684e-14 +inv_tolerance=0.001
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Consolas"> 150717.96 489970.61 0.00 inf<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Consolas"> 150857.90 509969.45 0.00 inf<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">However, it gives an error for any input with a Northing greater than 500000 for the inverse transformation:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Consolas">(echo 150717.96 489970.61 & echo 150857.90 509969.45) | cct -d 2 -z 0 +proj=horner +deg=3 +fwd_origin=0,0 +fwd_u=+120025.8705,+0.999986846,-0.0250515e-8,+0.1471965e-14,+0.006918172,-0.0860518e-8,-0.0398590e-14,+0.0011249e-8,+0.4556545e-14,+0.1991105e-14
+fwd_v=+390181.7849,+0.999993024,-0.0338879e-8,+0.1184845e-14,-0.006907254,-0.0492600e-8,-0.1945916e-14,+0.0222878e-8,+0.9244571e-14,-0.1091684e-14 +inv_tolerance=0.001 +inv<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Consolas"> 30000.00 100000.00 0.00 inf<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Consolas"># Record 1 TRANSFORMATION ERROR: 150857.90 509969.45<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Consolas">(Point outside of projection domain)<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Is this a problem with the original polynomial transformation or a bug in PROJ?<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Regards, Jochem<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></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>
[<a href="https://www.kadaster.nl/algemene-leveringsvoorwaarden" target="_blank">https://www.kadaster.nl/algemene-leveringsvoorwaarden</a>].<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>
[<a href="https://www.kadaster.com/general-terms-and-conditions" target="_blank">https://www.kadaster.com/general-terms-and-conditions</a>]. <br>
</font>
</div>
_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
</div></blockquote></div>