演算法說明

參考文獻(網址)：

http://desktop.arcgis.com/zh-cn/arcmap/10.3/tools/spatial-analyst-toolbox/how-aspect-works.htm

計算公式說明：

坡向分析是參考ArcGIS網站上所提供之演算法，用於計算每個像元高程值與鄰近像元方向變化率最大的下坡方向。

計算公式如下：

a	b	c
d	e	f
g	h	i

公式1：[dz/dx] = ((c + 2f + i) - (a + 2d + g)) / 8----------像元e在x方向的變化率計算
公式2：[dz/dy] = ((g + 2h + i) - (a + 2b + c)) / 8---------- 像元e在y方向的變化率計算
公式3：aspect = 57.29578 * atan2 ([dz/dy], -[dz/dx])---------- 像元e在y方向的變化率計算
公式4：
if aspect < 0
    cell = 90.0 - aspect
else if aspect > 90.0
    cell = 360.0 - aspect + 90.0
else
    cell = 90.0 - aspect---------- 再將坡向值轉換成為0-360度的方向