Difference between revisions of "Brazilian Coordinate Reference Systems"

From OSGeo
Jump to navigation Jump to search
 
(9 intermediate revisions by 2 users not shown)
Line 1: Line 1:
This page discusses the Coordinate Reference Systems used in Brazil. It is written in Brazilian Portuguese and started as an attempt to reach a consensus on the transformation parameters that should be submitted to the "upstream" of open source programs such as proj4. The discussion was held mostly on the brasil@lists.osgeo.org starting XXXX (put date here) (XXXX put link to list archives here XXXX).
+
This page discusses the Coordinate Reference Systems used in Brazil. It started as an attempt to reach a consensus on the transformation parameters that should be submitted to the "upstream" of open source programs such as Proj4, libgeotiff, GDAL and Quantum GIS. Part of this document is written in Brazilian Portuguese, in-text translations are provided.
  
 +
The discussion was held mostly on the brasil@lists.osgeo.org mailing list starting April 25, 2012 (http://lists.osgeo.org/pipermail/brasil/2012-April/000657.html).
  
= Sistemas de referência utilizados no Brasil =
+
The following references are a substancial source of information (in Portuguese), as well as others below:
  
O Brasil já utilizou diversos sistemas de referência. A lei XXXX determina que, a partir de XXXX, será adotado unicamente o sistema SIRGAS 2000.
+
http://www.ibge.gov.br/home/geociencias/geodesia/pmrg/faq.shtm
 +
 
 +
http://www.ibge.gov.br/home/geociencias/geodesia/param_transf/default_param_transf.shtm
 +
 
 +
= List of reference systems =
 +
 
 +
Brazil has used various reference systems. Brazilian law states that, starting from 2015, the SIRGAS 2000 system will be the only official coordinate reference system.
 +
 
 +
The following list provides strings for use with the PROJ4 library and compatible implementations using the long/lat projection. For adapting the string to UTM or other projections, see the appropriate section of this page.
 +
 
 +
 
 +
 
 +
''O Brasil já utilizou diversos sistemas de referência. A [ftp://geoftp.ibge.gov.br/documentos/geodesia/projeto_mudanca_referencial_geodesico/legislacao/rpr_01_25fev2005.pdf resolução 01/2005 do IBGE] determina que, a partir de 25/fevereiro/2015 será adotado unicamente o sistema SIRGAS 2000, havendo um período de transição não superior a 10 anos a partir da publicação da resolução (2005) em que também poderá ser utilizado o sistema SAD69.''
 +
 
 +
''Segue uma lista dos sistemas de coordenadas brasileiros (também denominados sistemas de referência espacial) contendo os parâmetros para utilização na biblioteca PROJ4 e implementações compatíveis. Os parâmetros são fornecidos utilizando a projeção long/lat. Para utilizar projeção UTM ou outras, é necessário alterar os parâmetros de acordo com o descrito mais abaixo nesta página.''
  
 
== SIRGAS 2000 ==
 
== SIRGAS 2000 ==
XXXX
+
Projection:        longlat
 +
EPSG code:          4989 (3D) or 4674 (2D)
 +
PROJ string:        '+proj=longlat +ellps=GRS80 +towgs84=0,0,0 +no_defs'
 +
Reference:          http://www.epsg-registry.org/ (unofficial reference)
 +
Note:              Considered identical to WGS84
 +
 
 +
== SAD69 from 1995 onwards ''(SAD69 a partir de 1995)''==
 +
Projection:        longlat
 +
EPSG code:          5527
 +
PROJ string:        '+proj=longlat +ellps=aust_SA +towgs84=-67.35,3.88,-38.22'
 +
PROJ string (ntv2): '+proj=longlat +ellps=aust_SA +nadgrids=sad96_003.gsb +nodefs'
 +
Reference:          [ftp://geoftp.ibge.gov.br/documentos/geodesia/projeto_mudanca_referencial_geodesico/legislacao/rpr_01_25fev2005.pdf ''Resolução 01/2005 de 25 de fevereiro de 2005 - IBGE'']
 +
Note:              Previously EPSG code 4618 was listed here, but that is actually the older version,
 +
                    acording to EPSG database (www.epsg.org).
 +
 
 +
== SAD69 before 1995 ''(SAD69 antes de 1995)''==
 +
Projection:        longlat
 +
EPSG code:          4618
 +
PROJ string:        '+proj=longlat +ellps=aust_SA +towgs84=-66.87,4.37,-38.52'
 +
PROJ string (ntv2): '+proj=longlat +ellps=aust_SA +nadgrids=sad69_003.gsb +nodefs'
 +
Reference:          [ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/rpr_2389.pdf ''Resolução 23/89 de 21 de fevereiro de 1989 - IBGE'']
 +
Note:              EPSG code 4618 was moved from the new SAD69 to this one, which is actually the correct
 +
                    designation according to EPSG database (www.epsg.org).
 +
Note 2:            Do not use EPSG 4291 (deprecated). EPSG codes 4291 and 4618 do not display any
 +
                    difference in EPSG database except for the comment that 4618 is valid in Brazil
 +
                    only and the datum for 4291 is missing the origin description.
 +
 
 +
==SICAD ''(Sistema Cartográfico do Distrito Federal)''==
 +
Projection:        longlat
 +
EPSG code:          ???
 +
PROJ string:        '+proj=longlat +ellps=intl +towgs84=-144.35,242.88,-33.22'
 +
Reference:          [http://www.terracap.df.gov.br/internet/arquivos/0073607515.doc ''Decreto 32.575/2010 (SICAD-SIRGAS) - GDF'']
 +
Note:              only valid in UTM 23s zone
 +
 
 +
==Corrego Alegre 1961 ''(Córrego Alegre 1961)''==
 +
Projection:        longlat
 +
EPSG code:          5524
 +
PROJ string:        ???
 +
PROJ string (ntv2): '+proj=longlat +ellps=intl +nadgrids=ca61_003.gsb +nodefs'
 +
Reference:          ???
 +
 
 +
==Corrego Alegre 1970-72 ''(Córrego Alegre 1970-72)''==
 +
Projection:        longlat
 +
EPSG code:          4225
 +
PROJ string:        '+proj=longlat +ellps=intl +towgs84=-205.57,168.77,-4.12'
 +
PROJ string (ntv2): '+proj=longlat +ellps=intl +nadgrids=ca7072_003.gsb +nodefs'
 +
Reference:          [http://www.mar.mil.br/dhn/chm/download/ita06.pdf ''Instrução Técnica A-06A - Marinha do Brasil''] (calculated by combining the transformations in the two following documents)
 +
                    [ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/bservico1602.pdf ''Resolução 22/83 de 21 de juho de 1983 - IBGE'']
 +
                    [ftp://geoftp.ibge.gov.br/documentos/geodesia/pdf/rpr_2389.pdf ''Resolução 23/89 de 21 de fevereiro de 1989 - IBGE'']
 +
 
 +
==Chua ''(Astro Chuá)''==
 +
Projection:        longlat
 +
EPSG code:          4224
 +
PROJ string:        '+proj=longlat +ellps=intl +towgs84=-143.87,243.37,-33.52'
 +
Reference:          http://www.topografia.ufsc.br/Geodesia.pdf (unofficial reference)
 +
Note:              Possibly the same as SICAD
 +
 
 +
==Summary Table==
 +
{| border="1" cellpadding="2" cellspacing="2"
 +
! Name
 +
! EPSG
 +
! towgs84
 +
! towgs84 currently in use by gdal (incorrect)
 +
! ntv2 grid from [ftp://geoftp.ibge.gov.br/aplicativos/transformacao_coordenadas/ IBGE]
 +
|-
 +
| SAD69 || 4618 || -67.35,3.88,-38.22 || -57,1,-41 || sad96_003.gsb
 +
|-
 +
| SAD69 before 1995 || 4291 || -66.87,4.37,-38.52 || -57,1,-41 || sad69_003.gsb
 +
|-
 +
| SICAD || ??? || -144.35,242.88,-33.22 || - || -
 +
|-
 +
| Corrego Alegre 1961 || 5524 || ??? || - || ca61_003.gsb
 +
|-
 +
| Corrego Alegre 1970-72 || 4225 || -205.57,168.77,-4.12 || -206,172,-6 || ca7072_003.gsb
 +
|-
 +
| Chua || 4224 || -143.87,243.37,-33.52 || -134,229,-29 || -
 +
|}
 +
 
 +
=Ellipsoids in use=
 +
 
 +
==aust_SA==
 +
PROJ name:        Australian Natl &  S. Amer. 1969
 +
PROJ string:      '+ellps=aust_SA'
 +
Semimajor radius  a=6378160.0
 +
Flattening        rf=298.25
 +
EPSG name:        GRS 1967 Modified
 +
EPSG code:        7050
 +
 
 +
==intl==
 +
PROJ name:        International 1909 (Hayford)
 +
PROJ string:      '+ellps=intl'
 +
Semimajor radius  a=6378388.0
 +
Flattening        rf=297.0
 +
EPSG name:        International 1924
 +
EPSG code:        7022
 +
 
 +
==GRS80==
 +
PROJ name:        GRS 1980(IUGG, 1980)
 +
PROJ string:      '+ellps=GRS80'
 +
Semimajor radius  a=6378137.0
 +
Flattening        rf=298.257222101
 +
EPSG name:        GRS 1980
 +
EPSG code:        7019
 +
 
 +
=Using UTM projected coordinates=
 +
To use the UTM projection, substitute "+proj=longlat" with the desired projection parameter, zone and units for the area of interest
  
==SAD69 a partir de 1995==
+
For exemple, take the SAD69 long/lat string:
Projeção:          longlat
+
  +proj=longlat +ellps=aust_SA +towgs84=-67.35,3.88,-38.22
Código EPSG:       4618
+
To use the UTM projection in zone 23s, the parameters become:
  String PROJ:        +proj=longlat +ellps=aust_SA +towgs84=-67.35,+3.88,-38.22
+
  +proj=utm +zone=23 +south +units=m +ellps=aust_SA +towgs84=-67.35,3.88,-38.22
  Versão grid (ntv2): +proj=longlat +ellps=aust_SA +nadgrids=sad96_003.gsb +nodefs
+
To use the UTM projection in zone 21n, the parameters become:
Referência:
+
+proj=utm +zone=21 +north +units=m +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==
+
NTV2 grids (*.gsb files) are distributed with the Progrid software package from IBGE (ftp://geoftp.ibge.gov.br/aplicativos/transformacao_coordenadas/).
Projeção:          longlat
 
Código EPSG:        4291
 
String PROJ:        +proj=longlat +ellps=aust_SA +towgs84=-66.87,+4.37,-38.52
 
Versão grid (ntv2): +proj=longlat +ellps=aust_SA +nadgrids=sad69_003.gsb +nodefs
 
Referência:
 
  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==
 
Sistema Cartográfico do Distrito Federal
 
(obs: só é válido na zona UTM 23s)
 
Projeção:          longlat
 
Código EPSG:        ???
 
String PROJ:        +proj=longlat +ellps=intl +towgs84=-144.35,+242.88,-33.22
 
Referência:
 
  GDF DECRETO 32.575-2010_SICAD-SIRGAS, http://www.terracap.df.gov.br/internet/arquivos/0073607515.doc
 
  
==Córrego Alegre 1961==
 
Projeção:          longlat
 
Código EPSG:      5524
 
Versão grid (ntv2): +proj=longlat +ellps=intl +nadgrids=ca61_003.gsb +nodefs
 
Referência:
 
  XXXX
 
  
==Córrego Alegre 1970-72==
 
Projeção:          longlat
 
Código EPSG:      4225
 
String PROJ:        +proj=longlat +ellps=intl +towgs84=-205.57,+168.77,-4.12
 
Versão grid (ntv2): +proj=longlat +ellps=intl +nadgrids=ca7072_003.gsb +nodefs
 
Referência:
 
  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á==
+
''Para utilizar a projeção UTM ao invés de long/lat, substitua "+proj=longlat" pelos parâmetros de projeção, zona e unidades para a zona de interesse.''
Projeção:          longlat
 
Código EPSG:      4224
 
String PROJ:        +proj=longlat +ellps=intl +towgs84=-143.87,+243.37,-33.52
 
Referência:
 
  http://www.topografia.ufsc.br/Geodesia.pdf (referência não oficial)
 
  
=Elipsóides utilizados=
+
''Por exemplo, os parâmetros long/lat do SAD69 são:''
aust_SA
+
+proj=longlat +ellps=aust_SA +towgs84=-67.35,3.88,-38.22
nome: Australian Natl&  S. Amer. 1969
+
''Para uso na zona UTM 23s os parâmetros modificados são:''
  a=6378160.0
+
  +proj=utm +zone=23 +south +units=m +ellps=aust_SA +towgs84=-67.35,3.88,-38.22
rf=298.25
+
''Para uso na zona UTM 21n os parâmetros modificados são:''
intl
+
  +proj=utm +zone=21 +north +units=m +ellps=aust_SA +towgs84=-67.35,3.88,-38.22
nome: International 1909 (Hayford)
 
  a=6378388.0
 
rf=297.0
 
  
=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
 
  
Comparação entre o método grid (ntv2) e parâmetros towgs84: [http://www.geobases.es.gov.br/portal/attachments/032_SAD69_SIRGAS2000_ArcGis_NTV2.pdf SAD69 para SIRGAS2000]
+
''Par uma comparação entre o método grid (ntv2) e parâmetros towgs84 na realização de conversões de coordenadas, veja:'' [http://www.geobases.es.gov.br/portal/attachments/032_SAD69_SIRGAS2000_ArcGis_NTV2.pdf ''SAD69 para SIRGAS2000'']. ''Normalmente utiliza-se apenas os parâmetros towgs84, que permitem inclusive a conversão de altitudes.''
  
Os arquivos *.gsb são distribuídos junto com o sofrware progrid do
+
''Os arquivos *.gsb necessários para a utilização do método grid (ntv2) são distribuídos junto com o software PROGRID do IBGE ([ftp://geoftp.ibge.gov.br/aplicativos/transformacao_coordenadas/ PROGRID]) e devem ser copiados para a pasta "/usr/share/proj" (linux) para esses parâmetros funcionarem com o proj.''
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.
 
  
 +
''AVISO: o proj não dá erro se você esquecer de copiar o arquivo ou errar o nome do arquivo na string de conversão.''
 +
''Cheque os resultados antes de utilizá-los.''
  
=Pendências=
+
''AVISO: utilize os parâmetros TOWGS84 positivos sem o sinal de "+" na biblioteca proj4as (versão do proj para FLEX da openscales).''
*testar as transformações do tipo grid (ntv2)
+
''Por exemplo: +towgs84=-67.35,3.88,-38.22 (correto); +towgs84=-67.35,+3.88,-38.22 (incorreto, sinal de + antes do 3.88)''
*preencher os campos marcados com XXXX
 
*completar a documentação dos sistemas de coordenadas
 
*atingir consenso na lista antes de mandar para os desenvolvedores upstream ou para a EPSG.
 

Latest revision as of 11:33, 13 September 2014

This page discusses the Coordinate Reference Systems used in Brazil. It started as an attempt to reach a consensus on the transformation parameters that should be submitted to the "upstream" of open source programs such as Proj4, libgeotiff, GDAL and Quantum GIS. Part of this document is written in Brazilian Portuguese, in-text translations are provided.

The discussion was held mostly on the brasil@lists.osgeo.org mailing list starting April 25, 2012 (http://lists.osgeo.org/pipermail/brasil/2012-April/000657.html).

The following references are a substancial source of information (in Portuguese), as well as others below:

http://www.ibge.gov.br/home/geociencias/geodesia/pmrg/faq.shtm

http://www.ibge.gov.br/home/geociencias/geodesia/param_transf/default_param_transf.shtm

List of reference systems

Brazil has used various reference systems. Brazilian law states that, starting from 2015, the SIRGAS 2000 system will be the only official coordinate reference system.

The following list provides strings for use with the PROJ4 library and compatible implementations using the long/lat projection. For adapting the string to UTM or other projections, see the appropriate section of this page.


O Brasil já utilizou diversos sistemas de referência. A resolução 01/2005 do IBGE determina que, a partir de 25/fevereiro/2015 será adotado unicamente o sistema SIRGAS 2000, havendo um período de transição não superior a 10 anos a partir da publicação da resolução (2005) em que também poderá ser utilizado o sistema SAD69.

Segue uma lista dos sistemas de coordenadas brasileiros (também denominados sistemas de referência espacial) contendo os parâmetros para utilização na biblioteca PROJ4 e implementações compatíveis. Os parâmetros são fornecidos utilizando a projeção long/lat. Para utilizar projeção UTM ou outras, é necessário alterar os parâmetros de acordo com o descrito mais abaixo nesta página.

SIRGAS 2000

Projection:         longlat
EPSG code:          4989 (3D) or 4674 (2D)
PROJ string:        '+proj=longlat +ellps=GRS80 +towgs84=0,0,0 +no_defs'
Reference:          http://www.epsg-registry.org/ (unofficial reference)
Note:               Considered identical to WGS84

SAD69 from 1995 onwards (SAD69 a partir de 1995)

Projection:         longlat
EPSG code:          5527
PROJ string:        '+proj=longlat +ellps=aust_SA +towgs84=-67.35,3.88,-38.22'
PROJ string (ntv2): '+proj=longlat +ellps=aust_SA +nadgrids=sad96_003.gsb +nodefs'
Reference:          Resolução 01/2005 de 25 de fevereiro de 2005 - IBGE
Note:               Previously EPSG code 4618 was listed here, but that is actually the older version,
                    acording to EPSG database (www.epsg.org).

SAD69 before 1995 (SAD69 antes de 1995)

Projection:         longlat
EPSG code:          4618
PROJ string:        '+proj=longlat +ellps=aust_SA +towgs84=-66.87,4.37,-38.52'
PROJ string (ntv2): '+proj=longlat +ellps=aust_SA +nadgrids=sad69_003.gsb +nodefs'
Reference:          Resolução 23/89 de 21 de fevereiro de 1989 - IBGE
Note:               EPSG code 4618 was moved from the new SAD69 to this one, which is actually the correct
                    designation according to EPSG database (www.epsg.org).
Note 2:             Do not use EPSG 4291 (deprecated). EPSG codes 4291 and 4618 do not display any
                    difference in EPSG database except for the comment that 4618 is valid in Brazil
                    only and the datum for 4291 is missing the origin description.

SICAD (Sistema Cartográfico do Distrito Federal)

Projection:         longlat
EPSG code:          ???
PROJ string:        '+proj=longlat +ellps=intl +towgs84=-144.35,242.88,-33.22'
Reference:          Decreto 32.575/2010 (SICAD-SIRGAS) - GDF
Note:               only valid in UTM 23s zone

Corrego Alegre 1961 (Córrego Alegre 1961)

Projection:         longlat
EPSG code:          5524
PROJ string:        ???
PROJ string (ntv2): '+proj=longlat +ellps=intl +nadgrids=ca61_003.gsb +nodefs'
Reference:          ???

Corrego Alegre 1970-72 (Córrego Alegre 1970-72)

Projection:         longlat
EPSG code:          4225
PROJ string:        '+proj=longlat +ellps=intl +towgs84=-205.57,168.77,-4.12'
PROJ string (ntv2): '+proj=longlat +ellps=intl +nadgrids=ca7072_003.gsb +nodefs'
Reference:          Instrução Técnica A-06A - Marinha do Brasil (calculated by combining the transformations in the two following documents)
                    Resolução 22/83 de 21 de juho de 1983 - IBGE
                    Resolução 23/89 de 21 de fevereiro de 1989 - IBGE

Chua (Astro Chuá)

Projection:         longlat
EPSG code:          4224
PROJ string:        '+proj=longlat +ellps=intl +towgs84=-143.87,243.37,-33.52'
Reference:          http://www.topografia.ufsc.br/Geodesia.pdf (unofficial reference)
Note:               Possibly the same as SICAD

Summary Table

Name EPSG towgs84 towgs84 currently in use by gdal (incorrect) ntv2 grid from IBGE
SAD69 4618 -67.35,3.88,-38.22 -57,1,-41 sad96_003.gsb
SAD69 before 1995 4291 -66.87,4.37,-38.52 -57,1,-41 sad69_003.gsb
SICAD ??? -144.35,242.88,-33.22 - -
Corrego Alegre 1961 5524 ??? - ca61_003.gsb
Corrego Alegre 1970-72 4225 -205.57,168.77,-4.12 -206,172,-6 ca7072_003.gsb
Chua 4224 -143.87,243.37,-33.52 -134,229,-29 -

Ellipsoids in use

aust_SA

PROJ name:         Australian Natl &  S. Amer. 1969
PROJ string:       '+ellps=aust_SA'
Semimajor radius   a=6378160.0
Flattening         rf=298.25
EPSG name:         GRS 1967 Modified 
EPSG code:         7050

intl

PROJ name:         International 1909 (Hayford)
PROJ string:       '+ellps=intl'
Semimajor radius   a=6378388.0
Flattening         rf=297.0
EPSG name:         International 1924
EPSG code:         7022

GRS80

PROJ name:         GRS 1980(IUGG, 1980)
PROJ string:       '+ellps=GRS80'
Semimajor radius   a=6378137.0
Flattening         rf=298.257222101
EPSG name:         GRS 1980
EPSG code:         7019

Using UTM projected coordinates

To use the UTM projection, substitute "+proj=longlat" with the desired projection parameter, zone and units for the area of interest

For exemple, take the SAD69 long/lat string:

+proj=longlat +ellps=aust_SA +towgs84=-67.35,3.88,-38.22

To use the UTM projection in zone 23s, the parameters become:

+proj=utm +zone=23 +south +units=m +ellps=aust_SA +towgs84=-67.35,3.88,-38.22

To use the UTM projection in zone 21n, the parameters become:

+proj=utm +zone=21 +north +units=m +ellps=aust_SA +towgs84=-67.35,3.88,-38.22

NTV2 grids (*.gsb files) are distributed with the Progrid software package from IBGE (ftp://geoftp.ibge.gov.br/aplicativos/transformacao_coordenadas/).



Para utilizar a projeção UTM ao invés de long/lat, substitua "+proj=longlat" pelos parâmetros de projeção, zona e unidades para a zona de interesse.

Por exemplo, os parâmetros long/lat do SAD69 são:

+proj=longlat +ellps=aust_SA +towgs84=-67.35,3.88,-38.22

Para uso na zona UTM 23s os parâmetros modificados são:

+proj=utm +zone=23 +south +units=m +ellps=aust_SA +towgs84=-67.35,3.88,-38.22

Para uso na zona UTM 21n os parâmetros modificados são:

+proj=utm +zone=21 +north +units=m +ellps=aust_SA +towgs84=-67.35,3.88,-38.22


Par uma comparação entre o método grid (ntv2) e parâmetros towgs84 na realização de conversões de coordenadas, veja: SAD69 para SIRGAS2000. Normalmente utiliza-se apenas os parâmetros towgs84, que permitem inclusive a conversão de altitudes.

Os arquivos *.gsb necessários para a utilização do método grid (ntv2) são distribuídos junto com o software PROGRID do IBGE (PROGRID) e devem ser copiados para a pasta "/usr/share/proj" (linux) 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 do arquivo na string de conversão.
Cheque os resultados antes de utilizá-los.
AVISO: utilize os parâmetros TOWGS84 positivos sem o sinal de "+" na biblioteca proj4as (versão do proj para FLEX da openscales).
Por exemplo: +towgs84=-67.35,3.88,-38.22 (correto); +towgs84=-67.35,+3.88,-38.22 (incorreto, sinal de + antes do 3.88)