劉 鑫, 馮 潔, 楊舒明

(云南師范大學(xué) 物理與電子信息學(xué)院, 云南 昆明 650500)

基于多光譜圖像選取馬鈴薯葉片的最佳彩色波段組合

劉 鑫, 馮 潔, 楊舒明

(云南師范大學(xué) 物理與電子信息學(xué)院, 云南 昆明 650500)

為了快速目測識別馬鈴薯葉片,利用Spectrocam多光譜相機(jī)獲取健康C-88馬鈴薯葉片的可見光及近紅外通道的多光譜圖像。采用波段指數(shù)法提取葉片的特征波段,通過真彩色原理及標(biāo)準(zhǔn)假彩色對所提取的特征波段進(jìn)行彩色合成,在此基礎(chǔ)上得到最佳彩色波段組合。實驗結(jié)果表明,用波段指數(shù)法提取多光譜圖像的特征波段進(jìn)行彩色合成,能快速獲取馬鈴薯葉片的最佳波段,并得到680,558,475 nm和800,680,558 nm為最佳彩色波段組合,為遙感圖像的目視解譯和更多有效信息的提取提供了可靠依據(jù)。

多光譜圖像; 最佳彩色波段組合; 馬鈴薯葉片; 波段指數(shù)法

引 言

馬鈴薯是全世界范圍內(nèi)一種重要的糧食作物,我國馬鈴薯產(chǎn)量居世界首位。目前許多學(xué)者對馬鈴薯進(jìn)行了研究,宋英博等[1]通過光譜儀診斷馬鈴薯葉片中對氮素敏感的波段,張小燕等[2]研究了馬鈴薯主要營養(yǎng)成分的近紅外光譜,徐明珠等[3]利用高光譜技術(shù)識別馬鈴薯葉片早疫病,并建立了識別模型。

作物的光譜特性包含著作物的很多重要信息,如何既能最大程度地有效利用光譜信息,又能較快地處理光譜數(shù)據(jù)成為一個研究熱點[4-5]。傳統(tǒng)光譜技術(shù)對植株的檢測相關(guān)性低,多光譜成像技術(shù)利用敏感波段能有效反映目標(biāo)光譜和空間特征的信息?；赟pectrocam多光譜相機(jī)和計算機(jī)軟件組成的多光譜成像系統(tǒng)可以同時從光譜維和空間維獲取被測目標(biāo)的信息,能夠解決傳統(tǒng)科學(xué)領(lǐng)域“成像無光譜”和“光譜不成像”的問題,其特點是實現(xiàn)簡單、圖譜合一[6-7]。

由于人眼對彩色敏感且分辨能力強(qiáng),多光譜數(shù)據(jù)往往需要進(jìn)行波段組合,合成彩色圖像以便于進(jìn)行目標(biāo)識別和圖像特征提取[8-10]。但多光譜數(shù)據(jù)波段較多,各波段之間既具有互補(bǔ)性,又存在著大量的數(shù)據(jù)冗余[11],因此對多光譜圖像進(jìn)行彩色合成時,有必要從眾多的波段中提取特征波段進(jìn)行最佳波段組合。

鑒于此,本文基于Spectrocam多光譜成像系統(tǒng)通過波段指數(shù)法提取馬鈴薯葉片的特征波段,并將特征波段進(jìn)行最佳波段組合合成假彩色圖像。

1 實驗原理

1.1波段指數(shù)法

波段指數(shù)[12]的定義為

(1)

式中:σi為第i個通道圖像的標(biāo)準(zhǔn)差,其值越大,表明此波段所含的信息量越豐富;Ri表示Rw和Ra之和,其中,Rw為第i波段與所在組內(nèi)其他波段相關(guān)系數(shù)的絕對值之和的平均值,Ra為第i波段與所在組以外的其他波段之間的相關(guān)系數(shù)的絕對值之和,其值越小,表明波段之間的獨立性越高,信息冗余度越低。設(shè)ρij為波段i與j之間的相關(guān)系數(shù),光譜數(shù)據(jù)分為k組,每組的波段數(shù)分別為n1,n2,…,nk,則

Ri=Rw+Ra

(2)

(3)

波段指數(shù)要求按照波段均方差越大,波段間相關(guān)系數(shù)的絕對值越小的特點進(jìn)行最佳波段選擇。波段指數(shù)綜合考慮了各波段的信息量及其相關(guān)性,因此,波段指數(shù)對特征波段的選取具有重要意義。

1.2彩色合成原理

彩色合成[13]有兩種方案:真彩色原理合成,即R、G、B分別賦予紅光、綠光、藍(lán)光波段;標(biāo)準(zhǔn)假彩色合成,即R、G、B分別賦予近紅外、紅光、綠光波段。彩色合成的關(guān)鍵是如何進(jìn)行波段組合,波段選擇要遵循以下原則:所選的三個波段的信息總量要大;所選的三個波段的相關(guān)性較弱;目標(biāo)地物類型要在所選的波段組合內(nèi)與其他地物具有很好的可區(qū)分性。

2 材料與方法

2.1多光譜成像系統(tǒng)概述

實驗采用多光譜成像系統(tǒng)進(jìn)行實時無損的馬鈴薯葉片多光譜圖像采集。實驗使用的儀器設(shè)備主要由海洋薄膜公司全新Spectrocam多光譜相機(jī)和計算機(jī)控制軟件組成。Spectrocam多光譜成像系統(tǒng)具有可互換式濾光片,選用濾光片中心波長分別為425(1),475(2),509(3),515(4),558(5),578(6),620(7),650(8),680(9),717(10),750(11),800(12),832(13),850(14),972(15)nm的15個通道,其對應(yīng)的半帶寬分別是100,100,20,10,5,10,10,10, 10,10,10,10,5,100,10 nm。

2.2材料制備及數(shù)據(jù)采集

實驗選取的健康馬鈴薯樣本來自云南師范大學(xué)馬鈴薯學(xué)院,馬鈴薯品種為C-88,生長時間為80 d,生長狀況良好。在溫室中分別采集盆栽培養(yǎng)的健康馬鈴薯葉片和參考白板的多光譜圖像,采集時間為12:00—14:00,多光譜相機(jī)與樣本之間的距離為1.2 m,參考白板用于測量時的系統(tǒng)標(biāo)定。圖1所示為15個通道下的多光譜圖像。

圖1 各通道下的多光譜圖像Fig.1 Multispectral images of each channel

如圖2所示,在所采集的多光譜圖像中選擇5個不同區(qū)域進(jìn)行灰度值提取,并同時提取參考白板的灰度值。

圖2 馬鈴薯葉片選取的區(qū)域Fig.2 The selected area of potato leaves

2.3實驗結(jié)果及分析

根據(jù)光譜反射率的定義得到健康馬鈴薯葉片的光譜反射比,通過MATLAB平臺擬合得到光譜反射率曲線,如圖3所示。

從圖3可以看出,在波長558 nm處出現(xiàn)一反射峰,在可見光波段與近紅外波段之間,即680～750 nm波段,反射率急劇上升且曲線陡,接近于直線的形狀,形成紅邊現(xiàn)象,這是典型綠色植物光譜曲線的明顯特征,故實驗數(shù)據(jù)滿足實驗分析要求。

根據(jù)波段指數(shù)法求出實驗中不同馬鈴薯葉片灰度值的標(biāo)準(zhǔn)差和不同通道之間的相關(guān)系數(shù),圖4為馬鈴薯葉片成像光譜數(shù)據(jù)的標(biāo)準(zhǔn)差,表1為不同通道之間的部分相關(guān)系數(shù)。

在利用多光譜數(shù)據(jù)進(jìn)行假彩色合成時,一般應(yīng)從近紅外、紅光、藍(lán)綠光波段中各選一個最優(yōu)波段進(jìn)行合成。為此,可將15個通道根據(jù)波長大小分為3組,即通道1～7(除紅光外的可見光組)、8～11(紅光組)、12～15(近紅外組)。在得到不同波段的標(biāo)準(zhǔn)差和相關(guān)系數(shù)后,利用式(1)計算出馬鈴薯葉片不同通道的波段指數(shù)值,表2為各組的波段指數(shù)值。

圖3 健康馬鈴薯葉片的光譜反射率曲線Fig.3 Spectral reflectance of healthy potato leaves

圖4 馬鈴薯葉片的各波段標(biāo)準(zhǔn)差曲線Fig.4 Standard deviation curve of potato leaves in every band

通道12345678910111213141511.0020.791.0030.880.851.0040.750.870.971.0050.15-0.060.420.431.0060.720.650.930.930.711.0070.710.730.950.970.600.971.0080.500.760.840.950.500.870.921.009-0.12-0.310.170.180.590.460.420.301.0010-0.11-0.070.350.340.670.520.560.390.891.00110.27-0.010.450.390.760.610.620.390.870.981.00120.15-0.090.380.340.770.560.570.370.880.990.991.00130.27-0.050.430.350.750.580.580.330.860.961.000.991.00140.20-0.080.410.350.800.600.590.370.910.971.000.990.991.00150.26-0.030.460.410.830.650.630.410.910.960.990.990.990.991.00

從表2得知可見光通道下馬鈴薯葉片在475,680,717,750 nm處波段指數(shù)較大,近紅外通道下馬鈴薯葉片在800 nm和850 nm處波段指數(shù)相對較大。根據(jù)植物光譜特征,在通道558 nm處存在一個反射峰,植物葉片在該波段附近信息相對豐富。綜合波段指數(shù)的特征波段選取原理和植物光譜特征,得出波段475,558,680,715,750,800,850 nm是可以獲得馬鈴薯葉片光譜信息的特征波段,即獲得馬鈴薯葉片特征光譜信息的最佳通道是2、5、9、10、11、12和14。

通過真彩色原理,即R、G、B分別賦予紅光、綠光、藍(lán)光波段,和標(biāo)準(zhǔn)假彩色即R、G、B分別賦予近紅外、紅光、綠光波段,分別對所提取的特征波段進(jìn)行彩色合成,合成方案如表3所示。

表2 馬鈴薯葉片各通道的波段指數(shù)值Tab.2 Band index of each channel of potato leaves

表3 合成方案Tab.3 Synthetic scheme

為了得到最佳波段組合,綜合考慮3個波段的信息總量及相關(guān)性,從表1可以看出相鄰?fù)ǖ?、10、11和通道5、2之間的相關(guān)性分別為-0.31、-0.07、-0.01和0.67、0.70、0.77,可見通道9的相關(guān)性較小,故選取通道9、5、2組合為最佳波段組合,同理得知標(biāo)準(zhǔn)假彩色方案中通道12、9、5為最佳波段組合。通過ENVI軟件將R、G、B分別賦予通道9、5、2和12、9、5合成假彩色圖像,如圖5所示。

從圖5可以看出,通道9、5、2合成的假彩色圖像與植物的真實顏色十分接近,通道12、9、5合成的假彩色圖像中植物顯示紅色。

圖5 假彩色圖像Fig.5 False color images

3 結(jié) 論

提出利用波段指數(shù)法提取健康馬鈴薯葉片的特征波段進(jìn)行假彩色合成并得到最佳波段組合。實驗結(jié)果得出通道475,558,680,717,750,800,850 nm是可以獲得馬鈴薯葉片光譜信息的特征波段,能快速獲取馬鈴薯葉片的信息,在此基礎(chǔ)上得出通道9(680 nm)、5(558 nm)、2(475 nm)和通道12(800 nm)、9(680 nm)、5(558 nm)為假彩色最佳波段組合。通道9、5、2合成的假彩色圖像與植物的真實顏色十分接近,由于人眼對彩色比較敏感且分辨能力強(qiáng),信息豐富的彩色合成圖像對有效進(jìn)行多光譜數(shù)據(jù)處理分析及信息提取至關(guān)重要。因為植物的輻射量差異在近紅外波段最敏感,通道12、9、5合成的假彩色圖像中植物顯示紅色,紅色植物非常容易和其他地物區(qū)分,在遙感植物研究中經(jīng)常用到。利用波段指數(shù)法提取特征波段選取假彩色最佳波段組合,此研究可為后續(xù)的圖像特征提取、目標(biāo)識別等工作提供參考。

[1] 宋英博.光譜診斷馬鈴薯葉片氮素敏感波段的研究[J].中國馬鈴薯,2010,24(3):176-178.

[2] 張小燕,楊炳南,劉威,等.馬鈴薯主要營養(yǎng)成分的近紅外光譜分析[J].食品科學(xué),2013,34(2):165-169.

[3] 徐明珠,李梅,白志鵬,等.馬鈴薯葉片早疫病的高光譜識別研究[J].農(nóng)機(jī)化研究,2016,38(6):205-209.

[4] 楊瑋,SIGRIMIS N,李民贊.基于多光譜圖像分析的溫室黃瓜葉片營養(yǎng)元素檢測與診斷[J].光譜學(xué)與光譜分析,2010,30(1):210-213.

[5] 馮潔,廖寧放,趙波,等.多光譜成像技術(shù)診斷植物病蟲害的人工神經(jīng)網(wǎng)絡(luò)模型[J].光學(xué)技術(shù),2008,34(5):717-720.

[6] 曹鵬飛,李宏寧,羅艷琳,等.基于多光譜成像選取四季豆葉片的特征波段[J].激光與光電子學(xué)進(jìn)展,2014,51(1):011101.

[7] 馮潔,李宏寧,楊衛(wèi)平,等.園藝作物病害的多光譜組合分類[J].光譜學(xué)與光譜分析,2010,30(2):426-429.

[8] 林文鵬,王長耀,儲德平,等.基于光譜特征分析的主要秋季作物類型提取研究[J].農(nóng)業(yè)工程學(xué)報,2006,22(9):128-132.

[9] 馮恩國,李淑香.居民地信息提取的最佳假彩色合成方法——以開封市區(qū)為例[J].平頂山學(xué)院學(xué)報,2006,21(2):91-94.

[10] 楊金紅,尹球,顧松山,等.城區(qū)高光譜遙感數(shù)據(jù)假彩色波段組合研究[J].南京氣象學(xué)院學(xué)報,2005,28(3):289-295.

[11] 張志敏,張?zhí)鞓?康冬舟,等.遙感圖像多光譜假彩色合成處理及其在環(huán)境地學(xué)中的應(yīng)用[J].影像技術(shù),2002(1):33-36.

[12] 姜小光,王長耀,王成.成像光譜數(shù)據(jù)的光譜信息特點及最佳波段選擇——以北京順義區(qū)為例[J].干旱區(qū)地理,2000,23(3):214-220.

[13] 吳云霞,史巍,史繼忠.多波段遙感圖像彩色合成處理解析——遙感圖像彩色合成處理教學(xué)設(shè)計[J].科技視界,2012(29):54-55.

Selectingbestcolorbandcombinationsforpotatoleavesbasedonmulti-spectralimages

LIU Xin, FENG Jie, YANG Shuming

(School of Physics and Electronic Information, Yunnan Normal University, Kunming 650500, China)

In order to identify the potato leaves rapidly,we use a Spectrocam multispectral camera to capture multispectral images of the healthy C-88 potato leaves in visible light and near infrared channels.The characteristics of the leaves are obtained by the band index method,and then color synthesis is carried out based on the principle of true color and false color.Using the above experiments,we obtain the best color band combinations.Experimental results show that extracting the characteristics of the leaves through the band index method and then finishing color synthesis can quickly yield the best color band combinations.The bands 680,558,475 nm and the bands 800,680,558 nm are the best color band combinations.The study provides reliable criteria for the visual interpretation of remote sensing images and the extraction of more effective information.

multispectral image; best color band combinations; potato leaf; band index method

1005－5630（2017）05－0035－06

2016-09-07

國家大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計劃(201510681005);云南省科技計劃(2016FB108)

劉 鑫(1993—),女,碩士研究生,主要從事新型光譜成像技術(shù)方面的研究。E-mail:lxin_66@126.com

馮 潔(1975—),女,教授,主要從事多光譜成像和顏色科學(xué)方面的研究。E-mail:fengjie_ynnu@126.com

O 439

A

10.3969/j.issn.1005-5630.2017.05.006

(編輯:張磊)