[Portugal] CRS 3763 no PostGIS

[Ricardo Sena]

Caro Luís Sousa,

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).

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.

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):


+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112
+k=1.000000 +x_0=180.598
+y_0=-86.98999999999999 +ellps=intl +units=m
+towgs84=-231.03,102.61,26.84,0.61,-0.20,0.88,1.79
Esta definição é a que está no campo *proj4text*. 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.

O sistema EPSG:3763 deve ter a seguinte definição:

+proj=tmerc +lat_0=39.6682583333333 +lon_0=-8.13310833333333 +ellps=GRS80
+towgs84=0,0,0 +units=m +no_defs

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:

1 - Fazer o download das grelhas (
http://www.fc.up.pt/pessoas/jagoncal/coordenadas/grelhas.zip) e instalar os
ficheiro na directoria proj do PostGIS
(C:\Programas\PostgreSQL\8.3\share\contrib\postgis\proj)

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:

EPSG:27492
+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112
+k=1.000000 +x_0=180.598
+y_0=-86.98999999999999 +ellps=intl +units=m +nadgrids=pt73_e89.gsb

EPSG:27493
+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112
+k=1.000000 +x_0=180.598
+y_0=-86.98999999999999 +ellps=intl +units=m +nadgrids=pt73_e89.gsb

 EPSG:20790
+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111 +k=1
+x_0=200000 +y_0=300000
+ellps=intl +units=m  +nadgrids=ptLX_e89.gsb

 EPSG:20791
+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111 +k=1
+x_0=200000 +y_0=300000
+ellps=intl +units=m  +nadgrids=ptlx_e89.gsb

Recomendo uma visita à página do Proj José Alberto Gonçalves (
http://www.fc.up.pt/pessoas/jagoncal/coordenadas) 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 (
http://spreadsheets.google.com/ccc?key=pqto3s91EMthHVKvbNX7MVQ&hl=en). Para
além de ser o elemento dinamizador deste grupo, julgo que o Ricardo Pinho
também já colocou estas definições na GISVM.

Cumprimentos,

Ricardo Sena


2009/7/10 Luís de Sousa <luis.a.de.sousa at gmail.com>

> Vivam,
>
> Tenho estado a usar o CRS 3763 por forma a por termo às projecções
> nacionais obsoletas. No postGIS este CRS não está incluído de origem e
> tem de ser inserido na tabela public.spatial_ref_sys. O sítio
> SpatialReferences.org disponibiliza o comando insert que faz esta
> operação:
>
> http://www.spatialreferences.org/ref/epsg/3763/
>
> Há no entanto algo de estranho com este comando, o campo srid é
> diferent do campo auth_srid. Em nenhum dos registos originais do
> postGIS isto é assim, em todos eles os dois campos são iguais.
>
> O postGIS usa o campo srid como chave primária, mas certas aplicações
> usam o campo auth_srid como referência (e.g. gvSIG, QGIS) e isso
> tem-me causado diversos problemas.
>
> Existe alguma razão para este o sítio SpatialReferences.org
> disponibilizar o comando desta forma? Ou isto é um erro que deve ser
> relatado?
>
> Obrigado,
>
> Luís
> _______________________________________________
> Portugal mailing list
> Portugal at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/portugal
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/portugal/attachments/20090710/21332d6c/attachment-0001.html