[Gdal-dev] sterea projection

[Jan Hartmann]

Hi,

I am trying to warp Geotiffs from epsg:28992 (Dutch national grid) to 
EPSG4326 (lat-lon WGS84). The PROJ epsg-file gives "stere" as projection 
for 28992, but "sterea" gives a better approximation in this case (see 
http://bugzilla.remotesensing.org/show_bug.cgi?id=980). Also, a 
"towgs84" parameter is needed for transforming to WGS84. I can use the 
complete, correct PROJ-string in cs2cs, mapserver and postgis by 
replacing the values in the epsg file/spatial_ref_sys table, but for 
gdalwarp this does not work. When I give the complete PROJ-string in 
gdalwarp, e.g.

gdalwarp -s_srs "+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 +to +proj=latlong +datum=WGS84" -t_srs input.tif output.tif

I get an error "Error 1: Translating source or target SRS failed." With 
"+proj=stere" everything works, but then the transformation is not exact.

Is there any way to get these parameters into gdalwarp?

Jan Hartmann
Department of Geography
University of Amsterdam