實驗日誌:
[實驗]
崩塌深度模型:
Soil depth =
(1 - a *G - b* Driver/Driver_max + c *Curvature + d* TWI/TWI_max)^e
G = slope gradient
【->直接用Slope工具做出的分級(degree)即可】
Driver = is the relative distance to the river channel on the slope
【->使用Euclidean工具,不過”搜尋半徑”要用不同尺度都測試看看,有什麼差異】
Curvature = Profile curvature
【->使用Curvature工具,產出的三種曲率中,選擇Profile curvature】
TWI = topographic wetness index
【-> "oxoxox"為圖層名
○1用DEM做"flow direction" (PS:沿用之前做子水系時的圖層)
○2用flow direction做"flow accumulation"(PS1:沿用之前做子水系時的圖層。PS2: flow accumulation_1000圖層是當初為了從flow accumulation圖層產製水系,所設定的下一步閾值圖層(1000),本次分析應使用原本的flow accumulation)
○3用DEM做"slope"(PS:沿用之前做的圖層)
○4使用raster calculator得到slope的弧度角(radian)
[公式] Tan("slope" * 1.570796) / 90 = "slope radian" (PS: 1.570796為tanθ值)
○5使用raster calculator得到tan(slope)
[公式] Con("slope radian"= =0 , 0.001 , "slope radian")) = "tan(slope)"
○6使用raster calculator得到圖層▲
[公式] ("flow accumulation" + 1) * "某一年度DEM(cell size(m))" = "▲" (PS1:cell size為分析時,使用哪一張尺度的DEM的網格大小。PS2:▲代表圖層名,無特別意義。)
○7使用raster calculator
[公式] Ln("▲" / "tan(slope)") = TWI
(所有PS1:圖層兩邊都要加上” “。PS2:數學函數Tan、Cos…等等,開頭都要大寫,GIS才會真正叫出函數)
】
a,b,c,d,e = 常數
【-> 需透過迭代方式找出最適合解,迭代的input為上述模型的output】
[論文]
*坡度與深度不會有直接關聯*
*坡度與深度不會有直接關聯*
*坡度與深度不會有直接關聯*