Hello,<br><br>I'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"><<a href="mailto:ragi@burhum.com">ragi@burhum.com</a>></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["ED50 / UTM zone 31N",<br>
GEOGCS["ED50",<br> DATUM["European_Datum_1950",<br> SPHEROID["International 1924",6378388,297.0000000000014,<br> AUTHORITY["EPSG","7022"]],<br>
AUTHORITY["EPSG","6230"]],<br> PRIMEM["Greenwich",0],<br> UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4230"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",0],<br> PARAMETER["central_meridian",3],<br> PARAMETER["scale_factor",0.9996],<br> PARAMETER["false_easting",500000],<br> PARAMETER["false_northing",0],<br>
UNIT["metre",1,<br> AUTHORITY["EPSG","9001"]],<br> AUTHORITY["EPSG","23031"]]<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'46.94"E, 42d42'34.10"N)<br>Lower Left ( 327720.000, 4728383.500) ( 0d53'49.56"E, 42d41'16.35"N)<br>
Upper Right ( 331193.000, 4730783.500) ( 0d56'19.50"E, 42d42'36.88"N)<br>Lower Right ( 331193.000, 4728383.500) ( 0d56'22.08"E, 42d41'19.12"N)<br>Center ( 329456.500, 4729583.500) ( 0d55'4.52"E, 42d41'56.62"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 "+proj=utm +zone=31 +ellps=intl +units=m +towgs84=-87,-101,-121 +no_defs" -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>