Name
ST_Slope — Gibt die Neigung (standardmäßig in Grad) eines Höhenrasterbandes zurück. Nützlich für Terrain-Analysen.
Synopsis
raster ST_Slope(raster rast, integer nband=1, text pixeltype=32BF, text units=DEGREES, double precision scale=1.0, boolean interpolate_nodata=FALSE);
raster ST_Slope(raster rast, integer nband, raster customextent, text pixeltype=32BF, text units=DEGREES, double precision scale=1.0, boolean interpolate_nodata=FALSE);
Beschreibung
Gibt die Neigung (standardmäßig in Grad) eines Höhenrasterbandes zurück. Verwendet Map Algebra und wendet die Neigungsgleichung auf benachbarte Pixel an.
units die Einheiten für die Neigung. Mögliche Werte sind: RADIANS, DEGREES (Standardeinstellung) und Prozent.
scale ist das Verhältnis zwischen vertikalen und horizontalen Einheiten. Für Feet:LatLon verwenden Sie bitte scale=370400, für Meter:LatLon scale=111120.
Wenn interpolate_nodata TRUE ist, dann werden die NODATA Pixel des Ausgangsraster mit ST_InvDistWeight4ma interpoliert, bevor die Neigung der Oberfläche berechnet wird.
![[Note]](../images/note.png) | |
Für weitere Information über Neigung, Exposition und Schummerung siehe ESRI - How hillshade works und ERDAS Field Guide - Slope Images. |
Verfügbarkeit: 2.0.0
Erweiterung: 2.1.0 Verwendet ST_MapAlgebra() und es wurden die optionalen Funktionsparameter units, scale und interpolate_nodata hinzugefügt
Änderung: 2.1.0 In Vorgängerversionen wurden die Werte in Radiant ausgegeben. Nun werden die Werte standardmäßig in Grad ausgegeben.
Beispiele: Variante 1
WITH foo AS (
SELECT ST_SetValues(
ST_AddBand(ST_MakeEmptyRaster(5, 5, 0, 0, 1, -1, 0, 0, 0), 1, '32BF', 0, -9999),
1, 1, 1, ARRAY[
[1, 1, 1, 1, 1],
[1, 2, 2, 2, 1],
[1, 2, 3, 2, 1],
[1, 2, 2, 2, 1],
[1, 1, 1, 1, 1]
]::double precision[][]
) AS rast
)
SELECT
ST_DumpValues(ST_Slope(rast, 1, '32BF'))
FROM foo
st_dumpvalues
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------
(1,"{{10.0249881744385,21.5681285858154,26.5650520324707,21.5681285858154,10.0249881744385},{21.5681285858154,35.2643890380859,36.8698959350586,35.2643890380859,21.5681285858154},
{26.5650520324707,36.8698959350586,0,36.8698959350586,26.5650520324707},{21.5681285858154,35.2643890380859,36.8698959350586,35.2643890380859,21.5681285858154},{10.0249881744385,21.
5681285858154,26.5650520324707,21.5681285858154,10.0249881744385}}")
(1 row)
Beispiele: Variante 2
Ein vollständiges Beispiel mit Coveragekacheln. Diese Abfrage funktioniert nur mit PostgreSQL 9.1 oder höher.
WITH foo AS (
SELECT ST_Tile(
ST_SetValues(
ST_AddBand(
ST_MakeEmptyRaster(6, 6, 0, 0, 1, -1, 0, 0, 0),
1, '32BF', 0, -9999
),
1, 1, 1, ARRAY[
[1, 1, 1, 1, 1, 1],
[1, 1, 1, 1, 2, 1],
[1, 2, 2, 3, 3, 1],
[1, 1, 3, 2, 1, 1],
[1, 2, 2, 1, 2, 1],
[1, 1, 1, 1, 1, 1]
]::double precision[]
),
2, 2
) AS rast
)
SELECT
t1.rast,
ST_Slope(ST_Union(t2.rast), 1, t1.rast)
FROM foo t1
CROSS JOIN foo t2
WHERE ST_Intersects(t1.rast, t2.rast)
GROUP BY t1.rast;