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 <[email protected]>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 > [email protected] > http://lists.osgeo.org/mailman/listinfo/portugal >
_______________________________________________ Portugal mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/portugal
