成長日記

4月 8日的日記[實驗記錄]
實驗日誌: [實驗] 崩塌深度模型: 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】 [論文] *坡度與深度不會有直接關聯* *坡度與深度不會有直接關聯* *坡度與深度不會有直接關聯*

成長日記