[gdal-dev] help with gdalwarp

[Alex Bernstein]

Hi all,

Sorry for a newbie question. I've been reading this list and searching
elsewhere for a few days now to try to figure this out, but I'm still
missing something. I need to reproject a map from Plate Caree
projection to Lambert Azimuthal Equal Area projection centered on
0N,15E with gdalwarp. I'm using a Natural Earth II map from
<http://www.shadedrelief.com/natural2/>. The map
<http://www.nacis.org/data/NE2/NE2_land_only.jpg> is 16200x8100
pixels, plate caree, WGS84. It also has a world file
<http://www.shadedrelief.com/natural2/images/world_16000.zip> with
following information.
          0.02222222222222
          0.00000000000000
          0.00000000000000
          -0.02222222222222
          -180.00000000000000
          90.00000000000000

I run gdalwarp with following parameters:
gdalwarp -s_srs "+proj=eqc +lat_ts=0 +lon_0=0 +datum=WGS84" -t_srs
"+proj=laea +lat_0=0 +lon_0=15 +x_0=0 +y_0=0 +datum=WGS84" -co
INTERLEAVE=PIXEL -co TFW=YES -of GTiff
c:\dgwarp\source_images\NE2_land_only.jpg
c:\dgwarp\NE2_land_only_laea.tif

But the resulting image looks basically the same as the input,
definitely not like LAEA should look. gdalinfo on the output image
shows strange corner coordinates:
Driver: GTiff/GeoTIFF
Files: NE2_land_only_laea.tif
Size is 16166, 8168
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235630016,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Lambert_Azimuthal_Equal_Area"],
    PARAMETER["latitude_of_center",0],
    PARAMETER["longitude_of_center",15],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-1665206.343060958000000,90.180212867784576)
Pixel Size = (0.022078666990796,-0.022078666990796)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-1665206.343,      90.180) (  0d 0'5.82"W,  0d 0'2.91"N)
Lower Left  (-1665206.343,     -90.158) (  0d 0'5.82"W,  0d 0'2.91"S)
Upper Right (-1664849.419,      90.180) (  0d 0'5.82"E,  0d 0'2.91"N)
Lower Right (-1664849.419,     -90.158) (  0d 0'5.82"E,  0d 0'2.91"S)
Center      (-1665027.881,       0.011) (  0d 0'0.00"W,  0d 0'0.00"N)
Band 1 Block=16166x1 Type=Byte, ColorInterp=Red
Band 2 Block=16166x1 Type=Byte, ColorInterp=Green
Band 3 Block=16166x1 Type=Byte, ColorInterp=Blue

What am I missing?