ST_Aspect
Prev  11.14. 래스터 공간 처리  Next

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π 라디안 사이의 값입니다.

units = DEGREES일 경우, 북쪽에서 시계방향으로 측정한 0도와 360도 사이의 값입니다.

픽셀의 경사가 0일 경우, 픽셀의 향은 -1입니다.

[Note]

경사(slope), 향(aspect), 음영기복(hillshade)에 대한 자세한 내용을 알고 싶다면, ESRI - How hillshade works 및 ERDAS Field Guide - Aspect Images 를 참조하십시오.

2.0.0 버전부터 사용할 수 있습니다.

개선 사항: 2.1.0 버전부터 ST_MapAlgebra()를 이용하며, 선택적인 interpolate_nodata 함수 파라미터가 추가됐습니다.

변경 사항: 2.1.0 미만 버전에서는 반환되는 값이 라디안 단위였습니다. 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;