Utilizzare le interrogazioni spaziali¶

Descrizione dell’esercizio¶

Lavoreremo con 2 dataset - un layer lineare che rappresenta i fiumi e un layer puntuale che descrive le città. L’obiettivo è quello di lanciare un’interrogazione spaziale per trovare tutte le città che non distano più di 10 chilometri da un fiume.

Ottenere i dati necessari¶

Useremo i file ne_10m_rivers_lake_centerlines e 10m_populated_places_simple . Entrambi dataset tratti da Natural Earth.

Download Rivers and Lake Centerlines

Scaricate Populated Places.

Fonte Dati: [NATURALEARTH]

Procedimento¶

1. Dopo che avete scaricato i dati, aprite QGIS. Andate su Layer ‣ Aggiungi Vettore....

2. Fate Click su Sfoglia e portatevi nella cartella dove avete scaricato i file .zip.

3. Trovate i file .zip scaricati. Tenete premuto il tasto Shift e fate click su entrambi i file zip per selezionarli. Click su Apri.

4. Apparirà una finestra di dialogo che vi chiederà di selezionare il layer da aprire tra i file zip. Selezionate ne_10m_rivers_lake_centerlines.shp e fate click su OK.

5. Visto che avete operato una selezione multipla, ripetete l’operazione per il file successivo. Selezionate 10m_populated_places_simple.shp e fate click su OK.

6. Potrete vedere ora entrambi gli shapefile caricati in QGIS.

7. Genereremo dei buffer sia nel layer puntuale che in quello lineare. Lo strumento di geoprecessing per il Buffer in QGIS usa la numerazione metrico decimale ordinaria per calcolare la misura del buffer. Ma i layer con cui stiamo lavorando sono in un Sistema di Riferimento (SR) basato su Coordinate Geografiche che, com’è noto, vengono calcolate in gradi. . Questo non soddisfa le esigenze delle nostre analisi che prevedono invece l’uso di unità di misura come i metri o i chilometri. Per risolvere questo problema dobbiamo allora ri-proiettare i nostri layer su un Sistema di Coordinate Proiettare (SR). Click sul tasto destro del mouse sul layer 10m_populated_places_simple layer e quindi scegliete Salva con Nome.

8. Nella finestra di dialogo Salva i vettori come... scegliete il tasto Sfoglia accanto alla voce Salva con Nome e assegnate al nuovo file il nome di populated_places_reprojected.shp . Ora fate click sul successivo tasto Sfoglia vicino alla voce Sist Rif.

9. Ora dobbiamo scegliere un Sistema di Riferimento (SR) che sia adeguato ai nostri scopi. Per creare dei buffer una proiezione Azimutale Equidistante dovrebbe essere la scelta migliore. Nel nostro caso, considerato che il dataset è globale, sceglieremo una proiezione planetaria. Nella finestra di dialogo Selettore di sistema di riferimento avviate la ricerca digitando world az.. e vedrete comparire immediatamente il risultato. Selezionate World_Azimuthal_Equidistant e fate click su OK.

10. Tornate nella finestra di dialogo Salva i vettori come ... e spuntate la casella Aggiungi il file salvato sulla mappa quindi fate click su OK.

11. Adesso ripetete l’operazione di ri-proiezione per il file ne_10m_rivers_lake_centerlines e salvate il nuovo layer di laghi e fiumi con il nome di rivers_lake_reprojected.shp.

12. Adesso avete 4 layer nel vostro:guilabel:Pannello dei Layer. Deselezionate i layer originali per vedere soltanto i layer sui quali sono state effettuate le riproiezioni. I layer riproiettati sono ancora rappresentati nei CRS di tipo Geografico per via delle impostazioni di default. Vediamo di annullarle. Facciamo click sul tasto Stato del SR in basso a destra. Questa finestra può essere raggiunta anche da menu Progetto ‣ Proprietà del Progetto.

13. Nella scheda Sistema di Riferimento (SR) della finestra di dialogo Proprietà del progetto disattivate la casella Abilita la riproiezione al volo. Click su OK.

14. Torniamo adesso nella finestra principale di QGIS, click su uno dei layer riproiettati e click su ” Zoom all’estensione del Layer.

15. Adesso vedremo finalmente nostri dati nel Sistema di Riferimento che abbiamo scelto per i nostri layer. Adesso dobbiamo creare dei buffer per entrambi i dataset. Click su Vettore ‣ Strumenti di Geoprocessing ‣ Buffer.

16. Nello strumento Buffer scegliete come file di input populated_places_reprojected . Inserite una distanza di buffer pari a 10000. Notate che vogliamo un buffer di 10 chilometri e, considerato che le unità del sistema di riferimento sono i metri, dobbiamo inserire 10.000. Inserite come nome del file in uscita populated_places_buffer.shp. E fate click su OK.

17. Una volta che il processo di buffering è finito, fate click su Yes quando vi verrà chiesto se volete aggiungere il nuovo vettore alla TOC.

18. Ripetete lo stesso processo di buffer per il file rivers_lake_reprojected e create un file di output che chiameremo rivers_lake_buffer.shp.

19. Il problema è che il nuovo layer rivers_lake_buffer contiene geometria che riguardano tanto i fiumi quanto i laghi. La nostra ricerca è interessata esclusivamente alle geometrie dei fiumi, Allora facciamo ricorso ad un’interrogazione che ci permetta di selezionare esclusivamente le geometrie dei fiumi. Click con il tasto destro su rivers_lake_buffer e selezioniamo Apri tabella attributi.

20. Come potete notare gli attributi della colonna “featurecla” contengono informazioni che possiamo usare per separare le geometrie dei fiumi da quelle dei laghi. Fate click sul tasto Seleziona elementi usando un espressione.

21. Inserite l’espressione “featurecla” = “River” fate click su Seleziona e infine click su Chiudi per tornare alla finestra principale di QGIS.

22. Ora siamo pronti per realizzare l’interrogazione spaziale. Dovete aver attivato il Plugin Interrogazione spaziale per utilizzare queste funzionalità. Vedetevi Usare i plugin per maggiori dettagli. Quando il plugin è attivo fate click su Vettore ‣ Interrogazione Spaziale ‣ Interrogazione Spaziale.

23. A noi interessa selezionare con la nostra interrogazione tutte le geometrie bufferizzate delle città che intersecano le linee bufferizzate dei fiumi. Assicuratevi che la casella vicino alla dicitura “geometrie selezionate” sia spuntata e indichi il numero dei fiumi selezionati. Questo ci assicura che la nostra interrogazione interesserà esclusivamente le geometrie dei fiumi che abbiamo precedentemente selezionato. Fate click su Apply.

24. Una volta che la nostra interrogazione è completa vedrete aggiungersi alla finestra di dialogo una nuova scheda chiamata Oggetti Selezionati. Fate click sul pulsante Crea layer con quello selezionato. Un nuovo layer si aggiungerà al pannello dei layer. Fate click su Chiudi.

25. Usate lo strumento ingrandisci (zoom in) in alcune aree per verificare i risultati. Noterete che il nuovo layer contiene solo le geometrie delle città che intersecano il buffer dei fiumi.

26. Dobbiamo sempre verificare che i nostri risultati siano validi. Un modo per verificare il risultato è esportare questo layer come un file KML e caricarlo in Google Earth. Così potrete vedere se le aree che abbiamo identificato sono realmente dentro i dieci chilometri dal fiume. Tasto destro del mouse sul layer e fate click su Salva con nome.

27. ella finestra Salva i vettori come...,scegliete WGS84 come SR. Questo perché il formato KML vuole che le coordinate siano in questo specifico sistema di riferimento. Chiamate il vostro file KML cities_near_river.kml.

28. Aprite Google Earth e verificate che le città rappresentate con questi buffer siano effettivamente vicine ai fiumi.