[Gdal-dev] Problems with warping

[Hugo Martins]

Hi gdal folks,

Well i have been experiencing some problems with gdalwarp. I downloaded some
imagery in jpg format. After i converted it to tif and no problem happened.
These images are in Lisboa Hayford Gauss IgeoE projection and i need to
reproject them to WGS84. Here are the gdalinfo from both files:

- before warping:

Driver: GTiff/GeoTIFF
Files: 41.tif
       41.aux
Size is 640, 480
Coordinate System is:
PROJCS["Lisboa_Hayford_Gauss_IGeoE",
    GEOGCS["GCS_Datum_Lisboa_Hayford",
        DATUM["unknown",
            SPHEROID["unnamed",6378388,297.0000000000014]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",39.66666666666666],
    PARAMETER["central_meridian",-8.131906111111112],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",200000],
    PARAMETER["false_northing",300000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (173002.500000000000000,282780.500000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  (  173002.500,  282780.500) (  8d26'44.94"W, 39d30'40.15"N)
Lower Left  (  173002.500,  282300.500) (  8d26'44.87"W, 39d30'24.59"N)
Upper Right (  173642.500,  282780.500) (  8d26'18.15"W, 39d30'40.23"N)
Lower Right (  173642.500,  282300.500) (  8d26'18.08"W, 39d30'24.66"N)
Center      (  173322.500,  282540.500) (  8d26'31.51"W, 39d30'32.41"N)
Band 1 Block=640x12 Type=Byte, ColorInterp=Red
Band 2 Block=640x12 Type=Byte, ColorInterp=Green
Band 3 Block=640x12 Type=Byte, ColorInterp=Blue


- after warping

Driver: GTiff/GeoTIFF
Files: 42.tif
       42.aux
Size is 693, 403
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.2572235630016,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (-8.445817277772575,39.510366486891535)
Pixel Size = (0.000010765710851,-0.000010765710851)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  (  -8.4458173,  39.5103665) (  8d26'44.94"W, 39d30'37.32"N)
Lower Left  (  -8.4458173,  39.5060279) (  8d26'44.94"W, 39d30'21.70"N)
Upper Right (  -8.4383566,  39.5103665) (  8d26'18.08"W, 39d30'37.32"N)
Lower Right (  -8.4383566,  39.5060279) (  8d26'18.08"W, 39d30'21.70"N)
Center      (  -8.4420870,  39.5081972) (  8d26'31.51"W, 39d30'29.51"N)
Band 1 Block=693x11 Type=Byte, ColorInterp=Red
Band 2 Block=693x11 Type=Byte, ColorInterp=Green
Band 3 Block=693x11 Type=Byte, ColorInterp=Blue

As you see, there is some shift with the WGS84 projection.
Could somebody help me on this?
Thanks,

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