説明

標高ラスタバンドの傾斜方向(デフォルトの単位は度)を返します。地図代数を利用て、傾斜方向方程式を隣接ピクセルに適用します。

unitsは、傾斜方向の単位を示します。取りえる値はRADIANS, DEGREES(デフォルト)です。

unitsがRADIANSの時、値は0から2πラジアンの間で、北から時計回りに計ります。

unitsがDEGREESの時、値は0から360度の間で、北から時計回りに計ります。

ピクセルの勾配が0の場合には、傾斜方向は-1とします。

	勾配、傾斜方向、陰影起伏に関する詳細情報については、ESRI - How hillshade worksとERDAS Field Guide - Aspect Imagesをご覧下さい。

初出: 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;
                                        

関連情報

ST_MapAlgebra, ST_TRI, ST_TPI, ST_Roughness, ST_HillShade, ST_Slope