[OSGeo-Brasil] Sistemas de referências incorretos no Proj, GDAL/OGR e PostGIS

Daniel Araujo Miranda miranda.dam em dpf.gov.br
Quinta Abril 26 17:36:30 EDT 2012 

Caros,

Para adicionar à lista, seguem novas transformações utilizando o método 
NTv2 (grids), referentes à dica do Etienne, que indicou onde achar a 
transformação no site epsg-registry.org:

 > O EPSG tem entradas de tipo "Coordinate Transformation", voce pode ver
 > no site http://www.epsg-registry.org/ e colocar "Coordinate
 > Transformation"  no campo "Type:" e "sad69" no campo "Name"

SAD69 a partir de 95 (EPSG 4618)
   +proj=longlat +ellps=aust_SA +nadgrids=sad96_003.gsb +nodefs

SAD69 original (EPSG 4291)
   +proj=longlat +ellps=aust_SA +nadgrids=sad69_003.gsb +nodefs

Córrego alegre 1961 (EPSG 5524)
   +proj=longlat +ellps=intl +nadgrids=ca61_003.gsb +nodefs

Córrego alegre 1970-72 (EPSG 4225)
   +proj=longlat +ellps=intl +nadgrids=ca7072_003.gsb +nodefs

Tem um texto em 
"http://www.geobases.es.gov.br/portal/attachments/032_SAD69_SIRGAS2000_ArcGis_NTV2.pdf" 
que compara os dois métodos.

--Miranda

P.S.1: Testei a primeira transformação, e nos pontos testados a 
diferença ficou entre 4cm (Distrito Federal) e 1m (Pará) entre a 
transformação utilizando grids e a utilizando parâmetros towgs84.

P.S.2: Os arquivos *.gsb são distribuídos junto com o sofrware progrid 
do IBGE e devem ser copiados para a pasta "/usr/share/proj" para esses 
parâmetros funcionarem com o proj. AVISO: o proj não dá erro se você 
esquecer de copiar o arquivo ou errar o nome. Cheque os resultados.


On 04/26/2012 12:31 PM, Daniel Araujo Miranda wrote:
> Segue a versão mais recente, após contribuição do Israel e errata do
> Astro chuá.
>
> *Sistemas de coordenadas
>
> Organização:
> Datum
> Projeção
> Código EPSG
> String PROJ
> Referência
>
> sad69 a partir de 95
> longlat
> 4618
> +proj=longlat +ellps=aust_SA +towgs84=-67.35,+3.88,-38.22
> RPR 01/2005 de 25/02/2005 - IBGE
> (ftp://geoftp.ibge.gov.br/documentos/geodesia/projeto_mudanca_referencial_geodesico/legislacao/rpr_01_25fev2005.pdf)
>
>
> sad69 antes de 95
> longlat
> 4291
> +proj=longlat +ellps=aust_SA +towgs84=-66.87,+4.37,-38.52
> Resolução da Presidência do IBGE n° 23, de 21/02/89 (R.PR 23/89)
> (ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/rpr_2389.pdf)
>
> sicad (obs: só aceita zona UTM 23s)
> longlat
>
> +proj=longlat +ellps=intl +towgs84=-144.35,+242.88,-33.22
> GDF DECRETO 32.575-2010_SICAD-SIRGAS,
> www.terracap.df.gov.br/internet/arquivos/0073607515.doc
>
> Córrego Alegre
> longlat
> 4225
> +proj=longlat +ellps=intl +towgs84=-205.57,+168.77,-4.12
> https://www.mar.mil.br/dhn/chm/download/ita06.pdf, composição de
> ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/bservico1602.pdf com
> ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/rpr_2389.pdf
>
> Astro Chuá
> longlat
> 4224
> +proj=longlat +ellps=intl +towgs84=-143.87,+243.37,-33.52
> http://www.topografia.ufsc.br/Geodesia.pdf (referência não oficial)
>
>
> *Elipsóides utilizados
> aust_SA
> a=6378160.0
> rf=298.25
> Australian Natl & S. Amer. 1969
>
> intl
> a=6378388.0
> rf=297.0
> International 1909 (Hayford)
>
>
> *Observações
> Para utilizar a projeção UTM, substitua
> "+proj=longlat"
> pelos parâmetros de projeção, zona e unidades para a zona de interesse.
> Por exemplo:
> "+proj=utm +zone=23 +south +units=m" para zona 23s
> ou
> "+proj=utm +zone=21 +north +units=m" para zona 21n
>
>
> --Miranda
>
> On 04/25/2012 04:45 PM, Etienne Tourigny wrote:
>> Prezados,
>>
>> os dados de http://spatialreference.org usam o GDAL/OGR para obter os
>> parametros. Eo não tenho certeza sobre a fonte que usam proj.4 e
>> postgis.
>>
>> A fonte de GDAL/OGR é a base EPSG (epgs.org / epsg-registry.org).
>> Qualquer erro tem ser ser encaminhado ao EPSG, e esse serà propagado
>> ao softwares citados.
>>
>> Etienne
>>
>> On Wed, Apr 25, 2012 at 4:24 PM, Luiz Motta<motta.luiz em gmail.com> wrote:
>>> Daniel,
>>>
>>> Vou ver com o Thiago (Terra Legal - MDA/INCRA) como ele está
>>> trabalhando com
>>> as projeções com Postgis e GDAL/OGR.
>>>
>>> Os FOSS4G trabalha seguindo a definição da proj4
>>> (http://spatialreference.org).
>>>
>>> Abs
>>> Luiz
>>>
>>> Em 25 de abril de 2012 16:02, Daniel Araujo
>>> Miranda<miranda.dam em dpf.gov.br>
>>> escreveu:
>>>
>>>> Caros, fiz uma pesquisa recente e descobri que os parâmetros de
>>>> conversão
>>>> para diversos sistemas de coordenadas (SRS, CRS) brasileiros estão
>>>> incorretos nos seguintes aplicativos:
>>>>
>>>> Proj4
>>>> GDAL/OGR
>>>> PostGIS
>>>> e talvez outros.
>>>>
>>>> Eles estão com as mesmas informações que o site
>>>> http://spatialreference.org/ (vide a opção proj4 para revelar os
>>>> parâmetros
>>>> de transformação)
>>>>
>>>> Levantei os seguintes parâmetros que acredito estarem corretos. Estou
>>>> divulgando esses dados aqui para vocês utilizarem se tiverem
>>>> interesse e
>>>> para criticarem se houver algo errado (Motta, me ajudaria muito se pelo
>>>> menos você desse uma olhada!).
>>>>
>>>> --Miranda
>>>>
>>>>
>>>> organização:
>>>> Datum
>>>> Projeção
>>>> EPSG
>>>> String PROJ
>>>> Referência
>>>>
>>>>
>>>>
>>>> sad69 apos 94
>>>> longlat
>>>> 4618
>>>> +proj=longlat +a=6378160.0 +rf=298.25 +towgs84=-67.35,+3.88,-38.22
>>>> RPR 01/2005 de 25/02/2005 - IBGE
>>>> (ftp://geoftp.ibge.gov.br/documentos/geodesia/projeto_mudanca_referencial_geodesico/legislacao/rpr_01_25fev2005.pdf)
>>>>
>>>>
>>>> sad69 antes de 94
>>>> longlat
>>>> 4291
>>>> +proj=longlat +a=6378160.0 +rf=298.25 +towgs84=-66.87,+4.37,-38.52
>>>> Resolução da Presidência do IBGE n° 23, de 21/02/89 (R.PR 23/89)
>>>> (ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/rpr_2389.pdf)
>>>>
>>>> sicad
>>>> longlat
>>>>
>>>> +proj=longlat +a=6378388.0 +rf=297.00
>>>> +towgs84=-144.35,+242.88,-33.22
>>>> GDF DECRETO 32.575-2010_SICAD-SIRGAS,
>>>> www.terracap.df.gov.br/internet/arquivos/0073607515.doc
>>>>
>>>> sicad
>>>> utm 23s
>>>>
>>>> +proj=utm +zone=23 +south +units=m +a=6378388.0 +rf=297.00
>>>> +towgs84=-144.35,+242.88,-33.22
>>>> GDF DECRETO 32.575-2010_SICAD-SIRGAS,
>>>> www.terracap.df.gov.br/internet/arquivos/0073607515.doc
>>>>
>>>> Córrego Alegre
>>>> longlat
>>>> 4225
>>>> +proj=longlat +a=6378388.0 +rf=297.00
>>>> +towgs84=-205.57,+168.77,-4.12
>>>> https://www.mar.mil.br/dhn/chm/download/ita06.pdf, composição de
>>>> ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/bservico1602.pdf com
>>>> ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/rpr_2389.pdf
>>>>
>>>> Astro Chuá
>>>> longlat
>>>> 4224
>>>> +proj=longlat +a=6378388.0 +rf=297.00
>>>> +towgs84=-143.87,+243.37,+33.52
>>>> http://www.topografia.ufsc.br/Geodesia.pdf (referência não oficial)
>>>> _______________________________________________
>>>> Brasil mailing list
>>>> Brasil em lists.osgeo.org
>>>> http://lists.osgeo.org/mailman/listinfo/brasil
>>>
>>>
>>>
>>> _______________________________________________
>>> Brasil mailing list
>>> Brasil em lists.osgeo.org
>>> http://lists.osgeo.org/mailman/listinfo/brasil
>>>
>>
> _______________________________________________
> Brasil mailing list
> Brasil em lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/brasil
>

Mais detalhes sobre a lista de discussão Brasil