Creazione di una trasformazione personalizzata basata su una trasformazione a 7 parametri (Bursa Wolf)

Se si vogliono utilizzare i comandi della libreria di FWTools per convertire dal datum catastale al datum WGS84,  il parametro +towgs84=...  che consente tale passaggio (nell'ambito del territorio provinciale di Agrigento, escluse le isole Pelagie) č:

+towgs84=763.901,183.452,519.534,-2.0388,-6.0089,-7.1995,-35.296542

Dove i parametri di rotazione della trasformazione Bursa-Wolf ricavati con ILWIS sono stati cambiati di segno.

Pertanto, per effettuare conversioni di coordinate utilizzando la trasformazione conforme 3D, si puņ utilizzare la seguente stringa
cs2cs  +proj=cass +lat_0=37.41457777777778 +lon_0=13.780637500000 +x_0=0.0 +y_0=0.0 +ellps=bessel  +towgs84=763.901,183.452,519.534,-2.0388,-6.0089,-7.1995,-35.296542 +units=m +to +proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m < NomeFileCASSINI.txt > NomeFileUTMWGS84.txt

Per la trasformazione inversa, utilizzare la stringa modificata come in appresso:
cs2cs -I +proj=cass +proj=cass +lat_0=37.41457777777778 +lon_0=13.780637500000 +x_0=0.0 +y_0=0.0 +ellps=bessel  +towgs84=763.901,183.452,519.534,-2.0388,-6.0089,-7.1995,-35.296542 +units=m +to +proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m < NomeFileUTMWGS84.txt > NomeFileCASSINI.txt
Da notare che sempre il nome del file in input viene inserito prima di quello del file in output

Per convertire un vettoriale (shapefile):

ogr2ogr -s_srs "+proj=cass +lat_0=37.41457777777778 +lon_0=13.780637500000 +x_0=0.0 +y_0=0.0 +ellps=bessel  +towgs84=763.901,183.452,519.534,-2.0388,-6.0089,-7.1995,-35.296542 +units=m" -t_srs "+proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m"  OutputUTMWGS84.shp InputCassini.shp

Per convertire un raster tiff:

gdalwarp -s_srs "+proj=cass +lat_0=37.41457777777778 +lon_0=13.780637500000 +x_0=0.0 +y_0=0.0 +ellps=bessel  +towgs84=763.901,183.452,519.534,-2.0388,-6.0089,-7.1995,-35.296542 +units=m" -t_srs "+proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m" NomeFileCassini.tif NomeFileUTMWGS84.tif