<div>Caro Luís Sousa,</div>
<div> </div>
<div>Uma vez que se trata de um CRS que não vem registado por defeito na instalação do PostGIS, julgo que utilizam o prefixo 9 para indicar que este sistema foi definido pelo utilizador. O mesmo se passa com a nova definição do CRS relativo ao Datum 73 (substitui a anterior definição do 27492), que está definido pelo EPSG com o código 27493 e é inserido com o SRID 927493. Recomendo que substitua os códigos com prefixo 9 pelos códigos originais do EPSG (27492 em vez de 927492). </div>
<div> </div>
<div>Chamo a atenção para o facto de as definições dos CRS portugueses feitas pelo EPSG não incluírem o parâmetro relativo à transformação do datum local para o datum global (Datum 73 para ETRS89, por exemplo), o que provoca a introdução de erros de várias dezenas de metros (mais de 100 metros) quando se faz transformação de coordenas entre sistemas com diferentes datuns. Recomendo por isso que introduza nas definições que estão na tabela Spatial_Ref_Sys os valores desse parâmetro. </div>
<div> </div>
<div>Para uma transformação Bursa-Wolf (7 parâmetros), a definição do sistema EPSG:27492 seria a seguinte (texto a acrescentar assinalado a amarelo):</div>
<div> </div>
<div>
<p class="western"><font size="2">+proj=tmerc </font><font size="2">+lat_0=39.66666666666666 </font><font size="2">+lon_0=-8.131906111111112 </font><font size="2">+k=1.000000 </font><font size="2">+x_0=180.598 </font><font size="2">+y_0=-86.98999999999999 </font><font size="2">+ellps=intl <font size="2">+units=m </font></font><font size="2"><span style="BACKGROUND: rgb(255,255,0)">+towgs84=-231.03,102.61,26.84,0.61,-0.20,0.88,1.79</span></font></p>
</div>
<div>Esta definição é a que está no campo <strong>proj4text</strong>. Uma vez que o PostGIS se baseia na biblioteca PROJ.4 para fazer a transformação de coordenadas, é a definição que está neste campo que é utilizada para fazer as transformações.</div>
<div> </div>
<div>O sistema EPSG:3763 deve ter a seguinte definição:</div>
<div> </div>
<div><font size="2">+proj=tmerc +lat_0=39.6682583333333 +lon_0=-8.13310833333333 +ellps=GRS80 +towgs84=0,0,0 +units=m +no_defs</font><font size="3"> </font></div>
<div> </div>
<div>A transformação de coordenadas também pode ser realizada através do método das grelhas. O Prof. José Alberto Gonçaves, da Universidade do Porto, produziu as grelhas para a passagem entre os sistemas mais utilizados em Portugal (datum Lisboa, datum 73, ED50, bessel) para ETRS89. De acordo com a informação disponibilizada sobre a utilização destes método, a transformação com grelhas tem erro médio quadrático de 5 cm nas duas coordenadas no caso do datum 73 e 8 cm no caso do datum Lisboa. Para configurar a utilização deste método no PostGIS, é necessário executar os seguintes procedimentos:</div>
<div> </div>
<div>1 - Fazer o download das grelhas (<a href="http://www.fc.up.pt/pessoas/jagoncal/coordenadas/grelhas.zip">http://www.fc.up.pt/pessoas/jagoncal/coordenadas/grelhas.zip</a>) e instalar os ficheiro na directoria proj do PostGIS (C:\Programas\PostgreSQL\8.3\share\contrib\postgis\proj)</div>
<div> </div>
<div>2 - Alterar as definições dos sistemas portugueses que estão no campo proj4text da tabela Spatial_Ref_Sys, de modo a incluir a referência para o respectivo ficheiro da grelha:</div>
<div> </div>
<div>EPSG:27492</div>
<div>+proj=tmerc <font size="2">+lat_0=39.66666666666666 </font><font size="2">+lon_0=-8.131906111111112 </font><font size="2">+k=1.000000 </font><font size="2">+x_0=180.598 </font><font size="2">+y_0=-86.98999999999999 </font><font size="2">+ellps=intl <font size="2">+units=m +nadgrids=pt73_e89.gsb</font></font></div>
<div> </div>
<div>EPSG:27493</div>
<div>+proj=tmerc <font size="2">+lat_0=39.66666666666666 </font><font size="2">+lon_0=-8.131906111111112 </font><font size="2">+k=1.000000 </font><font size="2">+x_0=180.598 </font><font size="2">+y_0=-86.98999999999999 </font><font size="2">+ellps=intl <font size="2">+units=m <font style="BACKGROUND-COLOR: #ffff33">+nadgrids=pt73_e89.gsb</font></font></font></div>
<div> </div>
<div>
<div>EPSG:20790</div>
<div><font size="2">+proj=tmerc </font><font size="2">+lat_0=39.66666666666666 </font><font size="2">+lon_0=-8.131906111 </font><font size="2">+k=1 </font><font size="2">+x_0=200000 </font><font size="2">+y_0=300000 </font><font size="2">+ellps=intl +units=m <font style="BACKGROUND-COLOR: #ffff33">+nadgrids=ptLX_e89.gsb</font></font></div>
<div> </div>
<div>
<div>EPSG:20791</div>
<div><font size="2">+proj=tmerc </font><font size="2">+lat_0=39.66666666666666 </font><font size="2">+lon_0=-8.131906111 </font><font size="2">+k=1 </font><font size="2">+x_0=200000 </font><font size="2">+y_0=300000 </font><font size="2">+ellps=intl +units=m <font style="BACKGROUND-COLOR: #ffff33">+nadgrids=ptlx_e89.gsb</font></font></div>
</div>
<div> </div>
<div>Recomendo uma visita à página do Proj José Alberto Gonçalves (<a href="http://www.fc.up.pt/pessoas/jagoncal/coordenadas">http://www.fc.up.pt/pessoas/jagoncal/coordenadas</a>) e do Grupo de Trabalho que foi consituído para compatibilizar as definições dos CRS portugueses do EPSG com a utilização de software livre para SIG (<a href="http://spreadsheets.google.com/ccc?key=pqto3s91EMthHVKvbNX7MVQ&hl=en">http://spreadsheets.google.com/ccc?key=pqto3s91EMthHVKvbNX7MVQ&hl=en</a>). Para além de ser o elemento dinamizador deste grupo, julgo que o Ricardo Pinho também já colocou estas definições na GISVM.</div>
<div> </div>
<div>Cumprimentos,</div>
<div> </div>
<div>Ricardo Sena</div>
<div> </div></div>
<div> </div>
<div class="gmail_quote">2009/7/10 Luís de Sousa <span dir="ltr"><<a href="mailto:luis.a.de.sousa@gmail.com">luis.a.de.sousa@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Vivam,<br><br>Tenho estado a usar o CRS 3763 por forma a por termo às projecções<br>nacionais obsoletas. No postGIS este CRS não está incluído de origem e<br>
tem de ser inserido na tabela public.spatial_ref_sys. O sítio<br>SpatialReferences.org disponibiliza o comando insert que faz esta<br>operação:<br><br><a href="http://www.spatialreferences.org/ref/epsg/3763/" target="_blank">http://www.spatialreferences.org/ref/epsg/3763/</a><br>
<br>Há no entanto algo de estranho com este comando, o campo srid é<br>diferent do campo auth_srid. Em nenhum dos registos originais do<br>postGIS isto é assim, em todos eles os dois campos são iguais.<br><br>O postGIS usa o campo srid como chave primária, mas certas aplicações<br>
usam o campo auth_srid como referência (e.g. gvSIG, QGIS) e isso<br>tem-me causado diversos problemas.<br><br>Existe alguma razão para este o sítio SpatialReferences.org<br>disponibilizar o comando desta forma? Ou isto é um erro que deve ser<br>
relatado?<br><br>Obrigado,<br><br>Luís<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>