[Dutch] Fwd: RD/ETRS89 Kennissessie 17 juni | presentaties en rapport beschikbaar

[Richard Duivenvoorde]

On 03-07-15 12:56, Sebastiaan Couwenberg wrote:
> On 07/03/2015 12:22 PM, Richard Duivenvoorde wrote:
>> Wat is de status van GDAL/28992+correctiegrid?
>> Want als ik het goed begrijp, heb je daarmee dan eigenlijk een goed
>> uitgangspunt om evt die transformatie met GDAL te maken?
>>
>> Zoiezo zou het natuulijk ideaal zijn als die binnen de GDAL-lib zou
>> kunnen worden opgelost (of een andere Open Source lib), dan kunnen
>> andere gispakketten daarop inhaken met extenties e.d.
> 
> Volgens mij was de bedoeling van de rdnaptrans datumgrids om deze
> transformatie te ondersteunen. Zie de eerdere discussie op deze lijst:
> 
> RD/NAP naar ETRS89 met NTv2 en VDatum
> http://lists.osgeo.org/pipermail/dutch/2014-October/000975.html

Wow, ook een must-read bijna :-)

Petje af voor jou en Lennart hoor dat jullie dat zo goed hebben
uitgezocht EN geregeld!

Ben nog wel benieuwd naar de cliff-hanger:

wordt (de ontvoerde) rdtrans2008.gsb nu wel of niet gevonden door
(inspecteur) Debian ;-)

Ik zag de laatste aflevering niet in de thread, maar ik denk dat het
goed afgelopen is, omdat ik blijkbaar [0] de juiste uitkomst krijg:

cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -E -f %.9f
155000 463000 0
155000 463000 0	52.155172899	5.387203658 43.255106923



Blijft voor mij praktisch nog even de vraag aan de mensen die er
verstand van hebben:

Kan ik nu mijn kaas.shp bestand die ik hier heb staan als zijnde een
epsg:28992 (2D) bestand nu direct omzetten naar epsg:4258 (3D)?

cs2cs -s +init=epsg:28992 +to +init=epsg:4258 -E -f %.9f
155000 463000
155000 463000	52.155172302	5.387203508 43.346941932
^^^ lijkt goed, ik kom hier die 43 tegen uit [1]

Of moet ik dat doen via epsg:28992 (2D) naar (met een aanname van de
hoogte blijkbaar) naar epsg:4258 (3D)?
Ik neem dan de hoogte even 43 zie Appendix II in [1]

cs2cs -s +init=epsg:28992 +to +init=epsg:4258 -E -f %.9f
155000 463000 43
155000 463000 43	52.155172308	5.387203509 86.347117045
^^^ mmm blijkbaar ook niet, ik lig nu 2x zo hoog als ik wilde?

Of moet ik net doen of het mijn 28992 bestand een rndap:rd bestand is en
dus van rdnap:rd naar epsg:4258 gaan?
cs2cs -s +init=rdnap:rd +to +init=epsg:4258 -E -f %.9f
155000 463000
155000 463000	52.155172899	5.387203658 0.000065284
^^^ mmm, zou kunnen

Of moet ik dat doen via epsg:28992 (2D) naar rdnap:rdnap (3D) en dan
naar epsg:4258 (3D)?

cs2cs -s +init=rdnap:rd +to +init=rdnap:rdnap -E -f %.9f
155000 463000
155000 463000	463000.000000000	155000.000000000 -43.255041641

cs2cs -s +init=rdnap:rdnap +to +init=epsg:4258 -E -f %.9f
155000.000000000  463000.000000000 -43.255041641
155000.000000000  463000.000000000 -43.255041641	52.155172899
5.387203658 0.000065284
^^^ hee, dezelfde als hierboven, maar 43 meter beneden het 'juiste' antwoord


Of mis ik toch nog iets?
Ik verwacht dus eigenlijk dat het antwoord iets zal zijn als: "je moet
je zelf afvragen wat voor data je had toen je 155000,463000 vastlegde
als een punt. Op ANP nivo == 0 of op de hoogte zoals die in het
hoogtegrid zit (V-Datum?) of ....

Groet,

Richard Duivenvoorde

ps, feel free to ignore, just thinking out loud

[0] https://lists.osgeo.org/pipermail/dutch/2014-October/000992.html
[1]
http://www.geonovum.nl/sites/default/files/rd_etrs89_resultaten_impactanalyse.pdf