In sintesi

La mappa stima il profilo della popolazione che risiede in un'area scelta liberamente, sommando i dati delle sezioni di censimento ISTAT 2021 che ricadono al suo interno. Tutti i calcoli avvengono nel browser, in tempo reale, su dati preparati in anticipo: non c'è un server che elabora le richieste e nessun dato dell'utente viene inviato altrove.

1. Preparazione: dati ISTAT, confini, edifici e modello del terreno vengono uniti, riproiettati in WGS84 e convertiti in formati leggeri per il web.
2. Selezione: l'area di analisi (cerchio o poligono) individua le sezioni il cui centroide cade al suo interno.
3. Aggregazione: i conteggi delle sezioni selezionate vengono sommati per ogni argomento.
4. Contesto: la zona viene collocata nella sua UPL, nel quartiere e nella circoscrizione, e descritta dal punto di vista morfologico con il DTM.

Fonti dei dati

Preparazione dei dati

Le operazioni più pesanti sono eseguite una volta sola, fuori dal browser, con strumenti GIS (GDAL/OGR, QGIS e tippecanoe). Il risultato sono file statici che l'applicazione scarica solo quando servono.

• Unione dati–geometrie: le variabili ISTAT sono agganciate alle sezioni tramite il codice univoco SEZ21_ID. A ogni sezione vengono associati la circoscrizione, il quartiere e l'UPL di appartenenza, la superficie e le coordinate del centroide (lon, lat).
• Indicatori: le 3.600 sezioni con tutti i campi usati dai grafici sono raccolte in un unico file JSON (circa 6 MB), caricato all'avvio.
• Tile vettoriali: sezioni, confini amministrativi ed edifici sono convertiti in PMTiles, un archivio di tile che il browser legge a pezzi con richieste HTTP parziali, senza bisogno di un tile server.
• Edifici: ogni edificio riceve il codice della sezione in cui ricade e ne eredita la densità di popolazione e la copertura edificata, così la mappa può colorarli senza calcoli al volo.
• Terreno: dal DTM a 5 m derivano tre prodotti: le tile di elevazione a colori, le tile di rilievo in codifica Terrarium per la vista 3D e una griglia di punti a passo 50 m con gli indici morfologici, distribuita come tile vettoriali.

Unità di analisi

La sezione di censimento è la più piccola ripartizione per cui ISTAT pubblica variabili su popolazione, famiglie, istruzione, occupazione e abitazioni: nel centro storico corrisponde spesso a un isolato. Il Comune di Palermo ne conta 3.600, con 635.439 residenti al 2021.

• 3.079 sezioni con dati di popolazione, di cui 112 con zero residenti (valore valido: aree produttive, verde, infrastrutture).
• 521 sezioni senza dati: il dataset ISTAT non riporta valori di popolazione per queste sezioni.

Come si forma la zona

Una sezione entra nella zona se il suo centroide cade all'interno dell'area di analisi. Non c'è ripartizione proporzionale per superficie: una sezione conta per intero oppure non conta.

• Area circolare: il test è la distanza geodetica (formula dell'emisenoverso, haversine) fra il centro e ciascun centroide, confrontata con il raggio. Il raggio iniziale è 250 m (diametro 500 m), il minimo 50 m.
• Poligono: le sezioni vengono prima filtrate con il rettangolo che contiene il poligono, poi verificate con il test punto-in-poligono (ray casting).
• Aggiornamento continuo: durante il trascinamento il calcolo viene ripetuto ogni 60 ms. Con 3.600 centroidi il confronto uno a uno è immediato e non serve un indice spaziale.

Con raggi piccoli (sotto i 150–200 m) il risultato dipende molto dalla posizione esatta del centro: spostandolo di pochi metri una sezione può entrare o uscire.

Aggregazione degli indicatori

I valori mostrati sono somme dei conteggi delle sezioni incluse. Le percentuali sono calcolate sui totali della zona, non come media delle percentuali di sezione: così ogni residente pesa allo stesso modo, qualunque sia la sezione in cui vive.

• Sezioni senza dati: se la zona ne include, vengono contate e segnalate con un avviso, ma non contribuiscono ai totali. Un avviso frequente indica che la stima è parziale.
• Filtro stranieri: se l'argomento Stranieri è attivo insieme ad altri, i grafici compatibili mostrano solo la popolazione straniera. La piramide passa a 3 fasce larghe (0–14, 15–64, 65+), le uniche che ISTAT incrocia con il sesso per gli stranieri.
• Scala comune: con più argomenti attivi, l'opzione Scala usa lo stesso massimo su tutte le barre, così i valori sono confrontabili in assoluto.
• Confronto A/B: la zona B segue esattamente le stesse regole della zona A, quindi i risultati delle due zone sono direttamente confrontabili.

Densità sulla mappa

Gli edifici sono colorati con l'indicatore della sezione in cui ricadono, calcolato una volta per sezione:

• Densità di popolazione = residenti della sezione ÷ superficie della sezione, in abitanti per ettaro. La scala va da 0 a 400 e oltre.
• Copertura edificata = superficie coperta da edifici ÷ superficie della sezione, in percentuale (0–100%).

Dentro la zona di analisi gli edifici mostrano sempre la densità di popolazione: un edificio è considerato interno se il suo centroide cade nella zona. Poiché tutti gli edifici di una sezione hanno lo stesso colore, la mappa non distingue un palazzo abitato da un magazzino nella stessa sezione.

Classifiche territoriali

La popolazione di 8 circoscrizioni, 25 quartieri e 55 UPL (Unità di Primo Livello) è la somma delle sezioni che vi appartengono. Gli stranieri sono il campo ST1; gli italiani sono ricavati per differenza (P1 − ST1). L'unità evidenziata è quella in cui cade il centro della zona A, individuata sul confine UPL.

Analisi DTM

Il pannello di destra descrive il terreno con i dati dell'HR-DTM-5m dell'IRPI-CNR, un modello del terreno a 5 m di risoluzione ottenuto da rilievi LiDAR integrati con TINITALY. Gli indici sono stati calcolati sul DTM a 5 m e poi campionati su una griglia a passo 50 m.

• Punto più vicino: l'applicazione carica le tile della griglia attorno al centro della zona A (la tile che lo contiene e le otto adiacenti) e mostra il punto con distanza minima. È un valore puntuale, non una media sulla zona.
• Indici: quota, pendenza, esposizione, forma del terreno, stabilità dei versanti, costruibilità, indici morfometrici (TRI, TPI, SRI), idrologici (TWI, SPI, accumulo di flusso), energetici (cielo visibile, potenziale fotovoltaico, ombreggiamento), velocità di cammino secondo Tobler ed erosione (fattore LS della RUSLE).
• Vista 3D: il rilievo è lo stesso DTM, esagerato 1,5 volte per rendere leggibili le pendenze.

Il dettaglio di ogni indice è nella mappa palermo_dtm5m.

Esportazione

Il CSV dei totali riporta, per la zona A (ed eventualmente B), tipo di area, coordinate del centro, superficie, numero di sezioni, sezioni senza dati e tutti gli argomenti aggregati. La superficie è π·r² per il cerchio; per il poligono è calcolata con la formula di Gauss su una proiezione locale. Il CSV delle sezioni elenca i valori ISTAT originali di ogni sezione inclusa. I file usano il punto e virgola come separatore e la virgola decimale, così si aprono direttamente in Excel e LibreOffice.

Limiti

• Data: i dati fotografano la popolazione residente al 2021 e non registrano i cambiamenti successivi.
• Inclusione per centroide: ai margini della zona la popolazione può essere sovrastimata o sottostimata. Le sezioni molto estese (periferia, aree verdi, porto) pesano per intero anche quando l'area di analisi ne copre solo una parte.
• Effetto della forma: aree diverse sullo stesso territorio possono dare risultati diversi (problema dell'unità areale modificabile, MAUP). Conviene provare più raggi e posizioni.
• Dati mancanti: le 521 sezioni senza dati rendono i totali una stima per difetto nelle zone in cui sono concentrate.
• Edifici: il layer combina fonti di anni diversi; forme e altezze sono indicative.

I risultati sono stime a fini esplorativi e divulgativi e non sostituiscono le statistiche ufficiali.