Caro José, Penso que não compreendeu os meus comentários, não estou a fazer qualquer confusão entre PC e PAF. Notei várias vezes que nem o WKT nem o proj4 permitem a definição do PAF, o que tal como escreve, é deixado ao critério do utilizador aquando da invocação da rotina de transformação. Tal como escrevi anteriormente até esse momento o PAF e o posicionamento do elipsóide em geral são irrelevantes.
O que me ensiranam na escola foi que um sistema de referência é composto por dois subsistemas, um geográfico, definido por um elipsóide e o seu posicionamento relativamente ao geóide e um cartográfico, definido por uma função matemática de projecção da superfície do elipsóide num plano. Ora, se o mesmo ponto no geóide (e.g Observatório Astronómico) tem coordenadas geográficas diferentes então estamos perante dois sistemas diferentes. Claro que podemos assumir que em ambos os sistemas o posicionamento do elipsóide é o mesmo, sendo a diferença apenas na origem das coordenadas geográficas. Mas isso é apenas uma assumpção, nada na definição dos dois sistemas o indica, pois nenhum dos formatos permite essas definições. Bom, e isso é a origem desta longo conversa, o WKT e o proj4 são omissos na definição do sub-sistema geográfico. Saudação, Luís 2010/5/20 Jose Gonçalves <[email protected]>: > Caro Luís Sousa > >> +pm=lisbon está em directa contradicção com o que escreve o prof. >> Casaca. Enquanto que o observatório astronómico tem coordenadas >> geográficas 39º40' N, 8º 07' 54.862'' W no SHGM no EPSG20790 será > > Está a confudir duas coisas diferentes, que são o ponto astronómico > fundamental (PAF) e o ponto central (PC) da projecção. > O ponto em Lisboa, no Castelo de São Jorge, é o PAF, isto é um ponto onde, > por observação astronómica se fez a medição da latitude e longitude do > lugar. Esses valores foram: > > LON = 09d07'54.862"W > LAT = 38d42'43.631"N (e não 39d40') > > Nesse ponto considera-se as coordenadas geográficas geodésicas iguais às > astronómicas, ou seja "fixa-se o elipsóide". > A partir deste ponto, por triangulação, determinaram-se as coordenadas > geográficas geodésicas, no elipsóide de Hayford, de todos os vértices da > rede geodésica de 1ª ordem. Este conjunto de vértices, com as respectivas > coordenadas geográficas, materializam o datum Lisboa. > > Para efeito de cálculo topográfico e produção de cartografia estabeleceu-se > a projecção cartográfica e só aí é que aparece o PC. Esse ponto foi > estabelecido como sendo um ponto exactamente à longitude 1d00'00.000"E do > meridiano de Lisboa (ou melhor, do vértice geodésico existente no Castelo de > S. Jorge) e à latitude 39d40'00.000"N. Este ponto cai algures no centro do > país, provavelmente no meio do mato, por acaso relativamente próximo do > vértice geodésico da Melriça (a uns 3 km). Esse 1 grau a este de Lisboa faz > com que a longitude do PC, em relação ao meridiano internacional, seja > 8d07'54.862"W. > A projecção utilizada é a de Gauss (ou Gauss-Krueger ou Transversa de > Mercator, ou que se lhe quiser chamar, mas isso é outra discussão...), que > com a translacção miltar constitui o SHGM. > > O PC é importante para a configuração do sistema de coordenadas, mas o PAF > não. Quando muito, o PAF poderá ser considerado como meridiano de longitude > zero, como está no sistema EPSG:20790, e como foi usado na nossa cartografia > durante algum tempo . Se reparar nas folhas mais antigas da carta militar > (pelo menos nas edições de 1949) só aparecem na carta longitudes em relação > ao meridiano de Lisboa. > Espero que isto esclareça as dúvidas sobre a projecção. > > Em relação aos parâmetros de transformação de datum, eles indicam o > posicionamento do elipsóide do datum local em relação ao elipsóide do datum > global (não ao geóide). Na base de dados EPSG muitos sistemas (quase todos > os europeus) têm diversos conjuntos de parâmetros. A opção no software > PROJ.4 foi de, nesses casos, extrair da base de dados EPSG apenas a > informação do elipsóide e da projecção e não os parâmetros de datum. Esses > ficam ao critério do utilizador. > No caso do software ESRI, também ocorre algo semelhante: quando uma mudança > de coordenadas involve mudança de datum (uma "geographic transformation") é > pedida ao utilizador a escolha numa lista. > > Cumprimentos > > José Alberto Gonçalves > > > Em 20 de maio de 2010 17:30, Luís de Sousa <[email protected]> > escreveu: >> >> Caro Luís, >> >> Por favor leia os meis comentários em seguida: >> >> > <20790> +proj=tmerc +lat_0=39.66666666666666 +lon_0=1 +k=1 +x_0=200000 >> > +y_0=300000 +ellps=intl +pm=lisbon +units=m +no_defs <> >> > >> > em que lat_0= 39.66(6) é equivalente a 39º 40' N referido pelo Prof. >> > Casaca. >> >> +pm=lisbon está em directa contradicção com o que escreve o prof. >> Casaca. Enquanto que o observatório astronómico tem coordenadas >> geográficas 39º40' N, 8º 07' 54.862'' W no SHGM no EPSG20790 será >> 39º40' N, 1º00'E (veja o e-mail inicial do prof. J.A. Gonçalves) isto >> assumindo que o elipsóide está posicionado no mesmo ponto. >> >> > Os parâmetros +lon_0=1 e +pm=lisbon definem 8º 07' 54.862'' WGrw também >> > referido no livro citado. >> >> Não é a mesma coisa, no livro do prof. Casaca é explicitamente >> referido que o meridiano de origem do SHGM é o de Greenwhich. >> >> > Note que para transformar,por exemplo, coordenadas de Datum Lisboa para >> > WGS84, não precisa de saber a definição do Datum Lisboa. Precisa apenas >> > de >> > ter conhecimento dos parâmetros de transformação e garantir que as >> > coordenadas que quer transformar estão efectivamente em Datum Lisboa. >> >> Compreendo, os paramêtros de transformação são uma forma indirecta de >> indicar a orientação dos elipsóides em relação ao geóide, portanto o >> software poderá ignorar indicações como a do meridiano de origem. >> >> Saudação, >> >> Luís >> _______________________________________________ >> Portugal mailing list >> [email protected] >> http://lists.osgeo.org/mailman/listinfo/portugal > > > _______________________________________________ > Portugal mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/portugal > > _______________________________________________ Portugal mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/portugal
