Svarbu

Vertimas yra bendruomenės pastangos, prie kurių jūs galite prisijungti. Šis puslapis šiuo metu išverstas 78.33%.

24.1.17. Rastro paviršiaus analizė

24.1.17.1. Aspektas

Skaičiuoja įvestyje pateikto skaitmeninio paviršiaus modelio aspektą. Gautame aspekto rastro sluoksnyje yra reikšmės nuo 0 iki 360, kurios išreiškia nuolydžio kryptį, pradedant nuo šiaurės (0°) ir skaičiuojant prieš laikrodžio rodyklę.

Fig. 24.30 Aspekto reikšmės

Žemiau pateikiama iliustracija, kur aspekto sluoksnis perklasifikuotas naudojant spalvų rampą:

Fig. 24.31 Perklasifikuotas aspekto sluoksnis

Parametrai

Užrašas	Pavadinimas	Tipas	Aprašymas
Aukščio sluoksnis	INPUT	[rastras]	Skaitmeninio paviršiaus modelio rastro sluoksnis
Z koeficientas	Z_FACTOR	[skaičius: slankaus kablelio] Numatytas: 1.0	Vertikalus didinimas. Šis parametras naudingas, kai Z vienetai skiriasi nuo X ir Y vienetų, pavyzdžiui pėdos ir metrai. Jūs galite naudoti šį parametrą, kad prisiderintumėte prie šito. Numatytoji reikšmė yra 1 (be padidinimo).
Aspektas	OUTPUT	[rastras] Numatytas: [Įrašyti į laikiną failą]	Nurodykite išvesties aspekto rastro sluoksnį. Vienas iš: Įrašyti į laikiną failą Įrašyti į failą…

Išvestys

Užrašas	Pavadinimas	Tipas	Aprašymas
Aspektas	OUTPUT	[rastras]	Išvesties aspekto rastro sluoksnis

Pythono kodas

Algoritmo ID: qgis:aspect

import processing
processing.run("algorithm_id", {parameter_dictionary})

algoritmo id rodomas, kai užvedate pelę virš algoritmo apdorojimo įrankinėje. parametrų žodynas teikia parametrų pavadinimus ir reikšmes. Daugiau informacijos apie tai, kaip vykdyti apdorojimo algoritmus Pythono konsolėje rasite skyriuje Using processing algorithms from the console.

24.1.17.2. DTM filtras (pagal nuolydį)

Gali būti naudojamas filtruojant skaitmeninį aukščio modelį klasifikuojant jo celes į žemės ir objektų (ne žemės) celes.

Įrankis naudoja koncepcijas, apibrėžtas Vosselman (2000) ir remiasi prielaida, kad dideli aukščio skirtumai tarp gretimų celių negali būti dėl paviršiaus nuolydžio. Tikimybė, kad aukštesnės celės gali būti ne žemė didėja mažėjant atstumui tarp dviejų celių. Todėl filtras apibrėžia maksimalų aukščio skirtumą (dz_max) tarp celių kaip funkciją atstumo (d) tarp celių (dz_max( d ) = d). Celė klasifikuojama kaip paviršiaus, jei branduolio spindulio atstumu nėra celės, kurios aukščio skirtumas būtų didesnis nei nurodytas aukščio skirtumų maksimumas nurodytu atstumu tarp dviejų celių.

Apytikslio paviršiaus nuolydžio (s) parametras naudojamas keisti filtro funkciją, kad ji atitiktų bendrą nuolydį studijuojamame plote (dz_max( d ) = d * s). 5 % užtikrintumo intervalas (ci = 1.65 * sqrt( 2 * stddev )) gali būti naudojamas keisti filtro funkciją dar labiau, arba atlaisvinant (dz_max( d ) = d * s + ci), arba padidinant (dz_max( d ) = d * s - ci) filtro kriterijų.

Šaltinis: Vosselman, G. (2000): Slope based filtering of laser altimetry data. IAPRS, Vol. XXXIII, Part B3, Amsterdam, The Netherlands, 935-942

Taip pat žiūrėkite

Šis įrankis yra perkeltas iš SAGA DTM Filtro (paremto nuolydžiu)

Parametrai

Baziniai parametrai

Užrašas	Pavadinimas	Tipas	Aprašymas
Įvesties sluoksnis	INPUT	[rastras]	Skaitmeninio paviršiaus modelio rastro sluoksnis
Juostos numeris	BAND	[rastro juosta]	Tiriama DEM juosta
Branduolio spindulys (pikseliais)	RADIUS	[skaičius: sveikas skaičius] Numatytas:	Filtro branduolio spindulys (pikseliais). Turi būti pakankamai didelis, kad pasiektų paviršiaus celes greta ne žemės objektų.
Paviršiaus nuolydis (%, pikselio dydis/vertikalūs vienetai)	TERRAIN_SLOPE	[skaičius: slankaus kablelio] Numatytas: 30.0	Apytikslis paviršiaus nuolydis %. Paviršiaus nuolydį reikia patikslinti, kad būtų atsižvelgta į aukščio vienetų ir rastro pikselių matmenų santykį. Naudojama atlaisvinti filtro kriterijų statesniame paviršiuje.
Filtro modifikacija	FILTER_MODIFICATION	[sąrašas] Numatytas: 0	Pasirinkite, ar taikyti filtro branduolį be pokyčių, ar naudoti užtikrintumo intervalą aukščio kriterijaus atlaisvinimui ar padidinimui. 0 - Jokio 1 - Atlaisvinti filtrą 2 - Padidinti
Standartinis nuokrypis	STANDARD_DEVIATION	[skaičius: slankaus kablelio] Numatytas: 0.1	Standartinis nuokrypis naudojamas skaičiuojant 5% užtikrintumo intervalą, taikomą aukščio slenksčiui.
Išvesties sluoksnis (žemė) Pasirinktinis	OUTPUT_GROUND	[rastras] Numatytas: [Įrašyti į laikiną failą]	Nurodykite filtruotą DEM, kuriame yra tik celės, klasifikuotos kaip žemė. Vienas iš: Praleisti išvestį Įrašyti į laikiną failą Įrašyti į failą…
Išvesties sluoksnis (ne žemės objektai) Pasirinktinis	OUTPUT_NONGROUND	[rastras] Numatytas: [Praleisti išvestį]	Nurodykite ne žemės objektus, pašalintus filtro. Vienas iš: Praleisti išvestį Įrašyti į laikiną failą Įrašyti į failą…

Išmanesni parametrai

Užrašas	Pavadinimas	Tipas	Aprašymas
Kūrimo parinktys Pasirinktinis	CREATE_OPTIONS	[tekstas] Numatytas: ‚‘	Norint pridėti vieną ar daugiau kūrimo parinkčių, valdančių rastro kūrimą (spalvos, bloko dydis, failo suspaudimas…). Patogumui jūs galite remtis iš anksto sukurtais profiliais (žiūrėkite GDAL tvarkyklės parinkčių skiltį). Paketinis apdorojimas ir Modelių kūrimas: atskirkite kelias parinktis statmenu brūkšniu (|).

Išvestys

Užrašas	Pavadinimas	Tipas	Aprašymas
Išvesties sluoksnis (žemė)	OUTPUT_GROUND	[rastras]	Filtruotas DEM, kuriame yra tik celės, klasifikuotos kaip žemė.
Išvesties sluoksnis (ne žemės objektai)	OUTPUT_NONGROUND	[rastras]	Ne žemės objektai, kuriuos pašalino filtras.

Pythono kodas

Algoritmo ID: native:dtmslopebasedfilter

import processing
processing.run("algorithm_id", {parameter_dictionary})

algoritmo id rodomas, kai užvedate pelę virš algoritmo apdorojimo įrankinėje. parametrų žodynas teikia parametrų pavadinimus ir reikšmes. Daugiau informacijos apie tai, kaip vykdyti apdorojimo algoritmus Pythono konsolėje rasite skyriuje Using processing algorithms from the console.

24.1.17.3. Hillshade

Calculates the hillshade raster layer given an input Digital Terrain Model.

The shading of the layer is calculated according to the sun position: you have the options to change both the horizontal angle (azimuth) and the vertical angle (sun elevation) of the sun.

Fig. 24.32 Azimuth and vertical angle

The hillshade layer contains values from 0 (complete shadow) to 255 (complete sun). Hillshade is used usually to better understand the relief of the area.

Fig. 24.33 Hillshade layer with azimuth 300 and vertical angle 45

Particularly interesting is to give the hillshade layer a transparency value and overlap it with the elevation raster:

Fig. 24.34 Overlapping the hillshade with the elevation layer

Parametrai

Užrašas	Pavadinimas	Tipas	Aprašymas
Aukščio sluoksnis	INPUT	[rastras]	Skaitmeninio paviršiaus modelio rastro sluoksnis
Z koeficientas	Z_FACTOR	[skaičius: slankaus kablelio] Numatytas: 1.0	Vertical exaggeration. This parameter is useful when the Z units differ from the X and Y units, for example feet and meters. You can use this parameter to adjust for this. Increasing the value of this parameter will exaggerate the final result (making it look more „hilly“). The default is 1 (no exaggeration).
Azimuth (horizontal angle)	AZIMUTH	[skaičius: slankaus kablelio] Default: 300.0	Set the horizontal angle (in degrees) of the sun (clockwise direction). Range: 0 to 360. 0 is north.
Vertical angle	V_ANGLE	[skaičius: slankaus kablelio] Default: 40.0	Set the vertical angle (in degrees) of the sun, that is the height of the sun. Values can go from 0 (minimum elevation) to 90 (maximum elevation).
Hillshade	OUTPUT	[rastras] Default: Save to temporary file	Specify the output hillshade raster layer. One of: Įrašyti į laikiną failą Įrašyti į failą…

Išvestys

Užrašas	Pavadinimas	Tipas	Aprašymas
Hillshade	OUTPUT	[rastras]	The output hillshade raster layer

Pythono kodas

Algorithm ID: qgis:hillshade

import processing
processing.run("algorithm_id", {parameter_dictionary})

algoritmo id rodomas, kai užvedate pelę virš algoritmo apdorojimo įrankinėje. parametrų žodynas teikia parametrų pavadinimus ir reikšmes. Daugiau informacijos apie tai, kaip vykdyti apdorojimo algoritmus Pythono konsolėje rasite skyriuje Using processing algorithms from the console.

24.1.17.4. Hypsometric curves

Calculates hypsometric curves for an input Digital Elevation Model. Curves are produced as CSV files in an output folder specified by the user.

A hypsometric curve is a cumulative histogram of elevation values in a geographical area.

You can use hypsometric curves to detect differences in the landscape due to the geomorphology of the territory.

Parametrai

Užrašas	Pavadinimas	Tipas	Aprašymas
DEM to analyze	INPUT_DEM	[rastras]	Digital Terrain Model raster layer to use for calculating altitudes
Boundary layer	BOUNDARY_LAYER	[vektorius: poligonas]	Polygon vector layer with boundaries of areas used to calculate hypsometric curves
Step	STEP	[skaičius: slankaus kablelio] Numatytas: 100.0	Vertical distance between curves
Use % of area instead of absolute value	USE_PERCENTAGE	[boolean] Numatytas: ne	Write area percentage to “Area” field of the CSV file instead of the absolute area
Hypsometric curves	OUTPUT_DIRECTORY	[aplankas]	Specify the output folder for the hypsometric curves. One of: Įrašyti į laikiną aplanką Įrašyti į aplanką

Išvestys

Užrašas	Pavadinimas	Tipas	Aprašymas
Hypsometric curves	OUTPUT_DIRECTORY	[aplankas]	Directory containing the files with the hypsometric curves. For each feature from the input vector layer, a CSV file with area and altitude values will be created. The file names start with histogram_, followed by layer name and feature ID.

Pythono kodas

Algorithm ID: qgis:hypsometriccurves

import processing
processing.run("algorithm_id", {parameter_dictionary})

algoritmo id rodomas, kai užvedate pelę virš algoritmo apdorojimo įrankinėje. parametrų žodynas teikia parametrų pavadinimus ir reikšmes. Daugiau informacijos apie tai, kaip vykdyti apdorojimo algoritmus Pythono konsolėje rasite skyriuje Using processing algorithms from the console.

24.1.17.5. Relief

Creates a shaded relief layer from digital elevation data. You can specify the relief color manually, or you can let the algorithm choose automatically all the relief classes.

Fig. 24.35 Relief layer

Parametrai

Užrašas	Pavadinimas	Tipas	Aprašymas
Aukščio sluoksnis	INPUT	[rastras]	Skaitmeninio paviršiaus modelio rastro sluoksnis
Z koeficientas	Z_FACTOR	[skaičius: slankaus kablelio] Numatytas: 1.0	Vertical exaggeration. This parameter is useful when the Z units differ from the X and Y units, for example feet and meters. You can use this parameter to adjust for this. Increasing the value of this parameter will exaggerate the final result (making it look more „hilly“). The default is 1 (no exaggeration).
Generate relief classes automatically	AUTO_COLORS	[boolean] Numatytas: ne	If you check this option the algorithm will create all the relief color classes automatically
Relief colors Pasirinktinis	COLORS	[table widget]	Use the table widget if you want to choose the relief colors manually. You can add as many color classes as you want: for each class you can choose the lower and upper bound and finally by clicking on the color row you can choose the color thanks to the color widget. Fig. 24.36 Manually setting of relief color classes The buttons in the right side panel give you the chance to: add or remove color classes, change the order of the color classes already defined, open an existing file with color classes and save the current classes as file.
Relief	OUTPUT	[rastras] Numatytas: [Įrašyti į laikiną failą]	Specify the output relief raster layer. One of: Įrašyti į laikiną failą Įrašyti į failą…
Frequency distribution Pasirinktinis	FREQUENCY_DISTRIBUTION	[vektorius: lentelė] Numatytas: [Praleisti išvestį]	Specify the CSV table for the output frequency distribution. One of: Praleisti išvestį Įrašyti į laikiną failą Įrašyti į failą…

Išvestys

Užrašas	Pavadinimas	Tipas	Aprašymas
Relief	OUTPUT	[rastras]	The output relief raster layer
Frequency distribution	OUTPUT	[vektorius: lentelė]	The output frequency distribution

Pythono kodas

Algorithm ID: qgis:relief

import processing
processing.run("algorithm_id", {parameter_dictionary})

algoritmo id rodomas, kai užvedate pelę virš algoritmo apdorojimo įrankinėje. parametrų žodynas teikia parametrų pavadinimus ir reikšmes. Daugiau informacijos apie tai, kaip vykdyti apdorojimo algoritmus Pythono konsolėje rasite skyriuje Using processing algorithms from the console.

24.1.17.6. Ruggedness index

Calculates the quantitative measurement of terrain heterogeneity described by Riley et al. (1999). It is calculated for every location, by summarizing the change in elevation within the 3x3 pixel grid.

Each pixel contains the difference in elevation from a center cell and the 8 cells surrounding it.

Fig. 24.37 Ruggedness layer from low (red) to high values (green)

Parametrai

Užrašas	Pavadinimas	Tipas	Aprašymas
Aukščio sluoksnis	INPUT	[rastras]	Skaitmeninio paviršiaus modelio rastro sluoksnis
Z koeficientas	Z_FACTOR	[skaičius: slankaus kablelio] Numatytas: 1.0	Vertical exaggeration. This parameter is useful when the Z units differ from the X and Y units, for example feet and meters. You can use this parameter to adjust for this. Increasing the value of this parameter will exaggerate the final result (making it look more rugged). The default is 1 (no exaggeration).
Ruggedness	OUTPUT	[rastras] Numatytas: [Įrašyti į laikiną failą]	Specify the output ruggedness raster layer. One of: Įrašyti į laikiną failą Įrašyti į failą…

Išvestys

Užrašas	Pavadinimas	Tipas	Aprašymas
Ruggedness	OUTPUT	[rastras]	The output ruggedness raster layer

Pythono kodas

Algorithm ID: qgis:ruggednessindex

import processing
processing.run("algorithm_id", {parameter_dictionary})

algoritmo id rodomas, kai užvedate pelę virš algoritmo apdorojimo įrankinėje. parametrų žodynas teikia parametrų pavadinimus ir reikšmes. Daugiau informacijos apie tai, kaip vykdyti apdorojimo algoritmus Pythono konsolėje rasite skyriuje Using processing algorithms from the console.

24.1.17.7. Slope

Calculates the slope from an input raster layer. The slope is the angle of inclination of the terrain and is expressed in degrees.

Fig. 24.38 Flat areas in red, steep areas in blue

Parametrai

Užrašas	Pavadinimas	Tipas	Aprašymas
Aukščio sluoksnis	INPUT	[rastras]	Skaitmeninio paviršiaus modelio rastro sluoksnis
Z koeficientas	Z_FACTOR	[skaičius: slankaus kablelio] Numatytas: 1.0	Vertical exaggeration. This parameter is useful when the Z units differ from the X and Y units, for example feet and meters. You can use this parameter to adjust for this. Increasing the value of this parameter will exaggerate the final result (making it steeper). The default is 1 (no exaggeration).
Slope	OUTPUT	[rastras] Numatytas: [Įrašyti į laikiną failą]	Specify the output slope raster layer. One of: Įrašyti į laikiną failą Įrašyti į failą…

Išvestys

Užrašas	Pavadinimas	Tipas	Aprašymas
Slope	OUTPUT	[rastras]	The output slope raster layer

Pythono kodas

Algorithm ID: qgis:slope

import processing
processing.run("algorithm_id", {parameter_dictionary})

algoritmo id rodomas, kai užvedate pelę virš algoritmo apdorojimo įrankinėje. parametrų žodynas teikia parametrų pavadinimus ir reikšmes. Daugiau informacijos apie tai, kaip vykdyti apdorojimo algoritmus Pythono konsolėje rasite skyriuje Using processing algorithms from the console.