ST_HillShade
上一个  11.14. 栅格处理：DEM（高程）  下一个

名称

ST_HillShade — 使用提供的方位角、高度、亮度和比例输入，返回高程栅格波段的假想照明。

概要

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);

描述

使用方位角、高度、亮度和比例输入，返回高程栅格波段的假想照明。 利用地图代数并将山体阴影方程应用于相邻像素。 返回的像素值介于 0 和 255 之间。

azimuth 是一个介于 0 到 360 度之间的值，从正北开始顺时针测量。

altitude 是一个介于 0 到 90 度之间的值，其中 0 度在地平线上，90 度在正上方。

max_bright 是一个介于 0 和 255 之间的值，其中 0 表示无亮度，255 表示最大亮度。

scale 是垂直单位与水平单位的比率。 对于英尺：经纬度，使用 scale=370400，对于米：经纬度，使用 scale=111120。

如果 interpolate_nodata 为 TRUE，则在计算山体阴影照明之前，将使用 ST_InvDistWeight4ma 插值输入栅格的 NODATA 像素的值。

[Note]

有关山体阴影的更多信息，请参阅 山体阴影的工作原理。

可用性：2.0.0

增强功能：2.1.0 使用 ST_MapAlgebra() 并添加了可选的 interpolate_nodata 函数参数

已更改：2.1.0 在早期版本中，方位角和高度以弧度表示。 现在，方位角和高度以度表示

示例：变体 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)

示例：变体 2

覆盖范围的瓦片的完整示例。 此查询仅适用于 PostgreSQL 9.1 或更高版本。

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;