成長日記

3月 14日的日記[實驗記錄]
實驗日誌:(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年

成長日記