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:
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.
| 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 |
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 |
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