Parametros de referenciamento

[André Brito]

Galera!

Alguém ai já teve aquela trabalheira danada pra transformar algum .PRJ
do ArcGIS pra sintaxe PROJ.4?
Uma pessoa muito inspirada resolveu fazer uma coordenada personalizada
e misturou SICAD(O sistema de Coordenadas do DF, q deveria ter sido
abandonado à anos) com UTM. Consegui encontrar o SICAD na sintaxe do
PROJ.4, mas logo q li o .PRJ, vi q tinha alguma coisa a mais.
Notei q toda vez q eu importava o arquivo, ele vinha com coordenadas
UTM, e não latlong.
Já tentei utilizar uma espécie de importador de sistema de
coordenadas, em Python, mas ele sempre me retorna um sistema UTM, na
zona23S. Mas sem especificar o elipsóide e algumas outras variavéis,
as geometrias não se encaixam direito com outros sistemas de
coordenadas.

O .PRJ original, encontra-se assim:

PROJCS["SICAD",GEOGCS["GCS_CHUA",DATUM["D_Chua",SPHEROID["International_1924",6378388.0,297.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",10000000.0],PARAMETER["Central_Meridian",-45.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]

O importador do python sempre me passa a seguinte configuração:

+proj=utm +zone=23 +south +ellps=intl +units=m +no_defs

Estou tentando ler a documentação do PROJ.4, mas tá um pouco tenso,
ainda não encontrei nenhum exemplo q possa me elucidar.

Alguém pode ajudar nessa tradução ai?!

Muito obrigado

--
André Filipe de Assunção e Brito

[Ander]

Explicando muito rapidamente, 

O seu prj original trata-se de um sistema projetada, logo, sempre retornará em coordenadas projetadas, e não em coordenadas geográficas como deseja.

Para ver em coordenadas geográficas, vc deve converter de um sistema para outro (Chuá para Wgs84 por exemplo), e depois para uma Wgs 84 projetada para wgs84 geográfica.

--
 
Anderson Roberto da Silva
Engenheiro Cartógrafo

[yjmenezes]

Andre,

Bom dia.

Em segunda-feira, 13 de agosto de 2012 18h44min02s UTC-3, decko escreveu:

Galera!

Alguém ai já teve aquela trabalheira danada pra transformar algum .PRJ
do ArcGIS pra sintaxe PROJ.4?
Uma pessoa muito inspirada resolveu fazer uma coordenada personalizada
e misturou SICAD(O sistema de Coordenadas do DF, q deveria ter sido
abandonado à anos) com UTM. Consegui encontrar o SICAD na sintaxe do
PROJ.4, mas logo q li o .PRJ, vi q tinha alguma coisa a mais.
Notei q toda vez q eu importava o arquivo, ele vinha com coordenadas
UTM, e não latlong.

O SICAD foi estabelecido localmente para a construcao de Brasilia.
Eu escrevi um codigo C para transformacao em SAD69.
Os parametros estao neste fragmento de codigo:

    case 8 :  MD_.Datum_1 = ASTROCHUA;
              MD_.Datum_2 = SAD69;
              MD_.a1 = 6378388.0;
              MD_.f1 = (1.0 / 297.0);
              MD_.a2 = 6378160.0;
              MD_.f2 = (1.0 / 298.25);
              MD_.dx =  -77.0;
              MD_.dy =  239.0;
              MD_.dz =    5.0;
            break; /* 8 */

Meu codigo esta aqui:
http://www.gpl2.com.br/gnu_scripts/c_qgis_spatia_pgis/xdatum_2012-01-06.zip

Teria de alterar e recompilar para o AstroChua.
Caso seus dados estejam num formato texto ASCII,  se me enviar um arquivo eu transformo para voce, basta me indicar os sistemas origem e destino.

[]s
julio menezes

[yjmenezes]

Andre,

So acrescentando,

1- Se nao me falha a memoria, a CODEPLAN usa fuso 23 expandido na projecao UTM, ja que Brasilia eh cortada pelo meridiano 48.
2- Os parametros que passei referem-se ao SAD69 original.

[]s
julio menezes

[André Brito]

Ander

Eu posso utilizar projeção "on the fly" do qgis, não? Ainda que as
coordenadas originais sejam projetadas, no processo ele faz a
transformação, não?

2012/8/13 Ander <bya...@gmail.com>:

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

--
André Filipe de Assunção e Brito

Coordenadoria de Articulação Intergovernamental
Secretaria de Estado de Governo
Governo do Distrito Federal
61 78129408
ID: 55*133*2058

[André Brito]

Galera

Acho q cheguei a algo
Lendo a documentação do Proj.4 e me baseando no .prj q li, desenvolvi
a seguinte linha.
Tb resolvi fazer alguns experimentos e já utilizar UTM como projeção.

Pelos testes bem simplórios q fiz aqui, chegou mais perto do q eu
estava esperando.
Habilitei a reprojeção "on-the-fly", para Sirgas2000, projeção latlong.

Com um arquivo de pontos, a aproximação está na casa dos segundos.
Se alguém precisar, ou quiser ajudar a melhorar, as linhas seguem abaixo:

SICAD DF
+ellps=intl +lat_0=0.0 +lon_0=-45.0 +k_0=0.9996 +x_0=500000.0
+y_0=10000000.0 +proj=tmerc +towgs84=-144.35,+242.88,-33.22

SICAD DF (Baseado em AstroChua)
+ellps=intl +lat_0=0.0 +lon_0=-45.0 +k_0=0.9996 +x_0=500000.0
+y_0=10000000.0 +proj=tmerc +towgs84=-143.87,+243.37,-33.52

SICAD DF (UTM)
+ellps=intl +k_0=0.9996 +proj=utm +zone=23 +south +unit=m
+towgs84=-143.87,+243.37,-33.52

2012/8/14 yjmenezes <yjme...@gmail.com>:

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

--
André Filipe de Assunção e Brito

[yjmenezes]

Ander,

Bom dia.

Estes parametros aqui batem com os que postei anteriormente de AstroChua para SAD69 original.
Estes levam de AstroChua para WGS84

SICAD DF (Baseado em AstroChua)
+ellps=intl +lat_0=0.0 +lon_0=-45.0 +k_0=0.9996 +x_0=500000.0
+y_0=10000000.0 +proj=tmerc +towgs84=-143.87,+243.37,-33.52

[]s
julio