赫曉慧， 馮坤， 郭恒亮， 田智慧

(1.鄭州大學(xué)地球科學(xué)與技術(shù)學(xué)院(智慧城市研究院)，河南 鄭州 450001；2.鄭州大學(xué)水利科學(xué)與工程學(xué)院，河南 鄭州 450001)

高產(chǎn)是大豆育種的主要目標(biāo)，高產(chǎn)育種的基礎(chǔ)是群體的光合效率[1]，但是如何快速、有效地篩選高產(chǎn)大豆種子卻是個(gè)難題。葉面積指數(shù)是常用的表征群體特征的參數(shù)[2]，在一定程度上，作物產(chǎn)量隨LAI的增大而提高，可以很好地反映農(nóng)作物生長(zhǎng)的狀況與趨勢(shì)。因此，有效地估測(cè)LAI對(duì)大豆育種具有重要意義。

目前，對(duì)農(nóng)作物生長(zhǎng)參數(shù)的反演方法主要有經(jīng)驗(yàn)?zāi)Ｐ头?、物理模型法和綜合反演法[3]。經(jīng)驗(yàn)?zāi)Ｐ头ㄍǔＪ窃诠庾V和相關(guān)生化參量之間建立回歸經(jīng)驗(yàn)?zāi)Ｐ?，該方法雖然簡(jiǎn)單、快速、易于操作，但該方法受限于飽和效應(yīng)、場(chǎng)景和遙感傳感器等方面的影響，反演結(jié)果的準(zhǔn)確性和普適性難以得到保證。物理模型法則是綜合考慮了葉片、冠層、土壤、觀測(cè)情況等因素，更能反映復(fù)雜的地表環(huán)境，反演精度相對(duì)較高，具有較強(qiáng)的普適性和外延性。目前，物理模型中應(yīng)用最為廣泛的物理模型則是PROSAIL模型[4]，如LI等[5-6]結(jié)合農(nóng)學(xué)知識(shí)，將PROSAIL模型應(yīng)用于作物監(jiān)測(cè)中，證明其在估測(cè)葉綠素含量等參數(shù)上的優(yōu)勢(shì)。楊貴軍等[7]利用PROSAIL模型反演春小麥的LAI時(shí)指出：實(shí)測(cè)值與觀測(cè)值誤差較小，可以準(zhǔn)確地反映出春小麥長(zhǎng)勢(shì)。張明政等[8]通過(guò)利用PROSAIL模型對(duì)玉米冠層葉面積指數(shù)和葉綠素含量反演，表明在三葉期時(shí)的LAI反演結(jié)果和拔節(jié)期的LCC反演結(jié)果總體精度最高，分別為0.62和0.59，為農(nóng)作物反演策略提功力研究思路和理論依據(jù)。綜合反演法包括神經(jīng)網(wǎng)絡(luò)、有影追蹤和回歸樹(shù)方法等，該類方法具有反演速度快、無(wú)需參數(shù)初始值和能有效考慮關(guān)鍵變量的變化等優(yōu)點(diǎn)。在綜合反演法中，BP神經(jīng)網(wǎng)絡(luò)模型因在處理復(fù)雜的非線性問(wèn)題，及模式識(shí)別方面有著無(wú)可比擬的優(yōu)勢(shì)而被廣泛采用[9]。陳士城等[10]在LAI的遙感反演中，從網(wǎng)絡(luò)隱含層層數(shù)、隱含層節(jié)點(diǎn)數(shù)和學(xué)習(xí)速率等3個(gè)方面進(jìn)行分析與討論，確定了最佳的BP網(wǎng)絡(luò)結(jié)構(gòu)，取得了較好的反演結(jié)果。劉俊等[11]通過(guò)構(gòu)建不同隱含層的BP神經(jīng)網(wǎng)絡(luò)模型，表明基于BP神經(jīng)網(wǎng)絡(luò)模型的夏玉米LAI反演精度整體較好，符合實(shí)際情況，基本吻合夏玉米實(shí)際生長(zhǎng)情況。BP神經(jīng)網(wǎng)絡(luò)模型是一種按照誤差反向傳播算法算法訓(xùn)練的多層前饋網(wǎng)絡(luò)，通常情況下容易陷入局部最優(yōu)，這時(shí)需要利用遺傳算法[12-14]優(yōu)化其權(quán)值和閾值，利用其全局搜索能力，逐步搜尋最優(yōu)解。

BP神經(jīng)網(wǎng)絡(luò)模型和PROSAIL模型在農(nóng)作物生長(zhǎng)參數(shù)反演中占據(jù)著極其重要的地位，但不足的是以上研究?jī)H針對(duì)單一農(nóng)作物類型而言，模型在同一農(nóng)作物的多個(gè)基因型種群中的優(yōu)劣并未提及。因此，本研究以人工定向培育的RIL大豆種群和黃河流域部分地區(qū)自然野生大豆種群為數(shù)據(jù)源，利用PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型，開(kāi)展大豆種群冠層葉片的LAI建模反演，并結(jié)合野外實(shí)測(cè)數(shù)據(jù)對(duì)反演結(jié)果進(jìn)行驗(yàn)證與評(píng)價(jià)，以期在不同的大豆種群生育期內(nèi)中找出最好的反演模型，為大豆長(zhǎng)勢(shì)監(jiān)測(cè)和規(guī)?；N提供更為簡(jiǎn)單有效的一種方法。

1 數(shù)據(jù)與方法

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

研究區(qū)位于河南省新鄉(xiāng)市河南農(nóng)業(yè)大學(xué)原陽(yáng)試驗(yàn)基地(113°56′27.95″E，35°6′27.18″N)，試驗(yàn)區(qū)總面積約為4 800 m2，地勢(shì)平坦，土壤肥沃。時(shí)間為2019-06-20至2019-09-26，供試品種為分別為重組自交系大豆和自然野生大豆，其中RIL種群大豆由親本為DN50和W82的子一代自交9次所得，自然野生種群大豆的種子來(lái)源則是黃河流域附近的野生大豆?？紤]到自然野生種群大豆品種多，基因型更為豐富，且2個(gè)種群需做對(duì)照分析的因素，本次試驗(yàn)分為2塊區(qū)域種植，其中RIL種群種植面積約1 200 m2，劃分為32個(gè)樣區(qū)，自然野生種群約3 600 m2，劃分為56個(gè)樣區(qū)。由于2個(gè)大豆種群的生育期時(shí)間存在差異，RIL種群比自然野生種群生育期晚8 d左右，因此分別對(duì)2個(gè)種群進(jìn)行了7次觀測(cè)，每次觀測(cè)的具體日期如表1所示。

1.2 光譜和LAI的測(cè)定與預(yù)處理

光譜測(cè)定采用的儀器為美國(guó)ASD公司(Analytical Spectral Devices，ASD)生產(chǎn)的FiledSpecHH2手持式地物光譜儀，波長(zhǎng)范圍為325～1 075 nm，精度±1 nm，采樣間隔1.5 nm。為了有效地減少外部環(huán)境對(duì)數(shù)據(jù)噪聲影響，盡可能地選擇晴朗無(wú)云或少云，無(wú)風(fēng)的天氣，每次光譜的測(cè)定時(shí)間都為9:30—12:00。光譜采集時(shí)，將儀器探頭置于樣本上方約15 cm處，采集10組光譜數(shù)據(jù)，取這10組光譜平均值作為該樣區(qū)的最終大豆冠層光譜反射率；經(jīng)統(tǒng)計(jì)，RIL種群共224組數(shù)據(jù)，自然野生種群共392組數(shù)據(jù)。

大豆葉面積指數(shù)LAI的采集使用儀器為英國(guó)Delta-T公司生產(chǎn)的SunScan無(wú)線型(SS1-COM-R4)植被冠層分析儀。LAI采集應(yīng)與光譜采集同步進(jìn)行，將探桿置于大豆冠層下方，并水平測(cè)量5次，取其平均值作為該樣區(qū)的最終大豆葉面積指數(shù)。

盡管同種群同生育期內(nèi)大豆冠層光譜為同一短時(shí)間段內(nèi)所測(cè)，但是各個(gè)冠層葉片的理化參數(shù)有所不同，數(shù)據(jù)采集時(shí)儀器與樣本的距離和角度不可避免的存在差異。為充分發(fā)掘出各參數(shù)的不同組合所模擬模的光譜信息，因此在全生育期模擬光譜中隨機(jī)選擇40組數(shù)據(jù)用于驗(yàn)證集，典型生育期模擬光譜則隨機(jī)選擇20組數(shù)據(jù)用于驗(yàn)證集，其余數(shù)據(jù)則分別作為各數(shù)據(jù)集的建模集。

1.3 遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)模型及其反演策略

BP神經(jīng)網(wǎng)絡(luò)模型是一種按照誤差逆向傳播算法訓(xùn)練的多層前饋網(wǎng)絡(luò)，是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)模型之一。本研究在BP神經(jīng)網(wǎng)絡(luò)模型的基礎(chǔ)上，引入遺傳算法對(duì)其進(jìn)行優(yōu)化，以期對(duì)LAI的估測(cè)達(dá)到更好的精度要求。其基本思路是將BP神經(jīng)網(wǎng)絡(luò)算法和遺傳算法相結(jié)合，當(dāng)BP算法訓(xùn)練網(wǎng)絡(luò)出現(xiàn)收斂速度緩慢時(shí)，將BP神經(jīng)網(wǎng)絡(luò)模型的各隱含層節(jié)點(diǎn)的閾值和權(quán)值作為遺傳算法的輸入信息，并將它們編碼成染色體，接著運(yùn)用遺傳算法的選擇、交叉、變異產(chǎn)生新的后代，最后將結(jié)果作為BP算法的初始值，繼續(xù)BP算法訓(xùn)練網(wǎng)絡(luò)，該操作重復(fù)運(yùn)行，直至達(dá)到問(wèn)題要求的誤差精度，輸出值為全局最優(yōu)。

本次研究中需要進(jìn)行遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)模型是一種單向傳播的多層前向網(wǎng)絡(luò)，包括3個(gè)層級(jí)，分別是輸入層、隱含層和輸出層。研究擬通過(guò)不同大豆種群實(shí)測(cè)冠層葉片反射率來(lái)反演LAI，故模型輸入層為實(shí)測(cè)反射率，輸出層為L(zhǎng)AI。如表2所示，考慮到可見(jiàn)光與近紅外區(qū)的光譜反射率差異較大，而以神經(jīng)網(wǎng)絡(luò)來(lái)擬合實(shí)測(cè)數(shù)據(jù)需要發(fā)揮不同波段間的互補(bǔ)優(yōu)勢(shì)，為了確保不遺漏有用的信息，因此在藍(lán)色、綠色、紅色和近紅外4個(gè)波段范圍內(nèi)，分別選擇與LAI相關(guān)性指數(shù)最大、在P=0.05水平上顯著相關(guān)的波長(zhǎng)，并將其作為各波段范圍內(nèi)的敏感波段。

表2 不同大豆種群敏感波段的選擇

PROSAIL模型是通過(guò)耦合PRSPECT葉片輻射傳輸模型和SAILH冠層結(jié)構(gòu)模型得到的整體模型。PROSPECT模型主要描述植物葉片在400～2 500 nm光譜范圍內(nèi)的光學(xué)特征，本研究采用常用波段為325～1 075 nm，并將波段進(jìn)行5點(diǎn)平滑處理及歸一化處理，通過(guò)輸入葉片結(jié)構(gòu)參數(shù)N、葉綠素含量Cab、干物質(zhì)含量Cm和等效水厚度Cw等，可以得到325～1 075 nm光譜范圍內(nèi)的葉片透射率和反射率。SAIL模型是在水平均勻的假設(shè)下主要考察冠層的垂直分層結(jié)構(gòu)和葉傾角分布，輸入?yún)?shù)主要包括葉片光譜信息，LAI、平均葉傾角、土壤參數(shù)、熱點(diǎn)參數(shù)、漫反射系數(shù)、太陽(yáng)天頂角、觀測(cè)相對(duì)方位角等，其中葉片光譜信息(反射率、透射率)采用PROSPECT模型模擬的輸出結(jié)果，其他參數(shù)則根據(jù)光譜獲取時(shí)的實(shí)際觀測(cè)信息確定。

本研究首先根據(jù)不同大豆種群冠層光譜采集過(guò)程的實(shí)際情況，得到查找表的各輸入?yún)?shù)及變化范圍；然后通過(guò)將各參數(shù)的不同組合帶入耦合之后的PROSAIL模型進(jìn)行冠層反射率模擬，建立不同大豆種群葉面積指數(shù)的查找表；接著在敏感波段范圍內(nèi)，將實(shí)測(cè)光譜發(fā)射率與查找表中反射率比較，利用代價(jià)函數(shù)確定最優(yōu)解；最后通過(guò)最優(yōu)解，便可利用實(shí)測(cè)光譜反射率對(duì)葉面積指數(shù)做出估測(cè)。

1.5 模型驗(yàn)證指標(biāo)

利用LAI野外實(shí)測(cè)數(shù)據(jù)從穩(wěn)定性、估測(cè)能力和估測(cè)精度3個(gè)方面對(duì)反演結(jié)果進(jìn)行驗(yàn)證，穩(wěn)定性使用確定系數(shù)R2來(lái)衡量，估測(cè)能力使用均方根誤差RMSE(root mean square error，RMSE)來(lái)衡量，估測(cè)精度使用EA來(lái)衡量，計(jì)算公式如下：

(1)

(2)

(3)

2 結(jié)果與分析

PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型所模擬的光譜反射率是成功反演植被生長(zhǎng)參數(shù)的前提條件，因此本研究首先分別從2個(gè)大豆種群各生育期內(nèi)的實(shí)測(cè)光譜中選擇與LAI相關(guān)性較好的數(shù)據(jù)進(jìn)行平均光譜計(jì)算，然后對(duì)2個(gè)模型模擬光譜進(jìn)行評(píng)價(jià)與模型結(jié)構(gòu)調(diào)整；接著將全生育期和典型生育期的驗(yàn)證集中實(shí)測(cè)光譜與PROSAIL模型所構(gòu)建的查找表中的模擬光譜作比較，并利用代價(jià)函數(shù)對(duì)LAI進(jìn)行合理性的估測(cè)，而遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型則是以實(shí)測(cè)光譜的敏感波段作為模型輸入信息，以LAI作為模型輸出信息，最后以各生育期內(nèi)實(shí)測(cè)點(diǎn)LAI分別作為2個(gè)模型精度驗(yàn)證的依據(jù)，對(duì)比分析不同模型在不同種群中的反演精度。

2.1 實(shí)測(cè)光譜與模擬光譜對(duì)比分析

圖1展示了研究區(qū)內(nèi)的全生育期綜合模擬光譜與實(shí)測(cè)光譜的分析結(jié)果。分析發(fā)現(xiàn)，PROSAIL模型相比于遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型，在對(duì)RIL種群大豆冠層葉片光譜的模擬中有著更好的效果，其模擬光譜與實(shí)測(cè)光譜擬合效果最佳，在整個(gè)可見(jiàn)光和近紅外波段內(nèi)與實(shí)測(cè)光譜的變化趨勢(shì)相同；而遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型對(duì)RIL種群的模擬光譜，在可見(jiàn)光波段內(nèi)與實(shí)測(cè)光譜擬合較好，同PROSAIL模型差別不大，在近紅外波段內(nèi)其擬合程度較PROSAIL模型差。自然野生種群中，2個(gè)模型的模擬光譜均高于實(shí)測(cè)光譜，PROSAIL模型的模擬光譜在可見(jiàn)光波段內(nèi)較遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型模擬光譜低，更加貼合實(shí)測(cè)光譜，而遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型在近紅外波段內(nèi)的模擬光譜卻優(yōu)于PROSAIL模型模擬光譜，擬合效果更好。

圖1 全生育期綜合模擬光譜與實(shí)測(cè)光譜

為進(jìn)一步發(fā)掘PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型在不同種群典型生育期內(nèi)的差異性，分別利用所構(gòu)建的模型對(duì)RIL和自然野生大豆種群始花期(R1)、始莢期(R3)和始粒期(R5)進(jìn)行光譜的模擬，以各生育期內(nèi)所有實(shí)測(cè)光譜平均值作為評(píng)價(jià)依據(jù)分析，結(jié)果如圖2所示。RIL種群中，遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型對(duì)R1期的模擬光譜可見(jiàn)光波段高于實(shí)測(cè)光譜，PROSAIL模型對(duì)R3期和R5期的模擬光譜在可見(jiàn)光波段內(nèi)低于實(shí)測(cè)光譜；近紅外波段內(nèi)，2個(gè)模型對(duì)R1期的模擬光譜均高于實(shí)測(cè)光譜，R3期內(nèi)的模擬光譜均低于實(shí)測(cè)光譜，遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型模擬效果較PROSAIL模型差；R5期時(shí)遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型模擬值低于實(shí)測(cè)光譜，PROSAIL模型模擬值高于實(shí)測(cè)光譜，后者誤差較小。自然野生種群中，可見(jiàn)光波段內(nèi)典型生育期的模擬光譜或高或低的分布在實(shí)測(cè)光譜上下，與其差異較??；近紅外波段內(nèi)PROSAIL模型對(duì)R3期的模擬光譜最為突出，與實(shí)測(cè)光譜誤差最大，而遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型對(duì)R1期和R3期的模擬光譜效果都較好，R5期擬合度最佳，模擬光譜與實(shí)測(cè)光譜變化趨勢(shì)一樣。

圖2 典型生育期模擬光譜與實(shí)測(cè)光譜

在RIL種群中，PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型的LAI反演結(jié)果與野外實(shí)測(cè)值擬合效果都較好，且回歸模型t檢驗(yàn)P值均小于0.05，在α=0.05的水平上顯著。其中，PROSAIL模型確定系數(shù)R2最大，RMSE最小(表2)，說(shuō)明對(duì)應(yīng)試驗(yàn)點(diǎn)的反演結(jié)果與實(shí)測(cè)值離散程度較小，較好的分布在一階擬合線上下；遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型雖在RIL種群中反演能力弱于PROSAIL模型，但是對(duì)比該模型在自然野生種群中的結(jié)果參數(shù)發(fā)現(xiàn)，其依然存在著一定優(yōu)勢(shì)，更適用于品種繁多的大豆種群。在自然野生種群中，PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型的反演結(jié)果同樣較為均勻地分布在一階擬合線兩側(cè)，但是對(duì)比2個(gè)模型在RIL種群中的擬合效果發(fā)現(xiàn)，PROSAIL模型模擬值同實(shí)測(cè)值差距較大，該模型在品種繁多的大豆種群中的性能體現(xiàn)較遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型差；對(duì)比RMSE發(fā)現(xiàn)，遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型為0.12，PROSAIL模型為0.15，說(shuō)明二者數(shù)值雖然相差不多，但是考慮到LAI大面積觀測(cè)所需要的精確性，PROSAIL模型對(duì)自然野生大豆種群LAI的反演存在明顯的弊端，反演效果弱于遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型，模擬值與實(shí)測(cè)值差距較大，產(chǎn)生的誤差不利于對(duì)品種繁多的大豆種群的LAI進(jìn)行估測(cè)。

如表3所示，RIL種群所構(gòu)建的模型中，PROSAIL模型性能最好，R2為0.85，RMSE為0.11，EA為96%，自然野生種群的建模集中，遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型性能最好，R2為0.83，RMSE為0.12，EA為94%，說(shuō)明PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型不同種群中反演LAI的能力各有優(yōu)勢(shì)和不足，人工定向培育的大豆具有一定的目的性，在生理性狀和表型特征上更為接近，通過(guò)PROSAIL模型使用其葉片理化參數(shù)更容易獲得LAI的反演值，且精度較高；而對(duì)種類繁多的自然野生種群而言，遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型在選擇、交叉和變異后更容易找到與實(shí)測(cè)值接近的最優(yōu)解，具有很好的LAI反演能力。

2.3 典型生育期反演結(jié)果分析

針對(duì)不同的大豆種群進(jìn)行LAI的反演工作，需根據(jù)實(shí)際情況來(lái)選擇合適的模型，以2個(gè)大豆種群R1期、R3期和R5期實(shí)測(cè)LAI為例，結(jié)果如表4所示。在R1期階段，PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型穩(wěn)定性較差，這是因?yàn)镽1期大豆冠層葉片較小，參與建模的光譜信息容易受到背景信息的干擾而降低模型精度。R3期是大豆生育最旺盛的時(shí)期，隨著冠層葉片逐漸增大，地面噪聲信息影響程度減弱，模型穩(wěn)定性、估測(cè)能力和估測(cè)精度都較R1期有了提升。R5期時(shí)，在RIL種群中，PROSAIL模型確定性系數(shù)為0.89，高于遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型，而在自然野生種群中，遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型性能優(yōu)于PROSAIL模型，典型生育期內(nèi)LAI反演結(jié)果與全生育期內(nèi)LAI 反演結(jié)果大致相同。

表4 典型生育期LAI反演模型估測(cè)精度

為更加直觀地觀察模型對(duì)未知樣本的適應(yīng)性，將PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型在不同種群中的估測(cè)數(shù)據(jù)與對(duì)應(yīng)實(shí)測(cè)值進(jìn)行回歸擬合，R5期結(jié)果如圖5、圖6所示。R5期時(shí)，PROSAIL模型在RIL種群中的估測(cè)值和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型在自然野生種群中的預(yù)估測(cè)值，二者都與實(shí)測(cè)值較為接近，估測(cè)結(jié)果沿直線y=x的分布狀態(tài)較優(yōu)，擬合效果較好。R5期是大豆種子形成的重要時(shí)期，對(duì)比分析各生育期內(nèi)模型反演趨勢(shì)線并結(jié)合表4發(fā)現(xiàn)，隨著大豆生育期的變化，模型確定系數(shù)也在不斷增加，均方根誤差不斷減小，估測(cè)精度逐步提高，越接近R5期，模型的可行性越強(qiáng)，遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型在自然野生種群中最大為85%，PROSAIL模型在RIL種群中最大為97%，通過(guò)利用2個(gè)不同的LAI反演模型，分別對(duì)RIL和自然野生大豆種群R5期LAI進(jìn)行模擬估測(cè)，可以在大豆種子形成之前對(duì)其長(zhǎng)勢(shì)做出精確估測(cè)，并施加一定外界干擾以達(dá)到規(guī)?；N的目的。

圖3 RIL種群LAI估測(cè)精度

圖4 自然野生種群LAI估測(cè)精度

圖5 RIL種群R5期LAI估測(cè)精度

圖6 自然野生種群R5期LAI估測(cè)精度

3 結(jié)論

本研究利用PROSAIL模型和遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)模型分別對(duì)RIL和自然野生大豆種群冠層葉片LAI進(jìn)行反演，利用實(shí)測(cè)數(shù)據(jù)對(duì)全生育期和典型生育期內(nèi)反演結(jié)果進(jìn)行精度驗(yàn)證。結(jié)果表明，全生育期內(nèi)，PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型對(duì)2個(gè)大豆種群冠層葉片LAI反演結(jié)果總體精度較高，與實(shí)測(cè)值的相關(guān)性較好，PROSAIL模型在RIL種群中的R2最大為0.85，RMSE為0.11，EA為96%，而遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型在自然野生大豆種群中的R2最大為0.83，RMSE為0.12，EA為94%。典型生育期內(nèi)，在對(duì)RIL種群大豆冠層葉片LAI反演中，PROSAIL模型表現(xiàn)出了更優(yōu)異的反演能力，而對(duì)自然野生大豆種群LAI反演中，遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型表現(xiàn)出了更好的適用性，并且上述模型在R5期時(shí)性能均為最佳，R2分別為0.89和0.85，RMSE分別為0.11和0.13，EA均為97%，2個(gè)模型在不同大豆種群R5期內(nèi)的反演性能均優(yōu)于全生育期綜合反演性能，并且該精度能夠滿足農(nóng)業(yè)區(qū)域植被理化參數(shù)信息反演的需求；通過(guò)對(duì)R5時(shí)期LAI的估測(cè)，可以較為精確的反映出大豆長(zhǎng)勢(shì)情況，在其種子形成前作出人工干預(yù)，以此來(lái)方便大豆的規(guī)模化育種工作。