Re: [OSGeo-Brasil] Sistemas de referências incorretos no Proj, GDAL/OGR e PostGIS
Etienne Tourigny
etourigny.dev em gmail.com
Quinta Abril 26 08:55:00 EDT 2012
On Wed, Apr 25, 2012 at 6:47 PM, Daniel Araujo Miranda
<miranda.dam em dpf.gov.br> wrote:
> Etienne,
> Salvo engano, o spatiareference.org e o epsg.org não têm os parâmetros de
> transformação (dX, dY, dZ, etc.). Você pode confirmar isso? Até baixei a
> base de dados em postgresql mas não consegui achar essa informação.
> Se for isso mesmo, não tem como mandar as correções para eles, pois o
> problema não é nos elipsóides, é principalmente nos parâmetros de
> transformação, e isso eles não têm.
>
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"
Os dados brutos do EPSG (epsg.org) são a principal fonte de dados para
o GDAL/OGR e GeoTIFF.
Não se porque spatialreference.org não tem os parametros TOWGS84 ...
Na pagina http://spatialreference.org/about/ esta escrito "GDAL and
Proj.4 are the backbone of this website", porém o GDAL tem os
parametros TOWGS84.
O GDAL/OGR tem varios parametros TOWGS84 (no arquivo datum_shift.csv)
para cada DATUM, mas usa o "preferred datum" nas consultas:
$ gdalsrsinfo EPSG:4618
PROJ.4 : '+proj=longlat +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +no_defs '
E no arquivo datum_shift.csv:
$ grep 4618 /home/soft/share/gdal/datum_shift.csv
13,1864,4618,4326,Derived at 84 stations.,"For military purposes only.
Accuracy 15m, 6m and 9m in X, Y and Z
axes.",1358,-45,12.51,-81.4,-29.03,1,0,9603,-57,1,-41,,,,,1
14,1865,4618,4326,"Derived at 10 stations. Note: SAD69 not adopted in
Argentina: see Campo Inchauspe (CRS code 4221).",For military purposes
only. Accuracy 5m in each
axis.,3215,-52.43,-21.78,-73.58,-53.65,1,0,9603,-62,-1,-37,,,,,0
15,1866,4618,4326,"Derived at 4 stations. Note: SAD69 not adopted in
Bolivia: see PSAD56 (CRS code 4248).",For military purposes. Accuracy
15m in each axis.,1049,-22.9,-9.68,-69.66,-57.52,1,0,9603,-61,2,-48,,,,,0
16,1867,4618,4326,Derived at 22 stations.,"For military purposes only.
Accuracy 3m, 5m and 5m in X, Y and Z
axes.",3845,-35.71,7.04,-60.57,-29.03,1,0,9603,-60,-2,-41,,,,,0
17,1868,4618,4326,"Derived at 9 stations. Note: SAD69 not adopted in
Chile.","For military purposes only. Accuracy 15m, 8m and 11m in X, Y
and Z axes.",3227,-45,-17.51,-75.22,-67,1,0,9603,-75,-1,-44,,,,,0
18,1869,4618,4326,"Derived at 7 stations. Note: SAD69 not adopted in
Colombia: see Bogota 1975 (CRS code 4218).","For military purposes
only. Accuracy 6m, 6m and 5m in X, Y and Z
axes.",3229,-4.24,12.51,-79.1,-66.87,1,0,9603,-44,6,-36,,,,,0
19,1870,4618,4326,"Derived at 11 stations. Note: SAD69 not adopted in
Ecuador: see PSAD56 (CRS code 4248).",For military purposes. Accuracy
3m in each axis.,3241,-5,1.45,-81.03,-75.22,1,0,9603,-48,3,-44,,,,,0
20,1871,4618,4326,"Derived at 1 station. Note: SAD69 not adopted in
Ecuador.",For military purposes. Accuracy 25m in each
axis.,2356,-1.41,0.17,-91.71,-89.2,1,0,9603,-47,26,-42,,,,,0
21,1872,4618,4326,"Derived at 5 stations. Note: SAD69 not adopted in
Guyana.","For military purposes only. Accuracy 9m, 5m and 5m in X, Y
and Z axes.",3259,1.19,8.57,-61.39,-56.47,1,0,9603,-53,3,-47,,,,,0
22,1873,4618,4326,"Derived at 4 stations. Note: SAD69 not adopted in
Paraguay.",For military purposes. Accuracy 15m in each
axis.,1188,-27.58,-19.3,-62.64,-54.24,1,0,9603,-61,2,-33,,,,,0
23,1874,4618,4326,"Derived at 6 stations. Note: SAD69 not adopted in
Peru: see PSAD56 (CRS code 4248).",For military purposes. Accuracy 5m
in each axis.,3292,-18.35,-0.04,-81.4,-68.67,1,0,9603,-58,0,-44,,,,,0
24,1875,4618,4326,"Derived at 1 station. Note: SAD69 not adopted in
Trinidad and Tobago.",For military purposes only. Accuracy 25m in each
axis.,3143,9.99,10.89,-61.97,-60.86,1,0,9603,-45,12,-33,,,,,0
25,1876,4618,4326,"Derived at 5 stations. Note: SAD69 not adopted in
Venezuela: see PSAD56 (CRS code 4248).","For military purposes only.
Accuracy 3m, 6m and 3m in X, Y and Z
axes.",3327,0.65,12.25,-73.38,-59.8,1,0,9603,-45,8,-33,,,,,0
26,1877,4618,4326,"Derived by Brazilian Institute of Geography and
Statistics (IBGE) in 1989. Used by ANP. (Note: for historic reasons
associated with one-time web url, tfm version uses initials IGBE, not
IBGE). Replaced by SAD69 to WGS 84 (15) (tfm code 5528).",Medium and
small scale mapping.,3845,-35.71,7.04,-60.57,-29.03,1,0,9603,-66.87,4.37,-38.52,,,,,0
543,1271,4293,4326,"Beware! Schwarzeck CRS uses German legal metres.
Example: Schwarzeck Lat 19d 35m 46.952s S Long 20d 41m 50.649s E;
X=5623409.40 Y=2124618.00 Z=-2125847.62 GLM; X=5623485.86 Y=2124646.89
Z=-2125876.53 m; WGS84 X=5624101.50 Y=2124748.97 Z=2126132.34
m.","?",1169,-31.2,-16.99,8.32,25.28,1,0,9603,615.64,102.08,-255.81,,,,,0
abs
Etienne
> O spatialreference.org tem os parâmetros de alguns. Por exemplo:
> http://spatialreference.org/ref/epsg/4224/
> Os parâmetros citados são:
> +proj=longlat +ellps=intl +towgs84=-134,229,-29,0,0,0,0 +no_defs
> Se isso for o Astro Chuá mesmo, os parâmetros estão errados.
>
> No proj4.7 (pacote proj-data):
> /usr/share/proj/esri:
> <4224> +proj=longlat +ellps=intl +towgs84=-134,229,-29,0,0,0,0 +no_defs
> no_defs <>
> /usr/share/proj/epsg:
> <4224> +proj=longlat +ellps=intl +no_defs <>
>
> Quanto às fontes do proj e postgis: O postgis não especifica (seção 4.3.1 do
> manual), e o proj usa diversas fontes. A lista de arquivos do proj dá uma
> dica:
>
> linhas arquivo
> 7631 /usr/share/proj/epsg
> 5937 /usr/share/proj/esri
> 5189 /usr/share/proj/ntv1_can.dat
> 2985 /usr/share/proj/alaska
>
> Por enquanto, estou expondo os parâmetros para o pessoal poder criticar.
> Quando tivermos confiança nessas transformações, podemos procurar as fontes
> desse dado para atualizá-lo.
>
>
>
> **********************
> Aproveito para enviar a primeira errata:
> 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)
>
> (o sinal do deltaZ estava trocado)
> **********************
>
>
>
> --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
>>>
>>
>
Mais detalhes sobre a lista de discussão Brasil