Rastersko mozaičko kreiranje i isecanje (QGIS3)

Ovaj tutorijal istražuje osnovne tehnike rada sa rasterima u QGIS-u, kao što su mozaikiranje i podešavanje.

Pregled zadatka

Preuzećemo podatke o nadmorskoj visini za Šri Lanku u obliku SRTM pločica, spojiti ih i rezultujući mozaik iseći na granicu države.

Druge veštine koje ćete naučiti

  • Korišćenje Hillshade renderera za vizuelizaciju podataka o nadmorskoj visini.

Dobijte podatke

Centar za distribuiranu aktivnu arhivu kopnenih procesa (LP DAAC) pruža skup podataka „NASA Shuttle Radar Topography Mission (SRTM) Global 1 lučna sekunda <https://lpdaac.usgs.gov/products/srtmgl1v003/>“ kao pločice elevacije.

Jednostavan interfejs za preuzimanje pločica za dato područje je 30-Meter SRTM Tile Downloader autora Dereka Votkinsa. Preuzmite pojedinačne SRTM pločice koje pokrivaju Šri Lanku. Imajte na umu da će vam biti potreban besplatan Earth Data nalog da biste preuzeli podatke.

../../_images/data120.png

Takođe će nam biti potreban shapefile Admin 0 - Countries iz programa Natural Earth.

Radi lakšeg snalaženja, možete direktno preuzeti kopiju skupova podataka sa donjih linkova:

N05E080.SRTMGL1.hgt.zip

N06E079.SRTMGL1.hgt.zip

N06E080.SRTMGL1.hgt.zip

N06E081.SRTMGL1.hgt.zip

N07E079.SRTMGL1.hgt.zip

N07E080.SRTMGL1.hgt.zip

N07E081.SRTMGL1.hgt.zip

N08E079.SRTMGL1.hgt.zip

N08E080.SRTMGL1.hgt.zip

N08E081.SRTMGL1.hgt.zip

N09E080.SRTMGL1.hgt.zip

ne_10m_admin_0_countries.zip

Izvor podataka [SRTM] , [NATURALEARTH]

Procedura

  1. Otvorite QGIS i pronađite preuzete datoteke u panelu Browser. Proširite pojedinačne zip datoteke da biste prikazali .hgt datoteke. Držite pritisnut taster Ctrl i izaberite sve pojedinačne datoteke. Kada ih izaberete, prevucite ih na platno.

../../_images/1148.png

  1. Videćete 11 pojedinačnih slojeva učitanih u panelu Slojevi i prikazanih na platnu. Spojićemo ove pojedinačne pločice u jedan mozaik. Idite na Obrada ‣ Alatke.

../../_images/2126.png

  1. Pretražite i pronađite alatku GDAL ‣ Raster razno ‣ Spajanje. Kliknite dvaput da biste je pokrenuli.

../../_images/364.png

  1. U dijalogu Spajanje kliknite na dugme … pored Unesi slojeve. Kliknite na Izaberi sve da biste izabrali sve pojedinačne slojeve.

../../_images/435.png

  1. Kao što je pomenuto u detaljima sloja skupa podataka, tip ulaznih podataka je 16-bitni označeni ceo broj. Da bismo održali integritet podataka, trebalo bi da zadržimo isti tip podataka za spojeni sloj. Izaberite Int16 kao Izlazni tip podataka. Takođe, podrazumevani format izlaznih podataka je GeoTiff. GeoTiff datoteke mogu postati veoma velike ako se ne kompresuju. Izaberite Visoka kompresija kao Profil. Kliknite na Pokreni.

../../_images/531.png

  1. Kada se obrada završi, novi sloj „Spojen“ biće dodat na panel Slojevi. Ukoliko sloj nije na vrhu steka, izaberite ga i prevucite na vrh panela Slojevi.

../../_images/630.png

  1. Videćete da sloj „Spojeno“ sadrži spojene podatke o visini iz pojedinačnih ulaznih pločica. Podrazumevana vizuelizacija prikazuje samo vrednosti piksela u opsegu od 0-255. Ali naši podaci sadrže piksele sa vrednostima od -14 do 2371, što rezultira renderovanjem sa niskim kontrastom. Hajde da to promenimo u bolju vizuelizaciju. Kliknite na dugme Otvori panel za stilizovanje sloja u panelu Slojevi.

../../_images/730.png

  1. U panelu Stilizovanje sloja, kliknite na padajući meni Tip renderovanja i izaberite renderer Hilsejd. Ova opcija renderovanja je posebno pogodna za podatke o nadmorskoj visini.

../../_images/829.png

  1. Još jedna uobičajena operacija pri radu sa rasterima je isecanje rastera na područje koje vas zanima. U ovom tutorijalu, iseckaćemo spojeni sloj na granicu države Šri Lanka. Pronađite preuzetu datoteku ne_10m_admin_0_countries.zip i proširite je. Prevucite datoteku ne_10m_admin_0_countries.shp na platno.

../../_images/930.png

  1. Izaberite novododate slojeve „ne_10m_admin_0_countries“ u panelu „Slojevi“. Kliknite na dugme „Izaberite obeležja po oblasti ili jednim klikom“ na traci sa alatkama „Atributi“. Kada ih izaberete, kliknite na poligon za Šri Lanku da biste ga izabrali.

../../_images/1037.png

  1. Zadržite selekciju kakva jeste i otvorite Processing ‣ Toolbox. Potražite i pronađite alatku GDAL ‣ Raster extraction ‣ Iseci raster pomoću maski sloja. Dvaput kliknite da biste je pokrenuli.

../../_images/1149.png

  1. U dijalogu Iseci raster pomoću maskiranog sloja, podesite Spojeno kao Ulazni sloj. Izaberite ne_10m_admin_0_countries kao Maskiraj sloj i označite polje za potvrdu Samo odabrane karakteristike. Unesite 0.0000 kao Dodeli određenu vrednost nodata izlaznim opsezima. Kao i ranije, izaberite Visoka kompresija kao Profil. Kliknite na Pokreni.

../../_images/1239.png

  1. Novi sloj „Odsečen (maska)“ biće dodat na panel Slojevi. U ovom trenutku može biti teško videti rezultat jer imamo previše preklapajućih slojeva. Kliknite na dugme Upravljaj temama mape na panelu Slojevi i izaberite Sakrij sve slojeve.

../../_images/1337.png

  1. Uključite samo poslednji sloj „Odsečeno (maska)“ i stilizujte ga pomoću renderera „Hilshade“ kao što je prethodno urađeno.

../../_images/1434.png

  1. Spojeni i podskupljeni izlazni sloj elevacije za Šri Lanku je spreman.

../../_images/1530.png

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