基于無(wú)人機(jī)遙感影像的玉米冠層溫度提取及作物水分脅迫監(jiān)測(cè)

張智韜，于廣多，吳天奎，張譽(yù)馨，白旭乾，楊 帥，周永財(cái)

（1. 西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，楊凌 712100；2. 西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室，楊凌 712100）

0 引 言

及時(shí)準(zhǔn)確地獲取作物水分脅迫狀況對(duì)發(fā)展節(jié)水灌溉農(nóng)業(yè)、提高水分利用效率有重要意義[1]。冠層溫度由土壤—植物—大氣連續(xù)體內(nèi)熱量的傳輸與轉(zhuǎn)換決定[2]，冠層溫度及其相關(guān)指數(shù)可以指示作物水分脅迫狀況[3-4]，其中以作物水分脅迫指數(shù)（Crop Water Stress Index，CWSI）應(yīng)用最為廣泛。Idso 等[5]提出基于冠層溫度的CWSI 經(jīng)驗(yàn)公式，而后又有學(xué)者提出了CWSI 的理論公式與簡(jiǎn)化公式[6-8]，通過(guò)蒸散模型計(jì)算CWSI 的理論公式考慮了諸多氣象因素[9]，理論基礎(chǔ)較強(qiáng)，監(jiān)測(cè)精確度較高。近年來(lái)，國(guó)內(nèi)外低空無(wú)人機(jī)遙感技術(shù)快速發(fā)展，可借助熱紅外成像儀快速、大面積地采集冠層溫度信息，為監(jiān)測(cè)作物水分提供了一種可靠路徑[10-12]。然而當(dāng)作物未達(dá)到全覆蓋時(shí)，無(wú)人機(jī)獲取的試驗(yàn)區(qū)域圖像中包含冠層、土壤等多種地物信息，并且不同地物的溫度在一定范圍內(nèi)存在交集，這對(duì)提取冠層溫度有較大的干擾，因此剔除無(wú)人機(jī)影像中的土壤背景可以提高作物水分脅迫的監(jiān)測(cè)精度[13-14]，快速準(zhǔn)確地對(duì)熱紅外圖像進(jìn)行地物分類(lèi)并提取其中的冠層溫度是十分必要的。

按照是否借助可見(jiàn)光圖像可將冠層溫度提取方法分為兩類(lèi)：一種是依據(jù)冠層與土壤的溫度信息差異，直接在熱紅外圖像上分類(lèi)，如經(jīng)驗(yàn)公式法[15]、邊緣檢測(cè)法等[16-17]，統(tǒng)稱(chēng)為直接法。經(jīng)驗(yàn)公式法十分依賴所測(cè)定的大氣溫度的準(zhǔn)確性，誤差較大[18]；邊緣檢測(cè)法要求圖像具有較高的分辨率、不同地物溫度差異較大且地物類(lèi)別不宜過(guò)多[19]，這些因素在一定程度上限制了此方法的應(yīng)用。另一種是借助可見(jiàn)光圖像的光譜信息實(shí)現(xiàn)熱紅外圖像地物分離的間接法[20-22]。雖然高分辨率的可見(jiàn)光相機(jī)能夠很好地區(qū)分出復(fù)雜地物信息，但在實(shí)際操作中，受相機(jī)自身參數(shù)、拍攝角度、氣象因素以及人為因素等多個(gè)因素的干擾，兩圖像像元無(wú)法準(zhǔn)確無(wú)誤地一一對(duì)應(yīng)[23]，這直接影響冠層溫度的提取精度。直接法與間接法各有優(yōu)缺點(diǎn)，而現(xiàn)階段對(duì)兩類(lèi)方法的分類(lèi)效果進(jìn)行比較，并將兩類(lèi)方法相結(jié)合來(lái)協(xié)同提取冠層溫度的研究較少。

鑒于此，本研究以4 種不同水分處理的拔節(jié)期夏玉米為研究對(duì)象，利用無(wú)人機(jī)獲取試驗(yàn)區(qū)域每日13:00 的熱紅外和可見(jiàn)光圖像資料，分別采用 Otsu 算法、EXG-Kmeans 算法和Otsu-EXG-Kmeans 算法獲取冠層區(qū)域圖像，并從冠層輪廓圖、冠層區(qū)域提取精度、溫度直方圖以及與實(shí)測(cè)冠層溫度相關(guān)情況4 個(gè)方面對(duì)提取結(jié)果進(jìn)行精度評(píng)價(jià)，以確定提取熱紅外圖像冠層溫度的最佳方法，而后求得相應(yīng)CWSI，分析CWSI 同土壤含水率相關(guān)關(guān)系以及CWSI 日平均變化趨勢(shì)，對(duì)玉米水分虧缺狀況進(jìn)行定性分析，以期為無(wú)人機(jī)遙感精準(zhǔn)監(jiān)測(cè)作物水分脅迫狀況提供參考。

1 材料與方法

1.1 研究區(qū)域概況

研究區(qū)域位于陜西省楊凌示范區(qū)的中國(guó)旱區(qū)節(jié)水農(nóng)業(yè)研究院，34°17′N(xiāo)，108°04′E，海拔525 m。該區(qū)域?qū)儆跍貛Т箨懶约撅L(fēng)氣候，多年平均降水量約640 mm，主要集中在7－9 月，試驗(yàn)期間大氣平均溫度為33 ℃，總輻射較強(qiáng)。試驗(yàn)區(qū)域的土壤類(lèi)型為中壤土，凋萎系數(shù)8.6%，0～1 m 平均田間持水量為23%～26%（均為質(zhì)量含水率），平均土壤干容重1.44 g/cm3，地下水埋藏較深，其向上的補(bǔ)給量忽略不計(jì)。

1.2 試驗(yàn)設(shè)計(jì)

試驗(yàn)設(shè)置嚴(yán)重水分脅迫（T1）、中度水分脅迫（T2）、輕度水分脅迫（T3）以及充分灌溉（T4）4 種水分處理梯度，土壤含水率分別維持在田間持水量的50%、65%、80%以及95%～100%[17]，每個(gè)水分處理設(shè)有3 個(gè)重復(fù)，共計(jì)12 個(gè)小區(qū)（圖1），小區(qū)為4 m×4 m，2 個(gè)小區(qū)之間布置保護(hù)行。玉米品種為“華農(nóng)138”，播種時(shí)間為2020年6 月14 日，出苗時(shí)間為6 月20 日，播種深度約5 cm，種植行距60 cm，株距30 cm，行沿南北走向。試驗(yàn)區(qū)域設(shè)置小型氣象站，并配有移動(dòng)折疊雨棚以防止降雨干擾。

1.3 數(shù)據(jù)采集

拔節(jié)期玉米水分脅迫程度直接影響玉米長(zhǎng)勢(shì)，且無(wú)人機(jī)相機(jī)在晴朗無(wú)風(fēng)天氣拍攝效果較好?？紤]以上因素，本研究選擇在2020 年7 月26 日、7 月27 日、7 月29 日、7 月31 日、8 月2 日的13:00 進(jìn)行無(wú)人機(jī)圖像以及地面數(shù)據(jù)的采集。

1）熱紅外與可見(jiàn)光圖像采集

使用大疆公司研發(fā)的Matrice600Pro 六旋翼無(wú)人機(jī)，并搭載ZENMUSE XT 熱紅外相機(jī)與ZENMUSE X3 可見(jiàn)光相機(jī)進(jìn)行圖像采集。設(shè)置無(wú)人機(jī)飛行高度為20 m，鏡頭垂直地面進(jìn)行拍攝，熱紅外相機(jī)工作波段為7.5～13.5μm，像素為640×512，水平視場(chǎng)角32°，垂直視場(chǎng)角26°，鏡頭焦距19 mm，溫度分辨率0.05 ℃，20 m 圖像分辨率為1.8 cm?？梢?jiàn)光相機(jī)像素為4 000×3 000，視場(chǎng)角94°，鏡頭焦距20 mm，20 m 圖像分辨率為1.0 cm。

2）地面數(shù)據(jù)采集

在各小區(qū)內(nèi)選取可以代表小區(qū)平均長(zhǎng)勢(shì)且生長(zhǎng)良好的3 株玉米植株進(jìn)行標(biāo)記[24]，在圖像采集完成后，及時(shí)利用測(cè)溫槍獲取標(biāo)記植株承光面的冠層溫度，同時(shí)利用土鉆取土，取土深度為20 cm[25]，并通過(guò)烘干法測(cè)定土壤含水率。

1.4 熱紅外圖像處理

1.4.1 無(wú)人機(jī)圖像預(yù)處理

首先采用PIX4D 軟件對(duì)可見(jiàn)光與熱紅外圖像進(jìn)行拼接處理，而后以熱紅外圖像為基準(zhǔn)，利用ENVI Classic軟件手動(dòng)選取熱紅外與可見(jiàn)光圖像中明顯的地物特征點(diǎn)進(jìn)行配準(zhǔn)（特征點(diǎn)選取30 個(gè)以上，且保證每個(gè)小區(qū)包含至少2 個(gè)特征點(diǎn)，特征點(diǎn)均方根誤差小于1），最后利用FLIR Tools 軟件將熱紅外圖像的灰度值轉(zhuǎn)換為溫度，并用實(shí)測(cè)水溫校準(zhǔn)溫度。

1.4.2 熱紅外圖像冠層溫度提取

1）Otsu 算法

熱紅外相機(jī)獲取的圖像可分為兩類(lèi)：1）冠層像元，2）土壤像元[13]。Otsu 算法是灰度圖像閾值分割的經(jīng)典自適應(yīng)閾值算法，該算法以像元類(lèi)間方差為分割標(biāo)準(zhǔn)，將兩種像元類(lèi)間方差最大時(shí)的閾值作為冠層與土壤溫度的分界線，其核心公式如下[26]：

式中σ2為類(lèi)間方差；k為閾值；σ2(k)取最大值時(shí)的k值即為最佳閾值；w0為冠層像元出現(xiàn)的概率；w1為土壤像元出現(xiàn)的概率；μ為圖像灰度均值；μ0為冠層灰度均值；μ1為土壤灰度均值。

2）EXG-Kmeans 算法

可見(jiàn)光圖像的空間分辨率較高，且玉米冠層對(duì)綠波段反射強(qiáng)烈而土壤對(duì)綠波段反射較弱[27]，故本研究采用EXG 指數(shù)[28]區(qū)分冠層與土壤信息，EXG 指數(shù)計(jì)算式如下：

式中R、G、B分別表示可見(jiàn)光圖像紅、綠、藍(lán)波段灰度值。

同時(shí)為避免分類(lèi)結(jié)果受主觀因素的影響，采取非監(jiān)督分類(lèi)中的Kmeans 算法對(duì)EXG 指數(shù)圖像進(jìn)行分類(lèi)，生成包含玉米冠層像素輪廓特征的矢量文件，而后對(duì)熱紅外圖像進(jìn)行掩膜處理獲取冠層溫度信息。

3）Otsu-EXG-Kmeans 算法

研究表明Otsu 算法可以有效剔除圖像中陽(yáng)光直射土壤，而對(duì)陰影覆蓋土壤剔除效果較差[29]；借助可見(jiàn)光圖像可以準(zhǔn)確確定冠層區(qū)域，對(duì)陰影土壤與陽(yáng)光直射土壤的剔除效果較好，但由于存在配準(zhǔn)偏差，導(dǎo)致少量土壤區(qū)域被誤分到冠層區(qū)域中[23]，影響EXG-Kmeans 算法提取冠層溫度的精確度。因此本研究首先利用Otsu 算法剔除圖像中陽(yáng)光直射土壤，而后利用EXG-Kmeans算法剔除陰影土壤，稱(chēng)為Otsu-EXG-Kmeans 算法，冠層溫度提取過(guò)程通過(guò)MATLAB 與ENVI 實(shí)現(xiàn)，計(jì)算流程如圖2 所示。

1.5 作物水分脅迫指數(shù)

采用Jones 等[7]提出并由張仁華[9]改良后的理論模型計(jì)算CWSI，公式如下：

式中Tc為圖像冠層溫度，℃；Ta為大氣溫度，℃；γ為濕度計(jì)常數(shù)，kPa/℃；Δ為飽和水汽壓與溫度曲線斜率，kPa/℃；d為空氣飽和差，kPa；Rn為太陽(yáng)凈輻射，W/m3；G為土壤熱通量，W/m3；ra為空氣動(dòng)力阻力，s/m；ρ為空氣密度，kg/m3；Cp為空氣比熱，J/(kg·℃)。

1.6 圖像分類(lèi)精度評(píng)價(jià)

本研究采用生產(chǎn)者精度（Producer's Accuracy，PA）和用戶精度（User's Accuracy，UA）對(duì)熱紅外圖像冠層區(qū)域提取結(jié)果進(jìn)行量化評(píng)價(jià)，兩者越接近100%，說(shuō)明冠層區(qū)域提取精度越高，具體公式如下：

式中TP是預(yù)測(cè)為冠層的冠層樣本個(gè)數(shù)；FP是預(yù)測(cè)為冠層的土壤樣本個(gè)數(shù)；FN是預(yù)測(cè)為土壤的冠層樣本個(gè)數(shù)。

1.7 評(píng)價(jià)指標(biāo)

本研究通過(guò)相關(guān)系數(shù)（r）評(píng)價(jià)兩者的相關(guān)關(guān)系，r越接近1 兩者相關(guān)性越高；通過(guò)均方根誤差（Root Mean Square Error，RMSE）評(píng)價(jià)預(yù)測(cè)值的誤差，RMSE 越接近0，預(yù)測(cè)值的誤差越小。

2 結(jié)果與分析

2.1 冠層區(qū)域提取結(jié)果分析

2.1.1 不同算法冠層區(qū)域提取結(jié)果分析

以8 月2 日T1-1 小區(qū)中心區(qū)域?yàn)槔?，分別在RGB 影像與熱紅外影像中繪制不同算法確定的冠層區(qū)域圖（圖3）。對(duì)比RGB 影像可知，Otsu 算法確定的冠層區(qū)域包含極大部分玉米像元以及部分溫度與冠層溫度相近的土壤像元，而EXG-Kmeans 算法可在RGB 影像中準(zhǔn)確剔除陰影覆蓋土壤與陽(yáng)光直射土壤。由于RGB 影像與熱紅外影像在配準(zhǔn)時(shí)存在細(xì)微偏差，導(dǎo)致少量高于冠層溫度的土壤像元被誤劃于冠層區(qū)域，對(duì)比熱紅外影像可知，EXG-Kmeans 算法無(wú)法避免此類(lèi)現(xiàn)象發(fā)生，而Otsu-EXG-Kmeans 算法由于加入溫度閾值條件，可更好地剔除高于冠層溫度的土壤像元。

2.1.2 不同水分處理冠層區(qū)域提取結(jié)果分析

分別運(yùn)用 Otsu 算法、 EXG-Kmeans 算法和Otsu-EXG-Kmeans 算法對(duì)熱紅外圖像進(jìn)行分類(lèi)處理，每種方法的生產(chǎn)者精度和用戶精度如表1 所示。

表1 不同提取方法分類(lèi)精度Table 1 Classification accuracy of different extraction methods

依據(jù)圖3 所示冠層區(qū)域并對(duì)比表1 中同一算法的生產(chǎn)者精度與用戶精度可知：由于Otsu 算法確定的冠層區(qū)域包含大量土壤樣本，而土壤區(qū)域只有極少數(shù)的冠層樣本， 此算法的生產(chǎn)者精度遠(yuǎn)遠(yuǎn)大于用戶精度（96.1%>66.6%）；利用EXG-Kmeans 算法進(jìn)行地物分類(lèi)，誤分到冠層區(qū)域的土壤樣本與誤分到土壤區(qū)域的冠層樣本數(shù)量相近，生產(chǎn)者精度與用戶精度相差較??；Otsu-EXG-Kmeans 算法可剔除誤分到冠層區(qū)域的土壤樣本，而無(wú)法將誤分到冠層區(qū)域的土壤樣本召回，故用戶精度大于生產(chǎn)者精度（95.9%>85.6%）。

對(duì)比Otsu 算法對(duì)不同水分處理的分類(lèi)效果，其中T1處理的用戶精度高于T4 處理而生產(chǎn)者精度低于T4 處理(68.4%>64.0%，94.5%<97.2%)，由此推斷Otsu 算法在嚴(yán)重水分脅迫處理區(qū)域（T1）的適用性優(yōu)于充分灌溉處理（T4），但是誤分到土壤區(qū)域的冠層像元數(shù)量也要大于充分灌溉處理。EXG-Kmeans 算法中T4 處理的用戶精度與生產(chǎn)者精度均高于T1 處理(89.4%>77.9%，88.7%>79.1%)，故此算法對(duì)充分灌溉處理冠層區(qū)域提取效果要優(yōu)于嚴(yán)重水分脅迫處理。而Otsu-EXG-Kmeans 算法的用戶精度和生產(chǎn)者精度同水分脅迫情況呈負(fù)相關(guān)，且用戶精度均值達(dá)95.9%，對(duì)冠層提取效果優(yōu)于其他2 種算法。

2.2 冠層溫度頻率直方圖分析

水分脅迫處理越嚴(yán)重的小區(qū)冠層溫度越高，溫度極差越大[30]，不同算法溫度統(tǒng)計(jì)結(jié)果越明顯。以8 月2 日T1-1 小區(qū)為例，繪制原始圖像與不同分類(lèi)方法獲得的冠層溫度頻率直方圖（圖4），并統(tǒng)計(jì)冠層溫度的最小值、最大值以及平均值（表2）。

表2 圖像冠層溫度統(tǒng)計(jì)Table 2 Image canopy temperature statistics ℃

地面實(shí)測(cè)數(shù)據(jù)表明，陽(yáng)光直射土壤溫度大于冠層溫度，而陰影覆蓋土壤溫度與冠層溫度接近，甚至低于冠層溫度，故作物半覆蓋區(qū)域的溫度頻率直方圖大致為雙峰形狀。如圖4a 所示，原始溫度頻率直方圖左側(cè)峰頻率較大且溫度較低，代表冠層集中區(qū)域，右側(cè)峰頻率較小且溫度較高，代表土壤集中區(qū)域。原始圖像與EXG-Kmeans 算法相比于Otsu 算法與Otsu-EXG-Kmeans算法，提取的冠層溫度分布范圍更廣，進(jìn)一步說(shuō)明EXG-Kmeans 算法由于存在配準(zhǔn)偏差的現(xiàn)象，少量高溫的光照土壤像元無(wú)法有效剔除。由表 2 可知，Otsu-EXG-Kmeans 算法取得的冠層溫度最小值與EXG-Kmeans 算法同為31.0℃，最大值與Otsu 算法同為46.4℃，故此方法有效地剔除了像元溫度低于冠層溫度的陰影覆蓋土壤以及像元溫度高于冠層溫度的陽(yáng)光直射土壤。對(duì)比冠層溫度平均值，3 種方法均可有效降低土壤溫度的干擾，且Otsu-EXG-Kmeans 算法效果最佳，Otsu 算法次之，EXG-Kmeans 算法最差（冠層溫度：38.2℃<38.9℃<39.9℃<42.2℃）。

2.3 冠層溫度提取值與實(shí)測(cè)值相關(guān)性分析

為進(jìn)一步評(píng)估不同方法的提取效果，圖像提取溫度與地面實(shí)測(cè)溫度進(jìn)行相關(guān)性分析（n=60），如圖5 所示，隨著土壤像元的減少，圖像冠層溫度逐漸降低。對(duì)比3種處理方法，Otsu-EXG-Kmeans 算法提取的冠層溫度與地面實(shí)測(cè)溫度相關(guān)性最佳（r=0.788），誤差最?。≧MSE=1.901 ℃），同時(shí)Otsu-EXG-Kmeans 算法的擬合直線相比較其他方法，截距更小且斜率（回歸系數(shù)）更接近1，整體趨勢(shì)更接近1:1 線，故利用Otsu-EXG-Kmeans算法得到的冠層溫度精確度更高。

2.4 玉米水分脅迫監(jiān)測(cè)分析

2.4.1 冠層溫度指標(biāo)與土壤含水率相關(guān)性分析

上述研究表明，Otsu-EXG-Kmeans 算法對(duì)不同水分處理的適用性均較強(qiáng)，得到的冠層溫度精確度也較高，故可借助Otsu-EXG-Kmeans 算法監(jiān)測(cè)玉米水分脅迫狀況。作物通過(guò)根系吸收土壤水分來(lái)進(jìn)行蒸騰作用，從而維持冠層溫度的平衡，因此土壤含水率與玉米缺水狀況呈負(fù)相關(guān)[31]，是監(jiān)測(cè)玉米缺水狀況的有效指標(biāo)。Otsu-EXG-Kmeans 算法提取的冠層溫度值代入式（3）計(jì)算出相應(yīng)的作物水分脅迫指數(shù)，將冠層溫度與作物水分脅迫指數(shù)分別同0～20 cm 土壤含水率進(jìn)行相關(guān)性分析（n=60），以間接監(jiān)測(cè)玉米缺水狀況。結(jié)果表明（圖6），冠層溫度與作物水分脅迫指數(shù)均與土壤含水率負(fù)相關(guān)，即同玉米缺水狀況呈正相關(guān)，相比于冠層溫度，作物水分脅迫指數(shù)可以更為準(zhǔn)確地指示作物的水分虧缺狀況（r= -0.738）。

2.4.2 冠層溫度與CWSI 日平均變化趨勢(shì)對(duì)比

分別取各水分梯度3 個(gè)重復(fù)小區(qū)的平均值作為本梯度的冠層溫度值與作物水分脅迫指數(shù)，由此繪制不同水分處理的冠層溫度與作物水分脅迫指數(shù)日平均變化趨勢(shì)線如圖7 所示，不同水分處理的冠層溫度與作物水分脅迫指數(shù)具有明顯差異，并且T1>T2>T3>T4。根據(jù)試驗(yàn)記錄顯示，7 月28 日對(duì)T3、T4 處理進(jìn)行了灌水處理，故2種處理在29 日出現(xiàn)拐點(diǎn)，而T1 與T2 梯度在試驗(yàn)期間一直未進(jìn)行灌水處理，隨著時(shí)間的推移，土壤水分不斷被作物消耗，玉米受水分脅迫程度逐漸增加，相比于冠層溫度，作物水分脅迫指數(shù)日平均變化趨勢(shì)更好地描述了這一現(xiàn)象。

2.4.3 冠層溫度與CWSI 分布圖

為更直觀地展示試驗(yàn)區(qū)域不同處理的作物水分脅迫程度，以8 月2 日為例，繪制Otsu-EXG-Kmeans 算法獲取的冠層溫度與作物水分脅迫指數(shù)分布圖，如圖8 所示，整體冠層溫度在28～46.5 ℃之間分布，作物水分脅迫指數(shù)在0～1 之間分布。

3 討 論

本研究結(jié)果表明，Otsu 算法對(duì)水分虧缺嚴(yán)重的區(qū)域分類(lèi)效果較好，對(duì)充分灌溉區(qū)域分類(lèi)效果較差，然而，EXG-Kmeans 算法對(duì)不同水分脅迫處理的分類(lèi)效果與Otsu 算法截然相反。通過(guò)探究可知，在供水充足的狀態(tài)下，植物通過(guò)蒸騰作用散失大量輻射熱，從而降低葉溫和植物體溫，避免高溫傷害，同樣地表也可通過(guò)蒸發(fā)作用延緩?fù)寥罍囟壬仙?，與此同時(shí)，作物體內(nèi)的水分為葉片光合作用提供原料，葉綠體可正常工作，熱紅外圖像中的冠層與土壤溫度相差較小，而EXG-Kmeans 算法主要依據(jù)可見(jiàn)光圖像中冠層與土壤像元色彩差異對(duì)圖像進(jìn)行分類(lèi)，充分灌溉條件下，可見(jiàn)光圖像中的冠層像元為綠色，與棕黃色土壤色差較大，配準(zhǔn)法較為實(shí)用。作物水分脅迫區(qū)域，作物蒸騰作用與土壤蒸發(fā)作用減弱，溫度均升高，而由于土壤的比熱容相對(duì)水較小，升溫比葉片速率快，同時(shí)缺水會(huì)導(dǎo)致葉片氣孔導(dǎo)度下降，進(jìn)入葉內(nèi)的二氧化碳減少，光合產(chǎn)物在葉片中積累對(duì)光合作用產(chǎn)生反饋抑制作用，也會(huì)影響細(xì)胞伸長(zhǎng)并抑制蛋白質(zhì)的合成，葉片易出現(xiàn)暫時(shí)性萎蔫，使光合面積降低[32]，嚴(yán)重缺水時(shí)甚至造成葉綠體類(lèi)囊體結(jié)構(gòu)被破壞，葉片對(duì)紅光與藍(lán)光吸收能力下降[33]，熱紅外圖像中的冠層與土壤溫度相差較大，可見(jiàn)光圖像中的冠層像元偏向黃色，與棕黃色土壤色差較小，而Otsu 算法主要依據(jù)熱紅外圖像中冠層與土壤溫度差異對(duì)圖像進(jìn)行分類(lèi)，故在水分虧缺時(shí)Otsu 算法分類(lèi)效果較好。上述情況是水分脅迫程度與冠層溫度呈正相關(guān)的主要原因， 也是Otsu-EXG-Kmeans 算法在不同水分脅迫處理中適用性均較強(qiáng)的主要原因。

通過(guò)冠層溫度提取值與實(shí)測(cè)值相關(guān)性分析發(fā)現(xiàn)，水分脅迫嚴(yán)重的小區(qū)，部分冠層溫度提取值大于實(shí)測(cè)值，而充分灌溉的小區(qū)，部分冠層溫度提取值小于實(shí)測(cè)值。這一現(xiàn)象出現(xiàn)的原因是實(shí)測(cè)值選取的是生長(zhǎng)良好的承光面冠層區(qū)域，然而水分脅迫過(guò)重的區(qū)域玉米長(zhǎng)勢(shì)相對(duì)較差，且容易發(fā)生葉片卷曲現(xiàn)象（T1 處理中心區(qū)域），在提取過(guò)程中這部分溫度較高的區(qū)域參與計(jì)算，導(dǎo)致冠層溫度提取值高于實(shí)測(cè)值。反之充分灌溉區(qū)域玉米長(zhǎng)勢(shì)較好，冠層覆蓋度較高，部分陰影投映在冠層區(qū)域，在提取過(guò)程中溫度較低的陰影冠層參與計(jì)算，導(dǎo)致冠層溫度提取值低于實(shí)測(cè)值。

本研究方法也存在不足之處：首先，在水分嚴(yán)重脅迫區(qū)域，部分玉米葉片變黃，其冠層溫度與土壤溫度接近，本研究提出的方法無(wú)法準(zhǔn)確確定此部分冠層區(qū)域；其次，通過(guò)上述分析可知，由于剔除了可見(jiàn)光影像誤分到冠層中的土壤樣本，本研究方法確定的冠層區(qū)域小于僅借助可見(jiàn)光影像確定的冠層區(qū)域，即采用本研究方法計(jì)算的作物覆蓋度略低；此外，不同作物的耐旱性不同導(dǎo)致溫度分布不同，不同地區(qū)下墊面不同導(dǎo)致光譜信息不同，因此本研究提出的方法是否適用于其他地區(qū)的水分脅迫監(jiān)測(cè)仍需進(jìn)一步研究與驗(yàn)證。

4 結(jié) 論

本研究通過(guò)不同方法提取玉米冠層溫度，計(jì)算了基于冠層溫度的作物水分脅迫指數(shù)，探究了其在水分脅迫監(jiān)測(cè)中的精確度， 得出以下結(jié)論： 1 ） 利用Otsu-EXG-Kmeans 算法獲取的冠層區(qū)域圖像最優(yōu)，EXG-Kmeans 算法次之，Otsu 算法最差（用戶精度：95.9%>83.2%>66.6%）；Otsu-EXG-Kmeans 算法提取的冠層溫度與實(shí)測(cè)溫度相關(guān)關(guān)系接近1:1 線，相關(guān)性相比于其他方法更好（r=0.788），因此Otsu-EXG-Kmeans 算法是準(zhǔn)確獲取圖像冠層溫度的有效方法。2）相比于冠層溫度，作物水分脅迫指數(shù)與土壤含水率的相關(guān)性更高，作物水分脅迫指數(shù)日平均變化趨勢(shì)更符合實(shí)際情況，可更加精確地監(jiān)測(cè)玉米缺水狀況。