Tutorial QGIS gratuito: mappa della densità di popolazione dei comuni toscani con i dati ISTAT (QGIS 3.40 LTR)

Risorsa educativa aperta. Questo tutorial è gratuito e rilasciato con licenza Creative Commons Attribuzione – Condividi allo stesso modo 4.0 (CC BY-SA 4.0): puoi usarlo, modificarlo e ridistribuirlo, anche nei tuoi corsi, citando SmartLands e mantenendo la stessa licenza. È pensato per chi ha appena iniziato con QGIS e vuole arrivare a una mappa tematica completa, pronta da stampare.

Cosa ottieni alla fine: una mappa in PDF della densità di popolazione (abitanti per km²) dei comuni della Toscana, costruita interamente con dati aperti ISTAT.

Cosa impari: caricare e filtrare un layer vettoriale, lavorare in un sistema di riferimento proiettato, salvare in GeoPackage, calcolare nuovi campi, unire una tabella ai confini con un join, creare una simbologia graduata e impaginare la mappa nel layout di stampa.

Requisiti: QGIS 3.40 LTR (versione a lungo termine, scaricabile gratis da qgis.org). Tempo stimato: 60–90 minuti.

1. Scaricare i dati

Servono due file, entrambi pubblicati da ISTAT come dati aperti:

1. Confini dei comuni al 1° gennaio 2025, versione generalizzata: Limiti01012025_g.zip. Dentro lo zip, la cartella Com01012025_g contiene lo shapefile dei comuni. Il sistema di riferimento è WGS84 / UTM fuso 32N (EPSG:32632).
2. Popolazione residente per comune al 1° gennaio, dal portale demografico demo.istat.it: scegli la tavola della popolazione residente, seleziona la Toscana con dettaglio comunale e scarica il file CSV.

Scompatta tutto in una cartella di lavoro, per esempio tutorial_toscana.

2. Creare il progetto e impostare il sistema di riferimento

1. Apri QGIS e crea un nuovo progetto (Progetto › Nuovo).
2. In basso a destra clicca sul codice del sistema di riferimento e scegli EPSG:32632 – WGS 84 / UTM zone 32N.
3. Salva subito il progetto nella cartella di lavoro (Progetto › Salva con nome).

Perché è importante: in un sistema proiettato in metri, come UTM 32N, le aree si calcolano direttamente in metri quadrati. In un sistema geografico (gradi), un calcolo di superficie fatto male dà numeri senza senso.

3. Caricare i comuni e tenere solo la Toscana

1. Trascina nella mappa il file Com01012025_g_WGS84.shp (il nome esatto può variare leggermente tra le edizioni).
2. Apri la tabella degli attributi (tasto destro sul layer › Apri tabella attributi) e osserva i campi: COD_REG (codice regione), PRO_COM e PRO_COM_T (codice ISTAT del comune, numerico e testuale a 6 cifre), COMUNE (nome).
3. Tasto destro sul layer › Filtro… e scrivi l’espressione:
   "COD_REG" = 9
   Il codice 9 corrisponde alla Toscana. Conferma: in mappa restano solo i comuni toscani (sono 273).

4. Salvare in GeoPackage

1. Tasto destro sul layer filtrato › Esporta › Salva elementi con nome…
2. Formato: GeoPackage; nome file: toscana.gpkg; nome layer: comuni; SR: EPSG:32632.
3. Rimuovi dal progetto lo shapefile originale e lavora solo sul nuovo layer comuni.

Il GeoPackage è il formato consigliato da QGIS: un solo file, nessun limite di 10 caratteri sui nomi dei campi, più layer nello stesso contenitore.

5. Calcolare la superficie in km²

1. Apri la tabella attributi del layer comuni, attiva la modifica (icona della matita) e apri il Calcolatore di campi (icona dell’abaco).
2. Crea un nuovo campo: nome area_kmq, tipo Numero decimale (real).
3. Espressione:
   area($geometry) / 1000000
4. Salva le modifiche e disattiva la matita.

Nota: area($geometry) calcola l’area planimetrica nel sistema del layer (qui metri quadrati UTM). L’alternativa $area usa l’ellissoide impostato nelle proprietà del progetto: i due valori differiscono di poco. Per questo esercizio va bene il primo, basta dichiararlo nella fonte della mappa.

6. Unire la tabella della popolazione (join)

1. Apri Layer › Gestore origine dati › Testo delimitato, scegli il CSV scaricato da demo.istat.it, imposta il delimitatore corretto (di solito il punto e virgola) e seleziona Nessuna geometria. Se il file ha righe di intestazione prima dei nomi di colonna, indica quante righe saltare.
2. Controlla nella tabella del CSV la colonna con il codice del comune e quella con la popolazione totale.
3. Apri le Proprietà del layer comuni › scheda Join › pulsante +.
4. Layer da unire: il CSV. Campo di join: la colonna del codice comune. Campo destinazione: PRO_COM se il codice nel CSV è un numero, PRO_COM_T se è un testo a 6 cifre con lo zero iniziale (per esempio 048017).
5. Spunta solo il campo della popolazione e imposta un prefisso breve, per esempio pop_.

Verifica: apri la tabella attributi. Se il campo della popolazione è vuoto per tutti i comuni, i codici non coincidono (numero contro testo, oppure zero iniziale mancante). È l’errore più comune nei join: controlla i tipi dei due campi prima di tutto.

7. Calcolare la densità

1. Riattiva la modifica sul layer comuni e apri il Calcolatore di campi.
2. Nuovo campo densita, tipo decimale, espressione (adatta il nome del campo a quello del tuo join):
   to_real("pop_Totale") / "area_kmq"
3. Salva le modifiche.

Il campo unito con il join è “virtuale”: se rimuovi il join, sparisce. Il campo densita invece è scritto nel GeoPackage e resta.

8. Simbologia graduata

1. Proprietà › Simbologia: scegli Graduato, valore densita.
2. Modo: Interruzioni naturali (Jenks), 5 classi. Scala di colori: una sequenza chiara-scura, per esempio YlOrRd.
3. Clicca Classifica e poi OK.
4. Arrotonda le etichette della legenda: nella scheda Istogramma o nelle impostazioni del formato delle etichette imposta 0 decimali.

Da osservare: i valori alti si concentrano nella piana Firenze–Prato–Pistoia e lungo la costa, quelli bassi nelle aree interne e montane. Prova anche il metodo Quantile e confronta: la stessa informazione può raccontare storie diverse a seconda delle classi scelte.

9. Layout di stampa e PDF

1. Progetto › Nuovo layout di stampa, nome Densità Toscana.
2. Aggiungi la mappa (Aggiungi elemento › Aggiungi mappa) disegnando un riquadro sulla pagina.
3. Aggiungi legenda, barra di scala (in km) e freccia del nord.
4. Aggiungi un titolo e un’etichetta con la fonte:
   Fonte: ISTAT, confini amministrativi al 1/1/2025 e popolazione residente, licenza CC BY 4.0. Elaborazione: [il tuo nome] con QGIS.
5. Layout › Esporta come PDF.

10. Esercizi di verifica

1. Con Seleziona elementi usando un’espressione, trova i comuni con densità superiore a 1.000 ab/km². Quanti sono?
2. Cambia il filtro del punto 3 per lavorare su un’altra regione (i codici COD_REG vanno da 1 a 20) e rifai la mappa senza ripetere i passaggi: basta riapplicare lo stile salvato (Proprietà › Stile › Salva stile).
3. Aggiungi le etichette solo ai comuni sopra i 50.000 abitanti usando una regola nella scheda Etichette.

Errori frequenti

• Aree enormi o minuscole: il progetto o il layer sono in un sistema geografico (EPSG:4326). Riproietta in EPSG:32632.
• Join vuoto: tipi di campo diversi o zero iniziale mancante nel codice ISTAT.
• Densità “NULL”: il campo della popolazione è stato importato come testo; usa to_real() o to_int() nell’espressione.
• Accenti illeggibili nei nomi dei comuni: imposta la codifica UTF-8 nelle proprietà della sorgente del layer.

Licenza e riuso

Testo del tutorial © SmartLands Srl, rilasciato con licenza CC BY-SA 4.0. I dati sono di ISTAT e restano soggetti alla loro licenza (CC BY 4.0, con obbligo di citare la fonte). QGIS è software libero distribuito con licenza GNU GPL.

Se usi questo materiale in un corso o in un’esercitazione, ci fa piacere saperlo. E se vuoi fare lo stesso esercizio sui dati del tuo ente, con un docente dal vivo, trovi i nostri corsi QGIS.