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

[Sebastiaan Couwenberg]

Beste Lennard,

Hartelijk dank voor de uitvoerige informatie.

On 10/03/2014 09:19 PM, Huisman, Lennard wrote:
> 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  <>

Ik was benieuwt waar ik de de rdtrans2008.gsb & naptrans2008.gtx
bestanden konden, om te zien of deze in het PROJ.4 project opgenomen
kunnen worden.

Na het aanvragen van de RDNAPTRANS2008 download via
https://www.kadaster.nl/web/Themas/Registraties/Rijksdriehoeksmeting/Transformatie-van-coordinaten.htm
vond ik de bestanden in de NTv2_2008.zip file.

Helaas staat de licentie van RDNAPTRANS2008 herdistributeren niet toe,
waardoor de C++ voorbeeld code en de grid files niet in open source
projecten opgenomen kunnen worden.

De licentie voor NTv2_2008 is ook problematisch voor open source
projecten. Het niet mogen aanpassen van de correctie waardes zorgt
ervoor dat de licentie niet voldoet aan Debian Free Software Guidelines
en daarvan afgeleide Open Source Definition.

Ik snap de gedachte achter het niet aanpassen wel, net als de OSG
Document Notice voor de standard documenten. Voor Debian gebaseerde
distributies zou een source package in de non-free sectie vereist zijn
die rdtrans2008.gsb & naptrans2008.gtx bestanden beschikbaar maakt voor
de proj package.

De voorgestelde aanpassingen aan de epsg file kunnen onder de huidige
licentie voorwaarden dus niet het proj project zelf opgenomen worden
omdat de MIT licentie expliciet het recht op aanpassingen toestaat.

In de documentatie een notitie opnemen dat de transformatie procedure
niet identiek is aan RDNAPTRANS™2008 is mogelijk, gezien het een
"should" requirement is is het waarschijnlijk geen harde eis.

De verplichting de grid correcties te testen zorgt er mogelijk voor dat
tijdens het build proces de test fase niet overgeslagen kan worden bij
een strikte interpretatie.

Is het mogelijk de NTv2_2008 licentie aan te passen naar een
gebruikelijke Free Software licentie (of zelf public domain/CC-0) zodat
het in open source projecten opgenomen kan worden?

> 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 EPSG codes 1751, 15739 en 4830 komen niet voor in de epsg file van
PROJ.4. Wel deze:

# Amersfoort
<4289> +proj=longlat +ellps=bessel
+towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725
+no_defs  <>

# Amersfoort / RD Old
<28991> +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889
+k=0.9999079 +x_0=0 +y_0=0 +ellps=bessel
+towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725
+units=m +no_defs  <>

# 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  <>

> 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

Deze +towgs84 waardes zijn net iets preciezer dan in de epsg file
momenteel, dat kan overgenomen worden.

> 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).

Zou het wenselijk zijn om de naam in de epsg file aan te passen naar
RD-projection?

Mvg,

Bas

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