基于無人機(jī)高光譜分?jǐn)?shù)階微分的馬鈴薯地上生物量估算

劉 楊 馮海寬 孫 乾 楊福芹 楊貴軍

(1.北京農(nóng)業(yè)信息技術(shù)研究中心農(nóng)業(yè)農(nóng)村部農(nóng)業(yè)遙感機(jī)理與定量遙感重點(diǎn)實(shí)驗(yàn)室， 北京 100097；2.山東科技大學(xué)測繪科學(xué)與工程學(xué)院， 青島 266590； 3.國家農(nóng)業(yè)信息化工程技術(shù)研究中心， 北京 100097；4.河南工程學(xué)院土木工程學(xué)院， 鄭州 451191； 5.北京市農(nóng)業(yè)物聯(lián)網(wǎng)工程技術(shù)研究中心， 北京 100097)

0 引言

及時掌握馬鈴薯的生長狀況對于指導(dǎo)田間管理和優(yōu)化種植格局、挖掘生產(chǎn)潛力具有重要意義。地上生物量(AGB)是反映作物生長狀況的重要指標(biāo)，AGB的變化能夠直接表征作物有機(jī)物的積累能力，反映作物生長狀況和營養(yǎng)狀態(tài)[1-2]?？焖贉?zhǔn)確地獲取AGB信息對于監(jiān)測馬鈴薯光合作用能力和生長狀態(tài)至關(guān)重要[3-6]。傳統(tǒng)的人工測量手段獲取AGB信息準(zhǔn)確，但破壞性極強(qiáng)，且難以滿足大規(guī)模種植的監(jiān)測需求。高光譜遙感技術(shù)具有較高的波譜分辨率和較強(qiáng)的波段連續(xù)性，是對地觀測的重大突破[7-10]。由于作物冠層對太陽光的吸收和反射形成特有的光譜曲線，故可通過分析作物的高光譜反射特性實(shí)現(xiàn)作物AGB的快速無損估測[11-12]。

無人機(jī)高光譜遙感具有機(jī)動性強(qiáng)、成本低、可云下獲取影像等優(yōu)勢，已成為當(dāng)前農(nóng)業(yè)領(lǐng)域所關(guān)注的焦點(diǎn)。國內(nèi)外學(xué)者基于無人機(jī)高光譜信息在作物地上生物量估算研究方面取得了一定的成果[13-19]。

現(xiàn)有研究多基于原始冠層光譜特征或利用其構(gòu)建的植被指數(shù)與地上生物量建立關(guān)系，構(gòu)建作物地上生物量估算模型。光譜微分技術(shù)如分?jǐn)?shù)階微分能夠細(xì)化光譜信息，消除環(huán)境背景的影響，深度挖掘光譜中潛在的信息。近年來，利用分?jǐn)?shù)階微分技術(shù)在各個領(lǐng)域得到了廣泛研究與應(yīng)用[20-23]。然而，分?jǐn)?shù)階微分技術(shù)在作物營養(yǎng)監(jiān)測方面的研究還較少。

本文對波段454～950 nm的馬鈴薯冠層無人機(jī)高光譜數(shù)據(jù)進(jìn)行0～2階(間隔0.2)微分處理，得到每個生育期11種冠層分?jǐn)?shù)階微分光譜數(shù)據(jù)。通過相關(guān)性分析，挑選出各生育期與馬鈴薯AGB相關(guān)性較好的各階分?jǐn)?shù)階微分光譜，計(jì)算其相關(guān)系數(shù)絕對值，優(yōu)選不同階微分波段，構(gòu)建馬鈴薯AGB的多元線性回歸模型、隨機(jī)森林模型和人工神經(jīng)網(wǎng)絡(luò)模型，并進(jìn)行模型精度驗(yàn)證，最后篩選出最優(yōu)估算模型。

1 試驗(yàn)設(shè)計(jì)與數(shù)據(jù)獲取

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

于2019年3—7月在北京市昌平區(qū)小湯山鎮(zhèn)國家精準(zhǔn)農(nóng)業(yè)研究示范基地的馬鈴薯試驗(yàn)田開展試驗(yàn)，位于北緯40°10′34″，東經(jīng)116°26′39″，平均海拔為36 m，氣候類型為暖溫帶半濕潤大陸性季風(fēng)氣候。如圖1所示，試驗(yàn)區(qū)共設(shè)密度試驗(yàn)(N區(qū))、氮素試驗(yàn)(S區(qū))、鉀肥試驗(yàn)(K區(qū))3個試驗(yàn)區(qū)，采用完全隨機(jī)試驗(yàn)設(shè)計(jì)，N區(qū)和S區(qū)均設(shè)2個試驗(yàn)品種，分別為中薯5(P1)和中薯3(P2)，均為早熟品種，K區(qū)設(shè)1個試驗(yàn)品種，為中薯3。密度試驗(yàn)設(shè)3個水平：60 000株/hm2(T1)、72 000株/hm2(T2)、84 000株/hm2(T3)，6個處理，每個處理重復(fù)3次，共18個試驗(yàn)小區(qū)；氮素試驗(yàn)設(shè)4個水平(以尿素計(jì)量)：0 kg/hm2(N0)、244.65 kg/hm2(N1)、489.15 kg/hm2(N2,正常處理，15 kg純氮)、733.5 kg/hm2(N3)，8個處理，每個處理重復(fù)3次，共24個試驗(yàn)小區(qū)；鉀肥試驗(yàn)設(shè)3個水平：0 kg/hm2(K0)、970.5 kg/hm2(K1，N區(qū)和S區(qū)均為K1處理)、1 941 kg/hm2(K2)，2個處理，每個處理重復(fù)3次，共6個試驗(yàn)小區(qū)；試驗(yàn)小區(qū)總計(jì)48個，試驗(yàn)小區(qū)面積為6.5 m×5 m。

圖1 馬鈴薯試驗(yàn)位置和試驗(yàn)設(shè)計(jì)Fig.1 Experimental area and design of potato

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

分別獲取馬鈴薯現(xiàn)蕾期(2019年5月13日)、塊莖形成期(2019年5月28日)、塊莖增長期(2019年6月10日)、淀粉積累期(2019年6月20日)、成熟期(2019年7月3日)5個關(guān)鍵時期的地上生物量數(shù)據(jù)。馬鈴薯地上生物量通過收獲法獲取，在每個小區(qū)中選取代表性的3顆植株，將其莖葉分離，隨后用清水洗凈，105℃殺青，80℃干燥48 h以上，直到質(zhì)量恒定再進(jìn)行稱量。將植株莖、葉的干質(zhì)量求和得到樣本干質(zhì)量，最后通過樣本干質(zhì)量與種植密度得到每個小區(qū)的馬鈴薯地上生物量[24-25]。

1.3 無人機(jī)高光譜數(shù)據(jù)獲取及處理

試驗(yàn)利用無人機(jī)搭載的德國Cubert公司生產(chǎn)的UHD185型機(jī)載成像光譜儀，該光譜儀尺寸為195 nm×67 nm×60 nm，質(zhì)量470 g，波譜范圍為450～950 nm，共有125個光譜通道，采樣間隔4 nm，光譜分辨率8 nm，數(shù)字分辨率12位。選擇晴朗無云的天氣進(jìn)行拍攝，此時太陽光照強(qiáng)度穩(wěn)定。飛行時間10:30—14:00，飛行高度為50 m。

無人機(jī)獲取的高光譜數(shù)據(jù)處理主要包括輻射校正、影像拼接、影像融合和光譜提取4部分。根據(jù)UHD185型成像光譜儀中心波長和波長半幅寬在Matlab環(huán)境設(shè)計(jì)的輻射定標(biāo)系統(tǒng)，完成由影像DN值到地表反射率的輻射校正。利用Agisoft PhotoScan 軟件進(jìn)行高光譜影像拼接，變?yōu)閏ue格式數(shù)據(jù)，再提取子波段為jpg格式，最后再對各個拼接的子波段影像進(jìn)行合并生成馬鈴薯試驗(yàn)田的高光譜影像數(shù)據(jù)?；贏rcmap10.2軟件，繪制馬鈴薯每個小區(qū)矢量數(shù)據(jù)，基于IDL語言統(tǒng)計(jì)每個小區(qū)平均光譜作為馬鈴薯冠層光譜，得到高光譜反射率數(shù)據(jù)。

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

分?jǐn)?shù)階微分是從整數(shù)階微分發(fā)展而來的，本質(zhì)上是任意階光譜反射率的斜率。目前，比較流行的分?jǐn)?shù)階微分定義為Grunwald-Letnikov、Riemann-Liouville和Caputo 3種形式[15]。3種定義中，Grunwald-Letnikov因形式簡單而被廣泛運(yùn)用，故本文運(yùn)用此形式完成對馬鈴薯冠層光譜信息的深度挖掘，該微分形式為

(1)

式中Γ——Gamma函數(shù)λ——對應(yīng)波長

n——微分上下限之差

a——任意階數(shù)f(λ)——光譜

當(dāng)a=0、1或2時，f(λ)為原始光譜、一階微分光譜或二階微分光譜；當(dāng)a是小數(shù)時，則為Grunwald-Letnikov微分分?jǐn)?shù)階形式。

1.5 精度評價

本文對每個生育期分別挑選2/3樣本數(shù)據(jù)作為建模集，1/3樣本數(shù)據(jù)作為驗(yàn)證集，以此構(gòu)建馬鈴薯AGB估算模型。采用決定系數(shù)(Coefficient of determination,R2)、均 方 根 誤 差(Root mean square error, RMSE)、標(biāo)準(zhǔn)均方根誤差(Normalized root mean square error, NRMSE)評價模型的精度。R2越接近于1，RMSE和NRMSE越低，其估測模型的精度就越高。

2 結(jié)果與分析

2.1 原始冠層光譜與生物量相關(guān)性

由圖2可知，馬鈴薯冠層光譜曲線呈現(xiàn)典型的綠色植物反射特征，在可見光波段范圍內(nèi)，存在2個吸收帶，分別位于480 nm和670 nm附近，在550 nm附近形成小的反射峰，這一特征變化主要由于馬鈴薯葉片中葉綠素對綠光反射作用強(qiáng)，對藍(lán)光和紅光吸收強(qiáng)造成的。在波段674～780 nm之間馬鈴薯冠層反射率急增，形成典型綠色植被特有的“紅邊”特征。

圖2 不同試驗(yàn)的馬鈴薯冠層光譜曲線Fig.2 Potato canopy spectral curves of different experiments

為篩選出與馬鈴薯AGB相關(guān)性較好的原始冠

層光譜波段，通過相關(guān)性分析，得到各個生育期原始冠層光譜與AGB的相關(guān)性，結(jié)果如圖3所示。從圖3可以看出：①現(xiàn)蕾期，冠層原始光譜在波段454～718 nm范圍內(nèi)與馬鈴薯AGB呈極顯著負(fù)相關(guān)(P<0.01)，在波段746～914 nm范圍內(nèi)與AGB呈極顯著正相關(guān)。由于與AGB相關(guān)的光譜波段主要是可見光波段，故選取相關(guān)性較大的波長478、674 nm，其相關(guān)系數(shù)分別為-0.64和-0.66。②塊莖形成期，冠層原始光譜分別在波段558～698 nm、726～950 nm與AGB呈極顯著負(fù)、正相關(guān)。相關(guān)性較高的波長為766、778 nm，其相關(guān)系數(shù)分別為0.63和0.63。③塊莖增長期，冠層原始光譜分別在波段454～702 nm、718～950 nm與AGB呈極顯著負(fù)、正相關(guān)。相關(guān)性較大的波長為610、770 nm，其相關(guān)系數(shù)分別為-0.72和0.71。④淀粉積累期，冠層原始光譜分別在波段454～710 nm、714～950 nm與AGB呈極顯著負(fù)、正相關(guān)。相關(guān)性較高的波長為482、618、746 nm，其相關(guān)系數(shù)分別為-0.70、-0.67和0.76。⑤成熟期，冠層原始光譜在波段914～950 nm與AGB呈極顯著負(fù)相關(guān)，未見波段與AGB呈極顯著正相關(guān)。在可見光波段范圍內(nèi)選擇相關(guān)系數(shù)絕對值較大的波長550、694 nm，其相關(guān)系數(shù)為-0.33和0.25。

圖3 馬鈴薯不同生育期原始冠層光譜與AGB相關(guān)系數(shù)Fig.3 Correlation between original canopy spectrum of potato and AGB at different growth stages

2.2 分?jǐn)?shù)階微分光譜與生物量相關(guān)性分析

為篩選出每個生育期與馬鈴薯AGB相關(guān)性較好的分?jǐn)?shù)階微分光譜波段，通過相關(guān)性分析，得到不同生育期馬鈴薯冠層分?jǐn)?shù)階微分光譜與AGB的相關(guān)性，0.4階、0.8階、1.2階、1.6階和2.0階的相關(guān)系數(shù)變化如圖4所示。

根據(jù)馬鈴薯不同生育期冠層分?jǐn)?shù)階微分光譜與AGB的相關(guān)性，得到不同分?jǐn)?shù)階下微分光譜與AGB相關(guān)系數(shù)絕對值最大值，繪制不同生育期下不同階微分光譜與AGB相關(guān)系數(shù)絕對值最大值的折線圖，結(jié)果如圖5所示。分析圖5可知，不同生育期馬鈴薯AGB與分?jǐn)?shù)階微分光譜間的相關(guān)系數(shù)絕對值最大值出現(xiàn)的階數(shù)不同?，F(xiàn)蕾期，相關(guān)系數(shù)絕對值最大值在0.8階微分(470 nm)；塊莖形成期，相關(guān)系數(shù)絕對值最大值在1.8階微分(710 nm)，塊莖增長期和淀粉積累期，相關(guān)系數(shù)絕對值最大值都在1.6階微分(718、722、766 nm)；成熟期，相關(guān)系數(shù)絕對值最大值在1.0階微分(622 nm)。因此表明常用的“綠邊(502～554 nm)”和“紅邊(678～758 nm)”雖對作物營養(yǎng)狀況能夠很好地進(jìn)行監(jiān)測，但僅采用二者范圍內(nèi)的波段進(jìn)行估算AGB時，并不能夠深度挖掘光譜中潛在的信息，導(dǎo)致AGB估算時產(chǎn)生飽和現(xiàn)象，然而利用分?jǐn)?shù)階微分光譜能夠凸顯波段中隱藏的信息，提高光譜信息的敏感度，增強(qiáng)“綠邊”和“紅邊”位置對AGB估算的精度。

圖4 馬鈴薯不同生育期的分?jǐn)?shù)階微分光譜與AGB的相關(guān)系數(shù)Fig.4 Correlation between fractional differential spectrum of potato at different periods and AGB

圖5 不同生育期相關(guān)系數(shù)絕對值最大值出現(xiàn)的階數(shù)Fig.5 Order of maximum absolute value of correlation coefficient in different growth periods

2.3 敏感波段確定

為了盡可能避免各個分?jǐn)?shù)階微分光譜之間的共線性，將與AGB相關(guān)系數(shù)絕對值按從大到小排列，挑選出各個生育期與馬鈴薯AGB相關(guān)的前9個分?jǐn)?shù)階微分冠層光譜波段，并將各生育期的相關(guān)系數(shù)絕對值繪制成矩陣圖，結(jié)果如圖6所示。由圖可知，現(xiàn)蕾期，0.2階微分在波長666 nm處，0.4階微分在波長634 nm處，0.6階微分在波長470、578、594 nm處，0.8階微分在波長458、470 nm處，1.2階微分在波長790 nm處，1.6階微分在波長706 nm處，與馬鈴薯AGB達(dá)到0.01顯著水平，相關(guān)系數(shù)絕對值在0.67～0.72之間；塊莖形成期，1.2階微分在波長714 nm和718 nm處、1.4階微分在波長710 nm和718 nm處、1.6階和1.8階微分都在波長710 nm和714 nm處、2.0階微分在波長710 nm處，與馬鈴薯AGB達(dá)到0.01顯著水平，相關(guān)系數(shù)絕對值在0.72～0.76之間；塊莖增長期，0.8階微分在波長766 nm處、1.0階微分在波長742 nm處、1.2階微分在波長730 nm和894 nm處、1.4階微分在波長726 nm處、1.6和1.8階微分都在波長718 nm和722 nm處，與馬鈴薯AGB達(dá)到0.01顯著水平，相關(guān)系數(shù)絕對值均在0.76以上；淀粉積累期，0階微分在波長746 nm處、0.2階微分在波長730 nm和742 nm處、0.4階微分在波長726 nm處、0.6階微分在波長730 nm處、0.8階微分在波長726 nm處、1.0階和1.2階微分都在波長722 nm處、1.6階微分在766 nm波長處，與馬鈴薯AGB達(dá)到0.01顯著水平，相關(guān)系數(shù)絕對值在0.75～0.78之間；成熟期，0.2階微分和0.4階微分都在922 nm波長處，0.6階微分在波長686 nm和690 nm處，0.8階微分在波長622、682、686 nm處，1.0階微分和1.2階微分都在波長622 nm處，與馬鈴薯AGB達(dá)到0.01顯著水平，相關(guān)系數(shù)絕對值在0.50～0.59之間。綜合每個生育期挑選的前9個微分光譜波長可知，僅采用整數(shù)階微分并不能充分表達(dá)與AGB之間的聯(lián)系，使得冠層光譜中隱藏的信息被遺漏，而通過分?jǐn)?shù)階微分能夠細(xì)化光譜信息，深度挖掘出與AGB相關(guān)的有效信息，使得分?jǐn)?shù)階微分光譜較整數(shù)階微分光譜與AGB的相關(guān)性更高，提高AGB估算模型精度。

圖6 不同生育期馬鈴薯AGB與分?jǐn)?shù)階微分光譜相關(guān)系數(shù)矩陣Fig.6 Correlation coefficient matrix of potato AGB and fractional differential spectrum in different periods

2.4 基于分?jǐn)?shù)階微分光譜的馬鈴薯生物量最優(yōu)模型篩選

以馬鈴薯地上生物量(AGB)為因變量，由0～2階分?jǐn)?shù)階微分(0.2間隔)確定的前9個微分波段為自變量，建立分?jǐn)?shù)階微分冠層光譜與馬鈴薯AGB的多元線性回歸(MLR)模型、隨機(jī)森林(RF)模型和人工神經(jīng)網(wǎng)絡(luò)(ANN)模型，并驗(yàn)證模型精度，進(jìn)而挑選出各個生育期最優(yōu)估算模型，結(jié)果見表1。由表1可知，利用MLR、RF和ANN方法估算馬鈴薯AGB具有較高的精度，從現(xiàn)蕾期到塊莖增長期這3個生育期估算效果逐漸變優(yōu)，驗(yàn)證R2也是不斷增大，RMSE和NRMSE逐漸降低，從淀粉積累期到成熟期建模和驗(yàn)證R2逐漸降低，RMSE和NRMSE逐漸增加，每種方法均在塊莖形成期估算效果最佳。每個生育期的AGB估算，通過MLR構(gòu)建的模型效果最優(yōu)，其次為RF模型，而ANN模型效果最差，且采用3種方法估算AGB得到的模型均是驗(yàn)證效果要優(yōu)于相應(yīng)的建模效果。塊莖增長期，利用3種方法建模R2分別為0.76、0.72和0.61，RMSE分別為143.36、135.06、160.44 kg/hm2，NRMSE分別為16.67%、17.63%、19.58%；驗(yàn)證R2分別為0.84、0.74和0.69，驗(yàn)證RMSE分別為94.72、101.44、159.09 kg/hm2，驗(yàn)證NRMSE分別為11.61%、14.02%和17.76%。

為了弄清不同因素對模型的作用，以不同品種、密度和施肥狀況下的馬鈴薯AGB為因變量，采用上述3種方法構(gòu)建AGB估算模型，并計(jì)算各個模型的精度評價指標(biāo)，其結(jié)果見表2～4。綜合表2～4可知，利用MLR、RF和ANN方法估算3種狀況下的馬鈴薯AGB，也是從現(xiàn)蕾期到塊莖增長期這3個生長期估算效果逐漸變優(yōu)，驗(yàn)證R2也是不斷增大，RMSE和NRMSE逐漸降低，從淀粉積累期到成熟期建模和驗(yàn)證R2逐漸降低，RMSE和NRMSE逐漸增加，每種方法均在塊莖形成期估算效果最佳。每個生育期不同狀況的AGB估算，也是通過MLR構(gòu)建的模型效果最好，其次為RF模型，而ANN模型效果最差。每種狀況下AGB估算，采用3種方法得到的模型均是驗(yàn)證效果要優(yōu)于相應(yīng)的建模效果，其中以不同品種的AGB模型，每種方法以中薯5(P1)構(gòu)建模型的建模和驗(yàn)證精度要優(yōu)于相應(yīng)的以中薯3(P2)所構(gòu)建的模型精度。綜上可知，整體上得到的模型效果與不同狀況下的模型效果基本一致，說明品種、密度和施肥對AGB估算具有等效的影響力。

表1 不同生育期的馬鈴薯AGB估測精度對比Tab.1 Comparison of accuracy of potato AGB estimation at different growth stages

表2 不同品種的馬鈴薯AGB估算精度對比Tab.2 Comparison of accuracy of AGB estimation for different varieties of potatoes

表3 不同密度的馬鈴薯AGB估算精度對比Tab.3 Comparison of accuracy of AGB estimation for different densities of potatoes

表4 不同施肥的馬鈴薯AGB估算精度對比Tab.4 Comparison of accuracy of potato AGB estimation under different fertilizations

綜合建模和驗(yàn)證集可知，每個生育期建模R2低于驗(yàn)證R2，所得R2不僅與精度有關(guān)，也與樣本數(shù)量有關(guān)，整體、不同品種、不同密度和不同施肥條件下建模樣本分別為32個、(14個、18個)、12個和20個，驗(yàn)證樣本分別為16個、(7個、9個)、6個和10個，驗(yàn)證樣本的數(shù)量小于建模數(shù)量，所以驗(yàn)證R2高于建模的R2；驗(yàn)證模型的RMSE和建模的RMSE接近，且二者得到的NRMSE基本處于20%以內(nèi)，說明模型的穩(wěn)定性較好，預(yù)測能力較高。

每個生育期采用3種方法都在塊莖增長期表現(xiàn)效果最好，主要因?yàn)轳R鈴薯從現(xiàn)蕾期開始生長，最初為營養(yǎng)生長和生殖生長，體現(xiàn)在莖節(jié)迅速伸長，葉片快速擴(kuò)大。到了塊莖增長期，莖節(jié)和葉片等營養(yǎng)器官生長最為旺盛，此時植被覆蓋度最大，冠層光譜信息的提取不易受到地面土壤的影響，使得冠層光譜信息能夠充分表達(dá)和AGB之間的聯(lián)系，進(jìn)而構(gòu)建的模型效果最優(yōu)。而淀粉積累期后，地上部有機(jī)物不斷向地下輸送，外加連續(xù)多天大雨，造成地上部葉片迅速枯黃脫落，植被覆蓋度明顯降低，冠層光譜信息的提取容易受到地面土壤的影響，使得冠層光譜信息與AGB的聯(lián)系變差，進(jìn)而構(gòu)建的模型效果不好。

每個生育期利用MLR方法構(gòu)建馬鈴薯AGB估算模型的建模精度和驗(yàn)證模型的穩(wěn)定性優(yōu)于RF和ANN模型，主要因?yàn)闄C(jī)器學(xué)習(xí)適用于處理較大數(shù)據(jù)集，而對較小數(shù)據(jù)集表現(xiàn)優(yōu)勢不明顯，本文建模樣本和驗(yàn)證樣本都屬于較小數(shù)據(jù)集，所以RF和ANN的表現(xiàn)能力較差一些。而ANN模型估算精度和驗(yàn)證模型的穩(wěn)定性最差，可能是利用ANN方法構(gòu)建模型時，在訓(xùn)練過程中不斷地反復(fù)學(xué)習(xí)，過多地學(xué)習(xí)了細(xì)節(jié)，造成估算能力不佳。每個生育期利用3種方法估算馬鈴薯AGB時，每種方法的建模精度由高到低依次為塊莖增長期、塊莖形成期、淀粉積累期、現(xiàn)蕾期、成熟期。

3 結(jié)論

(1)以馬鈴薯為研究對象，利用無人機(jī)冠層高光譜數(shù)據(jù)，采用分?jǐn)?shù)階微分方法以0～2階(間隔0.2)進(jìn)行光譜微分處理，篩選出前9個分?jǐn)?shù)階微分光譜，結(jié)合地上實(shí)測的48個生物量數(shù)據(jù)，建立了估算馬鈴薯生物量的MLR、RF和ANN模型。

(2)在不同生育期，與AGB相關(guān)性較高的原始冠層光譜波段不同?，F(xiàn)蕾期，在可見光范圍內(nèi)，與AGB極顯著相關(guān)的波長為478、674 nm，其相關(guān)系數(shù)分別為-0.64、-0.66；塊莖形成期，與AGB極顯著相關(guān)的波長為766、778 nm，其相關(guān)系數(shù)均為0.63；塊莖增長期，與AGB極顯著相關(guān)的波長為610、770 nm，其相關(guān)系數(shù)分別為-0.72、0.71；淀粉積累期，與AGB極顯著相關(guān)的波長為482、618、746 nm，其相關(guān)系數(shù)分別為-0.70、-0.67和0.76；成熟期，冠層原始光譜在可見光范圍內(nèi)未存在極顯著正負(fù)相關(guān)的波長，而相關(guān)性較好的波長為550、694 nm，其相關(guān)系數(shù)為-0.33、0.25。

(3)不同生育期分?jǐn)?shù)階微分光譜與AGB相關(guān)系數(shù)絕對值最大值出現(xiàn)的階數(shù)不同?，F(xiàn)蕾期，相關(guān)系數(shù)絕對值最大值在0.8階微分(470 nm)，為0.72；塊莖形成期，相關(guān)系數(shù)絕對值最大值在1.8階微分(710 nm)，為0.76；塊莖增長期和淀粉積累期，相關(guān)系數(shù)絕對值最大值都在1.6階微分(718、722、766 nm)，分別為0.77和0.78；成熟期，相關(guān)系數(shù)絕對值最大值在1.0階微分(622 nm)，為0.59。通過比較馬鈴薯成熟期的原始光譜和微分光譜與AGB的相關(guān)性可知，分?jǐn)?shù)階微分光譜更能挖掘出光譜數(shù)據(jù)中隱藏的信息。

(4)由每個生育期挑選的前9個微分光譜波段可知，相較于整數(shù)階微分，分?jǐn)?shù)階微分與AGB的相關(guān)性更高。基于分?jǐn)?shù)階微分光譜的AGB估算模型，不同生育期利用MLR方法，以整體、不同品種、不同密度和不同施肥條件下的AGB為自變量，所得模型效果最好，其次是RF模型，而ANN模型效果最差。不同生育期利用3種方法估算馬鈴薯AGB時，其建模精度由高到低依次為塊莖增長期、塊莖形成期、淀粉積累期、現(xiàn)蕾期、成熟期。