淺談基于GIS的區(qū)域土壤侵蝕遙感快速調(diào)查技術

1、準»（sl190 - 96）的規(guī)定，土壤侵蝕強度分級，必須以年平均侵 蝕模數(shù)為判別指標，只有缺少實測及調(diào)查侵蝕模數(shù)資料時，可 以在經(jīng)過分析后，運用有關侵蝕方式的指標進行分級。鑒于在 我國大部分地區(qū)暫時還難以獲取足夠的侵蝕模數(shù)資料,一般運 用右關侵蝕方式的指標進行分級，常用的方法是多因素綜合 法，即全而考慮影響土壤侵蝕的各個因素，在分析各個因素與 土壤侵蝕量關系的基礎上，確定土壤侵蝕的強度。下面以水力 侵蝕強度分級為例,簡述第二次全國土壤侵蝕遙感調(diào)查屮基于 gis的多因素綜合法操作技術。1 土壤侵蝕調(diào)查的技術路線在本次全國土壤侵蝕調(diào)查中，土壤侵蝕強度分級采用多因 素綜合法，操作中存在截

2、然相反的兩種過程，技術路線分別見 圖1和圖2。第一種方法,首先將區(qū)域劃分為若干基本評價單 元（或網(wǎng)格），然后按基本單元分析影響土壤侵蝕各因索的屬性 （因子值），再綜合各因子值判別基本單元的侵蝕強度;第二種 方法，首先分析影響土壤侵蝕的因子值，并分別作為一個層面, 然后將各個層面疊加，得到疊加圖斑和與其對應的各因子值, 再按照一定規(guī)則進行制圖綜臺，得到評價單元及與其對應的各 因子值，綜合評價各因子值確定侵蝕強度。實踐證明，兩種方 法都是行之有效的，其結果是客觀的和可靠的。2影響土壌侵蝕因子值的提取土壤侵蝕強度分級,需要得到植被、地形（地而坡度）和土 地利用類型等3個因子值。下而就植被信息和地形信

3、息的提 取技術進行闡述（土地利用類型信息已從其他單位得到）。 211植被覆蓋度獲取方法由于綠色植物的光譜特征是在近紅 外波段呈反射峰，在紅光波段呈反射谷，故一般都采用這兩個波段來研究植被。在這 兩個波段組成的平面上，土壤背景的光譜數(shù)據(jù)基本上沿著與這 兩個軸呈45的直線分布，這條線被稱作“土壤線”，沿著“土壤 線”移向坐標原點,土壤濕度増加，因此采用近紅外波段（nir）第二步:給評價單元賦因子值|土壤侵蝕因子狀況圖斑號面積被瀕植覆地面 坡度土地利 用類型 12第三步:綜合判定侵蝕強度i土壤侵蝕強度圖i第一種方法技術路線圖林草覆蓋度地面坡度第三步:確定土壤侵蝕強度圖斑號面積被度蓋植覆土地利 用類型

4、 土填侵 蝕強度1n圖2第二種方法技術路線圖被指數(shù)，刃表不為rvi = nir/ r ,其屮nir表示近紅外波段的亮度值，r表示紅波段的亮度值。 不同傳感器的nir和r的選取如表1。表1不同傳感器比值植被指數(shù)波段選擇傳感器nir波段r波段tm4,5 ,73mss6,75avhrr21雖然簡單比值形式能夠用于綠色植物的相對測量，但地點 和軌道不同的植被指數(shù)可能有相當大的差異。(2) nqw植被指數(shù)。ndvi被指數(shù)(歸一化差值植被指 數(shù))被定義為紅波段和近紅外波段亮度值的差與和的比值即ndyi = (nir - r)/(nir + r)。地點和軌道 不同的ndv/植被指數(shù)值相對穩(wěn)定，并且對于植被覆

5、蓋相對較好的地區(qū)(如植被覆蓋度大于60 %的地區(qū))有 著很好的指示。當植被覆蓋度frl對較小時，n"/可能變成負 值，為對這種情形進行修正，捉出轉(zhuǎn)換植被指數(shù)的概念，即將 nd刃值增加()15后取其平方根，其表達式為tvi= (ndvi + 0. 5丿 1/2 o(3) kauth和thomas植被指數(shù)。該植被指數(shù)應歸屬于波段 線性組合的植被指數(shù)，是基于地物光譜統(tǒng)計特征設計的一種方 法。對于mss圖像其定義為gvi = - 0.29mss4 - 0. 562 mss5 +0.6mss6 +0.491 mss1;對于tm圖像其gvi定義為gv1 = - 0. 2848 tm + 2. 2

6、435 tm2 + 0. 5436 tm3+ 0. 7243 tm4 + 0. 0840 tm5 - 0. 1800 tm1 o(4) 土壤糾正植被指數(shù)。該植被指數(shù)消除了“土壤線”方向 上的濕度變化。在w/r和r波段組成的平而上，土壤背景的 遙感光譜基本上分布在一個橢圓形的區(qū)域內(nèi)。橢惻形的長軸 就是所謂的“土壤線”，它代表土壤濕度的變化。該植被指數(shù)被 定義為sav1 =(nir r)/(nir + r + l) j x( + l),式中l(wèi)為修正因子。由于土壤糾正植被指數(shù)消除了“土壤線” 方向上的濕度變化，同時進行了等值線的偏轉(zhuǎn)矯正，所以其精 度相對要比其他指數(shù)高。在處理過程中，為了計算機處理方

7、便和考慮到可能的負 值，對于植被覆蓋度的提取,可以采用修正nd刃植被指數(shù)，形 式如下：(0r > n1rndvi ( nr r)/( nir + r + l) j xk r < nir 式中r為方便擴張值。212地面坡度獲取技術除利用傳統(tǒng)的坡度圖、坡度尺量算坡度外，在g1s軟件支 度值 撚丿m將二.角形網(wǎng)轉(zhuǎn)化為貝有拓撲天糸的多邊形數(shù)據(jù)， 后將矢量的多邊形數(shù)據(jù)轉(zhuǎn)化為柵格坡度數(shù)據(jù)。這就是進行 析處理所需要的數(shù)據(jù)。這個過程可在arc/info軟件屮，通 擴展模塊tin來完成。tin繪專門的地表模型生成、顯示和 析模型，可以計算坡度、坡向及其他數(shù)據(jù)。(2)由dem模型生成坡度。dem (地

8、面數(shù)字高程模型) 地形的數(shù)學描述,數(shù)據(jù)存儲多為等距網(wǎng)格模式(矩陣形式)。 dem可以派主坡度、坡向、地勢圖等。利用dem計算坡度 算法，按照所用高程的個數(shù)分為4點法、8點法和9點法，按 計算方法分為算數(shù)平均值法、加權平均法、絕對值最大法、曲 擬合法心。各種算法的合理性、可靠性和適應性各具特色， 根據(jù)應用目的來選用。也可通過dem數(shù)據(jù)光照模擬處理求 平均坡度，即用光照亮度的-次微分判斷坡度，-次微分絕 值大則坡度較陡，反之則坡度較緩，若一次微分為零則表示 底或山脊，即為坡向轉(zhuǎn)折點'。3結語(1) 基于g1s的區(qū)域土壤侵蝕遙感調(diào)查的快速性。快速 表現(xiàn)在兩個階段：其一是在獲得土壤侵蝕強度分級

9、本底值 (一次調(diào)查)階段°利用gis技術，既可以在遙感信息中提 取響土壤侵蝕的因子值,又能夠分析評價單元(圖斑)的幾何 屬(如位崟和而積)，這樣較手工作業(yè)更有效、快速地分類分 級計、查詢侵蝕面積等。其二是在數(shù)據(jù)更新階段。在gis支 下，對有變化或變化較大的局部進行分析和修改，即可得到 局的現(xiàn)狀數(shù)據(jù)。實際上，基于gis的區(qū)域土壤侵蝕遙感調(diào)査 術的快速性史充分地表現(xiàn)在數(shù)據(jù)史新階段。(2) 快速性的實質(zhì)和基礎。影響土壤侵蝕因子值是通過 析遙感影像數(shù)據(jù)、地形圖、dem或其他數(shù)據(jù)提取的°因此，快 調(diào)查的基礎是具有良好的基礎數(shù)據(jù)，實質(zhì)是挖掘數(shù)據(jù)得到賽 的信息。這就啟示我們，為了快速查清水土流失現(xiàn)狀、分析 土流失發(fā)展趨勢，必須長期不懈地采集和積累全面的而不是 面的、系統(tǒng)的而不是零碎的、準確的而不是錯誤的數(shù)據(jù)。(3) 根據(jù)侵蝕類型調(diào)整因子提取方法。本文主要敘述了 力侵蝕分級指標獲取技術和侵蝕強度判別的技術路線。在 別風蝕、凍融侵蝕強度時，應根據(jù)分級指標，在數(shù)據(jù)源、數(shù)捧參考文獻i卜兆宏，唐萬龍,李士鴻遙感技術與應用，1993 ,8(1) :