Ważne

Tłumaczenie jest wysiłkiem społeczności QGISa przyłącz się. Ta strona jest obecnie przetłumaczona w 62.50%.

24.1.17. Raster terrain analysis

24.1.17.1. Ekspozycja

Calculates the aspect of the Digital Terrain Model in input. The final aspect raster layer contains values from 0 to 360 that express the slope direction, starting from north (0°) and continuing clockwise.

../../../../_images/aspect.png

Rys. 24.30 Aspect values

The following picture shows the aspect layer reclassified with a color ramp:

../../../../_images/aspect_2.png

Rys. 24.31 Aspect layer reclassified

Parametry

Etykieta

Nazwa

Typ

Opis

Elevation layer

INPUT

[raster]

Digital Terrain Model raster layer

Z factor

Z_FACTOR

[numeric: double]

Domyślnie: 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. The default is 1 (no exaggeration).

Ekspozycja

OUTPUT

[raster]

Domyślnie: [Zapisz do pliku tymczasowego]

Specify the output aspect raster layer. One of:

  • Zapisz do pliku tymczasowego

  • Zapisz do pliku…

Wyniki

Etykieta

Nazwa

Typ

Opis

Ekspozycja

OUTPUT

[raster]

The output aspect raster layer

Kod pythona

ID algorytmu: qgis:aspect

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

The algorithm id is displayed when you hover over the algorithm in the Processing Toolbox. The parameter dictionary provides the parameter NAMEs and values. See Using processing algorithms from the console for details on how to run processing algorithms from the Python console.

24.1.17.2. DTM filter (slope-based)

Can be used to filter a digital elevation model in order to classify its cells into ground and object (non-ground) cells.

The tool uses concepts as described by Vosselman (2000) and is based on the assumption that a large height difference between two nearby cells is unlikely to be caused by a steep slope in the terrain. The probability that the higher cell might be non-ground increases when the distance between the two cells decreases. Therefore the filter defines a maximum height difference (dz_max) between two cells as a function of the distance (d) between the cells (dz_max( d ) = d). A cell is classified as terrain if there is no cell within the kernel radius to which the height difference is larger than the allowed maximum height difference at the distance between these two cells.

The approximate terrain slope (s) parameter is used to modify the filter function to match the overall slope in the study area (dz_max( d ) = d * s). A 5 % confidence interval (ci = 1.65 * sqrt( 2 * stddev )) may be used to modify the filter function even further by either relaxing (dz_max( d ) = d * s + ci) or amplifying (dz_max( d ) = d * s - ci) the filter criterium.

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

Zobacz także

This tool is a port of the SAGA DTM Filter (slope-based)

Parametry

Parametry podstawowe

Etykieta

Nazwa

Typ

Opis

Warstwa wejściowa

INPUT

[raster]

Digital Terrain Model raster layer

Band number

BAND

[raster band]

The band of the DEM to consider

Kernel radius (pixels)

RADIUS

[numeric: integer]

Domyślnie: 5

The radius of the filter kernel (in pixels). Must be large enough to reach ground cells next to non-ground objects.

Terrain slope (%, pixel size/vertical units)

TERRAIN_SLOPE

[numeric: double]

Default: 30.0

The approximate terrain slope in %. The terrain slope must be adjusted to account for the ratio of height units vs raster pixel dimensions. Used to relax the filter criterium in steeper terrain.

Filter modification

FILTER_MODIFICATION

[lista]

Domyślnie: 0

Choose whether to apply the filter kernel without modification or to use a confidence interval to relax or amplify the height criterium.

  • 0 - None

  • 1 - Relax filter

  • 2 - Amplify

Standard deviation

STANDARD_DEVIATION

[numeric: double]

Domyślnie: 0.1

The standard deviation used to calculate a 5% confidence interval applied to the height threshold.

Output layer (ground)

Opcjonalne

OUTPUT_GROUND

[raster]

Domyślnie: [Zapisz do pliku tymczasowego]

Specify the filtered DEM containing only cells classified as ground. One of:

  • Pomiń wyjście

  • Zapisz do pliku tymczasowego

  • Zapisz do pliku…

Output layer (non-ground objects)

Opcjonalne

OUTPUT_NONGROUND

[raster]

Default: [Skip output]

Specify the non-ground objects removed by the filter. One of:

  • Pomiń wyjście

  • Zapisz do pliku tymczasowego

  • Zapisz do pliku…

Parametry zaawansowane

Etykieta

Nazwa

Typ

Opis

Opcje generowania

Opcjonalne

CREATE_OPTIONS

[ciąg tekstowy]

Domyślnie: «»

For adding one or more creation options that control the raster to be created (colors, block size, file compression…). For convenience, you can rely on predefined profiles (see GDAL driver options section).

Batch Process and Model Designer: separate multiple options with a pipe character (|).

Wyniki

Etykieta

Nazwa

Typ

Opis

Output layer (ground)

OUTPUT_GROUND

[raster]

The filtered DEM containing only cells classified as ground.

Output layer (non-ground objects)

OUTPUT_NONGROUND

[raster]

The non-ground objects removed by the filter.

Kod pythona

ID algorytmu: native:dtmslopebasedfilter

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

The algorithm id is displayed when you hover over the algorithm in the Processing Toolbox. The parameter dictionary provides the parameter NAMEs and values. See Using processing algorithms from the console for details on how to run processing algorithms from the Python console.

24.1.17.3. Cieniowanie

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.

../../../../_images/azimuth.png

Rys. 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.

../../../../_images/hillshade1.png

Rys. 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:

../../../../_images/hillshade_2.png

Rys. 24.34 Overlapping the hillshade with the elevation layer

Parametry

Etykieta

Nazwa

Typ

Opis

Elevation layer

INPUT

[raster]

Digital Terrain Model raster layer

Z factor

Z_FACTOR

[numeric: double]

Domyślnie: 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

[numeric: double]

Domyślnie: 300.0

Set the horizontal angle (in degrees) of the sun (clockwise direction). Range: 0 to 360. 0 is north.

Vertical angle

V_ANGLE

[numeric: double]

Domyślnie: 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).

Cieniowanie

OUTPUT

[raster]

Domyślnie: [Zapisz do pliku tymczasowego]

Specify the output hillshade raster layer. One of:

  • Zapisz do pliku tymczasowego

  • Zapisz do pliku…

Wyniki

Etykieta

Nazwa

Typ

Opis

Cieniowanie

OUTPUT

[raster]

The output hillshade raster layer

Kod pythona

ID algorytmu: qgis:hillshade

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

The algorithm id is displayed when you hover over the algorithm in the Processing Toolbox. The parameter dictionary provides the parameter NAMEs and values. See Using processing algorithms from the console for details on how to run processing algorithms from the Python 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.

Parametry

Etykieta

Nazwa

Typ

Opis

DEM to analyze

INPUT_DEM

[raster]

Digital Terrain Model raster layer to use for calculating altitudes

Boundary layer

BOUNDARY_LAYER

[vector: polygon]

Polygon vector layer with boundaries of areas used to calculate hypsometric curves

Krok

STEP

[numeric: double]

Domyślnie: 100.0

Vertical distance between curves

Use % of area instead of absolute value

USE_PERCENTAGE

[wartość logiczna]

Domyślnie: False

Write area percentage to “Area” field of the CSV file instead of the absolute area

Hypsometric curves

OUTPUT_DIRECTORY

[folder]

Specify the output folder for the hypsometric curves. One of:

  • Zapisz do katalogu tymczasowego

  • Zapisz w katalogu

Wyniki

Etykieta

Nazwa

Typ

Opis

Hypsometric curves

OUTPUT_DIRECTORY

[folder]

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.

../../../../_images/hypsometric.png

Kod pythona

ID algorytmu: qgis:hypsometriccurves

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

The algorithm id is displayed when you hover over the algorithm in the Processing Toolbox. The parameter dictionary provides the parameter NAMEs and values. See Using processing algorithms from the console for details on how to run processing algorithms from the Python console.

24.1.17.5. Rzeźba terenu

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.

../../../../_images/relief.png

Rys. 24.35 Relief layer

Parametry

Etykieta

Nazwa

Typ

Opis

Elevation layer

INPUT

[raster]

Digital Terrain Model raster layer

Z factor

Z_FACTOR

[numeric: double]

Domyślnie: 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

[wartość logiczna]

Domyślnie: False

If you check this option the algorithm will create all the relief color classes automatically

Relief colors

Opcjonalne

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.

../../../../_images/relief_table.png

Rys. 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.

Rzeźba terenu

OUTPUT

[raster]

Domyślnie: [Zapisz do pliku tymczasowego]

Specify the output relief raster layer. One of:

  • Zapisz do pliku tymczasowego

  • Zapisz do pliku…

Frequency distribution

Opcjonalne

FREQUENCY_DISTRIBUTION

[vector: table]

Default: [Skip output]

Specify the CSV table for the output frequency distribution. One of:

  • Pomiń wyjście

  • Zapisz do pliku tymczasowego

  • Zapisz do pliku…

Wyniki

Etykieta

Nazwa

Typ

Opis

Rzeźba terenu

OUTPUT

[raster]

The output relief raster layer

Frequency distribution

OUTPUT

[vector: table]

The output frequency distribution

Kod pythona

ID algorytmu: qgis:relief

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

The algorithm id is displayed when you hover over the algorithm in the Processing Toolbox. The parameter dictionary provides the parameter NAMEs and values. See Using processing algorithms from the console for details on how to run processing algorithms from the Python 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.

../../../../_images/ruggedness.png

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

Parametry

Etykieta

Nazwa

Typ

Opis

Elevation layer

INPUT

[raster]

Digital Terrain Model raster layer

Z factor

Z_FACTOR

[numeric: double]

Domyślnie: 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).

Szorstkość

OUTPUT

[raster]

Domyślnie: [Zapisz do pliku tymczasowego]

Specify the output ruggedness raster layer. One of:

  • Zapisz do pliku tymczasowego

  • Zapisz do pliku…

Wyniki

Etykieta

Nazwa

Typ

Opis

Szorstkość

OUTPUT

[raster]

The output ruggedness raster layer

Kod pythona

ID algorytmu: qgis:ruggednessindex

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

The algorithm id is displayed when you hover over the algorithm in the Processing Toolbox. The parameter dictionary provides the parameter NAMEs and values. See Using processing algorithms from the console for details on how to run processing algorithms from the Python console.

24.1.17.7. Nachylenie

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

../../../../_images/slope3.png

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

Parametry

Etykieta

Nazwa

Typ

Opis

Elevation layer

INPUT

[raster]

Digital Terrain Model raster layer

Z factor

Z_FACTOR

[numeric: double]

Domyślnie: 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).

Nachylenie

OUTPUT

[raster]

Domyślnie: [Zapisz do pliku tymczasowego]

Specify the output slope raster layer. One of:

  • Zapisz do pliku tymczasowego

  • Zapisz do pliku…

Wyniki

Etykieta

Nazwa

Typ

Opis

Nachylenie

OUTPUT

[raster]

The output slope raster layer

Kod pythona

ID algorytmu: qgis:slope

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

The algorithm id is displayed when you hover over the algorithm in the Processing Toolbox. The parameter dictionary provides the parameter NAMEs and values. See Using processing algorithms from the console for details on how to run processing algorithms from the Python console.