孫浩然，趙志根，趙佳星，陳衛(wèi)衛(wèi)

(1.安徽理工大學 空間信息與測繪工程學院，安徽 淮南 232001；2.礦山采動災(zāi)害空天地協(xié)同監(jiān)測與預(yù)警安徽普通高校重點實驗室，安徽 淮南 232001；3.礦區(qū)環(huán)境與災(zāi)害協(xié)同監(jiān)測煤炭行業(yè)工程研究中心，安徽 淮南 232001)

0 引言

土壤的物理、化學、生物等許多屬性都直接或間接地與有機質(zhì)的存在有關(guān)。土壤有機質(zhì)是土壤的重要組成部分，是評價土壤肥力的核心指標之一。耕作土壤表層有機質(zhì)含量通常在5%以下，雖只占土壤總質(zhì)量的很小一部分，卻在土壤肥力、環(huán)境保護、植物生長和農(nóng)業(yè)生產(chǎn)等方面有著積極的作用和意義[1-3]。土壤有機質(zhì)在空間分布上具有異質(zhì)性且呈連續(xù)變化特點，掌握有機質(zhì)的空間格局分布有利于農(nóng)業(yè)的精細化管理。土壤有機質(zhì)含量影響土壤色調(diào)、土壤孔隙度等進而影響光譜反射特性[4]。隨著衛(wèi)星遙感、土壤近地傳感、數(shù)據(jù)挖掘技術(shù)和建模技術(shù)的發(fā)展，采用遙感技術(shù)實現(xiàn)表層土壤有機質(zhì)空間格局反演及碳庫的數(shù)字制圖成為農(nóng)業(yè)與環(huán)境遙感領(lǐng)域的熱點[5]。

以往土壤有機質(zhì)空間分布格局反演通常采用地統(tǒng)計方法，在野外規(guī)則化地采集大量樣本數(shù)據(jù)并在室內(nèi)分析，通過插值法實現(xiàn)尺度上推[6-7]。該方法取得了較好的效果，但在實際應(yīng)用中，大面積的規(guī)則采樣操作復(fù)雜且成本過高。隨著遙感技術(shù)在光譜分辨率、空間分辨率、時間分辨率等方面顯著提高使得遙感影像在反演表層土壤有機質(zhì)中得到廣泛應(yīng)用，遙感技術(shù)反演區(qū)域表層土壤有機質(zhì)含量已成為土壤調(diào)查、土地質(zhì)量評價的主要手段。文獻[1，8]利用多光譜Landsat遙感影像DN值與表層土壤有機質(zhì)的聯(lián)系，建立預(yù)測模型，對區(qū)域有機質(zhì)含量進行大尺度反演，并完成表層土壤有機質(zhì)空間分布格局制圖。Landsat遙感影像空間分辨率較低，在反演表層土壤有機質(zhì)時受混合像元影響；其光譜分辨率低，土壤有機質(zhì)含量的差異與反射率(或DN值)之間的響應(yīng)關(guān)系不強，難以精確反演區(qū)域表層土壤有機質(zhì)含量空間格局分布。高光譜精細的光譜曲線構(gòu)成了獨特的超多維光譜空間，通過識別地物光譜的診斷特征可以在表層土壤有機質(zhì)反演過程中取得更好的效果[9]。文獻[10-13]利用高光譜數(shù)據(jù)的光譜反射率及反射率的變換形式，采用線性回歸法、多元線性回歸法、神經(jīng)網(wǎng)絡(luò)法及最小二乘回歸法分析構(gòu)建土壤有機質(zhì)預(yù)測模型。此方法精度雖高，但多基于田間和實驗室操作層面，尺度較小，不利于區(qū)域尺度表層土壤有機質(zhì)反演。文獻[14-15]將小波分析應(yīng)用到土壤光譜數(shù)據(jù)的處理中，確定了原始特征光譜的最佳小波分解尺度，并在小波包分解基礎(chǔ)上建立土壤有機質(zhì)預(yù)測模型，提出了一種新的土壤有機質(zhì)反演方法。文獻[16-17]將小波包分析方法應(yīng)用于Hyperion遙感影像數(shù)據(jù)中，利用高光譜遙感影像反演表層土壤有機質(zhì)含量，取得良好的預(yù)測效果。高光譜遙感反演表層土壤有機質(zhì)往往受限于空間尺度小，難以完成區(qū)域尺度表層土壤有機質(zhì)反演。隨著中國自主發(fā)射的珠海一號(Zhuhai-1)衛(wèi)星高光譜影像投入使用，以其光譜分辨率和空間分辨率獨特的優(yōu)勢，能夠高精準、定量化地反演表層土壤有機質(zhì)。本文以珠海一號的高光譜影像為實驗遙感影像，通過幾何校正、輻射校正、降噪等手段處理遙感影像，探索高光譜影像反演鄉(xiāng)鎮(zhèn)尺度表層土壤有機質(zhì)含量的方法。引入小波包多尺度分析方法對高光譜遙感數(shù)據(jù)進行處理，分別建立多元線性回歸預(yù)測模型與基于高斯核的支持向量機回歸預(yù)測模型，對研究區(qū)表層土壤有機質(zhì)含量的預(yù)測取得了可靠的結(jié)果。結(jié)果表明，以小波包分解高光譜數(shù)據(jù)的低頻分量建立的基于高斯核變換的支持向量機預(yù)測模型，可大尺度、高精準地反演區(qū)域表層有機質(zhì)含量，為區(qū)域農(nóng)業(yè)施肥管理與精準農(nóng)業(yè)提供科學管理依據(jù)。

1 研究區(qū)概況及采樣點設(shè)計

實驗地區(qū)位于淮南市舜耕山以南的三和鎮(zhèn)、曹庵鎮(zhèn)以及大通區(qū)的部分農(nóng)村(32°36′55.04″N～32°26′22.29″N，117°4′19.81″E～116°55′49.59″E)。該區(qū)地處平原地帶，日照充足，全年降水量893.4 mm，是淮南市重要的農(nóng)業(yè)生產(chǎn)基地。研究區(qū)內(nèi)農(nóng)耕地主要以種植小麥、水稻、玉米、大豆、油菜等為主，以瓜果為輔，種植方式主要為小麥與水稻、玉米、大豆間作和輪作以及水稻的連作等。2019年5月初，通過手持式全球定位系統(tǒng)在研究區(qū)內(nèi)隨機但均勻地采集無農(nóng)作物覆蓋區(qū)域土樣34個樣點，在道路交叉口以及顯著地標處采集8個控制點以用于幾何校正。樣點均分布無植被覆蓋或幼苗期農(nóng)田中，在30 m×30 m內(nèi)采用五點取樣法獲得1個樣本，采樣深度為0～10 cm。將樣點烘干，過100目篩，送往安徽理工大學分析測試中心，采用元素分析儀法測定土壤中含碳量，通過SOM=SOC×1.724，將測定的有機碳含量轉(zhuǎn)化為有機質(zhì)含量。

2 數(shù)據(jù)預(yù)處理

2.1 遙感影像及其預(yù)處理

影像數(shù)據(jù)源為歐比特公司提供的高光譜遙感影像，成像時間為2019年4月5日。高光譜影像獲取較難，采樣時間與成像時間相差數(shù)日，土壤有機質(zhì)變化幾乎可以忽略不計。采集土樣為研究區(qū)內(nèi)已收割完畢和幼苗期農(nóng)田土樣，可以忽略植被對采樣點光譜信息的影響。珠海一號高光譜影像結(jié)合成像技術(shù)與光譜技術(shù)，獲取高分辨率的連續(xù)、寬波段的圖像數(shù)據(jù)，單顆高光譜衛(wèi)星的最大成像范圍150 km×2 500 km，彌補了以往高光譜影像光譜分辨率高而空間分辨率低和幅寬窄的缺點。同時，珠海一號衛(wèi)星具有良好的時間周期性，可在較短的時間內(nèi)實現(xiàn)同一地區(qū)的多遍重訪掃描；其波長范圍為400～1 000 nm，光譜分辨率可達2.5 nm，精細的光譜曲線可以檢測到光譜信息的細微差異，有利于地物特征的高精度反演。集成高光譜分辨、高空間分辨率及高時間分辨率等特點的國產(chǎn)高光譜影像在真正意義上促進了高光譜遙感的廣泛應(yīng)用。本研究所采用的高光譜影像具有32個波段，藍、綠、紅、近紅外波段對應(yīng)的中心波長分別是480(band2)、566(band7)、670(band14)和880 nm(band28)，空間分辨率為10 m，其空間尺度適用于城鎮(zhèn)區(qū)域表層土壤有機質(zhì)空間分布格局的反演與制圖。

為了消除輻射量誤差對表層土壤有機質(zhì)反演過程中的影響，對高光譜遙感圖像進行輻射校正，獲得高光譜影像大氣上層表觀反射率產(chǎn)品。為精確提取樣本采集點的光譜信息，須首先對遙感圖像進行幾何校正，糾正遙感圖像的位置信息誤差，使得采樣點位置準確定位于遙感影像上。利用圖像自帶的RPC文件對高光譜影像進行幾何校正，誤差控制在 0.5個像元。利用研究區(qū)矢量數(shù)據(jù)裁剪研究區(qū)的范圍，并利用5×5高斯低通濾波對裁剪后的遙感影像做圖像增強處理，抑制高頻信息，使圖像平滑，消除遙感影像在成像過程中產(chǎn)生的噪聲。Zhuhai-1高光譜遙感波段信息及對應(yīng)中心波長如表1所示。

表1 高光譜遙感影像波段信息

2.2 小波包分析

由于高光譜遙感影像的光譜反射率受大氣和地表多重誤差影響，光譜反射率與平原地區(qū)表層土壤有機質(zhì)之間的微弱關(guān)系不足以準確描述二者之間的預(yù)測模型。小波包分析能夠讓信息能量集中，在細節(jié)中尋找有序性，篩選出其中的規(guī)律，為信號提供一種更加精細的分析方法。根據(jù)文獻[14,16]將小波包分析應(yīng)用在土壤有機質(zhì)中的研究，本文選用Db5小波基函數(shù)對采樣點原始光譜進行三層小波包分解。采用小波包分析多尺度分解采樣點光譜數(shù)據(jù)，將原始光譜數(shù)據(jù)分解為高頻分量和低頻分量。高頻系數(shù)能夠反映原始光譜微小的吸收特征，低頻系數(shù)與土壤質(zhì)地、土壤微量元素含量、土壤含水量等因素對光譜反射率的影響有關(guān)。利用高頻系數(shù)、低頻系數(shù)分別進行光譜重構(gòu)能夠達到突出微小細節(jié)特征、增強近似特征信息的效果。小波包分解結(jié)果如圖1所示。

圖1 小波包多尺度分析結(jié)果

2.3 相關(guān)性分析

在ENVI中提取采樣點在遙感影像上的光譜反射曲線，將實測表層土壤有機質(zhì)值和提取的高光譜32個波段光譜反射率值、光譜經(jīng)過小波分解的低頻分量、小波分解的高頻分量導入SPSS 19.0中。對以上數(shù)據(jù)進行描述統(tǒng)計分析、相關(guān)分析和Pearson 相關(guān)分析，分析土壤有機質(zhì)值與高光譜影像各個波段反射率之間的相關(guān)關(guān)系。

3 回歸分析與精度評估

3.1 多元線性回歸模型

利用相關(guān)性分析結(jié)果中原始光譜、小波包分解低頻分量、小波包分解高頻分量與土壤有機質(zhì)含量相關(guān)性較高的波段作為自變量，有機質(zhì)含量作為因變量，建立多元線性回歸分析模型，探究多元線性回歸分析模型對土壤有機質(zhì)含量的預(yù)測效果。

3.2 基于高斯核函數(shù)的支持向量機回歸模型

利用支持向量機回歸(support vector machine regression，SVR)建立數(shù)據(jù)模型時，基于高斯核函數(shù)實現(xiàn)一個非線性變換，將輸入空間(歐式空間Rn或離散幾何)對應(yīng)于一個特征空間(希爾伯特空間)，使得在輸入空間Rn中的超曲面模型對應(yīng)于特征空間中的超平面模型[18]。通過在特征空間中求解線性支持向量機建立擬合模型，使數(shù)據(jù)點到該平面的距離最近。

3.3 精度評估

對于多元線性回歸分析，采用決定系數(shù)R2(式(1))對土壤有機質(zhì)預(yù)測模型的準確性進行評價。基于高斯核函數(shù)的支持向量機回歸模型采用預(yù)測數(shù)據(jù)的均方誤差MSE(式(2))和平均絕對誤差MAE(式(3))對其預(yù)測準確性進行評價。

(1)

(2)

(3)

4 結(jié)果與分析

4.1 相關(guān)性分析結(jié)果

通過對研究區(qū)土壤采樣點有機質(zhì)含量與高光譜影像各波段光譜反射率值、光譜經(jīng)過小波分解的低頻分量、小波分解的高頻分量的相關(guān)性分析，土壤光譜反射率與實測有機質(zhì)含量之間的相關(guān)性如圖2所示。

圖2 土壤光譜反射率與有機質(zhì)含量之間的相關(guān)系數(shù)

土壤有機質(zhì)含量與高光譜影像多個波段程負相關(guān)，在656 nm處負相關(guān)性最強(r=-0.680)。經(jīng)小波包分解后的低頻分量與高頻分量最大相關(guān)波段的相關(guān)性均有所增強，低頻分量與土壤有機質(zhì)的相關(guān)性和原始光譜與土壤的相關(guān)性較為一致，說明小波包低頻系數(shù)分解能夠有效突出光譜近似特征。不同光譜形式與有機質(zhì)含量之間的最大相關(guān)性如表2所示。

表2 不同光譜變換與有機質(zhì)含量最大相關(guān)波段分析結(jié)果

4.2 土壤有機質(zhì)模型預(yù)測結(jié)果

對原始光譜、小波分解低頻分量、小波分解高頻分量與土壤有機質(zhì)含量相關(guān)波段做多次多元線性回歸實驗，得到結(jié)果如圖3所示。

圖3 不同光譜形式有機質(zhì)含量預(yù)測模型

采用珠海一號高光譜遙感影像反演土壤有機質(zhì)含量時利用多元線性回歸分析方法預(yù)測效果良好，R2均大于0.74；采用小波包多尺度分解的高頻分量和低頻分量建立土壤有機質(zhì)預(yù)測模型，效果優(yōu)于原始光譜預(yù)測效果。表3為不同光譜形式的土壤有機質(zhì)含量預(yù)測方程。

表3 不同光譜變換與有機質(zhì)含量最大相關(guān)波段分析結(jié)果

由多元回歸分析結(jié)果可以發(fā)現(xiàn)，采用高光譜影像原始光譜反射率需要多個波段進行反演土壤表層有機質(zhì)。光譜經(jīng)小波包分解的低頻分量和高頻分量實現(xiàn)了信息的集中，能夠以較少的波段達到較好的預(yù)測效果，且根據(jù)小波包分解后的高頻分量相關(guān)波段與低頻分量相關(guān)波段建立回歸模型，R2分別為0.770、0.789。結(jié)果表明，根據(jù)小波包分解得到土壤光譜信息的低頻分量與高頻分量所建立的土壤有機質(zhì)模型可取得更佳的效果。

在采用遙感的方法反演表層土壤有機質(zhì)過程中，遙感影像土壤反射率與經(jīng)過精細處理的土壤所測定的土壤有機質(zhì)含量存在誤差。同時，土壤含水量與土壤中所含的微量元素也會對遙感影像反射率曲線產(chǎn)生影響，土壤有機質(zhì)含量與遙感影像的相關(guān)性受多種因素干擾，遙感影像反射率與土壤有機質(zhì)含量之間并非線性關(guān)系。簡單的多元線性回歸模型沒有考慮到在反演土壤表層有機質(zhì)過程中的多種因素干擾及反射率與土壤有機質(zhì)含量之間并非線性關(guān)系，難以對研究區(qū)內(nèi)海量土壤光譜信息數(shù)據(jù)進行更佳的預(yù)測。采用遙感方法進行土壤有機質(zhì)數(shù)字制圖的過程中，基于高斯核函數(shù)變換的支撐向量機回歸模型處理遙感大數(shù)據(jù)往往更具有優(yōu)勢?；诟咚购撕瘮?shù)變換的支撐向量機回歸模型擅于解決非線性問題，通過對已有數(shù)據(jù)的訓練，尋找出一條最佳的預(yù)測曲線，實現(xiàn)對土壤有機質(zhì)含量的預(yù)測。根據(jù)小波包分解得到的低頻分量和高頻分量，建立支持向量機回歸表層土壤有機質(zhì)預(yù)測模型，預(yù)測效果比多元線性回歸模型更加精準。不同預(yù)測模型精度評估對比如表4所示。

表4 預(yù)測模型精度評估對比

由于原始光譜反射率詮釋表層土壤有機質(zhì)含量的敏感信息不集中，且本次實驗采樣點數(shù)據(jù)較少，根據(jù)原始光譜反射率建立的基于高斯核變換的支持向量機預(yù)測模型極容易出現(xiàn)過擬合現(xiàn)象，不具備泛化能力。從表4可以得出，根據(jù)小波包分解光譜后的低頻分量和高頻分量建立高斯核變換的支持向量機回歸模型其預(yù)測值均方誤差(MSE分別為2.591、5.428)均小于多元線性回歸模型預(yù)測值均方誤差(MSE分別為9.302、8.765)，該方法預(yù)測能力優(yōu)于多元回歸模型；根據(jù)小波包分解光譜后的低頻分量和高頻分量建立高斯核變換的支持向量機回歸模型其預(yù)測值平均絕對誤差(MAE)分別為1.426、1.520，預(yù)測值相比于多元線性回歸模型的預(yù)測值更加接近真值。研究結(jié)果表明，在小波包對土壤光譜分解結(jié)果的基礎(chǔ)上建立高斯核變換的支持向量機模型預(yù)測效果優(yōu)于多元線性回歸模型，其中根據(jù)小波包分解的低頻分量建立的基于高斯核函數(shù)變換的支持向量機回歸模型在表層土壤有機質(zhì)的預(yù)測中效果最好，為區(qū)域表層土壤有機質(zhì)的精準反演提供了一種新方法。

4.3 結(jié)果分析

研究區(qū)表層土壤有機質(zhì)含量與高光譜影像原始光譜反射率最大相關(guān)波段為656 nm處，國內(nèi)學者何挺等[10]通過對174個土樣與ASD光譜反射率進行分析，得出反射率與有機質(zhì)含量在675 nm 附近的相關(guān)性最強。考慮到珠海一號高光譜影像每個波段間隔接近15 nm，且不同區(qū)域土壤有機質(zhì)含量的光譜響應(yīng)波段略有差異，本研究中土壤有機質(zhì)含量與光譜反射率的相關(guān)性與何挺等人的研究結(jié)果近似。本文所提出的預(yù)測模型精度與文獻[16-17]略有不同，原因可能是由于采樣樣本不足和不同研究區(qū)土壤含量差異性造成的。東北地區(qū)土壤有機質(zhì)含量較高，土壤有機質(zhì)預(yù)測模型精度普遍較高；Hyperion數(shù)據(jù)光譜分辨高，光譜范圍大，土壤有機質(zhì)含量的敏感響應(yīng)波段相關(guān)性大，往往在預(yù)測土壤有機質(zhì)含量時會取得很好的效果，但其空間尺度較差，不足以廣泛應(yīng)用。珠海一號高光譜遙感具有較好的空間尺度，采用珠海一號高光譜遙感影像反演表層土壤有機質(zhì)分布具有精度和空間尺度上的雙重優(yōu)勢，有利于推進遙感大尺度高精準反演的應(yīng)用；其節(jié)約成本、較好的時間周期性等優(yōu)點，有利于表層土壤有機質(zhì)含量空間分布制圖。

研究區(qū)以農(nóng)業(yè)為主，表層土壤有機質(zhì)含量范圍主要位于8～38 g·kg-1之間。農(nóng)田類型主要為平原地帶耕作農(nóng)田、開發(fā)區(qū)農(nóng)田、舜耕山腳少量丘陵農(nóng)田。實驗發(fā)現(xiàn)，研究區(qū)內(nèi)平原地帶土壤有機質(zhì)含量普遍較高，介于25～38 g·kg-1之間，少數(shù)土壤有機質(zhì)含量過高得益于農(nóng)民的精細管理及部分水源充足、土地肥沃的田塊。開發(fā)區(qū)農(nóng)田原屬于耕作農(nóng)田，由于開發(fā)建設(shè)導致開發(fā)區(qū)內(nèi)農(nóng)田存在土地荒蕪現(xiàn)象，靠近城市開發(fā)區(qū)的一些農(nóng)田與部分耕作農(nóng)田土壤有機質(zhì)含量介于15～25 g·kg-1之間，土壤有機質(zhì)含量有所下降。土壤有機質(zhì)含量較低的地塊多為山腳坡地，該部分農(nóng)田由于土壤母質(zhì)原因?qū)е峦寥婪柿ω汃ぃ又車_發(fā)區(qū)建設(shè)，農(nóng)田種植作物復(fù)雜，低矮經(jīng)濟樹木、糧食作物、蔬菜、瓜果等交錯種植且少于管理，土壤有機質(zhì)含量介于9～15 g·kg-1之間。合理地劃分研究區(qū)土壤有機質(zhì)含量等級，利用珠海一號高光譜遙感影像，采用小波包分解光譜的低頻分量建立基于高斯核變換的支持向量機回歸模型，可實現(xiàn)區(qū)域表層土壤的精準反演，為區(qū)域精準農(nóng)業(yè)提供科學決策。

5 結(jié)束語

本文探索采用珠海一號高光譜影像反演區(qū)域表層土壤有機質(zhì)含量，通過小波包分析的方法，可有效分離原始光譜的微小吸收特征和近似特征，采用低頻分量、高頻分量建立的多元線性回歸分析預(yù)測土壤有機質(zhì)含量效果優(yōu)于原始光譜預(yù)測效果，提高了土壤有機質(zhì)含量的預(yù)測精度。依據(jù)原始光譜，通過小波包分析分解出的低頻分量、高頻分量，建立基于高斯核變換的支持向量機回歸模型效果優(yōu)于多元線性回歸模型，其中以小波包分解的低頻分量建立的基于高斯核變換的支持向量機回歸模型預(yù)測效果最好。由于珠海一號高光譜遙感影像具有良好的空間尺度、光譜分辨率、時間周期等優(yōu)點，研究結(jié)果為區(qū)域土壤有機質(zhì)反演制圖提供了新方法、新思路，促進了遙感大尺度高精準反演區(qū)域表層土壤有機質(zhì)的應(yīng)用。