Name
ST_HillShade — 使用提供的方位角、海拔高度、亮度和比例输入返回高程栅格波段的假设照明。
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);
描述
使用方位角、海拔高度、亮度和比例输入返回高程栅格标注栏的假设照明。利用地图代数并将山体阴影方程应用于相邻像素。返回的像素值介于0和255之间。
azimuth 是从北向顺时针测量的介于0度和360度之间的值。
altitude 是介于0到90度之间的值，其中0度位于地平线上，90度位于正上方。
max_bright 是介于0和255之间的值，0表示没有亮度，255表示最大亮度。
scale 是垂直单位与水平单位的比率。对于英尺：拉特隆使用比例=370400，对于米：拉特龙使用比例=111120。
如果 interpolate_nodata 为真，则输入栅格中的NODATA像素值将使用 ST_InvDistWeight4ma 在计算山体阴影照度之前。
![[Note]](images/note.png) | |
有关山体阴影的更多信息，请参阅 山体阴影的工作原理 。 |
可用性：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
Coverage平铺的完整示例。此查询仅适用于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;