[Portugal] CRS 3763 no PostGIS

Mon Jul 13 14:19:37 EDT 2009

Viva,

De forma muito telegráfica (de férias) aproveito para informar que o GISVM Server vem com os CRS portugueses: 20790, 20791, 27492, 3763 redefinidos no PostGIS. Foram de acordo com a lista que publiquei aqui logo apos o SASIG 2008:
http://spreadsheets.google.com/pub?key=pqto3s91EMthHVKvbNX7MVQ&output=html

Basicamente foi executada o script: sqri_pt.sql (que é distribuido com o GISVM Server)
http://docs.google.com/View?id=dc3bfsd6_1gj8xg5qc

Tomei a liberdade de substituir os CRS originais do EPSG, por aquelas definições que suponho correctas.(oficiais do IGP)
Caso detectem erros ou tenham uma melhor sugestão, agradeço que me comuniquem.
Obrigado!

Entretanto ainda sobre CRS, acrescento que o PostGIS do GISVM Server é distribuido como "ZIGGIS ready"!
O que significa que vem com todas as definições de coordenadas (CRS) da ESRI (ArcGIS) e portanto pode ser usado com o ArcGIS + ZigGIS!!!

Cumprimentos,
Ricardo Pinho

________________________________
De: Luís de Sousa <luis.a.de.sousa  gmail.com>
Para: Ricardo Sena <ricardogsena  gmail.com>
Cc: portugal  lists.osgeo.org
Enviadas: Sexta-feira, 10 de Julho de 2009 16:34:19
Assunto: Re: [Portugal] CRS 3763 no PostGIS

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

      ____________________________________________________________________________________
Veja quais são os assuntos do momento no Yahoo! +Buscados
http://br.maisbuscados.yahoo.com
-------------- próxima parte ----------
Um anexo em HTML foi limpo...
URL: http://lists.osgeo.org/pipermail/portugal/attachments/20090713/caa4de63/attachment.html