Comune di Palermo: aggiornamento del catasto dei soprassuoli percorsi dal fuoco
andy
Buon pomeriggio,
sui gruppi Telegram e FB di opendatasicilia si è parlato del tema in oggetto.
Il punto di partenza è questa notizia (e uno dei PDF a fondo pagina):
https://www.comune.palermo.it/palermo-informa-dettaglio.php?tp=4&id=41439
Creo questo post qui per farlo diventare un nodo centrale dove raccogliere idee, problemi, ecc..
Inizio con un contributo sui dati di base, perché leggo qui che scaricare ad esempio le particelle catastali può essere un problema.
La Regione Siciliana esponi i dati catastali qui:
https://map.sitr.regione.sicilia.it/gis/rest/services/catasto/cartografia_catastale/MapServer
Non è definito quando siano aggiornati, ma dovrebbero essere abbastanza aggiornati.
Sono esposti tramite ArcGIS Server REST API, e molti software GIS hanno l’accesso diretto.
Il problema è che sono layer grandi, e il server non è rapido e l’utilizzo visuale è lento.
In queste cose, la riga di comando o un linguaggio di scripting possono essere una soluzione.
Le particelle sono il layer con ID=4:
https://map.sitr.regione.sicilia.it/gis/rest/services/catasto/cartografia_catastale/MapServer/4
Ogni layer è interrogabile e c’è una GUI per testare delle query:
https://map.sitr.regione.sicilia.it/gis/rest/services/catasto/cartografia_catastale/MapServer/4/query
Una query base, per fare test, sempre valida, è quella in cui 1=1.
L’output HTML mostra i primi 1000 record ed è utile per vedere quali sono gli attributi.
C’è il campo PRO_COM_T che è tipicamente il codice ISTAT del Comune. Quello di Palermo è 082053.
Allora si può impostare una query in cui PRO_COM_T='082053'. È molto più lenta, ma restituisce un’anteprima HTML dei primi 1000 record per le particelle del Comune di Palermo.
La query, semplificando al minimo, ha questo URL (metto a capo per leggibilità):
https://map.sitr.regione.sicilia.it/gis/rest/services/catasto/cartografia_catastale/MapServer/4/
query?where=PRO_COM_T%3D%27082053%27&outFields=*&f=json
Con where si imposta la query per Comune, con outFields=* si chiedono in output tutti i campi, e con f=json si imposta il formato di output.
Per questi accessi e questi download, io consiglio e uso sempre GDAL, una libreria disponibile anche a riga di comando su tutti i sistemi operativi.
Il comando (qui sotto con la sintassi Linux) è
ogr2ogr -f GPKG output.gpkg \
"https://map.sitr.regione.sicilia.it/gis/rest/services/catasto/cartografia_catastale/MapServer/4/query?where=PRO_COM_T%3D%27082053%27&outFields=*&f=json&orderByFields=OBJECTID+ASC" ESRIJSON \
-oo FEATURE_SERVER_PAGING=YES \
--debug on \
-nln particelle \
-append
Note:
- con
-f GPKG output.gpkg, si impostano formato e nome file - poi URL per eseguire la query e nome layer di input (qui è
ESRIJSON) - poi
-oo FEATURE_SERVER_PAGING=YESper paginare la richiesta di 1000 in 1000 - poi
--debug onper avere a schermo riscontro di quello che avviene (vedi esempio sotto) - poi
-nln particelleper fissare il nome del layer di output - poi
-appendper non andare in scrittura in “distruzione”, ma in append, ovvero a integrare.
Fatto partire il comando, potete scendere a comprare un gelato e mangiarvelo.
Al ritorno avrete le particelle catastali. Le allego in 082053.7z.
C’è da fare il JOIN con le tabelle del PDF (che allego). E prima bisogna estrarre le tabelle dal PDF (con tabula, excel pymupd e decine di altri modi). Allego il CSV delle tabelle estratte dal PDF.
Qualche nota:
- di alcune coppie foglio particella, non mi sembra ci sia traccia nei dati del catasto citato sopra (sono quelle del file
not_found.csvallegato) - il file
_06082024055736.gpkgè file spaziale che deriva dal JOIN dei dati del PDF con i dati catastali - Gianni Vitrano ha verificato che nei dati del SIF ci sono più particella percorse dal fuoco. Non appena possibile sarebbero da estrarre e segnalare al Comune e/o al SIF
- ci sono nel PDF alcuni casi in cui la stessa coppia foglio particella appare più di una volta. Vuol dire che c’è stato più di un incendio?
Per le query ho usato duckdb e l’estensione spaziale. Scrivo la query principale, ma non entro nei dettagli.
duckdb -c "COPY (SELECT ST_AsWKB(A.geom) geom,B.* FROM st_read('082053.gpkg',layer=particelle) A
JOIN read_csv('_06082024055736.csv') B
ON A.FOGLIO = B.n_foglio AND A.NUMERO = B.n_particella)
to '_06082024055736.gpkg'
WITH (FORMAT GDAL, DRIVER 'GPKG',SRS 'EPSG:25833')"
Qui sotto, in rosso, le particelle del catasto dei soprassuoli percorsi dal fuoco.
Non ho fatto molte verifiche, e ci potrebbe essere qualche bruttura.
Ma l’elemento chiave è il file con le particelle, scaricato dal SITR regionale.
Saluti
website: https://medium.com/tantotanto
38° 7' 48" N, 13° 21' 9" E, EPSG:4326
___________________
"cercare e saper riconoscere chi e cosa,
in mezzo all’inferno, non è inferno,
e farlo durare, e dargli spazio"
Italo Calvino
Martina Oddo
Ho fatto io le richieste su facebook, e ringrazio anche quì chi ha dato il proprio contributo.
Mi trovo costretta a dividere la risposta perché la mailing list mi dice che il messaggio che ho scritto è troppo lungo.
Introduco la questione:
Faccio parte dell'associazione Fenice Verde che si occupa di incendi boschivi, ed in particolare seguo l'applicazione della legge quadro del 2000 e successive modifiche.
Questa normativa richiede ai comuni la pubblicazione -in formato digitale dal 2021- delle particelle bruciate durante l'anno precedente, al fine dell'applicazione di una serie di vincoli di protezione ambientale.
Secondo le norme, i comuni POSSONO AVVALERSI dei dati forniti dal Corpo Forestale (SIF), che non è mai completo, in quanto riporta solo gli eventi nei quali è intervenuto il servizio AIB forestale.
Di fatti, nella nostra rubrica fb "catasto" abbiamo evidenziato molti casi in cui il SIF ha pubblicato poligoni tagliati "con l'accetta", ad esempio secondo i quali il fuoco avrebbe seguito linee rette su versanti frastagliati.
Ma poco male, la legge è chiara quei dati regionali sono solo indicativi.
I comuni devono quindi aggiungere tutti gli incendi spenti dai vigili del fuoco, dalla protezione civile, dai cittadini o dalla "natura" (termine del combustibile, pioggia ecc...) e disegnare nuovi poligoni.
Nei successivi 30 giorni dalla pubblicazione la legge dà alla cittadinanza la possibilità di fare osservazioni /correzioni al "catasto delle aree percorse da fuoco".
Situazione al 5 luglio:
Il comune di Palermo ha pubblicato un elenco testuale (senza cartografia digitale allegata) con 5 giorni di ritardo rispetto alla scadenza del 31 luglio sancita dalla legge regionale n.13 del 2022, e lo ha fatto sotto le nostre pressioni: abbiamo portato avanti una campagna di partecipazione della cittadinanza chiamata "scrivi al sindaco" e ci siamo presentati al consiglio comunale l'1 agosto per interloquire con i consiglieri.
Questo elenco testuale identifica 20 macroaree bruciate, che (a me) sembrano ritagliati usando i limiti comunali dallo SHP del SIF.
Spero che non sia così, perché sarebbe imbarazzante metterci mesi cinque mesi a fare un ritaglio (il database SIF è online da metà febbraio!!!).
(continua)
Martina Oddo
Situazione attuale:
Dopo la pubblicazione abbiamo coinvolto dei volontari (di cui nessuno esperto GIS e pochi abili al computer) per farci aiutare a stilare le osservazioni.
Sicuramente mancano alcuni eventi/macroaree, di cui abbiamo la certezza perché fanno parte di una relazione scritta dai Rangers d'Italia che gestiscono la RNO di Monte Pellegrino.
Inoltre aggiungeremo la tabella "not found" abilmente ottenuta da @aborruso.
Nel frattempo sto cercando di creare dei layout da stampare in A4 ed affidare ai vari volontari, in modo che possano segnalare eventuali errori "a matita".
Stiamo scaricando dal Copernicus Browser le immagini di Sentinel 2 scattate qualche giorno prima e dopo ogni evento, ritagliano per estensione poco più grande della macroarea identificata dal Comune.
Sto valutando se si vede meglio la true color, la false color o NDVI (questo indice è il migliore per identificare la vegetazione bruciata).
Vi allego quelle della macroarea Monte Grifone 2, se mi sapete dire come scaricarle a migliore risoluzione, o se è più affidabile un altro satellite, vi ringrazio.
Se si potessero confrontare automaticamente le due date sarebbe un grande passo avanti. Io purtroppo sono autodidatta e non so processare i raster.
(continua)
Martina Oddo
La tabella è l'incrocio tra le particelle AGEA e quelle elencate dal Comune il 5 luglio.
In questo shp visualizzeremo con simboli diversi le particelle bruciate con estensione parziale e totale, e la colonna f/p come etichetta.
Lateralmente alla mappa metterò la tabella degli attributi visibili, con le colonne f/p; data osservazione; Ext_tot/parz; note.
(continua)
Martina Oddo
Ad esempio se vede che f1p1099 vede aree bruciate, o ha delle prove del fatto che è bruciata, scriverà nella casella Ext_Tot/Parz corrispondente T o P (totale o parziale), più eventuali note.
Poi raccoglieremo i fogli compilati e mi metterò a trasporli a mano su QGIS.
A quel punto stamperò le mappe+tabelle delle osservazioni da mandare al Comune.
Ho già il modello che hanno utilizzato in altri comuni siciliani per la parte burocratica.
Spero che questa spiegazione sia abbastanza chiara.
Come dicevo sono autodidatta di QGIS, probabilmente ho sbagliato qualche termine, sicuramente uso procedure elementari e sono all'oscuro della maggior parte delle capacità del software.
D'altra parte ho cercato un metodo per coinvolgere attivamente i volontari che si sono offerti per fare le osservazioni del Catasto.
Sarò felice di imparare dai vostri consigli.
Martina Oddo
(fine)
Giovan Battista Vitrano
Martina Oddo
--
Sito: http://opendatasicilia.it
Facebook: https://www.facebook.com/groups/opendatasicilia/
twitter: http://twitter.com/opendatasicilia
Gruppo Telegram: https://t.me/opendatasicilia
---
Hai ricevuto questo messaggio perché sei iscritto al gruppo "opendatasicilia" di Google Gruppi.
Per annullare l'iscrizione a questo gruppo e non ricevere più le sue email, invia un'email a opendatasicil...@googlegroups.com.
Per visualizzare questa discussione sul Web, visita https://groups.google.com/d/msgid/opendatasicilia/96439c66-2b39-4290-a9ef-468831db3ddan%40googlegroups.com.