datum e sumarize
Cássio Marcelo Castro
1) nos datuns relacionados pelo programa (SRC) qual a diferença entre a SAD 69 e SAD 69 (deprecated)? seria sad 69 e sad 69/1996?
2) como realizar a função SUMARIZE, como há no arcgis?
3) no layout, o que signfica EXPORTAR como RASTER, no item GERAL?
Agradeço
Saudações
--
Cássio Castro
Urbanista
Salvador, BA BRASIL
yjmenezes
Observando e consultando, em meu caso, um projeto em SAD96, meridiano
-57 , zona 21S
29191 ou 29180, ambos apresentam -57. 1, -41 e mais 4 parametros
zerados trocando apenas elipsoide.
+proj=utm +zone=21 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0
+units=m +no_defs
# WGS-84 -> SAD-69 original, classico, parametros do IBGE
6378137.000 298.257223563 -> 6378160.000 298.250000000 dx 66.870 dy
-4.370 dz 38.520
ja para o ajuste de 1996, parametros do IBGE
# WGS-84 ->SAD-69/1996 Chua (UGGI-1967) -57 00'00.0000"
6378137.000 298.257223563 -> 6378160.000 298.250000000 dx 67.350 dy
-3.880 dz 38.220
inverte-se o sinal para a transformacao reversa.
Os 7 parametros constantes no QGIS acredito serem DX,DY,DZ, rotacoes
e ppm de escala.
-57, 1, -41, 0, 0,0, 0
Estes valores *nao* batem com os divulgados pelo IBGE para as formulas
simplificadas de Molodensky.
Eu trabalho muito GPS geodesico e referencial correto eh
fundamental.
Tambem estou intrigado e tenho de verificar isso pois pode comprometer
os dados.
Nao sei se o QGIS usa o modelo de Molodensky, de qq modo ha de se
verificar a transformacao entre sistemas.
julio
Aproveitando:
Devo postar algo sobre os SIGs e a OGC verem o mundo em 2D e sem SRID
( shapefile, WKT )
O PostGIS implementa o EWKT que adciona a terceira dimensao e mais o
SRID.
Nao postei ainda pois estou comecando com o QGIS e *nao* quero postar
sem antes dar uma estudada.
Agradeco se tiver alguma observacao a este respeito.
On Nov 28, 1:06 am, Cássio Marcelo Castro <cassiomarc...@gmail.com>
wrote:
Cássio Marcelo Castro
Obrigado pelas observações
--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
yjmenezes
wrote:
> Valeu
> Obrigado pelas observações
>
> Cassio,
Boa tarde.
Fiz uma alteracao num programa meu aqui para ver os parametros SAD x
WGS do QuantumGIS.
Provavelmente o QuantumGIS use as excelentes bibliotecas OGROGR de
Frank Warmerdam mas, acredito, usem parametros para a America Latina.
Tendo tempo vou ver nos programa fonte o modelo matematico e
comentarios sobre a transformacao. Acesso ao codigo fonte,
maravilha :-)
Esta meio baguncado aqui pois fui fazendo e copiando saida de
console.
Ve se da para entender pois ha diferencas significativas:
Eu ate gostaria de saber o comportamento num software privativo como
ArqGIS e outros.
[]s julio
# Coordenadas retiradas da monografia oficial do IBGE. SAT90009
# 90009U96IBGE 615273.779 8797212.968 176.88 U96HE20S -63 IBGE
SAD96
# 90009U2KIBGE 615215.609 8797173.383 172.62 U2KHE20S -63 IBGE
SIRGAS2000
1) transformando de SAD96 para WGS84
jmenezes@gnu-00:~/Desktop/wutm$ ./utmmudatum -m -63 -s SAD96 -d WGS84
< sat90009u96he.txt
# WGS-84 -63 00'00.0000" altitude elipsoidal
# SAD-69/1996 Chua (UGGI-1967) -> WGS-84
# 6378160.000 298.250000000 -> 6378137.000 298.257223563 dx -67.350
dy 3.880 dz -38.220
#Pto xE yN altitude elipsoidal
90009 615215.615 8797173.388 172.603
2) transformando de SAD96 para SIRGAS2000
jmenezes@gnu-00:~/Desktop/wutm$ ./utmmudatum -m -63 -s SAD96 -d
SIRGAS2000 < sat90009u96he.txt
# SIRGAS-2000 (GRS-1980) -63 00'00.0000" altitude elipsoidal
# SAD-69/1996 Chua (UGGI-1967) -> SIRGAS-2000 (GRS-1980)
# 6378160.000 298.250000000 -> 6378137.000 298.257222101 dx -67.350
dy 3.880 dz -38.220
#Pto xE yN altitude elipsoidal
90009 615215.615 8797173.388 172.603
3) transformando de SIRGAS2000 para SAD96
jmenezes@gnu-00:~/Desktop/wutm$ ./utmmudatum -m -63 -s SIRGAS2000 -d
SAD96 < sat90009u2khe.txt
# SAD-69/1996 Chua (UGGI-1967) -63 00'00.0000" altitude elipsoidal
# SIRGAS-2000 (GRS-1980) -> SAD-69/1996 Chua (UGGI-1967)
# 6378137.000 298.257222101 -> 6378160.000 298.250000000 dx 67.350
dy -3.880 dz 38.220
#Pto xE yN altitude elipsoidal
90009 615273.774 8797212.963 176.896
4)-------- Agora usando os parametros do QuantumGIS e transformado de
SAD_OGC->WGS
jmenezes@gnu-00:~/Desktop/wutm$ ./utmmudatum -m -63 -s SADOGC -d WGS84
< sat90009u96he.txt
# WGS-84 -63 00'00.0000" altitude elipsoidal
# SAD-69/OGC Chua (UGGI-1967) -> WGS-84
# 6378160.000 298.250000000 -> 6378137.000 298.257223563 dx -57.000
dy 1.000 dz -41.000
#Pto xE yN altitude elipsoidal
90009 615223.387 8797172.030 180.404
5) voltando de WGS para SADOGC
jmenezes@gnu-00:~/Desktop/wutm$ ./utmmudatum -m -63 -s SADOGC -d WGS84
< sat90009u96he.txt | ./utmmudatum -m -63 -d SADOGC -s WGS84
# SAD-69/OGC Chua (UGGI-1967) -63 00'00.0000" altitude elipsoidal
# WGS-84 -> SAD-69/OGC Chua (UGGI-1967)
# 6378137.000 298.257223563 -> 6378160.000 298.250000000 dx 57.000
dy -1.000 dz 41.000
#Pto xE yN altitude elipsoidal
90009 615273.779 8797212.969 176.879
6) diferenca entre 1 e IBGE:
echo 90009U2KIBGE 615215.609 8797173.383 172.62 U2KHE20S -63
IBGE 90009 615215.615 8797173.388 172.603 | awk '{ printf "%s %.3f %.
3f %.3f\n",$1,$9-$2, $10 - $3 , $11 - $4 }'
90009U2KIBGE 0.006 0.005 -0.017
7) diferenca entre 4 e IBGE
jmenezes@gnu-00:~/Desktop/wutm$ echo 90009U2KIBGE 615215.609
8797173.383 172.62 U2KHE20S -63 IBGE 90009 615223.387 8797172.030
180.404 | awk '{ printf "%s %.3f %.3f %.3f\n",$1,$9-$2, $10 - $3 , $11
- $4 }'
90009U2KIBGE 7.778 -1.353 7.784
A diferenca de mais de 7 metros eh muito significativa para minhas
aplicacoes.
>