Caro David<br><br>As diferenças nos bounding box são de poucas centenas de metros: 134.67 em X e -243.78 em Y.<br>Com essa ordem de grandeza o erro é muito provavelmente a falta de uma transformação de datum.<br>Na definição inicial (listada com o gdalinfo) não existe informação de transformação de datum, isto é, algo como:<br>
TOWGS84[-304.046,-60.576,103.64,0,0,0,0]]<br>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)<br>Essas dificuldades não existiriam se os dados fossem antes transformados para ETRS89/PTTM06 com uma transformação adequada.<br>
<br>Cumprimentos<br><br>José Alberto Gonçalves<br><br><br><br><br><div class="gmail_quote">Em 3 de outubro de 2012 11:58, cremat0rio <span dir="ltr"><<a href="mailto:davidsantospinheiro@gmail.com" target="_blank">davidsantospinheiro@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Caros,<br>
<br>
Tenho aqui um problemita (que até pode ser muito básico) que é o seguinte:<br>
Tenho uma Layer configurada no Geoserver com a BBOX:<br>
33.492,16 -290.883,298 ; 35.296,322 -289.465,47<br>
<br>
O Gdalinfo devolve-me:<br>
/Driver: GTiff/GeoTIFF<br>
Files: 081001_a-20791.tif<br>
Size is 7103, 5582<br>
Coordinate System is:<br>
PROJCS["Lisbon (Lisbon) / Portuguese Grid",<br>
GEOGCS["Lisbon (Lisbon)",<br>
DATUM["Lisbon_1937_Lisbon",<br>
SPHEROID["International 1924",6378388,297.000000000005,<br>
AUTHORITY["EPSG","7022"]],<br>
AUTHORITY["EPSG","6803"]],<br>
PRIMEM["Lisbon",-9.131906111111112],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4803"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",39.66666666666666],<br>
PARAMETER["central_meridian",1],<br>
PARAMETER["scale_factor",1],<br>
PARAMETER["false_easting",0],<br>
PARAMETER["false_northing",0],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","20791"]]<br>
Origin = (33492.160000000003492,-289465.469999999972060)<br>
Pixel Size = (0.254000000000000,-0.254000000000000)<br>
Metadata:<br>
TIFFTAG_XRESOLUTION=200<br>
TIFFTAG_YRESOLUTION=200<br>
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)<br>
AREA_OR_POINT=Area<br>
Image Structure Metadata:<br>
INTERLEAVE=BAND<br>
Corner Coordinates:<br>
Upper Left ( 33492.160, -289465.470) ( 1d22'35.54"E, 37d 3'30.32"N)<br>
Lower Left ( 33492.160, -290883.298) ( 1d22'35.31"E, 37d 2'44.33"N)<br>
Upper Right ( 35296.322, -289465.470) ( 1d23'48.55"E, 37d 3'30.08"N)<br>
Lower Right ( 35296.322, -290883.298) ( 1d23'48.32"E, 37d 2'44.09"N)<br>
Center ( 34394.241, -290174.384) ( 1d23'11.93"E, 37d 3'7.21"N)/<br>
<br>
O que aparentemente está ok.<br>
<br>
O QGis faz o pedido pelo GetCapabilities que devolve, dentre outras coisas:<br>
<br>
/ <Layer queryable="1"><br>
<Name>vectorizacaoCadastroRustico:081001_a-20791</Name><br>
<Title>081001_a-20791</Title><br>
<Abstract/><br>
<KeywordList><br>
<Keyword>WCS</Keyword><br>
<Keyword>GeoTIFF</Keyword><br>
<Keyword>081001_a-20791</Keyword><br>
</KeywordList><br>
<SRS>EPSG:20791</SRS><br>
<br>
<LatLonBoundingBox minx="-7.755" miny="37.045" maxx="-7.735"<br>
maxy="37.058"/><br>
<BoundingBox SRS="EPSG:20791" minx="33492.16" miny="-290883.298"<br>
maxx="35296.322" maxy="-289465.47"/><br>
<br>
</Layer>/<br>
<br>
O que aparentemente também está ok.<br>
<br>
Acontece que a imagem que é apresentada no QGis está deslocada (para cima e<br>
para a esquerda).<br>
E daqui para a frente o pedido que o QGis faz ao geoserver (por WMS) é<br>
qualquer coisa como:<br>
/<a href="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/" target="_blank">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/</a><br>
<br>
Ou seja de alguma forma ele define outra BBOX que começa em<br>
33626.832816,-291127.076530, alguns metros deslocado.<br>
<br>
Já experimentei fazer o pedido WMS em vários datum e também em várias<br>
versões do QGis (1.7.1, 1.8.0 e 1.9.0-Master fefb3d0) e acontece o mesmo.<br>
<br>
Alguém faz uma ideia do que poderá estar a acontecer?<br>
<br>
<br>
Muito obrigado.<br>
Cumprimentos,<br>
David<br>
<br>
<br>
<br>
<br>
--<br>
View this message in context: <a href="http://osgeo-org.1560.n6.nabble.com/Raster-deslocado-no-QGis-tp5006068.html" target="_blank">http://osgeo-org.1560.n6.nabble.com/Raster-deslocado-no-QGis-tp5006068.html</a><br>
Sent from the OSGeo Portuguese Local Chapter mailing list archive at Nabble.com.<br>
_______________________________________________<br>
Portugal mailing list<br>
<a href="mailto:Portugal@lists.osgeo.org">Portugal@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/portugal" target="_blank">http://lists.osgeo.org/mailman/listinfo/portugal</a><br>
</blockquote></div><br>