tarantorespira

verde · il metodo

dati, formule, pesi, sensibilità, limiti

La mappa del verde risponde a una domanda sola: in quale ettaro di Taranto un nuovo parco servirebbe di più? Non dice dove si può fare — proprietà, bonifiche e costi non sono nei dati — ma dove il bisogno è più alto, con numeri che chiunque può rifare. Questa pagina è il metodo per intero; la versione testuale è servita tale e quale in dossier/metodo-verde.md.

in breve, per tutti

Dividiamo il comune in quadrati di 100 metri di lato: un ettaro ciascuno. Per ogni ettaro misuriamo quattro cose: quanto è caldo il suolo d'estate rispetto alla media della città (satellite Landsat); quante persone vivono entro cinque minuti a piedi (censimento 2021), e quante di loro sono bambini o anziani; quanto verde manca, cioè se non c'è un parco pubblico entro 300 metri e se la vegetazione intorno è poca (satellite Sentinel-2). Poi togliamo gli ettari dove un parco non si può fare: acqua, terreno già occupato da case e strade, aree archeologiche e tutelate, zone già verdi, posti dove non abita nessuno. Agli ettari rimasti diamo un punteggio da 0 a 1 e li ordiniamo. I primi cinquanta sono spiegati uno per uno.

che cosa misuriamo

da dove vengono i dati

datofontedatarisoluzione
temperatura del suoloUSGS Landsat 8/9 Collection 2 Level-2 ST (via Microsoft Planetary Computer)estati 2024–2025–202630 m → 100 m
vegetazione (NDVI)ESA Sentinel-2 L2A Collection 1 (via Earth Search, AWS)estati 2025–2026 (23 passaggi)10 m
copertura del suoloESA WorldCover v200202110 m
residenti per etàISTAT, Censimento permanente, sezioni di censimento2021sezione
confine comunaleISTAT, confini amministrativi1 gennaio 2026—
parchi, edifici, strade, acqua, usi esclusiOpenStreetMap (Overpass)2026-11-09vettoriale
vincoli: PPTR, SIN, beni culturaliRegione Puglia, servizi SIT (PPTR approvato, SIN, Vincoli)2026-11-09vettoriale

come si calcola

chi è candidato

Esclusi, per motivo (un ettaro può contarne più d'uno): acqua 2.394 · nessun residente entro 300 m 20.766 · suolo già occupato 3.507 · vincolo storico-culturale 385 · già verde 8.249 · senza temperatura estiva 0. Candidati: 2.711.

il punteggio

Ogni criterio è normalizzato fra il 5° e il 95° percentile dei candidati (0–1, tagliato agli estremi, così un valore anomalo non schiaccia gli altri), poi combinato con somma pesata:

criteriopesop5 → p95 (valori grezzi)
caldo (°C sopra la media)35 %-3,38 → 3,81
persone entro 300 m25 %141,49 → 6727,61
bambini e anziani entro 300 m15 %48,55 → 2588,55
verde che manca (0–1)25 %0,27 → 0,89

I cinque colori della mappa sono i quintili del punteggio fra i candidati: ogni classe contiene un quinto degli ettari.

quanto regge se cambiano i pesi

I pesi sono scelte. Per questo il calcolo viene rifatto con tre pesature diverse e si misura quanto cambia la graduatoria:

pesaturapesi (caldo / persone / fragili / verde)Spearmanprimi 50 in comune
equal0,25 / 0,25 / 0,25 / 0,250,98146 su 50
heat_first0,5 / 0,2 / 0,1 / 0,20,96744 su 50
people_first0,2 / 0,4 / 0,2 / 0,20,9445 su 50

39 ettari restano fra i primi cinquanta con tutte e quattro le pesature.

i siti strategici

La mappa del bisogno dice dove il bisogno è più alto; i siti dicono dove un parco cambierebbe di più, tenendo conto che parchi vicini si contendono le stesse persone. È il problema della copertura massima (Church e ReVelle, 1974), risolto in modo trasparente: si sceglie il sito che porta più residenti non ancora serviti entro 300 m da un'area verde di almeno mezzo ettaro, si segnano come serviti, si ripete. Un sito è un ettaro candidato con almeno il 50 % di suolo libero (il parco è ≥ 0,5 ha); ci si ferma a 20 siti o quando il prossimo servirebbe meno di 300 persone nuove. A parità, vince il sito più caldo. Il calcolo viene rifatto massimizzando bambini e anziani invece di tutti i residenti: 10 dei primi 10 siti hanno un sito dell'altra lista entro 300 m.

Oggi 49.875 tarantini su 189.461 (26 %) hanno un parco pubblico di almeno mezzo ettaro entro 300 metri da casa. Con 10 parchi nei punti giusti sarebbero 126.316 (67 %). Con tutti i 20 siti: 146.089 (77 %). Ettari ammessi come sito: 2.489.

n.quartiereviapersone nuove a 5 mindi cui bambini e anzianiserviti, cumulato
1Tre CarrareVia D'Alo Alfieri14.5865.29534 %
2MontegranaroVia Lombardia12.5475.29341 %
3Borgo UmbertinoVia Felice Cavallotti9.4973.30846 %
4Solito-CorviseaVia Plateja9.4473.82551 %
5TamburiVia Torquato Tasso7.9422.90755 %
6ItaliaVia Umbria6.2962.55858 %
7TalsanoVia Sanguzza4.5241.64961 %
8MontegranaroVia Lago di Varano4.4011.85063 %
9TamburiVia Le Citrezze4.0521.43665 %
10Tre CarrareVicolo Vincenzo Cuoco3.1481.17067 %
11SalinellaVia Lago di Montepulciano2.9211.11568 %
12LamaVia Brigantini2.42976370 %
13TalsanoVia Michelangelo Caravaggio2.34482571 %
14SalinellaVia Ancona2.20694672 %
15Quartiere Paolo VIViale Pietro Nenni2.01873773 %
16Quartiere Paolo VIViale Venticinque Aprile1.82058574 %
17TalsanoVia Eugenio Montale1.55153375 %
18LamaVia Begonie1.50561776 %
19Quartiere Paolo VIvia Michele Delvecchio1.49854876 %
20TalsanoVia Giuseppe Abbati1.48151877 %

Un sito è un'ipotesi che dice «qui un parco raggiungerebbe tante persone»: non dice che il suolo sia disponibile. Il perimetro SIN e i vincoli paesaggistici sono segnalati sito per sito.

il tracciatore e il modello di raffrescamento

Sulla mappa si può disegnare un parco, scegliere quanto del suo suolo sarà vegetato a maturità, e leggere che cosa ne prevede un modello costruito sui dati di Taranto: di quanto scenderebbe la temperatura estiva del suolo dentro e intorno, con un intervallo, e quante persone avrebbero un parco a cinque minuti. Il modello non è preso dalla letteratura: è una regressione lineare della temperatura del suolo (mediana estiva Landsat) di ogni ettaro abitato su ciò che l'ettaro contiene e su ciò che ha intorno.

variabile (per unità, 0–1)effetto, °Cintervallo 90 %
quota vegetata dell'ettaro-2,74-2,91 → -2,58
quota vegetata entro 300 m-3,96-4,85 → -3,06
quota dell'ettaro dentro un parco-0,82-1,67 → -0,18
quota di parco entro 300 m-0,11-1,40 → +1,10
quota di edifici-4,88-6,12 → -3,80
quota di strade-4,11-4,96 → -3,32
quota d'acqua entro 300 m-12,54-13,96 → -10,87
log(1 + metri dal mare)+0,88+0,70 → +1,09

In parole: ogni 10 punti in più di suolo vegetato nell'ettaro valgono circa -0,274 °C, ogni 10 punti in più nei 300 metri intorno -0,396 °C; un ettaro interamente dentro un parco aggiunge -0,818 °C oltre il suo verde. Celle usate: 18.311. Adattamento: R² 60 % in campione, 58 % fuori campione (validazione incrociata spaziale a 5 blocchi da 1000 m), errore tipico 1,943 °C per ettaro. Gli intervalli vengono da un bootstrap a blocchi spaziali (500 rifit su blocchi di 1000 m), non dagli errori standard classici: i residui sono fortemente autocorrelati (I di Moran 0,801) e gli errori classici sarebbero da tre a nove volte troppo stretti. Un termine quadratico sul verde è stato provato e scartato (errore fuori campione 1,942 contro 1,943).

Verifica sui parchi esistenti (temperatura osservata contro prevista nelle loro celle; PCI = quanto il parco è più fresco dell'anello 200–500 m non vegetato):

parcohavegetatoosservataprevistascartoPCI
Oasi regionale WWF palude La Vela121,848 %38,7341,02-2,298,43
Parco Cimino10,699 %38,4539,49-1,049,42
Parco del Mirto8,775 %45,945,690,222,8
Parco Archeologico delle Mura Greche "Massimo Battista"8,542 %47,746,631,07-1,5
Giardini Peripato4,582 %39,8541,94-2,093,68
Giardini Antonio De Curtis1,696 %43,6445,84-2,192,61
Giardini Virgilio1,545 %44,7745,16-0,390,6
Giardini Cordelia Jannelli1,362 %45,745,290,41,73
(senza nome)1,36 %48,7847,581,2-0,35
Giardini Garibaldi1,280 %42,0844,49-2,41-0,66
Parco della Gioia di Mister Sorriso0,96 %46,8846,840,04-1,54
Parco Sant Egidio0,867 %45,2745,61-0,342,03
Giardini Cremonesi0,698 %44,2944,84-0,551,45

Scarto medio nei parchi: -0,64 °C (± 1,22): il modello resta prudente, i parchi veri sono un poco più freschi di quanto preveda.

Il tracciatore applica il modello così: il poligono disegnato viene proiettato in UTM 33 N e ritagliato su ogni ettaro; in ogni ettaro toccato la quota vegetata diventa la media pesata fra quella attuale e quella scelta, e la quota di parco cresce della frazione coperta; per ogni ettaro entro 300 metri si ricalcolano le quote del vicinato; la variazione di temperatura è la somma dei termini del modello, ripetuta su ogni rifit del bootstrap per dare il 5° e il 95° percentile. Le persone sono i residenti degli ettari il cui centro dista meno di 300 metri da un ettaro toccato; «oggi senza parco» sono quelli senza un'area verde pubblica di almeno mezzo ettaro a quella distanza, se il disegno raggiunge il mezzo ettaro.

Limiti: è la temperatura della superficie a metà mattina d'estate, non l'aria; il modello è correlazionale (dove c'è più verde fa più fresco) e trasferirlo a un parco nuovo è un'estrapolazione che i parchi esistenti sostengono ma non dimostrano; la quota vegetata «a maturità» di alberi nuovi arriva dopo dieci o vent'anni; il suolo deve essere davvero libero.

i parchi, misurati

Ogni area verde pubblica di almeno mezzo ettaro con la quota di suolo che il satellite vede vegetata d'estate (mediana 76 %):

parcotipoettarivegetato
Oasi regionale WWF palude La Velanature_reserve121,848 %
Pineta Comunale Frantellapark42,199 %
Gravina di Riggiopark23,979 %
Parco Ciminopark10,699 %
Parco del Mirtopark8,775 %
Parco Archeologico delle Mura Greche "Massimo Battista"park8,542 %
Giardini Peripatopark4,582 %
(senza nome)park3,285 %
Parco Madre Teresa di Calcuttapark2,395 %
Giardini Antonio De Curtispark1,696 %
Parco Rimembranzapark1,597 %
Giardini Virgiliopark1,545 %
Giardini Cordelia Jannellipark1,362 %
(senza nome)park1,36 %
Cacacepark1,2—
Giardini Garibaldipark1,280 %
Parco della Musica Giandomenico Caramiaplayground1,149 %
(senza nome)park1,180 %
Parco della Civiltàpark1,083 %
(senza nome)recreation_ground1,066 %
8KBGpark1,092 %
Parco della Gioia di Mister Sorrisopark0,96 %
Parco UNICEFpark0,956 %
(senza nome)park0,98 %
Parco Sant Egidiopark0,867 %
Parco Madre Teresa di Calcuttapark0,8—
Villa Comunalepark0,793 %
Giardini Enrico Berlinguerpark0,72 %
(senza nome)park0,779 %
(senza nome)garden0,777 %
(senza nome)park0,75 %
Parco comunale Schiavonipark0,6—
Giardini Cremonesipark0,698 %
Piazza Giubileopark0,612 %
Parco "Monte Della Croce"park0,643 %
(senza nome)garden0,691 %
(senza nome)park0,688 %
Giardino Maresciallo Francesco Mignozzipark0,511 %
Parco Aldo Moropark0,57 %

limiti, dichiarati

riferimenti

  1. Norton B.A., Coutts A.M., Livesley S.J., Harris R.J., Hunter A.M., Williams N.S.G. (2015). Planning for cooler cities: A framework to prioritise green infrastructure to mitigate high temperatures in urban landscapes. Landscape and Urban Planning 134, 127–138. doi:10.1016/j.landurbplan.2014.10.018
  2. Bowler D.E., Buyung-Ali L., Knight T.M., Pullin A.S. (2010). Urban greening to cool towns and cities: A systematic review of the empirical evidence. Landscape and Urban Planning 97, 147–155. doi:10.1016/j.landurbplan.2010.05.006
  3. Reid C.E. et al. (2009). Mapping Community Determinants of Heat Vulnerability. Environmental Health Perspectives 117(11), 1730–1736. doi:10.1289/ehp.0900683
  4. Konijnendijk C.C. (2023). Evidence-based guidelines for greener, healthier, more resilient neighbourhoods: Introducing the 3–30–300 rule. Journal of Forestry Research 34, 821–830. doi:10.1007/s11676-022-01523-z
  5. WHO Regional Office for Europe (2016). Urban green spaces and health: a review of evidence. Copenhagen.
  6. Malakar N.K., Hulley G.C., Hook S.J., Laraby K., Cook M., Schott J.R. (2018). An Operational Land Surface Temperature Product for Landsat Thermal Data: Methodology and Validation. IEEE Transactions on Geoscience and Remote Sensing 56(10), 5717–5735. doi:10.1109/TGRS.2018.2824828
  7. Cook M., Schott J.R., Mandel J., Raqueno N. (2014). Development of an Operational Calibration Methodology for the Landsat Thermal Data Archive and Initial Testing of the Atmospheric Compensation Component of a Land Surface Temperature (LST) Product. Remote Sensing 6(11), 11244–11266. doi:10.3390/rs61111244
  8. Tucker C.J. (1979). Red and photographic infrared linear combinations for monitoring vegetation. Remote Sensing of Environment 8(2), 127–150. doi:10.1016/0034-4257(79)90013-0
  9. Voogt J.A., Oke T.R. (2003). Thermal remote sensing of urban climates. Remote Sensing of Environment 86(3), 370–384. doi:10.1016/S0034-4257(03)00079-8
  10. Oke T.R. (1982). The energetic basis of the urban heat island. Quarterly Journal of the Royal Meteorological Society 108(455), 1–24. doi:10.1002/qj.49710845502
  11. Zardo L., Geneletti D., Pérez-Soba M., Van Eupen M. (2017). Estimating the cooling capacity of green infrastructures to support urban planning. Ecosystem Services 26, 225–235. doi:10.1016/j.ecoser.2017.06.016
  12. Aram F., Higueras García E., Solgi E., Mansournia S. (2019). Urban green space cooling effect in cities. Heliyon 5(4), e01339. doi:10.1016/j.heliyon.2019.e01339
  13. Malczewski J. (2006). GIS-based multicriteria decision analysis: a survey of the literature. International Journal of Geographical Information Science 20(7), 703–726. doi:10.1080/13658810600661508
  14. Roberts D.R. et al. (2017). Cross-validation strategies for data with temporal, spatial, hierarchical, or phylogenetic structure. Ecography 40(8), 913–929. doi:10.1111/ecog.02881
  15. Moran P.A.P. (1950). Notes on continuous stochastic phenomena. Biometrika 37(1–2), 17–23. doi:10.1093/biomet/37.1-2.17
  16. Snyder J.P. (1987). Map Projections: A Working Manual. USGS Professional Paper 1395. doi:10.3133/pp1395
  17. Church R.L., ReVelle C.S. (1974). The maximal covering location problem. Papers of the Regional Science Association 32, 101–118. doi:10.1007/BF01942293
  18. Chen Y., Yu J., Khan S. (2010). Spatial sensitivity analysis of multi-criteria weights in GIS-based land suitability evaluation. Environmental Modelling & Software 25(12), 1582–1591. doi:10.1016/j.envsoft.2010.06.001
  19. Zanaga D. et al. (2022). ESA WorldCover 10 m 2021 v200. doi:10.5281/zenodo.7254221
  20. D.M. 2 aprile 1968, n. 1444, art. 3: standard di 9 m² di verde pubblico per abitante.

aggiornamenti: numeri che non invecchiano

Numeri calcolati il 9 novembre 2026 da 55 passaggi Landsat e 23 Sentinel-2.

Il server li ricalcola da sé: prossimo controllo di nuovi passaggi dal 7 dicembre 2026.

  • Primo calcolo su questo server: da qui in poi ogni ricalcolo dirà che cosa è cambiato.

Questa mappa non è una fotografia scattata una volta: il server che la pubblica la ricalcola da solo. Ogni mattina alle 5:30 controlla se è il momento (ogni quattro settimane; ogni due da giugno a settembre, quando arrivano i passaggi che contano) e, se lo è, interroga i cataloghi satellitari. I passaggi Landsat mai visti vengono letti e aggregati, quelli già letti sono conservati; se la selezione Sentinel-2 è cambiata, la maschera di vegetazione a 10 m viene rifatta per intero; OpenStreetMap si riscarica ogni 30 giorni e i vincoli regionali ogni 90; il censimento e la copertura del suolo 2021 restano finché ISTAT ed ESA non pubblicano un'edizione nuova. Poi rifà sempre i passi finali: accesso ai parchi, punteggio, siti strategici, modello di raffrescamento, questo documento. La finestra estiva scorre con il calendario: a giugno la nuova estate entra e si riempie passaggio dopo passaggio, l'anno più vecchio esce solo quando un'estate intera l'ha sostituito.

I nuovi file sostituiscono i vecchi in un colpo solo, e solo se tutti i passi sono riusciti: un errore lascia in linea i numeri precedenti e il tentativo si ripete il giorno dopo. Ogni ricalcolo scrive che cosa è cambiato (aggiornamenti.json: nuovi passaggi, ettari che hanno cambiato classe, siti confermati, parchi aggiunti o tolti, coefficienti del modello con i loro intervalli) e una riga di storia (storico.json). Nessun numero viene mai corretto a mano.

calcoloLandsatSentinel-2ettari candidatiserviti oggidopo 10 parchiparchieffetto vegetazioneR² CVclassi cambiate
2026-11-0955232.71126 %67 %39−2,74 °C58 %—

api: la mappa come biblioteca

Tutto ciò che la pagina mostra si può leggere in JSON, senza chiavi né registrazione, con gli stessi numeri e le stesse frasi. I file portano un ETag e restano in cache un'ora.

rifarlo

Il codice è in dossier/analisi/verde/ (Python; ogni passo ha i suoi test). I numeri di ogni ettaro candidato: cells.csv, cells.json; i primi cinquanta con le frasi: top.json; parametri, soglie e sensibilità: priority.json; i parchi misurati: parchi.geojson. Ultimo calcolo: 2026-11-09.

navigazione