Iniziative dal Forum > Aste, Mercatino e Notizie Utili
Mappe catastali in dxf su mappe di google con qgis usando i punti fiduciali
SC:
visualizzare le mappe catastali sovrapponendole alle immagini satellitari può essere una cosa utile.
non sono un geometra ne un ingegnere non conosco dunque nulla di topografia, mi sono applicato solo alla di una soluzione pratica dell'uso di mappe catastali all'interno di un programma di gis, non sono responsabile in alcun modo ne dell'uso che si fa di queste informazioni ne della loro precisione
purtroppo il casto italiano usa un sistema di coordinate non geografico del secolo scorso il cosiddetto cassini-soldner ( EPSG 9806 le cui trasformazioni sono descritte qui paragrafo 1.3.4)
ogni sistema di riferimento ha bisogno di un punto di origine, il catasto ne ha 849 che coprono il territorio italiano 31 Grandi estensioni e 818 piccoli , i fogli catastali e i punti fiduciali devono riferirsi ad una origine specifiche per avere senso
per esempio che ho usato mi sono riferito al mio paese palazzolo acreide , ho usato un foglio di mappa 72b in formato dxf di qualche anno fa, ho usato anche il punto fiduciale 10 presente nella stesso foglio.
palazzolo fa riferimento come grande estensione a Mineo
i punti fiduciali e le grandi estensioni possono essere presi dal sito www.fiduciali.it le grandi estensioni sono simboleggiati da un triangolo rosso con un punto nero al centro, presenta delle informazioni + utili ai nostri scopi in particolare per esempio per mineo
Longitudine WGS84 14.692473
Latitudine WGS84 37.267029
questo ci da la possibilità di usare la grande estensione come punto di origine anche se come vedremo il risultato sarà meno preciso
il requisiti per fare questa procedura sono
connessione ad internet
qgis programma opensource gratuito (www.qgis.org)
il plugin openlayer presente nel menù plugin di qgis, se non presente bisogno scaricarlo e installarlo attraverso la prima voce "gestisci e installa plugin"
un foglio di mappa in formato dxf, ne allego uno che mi è stato dato da un amico ingegnere , so che si possono comprare e/o richiedere al catasto
procedura usando come origine la grande estensione (mineo in questo caso)
aprire qgis
1) impostare un SR personalizzato
• andare sul menù impostazioni-> sr personalizzato
• aggiungere un nuovo sr , dargli il nome che preferiamo
• scrivere nella apposito riquadro
+proj=cass +lat_0=37.267029 +lon_0=14.692473 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs
• dove i valori 37,26.. e 14,69.. indicati in lat_0 e lon_0 corrispondono ai valori di latitudine e longitudine wgs84 della grande estensione MINEO , punto di origine per i fogli del mio paese, i valori li ho ricavati dal sito www.fiduciali.it
2) menù layer-> aggiungi vettoriale
• scegliere il file dxf del foglio di mappa
• chiederà di scegliere il sistema di riferimento, fra quelli definiti dall'utente in basso scegliere quello che abbiamo creato prima facendo attenzione che i valori dei parametri siano quelli che avevamo dato.
• importiamo sia punti che linee (le linee conterranno la mappa, i punti le scritture con il numero) scegliendo sempre il sistema di riferimento che avevamo creato per ambedue
• vi comparirà la mappa catastale con i punti all'interno , cliccando sulla casellina a sinistra del livello dei punti essi spariranno dalla visualizzazione
3) abilitare proiezione a volo e cambiare il sistema di riferimento in wgs84 pseudo marcato (quello di google)
• andare sul menù file >proprietà progetto-> SR
• selezionare la casellina in alto abilita riproeizione al volo, fare applica
• scegliere tra i sistemi di riferimento wgs84 pseudo marcator (sfera wgs84 su piano), fare applica
4) aggiungere il livello di google map
• andare sul menù plugin->openalayer->add google hybrid (o quella che preferite)
• vi verrà visualizzata una mappa del mondo , mentre sulla parte sinistra comparirà la casellina del ivello se non lo è spostatela in basso della lista in maniera che funga da base alla mappa catastale
5) zoom all'estensione del layer (livello mappa foglio catastale)
• andate sul livello della mappa del foglio catastale se avete usato un dxf come quello allegato si dovrebbe chiamare entities linestring cliccate il tasto destro del mouse
• rivisualizzerà la mappa del foglio catastale e dovrebbe vedersi in sottofondo la corrispondente parte di google
6) ingrandire e rimpicciolire per ricaricare la mappa di google e vedere i particolari
passaggio 1 PROCEDURA PER USARE I PUNTI FIDUCIALI COME ORIGINE, più precisa
a questo punto la mappa catastale dovrebbe essere sovrapposta ma probabilmente ci sarà un grosso livello di imprecisione dovuto alla distanza dall'origine , una soluzione che ho utilizzato io è stata quella di usare i punti fiduciali come punto di origine
per ogni punto fiduciale l'agenzia delle entrate rilascia una monografia,
http://wwwt.agenziaentrate.gov.it/servizi/Monografie/ricerca.php
ossia un pdf con le coordinate cassini , la foto, la posizione nella mappa e altre notizie mancano almeno nella schede che ho trovato le coordinate geografiche in wgs84 importanti per creare il sistema di riferimento in qgis nel passaggio 1
allego quello da me usato
le possiamo ricavare se usiamo google earth e individuiamo il punto fiduciale li, la foto è parecchio utile in tal senso
in questo modo molto empirico io ho ricavato per il punto fiduciale 10 del foglio 72b di palazzolo acreide le seguenti coordinate in wgs 84 lat. 37.063174 lon 14.90616
ora non essendo l'origine effettiva del sistema di riferimento cassini dobbiamo effettuare una correzione nel comando di proiezione attraverso i paramenti del falso est e falso nord x_0 , y_0 , per effettuarla possiamo usare i dati delle coordinate presenti nella monografia
x-22607.2
y 19005.3
invertendo x con y perchè sembra che cassini funzioni invertito rispetto alle coordinate tradizionali in pratica il comando da inserire sarà
+proj=cass +lat_0=37.063174 +lon_0=14.90616 +x_0=19005.3 +y_0=-22607.2 +ellps=WGS84 +units=m +no_defs
usando questa otterrete un passo avanti notevole nel posizionamento della mappa catastale
non escludo che ci possano essere ulteriori raffinamenti per migliorare la precisione , l'uso di altri parametri del comando proj (per esempio ellps e datum forse potrebbero essere utili), fattori di correzione, uso di pf migliori sono aperto a qualsiasi suggerimento , spero anzi che questo post sproni a ulteriori miglioramenti
sono stati di ispirazioni siti come:
http://wiki.gfoss.it/index.php/Catasto
punti fiduciali documentazione
http://geomatica.como.polimi.it/corsi/catasto/pf.pdf
P.S.: ho notato che i dxf forniti dal catasto non sono immediatamente utilizzabili da QGIS in quanto mancano dell'header un trucco può essere riconvertirli sempre in dxf tramite questo programma
http://www.opendesign.com/guestfiles/teighafileconverter
PS2 *solo caso sicilia* se si vuole usare al posto di google le ortofosiciliane (menu layer voce aggiungi layer WMS) usare la seguente stringa:
http://map.sitr.regione.sicilia.it/ArcGIS/services/WGS84_F33/Ortofoto_ATA20072008_f33/MapServer/WMSServer
FI:
Scusa salvo, ma guarda che la planimetria dxf da te usata è una scansione, rasterizzata, del foglio cartaceo quindi le imprecisioni sono inevitabili.
SC:
--- Citazione da: FI - 16:10:27 pm, 06 Novembre 2013 ---Scusa salvo, ma guarda che la planimetria dxf da te usata è una scansione, rasterizzata, del foglio cartaceo quindi le imprecisioni sono inevitabili.
--- Termina citazione ---
il file dxf è uno di quelli del catasto non so come l'hanno ottenuto, ad ogni modo quell'imprecisione non mi sembra dovuta ad un eventuale scansione visto che mi sembra costante,
tu che dici?
FI:
Considera che l'originale di quelle mappe sono state disegnate a mano e che nella scala 1:2000 ogni mm corrisponde a 2 metri quindi se consideri gli strumenti da disegno degli anni 60/70 ti lascio immaginare quale errore possa venire fuori.
SC:
--- Citazione da: FI - 21:03:23 pm, 10 Novembre 2013 ---Considera che l'originale di quelle mappe sono state disegnate a mano e che nella scala 1:2000 ogni mm corrisponde a 2 metri quindi se consideri gli strumenti da disegno degli anni 60/70 ti lascio immaginare quale errore possa venire fuori.
--- Termina citazione ---
beh ok comunque quell'errore non lo posso correggere io posso cercare solo di migliorare l'adattamento, e per quanto mi piacerebbe dire che è tutta colpa del catasto penso che ancora la precisione della procedura sia migliorabile purtroppo con le mia attuali conoscenze non riesco ad andare avanti
pertanto ogni consiglio che mi aiuti in questo obiettivo è gradito
Navigazione
[0] Indice dei post
Vai alla versione completa