Premessa
Più di una volta mi sono imbattuto nel problema di dover gestire dati catastali. Questi dati sono spesso forniti con proiezione di Cassini basata su ellissoide di Bessel orientato a Genova. Sono presenti inoltre numerosi centri di emanazione (origini di sistemi di riferimento locali) distribuiti sul territorio.
Basi di partenza
La lettura di questo articolo mi ha invogliato a provare a definire una stringa valida, per definire sistemi di riferimento personalizzati, per Proj4 e quindi per QGIS, GRASS e tutti gli altri software che si basano su Proj4.
Inoltre la disponibilità delle coordinate dei centri di emanazione (http://www.fiduciali.it) mi ha definitivamente convinto a provare.
Il risultato, per una zona del territorio umbro-ternano, è mostrato in figura.
- Un esempio del risultato ottenuto. La base è l’ortofoto UTM WGS84 fuso 33 tratta dal Portale Cartografico Nazionale. In rosso il dato riproiettato con un software specifico messo a disposizione dalla Regione Umbria per il territorio regionale (Ferens), in verde il dato riproiettato utilizzando la stringa di proj4 sotto descritta.
Dettagli
Dalla lettura del testo si possono desumere i 7 parametri di Helmert (basta invertire i termini di rotazione di Bursa-Wolf) per l’orientazione dell’ellisoide di Bessel a Genova:
+towgs84=656.5,138.2,506.5,5.187,-2.540,5.256,–12.61
Sappiamo inoltre che, per la provincia di Terni il centro di emanazione è Monte Mario, le cui coordinate sull’ellissoide di Bessel possono essere ottenute sia su fiduciali.it che sull’articolo sopra citato. Sapendo inoltre che la proiezione da utilizzare è quella di Cassini, si arriva a scrivere la stringa Proj4 per il sistema di riferimento catastale per la provincia di Terni e più in generale per il territorio coperto dal centro di emanazione di M. Mario:
+proj=cass +lat_0=41.923444 +lon_0=12.451900 +x_0=0 +y_0=0 +ellps=bessel +towgs84=656.5,138.2,506.5,5.187,-2.540,5.256,–12.61 +units=m +no_defs
Note a margine
L’articolo riporta anche i parametri per una trasformazione a 3 parametri di Molodensky, che però, almeno per il territorio su cui è stato svolto il test, è risultata meno precisa.
Il test descritto sopra è stato effettuato solo nella zona citata, sono quindi graditi commenti, smentite e test per altre zone italiane (ovviamente cambiando le coordinate del centro di emanazione).
