基于固定翼無(wú)人機(jī)多光譜影像的水稻長(zhǎng)勢(shì)關(guān)鍵指標(biāo)無(wú)損監(jiān)測(cè)

王偉康，張嘉懿，汪慧，曹強(qiáng)，田永超，朱艷，曹衛(wèi)星，劉小軍

基于固定翼無(wú)人機(jī)多光譜影像的水稻長(zhǎng)勢(shì)關(guān)鍵指標(biāo)無(wú)損監(jiān)測(cè)

王偉康，張嘉懿，汪慧，曹強(qiáng)，田永超，朱艷，曹衛(wèi)星，劉小軍

南京農(nóng)業(yè)大學(xué)農(nóng)學(xué)院/國(guó)家信息農(nóng)業(yè)工程技術(shù)中心/智慧農(nóng)業(yè)教育部工程研究中心/農(nóng)業(yè)農(nóng)村部農(nóng)作物系統(tǒng)分析與決策重點(diǎn)實(shí)驗(yàn)室/江蘇省信息農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室，南京 210095

【背景】近年來(lái)隨著遙感技術(shù)的快速發(fā)展，實(shí)時(shí)無(wú)損監(jiān)測(cè)作物生長(zhǎng)狀況已成為當(dāng)前研究熱點(diǎn)，遙感獲取的農(nóng)情信息將為實(shí)現(xiàn)大面積作物精確管理提供指導(dǎo)。在眾多遙感監(jiān)測(cè)平臺(tái)里，無(wú)人機(jī)因其操作簡(jiǎn)單、使用成本低等特點(diǎn)而受到廣泛關(guān)注，無(wú)人機(jī)搭載多光譜相機(jī)可以快速獲取作物的長(zhǎng)勢(shì)信息?！灸康摹繃L試將固定翼無(wú)人機(jī)多光譜影像紋理信息與光譜信息結(jié)合，探究“圖譜”信息對(duì)水稻長(zhǎng)勢(shì)指標(biāo)的監(jiān)測(cè)效果?！痉椒ā客ㄟ^(guò)開(kāi)展兩年涉及不同播期、品種、播栽方式、施氮水平的水稻田間試驗(yàn)，在水稻關(guān)鍵生育期使用固定翼無(wú)人機(jī)搭載Sequoia多光譜相機(jī)獲取水稻冠層遙感影像，同步進(jìn)行地上部破壞性取樣以獲取水稻葉面積指數(shù)（LAI）、地上部生物量（AGB）和植株氮含量（PNC）等農(nóng)學(xué)指標(biāo)，采用簡(jiǎn)單線性回歸、偏最小二乘回歸和人工神經(jīng)網(wǎng)絡(luò)回歸算法，構(gòu)建基于固定翼無(wú)人機(jī)多光譜影像的水稻長(zhǎng)勢(shì)指標(biāo)監(jiān)測(cè)模型，比較分析光譜紋理信息在不同模型中的監(jiān)測(cè)效果?！窘Y(jié)果】利用簡(jiǎn)單線性回歸方法探究了植被指數(shù)（VI）、單波段紋理特征與水稻LAI、AGB和PNC間的定量關(guān)系，結(jié)果表明植被指數(shù)與LAI和AGB之間有較強(qiáng)的相關(guān)性，表現(xiàn)最好的植被指數(shù)為CIRE和NDRE，2分別為0.80和0.76，但對(duì)于PNC的監(jiān)測(cè)，植被指數(shù)并未達(dá)到理想的效果，表現(xiàn)最好的RESAVI和NDRE與PNC的決定系數(shù)僅為0.13。通過(guò)簡(jiǎn)單線性回歸進(jìn)一步發(fā)現(xiàn)單波段的紋理特征在對(duì)水稻生長(zhǎng)指標(biāo)的監(jiān)測(cè)中表現(xiàn)并不理想；為進(jìn)一步分析影像紋理對(duì)上述3個(gè)指標(biāo)的監(jiān)測(cè)效果，參照植被指數(shù)的構(gòu)建方法構(gòu)建了歸一化紋理指數(shù)（NDTI）、比值紋理指數(shù)（RTI）和差值紋理指數(shù)（DTI），通過(guò)相關(guān)性分析發(fā)現(xiàn)新構(gòu)建的紋理指數(shù)（TI）相較于單波段紋理特征對(duì)水稻生長(zhǎng)指標(biāo)的監(jiān)測(cè)精度有所提升，但效果并未好于植被指數(shù)。為實(shí)現(xiàn)光譜與紋理間的結(jié)合，采用偏最小二乘和人工神經(jīng)網(wǎng)絡(luò)的建模方法，以VI、VI+TI為不同的輸入?yún)?shù)組合進(jìn)行水稻LAI、AGB和PNC的監(jiān)測(cè)模型構(gòu)建，結(jié)果表明采用偏最小二乘和人工神經(jīng)網(wǎng)絡(luò)的建模方法與簡(jiǎn)單線性回歸相比模型的監(jiān)測(cè)精度均得到了大幅提升，以VI+TI為輸入變量，采用人工神經(jīng)網(wǎng)絡(luò)構(gòu)建的模型在模型驗(yàn)證中取得了最佳效果，LAI模型的驗(yàn)證2由0.75提升至0.86，AGB和PNC的模型驗(yàn)證2也分別由0.72和0.26提升至0.92和0.86，同時(shí)模型的均顯著降低?！窘Y(jié)論】利用固定翼無(wú)人機(jī)采集水稻冠層多光譜影像，通過(guò)人工神經(jīng)網(wǎng)絡(luò)算法融合光譜和紋理信息能夠有效提升水稻LAI、AGB和PNC的監(jiān)測(cè)精度，該研究結(jié)果將為快速大面積作物長(zhǎng)勢(shì)監(jiān)測(cè)提供理論依據(jù)。

無(wú)人機(jī)；植被指數(shù)；紋理特征；長(zhǎng)勢(shì)；水稻

0 引言

【研究意義】水稻作為重要的糧食作物之一，維持著全球一半以上人口的生存，準(zhǔn)確且有效的監(jiān)測(cè)水稻長(zhǎng)勢(shì)對(duì)于提升糧食總產(chǎn)量和保障糧食安全至關(guān)重要[1]。葉面積指數(shù)（leaf area index，LAI）、地上部生物量（aboveground biomass，AGB）和植株氮含量（plant nitrogen content，PNC）是水稻重要的生理生化指標(biāo)，其在水稻氮素營(yíng)養(yǎng)診斷[2-4]、生產(chǎn)管理[5]和產(chǎn)量預(yù)測(cè)[6]等方面具有重要意義。因此快速、無(wú)損地獲取水稻的長(zhǎng)勢(shì)信息對(duì)實(shí)現(xiàn)水稻精確管理意義重大。隨著農(nóng)業(yè)集約化進(jìn)程的推進(jìn)，傳統(tǒng)的田間測(cè)量和化學(xué)分析方法來(lái)獲取作物長(zhǎng)勢(shì)信息已難以滿足大面積的作物生產(chǎn)管理需求。近年來(lái)，隨著傳感器技術(shù)的飛速發(fā)展，多種類型的遙感監(jiān)測(cè)平臺(tái)被廣泛應(yīng)用于農(nóng)業(yè)研究領(lǐng)域，根據(jù)監(jiān)測(cè)尺度的不同，常用的遙感平臺(tái)分為地面平臺(tái)、無(wú)人機(jī)遙感平臺(tái)以及衛(wèi)星遙感平臺(tái)。無(wú)人機(jī)遙感平臺(tái)兼具時(shí)空分辨率高、覆蓋范圍大、操作成本低等優(yōu)勢(shì)，在農(nóng)情遙感監(jiān)測(cè)方面得到了更為廣泛的應(yīng)用?！厩叭搜芯窟M(jìn)展】當(dāng)前已有大量學(xué)者利用無(wú)人機(jī)實(shí)現(xiàn)了對(duì)作物生長(zhǎng)指標(biāo)的反演監(jiān)測(cè)，如Fu等[7]利用六旋翼無(wú)人機(jī)搭載Airphen多光譜相機(jī)獲取了小麥冠層多光譜影像，通過(guò)提取小麥不同生育時(shí)期多光譜影像的NDVI值，利用神經(jīng)網(wǎng)絡(luò)等多種機(jī)器學(xué)習(xí)方法構(gòu)建了小麥LAI和葉片干物重估算模型，實(shí)現(xiàn)了對(duì)小麥長(zhǎng)勢(shì)的實(shí)時(shí)估測(cè)；Tao等[8]利用多旋翼無(wú)人機(jī)獲取了小麥關(guān)鍵生育期的高光譜影像，通過(guò)影像提取了光譜指數(shù)和小麥株高等參數(shù)，結(jié)合小麥實(shí)測(cè)株高，利用偏最小二乘、人工神經(jīng)網(wǎng)絡(luò)和隨機(jī)森林3種建模方法構(gòu)建了小麥產(chǎn)量估測(cè)模型，有效提升了產(chǎn)量的估測(cè)精度；Lu等[9]利用無(wú)人機(jī)搭載RGB相機(jī)獲取小麥田間彩色影像，結(jié)合作物表面模型提高了小麥生物量的估測(cè)精度。上述研究多關(guān)注于影像的光譜信息，缺乏對(duì)遙感影像的深度挖掘及分析過(guò)程。陳鵬飛等[10]研究表明，隨著作物生育進(jìn)程的推進(jìn)，作物冠層覆蓋度趨于飽和，導(dǎo)致光譜反射率信息對(duì)冠層生物量變化不敏感，影響反演的精度和效果。隨著傳感器及圖像處理技術(shù)的進(jìn)步，高質(zhì)量影像的獲取及處理流程不斷簡(jiǎn)化。低空無(wú)人機(jī)平臺(tái)獲取的高分辨率多光譜影像蘊(yùn)含了豐富的地物反射率信息、紋理信息以及結(jié)構(gòu)信息等，以往大部分研究?jī)H利用影像的光譜信息進(jìn)行作物長(zhǎng)勢(shì)監(jiān)測(cè)[11-12]。而與單一遙感數(shù)據(jù)相比，多元遙感數(shù)據(jù)融合能夠提高作物長(zhǎng)勢(shì)指標(biāo)監(jiān)測(cè)模型的可靠性和魯棒性[13]。例如萬(wàn)亮等[14]利用可見(jiàn)光圖像的顏色和紋理特征實(shí)現(xiàn)了對(duì)稻穗的精準(zhǔn)識(shí)別，將此結(jié)果與多光譜圖像反射率信息融合，利用隨機(jī)森林算法有效提高了穗生物量的評(píng)估精度，建模2達(dá)到0.84；戴震[15]將玉米冠層光譜特征信息和紋理特征信息通過(guò)不同線性組合作為變量輸入，利用深度神經(jīng)網(wǎng)絡(luò)算法建立玉米葉片氮含量反演模型，模型決定系數(shù)達(dá)到0.85；Zheng等[16]利用無(wú)人機(jī)RGB影像和近地光譜儀獲取了小麥冠層的紋理特征和反射率信息，利用PLSR方法建立了光譜與紋理特征融合的小麥AGB估測(cè)模型，結(jié)果表明特征融合對(duì)小麥AGB的估測(cè)效果優(yōu)于單一特征。影像的紋理特征作為影像信息重要的組成成分，其分析過(guò)程實(shí)質(zhì)上是在自定義的窗口內(nèi)計(jì)算相鄰像素之間像素值的變化[17]，而紋理信息作為影像的一種固有屬性，能夠提供地物的空間幾何信息[18]，作物因營(yíng)養(yǎng)狀況的差異會(huì)導(dǎo)致生長(zhǎng)發(fā)育狀況不同，進(jìn)一步影響植株高度、葉片大小等，這些差異均會(huì)造成植株冠層紋理特征的差異，因此紋理可作為有效的補(bǔ)充信息用于作物生長(zhǎng)監(jiān)測(cè)?！颈狙芯壳腥朦c(diǎn)】紋理特征已在森林遙感監(jiān)測(cè)中得到廣泛應(yīng)用，研究者基于紋理特征進(jìn)行森林植被類型識(shí)別，生物量及覆蓋度估算等。紋理信息在作物生長(zhǎng)監(jiān)測(cè)方面的應(yīng)用多集中于對(duì)生物量的反演，而較少應(yīng)用于作物氮素營(yíng)養(yǎng)監(jiān)測(cè)，已有的研究中陳鵬飛等[19]基于無(wú)人機(jī)高空間分辨率影像，探討剔除土壤背景信息及增加紋理信息對(duì)棉花植株氮濃度反演的影響，為棉花氮素營(yíng)養(yǎng)精準(zhǔn)探測(cè)提供了新技術(shù)手段；陳鵬等[20]將植被指數(shù)與紋理特征通過(guò)主成分融合構(gòu)建了一種新的綜合指標(biāo)估算馬鈴薯葉片葉綠素含量，其模型估算效果優(yōu)于僅依賴植被指數(shù)或紋理特征建立的葉片葉綠素估算模型。光譜與紋理相結(jié)合在水稻氮素營(yíng)養(yǎng)監(jiān)測(cè)方面的效果尚未見(jiàn)詳細(xì)報(bào)道，同時(shí)在前人研究中，無(wú)人機(jī)平臺(tái)多選用旋翼無(wú)人機(jī)，而固定翼無(wú)人機(jī)對(duì)影像的采集頻率和影像質(zhì)量提出了新的要求，利用固定翼無(wú)人機(jī)多光譜影像信息進(jìn)行水稻長(zhǎng)勢(shì)監(jiān)測(cè)的效果有待探究。【擬解決的關(guān)鍵問(wèn)題】以固定翼無(wú)人機(jī)作為遙感平臺(tái)，搭載Sequoia多光譜相機(jī)，以水稻生長(zhǎng)指標(biāo)為研究對(duì)象，綜合分析水稻生長(zhǎng)指標(biāo)與無(wú)人機(jī)多光譜影像的光譜、紋理信息間的定量關(guān)系，利用簡(jiǎn)單線性回歸（simple liner regression，SLR）、偏最小二乘回歸（partial least squares regression，PLSR）和人工神經(jīng)網(wǎng)絡(luò)（artificial neural network，ANN）3種建模方法比較分析不同輸入變量組合下模型對(duì)水稻LAI、AGB和PNC的反演效果，為水稻生長(zhǎng)指標(biāo)的快速大面積監(jiān)測(cè)提供技術(shù)支撐。

1 材料與方法

1.1 試驗(yàn)區(qū)概況

于2018、2020年間在江蘇省興化市百萬(wàn)畝國(guó)家糧食生產(chǎn)功能區(qū)進(jìn)行不同品種、施氮水平、播栽方式的水稻田間試驗(yàn)（圖1），該地區(qū)年降雨量645 mm，年平均氣溫16.9 ℃，實(shí)行稻麥兩熟種植制度。

2018年試驗(yàn)設(shè)2個(gè)水稻品種（南粳9108、甬優(yōu)2640），采用2種播栽方式（缽苗、毯苗）；2020年試驗(yàn)設(shè)3個(gè)水稻品種（南粳9108、甬優(yōu)2640、武運(yùn)粳32），采用人工播栽。兩年試驗(yàn)均設(shè)置4個(gè)施氮水平，包括0（N0）、135 kg·hm-2（N1）、270 kg·hm-2（N2）、405 kg·hm-2（N3），栽插規(guī)格為粳稻15 cm×30 cm，雜交稻18 cm×30 cm，試驗(yàn)為小區(qū)試驗(yàn)，采用隨機(jī)分布，設(shè)3次重復(fù)。氮肥形態(tài)為尿素，施氮比例：基肥﹕分蘗肥﹕促花肥﹕?；ǚ?3﹕3﹕2﹕2。另配施磷肥：P2O590 kg·hm-2，作基肥一次性施入。鉀肥：K2O 150 kg·hm-2，基肥50%，倒4葉50%。其他管理措施同一般高產(chǎn)田。

氮水平N rate：N0（0）；N1（135 kg·hm-2）；N2（270 kg·hm-2）；N3（405 kg·hm-2）。品種Variety：V1（南粳9108 Nanjing 9108）；V2（甬優(yōu)2640 Yongyou 2640）。播栽方式Planting technique：P1（缽苗移栽Tray seeding transplanting）；P2（毯苗移栽Blanket seeding transplanting）

1.2 數(shù)據(jù)獲取

在水稻關(guān)鍵生育期，每個(gè)小區(qū)取代表性稻株3穴獲取LAI、AGB等農(nóng)學(xué)指標(biāo)，同步進(jìn)行無(wú)人機(jī)飛行測(cè)試獲取多光譜影像。

1.2.1 農(nóng)學(xué)數(shù)據(jù)獲取 根據(jù)植株器官發(fā)育情況，將樣品植株分離為葉、莖鞘和穗。將新鮮葉片清洗后通過(guò)LAI-3000臺(tái)式葉面積儀測(cè)定葉面積，隨后將分離的植株器官立刻放入烘箱在105 ℃殺青30 min，80 ℃烘干至恒重后稱量，計(jì)算各器官的干重。樣品粉碎后采用凱氏定氮法，使用AA3流動(dòng)分析儀（BRAN+ LUEBBE AA3；德國(guó)）測(cè)定植株各器官的含氮量，最后根據(jù)樣本的干重計(jì)算各個(gè)小區(qū)的PNC值。

1.2.2 無(wú)人機(jī)多光譜影像獲取 使用eBee-SQ固定翼無(wú)人機(jī)（圖2）搭載Parrot Sequoia多光譜相機(jī)來(lái)獲取水稻冠層多光譜影像數(shù)據(jù)，相機(jī)包含5個(gè)鏡頭，一個(gè)1 600萬(wàn)像素滾動(dòng)快門(mén)的RGB鏡頭，以及4個(gè)150萬(wàn)像素全局快門(mén)的單波段鏡頭，4個(gè)波段分別為綠光波段、紅光波段、紅邊波段和近紅外波段[21]。無(wú)人機(jī)飛行測(cè)試均在晴天和無(wú)風(fēng)條件下進(jìn)行，測(cè)試時(shí)間為10：00—14：00，每次飛行前對(duì)地面上的標(biāo)準(zhǔn)反射板進(jìn)行拍攝來(lái)完成相機(jī)的輻射校正。無(wú)人機(jī)飛行高度設(shè)為80 m，對(duì)應(yīng)地面像素分辨率為8 cm，旁向和航向重疊度均設(shè)置為70%。無(wú)人機(jī)影像使用Pix4Dmapper軟件進(jìn)行幾何標(biāo)定和拼接處理，并最終生成4幅單波段圖像[22]。

圖2 eBee-SQ固定翼無(wú)人機(jī)平臺(tái)

Fig. 2 eBee-SQ fixed-wing UAV platform

1.3 數(shù)據(jù)分析方法

1.3.1 無(wú)人機(jī)影像光譜分析 將拼接好的影像導(dǎo)入ENVI 5.3（Exelis Visual Information Solutions，美國(guó)），利用Region of Interest工具選取實(shí)驗(yàn)小區(qū)的影像處理區(qū)域，計(jì)算每個(gè)小區(qū)的反射率均值，參照前人研究，選取6個(gè)較為經(jīng)典的植被指數(shù)（vegetation index，VI）[23-24]對(duì)水稻LAI、AGB和PNC進(jìn)行簡(jiǎn)單線性回歸建模（表1）。

表1 本研究使用的植被指數(shù)

1.3.2 無(wú)人機(jī)影像紋理分析 使用ENVI 5.3軟件，利用灰度共生矩陣（CLCM）對(duì)每個(gè)波段的無(wú)人機(jī)影像進(jìn)行紋理特征提取，將紋理提取窗口設(shè)置為3*3大小，窗口移動(dòng)步長(zhǎng)設(shè)置為1，將4個(gè)角度（0°、45°、90°、135°）上紋理特征值取均值作為結(jié)果輸出。每個(gè)波段有8個(gè)紋理特征值：均值（Mea）、方差（Var）、同質(zhì)性（Hom）、對(duì)比度（Con）、異質(zhì)性（Dis）、熵（Ent）、二階矩（Sec）以及相關(guān)性（Cor）[25]。參照植被指數(shù)的構(gòu)建方式，將不同波段的紋理特征進(jìn)行兩兩組合[26]，分別構(gòu)建歸一化差值紋理指數(shù)（NDTI）、比值紋理指數(shù)（RTI）和差值紋理指數(shù)（DTI）（表2）。

表2 本研究所構(gòu)建的紋理指數(shù)

T1和T2為所選波段的紋理特征值

T1 and T2 are the texture feature values of the selected band

1.3.3 模型構(gòu)建 （1）簡(jiǎn)單線性模型：植被指數(shù)法是光學(xué)遙感估算植被長(zhǎng)勢(shì)指標(biāo)常用的經(jīng)驗(yàn)方法，在本研究中，首先利用簡(jiǎn)單線性回歸（SLR）方法分析全生育期植被指數(shù)以及單波段紋理特征與水稻LAI、AGB和PNC間潛在的相關(guān)性，之后為探究紋理組合特征對(duì)水稻長(zhǎng)勢(shì)的監(jiān)測(cè)效果，分別構(gòu)建3種類型的紋理指數(shù)（texture index，TI），通過(guò)與水稻長(zhǎng)勢(shì)指標(biāo)間的相關(guān)性分析，篩選出6個(gè)表現(xiàn)最好的TI進(jìn)行后續(xù)模型的構(gòu)建。

（2）偏最小二乘模型：偏最小二乘回歸（PLSR）是一種多元統(tǒng)計(jì)方法，其適用于分析多個(gè)自變量與多個(gè)因變量間的相關(guān)關(guān)系，它包括多元線性回歸和主成分分析的基本功能[27]，可以消除自變量之間多重共線性對(duì)模型精度帶來(lái)的干擾。本研究中通過(guò)留一法交叉驗(yàn)證的均方根誤差，確定PLSR中的最佳主成分?jǐn)?shù)量，通過(guò)獨(dú)立數(shù)據(jù)進(jìn)行模型的驗(yàn)證。PLSR建模在R Studio中用pls包進(jìn)行。

（3）人工神經(jīng)網(wǎng)絡(luò)模型：人工神經(jīng)網(wǎng)絡(luò)（ANN）是對(duì)生物神經(jīng)元的模擬和簡(jiǎn)化，可以實(shí)現(xiàn)輸入和輸出間的非線性映射關(guān)系[28]。多層感知器（multi-layer perceptron，MLP）是運(yùn)用最為廣泛的一種人工神經(jīng)網(wǎng)絡(luò)模型，特別適合于求解內(nèi)部機(jī)制復(fù)雜的問(wèn)題。一個(gè)典型的MLP模型由3部分組成，輸入層（input layer)、隱含層（hidden layer）和輸出層（output layer）。在多層感知器內(nèi)部，前一層的每個(gè)節(jié)點(diǎn)與相鄰后一層的每個(gè)節(jié)點(diǎn)相連接，利用非線性激活函數(shù)和反向傳播的監(jiān)督學(xué)習(xí)技術(shù)進(jìn)行訓(xùn)練，通過(guò)更新網(wǎng)絡(luò)權(quán)重參數(shù)，得到最優(yōu)結(jié)果。在本研究中，將隱含層數(shù)設(shè)置為1，隱含層激活函數(shù)設(shè)置為雙曲正切，輸出層的激活函數(shù)設(shè)置為恒等式，通過(guò)模型的決定系數(shù)（2）和均方根誤差（）來(lái)確定最佳的隱含層神經(jīng)元數(shù)量。ANN模型的構(gòu)建在SPSS中實(shí)現(xiàn)。

1.3.4 模型精度評(píng)價(jià) 試驗(yàn)數(shù)據(jù)取自2018年和2020年兩年水稻試驗(yàn)，共涉及分蘗、拔節(jié)、孕穗、抽穗4個(gè)關(guān)鍵生育期，其中70%作為建模數(shù)據(jù)，30%用于模型的檢驗(yàn)，建模精度與驗(yàn)證精度分別采用決定系數(shù)2和均方根誤差作為評(píng)估標(biāo)準(zhǔn)[29]，在所構(gòu)建的PLSR和ANN模型中，通過(guò)引入變量投影重要性（variable importance in the projection，VIP）指標(biāo)來(lái)衡量各輸入?yún)?shù)在模型中的貢獻(xiàn)度，通常情況下VIP值越高表明因子貢獻(xiàn)度越大[30]。

2 結(jié)果

2.1 水稻長(zhǎng)勢(shì)指標(biāo)的變異性分析

兩年水稻試驗(yàn)設(shè)置不同氮處理，樣本涵蓋范圍廣泛，使得模型普適性增強(qiáng)。如表3所示，不同品種、氮肥處理、播栽方式、播期水稻長(zhǎng)勢(shì)差異明顯，本研究所采用的數(shù)據(jù)集，LAI介于0.40—9.07，其變異系數(shù)為52.57%；AGB介于0.66—21.94，其變異系數(shù)為70.78%；PNC介于0.56—3.88，其變異系數(shù)為34.20%。整體來(lái)看，數(shù)據(jù)變化幅度較大，可解釋更多的可能性。因此，該數(shù)據(jù)可以為水稻長(zhǎng)勢(shì)監(jiān)測(cè)模型的建立提供可靠依據(jù)。

表3 不同處理下水稻葉面積指數(shù)、地上部生物量和植株氮含量的統(tǒng)計(jì)分析

2.2 基于植被指數(shù)（VI）的水稻長(zhǎng)勢(shì)監(jiān)測(cè)模型構(gòu)建與檢驗(yàn)

基于水稻冠層無(wú)人機(jī)多光譜影像數(shù)據(jù)，提取水稻冠層反射率，篩選并構(gòu)建了6個(gè)VI。表4顯示了全生育期VI與水稻LAI、AGB和PNC間的相關(guān)性?？梢钥闯鲈?個(gè)指標(biāo)的監(jiān)測(cè)中，VI與LAI的相關(guān)性最好，其次是AGB，對(duì)于PNC的反演效果最差。與LAI相關(guān)性較高的VI有CIRE、NDRE、RESAVI，2均＞0.70，其中表現(xiàn)最好的為CIRE，2=0.80；與AGB相關(guān)性較高的VI與LAI相同，其中NDRE的2達(dá)到0.76；植被指數(shù)RESAVIE和NDRE與PNC相關(guān)性高，2為0.13（圖3）。上述與LAI、AGB、PNC相關(guān)性較高的幾個(gè)VI均包含紅邊波段，由此可以確定在這兩個(gè)波段蘊(yùn)含豐富的水稻冠層植被信息。此外，在水稻孕穗期到抽穗期之間發(fā)現(xiàn)AGB和PNC均有一段跳躍式上升或下降過(guò)程，分析原因可能是在此階段水稻處于營(yíng)養(yǎng)生長(zhǎng)與生殖生長(zhǎng)并進(jìn)階段，生長(zhǎng)發(fā)育較為迅速，而取樣測(cè)試時(shí)期間隔較長(zhǎng)導(dǎo)致。將表現(xiàn)最好的VI進(jìn)行SLR建模，并使用獨(dú)立數(shù)據(jù)進(jìn)行驗(yàn)證，驗(yàn)證結(jié)果如圖3，其中LAI的模型驗(yàn)證2=0.75，=0.99，模型驗(yàn)證取得較好的結(jié)果。AGB的模型驗(yàn)證2=0.72，=3.12；PNC的模型驗(yàn)證2=0.26，=0.55。模型的驗(yàn)證結(jié)果表明所選取的VI在對(duì)LAI、AGB的監(jiān)測(cè)中均有較好的表現(xiàn)，而在PNC的驗(yàn)證中所選擇的6個(gè)VI表現(xiàn)均未達(dá)到理想效果。

表4 6個(gè)植被指數(shù)與水稻葉面積指數(shù)、地上部生物量和植株氮含量的相關(guān)性

2.3 基于紋理的水稻長(zhǎng)勢(shì)監(jiān)測(cè)模型構(gòu)建與檢驗(yàn)

通過(guò)分析多光譜影像各單波段的紋理特征與水稻LAI、AGB、PNC間的相關(guān)性（圖4）發(fā)現(xiàn)，單波段紋理特征與水稻長(zhǎng)勢(shì)指標(biāo)間的相關(guān)性并未好于植被指數(shù)。在對(duì)LAI的監(jiān)測(cè)中，不同波段紋理特征與LAI的相關(guān)性有較大差異，表現(xiàn)較好的波段有近紅外和紅邊波段，表現(xiàn)較好的紋理特征有均值、方差、同質(zhì)性等，其中與LAI相關(guān)性最高的紋理特征為近紅外波段的均值，2=0.48。與AGB相關(guān)性較高的紋理特征同樣多集中在近紅外和紅邊波段，這與LAI保持一致，但在本研究中與AGB相關(guān)性最好的紋理特征為綠光波段的均值，2=0.43。相較于LAI和AGB，各個(gè)波段的紋理特征與PNC的相關(guān)性均較低，其中表現(xiàn)最好的紅光波段的二階矩，2僅為0.09。

圖3 基于植被指數(shù)建立的水稻葉面積指數(shù)、地上部生物量、植株氮含量監(jiān)測(cè)模型和檢驗(yàn)

單一的紋理特征在水稻LAI、AGB和PNC的估測(cè)中并未取得較好的效果，為消除系統(tǒng)誤差，嘗試將4個(gè)波段的各個(gè)紋理特征（共32種）兩兩配對(duì)，分別構(gòu)建了歸一化紋理指數(shù)（NDTI）、比值紋理指數(shù)（RTI）和差值紋理指數(shù)（DTI）并分析了紋理指數(shù)與水稻LAI、AGB和PNC間的相關(guān)性，表5列舉了6個(gè)表現(xiàn)較好的TI，其中與LAI相關(guān)性最高的TI為NDTI（MEA550，MEA790），2=0.59；與AGB相關(guān)性最高的TI同樣為NDTI（MEA550，MEA790），2=0.48；與PNC相關(guān)性最高的TI為DTI（SEC735，ENT660），2=0.26。本研究發(fā)現(xiàn)與LAI和AGB相關(guān)性較高的TI多是綠光波段與近紅外波段的均值特征的組合，而與PNC相關(guān)性較高的TI多含有紅光波段的紋理特征。相較于單一紋理，通過(guò)紋理間的組合計(jì)算得到的TI與水稻LAI、AGB和PNC間的相關(guān)性有了一定程度的提升，但相較于VI來(lái)說(shuō)，TI除在PNC的監(jiān)測(cè)方面好于VI，對(duì)于LAI、AGB的監(jiān)測(cè)效果并未得到提升。上述結(jié)果表明，紋理特征作為除光譜信息外重要的影像信息可以在水稻生長(zhǎng)指標(biāo)監(jiān)測(cè)過(guò)程中發(fā)揮獨(dú)有的優(yōu)勢(shì)，尤其是在水稻氮素營(yíng)養(yǎng)監(jiān)測(cè)方面。

圖4 紋理特征與水稻葉面積指數(shù)、地上部生物量、植株氮含量的相關(guān)性

表5 6個(gè)紋理指數(shù)與水稻葉面積指數(shù)、地上部生物量和植株氮含量的相關(guān)性

圖5 偏最小二乘回歸模型中各植被指數(shù)的VIP值

2.4 基于偏最小二乘回歸的長(zhǎng)勢(shì)監(jiān)測(cè)模型構(gòu)建與檢驗(yàn)

上述研究利用SLR的建模方法進(jìn)行了水稻LAI、AGB和PNC監(jiān)測(cè)模型的構(gòu)建。為進(jìn)一步探究光譜紋理結(jié)合對(duì)水稻長(zhǎng)勢(shì)指標(biāo)監(jiān)測(cè)精度的提升效果，首先將本文選定的6個(gè)VI作為輸入?yún)?shù)進(jìn)行了水稻LAI、AGB和PNC的PLSR模型構(gòu)建。建模結(jié)果表明，LAI、AGB和PNC的PLSR模型在提取兩個(gè)主成分后對(duì)上述指標(biāo)達(dá)到了較好的預(yù)測(cè)效果。在變量因子重要性排序中，VI的VIP值均＞0.8，表明VI對(duì)3個(gè)指標(biāo)的解釋作用都較好，CIRE、NDRE、RESAVI在VIP值排序中均在前3位，表明在偏最小二乘的建模方法下，這3個(gè)VI對(duì)模型的估測(cè)效果影響最大（圖5）?；赩I的LAI偏最小二乘估測(cè)模型驗(yàn)證2=0.76，=0.97，相比基于植被指數(shù)CIRE建立的LAI簡(jiǎn)單線性回歸模型驗(yàn)證2=0.75，=0.99，2略有提升，有一定程度的降低?；赩I建立的水稻AGB和PNC偏最小二乘估測(cè)模型驗(yàn)證2分別為0.73和0.27，分別為3.04和0.55，相比基于植被指數(shù)NDRE建立的AGB簡(jiǎn)單線性回歸模型驗(yàn)證2=0.72，=3.12和基于植被指數(shù)RESAVI建立的PNC簡(jiǎn)單線性回歸模型驗(yàn)證2=0.26，=0.55，2和均未有顯著改善（圖6）。

上述研究發(fā)現(xiàn)基于VI的PLSR模型對(duì)LAI、AGB、PNC的預(yù)測(cè)精度提升效果不顯著，因此嘗試將TI納入PLSR模型的構(gòu)建中以提升模型預(yù)測(cè)效果。選擇上述與LAI、AGB和PNC相關(guān)性較好的6個(gè)TI，加入TI后，各指數(shù)的VIP值如圖7所示。由VIP值排序可以發(fā)現(xiàn)，VI對(duì)LAI和AGB有較好的解釋性，其中表現(xiàn)較好的VI為CIRE、NDRE和RESAVI，這與SLR建模時(shí)的結(jié)果保持一致，TI中表現(xiàn)較好的有RTI（MEA550，MEA790）、NDTI（MEA550，MEA790）。而對(duì)于PNC，紋理指數(shù)對(duì)其解釋效果要好于植被指數(shù)，其中表現(xiàn)較好的TI有DTI（COR550，SEC660）和RTI（VAR660，CON660）。與單獨(dú)使用VI進(jìn)行PLSR模型構(gòu)建相比，基于VI和TI的PLSR模型在對(duì)LAI的估測(cè)上，驗(yàn)證2=0.77提升1%，=0.96降低1.03%。在對(duì)AGB估測(cè)方面，與僅使用VI相比，模型的驗(yàn)證2=0.81提升7%，=2.52降低17.11%，有效改善了模型的預(yù)測(cè)精度。在對(duì)PNC估測(cè)方面，模型的驗(yàn)證2=0.36提升9%，=0.48降低9.40%（圖8）。通過(guò)比較發(fā)現(xiàn)，當(dāng)TI參與到PLSR模型的構(gòu)建后，模型對(duì)LAI、AGB和PNC的預(yù)測(cè)精度均得到一定程度的提升。

圖7 偏最小二乘回歸模型中各植被指數(shù)和紋理指數(shù)的VIP值

2.5 基于人工神經(jīng)網(wǎng)絡(luò)的水稻長(zhǎng)勢(shì)監(jiān)測(cè)模型構(gòu)建與檢驗(yàn)

為比較不同模型對(duì)水稻長(zhǎng)勢(shì)指標(biāo)的反演效果，使用ANN進(jìn)行水稻LAI、AGB和PNC的估測(cè)。為保證與PLSR的一致性，首先以6個(gè)VI作為輸入?yún)?shù)，建立ANN預(yù)測(cè)模型。模型的隱含層數(shù)設(shè)置為1，考慮到模型的運(yùn)算效率，將模型神經(jīng)元數(shù)量設(shè)置為3—30進(jìn)行測(cè)試，當(dāng)模型值最低時(shí)確定最佳的神經(jīng)元數(shù)量。圖9顯示在ANN模型中，植被指數(shù)的重要性排序。將VI作為唯一輸入?yún)?shù)進(jìn)行LAI預(yù)測(cè)模型構(gòu)建，當(dāng)神經(jīng)元數(shù)量設(shè)置為23個(gè)時(shí)，模型的取得最小值0.74，此時(shí)模型的驗(yàn)證2=0.85，與PLSR模型相比，模型2提升9%且降低23.71%。在AGB的ANN預(yù)測(cè)模型中，當(dāng)神經(jīng)元數(shù)量設(shè)置為13個(gè)時(shí)，模型的取得最小值1.86，此時(shí)模型驗(yàn)證2=0.88，相較于PLSR模型，模型精度提升約15%且降低38.82%，模型的預(yù)測(cè)精度大幅提升。當(dāng)神經(jīng)元設(shè)置為8個(gè)時(shí)，對(duì)于PNC的預(yù)測(cè)取得最小值0.46，此時(shí)模型預(yù)測(cè)精度為53%，相較于PLSR，模型的預(yù)測(cè)精度提升約48%，降低13.21%（圖10）。

圖8 基于植被指數(shù)和紋理指數(shù)的偏最小二乘回歸模型中LAI、AGB、PNC驗(yàn)證結(jié)果

圖9 人工神經(jīng)網(wǎng)絡(luò)模型中各植被指數(shù)的VIP值

為進(jìn)一步探索不同模型在多元參數(shù)輸入下對(duì)水稻LAI、AGB和PNC的預(yù)測(cè)精度，將TI納入到ANN預(yù)測(cè)模型中。同樣設(shè)置隱含層數(shù)為1，神經(jīng)元數(shù)量在3—30間進(jìn)行調(diào)試。圖11顯示當(dāng)水稻LAI、AGB和PNC預(yù)測(cè)模型達(dá)到最小時(shí)，各輸入變量的因子重要性排序。當(dāng)神經(jīng)元數(shù)量設(shè)置為25個(gè)時(shí)，模型對(duì)LAI的預(yù)測(cè)效果達(dá)到最優(yōu)，此時(shí)=0.72，驗(yàn)證2=0.86。相較于僅基于VI的ANN預(yù)測(cè)模型，加入TI后模型的精度有一定程度的提升。當(dāng)隱含層神經(jīng)元數(shù)量設(shè)置為29個(gè)時(shí)，模型對(duì)AGB的預(yù)測(cè)效果達(dá)到最優(yōu)，此時(shí)=1.56，驗(yàn)證2=0.92。相較于僅基于VI的ANN預(yù)測(cè)模型，模型的降低不明顯，但模型預(yù)測(cè)精度有較大的提升。對(duì)于PNC，當(dāng)神經(jīng)元數(shù)量設(shè)置為30個(gè)時(shí)，模型的取得最小值0.25，此時(shí)模型的驗(yàn)證2=0.86，相較于僅基于VI的ANN預(yù)測(cè)模型，降低45.65%，模型的預(yù)測(cè)誤差顯著降低（圖12）。綜上可以發(fā)現(xiàn)，當(dāng)在ANN模型中加入紋理信息后模型的精度和預(yù)測(cè)誤差均有一定程度的改善，因此紋理信息可以作為一種重要的信息應(yīng)用到作物生長(zhǎng)監(jiān)測(cè)中。

圖10 基于植被指數(shù)的人工神經(jīng)網(wǎng)絡(luò)模型中LAI、AGB、PNC驗(yàn)證結(jié)果

圖11 人工神經(jīng)網(wǎng)絡(luò)模型中各植被指數(shù)和紋理指數(shù)的VIP值

上述研究結(jié)果表明，當(dāng)VI或TI作為單一參數(shù)進(jìn)行SLR建模時(shí)，模型對(duì)水稻LAI、AGB和PNC的監(jiān)測(cè)效果最差，當(dāng)采用PLSR和ANN將多個(gè)VI進(jìn)行監(jiān)測(cè)模型構(gòu)建時(shí)，模型的精度和均方根誤差相較于SLR建模均有一定程度的提升，而當(dāng)TI與VI共同參與監(jiān)測(cè)模型的構(gòu)建時(shí)，模型對(duì)水稻LAI、AGB以及PNC的監(jiān)測(cè)效果最好，通過(guò)模型間的比較發(fā)現(xiàn)，以VI+TI作為輸入?yún)?shù)，采用ANN的模型構(gòu)建方法對(duì)水稻LAI、AGB和PNC的監(jiān)測(cè)取得了最優(yōu)的效果。

圖12 基于植被指數(shù)和紋理指數(shù)的人工神經(jīng)網(wǎng)絡(luò)模型LAI、AGB、PNC驗(yàn)證結(jié)果

Fig. 12 Validation results of LAI, AGB and PNC in ANN model based on vegetation index and texture index

3 討論

3.1 固定翼無(wú)人機(jī)在水稻長(zhǎng)勢(shì)監(jiān)測(cè)中的優(yōu)勢(shì)

與以往研究多使用旋翼無(wú)人機(jī)進(jìn)行影像數(shù)據(jù)采集不同，本研究中使用了固定翼無(wú)人機(jī)作為機(jī)載平臺(tái)，相較于旋翼無(wú)人機(jī)，固定翼無(wú)人機(jī)具有飛行速度快，續(xù)航時(shí)間久等優(yōu)點(diǎn)。在本研究中，當(dāng)飛行高度設(shè)為80 m時(shí)，固定翼無(wú)人機(jī)搭載多光譜相機(jī)單次平均作業(yè)時(shí)長(zhǎng)可達(dá)45 min，作業(yè)面積可達(dá)90 hm2，相較于大疆精靈4多光譜版27 min、大疆經(jīng)緯M600 Pro 38 min的作業(yè)時(shí)長(zhǎng)和作業(yè)效率有顯著提升。本研究分析了基于固定翼無(wú)人機(jī)多光譜影像信息在水稻長(zhǎng)勢(shì)監(jiān)測(cè)中的效果，結(jié)果表明基于光譜信息構(gòu)建的植被指數(shù)和基于紋理特征構(gòu)建的紋理指數(shù)與水稻長(zhǎng)勢(shì)指標(biāo)間均存在相關(guān)性，這與預(yù)期結(jié)果一致。水稻內(nèi)部所含的色素水平、水分含量及其他結(jié)構(gòu)控制著水稻特殊的光譜響應(yīng)。而隨著作物生長(zhǎng)發(fā)育的進(jìn)行，多光譜影像的紋理特征信息顯現(xiàn)出獨(dú)有的優(yōu)勢(shì)，紋理特征對(duì)作物發(fā)育后期植株的形態(tài)變化更為敏感[31]。光譜與空間維度信息的結(jié)合可以顯著提升水稻長(zhǎng)勢(shì)監(jiān)測(cè)的精度，應(yīng)用固定翼無(wú)人機(jī)搭載多光譜相機(jī)以較低的作業(yè)成本和高效的作業(yè)效率實(shí)現(xiàn)了對(duì)水稻長(zhǎng)勢(shì)指標(biāo)信息的準(zhǔn)確獲取，為實(shí)現(xiàn)大面積的農(nóng)田精確管理提供了方法。

3.2 單一影像信息在水稻長(zhǎng)勢(shì)指標(biāo)監(jiān)測(cè)中的表現(xiàn)

Zheng等[32]利用多旋翼無(wú)人機(jī)搭載六波段多光譜相機(jī)采集了水稻關(guān)鍵生育期冠層影像，通過(guò)構(gòu)建植被指數(shù)和紋理指數(shù)實(shí)現(xiàn)了對(duì)水稻AGB的估算，研究結(jié)果表明優(yōu)化的土壤調(diào)整植被指數(shù)（OSAVI）與水稻AGB的相關(guān)性最高（2=0.63），而歸一化紋理指數(shù)NDTI（MEA800，MEA550）與水稻AGB的建模2達(dá)到0.75，優(yōu)于所有植被指數(shù)。同時(shí)Zheng等[33]還發(fā)現(xiàn)基于紅邊條帶和近紅外條帶的紋理特征組成的NDTI在水稻葉片氮含量和植株氮含量的監(jiān)測(cè)中表現(xiàn)較好，表現(xiàn)最好的紋理指數(shù)為NDTI（DIS800，CON720）和NDTI（COR800，DIS720），建模2分別為0.31和0.41。本研究基于植被指數(shù)構(gòu)建的水稻LAI、AGB和PNC簡(jiǎn)單線性回歸模型表明與水稻LAI、AGB和PNC相關(guān)性最好的植被指數(shù)為CIRE、NDRE和RESAVI，決定系數(shù)分別為0.80、0.76和0.13。3個(gè)植被指數(shù)中均含有紅邊和近紅外波段，由此可以判斷這兩個(gè)波段區(qū)間在水稻長(zhǎng)勢(shì)監(jiān)測(cè)中具有獨(dú)特的作用，這一結(jié)果與前人研究一致[34]。本研究同時(shí)發(fā)現(xiàn)單波段紋理特征對(duì)3個(gè)指標(biāo)的監(jiān)測(cè)效果遠(yuǎn)低于植被指數(shù)，造成此結(jié)果的原因可能是紋理信息易受尺度縮放和視點(diǎn)變化的影響，導(dǎo)致提取的紋理信息過(guò)于冗余，難以區(qū)分有效信息和噪聲。在不同施氮水平下的水稻無(wú)人機(jī)圖像具有顯著的變化特征，而紋理特征均值表示了灰度共生矩陣中所有元素的平均值，均值越大，表示圖像中的紋理越明顯，反之則表示圖像中的紋理越平滑。因此均值紋理特征可以反映出不同處理下水稻冠層的結(jié)構(gòu)差異。雖然大部分紋理特征與水稻LAI和AGB的相關(guān)性較低，但是mean紋理特征單獨(dú)使用或者利用數(shù)學(xué)組合的形式結(jié)合在水稻LAI和AGB監(jiān)測(cè)上取得了較好效果。本文參照植被指數(shù)構(gòu)建方法構(gòu)建了紋理指數(shù)，相較于單波段紋理特征，紋理指數(shù)與水稻LAI、AGB和PNC間的相關(guān)性均有所提升，在LAI和AGB的監(jiān)測(cè)中表現(xiàn)最好的紋理指數(shù)為NDTI（MEA550，MEA790），建模2分別為0.59和0.48，表現(xiàn)較好的紋理指數(shù)均由綠光波段和近紅外波段組成，綠光波段主要反映了作物葉片色素的信息，而近紅外波段主要體現(xiàn)了植物的細(xì)胞構(gòu)造及內(nèi)部化學(xué)成分信息，兩個(gè)波段均值特征的組合一定程度上消除了土壤及其他背景因素的影響，增強(qiáng)了水稻植被信息，為水稻LAI和AGB的估算提供了新的參數(shù)，這一結(jié)論在前人研究中也有所體現(xiàn)[35]。在對(duì)水稻PNC的監(jiān)測(cè)中紋理指數(shù)的效果普遍好于單波段紋理特征和植被指數(shù)，Clark等[36]研究表明，紅邊區(qū)間與植被葉綠素、氮素、以及葉面積指數(shù)等結(jié)構(gòu)特征有重要相關(guān)性，而在本研究中與水稻PNC相關(guān)性較好的紋理指數(shù)也多由紅邊條帶組成，表現(xiàn)最好的紋理指數(shù)DTI（SEC735，ENT660）與PNC的相關(guān)性為0.26。實(shí)際上，目前對(duì)于紋理特征的篩選依然缺少一套標(biāo)準(zhǔn)，篩選出的紋理特征很難解釋其原理；另外不同紋理算法、窗口大小和光譜波段均會(huì)影響最佳紋理特征的選擇，有待今后進(jìn)一步探討。

3.3 光譜信息與紋理特征相結(jié)合在水稻長(zhǎng)勢(shì)指標(biāo)監(jiān)測(cè)中的效果

利用機(jī)器學(xué)習(xí)的方法將影像光譜、紋理、顏色等信息融合實(shí)現(xiàn)對(duì)作物長(zhǎng)勢(shì)的監(jiān)測(cè)已成為當(dāng)前研究熱點(diǎn)。本研究采用PLSR和ANN兩種機(jī)器學(xué)習(xí)方法，無(wú)論采用單一指標(biāo)（VI）或融合指標(biāo)（VI+TI）作為模型輸入?yún)?shù)，基于ANN構(gòu)建的模型監(jiān)測(cè)效果都要好于PLSR，這在前人研究中也有所體現(xiàn)[37]。紋理作為一種常見(jiàn)的視覺(jué)現(xiàn)象，反映了灰度圖像中含有強(qiáng)弱關(guān)系的相似圖形在圖像中的局部結(jié)構(gòu)或排列規(guī)則[38]，灰度共生矩陣提供了影像灰度方向、間隔、變化幅度與快慢的信息，利用灰度共生矩陣提取的水稻冠層紋理信息有效的反映了其空間分布特征，該信息能夠有效緩解水稻生長(zhǎng)后期冠層光譜飽和帶來(lái)的監(jiān)測(cè)精度下降問(wèn)題。通過(guò)將紋理特征值進(jìn)行歸一化、比值、差值處理，可以消除影像中因地形、陰影等帶來(lái)的影響[39]，進(jìn)一步突出了地物特征，因此多指標(biāo)結(jié)合的水稻長(zhǎng)勢(shì)監(jiān)測(cè)模型綜合了光譜信息與紋理特征對(duì)水稻LAI、AGB和PNC估算的共同貢獻(xiàn)，其模型預(yù)測(cè)精度要好于僅基于光譜所建立的模型，該結(jié)果也驗(yàn)證了前人關(guān)于光譜紋理融合指標(biāo)可以提升葉面積指數(shù)[17]、生物量[40]和植株氮含量[41]監(jiān)測(cè)精度的研究結(jié)論。在本研究中，植被指數(shù)和紋理指數(shù)對(duì)水稻LAI和AGB的監(jiān)測(cè)效果都要好于PNC，原因是在水稻生長(zhǎng)期間，LAI和AGB變化很大，這種變化遮蓋了葉綠素和氮素對(duì)冠層反射率和空間特征的貢獻(xiàn)，因此植被指數(shù)和紋理指數(shù)對(duì)水稻LAI和AGB的變化更為敏感。而在對(duì)水稻氮素含量遙感監(jiān)測(cè)的過(guò)程中，僅獲取了作物上層葉片光譜紋理信息，由此也導(dǎo)致了圖譜信息對(duì)植株整體氮素含量的反演效果并不理想。本研究中大部分紋理特征與植株氮含量的相關(guān)性很低，而通過(guò)構(gòu)建紋理指數(shù)的方式有效提升了對(duì)水稻PNC的監(jiān)測(cè)效果，但是相較于LAI和AGB，其監(jiān)測(cè)效果也并不理想。在本研究中，基于植被指數(shù)和紋理指數(shù)結(jié)合建立的人工神經(jīng)網(wǎng)絡(luò)模型在水稻LAI、AGB和PNC的監(jiān)測(cè)中取得了最好的效果，模型的驗(yàn)證2分別為0.86、0.92和0.86，同時(shí)與其他建模方法相比也達(dá)到最低值。通過(guò)對(duì)模型變量投影重要性分析，具有紅邊波段的植被指數(shù)對(duì)水稻LAI和AGB的反演具有較大貢獻(xiàn)，而對(duì)于水稻PNC，在不同建模方法中植被指數(shù)與紋理指數(shù)的重要性不同，因此尚無(wú)法得出統(tǒng)一結(jié)論。

本研究仍然存在一些不足之處，如已有研究表明在不同地面分辨率下，影像的紋理特征與作物生長(zhǎng)指標(biāo)間的相關(guān)性存在較大差異，本文并未討論在不同分辨率下紋理特征對(duì)模型精度的影響。同時(shí)，本研究尚未嘗試反射率與紋理特征直接的數(shù)學(xué)組合方式，“圖譜”指數(shù)對(duì)作物水稻長(zhǎng)勢(shì)指標(biāo)的監(jiān)測(cè)效果尚未明確。再則，由于取樣間隔時(shí)間長(zhǎng)，樣本數(shù)據(jù)較少，對(duì)于模型魯棒性有一定影響，后續(xù)可以增加樣本量、作物品種、不同年限以及更多地區(qū)等來(lái)驗(yàn)證模型，進(jìn)而提升模型普適性及魯棒性。

4 結(jié)論

基于固定翼無(wú)人機(jī)平臺(tái)搭載多光譜傳感器獲取水稻冠層關(guān)鍵生育期影像數(shù)據(jù)，通過(guò)無(wú)人機(jī)影像光譜和紋理信息進(jìn)行了植被指數(shù)與紋理指數(shù)的構(gòu)建。相較于單波段紋理特征，紋理指數(shù)與水稻LAI、AGB和PNC的相關(guān)性更好，但與植被指數(shù)相比，紋理指數(shù)并未取得更好的監(jiān)測(cè)效果，因此影像紋理只能作為光譜的輔助信息應(yīng)用到作物長(zhǎng)勢(shì)監(jiān)測(cè)中。進(jìn)一步采用SLR、PLSR和ANN 3種方法建立了水稻長(zhǎng)勢(shì)無(wú)損監(jiān)測(cè)模型，植被指數(shù)與紋理指數(shù)結(jié)合后的預(yù)測(cè)模型效果要好于單一變量。而基于PLSR和ANN兩種機(jī)器學(xué)習(xí)的建模效果要優(yōu)于SLR的建模效果，表現(xiàn)最好的建模方法為ANN。以植被指數(shù)和紋理指數(shù)為輸入?yún)?shù)，采用ANN的模型構(gòu)建方法，在對(duì)LAI的監(jiān)測(cè)中，模型的驗(yàn)證精度達(dá)到了86%，而在AGB和PNC的監(jiān)測(cè)中，模型的驗(yàn)證精度分別達(dá)到了92%和86%。研究結(jié)果證明了基于無(wú)人機(jī)多光譜影像數(shù)據(jù)實(shí)現(xiàn)對(duì)水稻長(zhǎng)勢(shì)和氮素狀況監(jiān)測(cè)的可行性，為作物長(zhǎng)勢(shì)快速監(jiān)測(cè)與精確管理提供了技術(shù)途徑。

[1] 陳明江, 劉貴富, 余泓, 王冰, 李家洋. 水稻高產(chǎn)優(yōu)質(zhì)的分子基礎(chǔ)與品種設(shè)計(jì). 科學(xué)通報(bào), 2018, 63(14): 1275-1289.

CHEN M J, LIU G F, YU H, WANG B, LI J Y. Towards molecular design of rice plant architecture and grain quality. Chinese Science Bulletin, 2018, 63(14): 1275-1289. (in Chinese)

[2] JIANG J, ATKINSON P, ZHANG J Y, LU R H, ZHOU Y Y, CAO Q, TIAN Y C, ZHU Y, CAO W X, LIU X J. Combining fixed-wing UAV multispectral imagery and machine learning to diagnose winter wheat nitrogen status at the farm scale. European Journal of Agronomy, 2022, 138: 126537.

[3] ZHAO B, ATA-UL-KARIM S T, DUAN A W, LIU Z D, WANG X L, XIAO J F, LIU Z G, QIN A Z, NING D F, ZHANG W Q, LIAN Y H. Determination of critical nitrogen concentration and dilution curve based on leaf area index for summer maize. Field Crops Research, 2018, 228: 195-203.

[4] ATA-UL-KARIM S T, ZHU Y, YAO X, CAO W X. Determination of critical nitrogen dilution curve based on leaf area index in rice. Field Crops Research, 2014, 167: 76-85.

[5] WANG X B, MIAO Y X, DONG R, ZHA H N, XIA T T, CHEN Z C, KUSNIEREK K, MI G H, SUN H, LI M Z. Machine learning-based in-season nitrogen status diagnosis and side-dress nitrogen recommendation for corn. European Journal of Agronomy, 2021, 123: 126193.

[6] BASCON M V, NAKATA T, SHIBATA S, TAKATA I, KOBAYASHI N, KATO Y, INOUE S, DOI K, MURASE J, NISHIUCHI S. Estimating yield-related traits using UAV-derived multispectral images to improve rice grain yield prediction. Agriculture, 2022, 12(8): 1141.

[7] FU Z P, JIANG J, GAO Y, KRIENKE B, WANG M, ZHONG K, CAO Q, TIAN Y C, ZHU Y, CAO W X, LIU X J. Wheat growth monitoring and yield estimation based on multi-rotor unmanned aerial vehicle. Remote Sensing, 2020, 12(3): 508.

[8] TAO H L, FENG H K, XU L J, MIAO M K, YANG G J, YANG X D, FAN L L. Estimation of the yield and plant height of winter wheat using UAV-based hyperspectral images. Sensors, 2020, 20(4): 1231.

[9] LU N, ZHOU J, HAN Z X, LI D, CAO Q, YAO X, TIAN Y C, ZHU Y, CAO W X, CHENG T. Improved estimation of aboveground biomass in wheat from RGB imagery and point cloud data acquired with a low-cost unmanned aerial vehicle system. Plant Methods, 2019, 15: 17.

[10] 陳鵬飛, TREMBLAY N, 王紀(jì)華, VIGNEAULT P, 黃文江, 李保國(guó). 估測(cè)作物冠層生物量的新植被指數(shù)的研究. 光譜學(xué)與光譜分析, 2010, 30(2): 512-517.

CHEN P F, TREMBLAY N, WANG J H, VIGNEAULT P, HUANG W J, LI B G. New index for crop canopy fresh biomass estimation. Spectroscopy and Spectral Analysis, 2010, 30(2): 512-517. (in Chinese)

[11] YANG G J, LIU J G, ZHAO C J, LI Z H, HUANG Y B, YU H Y, XU B, YANG X D, ZHU D M, ZHANG X Y, ZHANG R Y, FENG H K, ZHAO X Q, LI Z H, LI H L, YANG H. Unmanned aerial vehicle remote sensing for field-based crop phenotyping: current status and perspectives. Frontiers in Plant Science, 2017, 8: 1111.

[12] YANG K L, GONG Y, FANG S H, DUAN B, YUAN N G, PENG Y, WU X T, ZHU R S. Combining spectral and texture features of UAV images for the remote estimation of rice LAI throughout the entire growing season. Remote Sensing, 2021, 13(15): 3001.

[13] MAIMAITIJIANG M, SAGAN V, SIDIKE P, HARTLING S, ESPOSITO F, FRITSCHI F B. Soybean yield prediction from UAV using multimodal data fusion and deep learning. Remote Sensing of Environment, 2020, 237: 111599.

[14] 萬(wàn)亮, 杜曉月, 陳碩博, 于豐華, 朱姜蓬, 許童羽, 何勇, 岑海燕. 基于無(wú)人機(jī)多源圖譜融合的水稻稻穗表型監(jiān)測(cè). 農(nóng)業(yè)工程學(xué)報(bào), 2022, 38(9): 162-170.

WAN L, DU X Y, CHEN S B, YU F H, ZHU J P, XU T Y, HE Y, CEN H Y. Rice panicle phenotyping using UAV-based multi-source spectral image data fusion. Transactions of the Chinese Society of Agricultural Engineering, 2022, 38(9): 162-170. (in Chinese)

[15] 戴震. 基于無(wú)人機(jī)多源影像信息融合的玉米冠層氮含量遙感監(jiān)測(cè)研究[D]. 合肥: 安徽大學(xué), 2021.

DAI Z. Nitrogen content monitoring of corn canopy based on UAV multi-source data fusion[D]. Hefei: Anhui University, 2021. (in Chinese)

[16] ZHENG L, TAO J P, BAO Q, WENG S Z, ZHANG Y K, ZHAO J L, HUANG L S. Combining spectral and textures of digital imagery for wheat aboveground biomass estimation//International Conference on Electronic Information Technology (EIT 2022): SPIE, 2022: 304-307.

[17] WULDER M A, LEDREW E F, FRANKLIN S E, LAVIGNE M B. Aerial image texture information in the estimation of northern deciduous and mixed wood forest leaf area index (LAI). Remote Sensing of Environment, 1998, 64(1): 64-76.

[18] 劉俊, 畢華興, 朱沛林, 孫菁, 朱金兆, 陳濤. 基于ALOS遙感數(shù)據(jù)紋理及紋理指數(shù)的柞樹(shù)蓄積量估測(cè). 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2014, 45(7): 245-254.

LIU J, BI H X, ZHU P L, SUN J, ZHU J Z, CHEN T. Estimating stand volume offorest based on texture parameters and derivative texture indices of alos imagery. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(7): 245-254. (in Chinese)

[19] 陳鵬飛, 梁飛. 基于低空無(wú)人機(jī)影像光譜和紋理特征的棉花氮素營(yíng)養(yǎng)診斷研究. 中國(guó)農(nóng)業(yè)科學(xué), 2019, 52(13): 2220-2229.doi: 10.3864/j.issn.0571-1752.2019.13.003.

CHEN P F, LIANG F. Cotton nitrogen nutrition diagnosis based on spectrum and texture feature of images from low altitude unmanned aerial vehicle. Scientia Agricultura Sinica, 2019, 52(13): 2220-2229. doi: 10.3864/j.issn.0571-1752.2019.13.003. (in Chinese)

[20] 陳鵬, 馮海寬, 李長(zhǎng)春, 楊貴軍, 楊鈞森, 楊文攀, 劉帥兵. 無(wú)人機(jī)影像光譜和紋理融合信息估算馬鈴薯葉片葉綠素含量. 農(nóng)業(yè)工程學(xué)報(bào), 2019, 35(11): 63-74.

CHEN P, FENG H K, LI C C, YANG G J, YANG J S, YANG W P, LIU S B. Estimation of chlorophyll content in potato using fusion of texture and spectral features derived from UAV multispectral image. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(11): 63-74. (in Chinese)

[21] WANG Y Y, ZHANG K, TANG C L, CAO Q, TIAN Y C, ZHU Y, CAO W X, LIU X J. Estimation of rice growth parameters based on linear mixed-effect model using multispectral images from fixed-wing unmanned aerial vehicles. Remote Sensing, 2019, 11(11): 1371.

[22] LIU X H, WANG L. Feasibility of using consumer-grade unmanned aerial vehicles to estimate leaf area index in mangrove forest. Remote Sensing Letters, 2018, 9(11): 1040-1049.

[23] ZHANG J Y, QIU X L, WU Y T, ZHU Y, CAO Q, LIU X J, CAO W X. Combining texture, color, and vegetation indices from fixed-wing UAS imagery to estimate wheat growth parameters using multivariate regression methods. Computers and Electronics in Agriculture, 2021, 185: 106138.

[24] 田慶久, 閔祥軍. 植被指數(shù)研究進(jìn)展. 地球科學(xué)進(jìn)展, 1998, 13(4): 327-333.

TIAN Q J, MIN X J. Advances in study on vegetation indices. Advance in Earth Sciences, 1998, 13(4): 327-333. (in Chinese)

[25] 楊福芹, 馮海寬, 肖天豪, 李天馳, 郭向前. 融合無(wú)人機(jī)影像光譜與紋理特征的冬小麥氮營(yíng)養(yǎng)指數(shù)估算. 農(nóng)業(yè)現(xiàn)代化研究, 2020, 41(4): 718-726.

YANG F Q, FENG H K, XIAO T H, LI T C, GUO X Q. Nitrogen nutrition index estimation in winter wheat by UAV spectral information and texture feature fusion. Research of Agricultural Modernization, 2020, 41(4): 718-726. (in Chinese)

[26] 杭艷紅, 蘇歡, 于滋洋, 劉煥軍, 官海翔, 孔繁昌. 結(jié)合無(wú)人機(jī)光譜與紋理特征和覆蓋度的水稻葉面積指數(shù)估算. 農(nóng)業(yè)工程學(xué)報(bào), 2021, 37(9): 64-71.

HANG Y H, SU H, YU Z Y, LIU H J, GUAN H X, KONG F C. Estimation of rice leaf area index combining UAV spectrum, texture features and vegetation coverage. Transactions of the Chinese Society of Agricultural Engineering, 2021, 37(9): 64-71. (in Chinese)

[27] 李嵐?jié)? 郭宇龍, 韓鵬, 王丹丹, 王宜倫. 基于高光譜的冬小麥不同生育時(shí)期地上部生物量監(jiān)測(cè). 麥類作物學(xué)報(bào), 2021, 41(7): 904-913.

LI L T, GUO Y L, HAN P, WANG D D, WANG Y L. Estimation of shoot biomass at different growth stages of winter wheat based on hyperspectral reflectance. Journal of Triticeae Crops, 2021, 41(7): 904-913. (in Chinese)

[28] YUAN H H, YANG G J, LI C C, WANG Y J, LIU J G, YU H Y, FENG H K, XU B, ZHAO X Q, YANG X D. Retrieving soybean leaf area index from unmanned aerial vehicle hyperspectral remote sensing: Analysis of RF, ANN, and SVM regression models. Remote Sensing, 2017, 9(4): 309.

[29] 趙曉陽(yáng), 張建, 張東彥, 周新根, 劉小輝, 謝靜. 低空遙感平臺(tái)下可見(jiàn)光與多光譜傳感器在水稻紋枯病病害評(píng)估中的效果對(duì)比研究. 光譜學(xué)與光譜分析, 2019, 39(4): 1192-1198.

ZHAO X Y, ZHANG J, ZHANG D Y, ZHOU X G, LIU X H, XIE J. Comparison between the effects of visible light and multispectral sensor based on low-altitude remote sensing platform in the evaluation of rice sheath blight. Spectroscopy and Spectral Analysis, 2019, 39(4): 1192-1198. (in Chinese)

[30] LIU Y N, LIU S S, LI J, GUO X Y, WANG S Q, LU J W. Estimating biomass of winter oilseed rape using vegetation indices and texture metrics derived from UAV multispectral images. Computers and Electronics in Agriculture, 2019, 166: 105026.

[31] YUE J B, YANG G J, TIAN Q J, FENG H K, XU K J, ZHOU C Q. Estimate of winter-wheat above-ground biomass based on UAV ultrahigh-ground-resolution image textures and vegetation indices. ISPRS Journal of Photogrammetry and Remote Sensing, 2019, 150: 226-244.

[32] ZHENG H B, CHENG T, ZHOU M, LI D, YAO X, TIAN Y C, CAO W X, ZHU Y. Improved estimation of rice aboveground biomass combining textural and spectral analysis of UAV imagery. Precision Agriculture, 2019, 20(3): 611-629.

[33] ZHENG H B, ZHOU M, ZHU Y, CHENG T. Exploiting the textural information of UAV multispectral imagery to monitor nitrogen status in rice//IGARSS 2019-2019 IEEE International Geoscience and Remote Sensing Symposium, 2019: 7251-7253.

[34] ALI A, IMRAN M. Evaluating the potential of red edge position (REP) of hyperspectral remote sensing data for real time estimation of LAI & chlorophyll content of kinnow mandarin () fruit orchards. Scientia Horticulturae, 2020, 267: 109326.

[35] XU T Y, WANG F M, XIE L L, YAO X P, ZHENG J Y, LI J L, CHEN S T. Integrating the textural and spectral information of UAV hyperspectral images for the improved estimation of rice aboveground biomass. Remote Sensing, 2022, 14(11): 2534.

[36] Clark M L, Roberts D A, Clark D B. Hyperspectral discrimination of tropical rain forest tree species at leaf to crown scales. Remote sensing of environment, 2005, 96(3/4): 375-398.

[37] 劉秀英. 玉米生理參數(shù)及農(nóng)田土壤信息高光譜監(jiān)測(cè)模型研究[D]. 楊凌: 西北農(nóng)林科技大學(xué), 2016.

LIU X Y. Monitoring models of physiological parameters of corn and farmland soil information based on hyper-spectral reflectance[D]. Yangling: Northwest A&F University, 2016. (in Chinese)

[38] 龔紅菊, 姬長(zhǎng)英. 基于圖像處理技術(shù)的麥穗產(chǎn)量測(cè)量方法. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2007, 38(12): 116-119.

GONG H J, JI C Y. Estimating yield of wheat spike with its texture features based on image processing technology. Transactions of the Chinese Society for Agricultural Machinery, 2007, 38(12):116-119. (in Chinese)

[39] 鄭恒彪. 水稻生育期及生長(zhǎng)參數(shù)的近地面遙感監(jiān)測(cè)研究[D]. 南京: 南京農(nóng)業(yè)大學(xué), 2018.

ZHENG H B. Monitoring rice phenology and growth parameters using near-ground remote sensing platforms[D]. Nanjing: Nanjing Agricultural University, 2018. (in Chinese)

[40] LU D. Aboveground biomass estimation using Landsat TM data in the Brazilian Amazon. International Journal of Remote Sensing, 2005, 26(12): 2509-2525.

[41] 孫俊, 金夏明, 毛罕平, 武小紅, 朱文靜, 張曉東, 高洪燕. 基于高光譜圖像光譜與紋理信息的生菜氮素含量檢測(cè). 農(nóng)業(yè)工程學(xué)報(bào), 2014, 30(10): 167-173.

SUN J, JIN X M, MAO H P, WU X H, ZHU W J, ZHANG X D, GAO H Y. Detection of nitrogen content in lettuce leaves based on spectroscopy and texture using hyperspectral imaging technology. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(10): 167-173. (in Chinese)

Non-destructive monitoring of rice growth Key indicators based on fixed-wing UAV multispectral images

Wang WeiKang, Zhang JiaYi, Wang Hui, Cao Qiang, Tian YongChao, Zhu Yan, Cao WeiXing, Liu XiaoJun

College of Agriculture, Nanjing Agricultural University/National Engineering and Technology Center for Information Agriculture/ Engineering Research Center of Smart Agriculture, Ministry of Education/Key Laboratory for Crop System Analysis and Decision Making, Ministry of Agriculture and Rural Affairs/Jiangsu Key Laboratory for Information Agriculture, Nanjing 210095

【Background】In recent years, with the rapid development of remote sensing technology, real-time and non-destructive monitoring of crop growth status has become a research hotspot. Remote sensing-derived agricultural information will provide guidance for the precise management of large-scale crops. Among various remote sensing monitoring platforms, unmanned aerial vehicles (UAVs) have attracted wide attention due to their simple operation and low cost. UAVs equipped with multispectral cameras can quickly obtain crop growth conditions.【Objective】This study attempted to combine texture information and spectral information of multispectral images of fixed-wing UAVs to explore the monitoring effect of “atlas” information on rice growth indicators.【Method】A two-year rice field experiment involving different sowing dates, varieties, planting methods and nitrogen levels was conducted. During the key growth stages of rice, remote sensing images of the rice canopy were obtained using a Sequoia multispectral camera mounted on a fixed-wing UAV. Shoot destructive sampling was conducted simultaneously to obtain leaf area index (LAI), aboveground biomass (AGB), plant nitrogen content (PNC) and other agronomic indexes of rice. Simple regression, partial least squares regression and artificial neural network algorithms were used to construct rice growth index monitoring model based on multispectral images of fixed-wing UAV. The monitoring effects of spectral texture information in different models were compared and analyzed.【Result】The quantitative relationship between vegetation index (VI), single-band texture features and rice LAI, AGB, and PNC was explored using simple linear regression. The results showed that vegetation indexes had strong correlations with LAI and AGB, with the best-performing indexes being CIRE and NDRE, with2values of 0.80 and 0.76, respectively. However, for PNC monitoring, vegetation indexes did not achieve ideal results, with the best-performing RESAVI and NDRE having2values of only 0.13 with PNC. Further analysis using simple linear regression revealed that single-band texture features did not perform well in monitoring rice growth indicators. In order to further analyze the monitoring effect of image texture on the above three indexes, normalized texture indexes (NDTI), ratio texture indexes (RTI), and difference texture indexes (DTI) were constructed by referring to the construction method of VI. Correlation analysis showed that the newly constructed texture index (TI) improved the monitoring accuracy of rice growth indicators compared to single-band texture feature but did not perform better than vegetation indexes. To combine spectral and texture information, partial least squares and artificial neural network modeling methods were adopted in this paper. VI and VI+TI were used as different input parameter combinations to construct rice LAI, AGB and PNC monitoring models. The results showed that both partial least squares and artificial neural network modeling methods significantly improved the monitoring accuracy compared to simple linear regression. The best performance was achieved using VI+TI as input variables and an artificial neural network model for validation, with validation2values for LAI, AGB, and PNC models increasing from 0.75, 0.72, and 0.26 to 0.86, 0.92, and 0.86, respectively, while RMSE values were significantly reduced.【Conclusion】The monitoring accuracy of rice LAI, AGB and PNC can be effectively improved by using the fixed-wing UAV to collect multispectral images of rice canopy and using the texture features and reflectance information as input parameters of the model through the model construction method of artificial neural network. The research results will provide a theoretical basis for rapid monitoring of large area crop growth.

unmanned aerial vehicle (UAV); vegetation index (VI); texture feature; growth; rice

2023-02-23；

2023-05-10

國(guó)家自然科學(xué)基金（32071903）、國(guó)家重點(diǎn)研發(fā)計(jì)劃（2022YFD2301402）

王偉康，E-mail：2019801206@njau.edu.cn。通信作者劉小軍，E-mail：liuxj@njau.edu.cn

（責(zé)任編輯 岳梅）