[Dutch] RD/NAP naar ETRS89 met NTv2 en VDatum

[Huisman, Lennard]

Beste Sebastiaan,

Ik heb hetzelfde tesje gedaan met OpenSuSe 12.3 en cs2cs 4.8.0. De grids en het rdnap bestand heb ik gedownloaded vanuit git en in de map /usr/share/proj/epsg geplaatst:

~> echo "155000 463000 0" | PROJ_DEBUG=3 /usr/bin/cs2cs -v -s +init=rdnap:rdnap +to +init=epsg:4258 -E -f %.9f
pj_open_lib(rdnap): call fopen(/usr/share/proj/rdnap) - succeeded

pj_open_lib(epsg): call fopen(/usr/share/proj/epsg) - succeeded

# ---- From Coordinate System ----
#Oblique Stereographic Alternative
#	Azimuthal, Sph&Ell
# +init=rdnap:rdnap +proj=sterea +lat_0=52.15616055555555
# +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel
# +nadgrids=rdtrans2008.gsb +geoidgrids=naptrans2008.gtx +units=m +no_defs
# ---- To Coordinate System ----
#Lat/long (Geodetic alias)
#	
# +init=epsg:4258 +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
pj_open_lib(naptrans2008.gtx): call fopen(/usr/share/proj/naptrans2008.gtx) - succeeded

GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(/usr/share/proj/naptrans2008.gtx) - succeeded

pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(/usr/share/proj/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(/usr/share/proj/rdtrans2008.gsb) - succeeded

pj_apply_gridshift(): used NL_LAND 
155000 463000 0	52.155172899	5.387203658 43.255106923

ALs ik het rdtrans2008.gsb grid verwijder uit de map /usr/share/proj/epsg krijg ik het volgende:

~> echo "155000 463000 0" | PROJ_DEBUG=3 /usr/bin/cs2cs -v -s +init=rdnap:rdnap +to +init=epsg:4258 -E -f %.9f
pj_open_lib(rdnap): call fopen(/usr/share/proj/rdnap) - succeeded

pj_open_lib(epsg): call fopen(/usr/share/proj/epsg) - succeeded

# ---- From Coordinate System ----
#Oblique Stereographic Alternative
#	Azimuthal, Sph&Ell
# +init=rdnap:rdnap +proj=sterea +lat_0=52.15616055555555
# +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel
# +nadgrids=rdtrans2008.gsb +geoidgrids=naptrans2008.gtx +units=m +no_defs
# ---- To Coordinate System ----
#Lat/long (Geodetic alias)
#	
# +init=epsg:4258 +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
pj_open_lib(naptrans2008.gtx): call fopen(/usr/share/proj/naptrans2008.gtx) - succeeded

GTX 311x641: LL=(2.5,50.5) UR=(7.66666667,55.8333333)
pj_open_lib(naptrans2008.gtx): call fopen(/usr/share/proj/naptrans2008.gtx) - succeeded

pj_apply_gridshift(): used GTX Vertical Grid Shift File
pj_open_lib(rdtrans2008.gsb): call fopen(/usr/share/proj/rdtrans2008.gsb) - failed

Rel. 4.8.0, 6 March 2012
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): failed to load datum shift file
155000 463000 0	*	* 43.255041641

Wel een foutmelding dus.

Naast proj heb ik ook gdal en qgis geinstalleerd staan, ik weet niet of dat invloed heeft op de tests.

Gr,
Lennard
______________________________________
From: Sebastiaan Couwenberg [sebastic op xs4all.nl]
Sent: Friday, October 24, 2014 3:55 PM
To: Huisman, Lennard; dutch op lists.osgeo.org
Subject: Re: [Dutch] RD/NAP naar ETRS89 met NTv2 en VDatum

On 10/20/2014 11:37 PM, Huisman, Lennard wrote:
> Beste Sebastiaan,
>
>>> Mijn test resultaten heb ik online gezet voor geïnteresseerden:
>
> Ik heb de testresultaten bekeken. Voor x,y en phi,lambda zijn deze gelijk aan de resultaten die je krijgt als je optie +nadgrids weglaat uit de projectiestring. Het rdtrans2008.gsb correctiegrid wordt blijkbaar niet gebruikt/gevonden tijdens de test. Ik heb de indruk at het naptrans2008.gtx grid wel wordt gebruikt, de hoogte bij conversie van rdnap naar etrs89 klopt namelijk.
>
> Ter illustratie de output die ik krijg (proj 4.8.0 op windows):
>
> cs2cs -s +init=rdnap:nordnap +to +init=epsg:4258 -E -f %.9f
> 155000 463000 0
> 155000 463000 0 52.156160556    5.387638889 43.255041641
>
> Hierbij is in de 'nordnap' string de optie +nadgrids=rdtrans2008.gsb weggelaten. Het resultaat is gelijk aan die van het testscript
>
> cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -E -f %.9f
> 155000 463000 0
> 155000 463000 0 52.155172899    5.387203658 43.255106923
>
> Dit is de juiste uitkomst. De 'rdnap' string heb ik ge-copy-pasted uit het rdnap bestand in de git repository.

De rdtrans2008.gsb wordt inderdaad niet gebruikt, pj_open_lib word er
nooit voor aangeroepen:

$ echo "155000 463000 0" | PROJ_DEBUG=3 PROJ_LIB=. /usr/bin/cs2cs -v -s
+init=rdnap:rdnap +to +init=epsg:4258 -E -f %.9f
pj_open_lib(rdnap): call fopen(./rdnap) - succeeded

pj_open_lib(epsg): call fopen(./epsg) - succeeded

# ---- From Coordinate System ----
#Oblique Stereographic Alternative
#       Azimuthal, Sph&Ell
# +init=rdnap:rdnap +proj=sterea +lat_0=52.15616055555555
# +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel
# +nadgrids=rdtrans2008.gsb +geoidgrids=naptrans2008.gtx +units=m +no_defs
# ---- To Coordinate System ----
#Lat/long (Geodetic alias)
#
# +init=epsg:4258 +proj=longlat +ellps=GRS80 +no_defs
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
155000 463000 0 52.156160556    5.387638889 43.255041641


Als rdtrans2008.gsb hernoemd word en dus niet gevonden kan worden, word
er ook geen error geproduceerd voor het ontbreken van een niet optionele
datum grid in PROJ_LIB.

Interessant dat proj 4.8.0 op Windows het bij jou wel goed doet, maar
bij mij op Linux niet. Ik heb het ook op een Windows 7 64bit VM getest
en daar krijg ik hetzelfde resultaat als op 64bit Linux. Zowel met 32bit
als 64bit Cygwin i.c.m. respectievelijk OSGeo4W en OSGeo4W64 word de
rdtrans2008.gsb ook niet gebruikt.

Ik vermoed dat een van de voorgaande parameters ervoor zorgt dat de
nadgrids parameter word overgeslagen, de geoidgrids parameter word wel
gewoon gebruikt.

>>> In de git repository is de patch voor de test te vinden:
>
> In de tabel van het test script zit nog een kleine typo, ETS89 moet zijn ETRS89.

Good catch, fixed in git.

> Gr,
> Lennard

Mvg,

Bas

--
 GPG Key ID: 4096R/E88D4AF1
Fingerprint: 8182 DE41 7056 408D 6146  50D1 6750 F10A E88D 4AF1



Disclaimer:
De inhoud van dit bericht is uitsluitend bestemd voor geadresseerde.
Gebruik van de inhoud van dit bericht door anderen zonder toestemming van het Kadaster 
is onrechtmatig. Mocht dit bericht ten onrechte bij u terecht komen, dan verzoeken wij u 
dit direct te melden aan de verzender en het bericht te vernietigen. 
Aan de inhoud van dit bericht kunnen geen rechten worden ontleend.

Disclaimer:
The content of this message is meant to be received by the addressee only.
Use of the content of this message by anyone other than the addressee without the consent 
of the Kadaster is unlawful. If you have received this message, but are not the addressee, 
please contact the sender immediately and destroy the message.
No rights can be derived from the content of this message