[Portugal] Parâmetros transformação com .proj4

[Jose Gonçalves]

Olá

Com pontos de controlo com coordenadas conhecidas nos dois sistemas pode-se
obter uma relação de tranasformação directa entre os dois sistemas.
Em primeiro lugar, pela localização, dá para perceber que as coordenadas
militares (epsg:20790) estão com uma translação de 50 km em X e 200 km em
Y. Dois dos pontos de controlo corrigidos serão:

9105  71.23  40.69 100336.19   234735.28
9602  -0.79   9.21 100265.42   234701.08

Determinando uma transformação afim entre os dois sistemas obtém-se:

X_dlx_mil = 100266.56 + 0.99927416 * Xlocal - 0.03806220 * Ylocal
Y_dlx_mil = 234691.91 + 0.03806220 * Xlocal + 0.99927416 * Xlocal

Desta forma qualquer ponto no sistema local pode ser transformado e usado
na georreferenciação em coordenadas militares.

Havia a possibilidade de definir uma projeção a simular o sistema local.
Podia-se estabelecer por exemplo uma projecção transversa de Mercator com
um meridiano central tal que a convergência dos meridianos fosse igual à
diferença angular entre os dois sistemas, com um factor de escala de forma
a ter a mesma escala do sistema militar e uma translação de origem adequada
para se gerar as coordenadas do sistema local. Penso que não vale a pena
porque ainda daria algum trabalho.

Cumprimentos

José Alberto Gonçalves



Em 2 de agosto de 2012 11:03, Andre Mano <andre.s.mano  gmail.com> escreveu:

> Bom Dia,
>
> O sistema foi estabelecido por por triangulação topográfica, sim. Havia um
> ponto 0 (entretanto desaparecido) mas que não é um dos pontos de controlo
> que possuo.
>
> Preciso de desmontar o sistema antigo porque para a maioria dos desenhos -
> nalguns casos desenhos que apenas cobrem 4m2 não tenho forma de usar pontos
> de controlo no sistema 20790 (o actual). A minha ideia era, por isso, criar
> o meu próprio sistema de referência, georeferenciar e depois, na última
> fase transformar para 20790 e assim condensar toda a informação num único
> sistema.
>
> A minha ideia era a seguinte:
>
> 1 - aproveitar a string .proj4 que define o sistema 20790 e adaptar -
> Modificar a falsa origem de modo a que coincidisse com o ponto 0
>
> 2 - Usando a mesma projecção cartográfica não haveria distorções de maior
> (erro até 5cm é aceitável) no momento de transformar tudo para o sistema
> 20790
>
> 3 - Procurar uma constante de diferença nos 7 pontos de controlo (x60 -
> x90) que me permita aferir as coordenadas do ponto 0 expressas no sistema
> 20790 (sem deixar de ter em conta o efeito da rotação)
>
> 4 - Tenho um levantamento feito já nos anos 90, onde é possível medir a
> diferença angular entre o Norte dos anos 60 e o norte cartográfico do
> sistema 20790 (envio em pvt)
>
>
> Tudo isto porque o a rocha onde foi marcado o 0 da triangulação dos anos
> 60 desapareceu!
>
>
>
>
> 2012/8/2 Jose Gonçalves <jagoncal  gmail.com>
>
>> Caro André
>>
>> A biblioteca libproj4, desenvolvida originalmente por Gerald Evenden,
>> incluía um parâmetro +rot para rotações de eixos, mas penso que a
>> biblioteca PROJ.4 não considera esse parâmetro (
>> http://trac.osgeo.org/proj/wiki/GenParms).
>>
>> De qualquer forma a definição desse sistema de coordenadas antigo na
>> sintaxe do PROJ.4 obriga a conhecer todos os detalhes da projeção
>> cartográfica (elipsóide, projeção, parâmetros da projeção, etc.). É mais
>> provável que esse sistema não resulte de uma projeção cartográfica mas sim
>> de uma triangulação topográfica local. Em casos que vi de cartografia
>> relativamente antiga (por exemplo cartografia de escala 1:1000 de Coimbra,
>> de 1932) era isso que acontecia.
>> Qual a região dessa cartografia? Tem algum ponto com coordenadas (0,0)?
>> Dispondo de pontos de controlo nessa cartografia, com coordenadas num
>> sistema atual, o que há fazer é georreferenciar no sistema atual e esquecer
>> o antigo.
>>
>> Cumprimentos
>>
>> José Alberto Gonçalves
>>
>>
>> Em 1 de agosto de 2012 20:40, Andre Mano <andre.s.mano  gmail.com>escreveu:
>>
>>>
>>> Olá Boa noite,
>>>
>>> Estou a tentar criar os parâmetros de transformação de uma projecção
>>> local, usada nos anos 60, para uma actual (EPSG20790), de modo a poder
>>> georeferenciar desenhos antigos.
>>>
>>> Trata-se de uma área de 5ha para a qual tenho 7 pontos de controlo. Para
>>> cada um desses pontos tenho as coordenadas no sistema antigo e no sistema
>>> que é usado desde os anos 90 (20790).
>>>
>>> O problema é que o Norte usado nos anos 60 não foi o cartográfico  - tem
>>> uma declição de aproximadamente 2.º para a esquerda relativamente ao norte
>>> cartográfico usado nos desenhos feitos desde os anos 90, o que implica que
>>> não posso simplesmente aplicar uma falsa origem - tenho que ter em conta um
>>> parâmetro de rotação (digo eu!)
>>>
>>> É possível expressar esta diferença através da biblioteca .proj4? Alguém
>>> me pode dar exemplos?
>>>
>>> Qualquer dica é bem vinda!!
>>>
>>> cumprimentos
>>>
>>> --
>>> ..................................
>>> André Mano
>>> http://opussig.blogspot.com/
>>>
>>>
>>> _______________________________________________
>>> 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
>>
>>
>
>
> --
> ..................................
> André Mano
> http://opussig.blogspot.com/
>
>
> _______________________________________________
> Portugal mailing list
> Portugal  lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/portugal
>
>
-------------- próxima parte ----------
Um anexo em HTML foi limpo...
URL: <http://lists.osgeo.org/pipermail/portugal/attachments/20120803/8dd7de97/attachment.html>