Mozaičenje in obrezovanje rastrov (QGIS3)

V tem učbeniku so predstavljene osnovne tehnike dela z rastri v QGIS, kot sta mozaičenje in izrezovanje.

Pregled naloge

Za Šrilanko bomo prenesli višinske podatke v obliki ploščic SRTM, jih združili in nastali mozaik pripeli na državno mejo.

Česa se bomo še naučili

  • Uporabe upodobitvenega programa Hillshade za vizualizacijo podatkov o nadmorski višini.

Pridobivanje podatkov

Land Processes Distributed Active Archive Center (LP DAAC) zagotavlja NASA Shuttle Radar Topography Mission (SRTM) Global 1 arc second podatkovni niz v obliki višinskih delov (tile).

Enostaven vmesnik za prenos delov za določeno območje je 30-Meter SRTM Tile Downloader Dereka Watkinsa. Prenesite posamezne dele SRTM, ki pokrivajo Šrilanko. Upoštevajte, da boste za prenos podatkov potrebovali brezplačen račun Earth Data.

../../_images/data115.png

Potrebovali bomo tudi datoteko Admin 0 - Countries shapefile iz programa Natural Earth.

Kopijo podatkovnih nizov lahko neposredno prenesete s spodnjih povezav:

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

Vir podatkov [SRTM] , [NATURALEARTH]

Postopek

  1. Odprite QGIS in poiščite prenesene datoteke na plošči Browser. Razširite posamezne datoteke zip, da se prikažejo datoteke .hgt. Držite tipko Ctrl in izberite vse posamezne datoteke. Ko so izbrane, jih povlecite na platno.

../../_images/1108.png
  1. Na plošči Layers se naloži 11 slojev, ki so prikazani na platnu. Te posamezne plasti bomo združili v en sam mozaik. Pojdite na Processing ‣ Toolbox.

../../_images/2103.png
  1. Poiščite in poiščite orodje GDAL ‣ Raster miscellaneous ‣ Merge. Z dvoklikom ga zaženite.

../../_images/344.png
  1. V pogovornem oknu Merge kliknite gumb poleg Input layers. Kliknite Select All (izberi vse), da izberete vse posamezne plasti.

../../_images/423.png
  1. Kot je navedeno v podrobnostih o podatkovni plasti <https://lpdaac.usgs.gov/products/srtmgl1v003/>`_, je tip vhodnih podatkov 16-bit signed integer. Da bi ohranili celovitost podatkov, moramo za združeno plast ohraniti isti podatkovni tip. Izberite Int16 kot Output data type. Tudi privzeti format izhodnih podatkov je GeoTiff. Datoteke GeoTiff so lahko zelo velike, če niso stisnjene. Kot Profile izberite High Compression. Kliknite Run.

../../_images/522.png
  1. Ko se obdelava konča, se na ploščo Layers doda nova plast OUTPUT. Če plast ni na vrhu sklada, jo izberite in povlecite na vrh plošče Layers.

../../_images/621.png
  1. Videli boste, da sloj OUTPUT vsebuje združene višinske podatke iz posameznih vhodnih ploščic. Privzeta vizualizacija prikazuje samo vrednosti pikslov v razponu od 0-255. Naši podatki pa vsebujejo piksle z vrednostmi od -14 do 2371, zato je upodobitev nizkokontrastna. Spremenimo jo na boljšo vizualizacijo. Kliknite gumb Open the layer Styling panel (odpri ploščo za oblikovanje slojev) na plošči Layers (sloje).

../../_images/721.png
  1. Na plošči Layer Styling kliknite spustno okno Render type in izberite renderer Hillshade. Ta možnost upodabljanja je še posebej primerna za podatke o nadmorski višini.

../../_images/820.png
  1. Druga pogosta operacija pri delu z rastri je pripenjanje rastra na območje, ki vas zanima. V tem učbeniku bomo združeni sloj pripeli na državno mejo Šrilanke. Poiščite preneseno datoteko ne_10m_admin_0_countries.zip in jo razširite. Datoteko ne_10m_admin_0_countries.shp povlecite na platno.

../../_images/921.png
  1. Izberite novo dodane plasti ne_10m_admin_0_countries na plošči Layers. Kliknite gumb Select Features by area or single click (Izberi elemente po območju ali z enim klikom) v orodni vrstici Attributes Toolbar (Orodna vrstica Atributi). Ko je izbran, kliknite poligon za Šrilanko, da ga izberete.

../../_images/1028.png
  1. Ohranite izbor, kot je, in odprite Processing ‣ Toolbox. Poiščite orodje GDAL ‣ Raster extraction ‣ Clip raster by mask layer. Z dvojnim klikom ga zaženite.

../../_images/1132.png
  1. V pogovornem oknu Clip Raster by Mask Layer nastavite OUTPUT kot Input Layer. Izberite ne_10m_admin_0_countries kot Mask layer` in potrdite potrditveno polje Selected features only`. Vnesite 0,0000 kot Assign a specified nodata value to output bands. Kot prej izberite High compression kot Profile. Kliknite Run.

../../_images/1230.png
  1. Na ploščo Layers bo dodana nova plast OUTPUT. Na tej točki bo morda težko videti izhod, ker imamo vidnih preveč prekrivajočih se slojev. Kliknite gumb Manage Map Themes na plošči Layers in izberite Hide All Layers.

../../_images/1328.png
  1. Vključite samo zadnjo plast OUTPUT in jo oblikujte z prikazom Hilshade kot prej.

../../_images/1425.png
  1. Združeni in obrezani izhodni sloj višin za Šrilanko je pripravljen.

../../_images/1524.png
comments powered by Disqus