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@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@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@comune.grosseto.it
Comune di Grosseto jabber: bud@jabber.no
Via Ginori, 43 http://www.comune.grosseto.it/
58100 Grosseto (Tuscany, Italy)
http://www.comune.grosseto.it/interopEID/