[OSM-talk-nl] Verschuiving WGS1984 tov RD stelsel

[paul van den berg]

joris meijerink <meijerink at gmail.com> 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