[FWTools] FWtools 2.4.0 (gdalwarp - ERROR 1: Too many points 441 ...)

[Trent M Hare]

I've attempted to track down this issue but am not getting anywhere. Using 
FWTools 2.1.0 (GDAL 1.6.0dev), I was able to project from Equirectangular 
(aka Equidistant_Cylindrical) to Mercator (or other projections) using 
gdalwarp.  Switching to FWTool 2.4.0 (GDAL 1.7.0dev), and using the exact 
same gdalwarp command I now get:

>gdalwarp -t_srs "ESRI::Mars2000_Merc_clon180.prj" ESP_011909_1320_RED.JP2 
ESP_011909_1320_RED_merc180.tif
ERROR 1: Too many points (441 out of 441) failed to transform, unable to 
compute output bounds.

My suspicion is the source projection Equirectangular is being handled 
incorrectly  - potentially from the "latitude_of_origin" / 
"standard_parallel_1" keyword issue. That is the only difference that I 
can see. Any ideas? I know this is Mars but that is not the issue. This 
might be more suited to the GDAL list but I thought there might be some 
gdalwarp flags I can pass to help resolve the issue.

Source image used for testing (geoJpeg2000).
http://hirise-pds.lpl.arizona.edu/PDS/RDR/ESP/ORB_011900_011999/ESP_011909_1320/ESP_011909_1320_RED.JP2 
(272MB)
parent web site: http://hirise.lpl.arizona.edu/ESP_011909_1320

Contents of Mars2000_Merc_clon180.prj
PROJCS["Mars2000_Merc_clon180",GEOGCS["GCS_Mars_2000_Sphere",DATUM["D_Mars_2000_Sphere",SPHEROID["Mars_2000_Sphere_IAU_IAG",3396190.0,0.0]],PRIMEM["Reference_Meridian",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",180.0],PARAMETER["Standard_Parallel_1",0.0],UNIT["Meter",1.0]]

Thanks for any help or ideas.

Regards,
Trent




>gdalinfo ESP_011909_1320_RED.JP2
Driver: JP2ECW/ERMapper JPEG2000
Files: ESP_011909_1320_RED.JP2
       ESP_011909_1320_RED.J2w
       ESP_011909_1320_RED.JP2.aux.xml
Size is 15167, 36366
Coordinate System is:
PROJCS["Equirectangular MARS",
    GEOGCS[,
        DATUM["unknown",
            SPHEROID["unnamed",3386150.7470034,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Equirectangular"],
    PARAMETER["latitude_of_origin",-45],
    PARAMETER["central_meridian",180],
    PARAMETER["standard_parallel_1",0],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-6243330.250000000000000,-2815019.750000000000000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Corner Coordinates:
Upper Left  (-6243330.250,-2815019.750) ( 74d21'32.27"E, 92d37'54.79"S)
Lower Left  (-6243330.250,-2833202.750) ( 74d21'32.27"E, 92d56'22.40"S)
Upper Right (-6235746.750,-2815019.750) ( 74d29'14.21"E, 92d37'54.79"S)
Lower Right (-6235746.750,-2833202.750) ( 74d29'14.21"E, 92d56'22.40"S)
Center      (-6239538.500,-2824111.250) ( 74d25'23.24"E, 92d47'8.60"S)
Band 1 Block=15167x1 Type=UInt16, ColorInterp=Gray
  Min=0.000 Max=1023.000
  Minimum=0.000, Maximum=1023.000, Mean=196.147, StdDev=178.217
  Overviews: arbitrary
  Metadata:
    STATISTICS_MINIMUM=0
    STATISTICS_MAXIMUM=1023
    STATISTICS_MEAN=196.1465382488
    STATISTICS_STDDEV=178.21701451188
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/fwtools/attachments/20090602/ef1b4479/attachment.htm