ST_HillShade
Prev  11.14. Rasterverarbeitung: DEM (Elevation)  Next

Name

ST_HillShade — Gibt für gegebenen Horizontalwinkel, Höhenwinkel, Helligkeit und Maßstabsverhältnis die hypothetische Beleuchtung eines Höhenrasterbandes zurück.

Synopsis

raster ST_HillShade(raster rast, integer band=1, text pixeltype=32BF, double precision azimuth=315, double precision altitude=45, double precision max_bright=255, double precision scale=1.0, boolean interpolate_nodata=FALSE);

raster ST_HillShade(raster rast, integer band, raster customextent, text pixeltype=32BF, double precision azimuth=315, double precision altitude=45, double precision max_bright=255, double precision scale=1.0, boolean interpolate_nodata=FALSE);

Beschreibung

Gibt für gegebenen Horizontalwinkel, Höhenwinkel, Helligkeit und Maßstabsverhältnis die hypothetische Beleuchtung eines Höhenrasterbandes zurück. Verwendet Map Algebra und wendet die Schummerungsgleichung auf die Nachbarpixel an. Die zurückgegebenen Pixelwerte liegen zwischen 0 und 255.

azimuth der Richtungswinkel zwischen 0 und 360 Grad, im Uhrzeigersinn von Norden gemessen.

altitude der Höhenwinkel zwischen 0 und 90 Grad, wobei 0 Grad am Horizont und 90 Grad direkt über Kopf liegt.

max_bright ein Wert zwischen 0 und 255, wobei 0 keine Helligkeit und 255 maximale Helligkeit bedeutet.

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 Schummerung berechnet wird.

[Note]

Für weitere Information zur Schummerung siehe How hillshade works.

Verfügbarkeit: 2.0.0

Erweiterung: 2.1.0 Verwendet ST_MapAlgebra() und der optionale Funktionsparameter interpolate_nodata wurde hinzugefügt

Änderung: 2.1.0 In Vorgängerversionen wurden Richtungswinkel und Höhenwinkel in Radiant angegeben. Nun werden die Werte in Grad ausgedrückt.

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_Hillshade(rast, 1, '32BF'))
FROM foo

                                                                                                                       st_dumpvalues

------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
-----------------------------------------------------------------------
 (1,"{{NULL,NULL,NULL,NULL,NULL},{NULL,251.32763671875,220.749786376953,147.224319458008,NULL},{NULL,220.749786376953,180.312225341797,67.7497863769531,NULL},{NULL,147.224319458008
,67.7497863769531,43.1210060119629,NULL},{NULL,NULL,NULL,NULL,NULL}}")
(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_Hillshade(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;