2. Conversão de coordenadas com o programa CS2CS (PROJ.4)

A bibblioteca de projecções cartográficas PROJ.4 contém os programas PROJ e CS2CS. O programa PROJ efectua projecções cartográficas no sentido directo e no sentido inverso e é executado na linha de comando, definindo num string todas as características da projecção. A leitura e saída de dados faz-se pelo stdin/stdout (teclado/ecran), ou direccionando para ficheiros como no exemplo seguinte (projecção de Gauss do datum ETRS89 - sistema PT-TM06). O ficheiro de entrada deverá ter longitudes e latitudes ETRS89 e o de saída terá as coordenadas cartográficas resultantes.

        PROJ +ellps=WGS84 +proj=tmerc +lat_0=39d40'05.73"N +lon_0=08d07'59.19"W < file_in.txt > file_out.txt

O programa CS2CS converte datum e projecção entre dois sistemas de coordenadas geográficas ou cartográficas. Usa-se na linha de comando com dois strings, um para cada sistema de coordenadas.

        CS2CS string1 +to string2

Na página do PROJ.4 encontra-se documentação, manpages e FAQ sobre o uso dos dois programas.

Instalação

A forma mais simples de instalar este programa é através da instalação do programa FWTOOLS (versão Windows), da biblioteca GDAL. Este programa instala também o PROJ.4 e configura as variáveis de ambiente necessárias para correr os programas PROJ.EXE e CS2CS.EXE na linha de comando. Deverá ser corrida a shell (linha de comando) que ficou no ambiente de trabalho:

NOTA: Alternativamente pode instalar-se apenas o PROJ, de preferência na raíz do disco C. Os executáveis ficam em C:\PROJ\BIN enquanto que os ficheiros de configuração ficam em C:\PROJ\NAD. É necessário acrescentar a primeira ao PATH e criar uma variável de ambiente para a segunda (PROJ_LIB=C:\PROJ\NAD).

Configuração

Para configurar o programa para uso com as conversões com grelhas, colocar o conteúdo do ficheiro grelhas.zip (contém os 4 ficheiros de extensão GSB e o ficheiro PT que contém as definições dos sistemas nacionais) na pasta: C:\Programas\FWTools2.2.8\proj_lib. Os "strings" para os diferentes sistemas são os seguintes:

 
STRING            Sistema de coordenadas Código EPSG
etrs89 Geográficas ETRS89 4258
pttm06 PTTM06 do datum ETRS89 3763
pttm06mil  PTTM06 do datum ETRS89 (versão militar) N/A
etrs89utm  UTM-29N do datum ETRS89 25829
d73geo     Geográficas Datum 73 4274
d73hg      Hayford-Gauss datum 73 (IGP) 27493
d73hgmil   Hayford-Gauss datum 73 (IGP) + translaccao militar N/A
dlxgeo     Geográficas Datum Lisboa - Elips. Hayford 4207
dlxhg      Hayford-Gauss datum Lisboa 20791
dlxhgmil   Hayford-Gauss datum Lisboa (militar) 20790
ed50geo    Geográficas ED50 4666
ed50utm29  UTM 29-N Datum Europeu 1950 2963
dlxbgeo    Geográficas Datum Lisboa - Elipsoide de Bessel 4230
dlxbb      Sistema Bessel-Bonne (eixos invertidos) 23029

NOTA: Caso se pretenda uma transformação de coordenadas com transformação de Bursa-Wolf ou Molodensky pode usar-se os ficheiros PT7 ou PT3 respectivamente (também incluídos no ficheiro grelhas.zip). Os nomes dos sistemas de coordenadas são os mesmos. Os erros destas transformações são bastantes superiores aos das grelhas, especialmente no caso do datum Lisboa.

Exemplos:

Converter coordenadas militares Datum Lisboa para UTM datum europeu:

        cs2cs +init=pt:dlxhgmil +to +init=pt:ed50utm29

Converter datum 73 projectado para geográficas ETRS89 (input e output redireccionados para ficheiros):

      cs2cs +init=pt:d73hg +to +init=pt:etrs89 < in_file.txt > out_file.txt

Podem ser usados códigos EPSG (em alternativa às definições do ficheiro PT), combinados com a conversão de datum. Por exemplo, converter coordenadas HG datum 73 para geográficas WGS84 (tomado como igual ao ETRS89)

      cs2cs +init=epsg:27492 +nadgrids=pt73_e89.gsb +to +init=epsg:4326

 

Pontos de teste:

O ficheiro em formato Excel PONTOS_TESTE.XLS contém uma amostra de pontos da rede gódésica com coordenadas nos diferentes sistemas. As coordenadas foram obtidas na página do IGP. Todos os pontos da rede geodésica podem ser obtidos em:

Redes de 1ª e 2ª ordem (coordenadas ETRS89):   http://www.igeo.pt/produtos/geodesia/vg/rgn/rgn.asp 

Redes de 1ª, 2ª e 3ª ordem (sistemas tradicionais): http://www.igeo.pt/produtos/geodesia/vg/rgn/docs/osys.asp

As coordenadas obtidas por transformação diferirão por alguns centímetros ou decímetros das dadas.