[gdal-dev] GDAL 1.11 and EPSG:3857

[maw269]

Using GDAL 1.11

When I use the gdalwarp command:

gdalwarp -s_srs ESRI::EPG94.prj -t_srs epsg:*900913 *EPG94.ecw 900913.tif
I get legitimate gdalinfo with datum and spheroid defined.

Driver: GTiff/GeoTIFF
Files: 900913.tif
Size is 76318, 41566
Coordinate System is:
PROJCS["Google Maps Global 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"]]]
Origin = (13363025.433129203000000,-2674166.549594078200000)
Pixel Size = (0.163524580961778,-0.163524580961778)
Metadata:
  AREA_OR_POINT=Area
  COLORSPACE=RGB
  COMPRESSION_RATE_TARGET=9
  VERSION=2
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (13363025.433,-2674166.550) (120d 2'31.56"E, 23d29'18.24"S)
Lower Left  (13363025.433,-2680963.612) (120d 2'31.56"E, 23d32'40.94"S)
Upper Right (13375505.302,-2674166.550) (120d 9'15.15"E, 23d29'18.24"S)
Lower Right (13375505.302,-2680963.612) (120d 9'15.15"E, 23d32'40.94"S)
Center      (13369265.368,-2677565.081) (120d 5'53.36"E, 23d30'59.60"S)
Band 1 Block=76318x1 Type=Byte, ColorInterp=Red
  Description = Red
Band 2 Block=76318x1 Type=Byte, ColorInterp=Green
  Description = Green
Band 3 Block=76318x1 Type=Byte, ColorInterp=Blue
  Description = Blue
  
-----------------------------------------------------------------  
When I use the gdalwarp command:

gdalwarp -s_srs ESRI::EPG94.prj -t_srs epsg:*3857 *EPG94.ecw 3857.tif
I get illegitimate gdalinfo with datum and spheroid unknown and
unretrievable:

Driver: GTiff/GeoTIFF
Files: 3857.tif
Size is 76318, 41566
Coordinate System is:
LOCAL_CS["WGS 84 / Pseudo-Mercator",
    GEOGCS["WGS 84",
        *DATUM["unknown",*
            *SPHEROID["unretrievable - using*
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    AUTHORITY["EPSG","3857"],
    UNIT["metre",1]]
Origin = (13363025.433129203000000,-2674166
Pixel Size = (0.163524580961778,-0.16352458
Metadata:
  AREA_OR_POINT=Area
  COLORSPACE=RGB
  COMPRESSION_RATE_TARGET=9
  VERSION=2
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (13363025.433,-2674166.550)
Lower Left  (13363025.433,-2680963.612)
Upper Right (13375505.302,-2674166.550)
Lower Right (13375505.302,-2680963.612)
Center      (13369265.368,-2677565.081)
Band 1 Block=76318x1 Type=Byte, ColorInterp
  Description = Red
Band 2 Block=76318x1 Type=Byte, ColorInterp
  Description = Green
Band 3 Block=76318x1 Type=Byte, ColorInterp
  Description = Blue

 
This messes up my geoserver image pyramid when I try to use it.
---------------------------------

It seems epsg:900913 works but epsg:3857 does not.
The thing is, gdalwarp to 3857 worked in gdal 1.10 

Thanks,

Matt
 



--
View this message in context: http://osgeo-org.1560.x6.nabble.com/GDAL-1-11-and-EPSG-3857-tp5159753.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.