Analiza najbližjih sosedov (QGIS3)

GIS je zelo uporaben pri analizi prostorskih odnosov med značilnostmi. Ena takih analiz je ugotavljanje, katere značilnosti so najbližje določeni značilnosti. V QGIS lahko to analizo opravite na več načinov. Prostorsko združitev lahko izvedete z orodjem Join Attributes by Nearest ali pridobite razdalje do vseh elementov iz drugega sloja z orodjem Distance Matrix iz zbirke orodij za obdelavo. V tem vodniku bomo raziskali orodje z imenom Distance to nearest hub iz orodjarne Processing Toolbox, s katerim lahko ne le poiščete razdaljo do najbližjega elementa, temveč se mu za vizualizacijo rezultatov pridružite s črto.

Pregled naloge

Glede na lokacije vseh znanih potresov med letoma 1900 in 2000 poiščite najbližje naseljeno mesto za vsako lokacijo, kjer se je zgodil potres.

Pridobivanje podatkov

  1. V tej vadnici bomo prenesli nabor podatkov o potresih med letoma 1900 in 2000 iz Nacionalnega centra za geofizikalne podatke agencije NOAA, ki je pripravil odličen nabor podatkov o vseh pomembnih potresih od leta 2150 pred našim štetjem. obiščite portal NOAA NCEI <https://www.ngdc.noaa.gov/hazel/view/hazards/earthquake/search>`_ in vnesite Min kot 1900 in Max kot 2000. Tako boste dobili vse potrese, ki so se zgodili in jih je NOAA zabeležil med tema letoma. Za druge specifične rezultate lahko filtrirate z različnimi parametri. Kliknite :guilabel:Search.

../../_images/data011.png

  1. Posledično smo dobili 2585 potresov. Kliknite ikono Download TSV.

../../_images/data021.png

Natural Earth ima lepo zbirko podatkov Populated Places. Prenesite simple (less columns) dataset

Zaradi priročnosti lahko kopijo obeh zbirk podatkov neposredno prenesete s spodnjih povezav:

earthquakes_2021_11_25_14_31_59_+0530.tsv

ne_10m_populated_places_simple.zip

Viri podatkov: [NCEI] [NATURALEARTH]

Postopek

  1. Poiščite preneseno datoteko ne_10m_populated_places_simple.zip na plošči Browser in jo razširite. Datoteko ne_10m_populated_places_simple.shp povlecite na delovno površino.

../../_images/017.png

  1. Na plošči Layers boste videli naloženo nov sloj ne_10m_populated_places_simple. Ta sloj vsebuje točke, ki predstavljajo naseljena mesta. Zdaj bomo naložili sloj potresi. Ta sloj je na voljo kot Tab Serepated Values (TSV) besedilna datoteka. Če želite naložiti to datoteko, kliknite gumb Open Data Source Manager na plošči Data Source Toolbar. Uporabite lahko tudi Ctrl + L bližnjico na tipkovnici.

../../_images/027.png

  1. V pogovornem oknu Data Source Manager izberite Delimited Text.

../../_images/037.png

  1. Kliknite gumb … poleg File name in poiščite preneseno datoteko earthquakes-2021-11-25_13-39-30_+0530.tsv. Odvisno od operacijskega sistema morda datoteke ne boste videli v prenesenem imeniku. V tem primeru v pogovornem oknu All files (*; .) preklopite na Choose a Delimited Text File to Open. Po odprtju izberite Custom delimiters v razdelku File format in označite Tab. V razdelku Geometry definition (Opredelitev geometrije) izberite Point coordinates (Koordinate točke). Privzeto bosta vrednosti Polje X in Polje Y samodejno napolnjeni z ustreznimi polji v vnosu. V našem primeru sta to Longitude in Latitude. Polje Geometry CRS lahko nastavite na privzeto vrednost EPSG:4326 - WGS 84 CRS. Če vaša datoteka vsebuje koordinate v drugem CRS, lahko tukaj izberete ustrezen CRS. Kliknite Add in nato Close.

../../_images/047.png

  1. Povečajte in raziščite oba podatkovna niza. Vsaka rdeča točka predstavlja lokacijo potresa, vsaka zelena točka pa lokacijo naseljenega kraja. Naš cilj je, da za vsako točko v sloju potresov poiščemo najbližjo točko iz sloja naseljenih krajev. Preglejmo tabelo atributov sloja potresov. Izberite sloj in kliknite na ikono Open Attribute Table (Odpri atributno tabelo) v Toolbar (Orodna vrstica).

../../_images/057.png

  1. Obstaja 2586 funkcij, vendar podatki vsebujejo le nekaj vnosov brez podatkov o zemljepisni širini ali dolžini. To moramo odstraniti, preden nadaljujemo z delom. Zaprite tabelo atributov.

../../_images/067.png

  1. Pojdite na orodje Processing ‣ Toolbox ‣ Vector geometry ‣ Remove null geometries. Dvakrat kliknite, da ga odprete.

../../_images/077.png

  1. V pogovornem oknu Remove Null Geometries (Odstrani ničelne geometrije) izberite earthquakes-2021-11-25_13-39-30_+0530 kot Input layer (Vhodna plast) in potrdite polje Also remove empty geometries (Odstrani tudi prazne geometrije). Kliknite Run. Ko se obdelava konča, kliknite Close.

../../_images/086.png

  1. Na ploščo Layers bo dodan nov sloj Non null geometries. Za analizo bomo uporabili ta sloj namesto prvotnega sloja. Sliko earthquakes-2021-11-25_13-39-30_+0530 na plošči Layers` odkljukajte, da jo skrijete. Izberite plast Non null geometries in kliknite gumb Open Attribute Table (Odpri tabelo atributov) v orodni vrstici Attributes Toolbar (Orodna vrstica atributov).

../../_images/096.png

  1. Videlo se bo manjše število vseh elementov, saj so bile odstranjene vse vrstice s praznimi vrednostmi zemljepisne širine in dolžine. Zaprite tabelo atributov.

../../_images/1024.png

  1. Zdaj je čas za analizo najbližjega soseda. Poiščite in poiščite orodje Processing ‣ Toolbox ‣ Vector analysis ‣ Distance to nearest hub (line to hub). Z dvojnim klikom ga zaženite.

../../_images/1127.png

Opomba

Kot izhod lahko dodamo tudi plast točk, za kar uporabite orodje Distance to nearest hub (points).

  1. V pogovornem oknu Distance to Nearest Hub (Line to Hub) izberite Non null geometries kot Source points layer. Izberite ne_10m_populated_places_simple kot sloj Destination hubs layer. Izberite name kot atribut Hub name layer. Orodje bo izračunalo tudi neposredno razdaljo med naseljenim krajem in najbližjim potresom. Kot Measurement unit (merilna enota) nastavite Kilometers. Kliknite na ... v polju Hub Distance` in kliknite Save to File…`, da shranite datoteko kot earthquakes_with_nearest_city.gpkg . Kliknite Run. Ko se obdelava konča, kliknite Close.

../../_images/1225.png

  1. V glavnem oknu QGIS boste videli, da je na plošči Layers naložena nova vrstica z imenom earthquakes_with_nearest_city. Ta sloj ima linijske elemente, ki povezujejo vsako potresno točko z najbližjim naseljenim mestom. Izberite sloj earthquakes_with_nearest_city in kliknite ikono Open Attribute Tabel v Toolbar.

../../_images/1323.png

  1. Pomaknite se desno do zadnjih stolpcev in videli boste dva nova atributa, imenovana HubName in HubDist, ki sta dodana prvotnim značilnostim potresa. To je ime razdalje do najbližjega soseda iz sloja naseljenih krajev.

../../_images/1421.png

If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)