<div dir="ltr">J-S, from the looks of the coords, your original file may have already been in UTM to begin with. When you ran "gdal_translate -a_srs WGS84 file.tif file.wgs84.tif" you assigned a WGS84 projection to what appears to be a UTM tif file. Can you double-check this?<br>
<br>Glenn<br><br><div class="gmail_quote">On Mon, Sep 29, 2008 at 11:34 AM, Jean-Sébastien Guay <span dir="ltr"><<a href="mailto:jean-sebastien.guay@cm-labs.com">jean-sebastien.guay@cm-labs.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;">Hi all,<br>
<br>
This is probably just a sign of my inexperience in using gdal and GeoTiffs... I'm getting this error message when trying to convert a GeoTiff to UTM (in meters):<br>
<br>
========================================================================<br>
> gdalwarp -t_srs "+proj=utm +zone=32 +datum=WGS84 +units=m"<br>
  file.wgs84.tif file.utm.tif<br>
ERROR 1: Too many points (441 out of 441) failed to transform,<br>
unable to compute output bounds.<br>
========================================================================<br>
<br>
This is the info of the file:<br>
<br>
========================================================================<br>
> gdalinfo file.wgs84.tif<br>
Driver: GTiff/GeoTIFF<br>
Files: Gjoa-fledermaus.wgs84.tif<br>
Size is 6190, 7250<br>
Coordinate System is:<br>
GEOGCS["WGS 84",<br>
    DATUM["WGS_1984",<br>
        SPHEROID["WGS 84",6378137,298.2572235630016,<br>
            AUTHORITY["EPSG","7030"]],<br>
        AUTHORITY["EPSG","6326"]],<br>
    PRIMEM["Greenwich",0],<br>
    UNIT["degree",0.0174532925199433],<br>
    AUTHORITY["EPSG","4326"]]<br>
Origin = (548676.250000000000000,6807000.000000000000000)<br>
Pixel Size = (0.955613893376414,-0.956137931034483)<br>
Metadata:<br>
  AREA_OR_POINT=Area<br>
Image Structure Metadata:<br>
  INTERLEAVE=PIXEL<br>
Corner Coordinates:<br>
Upper Left  (  548676.250, 6807000.000)<br>
            (548676d15'0.00"E,6807000d 0'25769803776.00"N)<br>
Lower Left  (  548676.250, 6800068.000)<br>
            (548676d15'0.00"E,6800068d 0'25769803776.00"N)<br>
Upper Right (  554591.500, 6807000.000)<br>
            (554591d30'0.00"E,6807000d 0'25769803776.00"N)<br>
Lower Right (  554591.500, 6800068.000)<br>
            (554591d30'0.00"E,6800068d 0'25769803776.00"N)<br>
Center      (  551633.875, 6803534.000)<br>
            (551633d52'30.00"E,6803534d 0'25769803776.00"N)<br>
Band 1 Block=6190x1 Type=Byte, ColorInterp=Red<br>
  Mask Flags: PER_DATASET ALPHA<br>
Band 2 Block=6190x1 Type=Byte, ColorInterp=Green<br>
  Mask Flags: PER_DATASET ALPHA<br>
Band 3 Block=6190x1 Type=Byte, ColorInterp=Blue<br>
  Mask Flags: PER_DATASET ALPHA<br>
Band 4 Block=6190x1 Type=Byte, ColorInterp=Alpha<br>
========================================================================<br>
<br>
The first thing I noticed is that what is in the SPHEROID section doesn't seem to match the other coordinates in the file, but I'm not sure that's the problem, and even if my hunch is correct, I don't know how to fix it.<br>

<br>
For reference, the original file (which I got from importing point data into some third party software and then exporting a GeoTiff from it) had this info:<br>
<br>
========================================================================<br>
>gdalinfo file.tif<br>
Driver: GTiff/GeoTIFF<br>
Files: Gjoa-fledermaus.tif<br>
Size is 6190, 7250<br>
Coordinate System is:<br>
PROJCS["unnamed",<br>
    GEOGCS["unnamed",<br>
        DATUM["unknown",<br>
            SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],<br>
        PRIMEM["Greenwich",0],<br>
        UNIT[,0.0174532925199433]],<br>
    UNIT["unknown",1]]<br>
Origin = (548676.250000000000000,6807000.000000000000000)<br>
Pixel Size = (0.955613893376414,-0.956137931034483)<br>
Image Structure Metadata:<br>
  COMPRESSION=LZW<br>
  INTERLEAVE=PIXEL<br>
Corner Coordinates:<br>
Upper Left  (  548676.250, 6807000.000)<br>
Lower Left  (  548676.250, 6800068.000)<br>
Upper Right (  554591.500, 6807000.000)<br>
Lower Right (  554591.500, 6800068.000)<br>
Center      (  551633.875, 6803534.000)<br>
Band 1 Block=6190x1 Type=Byte, ColorInterp=Red<br>
  Mask Flags: PER_DATASET ALPHA<br>
Band 2 Block=6190x1 Type=Byte, ColorInterp=Green<br>
  Mask Flags: PER_DATASET ALPHA<br>
Band 3 Block=6190x1 Type=Byte, ColorInterp=Blue<br>
  Mask Flags: PER_DATASET ALPHA<br>
Band 4 Block=6190x1 Type=Byte, ColorInterp=Alpha<br>
========================================================================<br>
<br>
and I got file.wgs84.tif by running<br>
<br>
========================================================================<br>
> gdal_translate -a_srs WGS84 file.tif file.wgs84.tif<br>
========================================================================<br>
<br>
on it.<br>
<br>
Any hints would be appreciated. Thanks in advance,<br>
<br>
J-S<br>
-- <br>
______________________________________________________<br>
Jean-Sebastien Guay    <a href="mailto:jean-sebastien.guay@cm-labs.com" target="_blank">jean-sebastien.guay@cm-labs.com</a><br>
                               <a href="http://www.cm-labs.com/" target="_blank">http://www.cm-labs.com/</a><br>
                        <a href="http://whitestar02.webhop.org/" target="_blank">http://whitestar02.webhop.org/</a><br>
_______________________________________________<br>
osg-users mailing list<br>
<a href="mailto:osg-users@lists.openscenegraph.org" target="_blank">osg-users@lists.openscenegraph.org</a><br>
<a href="http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org" target="_blank">http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Glenn Waldron : Pelican Mapping : <a href="http://pelicanmapping.com">http://pelicanmapping.com</a> : +1.703.652.4791<br>
</div>