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
