[mapserver-users] reprojection

[Julien Cigar]

  Hello,

I'm trying to re-project a RASTER layer but without success so far.
I have the following TIFF file:

mage at mordor:/home/mage/gis % gdalinfo antarctica_sthn_ocean.tiff
Driver: GTiff/GeoTIFF
Files: antarctica_sthn_ocean.tiff
Size is 12216, 9094
Coordinate System is:
PROJCS[" Projection Name = WGS_1984_Stereographic_South_Pole Units = 
meters GeoTIFF Units = meters",
     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["Polar_Stereographic"],
     PARAMETER["latitude_of_origin",-71],
     PARAMETER["central_meridian",0],
     PARAMETER["scale_factor",1],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
Origin = (-8563840.575678950175643,6466755.550653258338571)
Pixel Size = (1420.390210742724321,-1420.394765779645695)
Metadata:
   TIFFTAG_XRESOLUTION=900
   TIFFTAG_YRESOLUTION=900
   TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-8563840.576, 6466755.551) ( 52d56'33.71"W,  8d 9'27.23"S)
Lower Left  (-8563840.576,-6450314.449) (126d59'13.87"W,  8d12'36.90"S)
Upper Right ( 8787646.239, 6466755.551) ( 53d39'3.45"E,  7d12'35.07"S)
Lower Right ( 8787646.239,-6450314.449) (126d16'45.99"E,  7d15'38.97"S)
Center      (  111902.832,    8220.551) ( 85d47'54.65"E, 88d58'2.41"S)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
   Overviews: 6108x4547, 3054x2274, 1527x1137, 764x569
   Metadata:
     LAYER_TYPE=athematic
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
   Overviews: 6108x4547, 3054x2274, 1527x1137, 764x569
   Metadata:
     LAYER_TYPE=athematic
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
   Overviews: 6108x4547, 3054x2274, 1527x1137, 764x569
   Metadata:
     LAYER_TYPE=athematic

If I understood well the output of gdalinfo, the projection seems to be 
EPSG:4326.
Now I want to re-project this to EPSG:3031.
I tried with the following MAP file :

MAP
     NAME "antarctic"
     STATUS ON
     SIZE 600 400
     EXTENT -8563840.576 -6450314.449 8787646.239 6466755.551
     IMAGECOLOR 255 255 255

     PROJECTION
         "init=epsg:3031"
     END

     OUTPUTFORMAT
         NAME png
         DRIVER "GD/PNG"
         MIMETYPE "image/png"
         IMAGEMODE RGBA
         EXTENSION "png"
     END

     LAYER
         NAME "ocean"
         TYPE RASTER
         STATUS DEFAULT
         DATA antarctica_sthn_ocean.tiff
         PROJECTION
             "init=epsg:4326"
         END
     END
END

... but all what I get is a gray 600x400 PNG image ...

Any idea what could be the problem ?

Thanks,
Julien
-------------- next part --------------
A non-text attachment was scrubbed...
Name: jcigar.vcf
Type: text/x-vcard
Size: 254 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/mapserver-users/attachments/20100921/5003182b/jcigar.vcf