[gdal-dev] gdalwarp extent issues

Jay L. jzl5325 at psu.edu
Tue Feb 21 19:54:23 EST 2012


List,

I am attempting to reproject 8 gtiffs, each in a local equirectangular
projection, back to gcs for use in geoserver.  I immediately went to
gdalwarp to perform this task and am having an issue.  The input files are
stored in positive longitude, while the output SRS is in negative
longitude.  The "shift" to negative longitude is a primary reason for
wishing to perform the warp.

The input images which cover 0-90 north and 0-90 south are warping
correctly.  The others are not.

The other warped output images have an extent of -180 180, when the should
cover 90 degrees.  The images which remain in positive longitude are
drawing correctly, but have the remained of the extent filled in black.
 Meaning the image which spans 90-180n and 90-180s renders, but has an
extent of -180 to 180.  The input image which are shifted to negative
longitude have a global extent and are completely black.

I have tried the following without success:

gdalwarp -t_srs moon2000.prf input output_warped.tif

gdalwarp -s_srs <copy of the gdalinfo srs information in a prf> -t_srs
moon2000.prf input.tif output.tif

gdalwarp -t_srs moon2000.prf -te <desired extent as xmin ymin xmax ymax>
input.tif output.tif

Where moon2000.prf is:

GEOGCS["Moon 2000",
    DATUM["D_Moon_2000",
        SPHEROID["Moon_2000_IAU_IAG",1737400.0,0.0]],
    PRIMEM["Greenwich",0],
    UNIT["Decimal_Degree",0.0174532925199433]]

Any insight would be greatly appreciated.
Jay

Here is the projection info as output by gdalinfo:

Coordinate System is:
PROJCS["EQUIRECTANGULAR MOON",
    GEOGCS["GCS_MOON",
        DATUM["D_MOON",
            SPHEROID["MOON_localRadius",1737400,0]],
        PRIMEM["Reference_Meridian",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Equirectangular"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",0],
    PARAMETER["standard_parallel_1",0],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (100.000000000000000,1819400.000000000000000)
Pixel Size = (100.000000000000000,-100.000000000000000)

The coordinate information:
Corner Coordinates:
Upper Left  (     100.000, 1819400.000) (  0d 0'11.87"E, 59d59'59.88"N)
Lower Left  (     100.000,    -100.000) (  0d 0'11.87"E,  0d 0'11.87"S)
Upper Right ( 2729200.000, 1819400.000) ( 90d 0'11.69"E, 59d59'59.88"N)
Lower Right ( 2729200.000,    -100.000) ( 90d 0'11.69"E,  0d 0'11.87"S)
Center      ( 1364650.000,  909650.000) ( 45d 0'11.78"E, 29d59'54.00"N)

Oddly the central meridian and standard_parallel_1 parameters do not change
with each input file.  Since I know that gdal can see the s_srs from within
the files I have tried using:
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20120221/8c067867/attachment.html


More information about the gdal-dev mailing list