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@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@lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/proj >
_______________________________________________ PROJ mailing list PROJ@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/proj
