[Portugal] Raster deslocado no QGis

[Jose Gonçalves]

Caro David

As diferenças nos bounding box são de poucas centenas de metros: 134.67 em
X e -243.78 em Y.
Com essa ordem de grandeza o erro é muito provavelmente a falta de uma
transformação de datum.
Na definição inicial (listada com o gdalinfo) não existe informação de
transformação de datum, isto é, algo como:
TOWGS84[-304.046,-60.576,103.64,0,0,0,0]]
Mesmo a ser feita essa transformação, ela terá sempre um erro de uns 2 ou 3
metros, e a imagem tem uma resolução superior (25.4 cm)
Essas dificuldades não existiriam se os dados fossem antes transformados
para ETRS89/PTTM06 com uma transformação adequada.

Cumprimentos

José Alberto Gonçalves




Em 3 de outubro de 2012 11:58, cremat0rio
<davidsantospinheiro  gmail.com>escreveu:

> 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.
> _______________________________________________
> Portugal mailing list
> Portugal  lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/portugal
>
-------------- próxima parte ----------
Um anexo em HTML foi limpo...
URL: <http://lists.osgeo.org/pipermail/portugal/attachments/20121003/3ecedf41/attachment.html>