ST_坡度
上一頁  11.14. 網格處理：DEM（高程）  下一頁

名稱

ST_坡度 — 返回高程網格圖層的坡度（預設單位為度）。適用於地形分析。

概要

raster ST_坡度(raster rast, integer nband=1, text pixeltype=32BF, text units=DEGREES, double precision scale=1.0, boolean interpolate_nodata=FALSE);

raster ST_坡度(raster rast, integer nband, raster customextent, text pixeltype=32BF, text units=DEGREES, double precision scale=1.0, boolean interpolate_nodata=FALSE);

說明

返回高程網格圖層的坡度（預設單位為度）。利用地圖代數並將坡度方程式應用於相鄰的像素。

units 表示坡度的單位。可能的值為：RADIANS（弧度）、DEGREES（度，預設值）、PERCENT（百分比）。

scale 是垂直單位與水平單位的比例。對於英尺：經緯度，使用 scale=370400；對於公尺：經緯度，使用 scale=111120。

如果 interpolate_nodata 為 TRUE，則輸入網格圖層中 NODATA 像素的值將在使用 ST_InvDistWeight4ma 計算表面坡度之前進行內插。

[Note]

有關坡度、坡向和陰影的更多資訊，請參閱 ESRI - 陰影如何運作 和 ERDAS 現場指南 - 坡度影像。

可用性：2.0.0

增強功能：2.1.0 使用 ST_MapAlgebra() 並新增了可選的 units、scale、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_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)

範例：變體 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_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;