[gdal-dev] issues with gdal2tiles.py

[Marcela Campo]

Hi,
  I am trying to use gdal2tiles with a set of geotiff and I get the
following error

ERROR 6: EPSG PCS/GCS code 900913 not found in EPSG support files.  Is this
a valid
EPSG coordinate system?
ERROR 6: No translation for an empty SRS to PROJ.4 format is known.


I have already transformed to projection EPSG:3857, here is the gdalinfo
output for one of the files

Driver: GTiff/GeoTIFF
Files: rec_230076_1994-12-31_f4_p.tif
       rec_230076_1994-12-31_f4_p.tfw
Size is 728, 741
Coordinate System is:
PROJCS["WGS 84 / Pseudo-Mercator",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
+lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext
 +no_defs"],
    AUTHORITY["EPSG","3857"]]
Origin = (-7130153.280253740958869,-2651286.133427014108747)
Pixel Size = (32.745668091799999,-32.745668092300001)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_IMAGEDESCRIPTION=ImageJ=1.48v

Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-7130153.280,-2651286.133) ( 64d 3' 4.52"W, 23d 9'33.23"S)
Lower Left  (-7130153.280,-2675550.673) ( 64d 3' 4.52"W, 23d21'34.15"S)
Upper Right (-7106314.434,-2651286.133) ( 63d50'13.59"W, 23d 9'33.23"S)
Lower Right (-7106314.434,-2675550.673) ( 63d50'13.59"W, 23d21'34.15"S)
Center      (-7118233.857,-2663418.403) ( 63d56'39.06"W, 23d15'33.82"S)
Band 1 Block=728x3 Type=Byte, ColorInterp=Red
Band 2 Block=728x3 Type=Byte, ColorInterp=Green
Band 3 Block=728x3 Type=Byte, ColorInterp=Blue


the command generates the tiles, but as the images appear a bit shifted I
wonder if it somewhat misinterpreted the projection.

I am using gdal 1.10.1 on ubuntu.

Thanks!
Marcela
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20140527/88c34442/attachment.html>