Verschuiving WGS1984 tov RD stelsel

classic Classic list List threaded Threaded
8 messages Options
Reply | Threaded
Open this post in threaded view
|

Verschuiving WGS1984 tov RD stelsel

Joris Meijerink-4
Ik heb van geofabrik.de de shapes van NL gedownload en in arcmap over een top25 kaart gezet.

De shapes van geofabrik zijn in WGS en de kaart in RD. Bij het importeren van de shapes geeft ArcMap aan dat hij de shapes moet omzetten van WGS naar RD. Na het omzetten vallen de top25 en de shapes echter niet samen. De shapes zijn ongeveer 130m naar het zuiden en 30m naar het westen verschoven. Bij het importeren kon ik zo geen nuttige opties vinden.

Iemand die een idee heeft waardoor dit komt?

gr Joris Meijerink


_______________________________________________
Talk-nl mailing list
[hidden email]
http://lists.openstreetmap.org/listinfo/talk-nl
Reply | Threaded
Open this post in threaded view
|

Re: Verschuiving WGS1984 tov RD stelsel

Stefan de Konink-3
On Tue, 10 Feb 2009, joris meijerink wrote:

> Iemand die een idee heeft waardoor dit komt?

Ik gok dat je de 'standaard' projectie hebt, en Milo heeft de 'echte'
projectie voor je ;)

Stefan


_______________________________________________
Talk-nl mailing list
[hidden email]
http://lists.openstreetmap.org/listinfo/talk-nl
Reply | Threaded
Open this post in threaded view
|

Re: Verschuiving WGS1984 tov RD stelsel

dbussche
In reply to this post by Joris Meijerink-4
Ik heb geen idee waarom dit zo is maar wij ervaren dit fenomeen ook.
Geen verschil of we de ESRI-projectie gebruiken of PROJ4.
Onze standaardoplossing sinds jaar en dag: een tooltje dat de verschuiving
automatisch ongedaan maakt en niet meer nadenken over het waarom...
Als er nog meer mensen dit ervaren kan ik dit tooltje (php/mapscript) ter
beschikking stellen om het op de OSM-server te zetten zodat iedereen eenvoudig
goede RD-shapes kan downloaden.
Maar eerst in de groep vragen of iemand het echte oorzaak kent, dat is dan
minder symptoombestrijding...



_______________________________________________
Talk-nl mailing list
[hidden email]
http://lists.openstreetmap.org/listinfo/talk-nl
Reply | Threaded
Open this post in threaded view
|

Re: Verschuiving WGS1984 tov RD stelsel

Stefan de Konink-3
Dirk Bussche wrote:
> Ik heb geen idee waarom dit zo is maar wij ervaren dit fenomeen ook.
> Geen verschil of we de ESRI-projectie gebruiken of PROJ4.

Geef je projectie string eens :)


Stefan

_______________________________________________
Talk-nl mailing list
[hidden email]
http://lists.openstreetmap.org/listinfo/talk-nl
Reply | Threaded
Open this post in threaded view
|

Re: Verschuiving WGS1984 tov RD stelsel

paul van den berg
In reply to this post by Joris Meijerink-4
joris meijerink <[hidden email]> wrote:

> Iemand die een idee heeft waardoor dit komt?

Ik heb hier ook een tijdje mee geworsteld.
Martijn van Excel heeft er een mooi verhaal over geschreven op
http://oegeo.wordpress.com/2008/05/20/note-to-self-the-one-and-only-rd-projection-string/

Zelf gebruik ik nu het volgende shell-script om 10x12.5 tiles in RD
te converteren naar wgs84 (met behulp van gdal en geotiff packages).

Gr, paulb


#!/bin/sh
# towgs.sh - convert 10x12.5 km RD-picture to wgs84 geotif
# $1 = left x
# $2 = lower y
# $3 = imagefile
# $4 = geotif image

# calculate RD corner coordinates for 10x12.5 images in meters
ulx=$(( $1 * 1000 ))
lry=$(( $2 * 1000 ))
uly=$(( $lry + 12500 ))
lrx=$(( $ulx + 10000 ))

# dump RD corner coordinates in temporary geotiff
gdal_translate -a_ullr $ulx $uly $lrx $lry $3 /tmp/1.tif

# gdal magic to translate RD to wgs84 in another geotif
gdalwarp -s_srs "+proj=sterea +lat_0=52.15616055555555\
 +lon_0=5.38763888888889 +k=0.9999079 +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"\
 -t_srs wgs84 /tmp/1.tif /tmp/2.tif

# produce high quality compressed jpeg inside resulting geotif image
geotifcp -c jpeg:80 /tmp/2.tif $4

# remove temporary files
rm /tmp/1.tif /tmp/2.tif


_______________________________________________
Talk-nl mailing list
[hidden email]
http://lists.openstreetmap.org/listinfo/talk-nl
Reply | Threaded
Open this post in threaded view
|

Re: Verschuiving WGS1984 tov RD stelsel

dbussche
In reply to this post by Stefan de Konink-3
Stefan, je hebt gelijk!
mijn projectie was
<28992> +proj=stere +lat_0=52.15616055555555 +lon_0=5.38763888888889
+k=0.999908 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs  <>
Dit is wat je vanzelf krijgt als je Proj4 installeert.
Hiermee ligt een projectie WGS84-RD ca. 100meter en een projectie Google
Mercator-RD ca. 300 meter ernaast.
De goede projectie is
<28992> +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 <>
Ik heb het net geprobeerd, met deze projectie lig je precies goed!
Ik ga dit op al onze computers aanpassen en ons verschuif-script weggooien.
Voor de meelezers met hetzelfde probleem: je moet dus het bestand
/usr/share/proj/epsg editen.
Achtergrondinfo
http://oegeo.wordpress.com/2008/05/20/
note-to-self-the-one-and-only-rd-projection-string/


_______________________________________________
Talk-nl mailing list
[hidden email]
http://lists.openstreetmap.org/listinfo/talk-nl
Reply | Threaded
Open this post in threaded view
|

Re: Verschuiving WGS1984 tov RD stelsel

Stefan de Konink-3
Dirk Bussche wrote:
> Stefan, je hebt gelijk!

Ik zal je quoten ;)


Stefan

_______________________________________________
Talk-nl mailing list
[hidden email]
http://lists.openstreetmap.org/listinfo/talk-nl
Reply | Threaded
Open this post in threaded view
|

Re: Verschuiving WGS1984 tov RD stelsel

osm-9
Stefan de Konink wrote:
> Dirk Bussche wrote:
>> Stefan, je hebt gelijk!
>
> Ik zal je quoten ;)
>
Heb je al vroeg dit jaar een keer gelijk gekregen. Goed werk :P

_______________________________________________
Talk-nl mailing list
[hidden email]
http://lists.openstreetmap.org/listinfo/talk-nl