[FWTools] Changes in epsg file not considered

[Fischer, Andreas]

Dear list,
I'm working with FWTools 2.4.6 (Windows Vista) to transform raster data files. In order to transform from EPSG:31467 to EPSG:25832 I use gdalwarp in the following way:

	gdalwarp -s_srs "+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +units=m +nadgrids=C:\Programme\FWTools2.4.6\proj_lib\BETA2007.gsb +wktext" -t_srs "+proj=utm +zone=32 +ellps=GRS80 +datum=WGS84 +units=m +no_defs" -co TFW=YES inputfile.tif outputfile.tif

The result ist fine, except of the GeoTiff-Tags in the outputfile, since they point at WGS84 as the following listing shows:

	Coordinate System is:
	PROJCS["WGS 84 / UTM zone 32N",
	    GEOGCS["WGS 84",
	        DATUM["WGS_1984",
	            SPHEROID["WGS 84",6378137,298.2572235630016,
	                AUTHORITY["EPSG","7030"]],
	            AUTHORITY["EPSG","6326"]],
	        PRIMEM["Greenwich",0],
	        UNIT["degree",0.0174532925199433],
	        AUTHORITY["EPSG","4326"]],
	    PROJECTION["Transverse_Mercator"],
	    PARAMETER["latitude_of_origin",0],
	    PARAMETER["central_meridian",9],
	    PARAMETER["scale_factor",0.9996],
	    PARAMETER["false_easting",500000],
	    PARAMETER["false_northing",0],
	    UNIT["metre",1,
	        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32632"]]

But for EPSG:25832 it should like this:

	Coordinate System is:
	PROJCS["ETRS89 / UTM zone 32N",
	    GEOGCS["ETRS89",
	        DATUM["European_Terrestrial_Reference_System_1989",
	            SPHEROID["GRS 1980",6378137,298.2572221010042,
	                AUTHORITY["EPSG","7019"]],
	            AUTHORITY["EPSG","6258"]],
	        PRIMEM["Greenwich",0],
	        UNIT["degree",0.0174532925199433],
	        AUTHORITY["EPSG","4258"]],
	    PROJECTION["Transverse_Mercator"],
	    PARAMETER["latitude_of_origin",0],
	    PARAMETER["central_meridian",9],
	    PARAMETER["scale_factor",0.9996],
	    PARAMETER["false_easting",500000],
	    PARAMETER["false_northing",0],
	    UNIT["metre",1,
	        AUTHORITY["EPSG","9001"]],
	    AUTHORITY["EPSG","25832"]]

So I changed the epsg file in FWTools2.4.6/proj_lib to ad the missing '+datum=WGS84-Parameter' and changed the commandline of gdalwarp like this:

gdalwarp -s_srs "+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +units=m +nadgrids=C:\Programme\FWTools2.4.6\proj_lib\BETA2007.gsb +wktext" -t_srs EPSG:25832 -co TFW=YES inputfile.tif outputfile.tif

Now the GeoTiff-Tags are fine but the result isn't any more. :-( 
I assume that the change in the epsg-file is not considered, because I even deleted this file and gdalwarp didn't bother at all. I tried this on a linux-Server with gdal (1.4.2) and proj (4.5.0) too. In this case may changes in epsg were considered! Is this a bug in FWTools, uses gdalwarp other definitions than in epsg or did I make a mistake?

Thank you so far!


Mit freundlichen Grüßen
Im Auftrag

Andreas Fischer

Kreis Unna - Der Landrat 
Zentrale Datenverarbeitung
Friedrich-Ebert-Straße 17
59425 Unna 

Fon 02303 27-4416
Fax 02303 27-2896 
andreas.fischer at kreis-unna.de
www.kreis-unna.de

Diese E-Mail wurde beim Ausgang auf Viren geprueft. Wegen der potentiellen Gefahr auf den Uebertragungswegen wird zu einer Vireneingangskontrolle geraten. Eine Haftung für Virenfreiheit wird ausgeschlossen.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/fwtools/attachments/20091201/1c6e51d6/attachment.htm