Olá Boa Noite,<br><br>Junto anexo a expressão .proj4 que (julgo) descreve o sistema que é usado actualmente - um 20790 transladado. É um bocado esquisito, mas só assim consegui que os pontos batessem onde era suposto (Castro Zambujal, Torres Vedras)<br>
<br>Se esta expressão estiver correcta, os parâmetros de transformação que o Prof. calculou continuam válidos? Uma transformação afim continua a poder ser utilizada?<br><br>Obrigado pela preciosa ajuda<br><br>cumprimentos,<br>
<br>André Mano<br><br><div class="gmail_quote">2012/8/3 Jose Gonçalves <span dir="ltr"><<a href="mailto:jagoncal@gmail.com" target="_blank">jagoncal@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Olá<br><br>Com pontos de controlo com coordenadas conhecidas nos dois sistemas pode-se obter uma relação de tranasformação directa entre os dois sistemas.<br>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:<br>
<br><span style="font-family:courier new,monospace">9105 71.23 40.69 100336.19 234735.28</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">9602 -0.79 9.21 100265.42 234701.08</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace"></span><br>Determinando uma transformação afim entre os dois sistemas obtém-se:<br><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">X_dlx_mil = 100266.56 + 0.99927416 * Xlocal - 0.03806220 * Ylocal</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">Y_dlx_mil = 234691.91 + 0.03806220 * Xlocal + 0.99927416 * Xlocal</span><br style="font-family:courier new,monospace"><br>Desta forma qualquer ponto no sistema local pode ser transformado e usado na georreferenciação em coordenadas militares.<br>
<br>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.<br>
<br>Cumprimentos<br><br>José Alberto Gonçalves<br><br><br><br><div class="gmail_quote">Em 2 de agosto de 2012 11:03, Andre Mano <span dir="ltr"><<a href="mailto:andre.s.mano@gmail.com" target="_blank">andre.s.mano@gmail.com</a>></span> escreveu:<div>
<div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Bom Dia,<br><br>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.<br>
<br>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.<br>
<br>A minha ideia era a seguinte:<br><br>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<br><br>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<br>
<br>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)<br><br>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)<br>
<br><br>Tudo isto porque o a rocha onde foi marcado o 0 da triangulação dos anos 60 desapareceu!<div><div><br><br><br><br><div class="gmail_quote">2012/8/2 Jose Gonçalves <span dir="ltr"><<a href="mailto:jagoncal@gmail.com" target="_blank">jagoncal@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Caro André<br><br>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 (<a href="http://trac.osgeo.org/proj/wiki/GenParms" target="_blank">http://trac.osgeo.org/proj/wiki/GenParms</a>).<br>
<br>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.<br>
Qual a região dessa cartografia? Tem algum ponto com coordenadas (0,0)?<br>Dispondo de pontos de controlo nessa cartografia, com coordenadas num sistema atual, o que há fazer é georreferenciar no sistema atual e esquecer o antigo.<br>
<br>Cumprimentos<br><br>José Alberto Gonçalves<br><br><br><div class="gmail_quote">Em 1 de agosto de 2012 20:40, Andre Mano <span dir="ltr"><<a href="mailto:andre.s.mano@gmail.com" target="_blank">andre.s.mano@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><br>Olá Boa noite,<br><br>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.<br>
<br>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).<br>
<br>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!)<br>
<br>É possível expressar esta diferença através da biblioteca .proj4? Alguém me pode dar exemplos?<br><br>Qualquer dica é bem vinda!!<br><br>cumprimentos<span><font color="#888888"><br><br>-- <br>..................................<br>
André Mano<br><a style="color:rgb(0,0,153)" href="http://opussig.blogspot.com/" target="_blank"><font size="1">http://opussig.blogspot.com/</font></a><br>
<br>
</font></span><br></div></div>_______________________________________________<br>
Portugal mailing list<br>
<a href="mailto:Portugal@lists.osgeo.org" target="_blank">Portugal@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/portugal" target="_blank">http://lists.osgeo.org/mailman/listinfo/portugal</a><br>
<br></blockquote></div><br>
<br>_______________________________________________<br>
Portugal mailing list<br>
<a href="mailto:Portugal@lists.osgeo.org" target="_blank">Portugal@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/portugal" target="_blank">http://lists.osgeo.org/mailman/listinfo/portugal</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>..................................<br>André Mano<br><a style="color:rgb(0,0,153)" href="http://opussig.blogspot.com/" target="_blank"><font size="1">http://opussig.blogspot.com/</font></a><br>
<br>
</div></div><br>_______________________________________________<br>
Portugal mailing list<br>
<a href="mailto:Portugal@lists.osgeo.org" target="_blank">Portugal@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/portugal" target="_blank">http://lists.osgeo.org/mailman/listinfo/portugal</a><br>
<br></blockquote></div></div></div><br>
<br>_______________________________________________<br>
Portugal mailing list<br>
<a href="mailto:Portugal@lists.osgeo.org">Portugal@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/portugal" target="_blank">http://lists.osgeo.org/mailman/listinfo/portugal</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>..................................<br>André Mano<br><a style="color:rgb(0,0,153)" href="http://opussig.blogspot.com/" target="_blank"><font size="1">http://opussig.blogspot.com/</font></a><br>
<br>