Analiza najbližeg suseda (QGIS3)¶
GIS je veoma koristan u analizi prostornog odnosa između objekata. Jedna takva analiza je pronalaženje koji su objekti najbliži datom objektu. Postoji više načina da se ova analiza izvrši u QGIS-u. Možete izvršiti prostorno spajanje koristeći alatku Spoji atribute po najbližem ili dobiti udaljenosti do svih objekata iz drugog sloja koristeći alatku Matrica udaljenosti iz okvira alata za obradu. U ovom tutorijalu ćemo istražiti alatku pod nazivom Udaljenost do najbližeg čvorišta iz okvira alata za obradu koji ne samo da može da pronađe udaljenost do najbližeg objekta već ga i da spoji linijom sa njim radi vizuelizacije rezultata.
Pregled zadataka¶
Na osnovu lokacija svih poznatih zemljotresa između 1900. i 2000. godine, pronađite najbliže naseljeno mesto za svaku lokaciju gde se zemljotres dogodio.
Dobijte podatke¶
Za ovaj tutorijal preuzećemo skup podataka o zemljotresima između 1900-2000. godine iz Nacionalnog centra za geofizičke podatke NOAA, koji proizvodi odličan skup podataka o svim značajnim zemljotresima od 2150. godine pre nove ere. Posetite NOAA NCEI portal i unesite Min kao
1900
i Max kao2000
. Ovo će vratiti sve zemljotrese koji su se dogodili i koje je NOAA zabeležila između tih godina. Za druge specifične rezultate, možete filtrirati pomoću različitih parametara. Kliknite na Search.
Kao rezultat toga, dobili smo 2585 incidenata zemljotresa. Kliknite na ikonu Preuzmi TSV.
Prirodna Zemlja ima lep skup podataka „Naseljena mesta <http://www.naturalearthdata.com/downloads/10m-cultural-vectors/10m-populated-places/>“. Preuzmite jednostavan skup podataka (sa manje kolona) <http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/10m/cultural/ne_10m_populated_places_simple.zip>`_
Radi lakšeg snalaženja, možete direktno preuzeti kopiju oba skupa podataka sa donjih linkova:
earthquakes_2021_11_25_14_31_59_+0530.tsv
ne_10m_populated_places_simple.zip
Izvori podataka: [NCEI] [NATURALEARTH]
Procedura¶
Pronađite preuzetu datoteku
ne_10m_populated_places_simple.zip
u panelu Pregledač i proširite je. Prevucite datotekune_10m_populated_places_simple.shp
na platno.
Videćete novi sloj
ne_10m_populated_places_simple
učitan u Layers panelu. Ovaj sloj sadrži tačke koje predstavljaju naseljena mesta. Sada ćemo učitati sloj zemljotresa. Ovaj sloj dolazi kao tekstualna datoteka Tab Serepated Values (TSV). Da biste učitali ovu datoteku, kliknite na dugme Otvori menadžer izvora podataka na Traci sa alatkama za izvore podataka. Takođe možete koristiti prečicu na tastaturi Ctrl + L.
U dijaloškom okviru Upravljač izvorima podataka izaberite Razgraničeni tekst.
Kliknite na dugme … pored Ime datoteke i pregledajte preuzetu datoteku
earthquakes-2021-11-25_13-39-30_+0530.tsv
. U zavisnosti od operativnog sistema, možda nećete videti datoteku u preuzetom direktorijumu. Ako je to slučaj, prebacite se na Sve datoteke (*; .) u dijalogu Izaberite datoteku sa razgraničenim tekstom za otvaranje. Kada se otvori, izaberite Prilagođeni razgraničnici u odeljku Format datoteke i označiteTab
. U odeljku Definicija geometrije izaberite Koordinate tačaka. Podrazumevano, vrednosti X polje i Y polje će biti automatski popunjene odgovarajućim poljima u unosu. U našem slučaju, to sugeografska dužina
igeografska širina
. Možete ostaviti Geometry CRS na podrazumevanomEPSG:4326 - WGS 84
CRS-u. Ako vaša datoteka sadrži koordinate u drugom CRS-u, možete ovde izabrati odgovarajući CRS. Kliknite na Dodaj, a zatim na Zatvori.
Zumirajte i istražite oba skupa podataka. Svaka crvena tačka predstavlja lokaciju pojave zemljotresa, a svaka zelena tačka predstavlja lokaciju naseljenog mesta. Naš cilj je da pronađemo najbližu tačku iz sloja naseljenih mesta za svaku od tačaka u sloju zemljotresa. Hajde da pregledamo tabelu atributa sloja zemljotresa. Izaberite sloj i kliknite na ikonu Otvori tabelu atributa u Traci sa alatkama.
Postoje objekti „2586“, ali podaci sadrže nekoliko unosa bez informacija o geografskoj širini ili dužini. Moramo to ukloniti pre nego što nastavimo dalje. Zatvorite tabelu atributa.
Idite na
. Dvaput kliknite da biste ga otvorili.
U dijaloškom okviru Ukloni nulte geometrije, izaberite
earthquakes-2021-11-25_13-39-30_+0530
kao Ulazni sloj i označite polje Takođe ukloni prazne geometrije. Kliknite na Pokreni. Kada se obrada završi, kliknite na Zatvori.
Novi sloj „Non null geometrije“ biće dodat na panel Slojevi. Za analizu ćemo koristiti ovaj sloj umesto originalnog sloja. Uklonite oznaku sa sloja
earthquakes-2021-11-25_13-39-30_+0530
na panelu Slojevi da biste ga sakrili. Izaberite slojNon null geometrije
i kliknite na dugme Otvori tabelu atributa iz Trake sa alatkama za atribute.
Videćete manji broj za ukupan broj obeležja jer su svi redovi sa praznim vrednostima geografske širine i dužine uklonjeni. Zatvorite tabelu atributa.
Sada je vreme da izvršite analizu najbližeg suseda. Potražite i pronađite alatku
. Dvaput kliknite da biste je pokrenuli.
Белешка
Takođe možemo dodati sloj tačaka kao izlaz, za to koristimo alatku Udaljenost do najbližeg čvorišta (tačke).
U dijaloškom okviru Udaljenost do najbližeg čvorišta (linija do čvorišta), izaberite
Nenulte geometrije
kao Sloj izvornih tačaka. Izaberitene_10m_populated_places_simple
kao Sloj odredišnih čvorišta. Izaberitenaziv
kao Atribut imena sloja čvorišta. Alat će takođe izračunati pravolinijsko rastojanje između naseljenog mesta i najbližeg zemljotresa. Postavitekilometar
kao jedinicu mere. Kliknite na...
u Udaljenost do čvorišta i kliknite na Sačuvaj u datoteku… da biste sačuvali datoteku kaoearthquakes_with_nearest_city.gpkg
. Kliknite na Pokreni. Kada se obrada završi, kliknite na Zatvori.
Nazad u glavnom QGIS prozoru, videćete novi linijski sloj pod nazivom „earthquakes_with_nearest_city”” učitan u Layers panelu. Ovaj sloj ima linijske karakteristike koje povezuju svaku tačku zemljotresa sa najbližim naseljenim mestom. Izaberite sloj „earthquakes_with_nearest_city” i kliknite na ikonu Otvori tabelu atributa u Toolbar.
Pomerite se udesno do poslednjih kolona i videćete 2 nova atributa pod nazivom HubName i HubDist dodata originalnim karakteristikama zemljotresa. Ovo je naziv udaljenosti do najbližeg suseda iz sloja naseljenih mesta.
If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)