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

Vr Okt 3 12:19:11 PDT 2014

Hallo Bas,

Er is nog geen EPSG-code voor deze datumtransformatie, maar er is wel een brainstorm geweest met EPSG over hoe de RD-defintie in het algemeen het beste kan worden opgenomen, hierbij is ook het NTv2-grid besproken. Voor het NTv2-grid is het een beetje een kip-ei probleem. Pas nu er een release is van de grids heeft EPSG iets om naar te refereren, terwijl gebruikers liever nu al een EPSG-code hadden gehad.

Overigens blijft de RD-projectie gewoon EPSG-28992. 

In het bestand: http://svn.osgeo.org/metacrs/proj/trunk/proj/nad/epsg (ik neem aan dat de versie is die voor PROJ.4 4.9.0 wordt gebruikt) kom ik de volgende string tegen voor epsg:28992
# Amersfoort / RD New
<28992> +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 +units=m +no_defs  <>

De +towgs84 optie in deze string heeft niets met de RD-projectie (en dus EPSG-28992) te maken, dit geldt ook voor de +nadgrids optie. Beide opties zijn bedoeld voor datumtransformaties, zie ook https://trac.osgeo.org/proj/wiki/GenParms punt 7. en 8.

De optie +nadgrids in de strings die ik heb gegeven kun je m.i. gewoon gebruiken in een nieuwe release van PROJ.4 als vervanging van de +towgs84 optie (deze twee opties kunnen sowieso niet gelijktijdig gebruikt in PROJ.4).

De string wordt dan:
# RD-projection
<28992> +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +nadgrids=rdtrans2008.gsb +units=m +no_defs  <>

Analoog kun je dit doen voor voor epsg:7415, daarvoor staat nu:
# Amersfoort / RD New + NAP height
<7415> +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 +units=m +vunits=m +no_defs  <>

Met NTv2 en VDatum wordt dit:
# RD-projection + NAP-height
<7415> +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 +vunits=m +no_defs  <>

Groeten,
Lennard

P.S.  nog wat opmerkingen over RD in het huidige bestand http://svn.osgeo.org/metacrs/proj/trunk/proj/nad/epsg
1)
De +towgs84 optie wordt dus gebruikt voor datumtransformaties, getallen die je in verschillende blogs, forums en het epsg bestand in PROJ.4 tegenkomt in combinatie met de RD-projectie (EPSG-28992) zijn bedoeld als de transformatieparameters tussen het ETRS89 en het RD. Voor deze transformatie heeft het Kadaster in 1996, 2005 en 2009 parameters gepubliceert. Deze verschillende versies van de parameters hebben ook een EPSG-code, respectievelijk en 1751, 15739 en 4830. EPSG vermeldt bij deze parameters overigens dat deze een nauwkeurigheid hebben van 50 centimeter, dat lees ik nergens terug in PROJ.4.

De huidige waarden bij de +towgs84 optie in zijn afgerond op basis van de 2009 waarden. De juiste waarden voor de +towgs84 optie bij gebruik van RD en ETRS89 zijn:

tot en met 2004:
+towgs84=565.040,49.910,465.840,-0.409394,0.359705,-1.868491,4.0772

van 2005 tot en met 2008:
+towgs84=565.2369,50.0087,465.6580,-0.406857,0.350733,-1.870347,4.0812

vanaf 2009:
+towgs84=565.4171,50.3319,465.5524,-0.398957,0.343988,-1.877402,4.0725

2)
De namen Amersfoort en RD New zijn weliswaar ingeburgerd bij GIS-gebruikers, maar komen niet in de defintie van RD voor. Het geheel is het stelsel van de Rijksdriehoeksmeting (kortweg RD-stelsel). Het RD-stelsel is gebaseerd op de Bessel-ellipsoide uit 1841 met als centraal punt Amersfoort, dit noemen we 'RD (Bessel)' (en niet Amersfoort). De kaartprojectie die wordt gebruik voor RD is de dubbelprojectie van Schreiber, dit noemen we de RD-projectie (en niet RD New).

________________________________________
From: dutch-bounces op lists.osgeo.org [dutch-bounces op lists.osgeo.org] on behalf of Sebastiaan Couwenberg [sebastic op xs4all.nl]
Sent: Friday, October 03, 2014 3:34 PM
To: dutch op lists.osgeo.org
Subject: Re: [Dutch] RD/NAP naar ETRS89 met NTv2 en VDatum

On 10/03/2014 11:43 AM, Huisman, Lennard wrote:
> In het tijdschrift Geo-Info van oktober 2014 staat een artikel over verschillende transformatie stussen RD en ETRS89. Gelijktijdig hiermee maken we een NTv2-grid en een VDatum-grid beschikbaar voor een nauwkeurige transformatie tussen RD/NAP en ETRS89, die op mm-niveau gelijk is aan de officiele RDNAPTRANS(TM)2008 procedure.
>
> De grids en wat achtergrond informatie zijn te vinden via de Kadaster website:
> https://www.kadaster.nl/web/Themas/Registraties/Rijksdriehoeksmeting/Transformatie-van-coordinaten.htm
>
> Als iemand de grids wil opnemen in het lijstje bij 'Datum shift grids' op http://trac.osgeo.org/proj/ juich ik dat toe.
>
> Met betrekking tot het gebruik:
> De string voor definiëren van RD/NAP met de NTv2- en VDatum-grid in proj4:
>
> # RDNAP with NTv2 and VDatum
> <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 <>
>
> # RD with NTv2 only
> <rd> +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +nadgrids=rdtrans2008.gsb +units=m +no_defs <>

Zijn er EPSG codes toegewezen aan deze projecties?

Recent is de ontwikkeling van PROJ.4 weer opgepikt en zijn er na bijna
een jaar sinds beta2 twee release candidates voor versie 4.9.0 uitgebracht.

Mogelijk kunnen de wijzigingen voor RD/NAP nog opgenomen worden voordat
4.9.0 final word uitgebracht.

Mvg,

Bas

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

_______________________________________________
Dutch mailing list
Dutch op lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/dutch

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