[gdal-dev] gdalwarp: utm to polarsterographic projection

[Anja Roesel]

Dear list,

i tried several options to reproject coordinates from utm to polar
stereographic but without success. It seems that I use some wrong
parameters.
Could you tell me what would be the best command to project it correctly?

These are the information from gdalinfo for the sourcefile:

> >gdalinfo utm.tif
Driver: GTiff/GeoTIFF
Size is 9101, 8561
Coordinate System is:
PROJCS["WGS 84 / UTM zone 10N",
    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["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-123],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32610"]]
Origin = (279000.000000,7853400.000000)
Pixel Size = (30.00000000,-30.00000000)
Metadata:
  AREA_OR_POINT=Point
Corner Coordinates:
Upper Left  (  279000.000, 7853400.000) (128d59'44.63"W, 70d41'12.77"N)
Lower Left  (  279000.000, 7596570.000) (128d23'2.66"W, 68d23'42.08"N)
Upper Right (  552030.000, 7853400.000) (121d35'1.70"W, 70d46'45.81"N)
Lower Right (  552030.000, 7596570.000) (121d43'44.87"W, 68d28'36.75"N)
Center      (  415515.000, 7724985.000) (125d10'26.99"W, 69d37'11.51"N)
Band 1 Block=9101x1 Type=Byte, ColorInterp=Gray


This is one of my commands:

gdalwarp -s_srs "+proj=utm  +zone=10 +datum=WGS84" -t_srs "+proj=stere
+lat_ts=70.0 +lat_0=90.0 +lon_0=45.0 +k_0=1.0 +x_0=0.0 +y_0=0.0" utm.tif
pol_ster.tif


Thank you very much for your comments and remarks.


-- Anja Rösel