joris meijerink <[email protected]> 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 [email protected] http://lists.openstreetmap.org/listinfo/talk-nl
