Claudio,
Bom dia.
On 4 jul, 09:17, Cláudiogeop <
claudiog...@gmail.com> wrote:
> Opa Julio, obrigado pela explicação, seria reprojeção simples onde tenho um
> arquivo shapefile no sistema geodésico WGS 84 e queria como saída Projeção
> UTM, Sistema de Coordenda SAD 69, Fuso 24S, Meridiano Central -39.
> Estarei testando essa ferramenta FrankW cs2cs.
Em *NIX, para testar o cs2cs pode fazer algo assim com entradas em
graus decimais e sinal algebrico:
cs2cs +proj=latlong +ellps=GRS67 +nadgrids=sad96_003.gsb +to
+proj=latlong +ellps=WGS84 +datum=WGS84 -E -v -I -w4 <<EOF
-53.80996650W -21.72888475S 720.470
-53.80996765W -21.72888557S 720.460
EOF
Repare que a altitude permanece quando deveria mudar tambem. O velho
problema do 3D que venho questionando faz tempo.
Eu implementei a transformacao de Molodensky em C para *NIX, o meu
fonte esta aqui, xdatum.c:
http://www.gpl2.com.br/gnu_scripts/c_qgis_spatia_pgis/
Vi que saiu o EPSG 7.11 em Abril de 2012, meus testes foram em janeiro
de 2012, e na versao 7.8 as definicoes, como falei, estavam erradas, /
nao/ sei se houveram correcoes para o Brasil, tenho ate de ver. Duro
eh arrumar um tempinho para isso :-(.
No PostGIS seria algo assim:
Se seu shp esta em WGS84=4326 e vai gerar SAD69=29194, pode fazer isso
no PostGIS.
Joga o shp para o PostGIS com shp2pgis.
Dentro do PostGIS seria algo assim:
SELECT ST_Transform(sua_geom_wgs, 29194) FROM sua_tabela;
ou na CLI *NIX:
$ psql seu_bando -c "SELECT ST_asEWKT(ST_Transform(sua_geom_wgs,
29194)) FROM sua_tabela;" > saida_utm_sad69.txt
Depois informe os resultados.
[]s
julio menezes