高洪燕，毛罕平，倪紀(jì)恒，公麗艷，張曉東

(1. 江蘇大學(xué)農(nóng)業(yè)工程學(xué)院，江蘇 鎮(zhèn)江212013； 2. 沈陽(yáng)師范大學(xué)糧食學(xué)院，遼寧 沈陽(yáng)110034)

生菜是重要的綠葉蔬菜之一，因營(yíng)養(yǎng)豐富、生熟食俱佳等優(yōu)點(diǎn)，栽培面積逐年擴(kuò)大.據(jù)統(tǒng)計(jì)，中國(guó)生菜產(chǎn)量占世界總產(chǎn)量的56%左右，而北京市的生菜年需求量達(dá)20萬(wàn)t.生菜屬于潛根系作物，生長(zhǎng)發(fā)育期短，一般40 d左右就可采收食用，因而對(duì)肥料的需求很敏感[1].實(shí)際生產(chǎn)中越來(lái)越多的農(nóng)民通過(guò)增加化肥的投入來(lái)提高生菜產(chǎn)量和經(jīng)濟(jì)效益，這不僅會(huì)造成土壤板結(jié)，有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)降低，酶活性減弱，也會(huì)引起生菜品質(zhì)下降，甚至減產(chǎn)[2].因此，研究不同施肥量下設(shè)施生菜養(yǎng)分質(zhì)量分?jǐn)?shù)、長(zhǎng)勢(shì)及產(chǎn)量三者的關(guān)系，是提高肥料利用率、有效發(fā)揮化肥在生菜生產(chǎn)中的作用的關(guān)鍵因素.

不同施肥制度對(duì)生菜長(zhǎng)勢(shì)及產(chǎn)量的影響有強(qiáng)弱和主次之分.通常采用多元回歸分析、多級(jí)通徑分析、主成分分析及相關(guān)分析等方法確定影響強(qiáng)度.MAH-LANGU等[3]研究了水培條件下不同施氮量對(duì)生菜長(zhǎng)勢(shì)及品質(zhì)的影響.GALIENI 等[4]分析了氮磷元素缺乏對(duì)生菜葉片形態(tài)和生理特性的影響.這些方法均基于“白色系統(tǒng)”，即需要大量數(shù)據(jù)，且要求數(shù)據(jù)具有典型的概率分布，但農(nóng)業(yè)生產(chǎn)中很難滿足這2個(gè)條件[5].農(nóng)業(yè)系統(tǒng)本身是一種灰色系統(tǒng)[6]，不同施肥制度與長(zhǎng)勢(shì)和產(chǎn)量間具有一定的灰色關(guān)聯(lián)性[7].在養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)關(guān)系研究方面，大多是針對(duì)整個(gè)生長(zhǎng)期建模，未考慮作物生長(zhǎng)規(guī)律與養(yǎng)分吸收規(guī)律.

因此，文中利用灰色關(guān)聯(lián)分析方法研究生菜養(yǎng)分質(zhì)量分?jǐn)?shù)、長(zhǎng)勢(shì)及產(chǎn)量三者間的關(guān)系，擬解決因數(shù)據(jù)量少及分布不典型的難題，明確影響生菜長(zhǎng)勢(shì)及產(chǎn)量的主次因素.在此基礎(chǔ)上，進(jìn)一步研究養(yǎng)分質(zhì)量分?jǐn)?shù)與產(chǎn)量關(guān)系，并根據(jù)生菜長(zhǎng)勢(shì)，建立不同時(shí)期養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)的關(guān)系模型，為設(shè)施生菜科學(xué)施肥提供科學(xué)依據(jù)和理論指導(dǎo).

1 材料與方法

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

試驗(yàn)在江蘇大學(xué)Venlo型溫室內(nèi)進(jìn)行，試驗(yàn)品種為意大利全年耐抽苔半結(jié)球生菜(南京沃蔬種業(yè)有限公司).采用營(yíng)養(yǎng)液加珍珠巖的無(wú)土栽培模式，將長(zhǎng)勢(shì)狀況相似的五葉一心幼苗定植到花盆中.研究選用山崎營(yíng)養(yǎng)液配方，共分為13個(gè)施肥處理，分別為標(biāo)準(zhǔn)配方中氮(N)、磷(P)和鉀(K)質(zhì)量分?jǐn)?shù)的25%，50%，75%，100%和125%.標(biāo)準(zhǔn)配方中含Ca(NO3)2·4H2O，236.00 mg/L；KNO3，404.00 mg/L；NH4H2PO4，57.00 mg/L；MgSO4·7H2O，123.00 mg/L；Fe-EDTA，16.00 mg/L；MnCl2·4H2O，1.20 mg/L；H3BO3，0.72 mg/L；ZnSO4·4H2O， 0.09 mg/L；CuSO4·5H2O，0.04 mg/L；(NO4)2MoO4，0.01 mg/L.其中NO3-濃度為6.0 mmol/L，NH4+濃度為0.5 mmol/L，K+濃度為4.0 mmol/L，PO43-濃度為0.5 mmol/L.在保證其他營(yíng)養(yǎng)元素濃度不變的情況下，缺失和過(guò)量處理按照標(biāo)準(zhǔn)配方中氮磷鉀質(zhì)量分?jǐn)?shù)進(jìn)行等比例減少或增加，不同水平氮磷鉀元素營(yíng)養(yǎng)液配方離子濃度c如表1所示.試驗(yàn)在春季進(jìn)行，共栽種生菜3批，每隔7 d栽種1批，每個(gè)處理6個(gè)樣本.根據(jù)天氣變化，通過(guò)控制風(fēng)機(jī)、濕簾、遮陽(yáng)網(wǎng)等設(shè)備調(diào)控溫室環(huán)境信息，使夜溫不低于15 ℃，晝溫不高于30 ℃，光強(qiáng)為200～400 μmol/(m2·s).

表1 不同水平氮磷鉀元素營(yíng)養(yǎng)液配方離子濃度

1.2 設(shè)施生菜長(zhǎng)勢(shì)信息獲取

自緩苗成功起，每隔3 d上午8：00開(kāi)始獲取生菜冠層俯視及主視圖像.在圖像采集前，對(duì)相機(jī)進(jìn)行白平衡標(biāo)定，并使用光圈優(yōu)先模式，光圈設(shè)定在F8，ISO為100，保證采集的圖像色彩還原準(zhǔn)確且有足夠的景深使生菜成像清晰，圖像分辨率為3 888×2 592.采用延遲拍攝模式，消除手動(dòng)操作帶來(lái)的微小抖動(dòng).采用標(biāo)準(zhǔn)坐標(biāo)紙為后期圖像處理的參考標(biāo)尺.采集主視圖像時(shí)，將花盆順時(shí)針旋轉(zhuǎn)90°，再次采集生菜的主視圖像，用以消除因個(gè)別葉片長(zhǎng)勢(shì)突出引起的偶然誤差.

1.3 設(shè)施生菜產(chǎn)量及氮磷鉀質(zhì)量分?jǐn)?shù)測(cè)定

采集的生菜去根后，立即用精度為0.1 mg的天平稱取質(zhì)量，單位為g.然后放入烘箱在105 °C下殺青1 h，再以75 °C烘干至恒重.利用H2SO4-H2O2消解0.1～0.2 g(過(guò)0.25～0.50 mm篩)磨細(xì)烘干后的生菜葉片樣本，氮和磷質(zhì)量分?jǐn)?shù)測(cè)定采用Auto Analyzer 3型連續(xù)流動(dòng)分析儀，鉀質(zhì)量分?jǐn)?shù)的測(cè)定采用火焰光度計(jì).營(yíng)養(yǎng)質(zhì)量分?jǐn)?shù)以單位干質(zhì)量的百分比(%)形式表示.

1.4 設(shè)施生菜長(zhǎng)勢(shì)特征提取

在不同施肥制度下，生菜長(zhǎng)勢(shì)會(huì)發(fā)生明顯變化.缺氮時(shí)，生菜生長(zhǎng)緩慢且外觀呈仿錐形；缺磷時(shí)，生菜葉片伸展度下降且葉片數(shù)減少；缺鉀時(shí)，生菜組織失去膨壓并萎蔫.STORY等[8]通過(guò)提取俯視圖像中冠幅投影面積ATPC(top projected canopy area)和冠幅周長(zhǎng)PTPC(top projected canopy perimeter)反映生菜缺肥情況.文中在此基礎(chǔ)上，結(jié)合營(yíng)養(yǎng)脅迫形態(tài)特征的變化，增加提取主視圖像中的株高HP(plant height)信息描述生菜的長(zhǎng)勢(shì)情況，并將2次主視圖像獲取HP的均值作為最終結(jié)果.

在Matlab軟件中運(yùn)用“2G-R-B”算子的分割冠幅區(qū)域，該方法計(jì)算量小且在較寬范圍光照條件下，能夠很好地提取生菜冠幅[9].圖1為生菜冠層俯視圖像和主視圖像二值化圖像，圖2為背景分割圖像.采用Roberts，Prewitt，Sobel，LoG，Zerocross和Canny算子對(duì)生菜進(jìn)行邊緣提取，通過(guò)分析不同算子提取結(jié)果可知，Roberts算子因定位不準(zhǔn)確，提取的邊緣較粗；Prewitt和Sobel算子提取的結(jié)果相近，在降低噪聲和偽邊緣的同時(shí)也平滑了真正邊緣；LoG算子和Zerocross算子提取結(jié)果相近，定位精度高，對(duì)噪聲較敏感，在一定的程度上模糊了目標(biāo)圖像邊緣[10]；而Canny算子對(duì)圖像進(jìn)行了“非極大值抑制”和形態(tài)學(xué)連接操作，所獲得邊緣特征輪廓清晰，且封閉性好.因此，文中采用Canny算子提取生菜周長(zhǎng)，提取結(jié)果如圖3所示，圖中i,j分別為長(zhǎng)度和高度方向像素個(gè)數(shù).最終得到ATPC,PTPC和HP隨時(shí)間變化曲線如圖4所示.

圖1 二值化圖

圖2 背景分割圖

圖3 邊緣提取

圖4 ATPC,PTPC和HP隨時(shí)間變化曲線

生菜ATPC，PTPC和HP的單位分別為cm2，cm和cm，計(jì)算公式為

(1)

(2)

(3)

式中：NL為俯視圖像中生菜區(qū)域的總像素?cái)?shù)；NP為俯視圖像中生菜外邊緣的總像素?cái)?shù)；NH為主視圖像中生菜高度方向的總像素?cái)?shù)；f1為比例因子，表示俯視圖像中1 cm2內(nèi)所含像素?cái)?shù)量；f2為比例因子，表示主視圖像中1 cm長(zhǎng)度內(nèi)包含的像素?cái)?shù)量.

2 結(jié)果與分析

2.1 養(yǎng)分質(zhì)量分?jǐn)?shù)、長(zhǎng)勢(shì)及產(chǎn)量間灰色關(guān)聯(lián)分析

根據(jù)灰色關(guān)聯(lián)分析理論，構(gòu)建養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)信息間灰色系統(tǒng)，其中比較數(shù)列為生菜冠層氮、磷、鉀質(zhì)量分?jǐn)?shù)，參考數(shù)列為各時(shí)刻長(zhǎng)勢(shì)信息，得到r1-r3為養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)信息的關(guān)聯(lián)度.通過(guò)關(guān)聯(lián)度大小確定影響強(qiáng)度.計(jì)算結(jié)果如表2所示.影響生菜長(zhǎng)勢(shì)變化的元素從大到小依次為氮，磷，鉀.這主要是由于氮元素參與生菜各種生命活動(dòng)，促進(jìn)蛋白質(zhì)和葉綠素等大分子及器官形成[11].磷元素在生菜體內(nèi)以有機(jī)和無(wú)機(jī)態(tài)的形式存在，有機(jī)磷是核酸和磷酸酯的必要組成部分，無(wú)機(jī)磷可激活酶的活性，幫助碳水化合物合成、分解和運(yùn)輸[12].鉀元素起到提高品質(zhì)的作用，增施鉀肥可以改善生菜的顏色及內(nèi)在品質(zhì)[13].因此，鉀元素對(duì)生菜長(zhǎng)勢(shì)影響要小于氮磷2種元素.對(duì)長(zhǎng)勢(shì)的影響從大到小依次為ATPC，PTPC，HP.這主要是由于當(dāng)缺肥時(shí)，生菜葉片組織失去膨壓，葉片伸展度下降，葉球緊實(shí)度降低，所以ATPC和PTPC比HP變化敏感.

表2 養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)間灰色關(guān)聯(lián)度

構(gòu)建長(zhǎng)勢(shì)與產(chǎn)量間灰色系統(tǒng)，其中比較數(shù)列為ATPC，PTPC，HP，參考數(shù)列為生菜產(chǎn)量，r1-r3分別為ATPC，PTPC，HP對(duì)產(chǎn)量的關(guān)聯(lián)度，計(jì)算后得到r1=0.748，r2=0.674，r3=0.707.得到ATPC，HP是影響生菜產(chǎn)量的主要因素.通過(guò)篩選影響作物產(chǎn)量的主要農(nóng)藝性狀，為作物高產(chǎn)培育提供了理論基礎(chǔ).

2.2 養(yǎng)分質(zhì)量分?jǐn)?shù)與產(chǎn)量關(guān)系研究

由上可知，與磷鉀元素相比，氮素是影響生菜長(zhǎng)勢(shì)和產(chǎn)量的主要因素.因此，文中著重分析磷鉀質(zhì)量分?jǐn)?shù)及產(chǎn)量在不同施氮條件下的變化規(guī)律.圖5為收獲期氮、鉀質(zhì)量分?jǐn)?shù)wK,wN對(duì)生菜產(chǎn)量Y影響等高線圖.結(jié)果表明，不同施氮量對(duì)生菜體內(nèi)鉀質(zhì)量分?jǐn)?shù)影響顯著，氮鉀元素間具有促進(jìn)吸收作用[14]，鉀質(zhì)量分?jǐn)?shù)隨著施氮量的增加而增加，但生菜產(chǎn)量最大值并非出現(xiàn)在氮鉀質(zhì)量分?jǐn)?shù)最大時(shí).

圖5 氮鉀質(zhì)量分?jǐn)?shù)對(duì)生菜產(chǎn)量影響等高線圖

圖6為收獲期不同氮、磷質(zhì)量分?jǐn)?shù)wP，wN對(duì)生菜產(chǎn)量影響等高線圖.研究結(jié)果表明，氮素質(zhì)量分?jǐn)?shù)較低時(shí)，生菜產(chǎn)量與含氮量存在正相關(guān)關(guān)系；氮素質(zhì)量分?jǐn)?shù)較高時(shí)，產(chǎn)量反而降低.隨著生菜體內(nèi)含氮量的增加，生菜的含磷量是逐漸升高，同時(shí)磷在作物體內(nèi)靠主動(dòng)吸收，環(huán)境因素溫度、光照等條件都會(huì)影響植物對(duì)磷的吸收.生菜最大產(chǎn)量達(dá)235.15 g，此時(shí)體內(nèi)氮磷鉀3種元素質(zhì)量分?jǐn)?shù)分別為4.83%，0.63%和2.82%.

圖6 氮磷質(zhì)量分?jǐn)?shù)對(duì)生菜產(chǎn)量影響等高線圖

2.3 養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)關(guān)系模型的建立

生菜長(zhǎng)勢(shì)特征與產(chǎn)量間具有一定的關(guān)聯(lián)性，而相較于產(chǎn)量，長(zhǎng)勢(shì)信息可利用機(jī)器視覺(jué)技術(shù)無(wú)損獲取.即通過(guò)長(zhǎng)勢(shì)特征間接反映生菜產(chǎn)量，實(shí)現(xiàn)對(duì)整個(gè)生長(zhǎng)期的跟蹤.因此，文中期望通過(guò)建立養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)的關(guān)系模型，探明養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)間的定量關(guān)系.由于環(huán)境溫度對(duì)作物的光合、呼吸、蒸騰、水肥的吸收影響顯著，故文中將環(huán)境溫度信息作為影響生菜生長(zhǎng)的主要環(huán)境因素.模型的一般形式為

(4)

式中：y為長(zhǎng)勢(shì)信息(ATPC，PTPC和HP)；β1—β11為回歸方程擬合系數(shù)；x1為環(huán)境溫度信息，℃；x2為生菜氮質(zhì)量分?jǐn)?shù)，%；x3為生菜磷質(zhì)量分?jǐn)?shù)(為使磷質(zhì)量分?jǐn)?shù)與氮鉀質(zhì)量分?jǐn)?shù)在同一數(shù)據(jù)尺度)，g/kg；x4為生菜鉀質(zhì)量分?jǐn)?shù)，%.

由圖4可知，作物在生長(zhǎng)過(guò)程中存在“慢-快-慢”S型生長(zhǎng)規(guī)律[15-16]，故文中將分3個(gè)階段建模，即生長(zhǎng)初期、速生期和緩慢生長(zhǎng)期.第一階段生菜生長(zhǎng)量較小，生物量累積速率不及養(yǎng)分吸收速率，養(yǎng)分質(zhì)量分?jǐn)?shù)逐漸上升；第三階段，生菜生長(zhǎng)緩慢，為抽薹開(kāi)花準(zhǔn)備，葉片養(yǎng)分吸收量逐漸下降.因此，這兩個(gè)階段養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)間存在線性相關(guān)，但各養(yǎng)分質(zhì)量分?jǐn)?shù)間存在交互作用，故對(duì)這2個(gè)階段的模型考慮自變量的一次項(xiàng)和交互項(xiàng).第二階段生物量迅速增加，養(yǎng)分吸收量在速生點(diǎn)附近出現(xiàn)拐點(diǎn)，故這個(gè)同時(shí)考慮了自變量的一次項(xiàng)、二次項(xiàng)和交互項(xiàng).養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)的關(guān)系模型如表3所示.

3 結(jié) 論

通過(guò)機(jī)器視覺(jué)技術(shù)連續(xù)監(jiān)測(cè)生菜長(zhǎng)勢(shì)情況，獲取冠層主視和俯視圖像，從中提取了能直接反應(yīng)生菜冠層長(zhǎng)勢(shì)狀況的ATPC，PTPC和HP這3個(gè)特征，實(shí)現(xiàn)了生菜長(zhǎng)勢(shì)信息的連續(xù)無(wú)損獲取.

通過(guò)研究生菜養(yǎng)分質(zhì)量分?jǐn)?shù)、長(zhǎng)勢(shì)及產(chǎn)量三者間的灰色關(guān)系，定性分析了施肥對(duì)長(zhǎng)勢(shì)及產(chǎn)量的影響，得到不同營(yíng)養(yǎng)對(duì)生菜長(zhǎng)勢(shì)影響從大到小依次為氮、磷、鉀，由此進(jìn)一步分析了不同施氮條件下，生菜體內(nèi)的磷鉀質(zhì)量分?jǐn)?shù)及產(chǎn)量的變化規(guī)律，當(dāng)產(chǎn)量達(dá)到最大時(shí)，氮磷鉀3種元素質(zhì)量分?jǐn)?shù)分別為4.83%，0.63%和2.82%.

由于長(zhǎng)勢(shì)與產(chǎn)量具有一定相關(guān)性，長(zhǎng)勢(shì)信息對(duì)產(chǎn)量影響從大到小依次為ATPC，HP，PTPC，其中ATPC和HP是決定生菜產(chǎn)量的主要因素，在實(shí)際生產(chǎn)中長(zhǎng)勢(shì)信息可利用圖像技術(shù)無(wú)損連續(xù)獲取，因此進(jìn)一步結(jié)合環(huán)境信息和生菜生長(zhǎng)規(guī)律，建立了不同生長(zhǎng)時(shí)期養(yǎng)分質(zhì)量分?jǐn)?shù)與長(zhǎng)勢(shì)的關(guān)系模型，為指導(dǎo)設(shè)施生菜科學(xué)施肥提供了理論基礎(chǔ).