Raster projection shift

[Alessandro Pasotti]

Hello,

I have a problem with a raster file ETM+, it is a GeoTIFF format, there is a 
shift (about 50m) with my shapefiles (ported to postgis), please see the 
attached image.

postgis srid is 26591

gdalinfo says:

-------8<------------------------------------

gdalinfo ETMPLUS.tif
Driver: GTiff/GeoTIFF
Size is 17868, 15764
Coordinate System is:
PROJCS["WGS 84 / UTM zone 32N",
    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",9],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32632"]]
Origin = (449380.875000,5051290.125000)
Pixel Size = (14.25000000,-14.25000000)
Metadata:
  TIFFTAG_XRESOLUTION=72
  TIFFTAG_YRESOLUTION=72
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Corner Coordinates:
Upper Left  (  449380.875, 5051290.125)
Lower Left  (  449380.875, 4826653.125)
Upper Right (  703999.875, 5051290.125)
Lower Right (  703999.875, 4826653.125)
Center      (  576690.375, 4938971.625)
Band 1 Block=17868x1 Type=Byte, ColorInterp=Gray

-------8<------------------------------------

Extract from the mapfile:

-------8<------------------------------------
  ################################################################
  #  Raster layer
  ################################################################
  LAYER
    NAME "etmplus"
    GROUP "Sfondo"
    DATA  "ETMPLUS.tif"
    
    TYPE         RASTER
    OFFSITE      0 0 0
    STATUS       OFF

    METADATA
      "wms_srs" "EPSG:32632"
      "wms_name" "etmplus"
      "wms_group" "Sfondo"
      "wms_server_version" "1.1.1"
      "wms_format" "image/jpeg"
      "typelabel" "etmplus"
    END
    PROJECTION
      "init=epsg:32632"
    END
  END
  #################################################################
  # Start of layer - [Comuni]
  #################################################################
  LAYER
    NAME "comuni"
    GROUP Cartografia
    STATUS ON
    DEBUG ON
    
    DATA "the_geom from comuni"
    CONNECTIONTYPE postgis
    CONNECTION "user=--- dbname=--- host=localhost password=---"
    TYPE POLYGON
    UNITS METERS
    DUMP TRUE
    METADATA
      "wms_title"                 "Comuni"
      "wms_srs"                   "EPSG:26591"
      "wfs_title"                 "Comuni"
    END
    TOLERANCE 3
    LABELITEM "nome"
    CLASS
      TEMPLATE 'ttt'
      NAME "Comune"
      SYMBOL "circle"
      SIZE 3
      OUTLINECOLOR 200 0 0
      LABEL
          TYPE TRUETYPE
          FONT "fritqat"
          SIZE 10
          BACKGROUNDCOLOR 255 255 255
          BACKGROUNDSHADOWCOLOR 112 112 112
          MINSIZE 4
          MAXSIZE 256
          POSITION AUTO
          OFFSET 0 0
          ANGLE AUTO
          BUFFER 0
          MINDISTANCE -1
          MINFEATURESIZE -1
          COLOR 0 0 0
          OUTLINECOLOR 255 255 255
          ANTIALIAS TRUE
          PARTIALS FALSE
          FORCE FALSE
      END
      # End of Label
    END
    # End of Class
  END
  # End of Layer comuni
-------8<------------------------------------

I really don't know where to find a solution.

Thanks for any help.

Alessandro Pasotti
--
itOpen - "Open Solutions for the Net Age"
w3: www.itopen.it
ICQ: 245871392
Linux User# 167502
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mapserv-shift.png
Type: image/png
Size: 13622 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20050908/bbe01443/attachment.png>