[Gdal-dev] Converting Polar Stereographic to WGS84 long/lat

[Ethan Alpert]

I needed to set the -te parameter for this to work

 

 

 

________________________________

From: gdal-dev-bounces at lists.maptools.org
[mailto:gdal-dev-bounces at lists.maptools.org] On Behalf Of Ethan Alpert
Sent: Wednesday, February 07, 2007 12:30 PM
To: gdal-dev at lists.maptools.org
Subject: [Gdal-dev] Converting Polar Stereographic to WGS84 long/lat

 

I'm use GDAL 1.4.0

 

I thought the following gdalwarp invocation would do the trick but it
creates a black image positioned completely wrong. My file N.tif covers
the northern hemisphere. The gdalinfo is below as well as the gdalinfo
of the output file. Do I need to do something else to convert this data?

 

gdalwarp -t_srs "EPSG:4326"  N.tif WGS.tif

 

 

-bash-2.05b$ gdalinfo N.tif

Driver: GTiff/GeoTIFF

Size is 1024, 1024

Coordinate System is:

PROJCS["unnamed",

    GEOGCS["Coordinate System imported from GRIB file",

        DATUM["unknown",

            SPHEROID["unnamed",6367470,0]],

        PRIMEM["Greenwich",0],

        UNIT["degree",0.0174532925199433]],

    PROJECTION["Polar_Stereographic"],

    PARAMETER["latitude_of_origin",60],

    PARAMETER["central_meridian",-80],

    PARAMETER["scale_factor",90],

    PARAMETER["false_easting",0],

    PARAMETER["false_northing",0],

    UNIT["metre",1,

        AUTHORITY["EPSG","9001"]]]

Origin = (-12184935.659713378176093,12184935.659713353961706)

Pixel Size = (23813.000000000000000,-23813.000000000000000)

Metadata:

  AREA_OR_POINT=Area

Corner Coordinates:

Upper Left  (-12184935.660,12184935.660) (145d 0'0.00"E, 20d49'33.60"S)

Lower Left  (-12184935.660,-12199576.340) (124d57'56.16"W,
20d51'29.41"S)

Upper Right (12199576.340,12184935.660) ( 54d57'56.16"E, 20d51'29.41"S)

Lower Right (12199576.340,-12199576.340) ( 35d 0'0.00"W, 20d53'25.05"S)

Center      (    7320.340,   -7320.340) ( 35d 0'0.00"W, 89d54'0.57"N

 

-bash-2.05b$ gdalinfo WGS.tif

Driver: GTiff/GeoTIFF

Size is 2851, 157

Coordinate System is:

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"]]

Origin = (-175.679918813105417,-1.452730998140064)

Pixel Size = (0.124296121889026,-0.124296121889026)

Metadata:

  AREA_OR_POINT=Area

Corner Coordinates:

Upper Left  (-175.6799188,  -1.4527310) (175d40'47.71"W,  1d27'9.83"S)

Lower Left  (-175.6799188, -20.9672221) (175d40'47.71"W, 20d58'2.00"S)

Upper Right ( 178.6883247,  -1.4527310) (178d41'17.97"E,  1d27'9.83"S)

Lower Right ( 178.6883247, -20.9672221) (178d41'17.97"E, 20d58'2.00"S)

Center      (   1.5042029, -11.2099766) (  1d30'15.13"E, 11d12'35.92"S)

Band 1 Block=2851x2 Type=Byte, ColorInterp=Gray

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20070207/b872dd78/attachment.html