[Gdal-dev] gdalwarp reprojection mercator to polar stereographic question

[Dori]

I am trying to reproject a Southern Hemisphere GeoTiff satellite picture
from Mercator to Polar Stereographic.  I want to obtain an image similar to
the one in : http://aviationweather.gov/obs/sat/intl/ for ICAO area J or K.
Pls. see input image info at the end of the message.
I tried following options to just reproject the map with not so promising
outputs:
gdalwarp -t_srs EPSG:32761 -wo SAMPLE_GRID=YES in.tif out.tif  => original
image diagonally on a white background but with southern latitudes on top of
northern latitudes like //
gdalwarp -t_srs "+proj=stere" in.tif out.tif => fully white image
gdalwarp -t_srs "+proj=stere +lat_ts=0.0 +lat+0=-90.0 +lon_0=-179.9999 +k_0=
1.0 +x_0=0 +y_0=0" -wo SAMPLE_GRID=YES in.tif out.tif => original images
diagonally shown on a white background like \\

I tried some other combinations with same bad results.  I have already
looked into the list of questions posted already with no much luck either.

Does someone know which are the right arguments for gdalwarp to obtain the
image I want to obtain?


==========================================================================================
The input image info:
Driver: GTiff/GeoTIFF
Size is 7482, 1871
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235629972,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Equirectangular"],
    PARAMETER["latitude_of_origin",30],
    PARAMETER["central_meridian",142],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-179.999900,0.029583)
Pixel Size = (0.04811252,-0.04811252)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  (-179.9999000,   0.0295831) (141d59'53.28"E,  0d 0'0.00"N)
Lower Left  (-179.9999000, -89.9889464) (141d59'53.28"E,  0d 0'2.93"S)
Upper Right ( 179.9779928,   0.0295831) (142d 0'6.72"E,  0d 0'0.00"N)
Lower Right ( 179.9779928, -89.9889464) (142d 0'6.72"E,  0d 0'2.93"S)
Center      (  -0.0109536, -44.9796817) (142d 0'0.00"E,  0d 0'1.46"S)
Band 1 Block=7482x1 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
    0: 255,255,255,255
    1: 0,0,0,255
    2: 0,0,85,255
    3: 0,0,170,255
    4: 0,0,255,255
    5: 0,85,0,255
    6: 0,85,85,255
    7: 0,85,170,255
.......... rest of color table
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20070205/28915841/attachment.html