[Gfoss] conversione da proj.4 a WKT (per SRID in PostGIS)

[Bud P. Bruegger]

Ecco la documentazione di un altro passo che ho fatto.

In postgis, gli SRID nella tabella spatial_reference richiedono la
definizione dello sistema di riferimento in due formati:  proj.4 e WKT.

Avevo fatto sperimenti usando il formato proj.4 (vedi mails precedenti)
e cercavo come trasformare da proj.4 a WKT.  In GRASS, g.proj farebbe
il lavoro ma (per un problema di istallazione/configurazione/ambiente?)
non mi funzionavo.  Cosi nel seguito la soluzione in python usando la
GDAL library.  (Su Debian: apt-get install python-gdal):

postgres at gis:~/py$ cat proj2wkt.py
#!/usr/bin/python
import osr
proj4Def='+proj=cass +lat_0=43.3174975 +lon_0=11.33244638 +x_0=-2.0
+y_0=-16.0 + ellps=intl +units=m +no_defs'
sr = osr.SpatialReference()
sr.ImportFromProj4(proj4Def)
print sr.ExportToPrettyWkt()

postgres at gis:~/py$ ./proj2wkt.py
PROJCS["unnamed",
    GEOGCS["International 1909 (Hayford)",
        DATUM["unknown",
            SPHEROID["intl",6378388,297]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Cassini_Soldner"],
    PARAMETER["latitude_of_origin",43.3174975],
    PARAMETER["central_meridian",11.33244638],
    PARAMETER["false_easting",-2],
    PARAMETER["false_northing",-16],
    UNIT["Meter",1]]

un altro passo fatto ;-)

-b

-- 
Bud P. Bruegger, Ph.D.          +39-0564-488577 (voice),  -21139 (fax) 
   European Chair, Global Collaboration Forum on eID
   Chair, Porvoo Subgroup on collab. govs/operating systems
   Leader of the Permanent eID Status Observatory (PESO) project
Servizio Elaborazione Dati       e-mail:  bud at comune.grosseto.it
Comune di Grosseto               jabber:  bud at jabber.no
Via Ginori, 43                   http://www.comune.grosseto.it/
58100 Grosseto (Tuscany, Italy)
http://www.comune.grosseto.it/interopEID/