*， 陳佰鴻

(1.甘肅農(nóng)業(yè)大學(xué)機(jī)電工程學(xué)院，甘肅 蘭州 730070；2.甘肅農(nóng)業(yè)大學(xué)園藝學(xué)院，甘肅 蘭州 730070)

草莓原產(chǎn)自南美，全世界共有20 000多個(gè)品種[1]，是一種喜溫涼和光的植物。草莓根系分布淺、蒸騰量大，對(duì)水分要求嚴(yán)格，但草莓在不同生長(zhǎng)期對(duì)水分的要求又稍有不同，在草莓果實(shí)生長(zhǎng)和成熟期對(duì)水分的需求最多，達(dá)80%以上。作物水分大部分由葉片承擔(dān)，缺水則會(huì)導(dǎo)致作物葉片萎縮、面積減小、顏色短期加深等情況[2]。目前草莓灌溉都以漫灌為主，導(dǎo)致大量的水資源浪費(fèi)，研究草莓按需灌溉是節(jié)約水資源的措施之一。

圖像處理技術(shù)簡(jiǎn)便、快速、無(wú)接觸、無(wú)破壞性，被認(rèn)為是現(xiàn)代無(wú)損檢測(cè)的重要技術(shù)之一[3-7]。一些學(xué)者利用圖像特征對(duì)部分作物水分檢測(cè)進(jìn)行了研究，主要集中在棉花、小麥、葡萄、玉米等作物的檢測(cè)上。王方永等[8]分析并發(fā)現(xiàn)棉花顏色參數(shù)G-R與水分含量及水分含量指數(shù)相關(guān)性最好，以G-R建立了棉花水分含量及水分含量指數(shù)預(yù)測(cè)模型，預(yù)測(cè)精度達(dá)到90.71%和91.02%。Zakaluk等[9]以馬鈴薯植株圖像預(yù)測(cè)植株水勢(shì)，利用主成分分析法發(fā)現(xiàn)，土壤氮與葉片反射率呈顯著的負(fù)線性關(guān)系(r=-0.71，p=0.003)。江朝暉等[10]采用偏最小二乘回歸建立了越冬期小麥冠層圖像特征與冠層含水率檢測(cè)模型，檢測(cè)得到冬小麥品種的相對(duì)誤差均值僅為1.29%。勞東青等[11]以溫室大棚中的葡萄植株為研究對(duì)象，采用多項(xiàng)式曲線擬合的方式建立了一種圖像灰度均值的含水率估算模型，模型的相關(guān)系數(shù)達(dá)到70.83%。韓文霆等[12]采用信息灰度共生矩陣和灰度直方圖分析抽穗期玉米葉片與含水率的關(guān)系，研究得到葉片顏色特征均值與含水率的關(guān)系模型決定系數(shù)為70.17%。Kacira[13]和Takayama等[14]建議把西紅柿葉子的垂直投影面積作為水分虧缺指數(shù)。

本文對(duì)圖像特征與結(jié)果期草莓葉片含水率的關(guān)系進(jìn)行了研究，提取葉片圖像的RGB、HSV和L*a*b*顏色模型中的多種顏色特征以及葉片投影面積，對(duì)32種顏色特征以及投影面積與草莓葉片含水率的關(guān)系進(jìn)行相關(guān)性回歸分析，從中篩選出相關(guān)性較高的顏色特征建立草莓葉片水分的預(yù)測(cè)模型，由于圖像傳感器的成本低，本文的研究結(jié)果可為草莓灌溉的水分監(jiān)測(cè)提供一種準(zhǔn)確率和實(shí)用性都較高的估算方法。

1 試驗(yàn)材料與方法

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

本次試驗(yàn)在甘肅農(nóng)業(yè)大學(xué)進(jìn)行，時(shí)間為2019年6月，供試對(duì)象為盆栽草莓，品種為寧玉，土壤質(zhì)地為草炭和珍珠巖接近一比一的混合。試驗(yàn)對(duì)象是生長(zhǎng)狀態(tài)相同的30盆正常生長(zhǎng)的草莓，其中20盆作為建模組，10盆作為驗(yàn)證組。對(duì)于建模組的草莓又隨機(jī)分為3組，第1組作為對(duì)比組，對(duì)這5盆草莓進(jìn)行每日定時(shí)澆水至飽和狀態(tài)，第2組也有5盆，該組試驗(yàn)期間不做澆水處理。第1、2組每天拍照，但不測(cè)量其葉片含水率，這樣做的目的是可以連續(xù)觀察草莓葉片隨水分減少時(shí)圖像特征的變化情況。由于本文葉片含水率測(cè)量時(shí)會(huì)剪掉一盆草莓的部分葉片，對(duì)于連續(xù)監(jiān)測(cè)葉片圖像特征的草莓無(wú)法測(cè)量其含水率，故本文采用了替代法，即對(duì)于建模組其余10盆作為第3組，該組不間斷地測(cè)量葉片含水率作為第2組對(duì)應(yīng)的含水率。對(duì)于該組的草莓，試驗(yàn)開始后前9盆每天測(cè)定一次，每次一盆，最后一盆是在試驗(yàn)開始后的第14天測(cè)量(此時(shí)草莓因缺水已整體倒伏)。對(duì)第2、3組草莓試驗(yàn)前一次性灌溉至飽和，隨后再不澆水，灌溉之后在這花盆土壤表面鋪一層稀碎紙片，以防水分流失太快不利于試驗(yàn)進(jìn)行。驗(yàn)證組水處理方法與建模組相同，每天測(cè)量一盆草莓的圖像特征與含水率，首先拍照，然后測(cè)量含水率，持續(xù)10天。

1.2 含水率測(cè)定

1.2.1 葉片含水率測(cè)定

在建模組的第3組中隨機(jī)挑選1盆草莓，均勻剪下10片葉子，并放在鋁箔中使用精度為0.000 1 g的分析天平稱重，記錄鮮重x1，然后將葉片平鋪于A4紙上，放置在型號(hào)為YHG-300-S的遠(yuǎn)紅外快速恒溫干燥箱內(nèi)進(jìn)行干燥，干燥溫度為80 ℃，干燥時(shí)間約為1 h，干燥過(guò)程中每15 min稱重一次，稱重到葉片重量誤差小于0.000 4 g時(shí)記錄此時(shí)的重量x2，則可通過(guò)計(jì)算得到葉片含水率y葉。葉片投影面積相對(duì)值、葉片含水率、土壤含水率變化趨勢(shì)如圖1所示，由圖1可以看到第4天的葉片含水率有上升趨勢(shì)，第5天到第7天緩慢遞減，隨后呈現(xiàn)出明顯遞減的現(xiàn)象。

圖1 葉片投影面積相對(duì)值、葉片含水率、土壤含水率變化趨勢(shì)

1.2.2 土壤含水率測(cè)定

在測(cè)完葉片含水率的草莓花盆中，從表面到根部均勻取土壤試樣，并放在鋁箔中使用分析天平稱重，記錄此時(shí)重量z1，然后放入預(yù)設(shè)溫度為105 ℃的干燥箱中干燥，干燥時(shí)間為4 h，記錄干燥后的重量為z2，則可得到土壤含水率y土。由圖1可知土壤含水率均勻遞減，并且第9天到第14天的水分流失最明顯。

1.3 圖像獲取

圖像采集使用D1000-IR-120/Color雙目攝像儀，分辨率為1 280*480，該儀器是一種可以全景曝光并應(yīng)用在Windows端的智能攝像儀器。拍攝在室內(nèi)進(jìn)行，使用太陽(yáng)自然光源，拍攝時(shí)間為早上9：00～12：00，采集的圖像以png格式儲(chǔ)存在計(jì)算機(jī)中。拍攝時(shí)相機(jī)鏡頭朝下，鏡頭光軸垂直于地面。拍攝過(guò)程中將48色比色板(型號(hào)為SpyderCHECKR48)平鋪于盆栽草莓旁邊，調(diào)整位置使比色板完整顯示在雙目鏡頭中，為后期圖像處理提供方便。固定拍攝時(shí)間與位置，連續(xù)拍照觀察直至草莓葉片全部萎蔫。

1.4 圖像處理

1.4.1 顏色校正

由于拍攝環(huán)境是室內(nèi)，各個(gè)花盆的拍攝時(shí)間也不同，并且實(shí)驗(yàn)進(jìn)行期間天氣不同，所以會(huì)因?yàn)楣饩€問(wèn)題導(dǎo)致拍攝的圖像顏色出現(xiàn)人眼所分辨不了的差異，因而使用圖像處理軟件Lightroom5.7和圖像校正軟件Spydercheckr對(duì)所有圖像進(jìn)行顏色校正。校正前后的圖像樣例如圖2和圖3所示。

1.4.2 圖像分割

圖像顏色校正后，對(duì)圖像進(jìn)行分割處理。由于草莓葉片邊緣有細(xì)小鋸齒，人工精細(xì)分割圖像非常耗時(shí)，本文采用相關(guān)文獻(xiàn)[15]介紹的算法進(jìn)行圖像分割，該算法先通過(guò)目標(biāo)檢測(cè)方法檢測(cè)到每張葉片，然后采用圖割方法優(yōu)化葉片顏色的能量函數(shù)從而實(shí)現(xiàn)葉片分割。實(shí)驗(yàn)發(fā)現(xiàn)該算法能細(xì)致地分割出葉片邊緣。對(duì)圖3分割后得到的圖像如圖4所示。

1.5 葉片投影面積的提取

用上述分割算法產(chǎn)生的分割模板自動(dòng)計(jì)算葉片面積，由于拍攝時(shí)鏡頭垂直朝下，故該面積為葉片的垂直投影面積S。圖1給出了第1組和第2組葉片投影面積的均值在試驗(yàn)期間的變化趨勢(shì)，為了作圖方便，該圖中分別將第一天5盆葉面積均值作為基準(zhǔn)，試驗(yàn)期間其他各天測(cè)得的投影面積均值除以基準(zhǔn)值作為相對(duì)投影面積，以下分析中所述面積均為相對(duì)面積?？梢钥闯龅?組由于每日定時(shí)澆水至飽和狀態(tài)，所以葉面積幾乎沒有變化。而第2組在剛開始的3天里，花盆有充足水分，葉面積幾乎沒有變化，其后隨著水分流失，葉面積開始逐漸減小。

圖2 葉片原圖

圖3 葉片校正圖

圖4 分割后葉片圖

1.6 顏色特征提取

在以往的研究文獻(xiàn)中，研究者使用了多種顏色特征建立不同植物葉片含水率模型[16-21]，這些特征來(lái)自于RGB、HSV和L*a*b*顏色模型。RGB是最基本的顏色模型，其中R、G和B分別是圖像的紅、綠和藍(lán)色分量，與人對(duì)顏色的主觀感知相一致，不足之處在于光照強(qiáng)度變化時(shí)各分量值也隨之變化。HSV是由色相(H)、飽和度(S)和亮度(V)組合的一種色彩模式，可將亮度與顏色分離。L*a*b*是一種人眼可以看到所有色彩的色彩模式，L*是光度分量，a*、b*為2個(gè)色度分量，a*表示從綠到紅，b*表示從藍(lán)到黃。前人研究表明與不同植物含水率關(guān)系最為密切的顏色特征是不同的，本試驗(yàn)收集32種與葉片含水率相關(guān)的顏色特征，分別對(duì)其與葉片含水率進(jìn)行回歸分析，尋找最適合表征草莓葉片含水率的顏色特征，見表1。本文對(duì)分割后的葉片圖像按像素求出平均的R、G、B值，然后利用MATLAB編程求出其他顏色特征的值。

2 數(shù)據(jù)分析

2.1 葉片含水率與顏色特征的相關(guān)性分析

為了全面研究葉片含水率與葉片圖像顏色特征的關(guān)系，提取了校正后葉片圖像的32種顏色特征，并且通過(guò)回歸分析發(fā)現(xiàn)每一個(gè)特征對(duì)目標(biāo)值都有不同程度的影響，結(jié)果見表2。

表1 顏色特征

續(xù)表1

表2 葉片顏色特征與葉片含水率的關(guān)系

由表2可知：在這32種顏色特征中，RGB顏色模型的(G-R)/(B-G)、G-R、R-B、r、R/(G+B)和L*a*b*顏色模型的L*與葉片含水率之間具有顯著的線性相關(guān)性，特征(G-R)/(B-G)在顯著性水平0.01時(shí)達(dá)到了極顯著性檢驗(yàn)水平，相關(guān)系數(shù)為0.769。(G-R)/(B-G)的校正測(cè)定系數(shù)為0.541，表明特征(G-R)/(B-G)能說(shuō)明葉片含水率的54.1%，而其余45.9%需要其他特征來(lái)解釋。

對(duì)與葉片含水率相關(guān)性較高的特征(G-R)/(B-G)、G-R、R-B、r、R/(G+B)、L做進(jìn)一步方差分析，見表3。

由表3分析結(jié)果可知，在這6種顏色特征中，特征(G-R)/(B-G)的F值最大，其顯著性p小于0.01，所以特征(G-R)/(B-G)對(duì)葉片含水率有非常顯著的影響，即可用特征(G-R)/(B-G)來(lái)建立葉片含水率預(yù)測(cè)模型。

2.2 葉片含水率與葉面積的相關(guān)性分析

本文主要分析第2組葉片投影面積S與葉片含水率之間的關(guān)系，變化趨勢(shì)如圖1所示，在花盆有充足水分的情況下，S幾乎沒有變化，其后隨著水分流失，S開始逐漸減小。S與葉片含水率的回歸結(jié)果見表4，p=0.01<0.05且相關(guān)系數(shù)達(dá)0.766，顯示S與葉片含水率之間具有顯著的相關(guān)性。

2.3 葉片含水率與葉片顏色特征和葉面積的雙變量相關(guān)性分析

使用同樣的方法對(duì)葉片顏色特征和S與葉片含水率之間的關(guān)系進(jìn)行雙變量回歸分析，發(fā)現(xiàn)這32種顏色特征和S與葉片含水率之間的關(guān)系相關(guān)，有23種組合特征與葉片含水率之間的關(guān)系密切相關(guān)，見表5。

表3 顏色特征與葉片含水率方差分析

表4 葉片面積和葉片含水率相關(guān)性分析

表5 顏色特征和葉面積與葉片含水率的關(guān)系

3 葉片含水率模型與驗(yàn)證

3.1 葉片含水率預(yù)測(cè)模型的建立

通過(guò)葉片含水率與32種顏色特征以及葉片投影面積S的相關(guān)性分析可知，在單變量回歸分析中，顏色特征(G-R)/(B-G)與葉片含水率的關(guān)系最為顯著，S次之，相關(guān)系數(shù)R分別為0.769、0.766。而在雙變量回歸分析中顯示，顏色特征G、B、H和S與葉片含水率相關(guān)性最為顯著，決定系數(shù)R2分別為0.839、0.831、0.824。本試驗(yàn)對(duì)顏色特征(G-R)/(B-G)和S分別建立單變量葉片水分預(yù)測(cè)模型，以及對(duì)G+S、B+S建立了雙變量葉片水分預(yù)測(cè)模型。

特征(G-R)/(B-G)、S與葉片含水率之間的散點(diǎn)圖分別如圖5、圖6所示，回歸方程分別見式(1)、式(2)。

圖5 特征(G-R)/(B-G)與葉片含水率線性擬合

圖6 葉面積與葉片含水率線性擬合

(1)

y2=0.272xS+0.468

(2)

由G+S、B+S與葉片含水率的回歸分析建立的雙變量葉片水分預(yù)測(cè)模型分別為式(3)和式(4)。

y3=0.273xS+0.155xG+0.311

(3)

y4=0.272xS+0.054xB+0.417

(4)

3.2 模型檢驗(yàn)

對(duì)建立的草莓葉片含水率預(yù)測(cè)模型進(jìn)行精度檢驗(yàn)，檢驗(yàn)樣本為驗(yàn)證組的10個(gè)數(shù)據(jù)，檢驗(yàn)結(jié)果見表6。

表6 回歸模型檢驗(yàn)

由表6 可知，雙變量預(yù)測(cè)模型相比于單變量預(yù)測(cè)模型預(yù)測(cè)效果好，其中特征B+S雙變量的葉片含水率預(yù)測(cè)模型的預(yù)測(cè)效果最好，擬合度達(dá)0.805，且平均相對(duì)誤差和均方根誤差值均小于10%。因此，可推薦以特征B+S作為草莓葉片水分含量的最佳預(yù)測(cè)模型。

4 結(jié)論

植物在生長(zhǎng)過(guò)程中普遍存在逆境生長(zhǎng)，而水分缺乏是反應(yīng)逆境生長(zhǎng)的一種形式，合理的含水量在植被功能等方面同樣發(fā)揮著重要作用，解決好水分缺乏的問(wèn)題，就是解決了植物光合作用等一系列問(wèn)題[22-24]。快速、無(wú)損、有效地檢測(cè)草莓的生長(zhǎng)狀況是現(xiàn)代精細(xì)管理的基礎(chǔ)，草莓葉片含水率能夠直觀地反映出草莓苗的生長(zhǎng)以及草莓果實(shí)的成長(zhǎng)情況。為快速檢測(cè)并診斷草莓的生長(zhǎng)情況，本文獲取草莓葉片圖像并測(cè)量草莓葉片含水率，分析了圖像RGB顏色模型、HSV顏色模型和L*a*b*顏色模型下的多種顏色特征以及葉片投影面積與葉片含水率之間的關(guān)系，得出如下結(jié)論：

(1)草莓葉片圖像經(jīng)過(guò)處理后得到的顏色特征與葉片含水率相關(guān)性顯著，RGB顏色模型的特征(G-R)/(B-G)與葉片含水率之間的關(guān)系尤為顯著，相關(guān)系數(shù)達(dá)0.769。草莓葉片投影面積與葉片含水率之間的關(guān)系顯著，相關(guān)系數(shù)達(dá)0.766。

(2)顏色特征和葉面積的雙變量與葉片含水率之間的關(guān)系最為顯著，選擇相關(guān)系數(shù)最大的兩組建立葉片含水率預(yù)測(cè)模型，經(jīng)檢驗(yàn)擬合效果較好。

(3)實(shí)驗(yàn)中土壤含水率隨時(shí)間遞減比較明顯，考慮下一步分析草莓土壤含水率與葉片顏色特征之間的關(guān)系。

本文通過(guò)對(duì)32種草莓葉片顏色特征、葉面積以及顏色特征和葉面積的組合特征與葉片含水率進(jìn)行相關(guān)性分析，并篩選相關(guān)性顯著的特征建立模型，確定了最佳含水率預(yù)測(cè)模型，為草莓干旱的無(wú)損診斷和生長(zhǎng)監(jiān)測(cè)提供了理論依據(jù)。