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

[Sebastiaan Couwenberg]

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