[Portugal] Transformação de coordenadas no PostGis com Grelhas

[Luís de Sousa]

Olá a todos,

Descarreguei o ficheiro das grelhas da página do prof JA Gonçalves e
instalei nos directórios do proj; aprentemente está tudo a funcionar.
No entanto quando invoco a transformação no PostGis o resultado tem um
erro de cerca de 100 metros, aparentemente não está a usar o ficheiros
das grelhas. Em anexo segue um log destas operações.

Necessito de fazer alguma operação adicional para que o PostGis use as grelhas?

Obrigado,

Luís

P.S.: O meu sistema é Ubuntu 9.10
-------------- next part --------------
$ cat in.txt
-74584.120 2459.620

$ cs2cs +init=pt:dlxhg +to +init=pt:pttm06 < in.txt
-74582.57	2460.86 0.00

$ psql -d pontosDTU -h localhost -p 5432 -U arhtappPassword for user arhtapp: 
psql (8.4.4)
SSL connection (cipher: DHE-RSA-AES256-SHA, bits: 256)
Type "help" for help.

pontosDTU=# SELECT PostGIS_Full_Version();
                                  postgis_full_version                                  
----------------------------------------------------------------------------------------
 POSTGIS="1.4.0" GEOS="3.2.2-CAPI-1.6.2" PROJ="Rel. 4.7.1, 23 September 2009" USE_STATS
(1 row)

pontosDTU=# UPDATE spatial_ref_sys
pontosDTU-#    SET proj4text = '+proj=tmerc +lat_0=39.66666666666666 +lon_0=1 +k=1 +x_0=0 +y_0=0 +ellps=intl +pm=lisbon +units=m +no_defs +nadgrids=ptLX_e89.gsb'
pontosDTU-#  WHERE srid = 20791;
UPDATE 1

pontosDTU=# SELECT X(ST_Transform(PointFromText('POINT(-74584.120 2459.620)', 20791), 3763)),
pontosDTU-#        Y(ST_Transform(PointFromText('POINT(-74584.120 2459.620)', 20791), 3763));
         x         |        y         
-------------------+------------------
 -74477.6196590331 | 2281.82581261726
(1 row)