<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7652.24">
<TITLE>Changes in epsg file not considered</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->
<P><FONT FACE="Arial">Dear list,</FONT>
<BR><FONT FACE="Arial">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:</FONT></P>
<UL>
<P><FONT FACE="Arial">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</FONT></P>
</UL>
<P><FONT FACE="Arial">The result ist fine, except of the GeoTiff-Tags in the outputfile, since they point at WGS84 as the following listing shows:</FONT></P>
<UL>
<P><FONT FACE="Courier New">Coordinate System is:</FONT>
<BR><FONT FACE="Courier New">PROJCS["WGS 84 / UTM zone 32N",</FONT>
<BR><FONT FACE="Courier New"> GEOGCS["WGS 84",</FONT>
<BR><FONT FACE="Courier New"> DATUM["WGS_1984",</FONT>
<BR><FONT FACE="Courier New"> SPHEROID["WGS 84",6378137,298.2572235630016,</FONT>
<BR><FONT FACE="Courier New"> <SPAN LANG="en-us"> AUTHORITY["EPSG","7030"]],</SPAN></FONT>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> AUTHORITY["EPSG","6326"]],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> PRIMEM["Greenwich",0],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> UNIT["degree",0.0174532925199433],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> AUTHORITY["EPSG","4326"]],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> PROJECTION["Transverse_Mercator"],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> PARAMETER["latitude_of_origin",0],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> PARAMETER["central_meridian",9],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> PARAMETER["scale_factor",0.9996],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> PARAMETER["false_easting",500000],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> PARAMETER["false_northing",0],</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> UNIT["metre",1,</FONT></SPAN>
<BR><SPAN LANG="en-us"><FONT FACE="Courier New"> </FONT></SPAN><SPAN LANG="de"> <FONT FACE="Courier New">AUTHORITY["EPSG","9001"]],</FONT></SPAN>
</UL>
<P><SPAN LANG="de"><FONT FACE="Courier New"> AUTHORITY["EPSG","32632"]]</FONT></SPAN>
</P>
<P><SPAN LANG="de"><FONT FACE="Arial">But for EPSG:25832 it should like this:</FONT></SPAN>
</P>
<UL>
<P><SPAN LANG="de"><FONT FACE="Courier New">Coordinate System is:</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New">PROJCS["ETRS89 / UTM zone 32N",</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> GEOGCS["ETRS89",</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> DATUM["European_Terrestrial_Reference_System_1989",</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> SPHEROID["GRS 1980",6378137,298.2572221010042,</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> AUTHORITY["EPSG","7019"]],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> AUTHORITY["EPSG","6258"]],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> PRIMEM["Greenwich",0],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> UNIT["degree",0.0174532925199433],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> AUTHORITY["EPSG","4258"]],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> PROJECTION["Transverse_Mercator"],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> PARAMETER["latitude_of_origin",0],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> PARAMETER["central_meridian",9],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> PARAMETER["scale_factor",0.9996],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> PARAMETER["false_easting",500000],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> PARAMETER["false_northing",0],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> UNIT["metre",1,</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> AUTHORITY["EPSG","9001"]],</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Courier New"> AUTHORITY["EPSG","25832"]]</FONT></SPAN>
</P>
</UL>
<P><SPAN LANG="de"><FONT FACE="Arial">So I changed the epsg file in FWTools2.4.6/proj_lib to ad the missing '</FONT><FONT FACE="Arial">+datum=WGS84-Parameter' and changed the commandline of gdalwarp like this:</FONT></SPAN></P>
<P><SPAN LANG="de"><FONT FACE="Arial">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"<B> -t_srs EPSG:25832</B> -co TFW=YES inputfile.tif outputfile.tif</FONT></SPAN></P>
<P><SPAN LANG="de"><FONT FACE="Arial">Now the GeoTiff-Tags are fine but the result isn't any more. :-( </FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Arial">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?</FONT></SPAN></P>
<P><SPAN LANG="de"><FONT FACE="Arial">Thank you so far!</FONT></SPAN>
</P>
<BR>
<P><SPAN LANG="de"><FONT FACE="Arial">Mit freundlichen Grüßen</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Arial">Im Auftrag</FONT></SPAN>
</P>
<P><SPAN LANG="de"><FONT FACE="Arial">Andreas Fischer</FONT></SPAN>
</P>
<P><SPAN LANG="de"><FONT FACE="Arial">Kreis Unna - Der Landrat<BR>
Zentrale Datenverarbeitung<BR>
Friedrich-Ebert-Straße 17<BR>
59425 Unna </FONT></SPAN>
</P>
<P><SPAN LANG="fr"><FONT FACE="Arial">Fon 02303 27-4416<BR>
Fax 02303 27-2896</FONT></SPAN><SPAN LANG="de"></SPAN><SPAN LANG="de"><BR>
<FONT FACE="Arial">andreas.fischer@kreis-unna.de</FONT></SPAN>
<BR><SPAN LANG="de"><FONT FACE="Arial">www.kreis-unna.de</FONT></SPAN>
</P>
<pre>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.
</pre></BODY>
</HTML>