高士 邱偉 彭佳紅

摘要：數(shù)字圖像處理技術(shù)已經(jīng)成為精準(zhǔn)農(nóng)業(yè)的重要手段之一。以一個季度苧麻盆栽葉片圖像數(shù)據(jù)為基礎(chǔ)，對苧麻幼葉、熟葉、老葉生長過程中苧麻葉片氮、磷肥不同營養(yǎng)元素變化量與生長過程中苧麻葉片圖像的顏色特征進(jìn)行提取研究。使用RGB顏色模型對苧麻生長過程葉片圖像顏色特征進(jìn)行提取，并將其特征值與氮磷含量變化進(jìn)行相關(guān)分析得到：G/（R+G+B）可以作為缺氮診斷標(biāo)準(zhǔn)，其值大于0.43時，推斷苧麻生長缺氮;R/（R+G+B）可以作為缺磷診斷標(biāo)準(zhǔn)，其值大于0.34時，診斷苧麻缺磷。

關(guān)鍵詞：苧麻;缺素;特征值

中圖分類號：TP391? ? ? ?文獻(xiàn)標(biāo)識碼：A

文章編號：1009-3044（2019）24-0199-02

開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

1 背景

數(shù)字圖像處理是通過計算機(jī)對圖像進(jìn)行去除噪聲、增強(qiáng)、復(fù)原、分割、提取特征等處理的方法和技術(shù)。其中，數(shù)字處理技術(shù)在苧麻缺素診斷過程中的應(yīng)用過程主要包括圖像預(yù)處理、特征值提取、診斷三個部分。圖像預(yù)處理常用分量法灰度化圖像與中值濾波法去噪聲?；叶然瘓D像的每個像素只需一個字節(jié)存放灰度值[1、2]。將其灰度化處理，使R、G、B的值歸一到某個特定的區(qū)間或相等的這個值就是灰度值，范圍一般在0-255之間[3]。將采集到的圖像中的R、G、B分量作為圖像的灰度值，公式如下：

[F1（i，j）=R（i，j）;F2（i，j）=G（i，j）;F3（i，j）=B（i，j）;]? ? ? ? ? ? ? ? ? ? ? （1）

其中，F(xiàn)1、F2、F3即計算后的灰度值，i、j為灰度坐標(biāo)系坐標(biāo)。

中值濾波是一種常用的、能有效抑制噪聲的非線性信號處理技術(shù)，常用于保護(hù)邊緣信息[4]，是一種經(jīng)典的平滑噪聲的方法，其輸出圖形為g（x，y）=med{f（x-k，y-l），（k，l∈W）}，其中，f（x，y），g（x，y）分別為原始圖像和處理后圖像?；驹硎前褦?shù)字圖像中一點(diǎn)的值用該點(diǎn)的一個鄰域中各點(diǎn)值的中值代替，使周圍的像素值接近的真實值，從而消除孤立的噪聲點(diǎn)。

苧麻葉片圖像的特征值是苧麻葉片圖像的物理特征，主要包括其結(jié)構(gòu)、紋理、顏色、形狀等。顏色模型是數(shù)據(jù)來描述圖像的典型代表，常見的模型是RGB顏色模型，提取其R、G、B特征值算法如圖1所示。

2 基于苧麻氮磷缺素葉片的RGB特征值診斷

通過提取到的RGB特征值與氮磷肥的含量變化相關(guān)性進(jìn)行診斷。

實驗組共分為七個小組，施肥過程為苧麻生長的幼葉、熟葉、老葉期三個時期，每個時期施肥兩次，每次約50和20KG/hm2。如：氮肥的施量在適度范圍內(nèi)呈線性增長，其他營養(yǎng)元素比例調(diào)配適當(dāng)，使氮元素為單一變量，對于其他元素同樣如此。這樣，保證苧麻葉片由缺素到營養(yǎng)元素足夠的水平范圍內(nèi)[5]，見表1，單位為KG/hm2。

對苧麻葉片圖像提取的21個顏色特征值，見表2。

2.1 基于苧麻葉片RGB特征值的缺氮診斷

在同一苧麻品種的不同生長時期，將提取到的苧麻葉片顏色特征數(shù)據(jù)表2產(chǎn)生的值與生長過程的含氮肥含量變化表1分別進(jìn)行相關(guān)分析，得到相關(guān)性關(guān)系見表3。

由表3可看出，在幼葉期間，所提取的顏色特征值與氮肥的含量相關(guān)性很高，且結(jié)果都呈顯著相關(guān)關(guān)系。且誤差范圍在極小范圍之內(nèi)，對結(jié)果不造成影響。得出結(jié)論：氮肥的含量跟計算機(jī)處理得到的顏色特征值是有相關(guān)性的，可以為后面的缺素診斷奠定基礎(chǔ)[6]。在整個苧麻葉片生長過程的氮肥含量變化過程中，提取到R、G、B分別跟氮肥的含量變化呈負(fù)相關(guān)性，G/（R+G+B）綠光標(biāo)準(zhǔn)值在整個生長過程中與氮肥的含量變化呈負(fù)相關(guān)關(guān)系，根據(jù)統(tǒng)計學(xué)原理，在數(shù)字圖像處理技術(shù)中，可以利用該綠光標(biāo)準(zhǔn)值的變化范圍來診斷苧麻生長過程中氮肥含量的情況[7-8]。為此，在苧麻的生長過程中記錄了G/（R+G+B）綠光標(biāo)準(zhǔn)值與氮肥含量的關(guān)系，得出散點(diǎn)圖2。

由圖2中看出，不同時期的N肥施量情況影響G/（R+G+B）值的變化。根據(jù)經(jīng)驗，在0-50kg氮肥的情況下是屬于缺素階段，G/（R+G+B）值的變化范圍在[0.43，0.45]區(qū)間，當(dāng)G/（R+G+B）值的變化范圍在小于0.43范圍時，可診斷苧麻的氮肥含量充足或有施肥過量。綜上，利用數(shù)字圖像處理技術(shù)在提取苧麻葉片RGB值后，計算不同時期的G/（R+G+B）綠光標(biāo)準(zhǔn)值，然后依據(jù)其值的范圍就可以推斷出苧麻的缺氮情況。當(dāng)綠光標(biāo)準(zhǔn)值在G/（R+G+B）>0.43時，診斷為缺氮;而G/（R+G+B）<0.43，供氮情況良好。

2.2 基于苧麻葉片RGB特征值的缺磷診斷

對不同時期磷肥含量的變化（見表1）數(shù)據(jù)值與苧麻葉片顏色特征值數(shù)據(jù)（見表2）進(jìn)行相關(guān)性分析，得到表4。

表4中，不同時期的磷肥施量與R/（R+G+B）標(biāo)準(zhǔn)值有明顯的負(fù)相關(guān)性，依據(jù)統(tǒng)計學(xué)原理，這表明R/（R+G+B）紅光標(biāo)準(zhǔn)比值可以作為磷肥含量的診斷指標(biāo)，同理，依據(jù)R/（R+G+B）標(biāo)準(zhǔn)比值隨磷肥的含量變化建立散點(diǎn)圖，推斷出磷肥供應(yīng)不足時R/（R+G+B）的范圍區(qū)間，進(jìn)行缺素診斷，散點(diǎn)圖見圖3。

由圖3中可以看到，不同時期的P肥施量情況影響R/（R+G+B）值的變化，在0-50kg磷肥的情況下，屬于缺素階段，此時R/（R+G+B）值的變化范圍在[0.34，0.37]區(qū)間，當(dāng)R/（R+G+B）值的在小于0.34范圍的時候，可以診斷苧麻的磷肥含量充足或有施肥過量的情況。綜上，利用數(shù)字圖像處理技術(shù)在提取苧麻葉片RGB值后，篩選出不同時期的R/（R+G+B）紅光標(biāo)準(zhǔn)值，計算后依據(jù)其值的范圍可以推斷出苧麻的缺磷情況如何[8]。當(dāng)紅光標(biāo)準(zhǔn)值在R/（R+G+B）>0.34時，可以診斷為供磷不足時期，R/（R+G+B）<0.34時，可以診斷為供磷情況良好。