基于無(wú)人機(jī)熱紅外遙感的夏玉米蒸散量估算及其影響因子

段晨斐，胡振華，魏 征，張寶忠，陳 鶴，李 榮

(1.山西農(nóng)業(yè)大學(xué)林學(xué)院,山西 晉中 030801; 2.中國(guó)水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100038；3.國(guó)家節(jié)水灌溉北京工程技術(shù)研究中心,北京 100048)

0 引 言

蒸散量(ET)是水文循環(huán)、水資源管理的重要環(huán)節(jié)，也是土壤-作物-大氣系統(tǒng)之間物質(zhì)和能量交換的關(guān)鍵因素[1,2]。加強(qiáng)農(nóng)田蒸散量的準(zhǔn)確模擬預(yù)測(cè)研究，提高蒸散模型的精度、適用性對(duì)于作物需水量、農(nóng)業(yè)用水效率、灌溉用水管理和水文模擬等具有重要意義[3-5]。

Jackson等(1977年)提出基于冠氣溫差(Tc-Ta)估算作物ETd的方法[6-8]，并利用Penman-Monteith方程和冠層能量平衡單層模型，推導(dǎo)出rc/ra的冠層溫度表達(dá)式[9-11]，已有學(xué)者對(duì)該模型進(jìn)行驗(yàn)證，并應(yīng)用模型進(jìn)行農(nóng)田日蒸散量的估算[11]。隨著科學(xué)進(jìn)步和計(jì)算機(jī)、GIS、遙感、遙測(cè)、遙控等技術(shù)的快速發(fā)展，推動(dòng)了作物ET的研究。無(wú)人機(jī)遙感采集圖像技術(shù)具有靈活、高效、高精度、低成本等優(yōu)點(diǎn)，在國(guó)內(nèi)外被廣泛應(yīng)用[12,13]。國(guó)外學(xué)者M(jìn).P.Gonzalez-Dugo等(2006年)利用無(wú)人機(jī)熱紅外研究探討了棉花田冠層溫度變化與作物水分脅迫的關(guān)系[14]，Kohei Arai(2013年)等應(yīng)用無(wú)人機(jī)近紅外監(jiān)測(cè)茶葉含氮量[15]。Hoffmann et al等(2016年)應(yīng)用無(wú)人機(jī)獲取地面高分辨率溫度分別用于TSEB和DTD算法計(jì)算地表蒸散量[16]。Nieto et al等(2018年)基于無(wú)人機(jī)獲取地面溫度，將其用于TSEB模型計(jì)算地表蒸散量[17]。國(guó)內(nèi)，張智韜等(2018年)以棉花為研究對(duì)象借助無(wú)人機(jī)熱紅外圖像計(jì)算冠層溫度對(duì)水分脅迫進(jìn)行了研究[13]。孫圣等(2018年)基于無(wú)人機(jī)熱紅外圖像對(duì)土壤水分進(jìn)行了預(yù)測(cè)模型的建立與應(yīng)用[18]。

利用無(wú)人機(jī)熱紅外影像對(duì)夏玉米的蒸散量研究相對(duì)較少。因此，本文基于無(wú)人機(jī)熱紅外傳感器反演夏玉米的冠層溫度，采用基于冠層溫度的作物蒸散模型估算ETd,t，并與渦度相關(guān)系統(tǒng)實(shí)測(cè)值ETd,e、FAO Penman-Monteith模型計(jì)算值ETd,f進(jìn)行驗(yàn)證。在此基礎(chǔ)上對(duì)比分析了不同影響因子對(duì)模型估算ETd,t的影響，并分析了相關(guān)影響因素，為無(wú)人機(jī)遙感影像反演冠層溫度進(jìn)行區(qū)域尺度估算蒸散量提供參考依據(jù)。

1 材料與方法

1.1 試驗(yàn)地概況

試驗(yàn)在國(guó)家節(jié)水灌溉工程技術(shù)中心(北京)大興實(shí)驗(yàn)基地完成，試驗(yàn)地位于北京市南部的大興區(qū)內(nèi)，距離市中心約30 km，北緯39°37.25′，東經(jīng)116°25.51′，海拔30 m，屬暖溫帶半濕潤(rùn)大陸季風(fēng)氣候，年平均氣溫為12.1 ℃，平均水面蒸發(fā)量為1 800 mm，年平均降水量540 mm，多集中在汛期6-9月，約占全年降雨量80%左右，冬季干旱少雨，夏季濕潤(rùn)多雨。實(shí)驗(yàn)基地占地200 m×200 m，土壤土層深厚且以沙壤土為主，0～1 m平均田間持水量為33.4%左右。

田間試驗(yàn)于2018年6-9月中旬進(jìn)行，種植的夏玉米品種為“紀(jì)元168”。試驗(yàn)區(qū)域分為8個(gè)小區(qū)，每個(gè)小區(qū)面積平均為25 m2，外圍設(shè)置保護(hù)行，減少外邊界的影響。灌水利用地下管道輸送到田間，出水口處裝有水表，用來(lái)計(jì)量灌水量。8個(gè)小區(qū)田間管理措施均一致，播種前進(jìn)行灌水，施底肥。試驗(yàn)田內(nèi)裝有渦度相關(guān)系統(tǒng)、氣象站等儀器，用來(lái)測(cè)取所需數(shù)據(jù)資料。

1.2 實(shí)驗(yàn)觀測(cè)數(shù)據(jù)

1.2.1 冠層熱紅外影像

冠層熱紅外圖像采用大疆M600無(wú)人機(jī)搭載(FLIR Camera Tau2，USA)熱紅外相機(jī)測(cè)量，采集圖像像素為640×512，熱紅外波長(zhǎng)范圍7.5～13.5 μm。在夏玉米每個(gè)生育期內(nèi)，選擇晴朗、無(wú)風(fēng)或微風(fēng)的天氣采集數(shù)據(jù)2～3 次，數(shù)據(jù)采集時(shí)間為上午10∶00-12∶00，無(wú)人機(jī)飛行高度60 m，速度3 m/s，飛行時(shí)間約為3 min[19,20]。

圖1 試驗(yàn)小區(qū)布置Fig.1 Test plot distribution

1.2.2 其他數(shù)據(jù)

(1)田間尺度日蒸散量(ETd,e)實(shí)測(cè)值的獲取。采用渦度相關(guān)系統(tǒng)(Campbell Scientific Inc，USA)自動(dòng)測(cè)取[9]。

(2)夏玉米株高的測(cè)量。在每個(gè)小區(qū)內(nèi)選取3株具有代表性的植株，平均5 d人工測(cè)量一次植株高度。

(3)多光譜數(shù)據(jù)的獲取。無(wú)人機(jī)熱紅外數(shù)據(jù)獲取的同時(shí)，通過(guò)無(wú)人機(jī)多光譜遙感平臺(tái)獲取多光譜影像數(shù)據(jù)。

(4)氣象數(shù)據(jù)采集。實(shí)驗(yàn)基地內(nèi)安裝有小型氣象站，每隔0.5 h自動(dòng)監(jiān)測(cè)一次、包括空氣溫度、空氣濕度、大氣壓、太陽(yáng)輻射、風(fēng)速、降雨量等常規(guī)氣象資料[21]。8月份降雨頻繁，補(bǔ)充了土壤水分，圖2為夏玉米生育期內(nèi)日均氣溫與降雨量。

圖2 夏玉米生育期日均氣溫與降雨量Fig.2 Average temperature and rainfall during the growth period of summer maize

2 研究方法

2.1 數(shù)學(xué)模型

2.1.1 基于冠層溫度的作物蒸散模型

根據(jù)冠層和空氣阻力計(jì)算蒸散發(fā)的Penman-Monteith方程表達(dá)式為[22,23]：

(1)

基于rc/ra的冠層溫度表達(dá)式為[24,25]：

rc/ra=f(Tc-Ta)=

(2)

聯(lián)立(1)、(2)式得[10,26]：

(3)

式中：ETc為作物蒸散量，mm/d；Rn為凈輻射，MJ/(m2·d)；G土壤熱通量，MJ/(m2·d)；es為飽和水汽壓，kPa；ea為實(shí)際水汽壓，kPa；es-ea為飽和水汽壓差，kPa；Δ飽和水汽壓曲線斜率，kPa/℃；γ為濕度計(jì)常數(shù)，kPa/℃；ρa(bǔ)為空氣密度，kg/m3；Cp為空氣定壓比熱容，1.013×10-3MJ/(kg·℃)；rc為冠層阻抗；ra為空氣動(dòng)力學(xué)阻抗，s/m；Tc為冠層溫度；Ta為空氣溫度，℃。

計(jì)算空氣動(dòng)力學(xué)阻抗表達(dá)式為[2,24]：

(4)

式中：ra為空氣動(dòng)力學(xué)阻抗，s/m；zm為測(cè)風(fēng)速的高度，m；zh為測(cè)濕度的高度，m；zom為控制動(dòng)量傳遞的粗糙長(zhǎng)度，m；zoh控制熱通量和水汽傳輸?shù)拇植陂L(zhǎng)度，m；K為von Karman常數(shù)，取0.41；uz為z高度處的風(fēng)速，m/s。

2.1.2 渦度相關(guān)法

計(jì)算冠層蒸散量ET的表達(dá)式為[9,25,27]：

(5)

式中：λETc為冠層尺度潛熱通量，單位W/m2；ω′為垂直風(fēng)速脈動(dòng)量，m/s；q′為水汽密度脈動(dòng)值，g/m3。

2.1.3 Penman-Monteith法

利用由聯(lián)合國(guó)糧食及農(nóng)業(yè)組織FAO Penman-Monteith方法計(jì)算ET0(Allen et al.1998)[2,28,29]：

(6)

式中：ET0為參考作物需水量[30]，mm/d；T為2 m高處的日平均氣溫，℃；u2為2 m高處的風(fēng)速，m/s。

計(jì)算蒸散量ETc作物系數(shù)法[2]：

ETc=KcET0

(7)

式中：Kc為作物系數(shù)，取值Kcini=0.45，Kcmid=1.28，Kcend=0.5，取值具體詳見(jiàn)參考文獻(xiàn)[31]。

2.2 冠層熱紅外影像預(yù)處理

采用手持式地面測(cè)溫儀(Raytek ST6)對(duì)無(wú)人機(jī)熱紅外圖像進(jìn)行校準(zhǔn)，保證熱紅外圖像數(shù)據(jù)的準(zhǔn)確度。圖3(a)為試驗(yàn)區(qū)RGB圖像，試驗(yàn)小區(qū)布置如圖所示。圖3(b)為熱紅外遙感影像，選取各個(gè)小區(qū)的冠層，剔除保護(hù)行與小區(qū)邊緣，對(duì)熱紅外圖像點(diǎn)云信息提取計(jì)算和統(tǒng)計(jì)分析。

圖3 夏玉米熱紅外圖像Fig.3 Infrared image of summer maize

3 結(jié)果與分析

3.1 生育期內(nèi)熱紅外冠層溫度分布情況

夏玉米試驗(yàn)期內(nèi)利用反演的冠層溫度影像，繪制夏玉米熱紅外溫度直方圖，圖4為夏玉米各主要生育期內(nèi)熱紅外溫度分布情況。從圖4(a)中可知，夏玉米冠層溫度變化范圍為25.8～38.6 ℃，極差12.8 ℃，冠層溫度像元集中在28.9 ℃，頻率為3.4%。由于冠層溫度直方圖的主要貢獻(xiàn)像元集中在溫度直方圖的中間段，而前端部分像元溫度遠(yuǎn)離主體部分，尾部像元稀少且零散，誤差較大，因此做剔除處理，以最大頻率分布為中心取95%的累計(jì)像元為小區(qū)有效溫度像元取值，其余部分為有效溫度像元取值，其有效溫度變化范圍27.6～31.4 ℃，極差為3.8 ℃。由圖4(b)可看出，其有效溫度變化范圍28.7～36.4 ℃，極差為7.7 ℃，像元最大頻率為1.9%，溫度值32.5 ℃。由圖4(c)可得，其有效溫度變化范圍25.9～34.7 ℃，極差為8.8 ℃，像元最大頻率為1.9%，溫度值32.5 ℃。由圖4(d)可得，其有效溫度變化范圍28.4～34.6 ℃，極差為6.2 ℃，像元最大頻率為2.3%，溫度值31.0 ℃。

夏玉米各主要生育期內(nèi)冠層溫度均值在30.845 ℃左右(10∶00-12∶00)，從均值來(lái)看差異不是很大，出苗期溫度均值最低29.4 ℃，拔節(jié)期溫度均值最高32.2 ℃，基本趨于穩(wěn)定。分別從熱紅外冠層溫度直方圖空間分布來(lái)看，圖4(b)、圖4(c)夏玉米拔節(jié)期和灌漿期溫度空間差異較明顯，是夏玉米生長(zhǎng)最活躍的時(shí)期，拔節(jié)期內(nèi)玉米個(gè)體生長(zhǎng)快慢有所差異，熱紅外冠層溫度值較分散，總體上看呈正態(tài)分布。夏玉米在灌漿期內(nèi)，玉米籽粒生長(zhǎng)差異最大，加之玉米個(gè)體株高不同，導(dǎo)致熱紅外溫度值最為分散。而在生長(zhǎng)初期和生長(zhǎng)末期圖4(a)、圖4(d)，夏玉米生長(zhǎng)緩慢或停滯，玉米個(gè)體形態(tài)長(zhǎng)勢(shì)接近，熱紅外冠層溫度較為集中。

圖4 不同生育期夏玉米冠層溫度直方圖Fig.4 Histogram of summer maize canopy temperature at different growth stages

3.2 蒸散量ETd,t的估算

為驗(yàn)證基于熱紅外溫度的作物蒸散模型估算ETd,t的準(zhǔn)確性，利用(1)、(2)和(5)式結(jié)合小型氣象站得到的典型日氣象要素，分別計(jì)算得到夏玉米在整個(gè)生育期內(nèi)不同方法估算的日蒸散量ETd。

為保障萬(wàn)無(wú)一失，試驗(yàn)前米瑛無(wú)數(shù)次到采氣平臺(tái)，收集生產(chǎn)數(shù)據(jù)，反復(fù)核算對(duì)比，才最終制定分析方案，并一次試驗(yàn)成功。目前，該平臺(tái)4口井日產(chǎn)量14萬(wàn)立方米，為氣田穩(wěn)產(chǎn)增添了“底氣”。2017年，米瑛還負(fù)責(zé)完成了“涪陵頁(yè)巖氣田微注入測(cè)試分析”等工作，參與了“涪陵頁(yè)巖氣田開(kāi)發(fā)規(guī)律及合理配產(chǎn)研究”等項(xiàng)目。

圖5分別給出熱紅外溫度的作物蒸散模型計(jì)算的ETd,t，F(xiàn)AO Penman-Monteith模型計(jì)算的ETd,f，基于渦度相關(guān)儀器實(shí)測(cè)數(shù)據(jù)計(jì)算的ETd,e，在夏玉米整個(gè)生育期內(nèi)的ETd變化過(guò)程。從圖5中可以看出，熱紅外冠層溫度作物蒸散模型計(jì)算的ETd,t最低值出現(xiàn)在幼苗期為3.42 mm/d，最高值出現(xiàn)在灌漿期為10.94 mm/d，從蒸散量空間分布來(lái)看試驗(yàn)區(qū)邊緣ETd,t高于內(nèi)部，熱紅外溫度的作物蒸散模型計(jì)算的ETd,t與夏玉米熱紅外溫度分布具有很好的一致性。在夏玉米出苗期、成熟期各個(gè)像元的ETd,t較為集中，在夏玉米拔節(jié)期、灌漿期各個(gè)像元的ETd,t較為分散，與熱紅外溫度直方圖空間分布有良好的一致性。由于夏玉米在幼苗期植株性狀差異較小，植被覆蓋度較低，熱紅外遙感影像大部分都為土壤，下墊面蒸散量(ETd)較為均勻。成熟期由于夏玉米生長(zhǎng)發(fā)育停滯，葉片枯黃脫落，株高等個(gè)體差異較小，蒸散量(ETd)分布較為一致。拔節(jié)期、灌漿期夏玉米處于生長(zhǎng)發(fā)育迅速的階段，個(gè)體差異較為明顯，導(dǎo)致ETd,t差異明顯。3種方法計(jì)算出的作物蒸散量整體趨勢(shì)一致，在有降雨時(shí)期，日蒸散量(ETd)降低，降雨過(guò)后日蒸散量(ETd)增加。7月下旬三種模型估算的ETd相差最低，ETd,t與ETd,e差值(0.292 mm/d)，小于ETd,t與ETd,f差值(2.051 mm/d)。8月初，三種模型估算的ETd差距最大，ETd,t與ETd,e差值(8.209 mm/d)，大于ETd,t與ETd,f差值(5.612 mm/d)。很大原因是由于夏玉米處于拔節(jié)期，生長(zhǎng)速度較快，空氣動(dòng)力學(xué)阻抗ra計(jì)算值差異較大。8月中旬到9月底夏玉米(ETd)出現(xiàn)相同的發(fā)展趨勢(shì)，抽穗期、灌漿期和成熟期的ETd值保持良好的一致性。FAO Penman-Monteith模型估算的日蒸散量(ETd,f)值總體最低，ETd,t與ETd,e總體上看兩者相差不大，ETd,t估算精度達(dá)到80%，8月中旬和9月上旬ETd,上下波動(dòng)較大，很可能是頻繁降雨的影響。

圖5 不同模型計(jì)算的日蒸散量ETdFig.5 Daily evapotranspiration ETd calculated by different models

圖6為ETd,t分別與ETd,e、ETd,f建立的一元線性回歸模型，ETd,t與ETd,f模擬效果良好，一元線性模型擬合線整體位于1∶1線附近(R2=0.739，RMSE=0.676 mm/d)。ETd,t與ETd,e一元回歸模型(R2=0.742，RMSE=0.109 mm/d)。表1中兩種不同回歸分析在P<0.01水平上呈極顯著性關(guān)系。

3.3 蒸散量ETd,t與影響因子

3.3.1 蒸散量ETd,t與株高、日凈輻射Rnd、風(fēng)速W、空氣動(dòng)力學(xué)阻抗ra

將夏玉米生育期內(nèi)蒸散量(ETd,t)與同時(shí)期株高測(cè)量值進(jìn)行線性擬合(圖7)，兩者存在線性關(guān)系(R2=0.652)，呈正相關(guān)關(guān)系。夏玉米生長(zhǎng)期內(nèi)蒸散量(ETd,t)與氣象要素日凈輻射(Rnd)、田間夏玉米冠層風(fēng)速(W)和空氣動(dòng)力學(xué)阻抗(ra)進(jìn)行線性擬合(圖7)。日蒸散量(ETd,t)與日凈輻射(Rnd)、風(fēng)速(W)相關(guān)性較低，而ETd,t與ra相關(guān)性較高(R2=0.696)呈負(fù)相關(guān)關(guān)系。結(jié)果表明，株高與空氣動(dòng)力學(xué)阻抗ra對(duì)夏玉米農(nóng)田蒸散量影響較大，而日凈輻射Rnd、風(fēng)速W的影響相對(duì)較小。

圖6 熱紅外分別與渦度相關(guān)系統(tǒng)、FAO彭曼模型計(jì)算ETd的相關(guān)關(guān)系Fig.6 The correlation between thermal infrared and vorticity correlation system and ETd calculated by FAO Penman-Monteith model is shown

表1 夏玉米熱紅外溫度蒸散模型ETd,t與渦度相關(guān)系統(tǒng)ETd,e、FAO Penman-Monteith模型ETd,f的一元線性模型Tab.1 The ETd,t of the thermal infrared temperature evapotranspiration model of summer maize and the ETd,e of the vorticity correlation system, and the unitary linear model of FAO penman-monteith model ETd,f

3.3.2 蒸散量ETd,t與葉面積指數(shù)LAI

夏玉米葉面積指數(shù)(LAI)[圖8(a)]從7月下旬到8月中旬，快速增加，這是由于夏玉米處于拔節(jié)、抽穗期，夏玉米生長(zhǎng)發(fā)育較快。隨著植被覆蓋度的增加，熱紅外估算的日蒸散量(ETd,t)顯著增加，這與CUI Mingyue等[32]研究結(jié)果一致，根據(jù)這一變化趨勢(shì)來(lái)看，熱紅外估算的(ETd,t)能夠準(zhǔn)確反映出夏玉米在拔節(jié)期和抽穗期實(shí)際生長(zhǎng)情況，LAI與ETd,t建立線性回歸模型(R2=0.700)呈現(xiàn)良好的相關(guān)性[圖8(b)]。8月中旬到9月下旬，夏玉米灌漿期、乳熟期，株高停止生長(zhǎng)，作物葉片隨著玉米的成熟變黃枯萎凋落，葉面積指數(shù)呈緩慢下降趨勢(shì)，ETd,t在這一時(shí)期也呈現(xiàn)相同緩慢下降趨勢(shì)。

3.3.3 蒸散量ETd,t與歸一化植被指數(shù)NDVI

夏玉米蒸散量(ETd,t)與基于無(wú)人機(jī)多光譜影像反演歸一化植被指數(shù)(NDVI)在生育期內(nèi)的變化如圖9(a)所示，在夏玉米生長(zhǎng)初期，NDVI值較低，作物耗水量相對(duì)較低，夏玉米蒸散量(ETd,t)與NDVI的變化趨勢(shì)一致并保持穩(wěn)定。夏玉米進(jìn)入拔節(jié)期，隨著作物生長(zhǎng)發(fā)育加快，NDVI值迅速升高，同時(shí)期ETd,t同步上升。夏玉米抽穗期、灌漿期夏玉米生長(zhǎng)發(fā)育趨于穩(wěn)定，成熟期葉黃枯萎，NDVI在中、后期變化相對(duì)減小，并呈緩慢下降趨勢(shì)，此時(shí)期內(nèi)ETd,t與NDVI變化趨勢(shì)相同。圖9(b)為NDVI與ETd,t一元線性模型(R2=0.656)，表明夏玉米NDVI與ETd,t之間有良好的正相關(guān)關(guān)系，夏玉米ETd,t隨著NDVI的變化而變化，與實(shí)際田間觀測(cè)保持一致。

圖7 夏玉米生育期內(nèi)日蒸散量與影響因子的關(guān)系Fig.7 Relationship between daily evapotranspiration and influencing factors in summer maize growth period

圖8 夏玉米生育期葉面積指數(shù)及相關(guān)性分析Fig.8 LAI and its correlation analysis in summer maize growing period

圖9 夏玉米生育期歸一化植被指數(shù)與ETd,t變化趨勢(shì)及相關(guān)性分析Fig.9 The change trend and correlation analysis of NDVI, ETd,t in summer maize growth stage

4 結(jié) 論

本文利用無(wú)人機(jī)搭載熱紅外遙感平臺(tái)，結(jié)合采集得到不同生育期熱紅外冠層溫度影像反演出夏玉米其日蒸散量(ETd,t)變化狀況，與渦度計(jì)算的實(shí)測(cè)值(ETd,e)、FAO Penman-Monteith模型計(jì)算值(ETd,f)進(jìn)行對(duì)比，并對(duì)其與影響因子之間的關(guān)系進(jìn)行分析，所得結(jié)論如下：

(1)基于無(wú)人機(jī)熱紅外冠層溫度的作物蒸散模型可以準(zhǔn)確快速的估算田間尺度作物日蒸散量(ETd,t)，與實(shí)測(cè)值(ETd,e)、FAO Penman-Monteith模型計(jì)算值(ETd,f)結(jié)果對(duì)比，結(jié)果顯示ETd,t與ETd,e、ETd,f都有良好的一致性。ETd,t分別與ETd,e、ETd,f建立一元線性模型均呈現(xiàn)良好的線性相關(guān)關(guān)系(R2=0.739、0.742)，得出ETd,t最低值出現(xiàn)在幼苗期為3.42 mm/d，最高值出現(xiàn)在灌漿期為10.94 mm/d，試驗(yàn)區(qū)邊緣ETd,t高于內(nèi)部，并與渦度相關(guān)實(shí)測(cè)值ETd,e、FAO Penman-Monteith模型計(jì)算值ETd,f進(jìn)行驗(yàn)證，在P<0.01水平上呈顯著線性關(guān)系(R2=0.739 、0.742，RMSE=0.676、0.109 mm/d)，ETd,t估算精度達(dá)到80%以上。

(2)反演蒸散量(ETd,t)與影響因子株高、Rnd、W、ra、LAI和NDVI進(jìn)一步分析，株高、LAI、NDVI是主要的生物因子，其中LAI為夏玉米農(nóng)田最主要的生物因子，ETd,t隨著LAI的增加而顯著增加，LAI與ETd,t呈線性正相關(guān)關(guān)系(R2=0.700)，這與CUI Mingyue等[32]研究結(jié)果一致。Rnd、W、ra中，空氣動(dòng)力學(xué)阻抗(ra)是最主要的環(huán)境驅(qū)動(dòng)因子，ra與ETd,t呈線性負(fù)相關(guān)關(guān)系(R2= 0.696)。

相比傳統(tǒng)計(jì)算作物蒸散量的方法，以無(wú)人機(jī)熱紅外遙感為平臺(tái)，具有較高的空間精度和應(yīng)用潛力。對(duì)小區(qū)域作物實(shí)現(xiàn)精準(zhǔn)灌溉具有重要意義，同時(shí)也為今后實(shí)現(xiàn)多區(qū)域、大尺度的估算作物蒸散量狀況提供了一條新途徑。