[Portugal] Raster deslocado no QGis

[cremat0rio]

Caros,

Tenho aqui um problemita (que até pode ser muito básico) que é o seguinte:
Tenho uma Layer configurada no Geoserver com a BBOX:
33.492,16  -290.883,298 ; 35.296,322  -289.465,47

O Gdalinfo devolve-me:
/Driver: GTiff/GeoTIFF
Files: 081001_a-20791.tif
Size is 7103, 5582
Coordinate System is:
PROJCS["Lisbon (Lisbon) / Portuguese Grid",
    GEOGCS["Lisbon (Lisbon)",
        DATUM["Lisbon_1937_Lisbon",
            SPHEROID["International 1924",6378388,297.000000000005,
                AUTHORITY["EPSG","7022"]],
            AUTHORITY["EPSG","6803"]],
        PRIMEM["Lisbon",-9.131906111111112],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4803"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",39.66666666666666],
    PARAMETER["central_meridian",1],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","20791"]]
Origin = (33492.160000000003492,-289465.469999999972060)
Pixel Size = (0.254000000000000,-0.254000000000000)
Metadata:
  TIFFTAG_XRESOLUTION=200
  TIFFTAG_YRESOLUTION=200
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (   33492.160, -289465.470) (  1d22'35.54"E, 37d 3'30.32"N)
Lower Left  (   33492.160, -290883.298) (  1d22'35.31"E, 37d 2'44.33"N)
Upper Right (   35296.322, -289465.470) (  1d23'48.55"E, 37d 3'30.08"N)
Lower Right (   35296.322, -290883.298) (  1d23'48.32"E, 37d 2'44.09"N)
Center      (   34394.241, -290174.384) (  1d23'11.93"E, 37d 3'7.21"N)/

O que aparentemente está ok.

O QGis faz o pedido pelo GetCapabilities que devolve, dentre outras coisas:

/      <Layer queryable="1">
        <Name>vectorizacaoCadastroRustico:081001_a-20791</Name>
        <Title>081001_a-20791</Title>
        <Abstract/>
        <KeywordList>
          <Keyword>WCS</Keyword>
          <Keyword>GeoTIFF</Keyword>
          <Keyword>081001_a-20791</Keyword>
        </KeywordList>
        <SRS>EPSG:20791</SRS>
        
        <LatLonBoundingBox minx="-7.755" miny="37.045" maxx="-7.735"
maxy="37.058"/>
        <BoundingBox SRS="EPSG:20791" minx="33492.16" miny="-290883.298"
maxx="35296.322" maxy="-289465.47"/>
        
      </Layer>/

O que aparentemente também está ok.

Acontece que a imagem que é apresentada no QGis está deslocada (para cima e
para a esquerda).
E daqui para a frente o pedido que o QGis faz ao geoserver (por WMS) é
qualquer coisa como:
/http://localhost/geoserver/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&BBOX=33626.832816,-291127.076530,35411.759255,-289676.181886&SRS=EPSG:20791&WIDTH=1024&HEIGHT=831&LAYERS=vectorizacaoCadastroRustico:081001_a-20791&STYLES=&FORMAT=image/jpeg&DPI=96/

Ou seja de alguma forma ele define outra BBOX que começa em
33626.832816,-291127.076530, alguns metros deslocado.

Já experimentei fazer o pedido WMS em vários datum e também em várias
versões do QGis (1.7.1, 1.8.0 e 1.9.0-Master fefb3d0) e acontece o mesmo.

Alguém faz uma ideia do que poderá estar a acontecer?


Muito obrigado.
Cumprimentos,
David




--
View this message in context: http://osgeo-org.1560.n6.nabble.com/Raster-deslocado-no-QGis-tp5006068.html
Sent from the OSGeo Portuguese Local Chapter mailing list archive at Nabble.com.