Difference between revisions of "Transformação de coordenadas cartográficas"

From OSGeo
Jump to navigation Jump to search
(New page: '''Transformação de coordenadas Cartográficas DE e PARA o novo sistema de projecção oficial de Portugal Continental: o PT-TM06''' Fonte: [http://wiki.osgeo.org/wiki/User:Rpinho Ricar...)
 
Line 1: Line 1:
−
'''Transformação de coordenadas Cartográficas DE e PARA o novo sistema de projecção oficial de Portugal Continental: o PT-TM06'''
+
= Transformação de coordenadas cartográficas de e para o novo sistema de projecção oficial de portugal continental: o pt-tm06 =
  
−
Fonte: [http://wiki.osgeo.org/wiki/User:Rpinho Ricardo Pinho] (http://gisvm.com/forum/index.php?topic=38.0)
+
'''Fonte: [http://wiki.osgeo.org/wiki/User:Rpinho Ricardo Pinho] (http://gisvm.com/forum/index.php?topic=38.0)'''
  
 +
* PT-MT06, baseado no ETRS89 e indicado pelo IGP para novo Sistema de Coordenadas Cartográficas Oficial Português:
 +
http://www.igeo.pt/produtos/Geodesia/Inf_tecnica/sistemas_referencia/Datum_ETRS89.htm
 +
 +
* Página do Prof. José Alberto Gonçalves, da Faculdade de Ciências da Universidade do Porto, que propõe a transformação de coordenadas pelo método das Grelhas.
 +
http://www.fc.up.pt/pessoas/jagoncal/coordenadas/
 +
 +
* Ferramentas utilizadas: [http://proj.maptools.org/man_cs2cs.html CS2CS] e [http://www.gdal.org/ogr2ogr.html OGR2OGR], que se podem encontrar no pacote [http://fwtools.maptools.org/ fwtools] ou directamente instaladas com as bibliotecas [http://trac.osgeo.org/proj/ Proj.4] e [http://www.gdal.org/ogr/ OGR]
 +
 +
== DESCRIÇÃO DOS UTILITÁRIOS CS2CS e OGR2OGR ==
 +
 +
* CS2CS, é um comando para usar em linhas de comandos (terminal/consola) que permite a conversão de coordenadas de um ponto ou conjunto de pontos (num ficheiro) entre dois sistemas de coordenadas indicados (Coordinate System).
 +
 +
* OGR2OGR, é um comando para usar em linhas de comandos (terminal/consola) que permite a conversão de uma diversidade grande de formatos vectoriais abrangidos pela biblioteca OGR. Vamos aqui detalhar a possibilidade de o utilizar para conversão de coordenadas no mesmo formato, nomeadamente o formato Shape (ESRI).
 +
 +
A transformação de coordenadas em ambos os comandos é feita com base na biblioteca PROJ.4, utilizada por grande maioria do software FOSS4G (Free Open Source Software for Geospatial).
 +
 +
== DEFINIÇÕES DE COORDENADAS CARTOGRÁFICA NO PROJ4 ==
 +
 +
Embora exista a possibilidade de utilizar para Portugal Continental as definições de coordenadas definidas na lista universal do EPSG no proj4, sugere-se que se utilizem os seguintes parâmetros específicos para as definições dessas coordenadas cartográficas:
 +
 +
* '''[D73HG] HAYFORD-GAUSS DATUM 73'''
 +
http://www.igeo.pt/produtos/Geodesia/Inf_tecnica/sistemas_referencia/Datum_73.htm
 +
 +
+proj=tmerc
 +
+lat_0=39.66666666666666
 +
+lon_0=-8.131906111111112
 +
+k=1
 +
+x_0=180.598
 +
+y_0=-86.98999999999999
 +
+ellps=intl
 +
+towgs84=-239.749,88.181,30.488,0.263,0.082,1.211,2.229
 +
 +
* '''[DLXHG] HAYFORD-GAUSS DATUM LISBOA (IGP)'''
 +
http://www.igeo.pt/produtos/Geodesia/Inf_tecnica/sistemas_referencia/Datum_Lx.htm
 +
 +
+proj=tmerc
 +
+lat_0=39.66666666666666
 +
+lon_0=-8.131906111111112
 +
+k=1
 +
+x_0=0
 +
+y_0=0
 +
+ellps=intl
 +
+towgs84=-288.885,-91.744,126.244,-1.691,0.410,-0.211,-4.598
 +
 +
* '''[DLXHGMIL] HAYFORD-GAUSS DATUM LISBOA (Militar)'''
 +
+proj=tmerc
 +
+lat_0=39.66666666666666
 +
+lon_0=-8.131906111111112
 +
+k=1
 +
+x_0=200000
 +
+y_0=300000
 +
+ellps=intl
 +
+towgs84=-288.885,-91.744,126.244,-1.691,0.410,-0.211,-4.598
 +
 +
* '''[PTTM06] PT-TM06 do datum ETRS89'''
 +
http://www.igeo.pt/produtos/Geodesia/Inf_tecnica/sistemas_referencia/Datum_ETRS89.htm
 +
 +
+proj=tmerc
 +
+lat_0=39.6682583333333
 +
+lon_0=-8.13310833333333
 +
+ellps=GRS80
 +
+towgs84=0,0,0
 +
 +
== DEFINIÇÕES DE COORDENADAS PELO METODO DAS GRELHAS (Prof. José Alberto Gonçalves) ==
 +
 +
As definições das transformações de coordenadas pelo método das grelhas obriga a utilização de um conjunto de ficheiros ([http://www.fc.up.pt/pessoas/jagoncal/coordenadas/grelhas.zip grelhas.zip]) disponibilizados pelo Prof. José Alberto Gonçalves na sua página:
 +
 +
http://www.fc.up.pt/pessoas/jagoncal/coordenadas/
 +
 +
Para aplicar esta metodologia é necessário descarregar o ficheiro [http://www.fc.up.pt/pessoas/jagoncal/coordenadas/grelhas.zip grelhas.zip] que contem a definição dos seguintes Sistemas de Coordenadas.
 +
 +
{| border="1" class="wikitable" style="margin: 1em 1em 1em 0; background: #f9f9f9; border: 1px #aaaaaa solid; border-collapse: collapse;"
 +
 +
|-
 +
 +
! bgcolor="#99CCCC" | Referencia
 +
! bgcolor="#99CCCC" | Sistema de coordenadas
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF" | etrs89
 +
| bgcolor="#FFFFFF" | Geográficas ETRS89
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |pttm06
 +
| bgcolor="#FFFFFF"  |PTTM06 do datum ETRS89
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |pttm06mil
 +
| bgcolor="#FFFFFF"  |PTTM06 do datum ETRS89 (versão militar)
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |etrs89utm
 +
| bgcolor="#FFFFFF"  |UTM-29N do datum ETRS89
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |d73geo
 +
| bgcolor="#FFFFFF"  |Geográficas Datum 73
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |d73hg
 +
| bgcolor="#FFFFFF"  |Hayford-Gauss datum 73 (IGP)
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |d73hgmil
 +
| bgcolor="#FFFFFF"  |Hayford-Gauss datum 73 (IGP) + translaccao militar
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |dlxgeo
 +
| bgcolor="#FFFFFF"  |Geográficas Datum Lisboa - Elips. Hayford
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |dlxhg
 +
| bgcolor="#FFFFFF"  |Hayford-Gauss datum Lisboa
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |dlxhgmil
 +
| bgcolor="#FFFFFF"  |Hayford-Gauss datum Lisboa (militar)
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |ed50geo
 +
| bgcolor="#FFFFFF"  |Geográficas ED50
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |ed50utm29
 +
| bgcolor="#FFFFFF"  |UTM 29-N Datum Europeu 1950
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |dlxbgeo  
 +
| bgcolor="#FFFFFF"  |Geográficas Datum Lisboa - Elipsoide de Bessel
 +
 +
|-
 +
 +
| bgcolor="#FFFFFF"  |dlxbb
 +
| bgcolor="#FFFFFF"  |Sistema Bessel-Bonne (eixos invertidos)
 +
 +
|}
 +
 +
conforme ao sistema operativo usado é necessário descompactar os ficheiros de definição na pasta com as definições padrão da biblioteca Proj, em um sistema GNU/linux esta pasta é tipicamente
 +
 +
<pre>/usr/share/proj</pre>
 +
 +
== CS2CS - CONVERSÃO DE PONTOS ==
 +
 +
A utilização deste comando é efectuada com a seguinte sintaxe simplificada (ponto a ponto):
 +
 +
<pre>cs2cs <Def. Sistema Coordenadas Dados> +to <Def. Sistema Coordenadas Resultados> [enter]
 +
 +
<x> <y> <z>[enter] (dados das coordenadas do ponto)
 +
 +
X Y Z (resultados da transformação de coordenadas)</pre>
 +
 +
Como exemplo vou converter as coordenadas do vértice Geodésico de PATAIS, cujas coordenadas oficiais são:
 +
 +
<pre># Vértice Geodésico Ordem X-HG73     Y-HG73    X-HGLX   Y-HGLX    X-TM06 Y-TM06
 +
88 PATAIAS           1 -74583,460  2462,139  -74584,120  2459,620  -74582,612 2460,852</pre>
 +
 +
=== OPÇÃO 1. Converter de DATUM 73 para PT-TM06 utilizando transformação BursaWolf de 7 parâmetros: ===
 +
 +
<pre># cs2cs +proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1
 +
+x_0=180.598 +y_0=-86.98999999999999 +ellps=intl +towgs84=-231.03,102.61,26.84,0.61,-0.20,0.88,1.79
 +
+to +proj=tmerc +lat_0=39.6682583333333 +lon_0=-8.13310833333333 +ellps=GRS80 +towgs84=0,0,0
 +
 +
# -74583.460  2462.139
 +
 +
# -74582.96  2460.61    54.41 (resultado).</pre>
 +
 +
Nota:
 +
Interessa-nos aqui apenas a transformação planimétrica pelo que deve ser ignorada a cota resultante (54.41m)
 +
 +
=== OPÇÃO 2. Converter de DATUM 73 para PT-TM06 utilizando transformação pelo método das grelhas:===
 +
 +
<pre># cs2cs +init=pt:d73hg +to +init=pt:pttm06
 +
 +
# -74583.460 2462.139
 +
 +
# -74582.61  2460.81 0.00 (resultado)</pre>
 +
 +
== OGR2OGR - CONVERSÃO DE FICHEIROS SHAPE ==
 +
 +
Tendo por exemplo como dados o ficheiro shape da CAOP versão 2008.1 fornecido pelo IGP exclusivamente no sistema de coordenadas PT-MT06. Vamos converte-lo para Datum 73 utilizado o OGR2OGR.
 +
O ficheiro shape pode ser obtido em:
 +
http://www.igeo.pt/produtos/cadastro/caop/download/CAOP2008_1_Shapes/Cont_Freg_CAOP20081_pol.zip
 +
 +
=== OPÇÃO 1. Converter de PT-TM06 para DATUM 73 utilizando transformação BursaWolf de 7 parâmetros: ===
 +
 +
<pre># ogr2ogr -f "ESRI Shapefile" -s_srs "+proj=tmerc +ellps=GRS80 +towgs84=0,0,0
 +
+lat_0=39.6682583333333 +lon_0=-8.1331083333333" -t_srs "+proj=tmerc +lat_0=39.66666666666666
 +
+lon_0=-8.131906111111112 +k=1 +x_0=180.598 +y_0=-86.98999999999999 +ellps=intl
 +
+towgs84=-239.749,88.181,30.488,0.263,0.082,1.211,2.229" Cont_Freg_CAOP20081_pol_d73hg_bw7.shp Cont_Freg_CAOP20081_pol.shp</pre>
 +
 +
E resulta numa shape file (4 ficheiros) com nome: “Cont_Freg_CAOP20081_pol_d73hg_bw7”
 +
 +
=== OPÇÃO 2. Converter de PT-TM06 para DATUM 73 utilizando transformação pelo método das grelhas: ===
 +
 +
<pre># ogr2ogr -f "ESRI Shapefile" -s_srs "+init=pt:pttm06 +wktext" -t_srs "+init=pt:d73hg +wktext" Cont_Freg_CAOP20081_pol_d73hg_grid.shp Cont_Freg_CAOP20081_pol.shp</pre>
 +
 +
E resulta numa shape file (4 ficheiros) com nome: “Cont_Freg_CAOP20081_pol_d73hg_grid”
  
 
[[Category:Portugal]]
 
[[Category:Portugal]]

Revision as of 12:23, 6 February 2009

Transformação de coordenadas cartográficas de e para o novo sistema de projecção oficial de portugal continental: o pt-tm06

Fonte: Ricardo Pinho (http://gisvm.com/forum/index.php?topic=38.0)

  • PT-MT06, baseado no ETRS89 e indicado pelo IGP para novo Sistema de Coordenadas Cartográficas Oficial Português:

http://www.igeo.pt/produtos/Geodesia/Inf_tecnica/sistemas_referencia/Datum_ETRS89.htm

  • Página do Prof. José Alberto Gonçalves, da Faculdade de Ciências da Universidade do Porto, que propõe a transformação de coordenadas pelo método das Grelhas.

http://www.fc.up.pt/pessoas/jagoncal/coordenadas/

  • Ferramentas utilizadas: CS2CS e OGR2OGR, que se podem encontrar no pacote fwtools ou directamente instaladas com as bibliotecas Proj.4 e OGR

DESCRIÇÃO DOS UTILITÁRIOS CS2CS e OGR2OGR

  • CS2CS, é um comando para usar em linhas de comandos (terminal/consola) que permite a conversão de coordenadas de um ponto ou conjunto de pontos (num ficheiro) entre dois sistemas de coordenadas indicados (Coordinate System).
  • OGR2OGR, é um comando para usar em linhas de comandos (terminal/consola) que permite a conversão de uma diversidade grande de formatos vectoriais abrangidos pela biblioteca OGR. Vamos aqui detalhar a possibilidade de o utilizar para conversão de coordenadas no mesmo formato, nomeadamente o formato Shape (ESRI).

A transformação de coordenadas em ambos os comandos é feita com base na biblioteca PROJ.4, utilizada por grande maioria do software FOSS4G (Free Open Source Software for Geospatial).

DEFINIÇÕES DE COORDENADAS CARTOGRÁFICA NO PROJ4

Embora exista a possibilidade de utilizar para Portugal Continental as definições de coordenadas definidas na lista universal do EPSG no proj4, sugere-se que se utilizem os seguintes parâmetros específicos para as definições dessas coordenadas cartográficas:

  • [D73HG] HAYFORD-GAUSS DATUM 73

http://www.igeo.pt/produtos/Geodesia/Inf_tecnica/sistemas_referencia/Datum_73.htm

+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=180.598 +y_0=-86.98999999999999 +ellps=intl +towgs84=-239.749,88.181,30.488,0.263,0.082,1.211,2.229

  • [DLXHG] HAYFORD-GAUSS DATUM LISBOA (IGP)

http://www.igeo.pt/produtos/Geodesia/Inf_tecnica/sistemas_referencia/Datum_Lx.htm

+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +towgs84=-288.885,-91.744,126.244,-1.691,0.410,-0.211,-4.598

  • [DLXHGMIL] HAYFORD-GAUSS DATUM LISBOA (Militar)

+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=200000 +y_0=300000 +ellps=intl +towgs84=-288.885,-91.744,126.244,-1.691,0.410,-0.211,-4.598

  • [PTTM06] PT-TM06 do datum ETRS89

http://www.igeo.pt/produtos/Geodesia/Inf_tecnica/sistemas_referencia/Datum_ETRS89.htm

+proj=tmerc +lat_0=39.6682583333333 +lon_0=-8.13310833333333 +ellps=GRS80 +towgs84=0,0,0

DEFINIÇÕES DE COORDENADAS PELO METODO DAS GRELHAS (Prof. José Alberto Gonçalves)

As definições das transformações de coordenadas pelo método das grelhas obriga a utilização de um conjunto de ficheiros (grelhas.zip) disponibilizados pelo Prof. José Alberto Gonçalves na sua página:

http://www.fc.up.pt/pessoas/jagoncal/coordenadas/

Para aplicar esta metodologia é necessário descarregar o ficheiro grelhas.zip que contem a definição dos seguintes Sistemas de Coordenadas.

Referencia Sistema de coordenadas
etrs89 Geográficas ETRS89
pttm06 PTTM06 do datum ETRS89
pttm06mil PTTM06 do datum ETRS89 (versão militar)
etrs89utm UTM-29N do datum ETRS89
d73geo Geográficas Datum 73
d73hg Hayford-Gauss datum 73 (IGP)
d73hgmil Hayford-Gauss datum 73 (IGP) + translaccao militar
dlxgeo Geográficas Datum Lisboa - Elips. Hayford
dlxhg Hayford-Gauss datum Lisboa
dlxhgmil Hayford-Gauss datum Lisboa (militar)
ed50geo Geográficas ED50
ed50utm29 UTM 29-N Datum Europeu 1950
dlxbgeo Geográficas Datum Lisboa - Elipsoide de Bessel
dlxbb Sistema Bessel-Bonne (eixos invertidos)

conforme ao sistema operativo usado é necessário descompactar os ficheiros de definição na pasta com as definições padrão da biblioteca Proj, em um sistema GNU/linux esta pasta é tipicamente 

/usr/share/proj

CS2CS - CONVERSÃO DE PONTOS

A utilização deste comando é efectuada com a seguinte sintaxe simplificada (ponto a ponto): 

cs2cs <Def. Sistema Coordenadas Dados> +to <Def. Sistema Coordenadas Resultados> [enter]

<x> <y> <z>[enter] (dados das coordenadas do ponto)

X Y Z (resultados da transformação de coordenadas)

Como exemplo vou converter as coordenadas do vértice Geodésico de PATAIS, cujas coordenadas oficiais são: 

#	Vértice Geodésico Ordem	X-HG73	    Y-HG73    X-HGLX	  Y-HGLX    X-TM06	Y-TM06
88	PATAIAS	          1	-74583,460  2462,139  -74584,120  2459,620  -74582,612	2460,852

OPÇÃO 1. Converter de DATUM 73 para PT-TM06 utilizando transformação BursaWolf de 7 parâmetros:

# cs2cs +proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 
+x_0=180.598 +y_0=-86.98999999999999 +ellps=intl +towgs84=-231.03,102.61,26.84,0.61,-0.20,0.88,1.79 
+to +proj=tmerc +lat_0=39.6682583333333 +lon_0=-8.13310833333333 +ellps=GRS80 +towgs84=0,0,0

# -74583.460  2462.139

# -74582.96   2460.61    54.41 (resultado).

Nota: Interessa-nos aqui apenas a transformação planimétrica pelo que deve ser ignorada a cota resultante (54.41m)

OPÇÃO 2. Converter de DATUM 73 para PT-TM06 utilizando transformação pelo método das grelhas:

# cs2cs +init=pt:d73hg +to +init=pt:pttm06

# -74583.460 2462.139

# -74582.61   2460.81 0.00 (resultado)

OGR2OGR - CONVERSÃO DE FICHEIROS SHAPE

Tendo por exemplo como dados o ficheiro shape da CAOP versão 2008.1 fornecido pelo IGP exclusivamente no sistema de coordenadas PT-MT06. Vamos converte-lo para Datum 73 utilizado o OGR2OGR. O ficheiro shape pode ser obtido em: http://www.igeo.pt/produtos/cadastro/caop/download/CAOP2008_1_Shapes/Cont_Freg_CAOP20081_pol.zip

OPÇÃO 1. Converter de PT-TM06 para DATUM 73 utilizando transformação BursaWolf de 7 parâmetros:

# ogr2ogr -f "ESRI Shapefile" -s_srs "+proj=tmerc +ellps=GRS80 +towgs84=0,0,0 
+lat_0=39.6682583333333 +lon_0=-8.1331083333333" -t_srs "+proj=tmerc +lat_0=39.66666666666666 
+lon_0=-8.131906111111112 +k=1 +x_0=180.598 +y_0=-86.98999999999999 +ellps=intl 
+towgs84=-239.749,88.181,30.488,0.263,0.082,1.211,2.229" Cont_Freg_CAOP20081_pol_d73hg_bw7.shp Cont_Freg_CAOP20081_pol.shp

E resulta numa shape file (4 ficheiros) com nome: “Cont_Freg_CAOP20081_pol_d73hg_bw7”

OPÇÃO 2. Converter de PT-TM06 para DATUM 73 utilizando transformação pelo método das grelhas:

# ogr2ogr -f "ESRI Shapefile" -s_srs "+init=pt:pttm06 +wktext" -t_srs "+init=pt:d73hg +wktext" Cont_Freg_CAOP20081_pol_d73hg_grid.shp Cont_Freg_CAOP20081_pol.shp

E resulta numa shape file (4 ficheiros) com nome: “Cont_Freg_CAOP20081_pol_d73hg_grid”

Retrieved from "http://wiki.osgeo.org/w/index.php?title=Transformação_de_coordenadas_cartográficas&oldid=35056"