[Portugal] CRS 3763 no PostGIS

[Luís de Sousa]

Ricardo,

Obrigado pela resposta detalhada, de facto tirando o prefixo 9 tudo
funciona correctamente.

A informação sobre o método das grelhas com o postGIS também vem a
calhar, obrigado.

Luís

2009/7/10 Ricardo Sena <ricardogsena at gmail.com>:
> 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
>
>