Raster projection shift

[Armin Burger]

Alessandro,

you could try the parameter +towgs84=-225,-65,9,0,0,0,0 for the 
projection definition of the Gauss-Boaga/Monte Mario data to force a 
datum shift (taken from 
http://dbwww.essc.psu.edu/lasdoc/overview/geomreg/dtm_list.html, "Datum 
: (210) ROME 1940") and see if the data match better then.

Armin

Alessandro Pasotti wrote:
> 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
> 
> 
> ------------------------------------------------------------------------
>