[PROJ] Unexpected error in +proj=horner +inv

[Javier Jimenez Shaw]

Hi Jochem

After looking at the code, this error happens when the coordinates are out
of range. The default value for that range is 500000
You can change it with
https://proj.org/en/stable/operations/transformations/horner.html#cmdoption-arg-range

I hope it helps

On Wed, 11 Dec 2024 at 11:54, Lesparre, Jochem via PROJ <
proj at lists.osgeo.org> wrote:

> Hi list,
>
>
>
> I tried a polynomial transformation published in 1925 and it works fine
> for the forward transformation:
>
>
>
> (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
>
>     150717.96      489970.61          0.00           inf
>
>     150857.90      509969.45          0.00           inf
>
>
>
> However, it gives an error for any input with a Northing greater than
> 500000 for the inverse transformation:
>
>
>
> (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
>
>      30000.00      100000.00          0.00           inf
>
> # Record 1 TRANSFORMATION ERROR: 150857.90 509969.45
>
> (Point outside of projection domain)
>
>
>
> Is this a problem with the original polynomial transformation or a bug in
> PROJ?
>
>
>
> Regards, Jochem
>
>
>
>
> Disclaimer:
> De inhoud van deze e-mail is vertrouwelijk en uitsluitend bestemd voor de
> geadresseerde(n).
> Gebruik, openbaarmaking, vermenigvuldiging, verspreiding en/of
> verstrekking van deze informatie aan derden is niet toegestaan.
> Op al onze producten en diensten zijn onze algemene leveringsvoorwaarden
> van toepassing
> [https://www.kadaster.nl/algemene-leveringsvoorwaarden].
>
> Disclaimer:
> 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.
> 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.
> Our general terms and conditions of delivery apply to all our products and
> services
> [https://www.kadaster.com/general-terms-and-conditions].
> _______________________________________________
> PROJ mailing list
> PROJ at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/proj
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20241211/1374b1c3/attachment.htm>