基于微波遙感的裸露地表土壤鹽分含量的反演

摘要：基于Radarsat-2雷達影像數(shù)據對裸露地表土壤鹽分含量進行反演研究，利用數(shù)據模擬分析方法建立土壤鹽分含量的反演模型，并對土壤鹽分含量的反演值與實測值進行相關分析。結果表明：模擬值與實測值有很好的相關性（r2=0.904 6） ，可以利用Radarsat-2雷達數(shù)據進行大面積地表土壤鹽分的監(jiān)測與提取，能取得較好的結果。

關鍵詞：Radarsat-2；后向散射系數(shù)；反演；土壤鹽分

中圖分類號： S151.9+3文獻標志碼： A文章編號：1002-1302（2015）11-0442-03

收稿日期：2014-12-07

基金項目：國家自然科學基金（編號：51169016）。

作者簡介：弓永利（1976—），女，內蒙古呼和浩特人，副教授，高級工程師，主要從事工程測量、遙感、地理信息應用方面的科研和教學工作。E-mail：fllyzz@126.com。 土壤鹽漬化（soil salinization）是指隨著水分上升到地表并蒸發(fā)后，使土壤底層的鹽分或者是地下水的鹽分在表層土壤積累的過程，是土地退化和土壤荒漠化的主要表現(xiàn)形式。在干旱、半干旱地區(qū)，土地逐漸退化主要是由人類不合理的活動和脆弱的生態(tài)環(huán)境共同導致的，土壤鹽漬化已經成為現(xiàn)在世界上土地荒漠化和土地退化的主要類型之一[1-4]。鹽堿土不僅會引起土壤板結、肥力下降，加重土壤侵蝕，還會造成資源破壞，引發(fā)農業(yè)生產損失，從而嚴重威脅生物圈和生態(tài)環(huán)境安全[5]。

遙感及其相關技術在獲取地表信息方面具有宏觀性、實時性、動態(tài)性等特點，這使得土壤鹽漬化的監(jiān)測更加宏觀、快速、準確、有效，有利于對土地鹽堿化發(fā)展和演變機理的認識[4，6-7]。因此，對河套平原鹽漬化程度進行遙感監(jiān)測、反演可以作為該地區(qū)灌溉制度制定、鹽漬化土壤治理方案制定的基礎資料。

由于微波遙感技術具有全天時、全天候的工作能力，而且對云霧、雨雪和地物都有一定的穿透能力，可以更好地識別、分析地物，提取所需的信息。土壤鹽分含量變化趨勢與土壤水分含量有著密不可分的關系，再加上微波遙感對土壤水分十分敏感[8-10]，因此微波遙感技術已經成為當前監(jiān)測土壤鹽分的主要方法。

本研究基于Radarsat-2雷達影像數(shù)據的處理，對裸露地表土壤鹽分含量進行反演研究，以期利用SPSS軟件進行數(shù)據模擬分析，建立地表土壤鹽分含量后向散射模型，從而為大面積獲得地表土壤鹽分時空分布信息提供理論依據和技術支撐。

1材料與方法

1.1研究區(qū)概況

研究區(qū)位于內蒙古自治區(qū)巴彥淖爾市杭錦旗境內，地理坐標范圍為40°44′～41°01′ N、106°52′～107°13′ E，海拔1 032～1 050 m ，面積約為625 km2。該地區(qū)地貌主要由沖擊平原、洪積平原、河漫灘3種地形構成，地勢從西南到東北逐漸降低，屬于河套平原。氣候特點為溫帶大陸性氣候，干燥、晝夜溫差大，年均氣溫8.5 ℃，年均降水量138.2 mm左右，無霜期130 d 左右，日照時間較長，年均日照時間3 181 h。地面特征以耕地與鹽堿化裸地為主，“種植類型”主要有玉米、葵花、小麥等。

1.2雷達影像數(shù)據

采用的雷達數(shù)據為2013年3月13日的一景Radarsat-2標準四極化模式雷達數(shù)據，此時研究區(qū)地表狀態(tài)為春耕前期時的裸露地表。雷達影像的工作波段為C，中心波長5.7 cm，頻率5.3 GHz，極化方式為HH、HV、VH、VV，數(shù)據類型為 Single Look Comple（SLC）。

1.3樣品數(shù)據分析

本試驗在視野開闊并無遮擋物的裸露地表取樣，所采集土壤樣品深度限定在0～10 cm，且采樣時間與衛(wèi)星飛過試驗區(qū)時間同步。本試驗在試驗區(qū)內共均勻采集60個土壤樣品，同時使用GPS獲取樣點坐標位置信息。然后將所采集樣品分別放置鋁盒和塑封袋內，對鋁盒現(xiàn)場稱重，待到實驗室用EDTA滴定法獲取土壤含鹽量，烘干法獲得土壤含水量。

1.4雷達影像處理

與可見光、紅外等光學遙感相比，Radarsat-2雷達影像數(shù)據處理起來比較復雜。對于SLC類型的雷達影像，先對其進行斜距到地距的轉換后才能獲得真正的地距影像，然后利用ENVI 4.8軟件進行輻射定標、幾何校正、濾波去噪等處理后得到的影像，才能更好地進行后向散射系數(shù)的提取。

利用ARCGIS 9.3軟件把GPS 獲取采樣點的坐標位置制作成矢量文件，再通過ENVI 4.8軟件在雷達影像上顯示采樣點的位置信息，從而獲得樣點的后向散射系數(shù)。由于雷達影像數(shù)據存在嚴重的斑點和噪聲，既便是進行濾波處理后，影像還是有很多斑點，這樣就給樣點后向散射系數(shù)的提取造成很大的干擾。為了獲取更加準確的樣點后向散射系數(shù)，本研究將取樣點象元周圍的象元后向散射系數(shù)的平均值作為該中心樣點最后的后向散射系數(shù)。

2結果與分析

2.1模型的建立

SPSS是一種非常經典、強大的統(tǒng)計分析軟件，它有比較成熟的理論基礎，可以確定變量之間的定量關系，并可以進行相應的預測，且能反映出統(tǒng)計變量之間的統(tǒng)計變化規(guī)律。在獲取四極化地表后向散射系數(shù)和實驗室數(shù)據等基礎數(shù)據后，通過SPSS軟件“控制變量”的思想，逐步加入和剔除變量，觀察變量加入前后模型的變化，確定該變量對本研究結果是否有積極影響，來探究各變量之間的關系。最后利用SPSS軟件進行數(shù)據模擬回歸分析，建立回歸模型，并根據已有實測數(shù)據和圖像數(shù)據來估計模型的各個參數(shù)。

選取HH、HV、VH、VV、HH2、HV2、VH2、VV2以及含水量為自變量，采用向前和向后的建模方法，對變量進行加入和剔除，從而建立線性回歸模型。選取45個樣點數(shù)據進行軟件模型模擬，通過向前和向后方法的模型分析，最終建立以HH、VV、HH2、HV2、VH2、VV2以及含水量為自變量的土壤含鹽量反演公式。具體模型如下：

Y=1.362+1.128×a-0.48×b×10.557×c-0.981×d+ 0.747 ×e-2.906×f+1.925×g。 （1）

式中：Y為土壤含鹽量；a為考慮入射角影響的HH極化的后向散射系數(shù)；b為考慮入射角影響的VV極化的后向散射系數(shù)；c為考慮入射角影響的HH極化的后向散射系數(shù)的平方值除以100；d為考慮入射角影響的HV極化的后向散射系數(shù)的平方值除以100；e為考慮入射角影響的VH極化的后向散射系數(shù)的平方值除以100；f為考慮入射角影響的VV極化的后向散射系數(shù)的平方值除以100；g為土壤含水量。

2.2模型驗證

為了進一步驗證模型的正確性，使用公式（1）對其余的15個樣點進行驗證，通過軟件Envi4.8中Band Math 波段運算功能進行土壤鹽分含量的計算，結果如表1所示。同時，對反演值與實測值進行相關性分析，結果如圖1所示。

在表1中，26號、44號樣本的反演結果為負值，則表明反演出錯，其余反演結果誤差的絕對值最大為9號樣本，值為0.296 g/kg，最小的為5號樣本，值為0.012 g/kg。誤差的絕對值<0.1 g/kg的有8個，≥0.1 g/kg的有7個，≥0.2 g/kg的僅有3個。表明模型計算的土壤鹽分含量與實測值的誤差較小。利用該模型對所有野外觀測數(shù)據進行模擬，樣本實測值與反演值的統(tǒng)計結果如表3所示。

從表2可以看出，土壤鹽分含量的實測值與反演值之間最大值、最小值之差分別約為0.25、0.33 g/kg。為了研究不同鹽漬化程度的地表后向散射響應特征，根據采樣點表層的不同土壤鹽分含量分成4組，即非鹽漬地、輕度鹽漬地、中度鹽漬地和重度鹽漬地[11]，具體情況如表3所示。

利用ENVI4.8軟件和模型公式反演整個研究區(qū)域地表的土壤鹽分，最終獲得該研究區(qū)域土壤鹽分的空間分布圖，然后利用表4的分類方法進行統(tǒng)計，結果如圖2所示。

3結論

本研究基于Radarsat-2雷達影像數(shù)據和實測數(shù)據，利用SPSS軟件進行數(shù)據模擬分析，建立地表土壤鹽分含量的反演模型，反演獲得的土壤鹽分含量與實測值之間最大值、最小值、之差均為0.25、0.33g/kg，且相關系數(shù)達到了0.9046。

這表明此模型在不考慮其他更多地表因素情況下，可以獲得較好的反演結果。因此，Radarsat-2雷達微波遙感可以利用該模型進行大面積的裸露地表土壤鹽分含量監(jiān)測，并且有較好的應用前景。

參考文獻：

[1]劉興土. 松嫩平原退化土地整治與農業(yè)發(fā)展[M]. 北京：科學出版社，2001.

[2]李鳳全，吳樟榮. 半干旱地區(qū)土地鹽堿化預警研究——以吉林省西部土地鹽堿化預警為例[J]. 水土保持通報，2002，22（1）：57-59.

[3]韋建波. RADARSAT-2數(shù)據在干旱區(qū)鹽漬地信息提取中的應用[D]. 烏魯木齊：新疆大學，2009：1-67.

[4]董新光，周金龍，陳躍濱. 干旱內陸區(qū)水鹽監(jiān)測與模型研究及其應用[M]. 北京：科學出版社，2007：1-27.

[5]Dehaan R L，Taylor G R. Field-derived spectra of salinized soils and vegetation as indicators of irrigation-induced soil satinization[J]. Remote Sensing of Environment，2002（80）：406-418

[6]關元秀，劉高煥. 區(qū)域土壤鹽漬化遙感監(jiān)測研究綜述[J]. 遙感技術與應用，2001，1（1）：40-44.

[7]徐東瑞，高廣惠. 鹽漬化土壤遙感動態(tài)監(jiān)測及人工調控模式——以河北省滄州市為例 [J]. 華北農業(yè)學報，1995，10：175-180.

[8]余凡，趙英時. 合成孔徑雷達反演裸露地表土壤水分的新方法[J]. 武漢大學學報：信息科學版，2010，35（3）：317-321.

[9]陳晶，賈毅，余凡. 雙極化雷達反演裸露地表土壤水分[J]. 農業(yè)工程學報，2013，29（10）：109-115.

[10]鄭磊，王耀強，郭曉靜，等. 四極化Radarsat-2數(shù)據對裸露地表土壤水分的反演[J]. 江蘇農業(yè)科學，2014，42（5）：250-252.

[11]李和平. 新疆干旱土綱基層分類探討[J]. 干旱區(qū)研究，2001（2）：56-60.黃艷，徐祥明，冷雄，等. 贛南不同植被覆蓋下紅壤有機碳剖面分布特征[J]. 江蘇農業(yè)科學，2015，43（11：445-447.