Hello,<br><br>I&#39;m from Spain, and I can confirm your parameters. These worked for me as well.<br><br>Regards<br>Jorge<br><br><br><div class="gmail_quote">On Tue, May 12, 2009 at 1:03 AM, Ragi Y. Burhum <span dir="ltr">&lt;<a href="mailto:ragi@burhum.com">ragi@burhum.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">I was just sending this e-mail so that it can be searched through the archives. Thanks to Even Rouault for walking me through something that should have been obvious :)<br>

<br>Originally, I was having trouble finding the right parameters for gdalwarp to be able to warp some imagery that I got from Spain. These were in ED50 / UTM zone 31N:<br>
<br>rburhum@cuzconix:~/p_bug$ gdalinfo exported_from_original_sid_file.tif <br>Driver: GTiff/GeoTIFF<br>Files: exported_from_original_sid_file.tif<br>Size is 6946, 4800<br>Coordinate System is:<br>PROJCS[&quot;ED50 / UTM zone 31N&quot;,<br>


    GEOGCS[&quot;ED50&quot;,<br>        DATUM[&quot;European_Datum_1950&quot;,<br>            SPHEROID[&quot;International 1924&quot;,6378388,297.0000000000014,<br>                AUTHORITY[&quot;EPSG&quot;,&quot;7022&quot;]],<br>


            AUTHORITY[&quot;EPSG&quot;,&quot;6230&quot;]],<br>        PRIMEM[&quot;Greenwich&quot;,0],<br>        UNIT[&quot;degree&quot;,0.0174532925199433],<br>        AUTHORITY[&quot;EPSG&quot;,&quot;4230&quot;]],<br>

    PROJECTION[&quot;Transverse_Mercator&quot;],<br>
    PARAMETER[&quot;latitude_of_origin&quot;,0],<br>    PARAMETER[&quot;central_meridian&quot;,3],<br>    PARAMETER[&quot;scale_factor&quot;,0.9996],<br>    PARAMETER[&quot;false_easting&quot;,500000],<br>    PARAMETER[&quot;false_northing&quot;,0],<br>


    UNIT[&quot;metre&quot;,1,<br>        AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],<br>    AUTHORITY[&quot;EPSG&quot;,&quot;23031&quot;]]<br>Origin = (327720.000000000000000,4730783.500000000000000)<br>Pixel Size = (0.500000000000000,-0.500000000000000)<br>


Metadata:<br>  AREA_OR_POINT=Area<br>Image Structure Metadata:<br>  INTERLEAVE=PIXEL<br>Corner Coordinates:<br>Upper Left  (  327720.000, 4730783.500) (  0d53&#39;46.94&quot;E, 42d42&#39;34.10&quot;N)<br>Lower Left  (  327720.000, 4728383.500) (  0d53&#39;49.56&quot;E, 42d41&#39;16.35&quot;N)<br>


Upper Right (  331193.000, 4730783.500) (  0d56&#39;19.50&quot;E, 42d42&#39;36.88&quot;N)<br>Lower Right (  331193.000, 4728383.500) (  0d56&#39;22.08&quot;E, 42d41&#39;19.12&quot;N)<br>Center      (  329456.500, 4729583.500) (  0d55&#39;4.52&quot;E, 42d41&#39;56.62&quot;N)<br>


Band 1 Block=6946x1 Type=Byte, ColorInterp=Red<br>Band 2 Block=6946x1 Type=Byte, ColorInterp=Green<br>Band 3 Block=6946x1 Type=Byte, ColorInterp=Blue<br><br><br>In order to get it to work correctly, I had to use the right +towgs84 parameters from <a href="http://earth-info.nga.mil/GandG/coordsys/onlinedatum/CountryEuropeTable.html" target="_blank">http://earth-info.nga.mil/GandG/coordsys/onlinedatum/CountryEuropeTable.html</a> . Although the table showed  -84 (+/-5), -107(+/-6), -120 (+/-3), for my own purposes -87, -98, -121 seemed to work best based on my source data.<br>


<br>so the command line would be:<br><br> gdalwarp -s_srs &quot;+proj=utm +zone=31 +ellps=intl +units=m +towgs84=-87,-101,-121 +no_defs&quot; -t_srs EPSG:4326 exported_from_original_sid_file.tif testout.tif<br><br>Hope this helps someone else,<br>

<font color="#888888">
<br>- Ragi<br><br><br><br><br><br>
</font><br>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></blockquote></div><br>