利用数值地型模型进行自动化集水区河川网路分析.ppt

利用數值地型模型進行自動化集水區河川網路分析與林地分級之適切度進度報告 授課老師:朱子豪 教授 學生:莊永忠 2004.1.5 研究動機 本研究旨在探討數值地形資料應用於自動化萃取集水區區劃與河川網路分析之可行性 。 將不同網格解析度與門檻值套用於數值地形模型進行河川演算時，會由於不同地形或是地表而產生不同表現，而各種不同的網格解析度與門檻值設定並未實際套用於森林密度高、地形起伏大之台灣山區森林集水區， 故必須進行檢驗。 計算出不同網格解析度與門檻值設定下演算出之集水區分級分區結果，並套用於林地分級系統，並與人工萃取之集水區套用結果做比較， 如此就可以找出最接近人工萃取集水區之自動化演算方式，未來便可找出替代人工區劃與萃取之工作流程。 研究目的 依據網格式水文演算法基礎理論（D8 法）在ArcInfo下撰寫AML程式，以實際樣區研究比較在不同網格解析度、不同之門檻值設定下演算集水區水流網路與區劃之表現。 進行自動化演算後求得河川網路與集水區區劃，首先利用測量認可度（Cohen’s kappa 統計量數）來檢驗與人工化作業所得河川網路之比較誤差，找出與人工進行區劃萃取最接近之網格解析度與最適合之門檻值。其次將自動化與人工萃取之結果套用至林地分級系統中，並配合相同之地形因子，最後比較其分級結果差異。 數值地型模型之介紹 數值高程網格之介紹 在座標系統上用網格的格式來記錄地表的高程，每一網格的值即代表該網格中心點之高程。 數值地形模型之生產方式 航空照片立體像對測量（航照立體測圖儀）。 野外現場直接測量 。 從現有的等高線地形圖數化 。 網格資料的優劣分析 運算上簡易有效率。 取樣上之複雜度與資料生產困難度亦較低。 一般研究所取用之資料來源多為網格式數值地形模型，資料取得較易，在應用發展上用途較廣。 運算之結果受制於固定網格寬度之大小，無法根據不同之地形地貌特徵進行修正。 應用於水文網路分析時，其演算結果為正方形網格組成之鋸齒狀結構，與真實之河流網路不全符合。 由於網格解析度之限制，無法輕易地掌握高度突然變化之地形以及一些微小之地物。 數值地形模型應用於水文模擬研究發展進程 以網格式數值地形模型作為資料來源，不僅可利用電腦程式自動進行河流網路粹取，亦可減少傳統人工在等高線地形圖上手繪河流水系時所產生之人為主觀判斷誤差，在資料儲存方面則可有效保存河川網路演算結果。 由於網格式數值地形模型資料特徵獨特，促使研究人員必須針對不同地形地貌使用不同解析度之網格資料，並選用不同區間之門檻值進行試誤模擬，以求得與真實世界適切度最高之河川網路。 研究試區、研究材料與儀器設備 台北縣市陽明山區 試區總面積約882公頃(原始區劃登記面積約835公頃) 資料格式為4M*4M之數值地形模型(經過smooth處理) 儀器設備主要為裝有Arc/Info地理資訊系統與EMDS林地分級系統之個人電腦 利用數值地形模型自動粹取河川網路與集水區之流程與方法 自動化水文模擬演算法（D8法）之原理介紹 由O’Callahan and Mark (1984)所提出 。 水流沿著地表面八個方向（0、π/2、π、3π/2、π/4、3π/4、5π/4、7π/4）中最陡的坡度流動的假設 ，每一網點流向的考慮是根據其與八個臨點的高度差來決定 。 根據水流方向，進而求得每一網點上游所累積的入流網點的數量（累積流量） ，再由累積流量的大小進一步測定河谷的位置。 自動化水文模擬演算法工作流程 窪地填補 FILL in_grid out_grid {SINK | PEAK} {z_limit} {out_dir_grid} 流向分析 FLOWDIRECTION(surface_grid,{o_drop_grid}, {NORMAL | FORCE}) 自動化水文模擬演算法之工作流程 累積流量計算 ( 設定網格權重為1 ) FLOWACCUMULATION(dir_grid, {weight_grid}) 門檻值設定 自動化水文模擬演算法之工作流程 自動化水文模擬演算法之工作流程 自動化水文模擬演算法之工作流程 自動化水文模擬演算法之工作流程 ARC/INFO GRID 河道網路擷取作業流程 目前之進度報告 目前已完成數值地形資料之過濾處理(利用加權平均法進行3*3視窗平滑化) 與自動水文模擬演算，整個流程約為50分鐘。 對照保安林資料庫中河流線的結果，發現對照40M*40M網格式數值地形模型，門檻值設定約80左右最符合實際的水系分布。 在進行自動萃取時，發現會產生數個問題如下 : 1.河道失真的問題 2.未歸類集水區之問題(細小碎塊) 3.網格與向量