QGIS & Transformação de Datum (parâmetros oficiais)
Danielle Cintra
Hebert Guilherme de Azevedo
A reprojecao on the fly deve estar ligada.
--
Comunidade QGIS Brasil: http://qgisbrasil.org/
---
Você recebeu essa mensagem porque está inscrito no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para postar nesse grupo, envie um e-mail para qgisb...@googlegroups.com.
Acesse esse grupo em https://groups.google.com/group/qgisbrasil.
Para mais opções, acesse https://groups.google.com/d/optout.
yjmenezes
Bom dia.
Concordo que eh um topico importante e que nao tem merecido a atencao.
a) Os deslocamentos podem atingir algo em torno de 60 metros, dependendo da regiao do BraSil.
Tem um agravante, as diferentes realizacoes do SRC. Temos de cuidar com as datas para saber qual realizacao do SAD69 ou Corrego Alegre.
b) nunca usei ArqGIS mas eu creio que ele use os parametros da EPSG, assim como o QGIS. Ja comentei aqui na lista sobre erros no EPSG para o BR.
c) Veja esta monografia do IBGE, pode testar no ArqGIS e confirmar se ele faz corretamente.
d) eu desaconselho o o uso de SRC personalizado, a nao ser em casos muito especias, como Curitiba, que passui um sistema proprio, parecido com SAD69.
Parametros locais de Curitiba para uso na transformacao de Molodensky com ./xdatum CuritibaSAD69 para SIRGAS2000 dx=-66.163 dy=2.028 dz=-33.718 SIRGAS2000 para CuritibaSAD69 dx=66.163 dy=-2.028 dz=33.718
http://cartognu.org/geoshell/2014-12/
SAT91242.pdf
[]s
julio menezes
cartognu.org
Felipe Dantas
Danielle Cintra
--
Comunidade QGIS Brasil: http://qgisbrasil.org/
---
Você recebeu essa mensagem porque está inscrito no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para postar nesse grupo, envie um e-mail para qgisb...@googlegroups.com.
Acesse esse grupo em https://groups.google.com/group/qgisbrasil.
Para mais opções, acesse https://groups.google.com/d/optout.
Processamento Digital
Analista de Geoprocessamento
----------------------------------------------
yjmenezes
Danielle
app1: esta informacao deveria ser o proprio EPSG, com codigos e parametros especificos, e corretos, para cada realizacao.
tem ainda a questao dos modelos matematicos. Molodensky, NTv2
http://www.ibge.gov.br/english/geociencias/geodesia/param_transf/default_param_transf.shtm
app2: porque a persolalizacao do SRC tera de acompanhar os dados na distribuicao destes ou havera deslocamento dos mesmos.
em minha opniao, o IBGE deveria, como orgao oficial, fornecer as correcoes ao EPSG.
A transformacao de dados antigos em SIRGAS 2000,4 tem de ser feita com muito criterio, especialmente para escalas pequenas.
Saber claramente o SRC dos dados para adotar os parametros corretos.
[]s
julio menezes
cartognu.org
Carlos P
Bom dia pessoal,
Dando continuidade sobre as dúvidas apresentadas neste tópico...
Estou tendo várias dificuldades com relação a transformação de data SAD69 para SIRGAS2000 e SIRGAS2000 para SAD69.
Na semana passada encontrei diversos relatos sobre os problemas que o QGIS tem por não efetuar a correta transformação em função de não conter os parâmetros corretos conforme a resolução do IBGE, n° 01/2005.
Vi estes links:
http://qgisbrasil.org/blog/2015/02/25/qgis-e-transformacao-de-datum-parametros-oficiais-diagnostico/
http://trac.osgeo.org/proj/ticket/241
Aí, fala sobre os erros apresentados e sugestões de como corrigi-los, como a criação de um novo SRC para corrigir os parâmetros que estão errados ou não estão listados no QGIS. Nesta última versão, 2.14, até já estão listados alguns parâmetros (diferente das versões anteriores, onde os mesmos não eram nem listados). Porém, percebi que apenas está listado o parâmetro de transformação de SAD 69 para SIRGAS2000 e não o contrário. Como o que muda neste caso, é apenas o sinal, saberiam me dizer se é necessário criar um SRC personalizado para resolver isto ou o QGIS efetua isso de forma automática?
Outra questão, seria muito importante que a comunidade QGIS atentasse para estas incorreções, já que podem causar um grande transtorno, visto que é um erro do programa. Além disso, causa uma impressão negativa principalmente para novos usuários. Porém, como isto é uma questão que envolve conhecimentos de programação, deve chegar até quem tem condições de corrigir estes problemas.
Se alguém puder me ajudar, fico muito agradecido.
Abraços
Saberia me ajudar?
Muito obrigado, desde já, pela ajuda!
Abraços
Em sexta-feira, 24 de junho de 2016 07:57:42 UTC-3, Danielle Cintra escreveu:
yjmenezes
Em segunda-feira, 4 de julho de 2016 12:04:46 UTC-3, Carlos P escreveu:
Bom dia pessoal,
Dando continuidade sobre as dúvidas apresentadas neste tópico...
Estou tendo várias dificuldades com relação a transformação de data SAD69 para SIRGAS2000 e SIRGAS2000 para SAD69.
Na semana passada encontrei diversos relatos sobre os problemas que o QGIS tem por não efetuar a correta transformação em função de não conter os parâmetros corretos conforme a resolução do IBGE, n° 01/2005.
Aí, fala sobre os erros apresentados e sugestões de como corrigi-los, como a criação de um novo SRC para corrigir os parâmetros que estão errados ou não estão listados no QGIS. Nesta última versão, 2.14, até já estão listados alguns parâmetros (diferente das versões anteriores, onde os mesmos não eram nem listados). Porém, percebi que apenas está listado o parâmetro de transformação de SAD 69 para SIRGAS2000 e não o contrário. Como o que muda neste caso, é apenas o sinal, saberiam me dizer se é necessário criar um SRC personalizado para resolver isto ou o QGIS efetua isso de forma automática?
Outra questão, seria muito importante que a comunidade QGIS atentasse para estas incorreções, já que podem causar um grande transtorno, visto que é um erro do programa. Além disso, causa uma impressão negativa principalmente para novos usuários. Porém, como isto é uma questão que envolve conhecimentos de programação, deve chegar até quem tem condições de corrigir estes problemas.
Se alguém puder me ajudar, fico muito agradecido.
Abraços
Carlos,
O problema nao esta no QGIS, mas sim os arquivos que contem os parametros, que sao oriundos do: hhttp://www.epsg.org/.
Se usa algum outro GIS, transforme esses dados, (monografia em anexo), do IBGE e repasse o resultado.
Por favor, reporte o SW e os resultados.
[]s
julio menezes
PS: existe uma maxima em programacao: "garbage in garbage out". Dados ruins geram resultados ruins.
cartognu.org
Processamento Digital
Analista de Geoprocessamento
----------------------------------------------
--
Carlos P
Você recebeu essa mensagem porque está inscrito em um tópico no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse tópico, acesse https://groups.google.com/d/topic/qgisbrasil/SKocT3c1WqQ/unsubscribe.
Para cancelar inscrição nesse grupo e todos os seus tópicos, envie um e-mail para qgisbrasil+...@googlegroups.com.
yjmenezes
Anexo
arquivo baixado do EPSG. Versao 8.9
http://www.epsg.org/EPSGDataset/DownloadDataset.aspx
de uma olhada.
[]s
julio menezes
Hebert Guilherme de Azevedo
TRANSFORMAÇÃO DE COORDENADAS
11. Existem parâmetros de transformação entre WGS 84 e SIRGAS2000?
Atualmente não existem parâmetros de transformação entre SIRGAS2000 e WGS 84 porque eles são praticamente iguais, ou seja, DX = 0, DY = 0 e DZ = 0.
Desde o estabelecimento do sistema GPS (Global Positioning System), o seu Sistema Geodésico de Referência (WGS 84) já passou por quatro refinamentos. Nestas quatro atualizações o objetivo sempre foi aproximá-lo ao ITRF (International Terrestrial Reference Frame), materialização mais precisa do ITRS (International Terrestrial Reference System), desenvolvida pelo IERS (International Earth Rotation and Reference Systems Service). A mais recente atualização recebeu a denominação de WGS 84 (G1674), adotado no Sistema GPS a partir de 08 de fevereiro de 2012. Os parâmetros de transformação WGS 84/SAD 69, divulgados através da Resolução do Presidente do IBGE n° 23, de 21/02/89 (R.PR 23/89), permanecem válidos para transformar coordenadas determinadas por posicionamentos GPS realizados no período de 01/01/1987 a 01/01/1994 - quando a versão correspondente do WGS 84 se denominava WGS 84 (Doppler).
Parâmetros WGS 84 (Doppler) para SAD69:
DX = +66,87 m
DY = -4,37 m
DZ = +38,52 m
Os parâmetros SAD 69/SIRGAS2000 utilizados no ProGriD (opção: SAD 69 Técnica Doppler ou GPS) e divulgados através da Resolução do Presidente do IBGE n° 1, de 25/02/2005 (R.PR 01/05), são válidos para transformar coordenadas entre SAD 69/WGS 84 e SAD 69/SIRGAS2000 determinadas por posicionamentos GNSS realizados após 01/01/1994.
SAD 69 para SIRGAS2000 (≡ WGS 84 (G1150)):
DX = -67,35 m
DY = +3,88 m
DZ = -38,22 m
E posteriormente faço a transformação normalmente. Existe algum erro em adotar este procedimento?
Att
(31) 7531-2575
Skype: hebert.guilherme.de.azevedo
lattes: http://lattes.cnpq.br/6582724988395217
yjmenezes
fui dar ima olhada e vejam o resultado:
jmenezes@gnu001:/tmp$ cat -n /usr/share/proj/epsg | awk 'NR >7089 && NR < 7105 { print }'
7090 <29187> +proj=utm +zone=17 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs <>
7091 # SAD69 / UTM zone 18S
7092 <29188> +proj=utm +zone=18 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs <>
7093 # SAD69 / UTM zone 19S
7094 <29189> +proj=utm +zone=19 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs <>
7095 # SAD69 / UTM zone 20S
7096 <29190> +proj=utm +zone=20 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs <>
7097 # SAD69 / UTM zone 21S
7098 <29191> +proj=utm +zone=21 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs <>
7099 # SAD69 / UTM zone 22S
7100 <29192> +proj=utm +zone=22 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs <>
7101 # SAD69 / UTM zone 23S
7102 <29193> +proj=utm +zone=23 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs <>
7103 # SAD69 / UTM zone 24S
7104 <29194> +proj=utm +zone=24 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs <>
O problema continua, erros nos dados do grupo EPSG e nao no QGIS ou PROJ4.
[]s
julio menezes
cartognu.org
Hebert Guilherme de Azevedo
yjmenezes
segue o conteudo de arquivo epsg que acompanha o QGIS e 2 screenshots.
So lembro que o problema nao esta nas implementacoes do Frank W mas sim nos arquivos da EPSG.
[]s
julio menezes
cartognu.org
yjmenezes
QGIS2.14.3 vs QGIS2.16
a porcaria dos dados dos petroleiros EPSG.
No 2.14 esta ok, ja no 2.16 :(
Tenho esperanca que algum dia o IBGE informe ao EPSG os valores a serem usados nesta Torre de Babel.
[]s
julio menezes
Em segunda-feira, 11 de julho de 2016 17:34:16 UTC-3, Hebert Guilherme de Azevedo escreveu:
Setor Geo
Pessoal, me tirem uma dúvida,
O que contém os parâmetros de transformação é a EPSG? Esta não é apenas um sistema de codificação? Os parâmetros não estão na biblioteca, PROJ4, que por sua vez é utilizada pelo QGis?
Me esclareçam por favor, estou en dúvida.
Se não for isso, me expliquem o que é mais a EPSG.
Obrigado
yjmenezes
Os arquivos do EPSG contem parametros, inclusive de DATUNs locais da Petrobras, como Aratu, por exemplo.
A ESRI tambem tem seus parametros.
Espero que este projeto ponha ordem nesta baderna.
https://trac.osgeo.org/metacrs/wiki/WikiStart
@Guilherme,
jmenezes@gnu001:/tmp$ gdalsrsinfo "EPSG:5533"
PROJ.4 : '+proj=utm +zone=23 +south +ellps=aust_SA +towgs84=-67.35,3.88,-38.22,0,0,0,0 +units=m +no_defs '
veja:
AUTHORITY["EPSG","5533"]]
[]s
julio menezes
cartognu.org
PROJ.4 : '+proj=utm +zone=23 +south +ellps=aust_SA +towgs84=-67.35,3.88,-38.22,0,0,0,0 +units=m +no_defs '
OGC WKT :
PROJCS["SAD69(96) / UTM zone 23S",
GEOGCS["SAD69(96)",
DATUM["South_American_Datum_1969_96",
SPHEROID["GRS 1967 Modified",6378160,298.25,
AUTHORITY["EPSG","7050"]],
TOWGS84[-67.35,3.88,-38.22,0,0,0,0],
AUTHORITY["EPSG","1075"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","5527"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-45],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",10000000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","5533"]]
yjmenezes
hoje entrei no wiki do Proj4.
O Frank esta atento, ha um link apontando para o IBGE
https://github.com/OSGeo/proj.4/wiki
Ha referencias a grids NTv2 para duas realizacoes do Corrego Alegre e duas realizacoes do SAD69.
Brazilian grids for datums Corrego Alegre 1961, Corrego Alegre 1970-72, SAD69 and SAD69(96)
A mapping file (
epsg) has also been introduced mapping most EPSG coordinate systems to PROJ.4 format.VerticalDatums. Lembrando que da mesma forma que a rede horizontal, a rede de nivelamento do IBGE sofreu reajustamentos e uma mesma RN fisica recebeu novas altitudes.
Isto tem consequencias importantes em projetos de hidroeletricas e estudos hidraulicos. ( curva de inundacao, remanso... )[]s
julio menezes
cartognu.org
Germano Pirajá Martins
Ezequiel Faria
Dificuldades nos parâmetros de transformações:
Quando transformo de córrego para sad69 utilizando src personalizada (com os parâmetros ibge), observei que o shp cai no mesmo lugar sad 69(96) epsg 5533 caso eu fizesse esta transformação de córrego para sad69(96) epsg 5533 sem parâmetros personalizados ibge.
Quando transformo de córrego para sad69 epsg 29193 com parâmetros personalizados do ibge no arcgis, observei que o novo shp cai onde estão os arquivos sad69 epsg 29193 retirado do qgis sem src personalizado (apenas salvando com nova epsg um arquivo original em córrego para sad69 epgs 29193).
Resultado: quando pego o shp proveniente de uma transformação personalizada feita no qgis (com os parâmetros ibge), cujo o shp personalizado sai com o mesmo posicionamento da epsg 5533 de sad69(96) já existente no qgis. E posteriormente pego um shp proveniente de uma transformação personalizada córrego para sad69 (96) com epsg 5533 (com os parâmetros ibge) feita no arcgis e coloco os dois shp juntos com parâmetros personalizados pelo ibge, ambos não caem no mesmo lugar.
No final das contas fazendo uma transformação de shp utilizando parâmetros ibge de córrego para sad69 os shp não estão batendo, quando coloco arquivos no arcgis ou no qgis.
Notas:
Valores dos parâmetros de transformações SAD69 (material em anexo) coincide com epsg 5533 (no caso qgis)
Fontes:
Parâmetros do qgis:
Parâmetros do Arcgis:
Em sexta-feira, 24 de junho de 2016 07:57:42 UTC-3, Danielle Cintra escreveu:
Aldisio Carlos
yjmenezes
**** SAD69 classico:
jmenezes@gnu000:~$ gdalsrsinfo epsg:29195 -o proj4
'+proj=utm +zone=25 +south +ellps=aust_SA +towgs84=-66.87,4.37,-38.52,0,0,0,0 +units=m +no_defs '
**** SAD69 ajuste 1996:
jmenezes@gnu000:~$ gdalsrsinfo epsg:5533 -o proj4
'+proj=utm +zone=23 +south +ellps=aust_SA +towgs84=-67.35,3.88,-38.22,0,0,0,0 +units=m +no_defs '
julio
Ezequiel Faria
yjmenezes
Nao eh apenas uma questao de datas, mas da origem usada.
Tem de saber exatamente o referencial usado onde foram feitas as amarracoes, para poder transformar.
A outra alternativa disponivel no Progrid eh um grid, NTv2, onde as correcoes sao modeladas em funcao de desvios locais.
Se nao tiver informacao sobre o referencial das amaraccoes, eh tentar validar rastreando alguns pontos e amarrando a rede RBMC ou a vertices do IBGE e verificar qual das transformacoes se adapta melhor.
*** importante, a rede altimetrica do IBGE tambem foi reajustada.
GPS de navegacao, *todos que conheco, tem o SAD69 com parametros errados.
WGS<-->SIRGAS2000, para efeitos de engenharia, podem ser considerados iguais.
Ezequiel Faria
Tranformações pelo QGIS:
Sirgas 2000 para SAD69 (epgs 4618):
Extent: -48.00, -33.50, -42.00, 5.13
Proj4: +proj=utm +zone=23 +south +ellps=aust_SA +towgs84=-66.87,4.37,-38.52,0,0,0,0 +units=m +no_defs
Sirgas 2000 para SAD 69 (epsg 4291)
Extent: -81.41, -55.96, -25.28, 12.52
Proj4: +proj=longlat +ellps=GRS67 +towgs84=-57,1,-41,0,0,0,0 +no_defs
Sirgas 2000 para SAD69(96) (epgs 5527):
Extent: -48.00, -33.50, -42.00, 5.13
Proj4: +proj=utm +zone=23 +south +ellps=aust_SA +towgs84=-67.35,3.88,-38.22,0,0,0,0 +units=m +no_defs
Tranformações pelo ArcGis:
Sirgas 2000 para SAD69 (epgs 4618)
Paramentro 2 e o Progrid são os mesmos (NTv2 Sad69_003)
Paramentro 1 (dx 67,35 dy 3,88 e dz 38,22)
Sirgas 2000 para SAD69(96) (epgs 5527):
Paramentro 1 e o Progrid são os mesmos (NTv2 Sad96_003)
Paramentro 2 (dx 67,35 dy 3,88 e dz 38,22)
Resultado das transformações entre Sirgas e Sad pelo QGIS e ArcGis:
- Sirgas 2000 to Sad69(96) parâmetro 2 Arcgis, Sirgas 2000 to Sad69(96) Qgis, Sirgas 2000 to Sad69 parâmetro 1 Arcgis, são coincidentes.
- Sirgas 2000 to Sad69(96) parâmetro 1 Arcgis, Sirgas 2000 to Sad69(96) parâmetro progrid Arcgis, são coincidentes.
- Sirgas 2000 to Sad69 parâmetro 2 Arcgis, Sirgas 2000 to Sad69 parâmetro progrid Arcgis, são coincidentes.
- Sirgas 2000 to Sad69 Qgis, não coincide com nenhum
Fonte: https://wiki.osgeo.org/wiki/Brazilian_Coordinate_Reference_Systems
Fonte: http://qgisbrasil.org/blog/2015/02/25/qgis-e-transformacao-de-datum-parametros-oficiais-diagnostico/
11. Existem parâmetros de transformação entre WGS 84 e SIRGAS2000?
Atualmente
não existem parâmetros de transformação entre SIRGAS2000 e WGS 84 porque eles
são praticamente iguais, ou seja, DX = 0, DY = 0 e DZ = 0.
Desde o estabelecimento do sistema GPS (Global Positioning System),
o seu Sistema Geodésico de Referência (WGS 84) já passou por quatro
refinamentos. Nestas quatro atualizações o objetivo sempre foi aproximá-lo ao
ITRF (International Terrestrial Reference Frame),
materialização mais precisa do ITRS (International Terrestrial Reference
System), desenvolvida pelo IERS (International Earth Rotation
and Reference Systems Service). A mais recente atualização recebeu a
denominação de WGS 84 (G1674), adotado no Sistema GPS a partir de 08 de
fevereiro de 2012. Os parâmetros de transformação WGS 84/SAD 69, divulgados
através da Resolução do Presidente do IBGE n° 23, de 21/02/89 (R.PR 23/89),
permanecem válidos para transformar coordenadas determinadas por
posicionamentos GPS realizados no período de 01/01/1987 a 01/01/1994 - quando a versão
correspondente do WGS 84 se denominava WGS 84 (Doppler).
Parâmetros
WGS 84 (Doppler) para SAD69:
DX = +66,87 m
DY = -4,37 m
DZ = +38,52 m
Os parâmetros SAD 69/SIRGAS2000 utilizados no ProGriD (opção: SAD 69 Técnica Doppler ou GPS) e divulgados através da Resolução do Presidente do IBGE n° 1, de 25/02/2005 (R.PR 01/05), são válidos para transformar coordenadas entre SAD 69/WGS 84 e SAD 69/SIRGAS2000 determinadas por posicionamentos GNSS realizados após 01/01/1994.
SAD
69 para SIRGAS2000 (≡ WGS 84 (G1150)):
DX = -67,35 m
DY = +3,88 m
DZ = -38,22 m
Questionamento:
Após ver as transformações e alguns tutoriais na internet de paramentros transformações no QGIS e ArcGis (http://www.processamentodigital.com.br/wp-content/uploads/2015/06/ArcGIS103_Reprojecao_de_SAD1969_para_SIRGAS20001.pdf, https://geobases.es.gov.br/Media/Geobases/documentos/manuais/032_SAD69_SIRGAS2000_ArcGis_NTV2.pdf), que uso chego a algumas duvidas:
- Transformação de Sad 69 (epsg 4618) para Sirgas 2000, em ambos tutoriais utilizam os valores (dx 67,35 dy 3,88 e dz 38,22) ou o arquivo Sad69_003.gsb (progrid). Ao consultar está tabela:
Observa que o parâmetro para o Sad 69 (epsg 4618) para Sirgas 2000 já é outro arquivo Sad69_003.gsb.
Outra observação seria em relação aos datum e suas epsg:
- no software ArcGis quando digitamos a epsg 4291 ele submete a epsg 4618.
- se olharmos os valores da epsg 4618 do QGis (dx -66.87, dy 4.37,dz -38.52) não coincide com os valores da epsg 4618 do ArcGis.
- os valores da epsg 5527 do Qgis (dx -67.35, dy 3.88, dz -38.22) coincide que com os valores da epsg 5527 e 4618 do ArcGis.
Referencias:
http://qgisbrasil.org/blog/2015/02/25/qgis-e-transformacao-de-datum-parametros-oficiais-diagnostico/
https://geobases.es.gov.br/Media/Geobases/documentos/manuais/032_SAD69_SIRGAS2000_ArcGis_NTV2.pdf
https://wiki.osgeo.org/wiki/Brazilian_Coordinate_Reference_Systems