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 <[email protected]> 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 <[email protected]> > >> 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 <[email protected]>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 >>> [email protected] >>> http://lists.osgeo.org/mailman/listinfo/portugal >>> >>> >> >> _______________________________________________ >> Portugal mailing list >> [email protected] >> http://lists.osgeo.org/mailman/listinfo/portugal >> >> > > > -- > .................................. > André Mano > http://opussig.blogspot.com/ > > > _______________________________________________ > Portugal mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/portugal > >
_______________________________________________ Portugal mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/portugal
