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

Za Okt 4 00:31:58 PDT 2014

Hoi Sebastian,

Het NTV2 grid- en vdatum-bestand mag je volgens mij wel distribueren, 
zie deze Engelstalige quote uit de beschrijving:

"License
RDNAP makes the rdtrans2008.gsb and naptrans.gtx grid correction files 
freely available.
Implementation and distribution of the rdtrans2008 NTv2-grid and 
naptrans2008 VDatum-grid in
combination with other tools/software that require these grids for 
transformation between ETRS89 and
RD/NAP is allowed, under the following conditions:
1 . It is not allowed to modify the correction values of the grid 
correction files. (It is allowed to
convert the grid data to a different file format or data structure as 
long as this does not affect
the results of the calculations.)
2. It should be clear from the documentation of the software using the 
grid correction files that
the transformation procedure using the grid correction files is not 
identical to
RDNAPTRANSTM2008.
3. The correct implementation of the grid corrections has to be tested 
as specified in Appendix 2"

Voor de gridcorrectie kun je gebruik maken van bestaande software (in 
Java kun je bijv. JGridShift via GeoTools gebruiken) die wel een 
geschikte licentie heeft.

Groeten,

Jan

On 4-10-2014 1:29, Sebastiaan Couwenberg wrote:
> 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
>

-- 
----------------------------------------
Geodan
Buitenhaven 27-A
5211 TP 's-Hertogenbosch (NL)

T +31 (0)73 - 692 5151
M +31 (0)6 - 2395 1191
E jan.boonen at geodan.nl
www.geodan.nl | www.geodan.nl/disclaimer
----------------------------------------