數(shù)值高度模型之地形復(fù)雜度量度指標(biāo)研究-以蝕溝等級為例

1、1數(shù)值高度模型之地形複雜度量度指標(biāo)研究以蝕溝等級為例 論文口試 報告人：許秋玲指導(dǎo)老師：朱子豪教授2002/5/81101001000110100Log (Distance of two points)Log (Variance)y = 1.6615x + 0.508R2 = 0.9911. 81. 922. 12. 22. 32. 490603003303002702402101801501201.9912.1522.1632.2062.1912.081360=2.1542報告大綱1.緒論 1.1 1.1 前言前言 1.2 1.2 研究動機(jī)與目的研究動機(jī)與目的1.3 1.3 研究定位與範(fàn)疇研

2、究定位與範(fàn)疇2.文獻(xiàn)回顧 2.1 地形複雜度 2.2 簡確法之蝕溝等級量度 2.2 數(shù)值高度模型量度 2.3 地形特徵值之量度2.3.2 2.3.2 數(shù)值地形的資料解析數(shù)值地形的資料解析度、度、 抽樣間距抽樣間距 2.4 2.4 碎形參數(shù)碎形參數(shù)3.研究方法 3.1 3.1 研究架構(gòu)研究架構(gòu) 3.2 研究流程4.DEM地形複雜度量度實驗設(shè)計 4.1 4.1 地形複雜度實驗設(shè)計流地形複雜度實驗設(shè)計流程程5.DEM地形複雜度量度分析 5.1 5.1 分析結(jié)果分析結(jié)果 5.1.1 5.1.1 數(shù)學(xué)模擬表面分?jǐn)?shù)學(xué)模擬表面分析析 5.1.2 5.1.2 真實地表複雜度真實地表複雜度分析分析 不規(guī)則土地單

3、元不規(guī)則土地單元 矩形單元矩形單元 5.2 碎形參數(shù)與傳統(tǒng)地形參數(shù)之間的 關(guān)係6.結(jié)論與討論 6.1 6.1 結(jié)論結(jié)論 6.2 6.2 未來研究方向未來研究方向3 張石角（1987）簡確法之蝕溝等級1.1 前言統(tǒng)計型態(tài)量度統(tǒng)計型態(tài)量度地形量度簡確法蝕溝等級幾何型態(tài)量度幾何型態(tài)量度碎形分析Low DLow RHigh DLow RHigh RLow DHigh DHigh R 14801/224WavelengthAmplitude45454514727638376 碎形的興起主要是源自地圖概括化，持續(xù)應(yīng)用在其他的空間現(xiàn)象上，碎形維度提供地圖概括化的工具以及預(yù)測與抽樣空間有關(guān)的資料分割特性。簡確

4、法中的蝕溝等級之判定方式，係依照等高線彎曲的角度，這種判定方式也是一種等高線地形圖概括化的過程，要連結(jié)真實地表的複雜現(xiàn)象以及概括化的等級判定方式，可用碎形的概念表示地表的複雜程度。 目前碎形的研究多利用特別的演算法描述個別物件的特性，但是在描述地形這樣複雜的結(jié)構(gòu)時，無法找到特定的數(shù)學(xué)演算法來描述地形的形貌，也無法完全控制地形演變的因素，僅能利用統(tǒng)計的方法來計算碎形維度，以及其地形的碎形特性。 4 1.2 研究 與目的 幾何型態(tài)分析 數(shù)值模式的分析 解決傳統(tǒng)以等高線地形圖從事研究時所遭遇的問題，並且以一套可操作的程序所得到的結(jié)果，在不同的時間、不同的操作者的情況下得以重複驗證。 對於地形複雜度這

5、種抽象的概念如能夠以量化分析的方式，取代簡確法中依據(jù)等高線彎曲角度來判定蝕溝等級之方式。 碎形模擬地形碎形參數(shù)所代表的地形意義地形地貌的形態(tài)十分複雜，用傳統(tǒng)歐氏幾何學(xué)去描述地形則過於簡化，但如果利用碎形的函數(shù)則可創(chuàng)造出十分逼真的自然形態(tài)。用人工所模擬的地形看起來與真實地形相似， 卻無法進(jìn)一步解釋其意義。SlopeAspect5 1.2 研究 與目的本研究之範(fàn)疇屬於地形量測的一環(huán)，探討DEM中的地形複雜度量度，以簡確法之蝕溝等級所呈現(xiàn)的複雜度為量度的對象。對地表的形態(tài)（surface morphology）作一定量分析，並且引用碎形分析方式中的空間自相關(guān)統(tǒng)計概念，探討DEM之高程變化的空間分布特

6、性與傳統(tǒng)等高線地形圖中所呈現(xiàn)地形起伏資訊之異同。 為了了解如何在DEM的有限樣本點(diǎn)紀(jì)錄下，傳遞實際地表連續(xù)起伏變化的趨勢，並且與等高線的紀(jì)錄方式作一對應(yīng)， 本研究採用4m DEM與1：5,000等高線地形圖， 希冀達(dá)到下列具體的研究目標(biāo)： SlopeAspect6找出地形崎嶇度的有效量度方式探討碎形參數(shù)碎形維度值與水平的地形起伏頻率地形崎嶇度之關(guān)聯(lián)。 找出地形起伏度的有效量度方式探討碎形參數(shù)碎形截距值與水平的地形起伏幅度地形起伏度之關(guān)聯(lián)。 找出地形崎嶇度與起伏度的組合，量度蝕溝等級所呈現(xiàn)的地形複雜度。透過量度DEM之地形複雜度特徵值，找出DEM的碎形參數(shù)與地形圖上等高線的彎曲程度之關(guān)聯(lián)。 1.

7、2 研究動機(jī)與Curvature ProfileCurvature Plane71.3 研究 與範(fàn)疇創(chuàng)意層次引入數(shù)學(xué)與物理學(xué)領(lǐng)域中的碎形布朗表面來詮釋真實地表的空間自相關(guān)程度，並且應(yīng)用於蝕溝等級問題與考量地形複雜度因子組合，屬於組合運(yùn)用。 主題定位主題在於空間量度中地形複雜度的量度，以地形學(xué)為其應(yīng)用範(fàn)疇，因為傳統(tǒng)對於地表形態(tài)之描述多以抽象符號化的狀態(tài)陳述，本研究以定量研究方式尋找一綜合性指標(biāo)描述。 上下游定位屬於中游的技術(shù)研發(fā)層級，其下游為簡確法之評估應(yīng)用，上游則為數(shù)學(xué)之空間量度。 全面性定位本研究探討的地形複雜度量度方法的對象不限於簡確法中的土地單元，其產(chǎn)出的指標(biāo)其目的為描述各種地形單元的地

8、形特徵值描述，故此種量度方法屬於通論的研究，選取簡確法中的蝕溝級序之量度為代表性驗證的對象。 產(chǎn)出形式產(chǎn)出的形式為一組狀態(tài)描述的指標(biāo)，這些指標(biāo)有數(shù)學(xué)統(tǒng)計上以及對應(yīng)到真實世界的意義。 主要應(yīng)用層次產(chǎn)出的實用性為地形特徵值的理論補(bǔ)強(qiáng)，可應(yīng)用至水文分析或是土壤沖蝕等研究 。 8 1.3 研究定位與空間量度地形複雜度量度理論面方法面碎形幾何變異元空間統(tǒng)計拓樸幾何地形計測參數(shù)統(tǒng)計意義蝕溝等級地形特徵應(yīng)用面92. 文獻(xiàn)回顧 指標(biāo)碎形參數(shù)起伏度102.3.2 數(shù)值地形的資料解析度、抽樣間距 Zhang & Montgometry（1994）使用2, 4, 10, 30, 90公尺解析度的DEM進(jìn)行地

9、形與水文應(yīng)用，得到10公尺解析度較30以及90公尺好，而與2以及4公尺的分析結(jié)果差異不大。 黃誌川、徐美玲（2002）探討不同的坡度計算方式在不同解析度（5, 10, 20, 40公尺）的DEM下的量度穩(wěn)定性，得到由DEM所計算出來的結(jié)果，各種演算法在解析度越粗的條件下，有低估坡度的趨勢而且有坡度平滑化的現(xiàn)象。不同的坡度計算方法中，緩坡如臺地及河谷，由DEM所計算出來的坡度變化較陡坡大，陡坡部分的坡度變化則較小。解析度20公尺與1/5,000等高線地形圖的地勢分析結(jié)果最吻合。 DEM重視高度變化的解像力與強(qiáng)調(diào)地表平面解像力的等高線不同，兩個抽樣方式的誤差依據(jù)地形的平坦或崎嶇而定。11碎形(fr

10、actals)理論 碎形維度是一種以自我相似為基礎(chǔ)的量測指標(biāo)，它可以適當(dāng)?shù)孛枋鲎匀唤缇拔锏难}雜程度，而被視為一種可以表達(dá)地表特性的量化參數(shù)（Ouchi & Matsushita, 1992; 黃金聰，1998）。碎形幾何學(xué)與傳統(tǒng)歐氏幾何學(xué)的比較(李宗仰，1999):2.5 碎形參數(shù)比較項目 歐氏幾何 碎形幾何 附註 歷史 傳統(tǒng)的（超過2000年） 近代的（約20年） 量度單位 使用特徵尺度或比例 無需特別尺度或大小 特徵尺度 適用範(fàn)圍 人為物件 自然形狀 描述語言 代數(shù)公式 疊代運(yùn)算 描述性 定量 定性 維度值 整數(shù)（不連續(xù)） 可為分?jǐn)?shù)（連續(xù)） 碎形維度細(xì)部結(jié)構(gòu)性 有限 無窮 自相似性

11、 D=1D=2D=3r=1r=2r=3 AB12尺度不變區(qū)定率碎形尺度不變區(qū)碎形特性定率碎形（定率碎形（Deterministic Fractal ）自我相似性（self-similarity）：不規(guī)則曲線的特徵外形在縮小尺度下會重複出現(xiàn)，而且尺度不斷縮小，長度會無限增加，但是碎形維度不變，滿足尺度不變性（scale-invariant）。其函數(shù)f（x, y）與改變尺度之f（rx, ry）是相似的，碎形維度不因尺度的改變而變化（包含x軸與y軸方向） 。 滿足連續(xù)而不可微（non-differentiable）：該曲線是連續(xù)不規(guī)則的曲線，在轉(zhuǎn)折點(diǎn)無法微分。隨機(jī)碎形隨機(jī)碎形（Random Frac

12、tal ）自然界的複雜現(xiàn)象建構(gòu)在統(tǒng)計中的隨機(jī)所產(chǎn)生的自我類似性（self-affine）的特性上，仍會受到尺度的影響。自類似性的碎形是由碎形布朗運(yùn)動（fractal Brownian motion, fBm）發(fā)展而來（Mandelbrot, 1982） 。其函數(shù)f（x, y）與改變尺度之f（rx, rHy）是相似的，碎形維度是沿著x軸的方向上，不因為垂直於x軸方向上的尺度之變化而改變，應(yīng)用在地形上其分析的尺度是地表水平的距離（ Mandelbrot, 1980; Turcotte, 1992 ）。 自類似性的碎形與物件表面起伏的頻率有關(guān)。 尺度不變區(qū)非碎形體系非碎形體系隨機(jī)碎形尺度不變區(qū)自相似

13、性的概念(李宗仰，1999 ) 133.1 研究架構(gòu)變異元法矩形單元蝕溝等級地形複雜度數(shù)值地形資料特性碎形參數(shù)地形崎嶇度空間自相關(guān)碎形截距碎形維度地形起伏度地形複雜度指標(biāo)人工判釋傳統(tǒng)地形計測(Geomorphometry)GIS數(shù)值高度模型量度碎形分析(Fractals)空間統(tǒng)計(Geostatistics)143.1.3 DEM 碎形維度方法的比較計算方法 公 式 碎形維度的計算說明 優(yōu)缺點(diǎn)盒計法(Box Counting) N(e)：非空集合的方盒總數(shù)e：方盒的邊長 繪製log e 對應(yīng)log N(e) 的座標(biāo)圖；迴歸直線斜率：1-D 由DEM粹取河川網(wǎng)路資料；封閉的曲線資料（海島的海岸線

14、、等高線等）餘長問題，需大量數(shù)據(jù)才可以作為合理的結(jié)果。能譜法(Power Spectrum) P(f)：頻譜密度f：頻率 繪製log f 對應(yīng)log P(f)的座標(biāo)圖；迴歸直線斜率：-(5-2D) 前期處理較為複雜；Herzfeld & Overbeck（1999 ）認(rèn)為頻譜法會誇大D值。 變異元法(Variogram) Zi,Zj：位於i點(diǎn)與j點(diǎn)的高程dij：i與j點(diǎn)間的距離 繪製log dij 對應(yīng)log Var (Zi-Zj) 的座標(biāo)圖；迴歸直線的斜率2：3-D 可沿特定方向取樣，來獲得各方向空間變異情形。 DEM常用計算碎形維度的方法 DeeN)3(2DijjidZZVar D

15、ffP2515碎形維度變異元法選定起始點(diǎn)取樣距離 du （通常設(shè)定為影像的解析度）計算兩點(diǎn)間的距離等於取樣矩離之兩點(diǎn)對的高程差 計算高程差的變方: Var（z（u+du）z（u） 逐漸增大取樣距離，重複上述步驟。 繪製以Log 取樣距離 為橫軸，Log 高程差的變方 為縱軸的座標(biāo)圖。 選擇最佳密合數(shù)據(jù)的直線，並以線性迴歸方式計算此直線的斜率。碎形維度等於3減去斜率的二分之一。 兩邊取log，此時的變異元的斜率（坡度）相當(dāng)於2H，因此碎形維度可以由下式求得:變異元法是應(yīng)用空間統(tǒng)計分析（geostatistical analysis）的方法來分析空間資料，以處理DEM資料為例，其演算的步驟如下 :

16、 HjijiuuuZuZ HjijiuuuZuZE22 2jiijuZuZEh Hijijhh2HED) 1(21slopeEDjiuu 碎形布朗運(yùn)動具有位置增量的平均值為零，以及增量的變方與時間差的2H次方成比例的關(guān)係，亦即具有增量不變性以及等向性的性質(zhì)（Voss, 1985） 。布朗粒子的運(yùn)動方式是假設(shè)布朗粒子的位置（u） 。H為常數(shù)且其值介於0與1之間，布朗運(yùn)動就是H1/2時的情形。設(shè)u位於E維的歐幾里德空間RE ，則D與H的關(guān)係D=E+1-H 。16碎形維度變異元法原始變異元實驗變異元（分組方式） PairsLag Distance(h)Variogram ()Scale(C)Nugg

17、et EffectSillExperimental Variogram CurveVarianceModel CurveRange21slopeED量度蝕溝是在分析地形的紋理性(texture)特性，故採用第一段迴歸線在R2 0.990的碎形參數(shù)值17真實地表的DEM地形複雜度量度各種量度方法的整理與比較等高線地形圖等高線彎曲角度地形特徵值數(shù)學(xué)幾何處理:一次微分;二次微分統(tǒng)計總和崎嶇度起伏度單點(diǎn)屬性特徵值坡度坡向曲率平均值 標(biāo)準(zhǔn)差斜率峰度地形複雜度特徵值影響地形複雜度量度因子面積地質(zhì)構(gòu)造坡度拓樸幾何碎形幾何變異元分析碎形參數(shù)高度Log-Log Variogram影響範(fàn)圍(R)碎形維度(D)碎形

18、截距()單元區(qū)界定組合蝕溝分級蝕溝分級R-squared值3.2 研究流程地形複雜度量度研究流程184. DEM 地形複雜度量度實驗設(shè)計變異元法蝕溝等級(真實地表)地形起伏度地形崎嶇度碎形截距碎形維度數(shù)學(xué)模型(人造地形)?樣區(qū)大小 掃描方向多蝕溝等級混合樣本過濾雜訊蝕溝、地形起伏與碎形參數(shù)關(guān)聯(lián)圖19混合等級部分與整體系統(tǒng)性偵測最大侵蝕 方向出現(xiàn)頻率垂直水平 形狀Relief RoughnessD理論數(shù)學(xué)模型規(guī)則單元不規(guī)則土地單元面積 坡度實驗設(shè)計階層圖20蝕溝等級(真實地表)起伏度崎嶇度?碎形截距碎形維度樣本變異元法過濾雜訊樣區(qū)大小掃描方向性多蝕溝等級混合數(shù)學(xué)模型(人造地形)根據(jù)蔡宗勳（199

19、4）所使用的曲面 模擬地形起伏，使用的數(shù)學(xué)方程式z=f(x,y)的方程式如下: H1 , H2 :改變數(shù)學(xué)模擬地形起伏度的參數(shù)L1 , L2 :改變數(shù)學(xué)模擬地形崎嶇度（地形頻率）的參數(shù)C :控制地形高度（ Z值）基準(zhǔn)面的參數(shù) 利用此正（餘）弦函數(shù)，控制地形波峰波谷出現(xiàn)的頻率，以及表現(xiàn)出凹坡、 平坡、與凸坡，產(chǎn)生各種不同的地形。CLYHLXHZ2211cossin5.1.1.數(shù)學(xué)模擬表面分析2 LHCH2LC21Reliefz=100*sin(x/400)+100*cos(y/400)+1800110100100010000100000100000010000000110100100010000

20、log(distance of two points)log(variance)z=500*sin(x/400)+500*cos(y/400)+1800110100100010000100000100000010000000110100100010000log(distance of two points)log(variance)z=2500*sin(x/400)+2500*cos(y/400)+1800110100100010000100000100000010000000110100100010000log(distance of two points)log(variance)H=10

21、0, L=400H=500, L=400H=2500, L=40018002500500022Roughnessz=100*sin(x/400)+100*cos(y/400)+1800110100100010000100000100000010000000110100100010000log(distance of two points)log(variance)z=50*sin(x/200)+50*cos(y/200)+1800110100100010000100000100000010000000110100100010000log(distance of two points)log(v

22、ar iance)z=25*sin(x/100)+25*cos(y/100)+1800110100100010000100000100000010000000110100100010000log(distance of two points)log(variance)H=100, L=400H=50, L=200H=25, L=100控制地形起伏度 控制地形崎嶇度 HLD HLD 1004002.1186 (0.991)1-20 0.0529 1004002.1186 (0.991)1-20 0.0529 2004002.1186 (0.991)1-20 0.2117 1002002.1350

23、 (0.991)1-10 0.2190 4004002.1186 (0.991)1-20 0.8470 1001002.1652 (0.993)1-5 1.0023 5004002.1186 (0.991)1-20 1.3234 100502.1836 (1.000)1-2 4.0504 8004002.1186 (0.991)1-20 3.3877 502002.1350 (0.991)1-10 0.0548 25004002.1186 (0.991)1-20 33.083 251002.1652 (0.993)1-5 0.0626 數(shù)學(xué)模型證明： 地形起伏度碎形截距地形崎嶇度碎形維度23問

24、題核心：研究區(qū)與使用資料方法採用：資料前處理、碎形之變異元計算樣本設(shè)定：單元區(qū)的界定、變異元參數(shù)的設(shè)定 蝕溝等級(真實地表)起伏度崎嶇度?碎形截距碎形維度樣本變異元法過濾雜訊樣區(qū)大小掃描方向性多蝕溝等級混合數(shù)學(xué)模型(人造地形)研究範(fàn)圍6464 (pixel)3232 (pixel)50個樣本 0306090120150180掃描方向研究範(fàn)圍1/5,000數(shù)化河系4m DEM 696集流閾值244公尺解析度的DEM資料輸入5.1.2 實驗設(shè)計流程影響地形複雜度量度因子坡度碎形參數(shù)碎形截距() 碎形維度(D) 影響範(fàn)圍(R)崎嶇度起伏度解析度地形參數(shù)蝕溝等級的量度是依照單元內(nèi)等高線彎曲的角度 人工

25、判釋空間統(tǒng)計變異元分析實驗變異元分組方式取樣點(diǎn)對數(shù)Log-Log variogramR-squared0.99最大點(diǎn)對距離角度分隔半徑分隔不同掃描方向(每隔30 )多蝕溝等級混合 部分與整體隨機(jī)掃描(360 )矩形單元251.81.922.12.22.32.490603003303002702402101801501202.03442.17892.10491.99272.00782.0329360=2.0044方位角 0 30 60 90 120 150 360 值 1.91601-25 1.03111-25 0.67551-25 1.53571-25 2.99851-25 2.99361-2

26、5 1.71041-25 SD(D)=0.6714、SD()=0.9322 5.1.2 真實地表複雜度分析不規(guī)則土地單元261.81.922.12.22.32.490603003303002702402101801501202.08242.09052.10262.172.15562.1135360=2.1402方位角 0 30 60 90 120 150 360 值 1.57111-17 1.79061-25 3.55711-25 3.95001-25 3.11531-25 1.23851-25 1.57111-17 SD(D)=0.0341、SD()=1.0911 5.1.2 真實地表複雜度

27、分析不規(guī)則土地單元271.81.922.12.22.32.490603003303002702402101801501202.11942.06742.09372.16332.18762.1963360=2.1499方位角 0 30 60 90 120 150 360 值 3.43721-21 3.54321-25 3.24491-25 3.08041-25 3.34891-11 3.33731-10 3.36671-25 SD(D)=0.0501、SD()=0.1519 5.1.2 真實地表複雜度分析不規(guī)則土地單元281.81.922.12.22.32.49060300330300270240

28、2101801501202.06922.11812.10532.1542.18062.1153360=2.1233方位角 0 30 60 90 120 150 360 值 3.35891-18 3.81941-25 3.92371-25 2.68471-25 2.70401-25 3.19371-13 3.17181-25 SD(D)=0.0371、SD()=0.5055 ItemGully ScaleD(360) (360) Max_D (方向)Max_(方向)SD (D)SD ()無蝕溝2.00441-251.71041-252.1789 (60)1-252.9985 (120)1-250

29、.67140.9322小蝕溝2.14031-171.57111-172.1700 (0)1-173.9500 (90)1-250.03411.0911中蝕溝2.14991-253.36671-252.1963 (120)1-113.5432 (30)1-250.05010.1519大蝕溝2.12331-253.17181-252.1806 (150)1-133.9237 (60)1-250.03710.5055 5.1.2 真實地表複雜度分析不規(guī)則土地單元/ D0.67550.31001.57110.72403.34891.52483.19371.4646蝕溝等級與Max_D方向之Max_D以

30、及蝕溝的頻率成正比Max_D的方向與最大地形起伏頻率的方向相同291.81.922.12.22.32.490603003303002702402101801501202.01952.01442.01352.06582.17722.033360=2.03181.81.922.12.22.32.490603003303002702402101801501202.03872.01951.99892.06912.28412.0735360=2.0478A1A3A2A430A5A6A7IDMax_D (direction)Max_ (direction)IDMax_D (direction)Max_ (

31、direction)A12.1772 (150) 1-25 0.8080 (60) 1-25 A52.2968 (180) 1-10 0.5953 (90) 1-25 A22.2169 (150)1-25 0.6374 (60)1-25 A62.3363 (120) 1-16 2.5536 (30) 1-25 A32.2841 (150) 1-25 0.5438 (60) 1-25 A72.2227 (150) 1-16 0.9473 (60) 1-25 A42.3218 (120) 1-10 0.5231 (30) 1-25 A82.3523 (120) 1-22 0.5617 (30) 1

32、-25 / D0.07800.03580.09900.04470.12170.05330.17600.0769/ D0.17660.07690.21080.09020.23600.10620.35480.1508A8311.81.922.12.22.32.490603003303002702402101801501202.06351.98531.97822.05962.14432.1442360=2.0557B1B21.81.922.12.22.32.490603003303002702402101801501202.08162.01471.98962.06012.19942.3304360=

33、2.0573B3B432B5B6B71.81.922.12.22.32.490603003303002702402101801501202.03762.05782.1128222952.32382.1447360=2.1097B8IDMax_D (direction)Max_ (direction)IDMax_D (direction)Max_ (direction)B12.1443 (150) 1-12 0.9333 (90) 1-25 B52.2211 (60) 1-12 0.9330(160) 1-25 B22.2739 (120) 1-25 1.1569(0) 1-25 B62.254

34、6 (30) 1-9 2.7939 (120) 1-25 B32.3304 (120) 1-11 0.8370 (0) 1-25 B72.3884 (120) 1-23 0.6973 (30)1-25 B42.2996 (150) 1-9 0.5014 (90) 1-25 B82.3238 (150) 1-6 1.0051 (30) 1-25 / D0.32870.15330.37750.16600.39360.16890.47950.2085/ D0.51320.23110.52890.23460.57410.24040.63800.2746331.81.922.12.22.32.49060

35、3003303002702402101801501202.01692.09872.11472.15852.09741.9911360=2.0818C1C2C31.81.922.12.22.32.42.590603003303002702402101801501202.08062.00892.02732.13162.40292.2397360=2.0987C4341. 81. 922. 12. 22. 32. 490603003303002702402101801501202.15242.16182.04122.03192.11322.1716360=2.1583C51.81.922.12.22

36、.32.490603003303002702402101801501202.19662.24462.07432.01582.06072.1879360=2.1246C61.81.922.12.22.32.490603003303002702402101801501202.08562.04042.02122.15562.25262.2389360=2.1180C7C8IDMax_D (direction)Max_ (direction)IDMax_D (direction)Max_ (direction)C12.1585 (180) 1-9 1.8030 (60) 1-25 C52.1716 (

37、120) 1-10 1.9158 (30) 1-25 C20.2525 (60) 1-14 2.2946 (120) 1-25 C62.2446 (60) 1-12 1.6255 (150) 1-25 C32.3665 (120) 1-9 1.2388 (150) 1-12 C72.2526 (150) 1-8 1.7124(120) 1-25 C42.4029 (150) 1-25 1.0120 (150) 1-25 C82.2137 (60) 1-17 2.6984 (0) 1-25 / D0.85190.39470.91750.40730.98420.41591.01200.4212/

38、D1.03710.47761.08120.48171.16650.51781.32130.5969351.81.922.12.22.32.490603003303002702402101801501202.17322.14682.03982.07322.18282.1873360=2.1474D11.81.922.12.22.32.490603003303002702402101801501202.20632.12382.08662.09692.16902.2561360=2.1465D2D31.81.922.12.22.32.490603003303002702402101801501202

39、.02882.18242.28442.23662.21612.0589360=2.1558D436D5D61.81.922.12.22.32.490603003303002702402101801501202.17532.192.08972.02862.08572.2007360=2.1781D7D8IDMax_D (direction)Max_ (direction)IDMax_D (direction)Max_ (direction)D12.1873 (120) 1-21 2.8249 (60) 1-25 D52.2642 (120) 1-23 1.9911 (120) 1-23 D22.

40、2561 (120) 1-7 1.8591(150) 1-25 D62.3178 (120) 1-25 2.5978 (120) 1-25 D32.2273 (150) 1-4 1.8442 (30) 1-25 D72.2007 (120) 1-13 2.6104 (60) 1-18 D42.2844 (30) 1-9 3.4962 (150) 1-25 D82.3289(150) 1-25 3.1703 (150) / D1.44180.65921.52310.67511.71710.77091.83650.8039/ D1.99110.87942.59781.12081.25001.193

41、63.17031.361337IDMax_D (direction)Max_ (direction)IDMax_D (direction)Max_ (direction)無蝕溝A12.1772 (150) 1-25 0.8080 (60) 1-25 A52.2968 (180) 1-10 0.5953 (90) 1-25 A22.2169 (150)1-25 0.6374 (60)1-25 A62.3363 (120) 1-16 2.5536 (30) 1-25 A32.2841 (150) 1-25 0.5438 (60) 1-25 A72.2227 (150) 1-16 0.9473 (6

42、0) 1-25 A42.3218 (120) 1-10 0.5231 (30) 1-25 A82.3523 (120) 1-22 0.5617 (30) 1-25 小蝕溝B12.1443 (150) 1-12 0.9333 (90) 1-25 B52.2211 (60) 1-12 0.9330(160) 1-25 B22.2739 (120) 1-25 1.1569(0) 1-25 B62.2546 (30) 1-9 2.7939 (120) 1-25 B32.3304 (120) 1-11 0.8370 (0) 1-25 B72.3884 (120) 1-23 0.6973 (30)1-25

43、 B42.2996 (150) 1-9 0.5014 (90) 1-25 B82.3238 (150) 1-6 1.0051 (30) 1-25 中蝕溝C12.1585 (180) 1-9 1.8030 (60) 1-25 C52.1716 (120) 1-10 1.9158 (30) 1-25 C20.2525 (60) 1-14 2.2946 (120) 1-25 C62.2446 (60) 1-12 1.6255 (150) 1-25 C32.3665 (120) 1-9 1.2388 (150) 1-12 C72.2526 (150) 1-8 1.7124(120) 1-25 C42.

44、4029 (150) 1-25 1.0120 (150) 1-25 C82.2137 (60) 1-17 2.6984 (0) 1-25 D12.1873 (120) 1-21 2.8249 (60) 1-25 D52.2642 (120) 1-23 1.9911 (120) 1-23 5.1.2 真實地表複雜度分析矩形單元/ D0.07800.03580.09900.04470.12170.05330.17600.07690.32870.15330.37750.16600.39360.16890.47950.20850.85190.39470.91750.40730.98420.41591.

45、01200.42121.44180.65921.52310.67511.71710.77091.83650.8039/ D0.17660.07690.21080.09020.23600.10620.35480.15080.51320.23110.52890.23460.57410.24040.63800.27461.03710.47761.08120.48171.16650.51781.32130.59691.99110.87942.59781.12081.25001.19363.17031.3613蝕溝等級與Max_D方向之Max_D成正比蝕溝軸線相同或是無蝕溝的Max_D與Max_ 方向垂

46、直381.81.922.12.22.32.490603003303002702402101801501202.19332.17332.05142.03852.05272.0882360=2.0568 1.81.922.12.22.32.490603003303002702402101801501202.22762.16732.04392.02432.04502.0901360=2.05391.81.922.12.22.32.490603003303002702402101801501202.17982.23132.15072.12092.11742.0804360=2.1057 D (360) (360)Max_D （方向）Max_（方向） Max_D 方向之 DSD (D)SD ()左2.0541-25 2.337 2.228(90)1-21 3.946(150)1-25 0.4511 0.0772 1.1700 中2.0571-25 2.930 2.193(90)1-22 5.143(150)1-25 0.5399 0.0641 1.5633 右2.1061-25 3.037 2.231(60)1-25 5