陳燕瓊

（福州大學福建省空間信息工程研究中心空間數(shù)據(jù)挖掘與信息共享教育部重點實驗室，福建 福州350002）

0 引言

土壤侵蝕是全球性的環(huán)境問題之一，制約著經(jīng)濟、社會、環(huán)境的可持續(xù)發(fā)展。土壤侵蝕降低了土地生產(chǎn)力，直接影響了農(nóng)業(yè)的生產(chǎn)發(fā)展，開展土壤侵蝕強度研究并提出有效的治理措施具有一定的現(xiàn)實意義。然而，大范圍的土壤侵蝕強度監(jiān)測需要耗費大量的時間、人力和經(jīng)費。因此，在一定精度范圍內(nèi)，得到滿足科學計算要求的土壤侵蝕估算已成為一種趨勢。長期以來，國內(nèi)外學者[1-4]對土壤侵蝕規(guī)律及估算模型的研究，促進了土壤侵蝕強度估算和水土保持工作的推進。李忠佩等[5]分析了紅壤丘陵區(qū)水土流失的過程，并提出了相應的技術措施，為水土流失區(qū)綜合治理獻策。土壤侵蝕定量評價方法主要以土壤侵蝕估算模型為核心，其中，最為著名、應用最廣泛的是美國修正通用土壤流失方程（RUSLE，Revised Universal Soil Loss Equation）。國內(nèi)一些學者通過對RUSLE的影響因子，如降雨侵蝕力值[6]、土壤可蝕性[7]、植被覆蓋度[8]等的研究，為RUSLE 在國內(nèi)相應區(qū)域的應用奠定了理論基礎。同時，GIS 和RS技術的發(fā)展，為土壤侵蝕的定量估算提供了強有力的支持。Arekhi[9]、李輝[10]、徐涵秋[11-12]、林惠花[13]等運用地理信息系統(tǒng)和遙感技術，對研究區(qū)土壤侵蝕強度進行了估算，以了解土壤的侵蝕現(xiàn)狀，并提出相應的治理措施。以上的研究都是在已有軟件工具的基礎上運用土壤侵蝕估算模型進行分析研究，這對于非專業(yè)用戶來說尚有困難。計算結(jié)果以平面圖的形式展示，未充分發(fā)揮三維可視化的優(yōu)勢。因此，根據(jù)用戶需求開發(fā)數(shù)據(jù)處理工具，簡化土壤侵蝕評價的計算邏輯，可以有效提高工作效率。充分發(fā)揮三維可視化功能的優(yōu)勢，可以方便用戶更好地理解土壤侵蝕發(fā)生機理。

鑒于上述研究背景，結(jié)合COM 技術，在自主研發(fā)的虛擬森林環(huán)境平臺VisForest[14]上，集成修正通用土壤流失方程模型，構(gòu)建具有良好可操作性、能夠服務于土壤侵蝕評價的子系統(tǒng)。

1 系統(tǒng)設計

1.1 總體設計思想

土壤侵蝕估算涉及到大量空間數(shù)據(jù)和非空間數(shù)據(jù)的采集、處理與利用，土壤侵蝕評價與可視化子系統(tǒng)涉及數(shù)據(jù)輸入、因子計算分析、數(shù)據(jù)輸出、結(jié)果可視化4個主要功能。數(shù)據(jù)輸入主要包括矢量數(shù)據(jù)、遙感影像數(shù)據(jù)和DEM 數(shù)據(jù)的輸入。因子計算分析是子系統(tǒng)的核心模塊，為保證計算結(jié)果的準確性，需要選擇合理的參數(shù)計算公式。各個評價因子計算完成后，通過疊加等空間分析運算輸出。由于土壤侵蝕評價因子只是一個半定量的計算結(jié)果，為了定量反映土壤侵蝕強度，系統(tǒng)需要基于RUSLE 土壤侵蝕預報模型進行設計，估算土壤侵蝕量?？紤]到空間分析的方便性，數(shù)據(jù)輸出格式主要為柵格數(shù)據(jù)。結(jié)果可視化主要包括二維和三維兩種形式，將計算得到的二維結(jié)果圖疊加到三維地形上，提高可讀性。土壤侵蝕量的計算和土壤侵蝕結(jié)果可視化能夠初步反映土壤流失過程的特征，用于土地利用和土壤侵蝕的關系研究，為土地利用格局的優(yōu)化設計和水土流失綜合治理提供技術支持。

1.2 結(jié)構(gòu)設計

根據(jù)土壤侵蝕評價模塊的設計思想、計算邏輯和ArcGIS Engine 的技術特點，模塊結(jié)構(gòu)設計（見圖1）采用三層結(jié)構(gòu)模式：①數(shù)據(jù)層主要用來提供和管理所需的空間和非空間數(shù)據(jù)；②邏輯層是模塊功能實現(xiàn)的技術支持，通過VC++語言調(diào)用ArcGIS Engine 中的組件，構(gòu)建土壤侵蝕評價模塊的結(jié)構(gòu)框架；③表現(xiàn)層主要包括土壤侵蝕評價模塊的土壤侵蝕量計算和土壤侵蝕結(jié)果可視化兩個模塊。

圖1 土壤侵蝕評價與可視化系統(tǒng)的結(jié)構(gòu)框架

2 功能實現(xiàn)

基于ArcGIS Engine 的開發(fā)，在VC++開發(fā)環(huán)境中，實現(xiàn)土壤侵蝕評價與可視化的解決方案。其核心是在自主研發(fā)的虛擬森林環(huán)境平臺VisForest 上，定義插件遵循的接口，以及自定義組件接口的實現(xiàn)。模塊建設的技術流程如下：在工程中引用ArcObjects 對象庫；采用智能指針聲明所需接口，在類模塊中寫入實現(xiàn)特定ArcObject 接口的代碼；定義過程函數(shù)；添加初始化事件Oncreate、Click 等代碼；編譯成DLL；加載到VisForest 平臺中，根據(jù)這一過程構(gòu)建土壤侵蝕評價模塊。

2.1 土壤侵蝕量計算子模塊

本模塊主要是依據(jù)RUSLE 模型建立，并根據(jù)中國水土保持的特點和前人研究的結(jié)果對參數(shù)因子計算公式進行了修正。RUSLE 的數(shù)學表達公式為：式中：A為土壤侵蝕模數(shù)；R為降雨量侵蝕力指標；K為土壤可蝕性因子；LS、C、P分別為地形因子、水土保持因子、植被覆蓋與管理因子。

其主要功能是通過RUSLE 土壤侵蝕預報模型估算區(qū)域的土壤侵蝕量。根據(jù)RUSLE 模型涉及的參數(shù)因子設計了6 個模塊，分別為降雨侵蝕因子R、土壤可蝕性因子K、坡度因子S、坡長因子L、作物覆蓋與管理因子C、水土保持措施因子P、土壤侵蝕量計算。通過這些子模塊，可以分別獲得各計算因子的結(jié)果圖層，最后通過土壤侵蝕量計算子模塊的疊加，計算輸出土壤侵蝕量圖層。

2.2 土壤侵蝕量結(jié)果可視化子模塊

在自主研發(fā)的虛擬森林環(huán)境系統(tǒng)（VisForest）平臺基礎上，針對土壤侵蝕分析的業(yè)務需求，集成修正通用土壤流失方程（RUSLE），形成土壤侵蝕評價與可視化子系統(tǒng)。將遙感影像和計算得到的土壤侵蝕量圖層與三維地形進行疊加，形成土壤侵蝕強度三維場景如圖2 所示，在圖2中，可定位到關心的區(qū)域，查看該區(qū)域的地形地貌、土地覆蓋/利用等，結(jié)合土壤侵蝕強度圖進行分析。

圖2 土壤侵蝕三維場景

3 應用實例

本文主要數(shù)據(jù)源為：長汀縣周邊20 個氣象觀測站1971—2012 年的各月總降雨數(shù)據(jù)；長汀縣DEM數(shù)據(jù)（分辨率為30m）；長汀縣已有的2010年的土地利用類型圖和1∶5萬土壤類型矢量圖；2010年歸一化植被覆蓋指數(shù)分布圖。

首先，使用上述因子的取值及計算方法對數(shù)據(jù)進行預處理，然后，將各因子圖進行柵格連乘計算，即可獲得每個柵格的土壤侵蝕模數(shù)。按照《土壤侵蝕強度分級標準》對計算所得的長汀縣土壤侵蝕強度進行歸并，得到長汀縣2010 年土壤侵蝕強度的面積、侵蝕總量分布及不同侵蝕強度侵蝕量對侵蝕總量的貢獻率，見表1。

表1 2010年土壤侵蝕級別統(tǒng)計表

由表1可知：

（1）2010年水土流失的侵蝕總量為369.37萬t，平均土壤侵蝕模數(shù)為11.92t·hm-2·a-1，屬于輕度侵蝕；其與2003 年水土流失平均土壤侵蝕模數(shù)27.03t·hm-2·a-1相比[15]，年平均土壤侵蝕強度下降明顯；

（2）長汀全縣侵蝕面積為39 995hm2，估算結(jié)果符合當?shù)厍闆r[16]，其與1985 年（97 469hm2）和2000 年（70 364 hm2）相比[17]，土壤侵蝕總面積大大減少，其中，微度侵蝕面積24 937hm2，占侵蝕總面積的2.35%；輕度侵蝕面積13 217hm2，占侵蝕總面積的33.08%；中、強和極強度分別占侵蝕總面積的3.56%、0.77%、0.15%；劇烈侵蝕面積達37hm2，占侵蝕總面積的0.01%；治理成效卓越與長汀縣多年來大力扶持的政策及一系列有效措施的實施有著密切聯(lián)系。

地形的起伏變化是造成長汀縣內(nèi)土壤侵蝕空間差異性大的關鍵因素。統(tǒng)計長汀縣不同坡度下的土壤侵蝕模數(shù)（見表2）發(fā)現(xiàn)，隨著坡度的增加，土壤侵蝕模數(shù)的變化總體呈現(xiàn)遞增趨勢，坡度大于25°的侵蝕模數(shù)最大，其次是20～25°和15～20°，15°是土壤侵蝕模數(shù)的相對質(zhì)變點。將土壤侵蝕強度圖和遙感影像疊加在三維地形上形成三維場景圖（見圖3），圖3（a）、圖3（b）內(nèi)的矩形框區(qū)域在實地為坡度大于15°的果園，其植被覆蓋分布均勻，但矩形框內(nèi)的土壤侵蝕強度相比周邊坡度小于15°的區(qū)域，其土壤侵蝕明顯更為嚴重。

表2 不同坡度的土壤侵蝕

圖3 坡度對土壤侵蝕影響的三維表達

4 結(jié)論

本文借助組件式GIS技術，在自主研發(fā)的虛擬森林環(huán)境系統(tǒng)（VisForest）平臺基礎上，集成修正通用土壤流失方程（RUSLE），以可視化形式表達土壤侵蝕空間分布特征。本文的創(chuàng)新點在于自主研發(fā)的虛擬森林環(huán)境系統(tǒng)（VisForest）平臺集成了GIS和RUSLE，實現(xiàn)了土壤侵蝕的評價功能，簡化土壤侵蝕估算步驟。將計算得到的土壤侵蝕結(jié)果疊加在三維地形上，充分發(fā)揮三維可視化表達的優(yōu)勢，方便用戶分析導致土壤侵蝕嚴重的可能因素。

[1] Yu Guoqiang,Li Zhanbin.Study on Mechanism of Regulation Measures on Water Erosion Process in the Slope System[C]// Electric Technology and Civil Engineering(ICETCE).Lushan：2011 International Conference on.ed., 2011：2636-2639.

[2] Shi Z H, Fang N F, Wang L, et al. Soil Erosion Processes and Sediment Sorting Associated with Transport Mechanisms on Steep Slopes[J]. Journal of Hydrology, 2012，(454)：123-130.

[3] Yuichi Onda, Hiroaki Kato, Maki Tsujimura, et al. Analysis of Runoff Generation and Soil Erosion Processes by Using Environmental Radionuclides in Semiarid Areas of Mongolia[J].Journal of Hydrology,2007,333(1)：124-132.

[4] 李銳.中國主要水蝕區(qū)土壤侵蝕過程與調(diào)控研究[J].水土保持通報，2011，31（5）：1-6.

[5] 李忠佩，張?zhí)伊郑瑮钇G生，等.紅壤丘陵區(qū)水土流失過程及綜合治理技術[J]. 水土保持通報，2001，21（2）：12-17.

[6] 周伏建，黃炎和.福建省降雨侵蝕力指標R值[J].水土保持學報，1995，9（1）：13-18.

[7] 方綱清，阮伏水，吳雄海，等.福建省主要土壤可蝕性特征初探[J].福建水土保持，1997（2）：19-23.

[8] Wen Zhongming, BrianG Lees, Jiao Feng, et al. Stratified Vegetation Cover Index： A New Way to Assess Vegetation Impact on Soil Erosion[J].Catena,2010,(83)：87-93.

[9] Saleh Arekhi, Yaghoub Niazi, Aman Mohammad Kalteh.Soil Erosion and Sediment Yield Modeling Using RS and GIS Techniques： A Case Study, Iran[J]. Arabian Journal of Geosciences,2012,5(2)：285-296.

[10] Li Hui, Chen Xiaoling, KyoungJae Lim, et al. Assessment of Soil Erosion and Sediment Yield in Liao Watershed,Jiangxi Province,China,Using USLE,GIS,and RS[J].Journal of Earth Science,2010,21(6)：941-953.

[11] 徐涵秋.南方典型紅壤水土流失區(qū)地表裸土動態(tài)變化分析——以福建省長汀縣為例[J]. 地理科學，2013，33（4）：489-496.

[12] 楊冉冉，徐涵秋，林娜，等.基于RUSLE的福建省長汀縣河田盆地區(qū)土壤侵蝕定量研究[J]. 生態(tài)學報，2013，33（10）：2974-2982.

[13] 林惠花，武國勝，朱鶴健，等.基于GIS和USLE的朱溪河小流域土壤侵蝕經(jīng)濟損失評估[J].江西農(nóng)業(yè)大學學報，2011，33（1）：183-188.

[14] Tang LY, Chen CC, Zou J, et al. OntoPlant：An integrated virtual plant software package for different scale application[C]//Proceeding of 2011 IEEE International Conference on Spatial Data Mining and Geographic Knowledge Services (ICSDM2011) in conjunction with Eighth Beijing International Workshop on Georaphical Informational Science (BJ_IWGIS 2011). Fuzhou： IEEE, 2011： 308-314.

[15] 江洪，王欽敏，汪小欽，等.區(qū)域性水土流失遙感動態(tài)監(jiān)測技術研究與實現(xiàn)[J].中國水土流失與生態(tài)安全綜合科學考察總結(jié)及水土保持發(fā)展戰(zhàn)略研討會論文集，2008：257-263.

[16] 武國勝，林惠花，朱鶴健.基于馬爾柯夫模型的福建長汀土壤侵蝕動態(tài)預測[J].福建師范大學學報：自然科學版，2011，27（1）：93-98，104.

[17] 林晨，周生路，吳紹華.基于優(yōu)化GeoCA 模型的土壤侵蝕時空變化模擬——以福建省長汀縣為例[J].地理科學進展，2012，31（7）：911-920.