Name
ST_Aspect — 返回高程栅格波段的坡向（默认以度为单位）。 对于分析地形很有用。
Synopsis
raster ST_Aspect(raster rast, integer band=1, text pixeltype=32BF, text units=DEGREES, boolean interpolate_nodata=FALSE);
raster ST_Aspect(raster rast, integer band, raster customextent, text pixeltype=32BF, text units=DEGREES, boolean interpolate_nodata=FALSE);
描述
返回高程栅格波段的坡向（默认以度为单位）。 利用地图代数并将长宽比方程应用于相邻像素。
units 表示坡度的单位。可能的取值有：RADIANS（弧度），DEGREES（度，默认值）。
当units = RADIANS 时，值介于 0 到 2 * pi 弧度之间，从北向顺时针测量。
当units = DEGREES 时，值介于 0 到 360 度之间（从北向顺时针测量）。
如果像素的斜率为零，则像素的纵横比为-1。
![[Note]](../images/note.png)
|
|
|
有关坡度、坡向和山体阴影的更多信息，请参阅 ESRI - 山体阴影的工作原理和 ERDAS 现场指南 - 坡向图像。 |
可用性: 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_Aspect(rast, 1, '32BF'))
FROM foo
st_dumpvalues
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
----------------------------------
(1,"{{315,341.565063476562,0,18.4349479675293,45},{288.434936523438,315,0,45,71.5650482177734},{270,270,-1,90,90},{251.565048217773,225,180,135,108.434951782227},{225,198.43495178
2227,180,161.565048217773,135}}")
(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_Aspect(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;