實驗日誌:(03.14~03.21)
[實驗] 2010~2013 崩塌面積資料計算
[論文]看兩篇崩塌面積公式文章
[操作流程]:
因為將每一年的面積代入已設定好的迴歸式,再回推關係式參數,這邏輯不太對。因為原本就是回歸式的係數了,再回推不就一模一樣 = = ?
先嘗試用2012、2016兩年的DEM乘上面積,然後看求出來的2012、2016體積與用既定迴歸式得出來的數值差多少?
使用Calculate Geometry算出來的面積與Shape_Area面積最一致,所以捨棄Area的數據,各年面積都用Shape_Area
polygon to raster應該不能配合DEM的最小size(20 m),這樣轉起來網格面積會太大,得到的面積與高度相乘結果,應該很不準。 *LS的面積我都轉成1m,再乘上DEM (下周問看看老師這樣對不對?) *
將做完的LS_2012年_PolygonToRaster_1m做reclassify,把網格數值都改成1
把2012、2016年所有DEM做Mosaic,共選取29張DEM。兩年度的DEM圖幅編碼都一樣
(編號:972230XX )
(34,35,36,42,43,44,45,46,52,53,54,55,56,57,63,64,65,66,67,73,74,75,76,77,83,84,85,93,94) 計29張
20mDEM resample成1公尺網格,使用Nearest分類方法
(方法無法選擇Majority,因為Majority不支援float格式的資料)
計算2012年的崩塌體積,做Raster Calculator。
"LS_2012年_PolygonToRaster_1m_Reclass" * "Mosaic_DEM_2012年_Resample_1m" = 2012年崩塌體積
*錯,上面這個計算是錯的。用崩塌地直接乘DEM,會直接得到超大數值,那是從水平面開始算,乾 差點變成智障*
看來得採用DEM轉坡度->深度法,查表
(驗證過Mosaic_20米與Resample_1米,網格數值都相同,那可以直接用1米的做Slope)
( *用1m的DEM做slope會崩圖,不知道原因為何,下週問一下老師* )
做Slope, 用20米DEM做
將Slope20米做Resample成1米
再做一次Reclassify,將slope按照文獻對應的崩塌深度,分成四個組別
(<30, 30~40, 40~60, >60)
先測試NCDR(2004)的版本
計算2012年的崩塌體積,做Raster Calculator。
"Slope_2012年_1m_Reclass_崩塌深度_NCDR" * "LS_2012年_PolygonToRaster_1m_Reclass" = 2012年崩塌體積
用既定迴歸式,計算2012年的崩塌體積(Excel表格)
製作表格,整理各崩塌深度方法與迴歸式方法所得的體積值
https://imgur.com/KJbIsXV
將上述步驟重複乙次,製作完2016年