羅亞東,許 齊,郭宇宏,2,孫智鵬,金煜龍,陳俊英,2,余衛(wèi)華,2
(1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院,陜西 楊凌 712100;2.西北農(nóng)林科技大學(xué) 旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室,陜西 楊凌 712100)
土壤含水率是分析冬小麥生育情況的主要依據(jù),對于理解和解決農(nóng)業(yè)生產(chǎn)、生態(tài)規(guī)劃以及水資源管理中的科學(xué)與實際問題至關(guān)重要[1]。傳統(tǒng)的土壤含水率監(jiān)測方法操作復(fù)雜、且數(shù)據(jù)滯后。遙感技術(shù)因其快速、簡便以及無損的特點,已廣泛應(yīng)用于農(nóng)田農(nóng)情的監(jiān)測[2-7]?;谶b感的墑情監(jiān)測是農(nóng)村水利科技工作的重點領(lǐng)域[8]。但是在無人機多光譜遙感實踐中,存在尺度效應(yīng)影響監(jiān)測精度的問題。
近年來,遙感監(jiān)測已廣泛應(yīng)用于土壤含水率的研究,楊珺博[9]等人利用無人機遙感監(jiān)測數(shù)據(jù)反演土壤含水率,發(fā)現(xiàn)反演精度可達80%以上,且反演效果最好的土壤深度為10~20 cm。Cheng Minghan[10]等利用無人機多模態(tài)數(shù)據(jù)對玉米田間含水率進行4種模型建立與反演,得出多光譜數(shù)據(jù)在不同模型下均提供最精確的土壤水分估算。
尺度效應(yīng)是指在無人機遙感觀測中,隨著遙感分辨率的變化,遙感反演得到的數(shù)據(jù)與實際表現(xiàn)不一致的現(xiàn)象,作為遙感科學(xué)的核心問題,近年來受到了學(xué)者的廣泛關(guān)注[11,12]。歐盛華[13]等利用機載航空WIDAS 數(shù)據(jù),發(fā)現(xiàn)尺度效應(yīng)是空間異質(zhì)性和反演函數(shù)非線性的綜合體現(xiàn),空間異質(zhì)性和反演函數(shù)非線性程度的增加均會造成尺度效應(yīng)的增大。徐凱健[14]通過對多光譜影像的森林樹種識別及其空間尺度響應(yīng)的研究。結(jié)果表明,影像空間分辨率對區(qū)域樹種識別結(jié)果具有顯著影響。Zhang[15]等利用偏振遙感器,根據(jù)地物目標(biāo)的空間異質(zhì)性和成像分辨率帶來的限制,發(fā)現(xiàn)了尺度效應(yīng)對地物偏振圖像的影響。Tang[16]等對SEBAL 模型的空間尺度效應(yīng)進行探究,得出高低分辨率對觀測熱通量的影響??梢?,尺度效應(yīng)對于無人機遙感觀測有一定影響,但目前尺度效應(yīng)對于無人機遙感監(jiān)測土壤水分的研究還鮮有報道,在實踐中缺乏最佳的監(jiān)測尺度。在提取作物冠層信息過程中,土壤、陰影背景均會對冠層信息準(zhǔn)確性造成影響[17,18]。楊帥[19]等基于RGRI指數(shù)法獲取的玉米冠層溫度與實測冠層溫度的相關(guān)性,得出不剔除土壤背景效果最差。王佳兒[20]等通過剔除土壤背景反演玉米田間含水率,發(fā)現(xiàn)數(shù)據(jù)的相關(guān)性和回歸模型的精度始終比未剔除土壤背景數(shù)據(jù)的相關(guān)性和模型精度差。紋理特征在冠層信息提取中中是一種相當(dāng)重要的特征,它反映了圖像或物體本身的屬性[21]。王思宇[22]利用紋理特征為基礎(chǔ),基于多光譜影像并構(gòu)建玉米冠層FAPAR 估算模型,發(fā)現(xiàn)其精度較高。目前仍缺乏通過去土壤去陰影掩膜方法提取作物冠層信息。
本文通過無人機搭載六波段多光譜相機在不同高度下獲取冬小麥紋理特征,剔除土壤背景、陰影背景,在此基礎(chǔ)上研究土壤含水率和植被紋理特征之間的關(guān)系。探究在不同分辨率下無人機遙感數(shù)據(jù)與土壤含水率的相關(guān)性并分析尺度效應(yīng)對于監(jiān)測數(shù)據(jù)的影響。最后得出尺度效應(yīng)對無人機多光譜遙感監(jiān)測土壤含水率的影響,獲取無人機監(jiān)測土壤水分的最佳飛行高度。
研究區(qū)域位于陜西省楊陵區(qū)陜西楊凌西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室的灌溉試驗站(34°17'N,108°04'E,海拔521 m)。該地區(qū)為典型的溫帶大陸性季風(fēng)季候,夏季炎熱,冬季寒冷,蒸發(fā)強度較大。年平均氣溫在12 ℃左右,多年平均降水量為630 mm,年平均蒸發(fā)量為884 mm,土壤干容重為1.38 g/cm3。
本試驗冬小麥選用小偃22 號,于2020年10月20日進行人工播種,共播種7.68 kg種子,并且施用20 kg的復(fù)合肥料作為底肥。同時在其生長過程中嚴(yán)格控制水分條件,試驗田灌水上限分別采用田間持水量的50%、65%、80%和95%,設(shè)置4 個水分處理,3 個重復(fù),共12 個小區(qū),以獲得不同田持下的土壤含水率數(shù)據(jù),使數(shù)據(jù)更有代表性。每個小區(qū)的面積為4 m×4 m=16 m2,每個小區(qū)間隔150 cm,相鄰兩個小區(qū)之間設(shè)置保護行。小區(qū)內(nèi)種植的冬小麥行距為25 cm,每一行播種40 g種子。為減少外界環(huán)境的干擾,每個小區(qū)都配有擋雨棚,圖1為試驗小區(qū)布置。
圖1 實驗小區(qū)布置Fig.1 Layout of experimental plot
1.3.1 土壤水分獲取
分別在3 個生育期在在每個實驗小區(qū)取3 個采樣點A、B、C,在采樣點處用土鉆采集10、20、30、40、60 cm 深度的土樣,每個土樣30 g 左右,土樣采出后立即裝入鋁盒中。將采集完成的鋁盒放入干燥箱,在105 ℃干燥后稱重并計算土壤含水率。烘干法測定土壤質(zhì)量含水率的計算式為:
式中:W1為鋁盒加濕土壤質(zhì)量;W2為鋁盒加干土質(zhì)量;W3為空鋁盒質(zhì)量。
1.3.2 無人機平臺與冠層數(shù)據(jù)獲取
在試驗田內(nèi)選取晴朗天氣下的土壤樣本,以配合無人機多光譜遙感圖像的獲取。采集時間為冬小麥分蘗期(2020年12月24日)、拔節(jié)期(2021年3月23日)、抽穗期(2021年4月19日),每日11∶00-14∶00。
本試驗采用由深圳市大疆創(chuàng)新科技有限公司生產(chǎn)的經(jīng)緯Matrice 600 六旋翼無人機,搭載Micro-MCA 多光譜相機進行多光譜遙感圖像采集。Micro-MCA 相機共有6個鏡頭,分別對應(yīng)6 個波段,波長分別為900 nm(近紅外)、800 nm(近紅外)、720 nm(紅邊)、680 nm(紅光)、550 nm(綠光)、490 nm(藍光)。無人機飛行高度設(shè)置為19、37、55、74 和92 m,對應(yīng)的分辨率分別為10、20、30、40 和50 mm,鏡頭垂直向下采集光譜信息。拍攝之前在實驗區(qū)附近放置白板和黑板,用于研究區(qū)內(nèi)進行標(biāo)定,從而獲取冬小麥光譜反射率。
使用多光譜相機完成圖像采集后,使用PixelWrench2[23]軟件對分別對分蘗期、拔節(jié)期、抽穗期無人機多光譜遙感原始圖像進行預(yù)處理,將預(yù)處理后遙感圖像導(dǎo)入Pix4Dmapper 軟件,進行幾何校正、輻射定標(biāo)、試驗區(qū)拼接處理,得到試驗區(qū)多光譜遙感拼接圖像。
將拼接處理后的圖像導(dǎo)入ENVI5.3(64 bit)軟件中,剔除土壤和陰影背景。使用ENVI5.3(64 bit)軟件中Co-occurrence measure 工具,獲取多光譜圖像6 個波段基于二階概率統(tǒng)計濾波的均值、方差、協(xié)同性、對比度、相異性、信息熵、二階矩、相關(guān)性等八種紋理特征。即從每個多光譜圖像中提取48組紋理特征。
不同冠層光譜反射率和紋理特征之間差異較大,存在較強的共線性問題[23]。從理論上來[24]說,預(yù)測指標(biāo)與參考變量之間的相關(guān)性越強,其預(yù)測結(jié)果越可靠[25]。采用全變量建模時,會導(dǎo)致模型預(yù)測不準(zhǔn)確。
基于灰色系統(tǒng)理論的灰色關(guān)聯(lián)分析是一種利用比較序列和參考序列之間的相似度來判別各變量相關(guān)性強弱的方法,其模型準(zhǔn)確度高、數(shù)據(jù)流失性小[26,27]。
本文使用灰色關(guān)聯(lián)法研究步驟如下:
構(gòu)建參考序列:
構(gòu)建比較序列:
計算關(guān)聯(lián)系數(shù):
計算關(guān)聯(lián)度:
式中:y表示土壤含水率;xa、xb表示一個生育期下的紋理特征;ξi代表關(guān)聯(lián)系數(shù);si代表關(guān)聯(lián)度;n表示3 個生育期下12塊小區(qū)各自對應(yīng)的土壤含水率,n=36;k為6 個波段下的的紋理特征總數(shù),k=48;n為樣本組數(shù),n=36;i表示3 組掩膜處理,i=3。
本研究每組樣本經(jīng)灰色關(guān)聯(lián)法優(yōu)選后為8個輸入變量,采用多元線性回歸(MLR)、前饋神經(jīng)網(wǎng)絡(luò)(BPNN)和隨機森林(RF)3種回歸方法在5 個高度和3 種掩膜方法下分別建立土壤水分反演模型,共45個土壤水分模型。將36組樣本按2∶1的比例隨機取樣作為建模集和驗證集。即選取24 個樣本用于建模,其余12個用于驗證。
MLR 利用線性模型來擬合多個自變量和因變量的關(guān)系,從而確定多元線性回歸模型的參數(shù),回歸至原假設(shè)方程中,通過回歸方程來預(yù)測因變量的趨勢[28-30]。
BPNN 作為一種根據(jù)誤差逆?zhèn)鞑ニ惴ㄓ?xùn)練的多次前饋網(wǎng)絡(luò),它通過對人類神經(jīng)元的功能進行模擬,可以儲存及學(xué)習(xí)大量的輸入數(shù)據(jù)和輸出數(shù)據(jù),且不需對變量的映射關(guān)系進行描述,利用輸入和輸出數(shù)據(jù)建模,其對非線性系統(tǒng)具有很強的模擬能力[31,32]。本文搭建4 層BP 神經(jīng)網(wǎng)絡(luò),采用兩個隱含層,傳遞函數(shù)分別為雙曲正切S型函數(shù)和線性函數(shù),節(jié)點數(shù)分別為8和2。
RF 是一種基于分類樹的算法,通過對大量分類樹的匯總提高了模型的預(yù)測精度,是取代神經(jīng)網(wǎng)絡(luò)等傳統(tǒng)機器學(xué)習(xí)方法的新的模型。當(dāng)林中的樹數(shù)變大時,林的泛化誤差會收斂到極限??梢泽w現(xiàn)變量間的交互作用,具有很高的預(yù)測準(zhǔn)確率,對異常值和噪聲具有很好的容忍度,且不容易出現(xiàn)過擬合[33,34]。
構(gòu)建反演模型,比較3種掩膜方法下各反演模型的擬合優(yōu)度R2和均方根誤差RMSE,分析尺度效應(yīng)對于無人機遙感監(jiān)測土壤水分的影響。R2越接近1,RMSE越接近0,則反演模型的精度和準(zhǔn)確性越好[23]。
通過取土烘干法獲取的不同生育期下每個小區(qū)土壤含水率梯度值,結(jié)果如圖2所示。
從圖2可以看出,不同生育期的土壤含水率有所不同,各小區(qū)分別在分蘗期、拔節(jié)期和抽穗期的土壤含水率依次降低,土壤含水率的覆蓋范圍較大且具有代表性。
圖2 各生育期各小區(qū)土壤含水率Fig.2 SWC of each plot at each fertility stage
獲得3組不同掩膜方法下的紋理特征,每組數(shù)據(jù)統(tǒng)計各自48 個紋理特征。本文通過灰色關(guān)聯(lián)法結(jié)果顯示紋理特征與土壤含水率的相關(guān)度,進而對各分辨率與掩膜方法進行初步評價與分析。
使用灰色關(guān)聯(lián)法計算紋理特征與土壤含水率相關(guān)度,結(jié)果如圖3所示。
圖3 灰色關(guān)聯(lián)法處理結(jié)果Fig.3 Results of the GRA method
無掩膜(NM)處理中,10 mm 和30 mm 下紋理特征與土壤含水率相關(guān)性最顯著,20 mm和40 mm下紋理特征與土壤含水率相關(guān)性總體最佳,關(guān)聯(lián)度分布在0.7~0.8。去土壤掩膜(SRM)處理中,紋理特征與土壤含水率的相關(guān)性變化不大。去土壤去陰影掩膜(SSRM)處理中,總體相關(guān)性較差,分辨率為10 mm時紋理特征與土壤含水率相關(guān)性最顯著,分辨率為40 mm 時次之。
從圖3可以看出,在10 mm 下,紋理特征與土壤含水率呈現(xiàn)最高相關(guān)性,其次則為30 mm、20 mm、40 mm、50 mm。這是由于飛行高度提高,分辨率隨之降低,混合像元在總像元中的比例上升,對冬小麥紋理特征造成影響,紋理特征與土壤含水率的相關(guān)性下降。
在掩膜分析結(jié)果中,無掩膜與去土壤掩膜的結(jié)果差異較小,且優(yōu)于去土壤去陰影掩膜。本文認(rèn)為在進行去土壤掩膜處理時,混合像元容易分類,而陰影和冠層的差異較??;而在進行去土壤去陰影處理時,混合像元難以分類,導(dǎo)致過多地保留陰影或刪除必要的冠層像元導(dǎo)致紋理特征不能準(zhǔn)確反映冠層情況。
2.3.1 無掩膜處理下反演結(jié)果
圖4可以看出,多元線性回歸模型在分辨率為10 mm 時擬合優(yōu)度(R2)最高,其次是40 mm。BP 神經(jīng)網(wǎng)絡(luò)模型在分辨率為10 mm 時R2最高,30 mm 次之。隨機森林模型在分辨率為10 mm時R2較高,其次是30 mm。
圖4 無掩膜處理反演結(jié)果Fig.4 Regression results without mask
在無掩膜處理中,3 種回歸方法的最佳反演結(jié)果均在10 mm 時取得較好結(jié)果,擬合優(yōu)度在0.7 以上,其中BP 神經(jīng)網(wǎng)絡(luò)的擬合優(yōu)度為0.86,RMSE為0.011。值得注意的是,BP 神經(jīng)網(wǎng)絡(luò)和隨機森林兩種方法在30 mm 分辨率時R2超過了10 mm,因此在10 mm和30 mm分辨率下,冠層情況在紋理特征中得到準(zhǔn)確表達。
同時,3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出先降低再升高再降低的趨勢,在分辨率降低到30 mm 至40 mm 時R2出現(xiàn)升高趨勢。在此掩膜處理下,多元線性回歸模型與BP 神經(jīng)網(wǎng)絡(luò)的擬合優(yōu)度低于隨機森林模型。
2.3.2 去土壤掩膜處理下反演結(jié)果
將篩選后去土壤掩膜處理的36 組濾波信息分別使用多元線性回歸、BP 神經(jīng)網(wǎng)絡(luò)、隨機森林進行反演處理,結(jié)果如圖5所示。
由圖5得出,多元線性回歸模型在分辨率為40 mm時R2最高,其次是30 mm。BP 神經(jīng)網(wǎng)絡(luò)模型在分辨率為30 mm 時R2最高,其次是40 mm。隨機森林模型在分辨率為10 mm,40 mm次之。
圖5 去土壤掩膜處理反演結(jié)果Fig.5 Regression results with soil background removed by masking
在去土壤掩膜處理中,3種回歸方法的最佳反演結(jié)果不盡相同,總體分布顯示在該掩膜處理下,40 mm分辨率下冠層情況在紋理特征中得到準(zhǔn)確表達。
3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出先降低再升高再降低的趨勢,除10 mm 分辨率的結(jié)果較好之外,30 mm 與40 mm分辨率下的反演結(jié)果擬合優(yōu)度較高。在去土壤掩膜處理之下,3種模型的擬合優(yōu)度均在0.6~0.7之間。
2.3.3 去土壤去陰影掩膜處理下反演結(jié)果
將篩選后去土壤去陰影掩膜處理的36 組濾波信息分別使用多元線性回歸、BP 神經(jīng)網(wǎng)絡(luò)、隨機森林進行反演處理,結(jié)果如圖6所示。
圖6 去土壤去陰影掩膜處理反演結(jié)果Fig.6 Regression results with soil and shadow background removed by masking
圖6可以看出,在去土壤去陰影掩膜處理中,多元線性回歸模型在分辨率為10 mm 時R2最高,其次是20 mm。BP 神經(jīng)網(wǎng)絡(luò)模型在分辨率為10 mm時R2最高,其次是40 mm。隨機森林模型在分辨率為10 mm,40 mm次之。
在此種掩膜處理中,3 種回歸方法的擬合結(jié)果顯示,在10 mm與40 mm分辨率下取得相對較高的反演精度,見圖7。
圖7 各掩膜方法下擬合優(yōu)度對比Fig.7 Comparison of the goodness of fit under different masking method
3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出逐漸降低趨勢,值得注意的是,在去土壤去陰影掩膜處理之下,除10 mm 的擬合優(yōu)度高于0.7,其余分辨率在各自模型下的擬合優(yōu)度均在0.5左右。
分析反演結(jié)果可以得出在分辨率為10 mm和40 mm時,模型反演效果較好。本文認(rèn)為飛行高度較低時,遙感圖像分辨率較高,能較為精準(zhǔn)地表現(xiàn)冠層狀況,同樣使得冠層情況也能在紋理特征中得到準(zhǔn)確表達;分辨率為40 mm 時,反演效果僅次于最高分辨率10 mm,本文認(rèn)為該現(xiàn)象由于是混合像元在總像元中占比降低,當(dāng)分辨率由10 mm下降至30 mm時,混合像元在在總像元中的比重增加,在用閾值進行背景分割時,會保留部分混合像元中的背景,或?qū)⒉糠只旌舷裨械淖魑锾蕹?,?dǎo)致紋理特征不能準(zhǔn)確表達冠層信息。當(dāng)分辨率下降至40 mm 時,分辨率與葉面寬度相近,此時混合像元占比反而下降,從而獲得較為準(zhǔn)確的反演精度。由于遙感圖像中存在土壤、陰影和田間雜草等背景,在用紋理特征進行建模時作為無關(guān)變量干擾了模型精度。導(dǎo)致回歸模型的魯棒性較差。當(dāng)高度增加,分辨率降低至與作物葉面尺度相近時,混合像元中背景為優(yōu)勢類別的像元被完全剔除,冠層為優(yōu)勢類別的像元被保留,并且較少的背景像元與冠層像元混合后對紋理特征的影響降低,能一定程度上減小背景對回歸模型的干擾。
最后,在各掩膜處理下,3種建模方法的反演結(jié)果不盡相同。在無掩膜處理下,多元線性回歸的擬合精度最差;在去土壤掩膜處理下,隨機森林的擬合精度最差;而在去土壤陰影掩膜處理中,BP 神經(jīng)網(wǎng)絡(luò)的擬合精度最差。本文認(rèn)為原因如下:多元線性回歸對數(shù)據(jù)有比較嚴(yán)格的要求如:①隨機誤差項是一個期望值為0 的隨機變量;②解釋變量的所有觀測值,隨機誤差項有相同的方差;③隨機誤差項彼此不相關(guān)等。在無掩膜處理的紋理特征中,干擾值很多隨機誤差不可控,所以這些條件在無掩膜的情況下很難滿足。隨機森林在回歸分析中不能連順序輸出,也不能預(yù)測超出訓(xùn)練集范圍的數(shù)據(jù),導(dǎo)致在出現(xiàn)特定噪聲的數(shù)據(jù)集中建模時容易出現(xiàn)過擬合,在去土壤掩膜處理的紋理特征數(shù)據(jù)中,出現(xiàn)了陰影的噪聲數(shù)據(jù),對隨機森林回歸造成了干擾,并且在用支持向量機進行監(jiān)督分類時,部分土壤并沒有被掩膜處理,也對隨機森林回歸造成了干擾。而BP 神經(jīng)網(wǎng)絡(luò)存在局部極小化的問題,采取局部搜索的方法的容易陷入局部最小值, 并且BP神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)自由度很大,目前并沒有完整的理論指導(dǎo),因此本文采用的隱含層結(jié)構(gòu)可能較適合于有完整信息的數(shù)據(jù)集,而對噪聲的容忍能力較高,去土壤去陰影掩膜刪去過多有效信息,所以在去土壤去陰影掩膜處理的紋理特征中效果較差。
(1)本文通過對不同分辨率下獲得的紋理特征進行土壤含水率的反演,發(fā)現(xiàn)在拍攝分辨率為10、40 mm 時,反演模型的擬合優(yōu)度較好。但分辨率在40 mm 時對不同反演方法表現(xiàn)更穩(wěn)定。在實際應(yīng)用中,過高分辨率下工作效率低,不能大范圍進行遙感監(jiān)測,為工作帶來不便。故本文認(rèn)為無人機遙感監(jiān)測土壤水分的最佳分辨率為40 mm。
(2)30 mm-NM-RF 模型為改組數(shù)據(jù)反演冬小麥土壤含水率的最佳反演模型,但在30 mm 分辨率下其他模型反演效果不佳。
(3)對于3 種反演模型,表現(xiàn)最好的是NM-RF 模型,擬合優(yōu)度較高且在各分辨率下表現(xiàn)穩(wěn)定。綜合各掩膜處理方法表現(xiàn)最好反演模型是BPNN 模型,在大部分情況下表現(xiàn)穩(wěn)定有較好的反演效果,為機器學(xué)習(xí)在遙感監(jiān)測的應(yīng)用驗證了可行性。
(4)無掩膜處理和去土壤掩膜處理的結(jié)果優(yōu)于去土壤去陰影掩膜,且去土壤去陰影掩膜處理下分辨率較低時,反演效果較差。