PROJ.4 4.9.2 and points outside datum shift grid
Debian Bug #813066 [0] revealed a bug in the test handling of the
conversion between RD/NAP and ETRS89 when PROJ.4 4.9.2 is used.
Previous versions of PROJ.4 were unable to transform points outside the
available datum shift grid, as can be seen in the output of the test
script [1] in the attachment testrdtrans2008-4.9.1.
With version 4.9.2 some points can now be calculated as seen in the
testrdtrans2008-4.9.2 attachment, but for tests 07 to 10 (whose points
"are outside the region where interpolation between either the NLGEO2004
geoid or the RD correction grid points is possible") in the direction
from RD/NAP to ETRS89 the results are entirely incorrect, whereas the
results for those same tests in opposite direction are mostly correct
except test 10 (edge_rd).
The documentation accompanying the rdtrans2008.gsb NTv2-grid and
naptrans2008.gtx VDatum-grid [2] mentions that the points for test 07 to
10 are outside the correction grid and that their "RD coordinates can be
computed, but the output should be handled with care".
With earlier versions there was no usable output, 4.9.2 reports
incorrect values (for some tests).
It looks like the fix for #292 [3][4] is the cause of this change. I'm
not sure whether this change is a bug or a feature.
In de debug logging for 4.9.2 we can see that for tests 07 to 10
pj_apply_vgridshift is no longer called, instead pj_apply_gridshift is
used. For tests 07 to 10 no z value is specified, so this looks like the
intended course of action to not apply the vertical grid shift.
On the other hand, the previous behaviour caused pj_apply_vgridshift to
log correctly that it failed to find a grid shift table for the location.
How should points outside outside datum shift grids be handled by PROJ.4?
The previous behaviour seems more correct, but with 4.9.2 some
transformations can be done that weren't before.
[0] https://bugs.debian.org/813066
[1]
https://sources.debian.net/src/proj-rdnap/2008-4/debian/patches/add-test.patch/
[2]
https://sources.debian.net/data/non-free/p/proj-rdnap/2008-4/Use%20of%20RDTRANS2008%20and%20NAPTRANS2008.pdf
[3] https://github.com/OSGeo/proj.4/issues/292
[4]
https://github.com/OSGeo/proj.4/commit/5d2af3a89be0898458e4f0fe272affc36642d304
Kind Regards,
Bas
--
GPG Key ID: 4096R/6750F10AE88D4AF1
Fingerprint: 8182 DE41 7056 408D 6146 50D1 6750 F10A E88D 4AF1
Direction: From ETRS89 to RD/NAP
Test: 01 Texel
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 53.160753042 4.824761912 42.8614
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 117380.1203 575040.3398 1.0000
Expect: 117380.1200 575040.3400 1.0000
Test OK: From ETRS89 to RD/NAP - 01 Texel (Not identical, but within margin)
Test: 02 Noord-Groningen
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 53.419482050 6.776726674 42.3586
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 247380.5600 604580.7798 1.9999
Expect: 247380.5600 604580.7800 2.0000
Test OK: From ETRS89 to RD/NAP - 02 Noord-Groningen (Not identical, but within margin)
Test: 03 Amersfoort
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 52.155172897 5.387203657 43.2551
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 155000.0000 462999.9998 -0.0000
Expect: 155000.0000 463000.0000 0.0000
Test OK: From ETRS89 to RD/NAP - 03 Amersfoort (Not identical, but within margin)
Test: 04 Zeeuws-Vlaanderen
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 51.368607152 3.397588595 47.4024
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 16460.9098 377380.2298 3.0000
Expect: 16460.9100 377380.2300 3.0000
Test OK: From ETRS89 to RD/NAP - 04 Zeeuws-Vlaanderen (Not identical, but within margin)
Test: 05 Zuid-Limburg
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 50.792584908 5.773795547 174.9478
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 182260.4500 311480.6701 129.0000
Expect: 182260.4500 311480.6700 129.000
Test OK: From ETRS89 to RD/NAP - 05 Zuid-Limburg (Not identical, but within margin)
Test: 06 Maasvlakte
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 51.947393898 4.072887101 47.5968
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 64640.8899 440700.0099 3.9998
Expect: 64640.8900 440700.0101 4.0000
Test OK: From ETRS89 to RD/NAP - 06 Maasvlakte (Not identical, but within margin)
Test: 07 No_rd&geoid
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 50.687420405 4.608971812
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_vgridshift(): failed to find a grid shift table for
location (4.6093016dW,50.6882354dN)
tried: naptrans2008.gtx
Rel. 4.9.1, 04 March 2015
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
Output: * * -0.0001
Expect: 100000.6700 300000.8900
Test OK: From ETRS89 to RD/NAP - 07 No_rd&geoid (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 08 No_geoid
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 51.136825197 4.601375361
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_vgridshift(): failed to find a grid shift table for
location (4.6017018dW,51.1376928dN)
tried: naptrans2008.gtx
Rel. 4.9.1, 04 March 2015
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
Output: * * -0.0001
Expect: 100000.6700 350000.8900
Test OK: From ETRS89 to RD/NAP - 08 No_geoid (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 09 No_rd
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 52.482440839 4.268403889
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_vgridshift(): failed to find a grid shift table for
location (4.2686722dW,52.4834651dN)
tried: naptrans2008.gtx
Rel. 4.9.1, 04 March 2015
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
Output: * * -0.0001
Expect: 79000.0100 500000.2300
Test OK: From ETRS89 to RD/NAP - 09 No_rd (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 10 edge_rd
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 51.003976532 3.891247830
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_vgridshift(): failed to find a grid shift table for
location (3.8914723dW,51.0048262dN)
tried: naptrans2008.gtx
Rel. 4.9.1, 04 March 2015
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
Output: * * -0.0001
Expect: 50000.4500 335999.6700
Test OK: From ETRS89 to RD/NAP - 10 edge_rd (Expected output overriden: * * ^-?\d+\.\d+$)
Direction: From RD/NAP to ETRS89
Test: 01 Texel
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 117380.1200 575040.3400 1.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 53.160753044 4.824761908 42.861397590
Expect: 53.160753042 4.824761912 42.8614
Test OK: From RD/NAP to ETRS89 - 01 Texel (Not identical, but within margin)
Test: 02 Noord-Groningen
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 247380.5600 604580.7800 2.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 53.419482052 6.776726674 42.358705049
Expect: 53.419482050 6.776726674 42.3586
Test OK: From RD/NAP to ETRS89 - 02 Noord-Groningen (Not identical, but within margin)
Test: 03 Amersfoort
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 155000.0000 463000.0000 0.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 52.155172899 5.387203658 43.255106923
Expect: 52.155172897 5.387203657 43.2551
Test OK: From RD/NAP to ETRS89 - 03 Amersfoort (Not identical, but within margin)
Test: 04 Zeeuws-Vlaanderen
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 16460.9100 377380.2300 3.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 51.368607153 3.397588597 47.402397195
Expect: 51.368607152 3.397588595 47.4024
Test OK: From RD/NAP to ETRS89 - 04 Zeeuws-Vlaanderen (Not identical, but within margin)
Test: 05 Zuid-Limburg
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 182260.4500 311480.6700 129.000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 50.792584907 5.773795547 174.947778837
Expect: 50.792584908 5.773795547 174.9478
Test OK: From RD/NAP to ETRS89 - 05 Zuid-Limburg (Not identical, but within margin)
Test: 06 Maasvlakte
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 64640.8900 440700.0101 4.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 51.947393899 4.072887103 47.596976393
Expect: 51.947393898 4.072887101 47.5968
Test OK: From RD/NAP to ETRS89 - 06 Maasvlakte (Not identical, but within margin)
Test: 07 No_rd&geoid
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 100000.6700 300000.8900
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_vgridshift(): failed to find a grid shift table for
location (4.6093016dW,50.6882354dN)
tried: naptrans2008.gtx
Rel. 4.9.1, 04 March 2015
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
Output: * * 0.000000000
Expect: 50.687420405 4.608971812
Test OK: From RD/NAP to ETRS89 - 07 No_rd&geoid (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 08 No_geoid
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 100000.6700 350000.8900
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_vgridshift(): failed to find a grid shift table for
location (4.6017018dW,51.1376928dN)
tried: naptrans2008.gtx
Rel. 4.9.1, 04 March 2015
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
Output: * * 0.000000000
Expect: 51.136825197 4.601375361
Test OK: From RD/NAP to ETRS89 - 08 No_geoid (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 09 No_rd
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 79000.0100 500000.2300
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_vgridshift(): failed to find a grid shift table for
location (4.2686722dW,52.4834651dN)
tried: naptrans2008.gtx
Rel. 4.9.1, 04 March 2015
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
Output: * * 0.000000000
Expect: 52.482440839 4.268403889
Test OK: From RD/NAP to ETRS89 - 09 No_rd (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 10 edge_rd
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 50000.4500 335999.6700
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_vgridshift(): failed to find a grid shift table for
location (3.8914729dW,51.0048265dN)
tried: naptrans2008.gtx
Rel. 4.9.1, 04 March 2015
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
Output: * * 0.000000000
Expect: 51.003976532 3.891247830
Test OK: From RD/NAP to ETRS89 - 10 edge_rd (Expected output overriden: * * ^-?\d+\.\d+$)
Direction: From ETRS89 to RD/NAP
Test: 01 Texel
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 53.160753042 4.824761912 42.8614
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 117380.1203 575040.3398 1.0000
Expect: 117380.1200 575040.3400 1.0000
Test OK: From ETRS89 to RD/NAP - 01 Texel (Not identical, but within margin)
Test: 02 Noord-Groningen
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 53.419482050 6.776726674 42.3586
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 247380.5600 604580.7798 1.9999
Expect: 247380.5600 604580.7800 2.0000
Test OK: From ETRS89 to RD/NAP - 02 Noord-Groningen (Not identical, but within margin)
Test: 03 Amersfoort
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 52.155172897 5.387203657 43.2551
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 155000.0000 462999.9998 -0.0000
Expect: 155000.0000 463000.0000 0.0000
Test OK: From ETRS89 to RD/NAP - 03 Amersfoort (Not identical, but within margin)
Test: 04 Zeeuws-Vlaanderen
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 51.368607152 3.397588595 47.4024
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 16460.9098 377380.2298 3.0000
Expect: 16460.9100 377380.2300 3.0000
Test OK: From ETRS89 to RD/NAP - 04 Zeeuws-Vlaanderen (Not identical, but within margin)
Test: 05 Zuid-Limburg
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 50.792584908 5.773795547 174.9478
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 182260.4500 311480.6701 129.0000
Expect: 182260.4500 311480.6700 129.000
Test OK: From ETRS89 to RD/NAP - 05 Zuid-Limburg (Not identical, but within margin)
Test: 06 Maasvlakte
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 51.947393898 4.072887101 47.5968
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 64640.8899 440700.0099 3.9998
Expect: 64640.8900 440700.0101 4.0000
Test OK: From ETRS89 to RD/NAP - 06 Maasvlakte (Not identical, but within margin)
Test: 07 No_rd&geoid
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 50.687420405 4.608971812
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 100000.6700 300000.8899 2147479807.9999
Expect: 100000.6700 300000.8900
Test OK: From ETRS89 to RD/NAP - 07 No_rd&geoid (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 08 No_geoid
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 51.136825197 4.601375361
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 100000.6699 350000.8898 2147479935.9999
Expect: 100000.6700 350000.8900
Test OK: From ETRS89 to RD/NAP - 08 No_geoid (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 09 No_rd
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 52.482440839 4.268403889
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 79000.0100 500000.2300 2147479935.9999
Expect: 79000.0100 500000.2300
Test OK: From ETRS89 to RD/NAP - 09 No_rd (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 10 edge_rd
Exec: cs2cs -r +init=epsg:4258 +to +init=rdnap:rdnap -f %.4f
Input: 51.003976532 3.891247830
STDERR: pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
Output: 50000.4133 335999.6363 903779406.6674
Expect: 50000.4500 335999.6700
x coordinate exceeds 0.001 meters: 0.0366999999969266
y coordinate exceeds 0.001 meters: 0.0336999999708496
Test FAILED: From ETRS89 to RD/NAP - 10 edge_rd (Expected output overriden: * * ^-?\d+\.\d+$)
Direction: From RD/NAP to ETRS89
Test: 01 Texel
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 117380.1200 575040.3400 1.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 53.160753044 4.824761908 42.861397590
Expect: 53.160753042 4.824761912 42.8614
Test OK: From RD/NAP to ETRS89 - 01 Texel (Not identical, but within margin)
Test: 02 Noord-Groningen
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 247380.5600 604580.7800 2.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 53.419482052 6.776726674 42.358705049
Expect: 53.419482050 6.776726674 42.3586
Test OK: From RD/NAP to ETRS89 - 02 Noord-Groningen (Not identical, but within margin)
Test: 03 Amersfoort
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 155000.0000 463000.0000 0.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 52.155172899 5.387203658 43.255106923
Expect: 52.155172897 5.387203657 43.2551
Test OK: From RD/NAP to ETRS89 - 03 Amersfoort (Not identical, but within margin)
Test: 04 Zeeuws-Vlaanderen
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 16460.9100 377380.2300 3.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 51.368607153 3.397588597 47.402397195
Expect: 51.368607152 3.397588595 47.4024
Test OK: From RD/NAP to ETRS89 - 04 Zeeuws-Vlaanderen (Not identical, but within margin)
Test: 05 Zuid-Limburg
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 182260.4500 311480.6700 129.000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 50.792584907 5.773795547 174.947778837
Expect: 50.792584908 5.773795547 174.9478
Test OK: From RD/NAP to ETRS89 - 05 Zuid-Limburg (Not identical, but within margin)
Test: 06 Maasvlakte
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 64640.8900 440700.0101 4.0000
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: 51.947393899 4.072887103 47.596976393
Expect: 51.947393898 4.072887101 47.5968
Test OK: From RD/NAP to ETRS89 - 06 Maasvlakte (Not identical, but within margin)
Test: 07 No_rd&geoid
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 100000.6700 300000.8900
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: -50.688542790 -175.391028187 2134749119.469262600
Expect: 50.687420405 4.608971812
Latitude not expected: -50.688542790, expected: *
Longitude not expected: -175.391028187, expected: *
Test FAILED: From RD/NAP to ETRS89 - 07 No_rd&geoid (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 08 No_geoid
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 100000.6700 350000.8900
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: -51.137943931 -175.398624638 2134749575.923032999
Expect: 51.136825197 4.601375361
Latitude not expected: -51.137943931, expected: *
Longitude not expected: -175.398624638, expected: *
Test FAILED: From RD/NAP to ETRS89 - 08 No_geoid (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 09 No_rd
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 79000.0100 500000.2300
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: -52.483546993 -175.731596111 2134750552.422238827
Expect: 52.482440839 4.268403889
Latitude not expected: -52.483546993, expected: *
Longitude not expected: -175.731596111, expected: *
Test FAILED: From RD/NAP to ETRS89 - 09 No_rd (Expected output overriden: * * ^-?\d+\.\d+$)
Test: 10 edge_rd
Exec: cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -f %.9f
Input: 50000.4500 335999.6700
STDERR: pj_open_lib(rdnap): call fopen(./rdnap) - succeeded
pj_open_lib(epsg): call fopen(./epsg) - succeeded
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(./naptrans2008.gtx) - succeeded
pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
NTv2 NL_ALL 63x65: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
NTv2 NL_LAND 561x421: LL=(3,50.5) UR=(7.66666667,54)
NTv2 - loading grid NL_LAND
pj_open_lib(rdtrans2008.gsb): call fopen(./rdtrans2008.gsb) - succeeded
pj_apply_gridshift(): used NL_LAND
Output: -51.006648861 -176.108751657 890969128.758312702
Expect: 51.003976532 3.891247830
Latitude not expected: -51.006648861, expected: *
Longitude not expected: -176.108751657, expected: *
Test FAILED: From RD/NAP to ETRS89 - 10 edge_rd (Expected output overriden: * * ^-?\d+\.\d+$)
Reply to: