Delo s podatki o reliefu (QGIS3)

Podatki o reliefu so pomembni pri velikem številu prostorskih analiz, pogosto pa jih uporabljamo tudi pri kartiranju. QGIS3 ima vgrajene številna zmogljiva orodja za obdelavo reliefnih podatkov. V tej vadnici bomo ustvarili različne izdelke iz podatkov o višinah, na primer plastnice, senčen relief in podobno.

Pregled naloge

Cilj naloge je ustvariti karto plastnic in senčen relief okolice Everesta.

Česa se bomo še naučili

  • Iskanja in prenašanja brezplačnih podatkov o površju.

  • Izvažanja vektorskih podatkovnih slojev v zapis KML in prikaza s programom Google Zemlja (Google Earth).

  • Glajenja plastnic

Pridobivanje podatkov

Delali bomo s podatkovno zbirko GMTED2010, ki jo je pripravila ameriška geološka služba USGS. GMTED (Global Multi-resolution Terrain Elevation Data) je globalni nabor podatkov o terenu, ki je novejša različica nabora podatkov GTOPO30.

V naslednjih korakih je opisano iskanje in pobiranje podatkov s strani USGS EarthExplorer.

  1. Pojdite na USGS Earthexplorer . V zavihku Search Criteria kliknite World Features. V polje Feature Name vpišite Everest, v polje Country vpišite NEPAL in kliknite na Show. Prikazala se bo tabela z informacijami o lokaciji. V razdelku Placename izberite Everest.

../../_images/1168.png

  1. Zdaj se bo delovna površina premaknila na lokacijo Mt. Everest. Kliknite na Data Sets.

../../_images/2138.png

  1. Razširite skupino Digital Elevation in preverite GMTED2010. Kliknite na Results.

../../_images/377.png

  1. Kliknite gumb Download Options.

Opozorilo

Na tej točki se boste morali prijaviti na spletno mesto. Če računa še nimate, lahko ustvarite brezplačen račun.

../../_images/447.png

  1. Izberite možnost 30 ARC SEC in kliknite Download.

../../_images/541.png

Zdaj boste dobili datoteko z imenom GMTED2010N10E060_300.zip. Podatki o nadmorski višini so na voljo v različnih rastrskih formatih, kot so ASC, BIL, GeoTiff itd. QGIS podpira veliko različnih rastrskih formatov prek knjižnice GDAL. Podatki GMTED so na voljo v obliki datotek GeoTiff, ki so v tem arhivu zip.

Zaradi priročnosti lahko kopijo podatkov prenesete neposredno s spodnje strani.

GMTED2010N10E060_300.zip

Vir podatkov: [GMTED2010]

Postopek

  1. Odprite Layer ‣ Add Layer ‣ Add Raster Layer.

../../_images/640.png

  1. Kliknite na … pod Source, poiščite in izberite datoteko z imenom 10n060e_20101117_gmted_mea300.tif.

../../_images/740.png

  1. V QGIS delovni površini bodo prikazani podatki o terenu. Vsak piksel v rastru terena predstavlja povprečno nadmorsko višino v metrih na tem mestu. Temne pike predstavljajo območja z nizko nadmorsko višino, svetlejše pike pa območja z visoko nadmorsko višino.

../../_images/839.png

  1. Let’s find our area of interest. From Wikipedia, we find that the coordinates for our area of interest - Mt. Everest - is located at the coordinates 27.9881° N, 86.9253° E. Note that QGIS uses the coordinates in (latitude, longitude) format when the CRS is set to EPSG:4326. Therefore, paste the value 27.9881, 86.9253 these at the bottom of the QGIS window where it says Coordinate and press guilabel:Enter. The viewport will be centered at this coordinate. To zoom in, Enter 1:1000000 in the Scale field and press Enter. You will see the viewport zoom to the area around the Himalayas.

../../_images/940.png

  1. Zdaj bomo raster obrezali na to območje zanimanja. Poiščite izrez v Processing Toolbox. Izberite možnost Clip Raster by extent pod algoritmi GDAL.

../../_images/1047.png

  1. V oknu Clip Raster by Extent` izberite 10n060e_20101117_gmted_mea300 kot Input Layer`, kliknite ... v Clipping extent` in izberite Use Map canvas extent, kliknite ... v Clipped (extent)` in vnesite ime mt_everest.tif. Kliknite Run.

../../_images/1169.png

  1. Na delovni površini se pojavi nova plast mt_everest. Hill poiščite v Processing Toolbox. Izberite algoritem Hillshade pod algoritmi GDAL.

../../_images/1250.png

  1. V oknu Hillshade, izberite mt_everest kot Elevation Layer, v Azimuth (horizontal angle) vnesite 315.000, v Vertical angle vnesite 45.000. Kliknite ...` v Hillshade` in vnesite ime mt_everest_hillshade.tif. Kliknite Run.

../../_images/1348.png

  1. Na delovni površini se pojavi nova plast mt_everest_hillshade.

../../_images/1444.png

  1. Obris poiščite v Processing Toolbox. Izberite algoritem Contour pod GDAL algoritmi.

../../_images/1538.png

  1. In the Contour window, select mt_everest as Input Layer, enter 250 in Interval between contour lines. Click ... in Contours and enter the name as mt_everest_contour.gpkg. Click Run.

../../_images/1636.png

  1. Na delovni površini se pojavi nova plast mt_everest_contour. Z desno tipko miške kliknite na sloj in kliknite Open Attribute Table.

../../_images/1735.png

  1. Videli boste, da ima vsaka vrstica atribut z imenom ELEV. To je višina v metrih, ki jo predstavlja vsaka črta. Nekajkrat kliknite na glavo stolpca, da razvrstite vrednosti v padajočem vrstnem redu. Tu boste našli vrstico, ki predstavlja najvišjo nadmorsko višino v naših podatkih, tj. goro Everest.

../../_images/1830.png

  1. Izberite zgornjo vrstico in kliknite gumb Zoom to selection.

../../_images/1924.png

  1. Preklopite na glavno okno QGIS. Izbrana obrisna črta bo označena z rumeno barvo. To je območje najvišje nadmorske višine v našem podatkovnem nizu.

../../_images/2022.png

  1. Gladko iskanje v Processing Toolbox. Izberite Smooth pod Vector geometry.

../../_images/2139.png

  1. V oknu Smooth izberite mt_everest_contour kot Input Layer, v Iterations vnesite 5. Kliknite Run.

Opozorilo

Algoritem za glajenje deluje tako, da vzdolž črt dodaja dodatne vrhove. Z večanjem števila iteracij se število vrhov v obrisnih črtah močno poveča. Zato bodite previdni pri uporabi večjega števila iteracij. Velikost izhodne datoteke lahko zmanjšate tako, da jo izvozite kot datoteko shapefile in poenostavite rezultate z uporabo programa Mapshaper.

../../_images/2223.png

  1. Na delovni površini se bo pojavila nova plast Smoothed`. Ta plast bo imela bolj gladke robove v primerjavi s plastjo mt_everest_contour.

../../_images/2320.png

  1. Svoj sloj obrisov lahko vizualizirate in preverite svojo analizo tudi tako, da sloj obrisov izvozite kot KML in si ga ogledate v programu Google Earth. Z desnim gumbom miške kliknite na zglajeni sloj in izberite Export ‣ Save Feature As….

../../_images/2420.png

  1. Izberite Keyhole Markup Language [KML] kot Format. Kliknite ...` v polju File name` in vnesite ime contour_smoothed.kml. Kliknite OK.

../../_images/2519.png

  1. Poiščite izhodno datoteko na disku in jo dvakrat kliknite, da odprete program Google Earth Pro.

../../_images/2616.png

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