Il 24 settembre 2020 è una data storica e rimarrà impressa nella mente di molte persone in quanto le mappe catastali, usufruibili tramite servizio WMS del AdE, diventano utilizzabili grazie alla modifica della licenza da CC BY-NC-ND 4.0 a CC BY 4.0.
Per scaricare i dati dal WMS esistono vari metodi ma, in questo blog post, descriverò solo quello utilizzabile in QGIS, la logica è presa in prestito dal blog post di Andrea Borruso (grazie mille):
QGIS
Metodo GetFeatureInfo
e output in HTML:
Funzione personalizzata da usare nel Calcolatore di campi di QGIS:
il file zippato è scaricabile cliccando qui.
Dove mettere il file
Scaricare lo script WMS_catasto.py
e metterlo nella seguente cartella:
Per Windows:
C:\Users\nomeUtente\AppData\Roaming\QGIS\QGIS3\profiles\default\python\expressions
Per OSX:
/Users/[NAME USER]/Library/Application Support/QGIS/QGIS3/profiles/default/python/expressions
al riavvio di QGIS, trovate la funzione nel calcolatore di campi di QGIS nel gruppo Custom
:
Come usare la funzione personalizzata
La funzione è get_parcel_info()
, come da esempio nella guida in linea, richiede due argomenti (le coordinate X e Y) ovvero è possibile dargli in pasto due campi che contengono le coordinate X e Y.
In questo video è spiegato tutto: nella InfoBox trovate i marcatori di capitolo:
NOTE FINALI: Per favore non abusare di questo servizio, con questo metodo descritto sopra potrai mandare in tilt i server dell AdE e rischiare di far chiudere il servizio.
Un ringraziamento particolare ad Andrea Borruso, che ha pensato, realizzato e condiviso questi metodi.
Nel gist trovate le espressioni utilizzate per i campi virtuali.
RINGRAZIAMENTI
- Agenzia delle Entrate : https://www.agenziaentrate.gov.it/portale/it/web/guest/schede/fabbricatiterreni/consultazione-cartografia-catastale/servizio-consultazione-cartografia
- Andrea Borruso : https://twitter.com/aborruso
- Giulio Fattori : https://github.com/Korto19
RIFERIMENTI
- WMS Catasto : https://wms.cartografia.agenziaentrate.gov.it/inspire/wms/ows01.php
- Blog Post di Andrea Borruso : https://medium.com/tantotanto/le-mappa-castali-diventano-finalmente-utilizzabili-821db2f84533
- Gist : https://gist.github.com/pigreco/86589dddf5a59b3a7650267d5af237bd
Se il blog post vi è piaciuto cliccate su ‘Mi piace’, grazie!!!
if you liked the blog post click on ‘Like’, thank you !!!
Grazie Totò informazioni utilissime. ciao
"Mi piace"Piace a 1 persona
Provato…funziona perfettamente….grazie
Una domanda: come si potrebbe modificare la procedura per estrarre dati di un gruppo di particella senza doverle cliccare una ad una??
"Mi piace"Piace a 1 persona
Ciao,
Basta creare un file con due campi X e Y, li popoli con le coordinate e poi usi la funzione personalizzata.
"Mi piace""Mi piace"
Intanto grazie per il tuo lavoro preziosissimo e fatto veramente bene, anche perchè dici quei piccoli dettagli che ad un neofita sono molto ma molto utili.
Ti scrivo, da neofita per capire 2 cose:
– creare un file, file in excel con 2 colonne e relative coordinate (X,Y o Y,X)?
– funzione personalizzata: quale è, dove si trova?
Grazie ancora 🙂
"Mi piace""Mi piace"
Ciao,
non capisco la prima domanda, la seconda: la funzione personalizzata è get_parcel_info e devi costruirla usando quel codice presente nel blog post.
"Mi piace"Piace a 1 persona
Hai scritto: “Basta creare un file con due campi X e Y, li popoli con le coordinate e poi usi la funzione personalizzata.”
Il fatto e che le coordinate le devo estrapolare da un file KML quindi vorrei sapere se il file da creare posso farlo esterno a qgis e poi importarlo, se cio è possibile, lo creo in excel e lo importo ma non saprei come fare a inserire le coordinate presenti nel file excel dentro la tabella attributi.
grazie
"Mi piace""Mi piace"
Il post lo ho scritto più di un anno fa, ma non trovo scritto nulla di quello che virgoletti.
Se hai un file KML puoi importarlo in QGIS e estrapolare le coordinate.
Capisco che sei neofita, ma questi problemi vanno oltre il post.
Se hai bisogno di una consulenza contattami tramite il mio canale Telegram https://t.me/pigrecoinfinito
"Mi piace"Piace a 1 persona
Grazie delle informazioni chiare e puntuali che ci concedi…
Ho provato la funzione sopra indicata, ed opera egregiamente…DOMANDA: sarebbe possibile ottimizzare la funzione per un estrazione del dato catastale su un ingombro poligonale o lineare?
"Mi piace""Mi piace"
Ciao Mauro, scusa ma non capisco la domanda, potresti fare un esempio?
"Mi piace""Mi piace"
Ciao Totò,
se abbiamo ad esempio l’ingombro di una strada di 3 mt di larghezza e mt 100 di lunghezza che attraversa una serie di particelle, sarebbe possibile estrarre le informazioni dal wms catastale relativamente alle particelle intersecate dal tracciato? non avremmo più quindi dei punti con delle coordinate x-y ma uno shape poligonale o lineare georeferenziato.
pensavo ad una possibile applicazione per le procedure di esproprio ricavando i dati foglio e particella dal wms catastale incrociandoli con i dati .ter ecc. dell’agenzia delle entrate automatizzando in questo modo la ricerca delle proprietà, tipologia della coltura etc.
(ps. attualmente procedo con la rototraslazione dei fogli .cxf convertendoli poi in shape ed intersecandoli con un eventuale tracciato, ma non sempre i fogli .cxf sono disponibili o aggiornati)
"Mi piace"Piace a 1 persona
Ciao, ora è chiaro.
La risposta è sì, laboriosa ma fattibile.
"Mi piace""Mi piace"
Ciao! Stupendo come sempre! Ho un dubbio: si dice di usare EPSG 3045, ma QGIS lo segnala per zona nordest Europa, non per Italia. é un errore nel pannello di QGIS?
"Mi piace""Mi piace"
Ciao Luca,
molto probabilmente sì
"Mi piace""Mi piace"
Buongiorno sono un GeoDataAnalyst e mi sto approcciando da poco alle tematiche catastali. Volevo sapere se c’è modo di estrarre dal WMS le bounding box dei fogli e delle particelle in automatico presenti in una regione o in tutta Italia direttamente? Senza dover mettere io dei punti a definire quali voglio scaricare. Grazie in anticipo
"Mi piace""Mi piace"
Buongiorno Mattew,
che io sappia dal WMS è possibile estrarre i BBOX delle particelle, per i foglio non lo so.
Il WMS però è sottoposto a token con validità di circa 10 minuti e valido per provincia.
saluti
"Mi piace""Mi piace"
Ciao , com’è la sintassi dello script per interrogare il wms ed estrarre la bbox di una particella?
"Mi piace"Piace a 1 persona
Ciao, se sei interessato a questo argomento, contattami in privato.
"Mi piace""Mi piace"
Buongiorno, è possibile tramite una tabella in input, contenente i dati di comune, foglio e mappale risalire alla posizione in modo automatico della particella sulla mappa wms e contrassegnarla con un elemento vettoriale. Mi si pone il problema, in quanto dovrei individuare la posizione di parecchie decine di particelle sulla mappa. Grazie
"Mi piace"Piace a 1 persona
Ciao, non puoi farlo in modo diretto ma esiste altra via; se ti interessa contattami dalla sezione Dove mi trovo
"Mi piace""Mi piace"
ciao, ho provato e riprovato ma non sembra funzionare. il campo catasto non vuol saperne di popolarsi. ho provato con diversi EPSG calcolando il campo con get_parcel_info2( “x” , “y” , ‘EPGS:6706’ ), get_parcel_info2( “x” , “y” , ‘EPGS:3045’ ) oppure get_parcel_info2( “x” , “y” , @project_crs ) (la funzione è copiata tal quale dal tuo tutorial) ma il campo resta vuoto….
"Mi piace""Mi piace"
Ciao,
puoi mandarmi il progetto e i dati che stai usando? cosi faccio un test.
grazie
"Mi piace""Mi piace"
ciao. in questo file zip ho messo lo shape con i punti, il progetto qgis e il file .py con la funzione https://we.tl/t-Q5gDiFVY4M
grazie
"Mi piace""Mi piace"
Eccomi,
ho capito il problema, ma non posso risolverlo per adesso.
La soluzione consigliata è utilizza SR proiettati e non geografici (usa, per esempio EPSG:3045). 🙂
"Mi piace""Mi piace"
ciao, grazie per la risposta. ho testato il tutto con EPGS 3045 come suggerito, impostandolo come sr del progetto e del layer ma purtroppo non ne vuol sapere di funzionare…
"Mi piace""Mi piace"
Scusami, ma EPSG:3045 è quello usato nei test e nel video, quindi funziona, guarda bene i valori che scrivi, per esempi vedo che tu usi EPGS al posto di EPSG, forse è questo uno dei problemi.
"Mi piace""Mi piace"
hai ragione scusami, c’era un refuso, con EPSG 3045 funziona
"Mi piace"Piace a 1 persona
ottimo, buon lavoro
"Mi piace""Mi piace"
Buon pomeriggio,
è possibile partire dal numero delle particella e ottenere la sua estensione senza utilizzare le coordinate di un punto?
"Mi piace"Piace a 1 persona
Ciao Sonia,
spiegati meglio, cosa intendi per estensione?
"Mi piace""Mi piace"
Ciao,
intendo i vertici che delimitano l’area in cui è inserito il poligono (l’area rossa che si crea selezionando la particella catastale dalla geometria).
Grazie
"Mi piace""Mi piace"
Ciao, è possibile farlo ma non direttamente da QGIS ma tramite script.
"Mi piace""Mi piace"
Domanda banale forse, perdono: collegandomi al wms catasto AdE ottengo “vestizioni, particelle, etc” ma non uno shape puntuale “catasto”. Mi manca l’inizio quindi per seguire il ragionamento. Mi dareste una spinta ? Grazie mille !
"Mi piace""Mi piace"
Ciao,
rivedi il workshop da qui:
https://fb.watch/60JOKALWuj/
qui il repository:
https://github.com/pigreco/workshop-estate-gis-2021
"Mi piace""Mi piace"
Buongiorno,
sono una studentessa di GIS al politecnico di Barcellona e non essendo abituata a ricercare dati italiani mi trovo un po’ in difficoltà. Vedo che l’agenzia delle entrate rende ormai facilmente disponibili le geometrie ma non i dati di catasto (oltre a quelli che appaiono nel tutorial). Sarebbe possibile adattare l’esempio usato per le particelle per ottenere informazioni sui ‘buildings’, come ad esempio per ottenere la tipologia?
"Mi piace"Piace a 1 persona
Buongiorno Giulia,
purtroppo il WMS del catasto espone pochi dati, e in questo blog post li estraggo tutti.
ciao
"Mi piace""Mi piace"
Ciao sto cercando di realizzare un progetto per calcolare l’area e il volume di alcuni laghi segnati da punti sulla mappa,
Ho seguito la tua guida che mi ha aiutato moltissimo solo che adesso vorrei inserire all’interno della “Tabella attributi” una voce con l’area intorno ad ogni punto segnato, si intende che l’area la selezioni con la funzione “misura area”, cosa mi consigli di fare?
Grazie
"Mi piace""Mi piace"
Ciao è una domanda fuori tema,
se hai bisogno di consulenza contattami tramite il mio canale Telegram https://t.me/pigrecoinfinito
"Mi piace""Mi piace"
Ciao! Io avrei bisogno di ricercare p.lle catastali come da te indicato nel video presente su Youtube e datato 04/02/2022. Domanda: La procedura è sempre la stessa da te indicata nel video su questo blog?
Grazie mille
"Mi piace""Mi piace"