Rijksdriehoekstelsel (epsg:28992) verschuiving?

Nog steeds is de projectie informatie die wordt meegeleverd met veel Open Source GIS tools niet de juiste. Althans het is wel de juiste volgens de epsg-database op basis waarvan de proj-definities worden gemaakt, maar in de praktijk werkt een andere definitie beter.

Gebruik van de standaard proj-definities heeft bijvoorbeeld tot gevolg dat er een verschuiving optreedt bij het laden van een shape-bestand welke je over GoogleMaps wilt draperen. Dit komt door het ‘herprojecteren’ van je data van epsg:28992 naar epsg:900913 (of epsg:3857 zoals die nu officieel heet).

Martijn van Exel schreef er al een blog over

Snelle oplossing: als je nog steeds een van de volgende proj definities ziet in qgis bij ‘Amersfoort / RD new’ (je mist dus het stukje ‘+towgs84′) :


+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defsa

of

+proj=sterea +lat_0=52.156161 +lon_0=5.387639 +k=0.999908 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs

Maak dan een ‘custom’ projectie definitie aan met een van de definitie:

+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.999908 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 +no_defs no_defs

Voorbeeld van een qgis project, met daarin GoogleMaps en een OpenStreetMap met daarin NIET de juiste definitie:

28992verschuiving1

En nu een met daarin de juiste proj definitie gebruikt:

Opmerking: als je dit nog nooit gezien hebt, en toch de ‘verkeerde’ definitie gebruikt, dan komt dit waarschijnlijk omdat je nog nooit gebruik hebt gemaakt van de ‘on the fly projection’ van Qgis in combinatie met een andere projectie.

4 thoughts on “Rijksdriehoekstelsel (epsg:28992) verschuiving?”

  1. Hi Tammo Jan,

    er zijn volgens mij twee dingen die je kunt doen:

    1) een ‘custom crs’ aanmaken die je je eigen code geeft. En die dan voor je data gebruiken

    2) of ‘harder’ is het aanpassen van de parameters in de sqlite database die QGIS gebruikt voor alle crs definities. Die zitten in een bestandje ‘qgis.db’ welke je met bijvoorbeeld de firefox extensie ‘sqlitemanager’ kunt openen en bewerken. Bij mij staat die in .qgis directory, op windows waarschijnlijk in je Settings directory in application data oid.
    Ik meen me zelfs te herinneren dat QGIS er zelf een heeft, en dat je een persoonlijke kopie hebt. Bij mij zie ik er bijvoorbeeld ook mijn eigen ‘ruimtelijke bookmarks’ in staan.

  2. Ik had dit probleem ook in R met een bestand (gem2014) van het CBS. Gelukkig zag ik hier op de site het antwoord van Richard.
    Heb het met de onderstaande code opgelost: (moest nog wel een + zetten voor laatste no_defs)
    proj4string(gem2014) <- "+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.999908 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 +no_defs +no_defs"

Geef een reactie

Je e-mailadres wordt niet gepubliceerd. Verplichte velden zijn gemarkeerd met *

De volgende HTML tags en attributen zijn toegestaan: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>