董 明,张 洋,张晓晓,朱嘉玮,杨卫平
(云南师范大学 物理与电子信息学院,云南 昆明 650500)
基于LCTF成像和ABS算法的小白菜光谱特征波段提取研究
董明,张洋,张晓晓,朱嘉玮,杨卫平
(云南师范大学 物理与电子信息学院,云南 昆明650500)
基于液晶可调滤波器(LCTF)和CMOS组合的多光谱成像系统,在435~720 nm波段范围内,以每隔5 nm波段对小白菜叶片进行灰度值信息的提取,然后求出各个波段的灰度值平均值、标准差以及相关系数,并采用自适应波段选择法(ABS)提取出小白菜叶片的波段指数,最后通过波段指数的排序选取出小白菜叶片的有效特征波段。实验结果表明,用ABS的特征波段提取的算法,能够快速有效地获取小白菜叶片的光谱信息,在445 nm、450 nm、455 nm、680 nm、685 nm、690 nm、695 nm和710 nm波段具有较理想的波段指数值,有较大的光谱信息量。因此,这些波段可以很好地作为识别小白菜叶片的有效特征信息波段。
液晶可调滤波器 (LCTF); ABS算法; 特征波段; 小白菜叶片; 波段指数
我国是农业大国,大面积种植着各类园艺作物[1],地域广阔、种类繁多,同时作物的产量和质量还会受到病虫害的影响[2-4]。因而如何快速、有效地对农作物进行准确的区分和识别[5],对我国农业的发展有着指导性的意义[6]。以往大多数是通过人眼视觉和主观经验,利用已知农作物的形状、颜色等来获取作物的识别信息,具有较大的主观性,往往带来了很大的缺陷和不足,造成大量的错误识别等。近年来光谱成像技术在农业中得到了广泛应用[7],而基于液晶可调滤波器(liquid crystal tunable filter,LCTF)[8-9]和单色CCD组成的多光谱成像系统,可以快速有效地提取作物的图像光谱信息,具有较高的光谱分辨率和空间分辨率[10],可以较好地对多光谱图像进行数据采集[11]。但由于多光谱成像系统具有较多的波段,会使得数据的冗余度大,数据处理、计算耗时长,造成数据处理工作的繁琐和准确率的降低。
因此,本研究采用一种基于LCTF光谱成像系统的自适应波段选择(ABS)算法,对我国大部分居民食用的十字花科类蔬菜中具有抗癌、护眼、摄取维生素作用的小白菜叶片进行特征波段的提取,从大量丰富的波谱信息中提取具有代表性的特征波段[12-13],主要用于对小白菜叶片的正确识别、聚类分析[14]和数据库的建立[15]等。这种方法对数据的处理达到了快速、准确的作用,是一种非接触、安全、实用、有效的方法。
ABS算法[16]对各个波段之间的空间相关性和谱间相关性大小做出了充分的考虑,并构造了相应的数学模型,对已求得的各个波段指数进行从大到小的顺序排列,并根据设定的阈值自适应地选择所需波段。该数学模型中的Indexi指数表示为:
(1)
式中:Indexi为第i幅图像的指数大小;σi为第i个波段的标准差;Ri-1,i为第i个波段与前一个波段的相关系数;Ri,i+1为第i个波段与后一个波段的相关系数,相关系数越小,两个波段数据之间的独立性越高。
(2)
式中:M、N分别为图像的行、列像素;fi为第i幅图像;μi为第i幅图像的像素平均值。
(3)
式中:Ri,j为相关系数;E{ }为数学期望。
ABS算法获得的波段指数充分考虑了每幅图像的信息量与相邻波段的相似性,Indexi越大,相应图像的信息量就越大,越具有代表性。将ABS算法得到的波段指数由大到小依次排序,排序靠前的7~8个波段即可作为该算法的有效特征波段。
2.1实验仪器
本实验使用美国CRI公司生产的VariSpecTM液晶可调滤波器,与单色CCD相机[17]和计算机控制软件组成的多光谱成像系统,其采集的波段范围为430~720 nm,带宽(FWHM)为10 nm,半角可视范围为7.5°,响应时间为50 ms,每间隔5 nm连续改变光谱透过中心波长,同时将采集到的小白菜叶片和参考白板的灰度图像传至计算机存储。本实验采用自行搭建的LCTF光谱成像系统,用于采集小白菜叶片的灰度图像,采用LCTF代替传统的机械滤光片转轮,以实现快速波段调谐。图1为基于LCTF光谱成像系统选取的有效特征波段原理图。本实验用于数据处理的软件为Omnic 8.0和MATLAB 2010b软件。
图1 基于LCTF光谱成像系统选取有效特征波段的原理图
2.2实验样本
实验选取广东省农科院的十字花科类-芸苔属健康小白菜(品种为“17号白菜”)作为实验样本,实验样本采集的环境是在温室大棚内且天气晴朗无风,光照强度均匀,平均气温为23.5 ℃,平均湿度为68%,小白菜生长时间为20~40 d,叶片成深绿色,叶脉清晰,生长状况良好。
2.3数据采集
新鲜生长的小白菜在特定的光源下和参考白板同时进行多光谱成像的图像采集,在435~720 nm波段范围内,每隔5 nm采集一幅灰度图像,共采集6组,每组58个波段,共计348组数据。以波段490 nm处的灰度图像为例说明选取小白菜叶片和参考白板区域,如图2所示。
图2 小白菜叶片和参考白板选取的区域
利用C#语言编写程序,获取样本图像的平均灰度值,在各波段选定的区域获取小白菜叶片和参考白板图像的灰度值。对每个波段样本的5个不同叶片表面选取1个矩形区域进行灰度值提取,相当于对所选取的每一个叶片的特定区域进行了灰度平均,消除了叶片表面亮度不均匀等特征因素的影响。本实验在暗室内进行,采用D65光源,由于光源和实验样本的影响导致参考白板有个较小的峰值,其他波段是较为理想的平滑曲线,属正常情况。实验室光源如图3所示,参考白板和小白菜叶片选取的区域的平均灰度值曲线如图4所示。
图3 光源
图4 平均灰度值曲线。
图5为小白菜叶片在435 nm~720 nm范围内的反射比曲线,由图可以清晰地看到小白菜叶片的光谱曲线在460 nm、550 nm和680 nm附近存在两个吸收谷和一个反射峰。从图中可以得出所采集到小白菜叶片的光谱数据符合一般绿色植物的曲线走势图,因此本实验数据满足实验分析的要求。
图5 小白菜叶片的反射比曲线
利用ABS算法从实验中获取小白菜叶片不同波段的成像光谱数据(相关系数和标准差),表1给出不同波段区间具有代表性的部分相关系数。在波段465 nm处相对于450 nm、565 nm、710 nm处的相对系数较小,而在635 nm和465 nm波段的相关系数达到了负值,其他波段之间的相关系数绝对值均在0.014-1.000之间。
表1 小白菜叶片的部分相关系数
利用式(1)~(3)计算出小白菜叶片不同波段的Index值,得到表2所示的指数排序表及图6所示的Index指数图。根据表2波段指数大小排序可以得出:在445 nm波段处具有最大的波段指数值1.355,其次分别是695 nm、680 nm、690 nm、710 nm、450 nm、455 nm和685 nm处。综合计算,上述波段中所含的光谱信息量较为丰富,故选择445 nm、450 nm、455 nm、680 nm、685 nm、690 nm、695 nm、710 nm波段作为本实验中小白菜叶片的整体最优特征波段。
表2 小白菜叶片成像光谱数据波段指数排序
图6 小白菜叶片Index指数杆状图
本文以健康小白菜叶片的多光谱成像为例,利用ABS算法计算实验数据并讨论光谱特征,得出445 nm、450 nm、455 nm、680 nm、685 nm、690 nm、695 nm、710 nm波段能够更好地反映健康小白菜叶片的特征光谱信息,可以作为小白菜叶片的特征波段。实验结果表明,基于LCTF组成的多光谱成像系统和ABS算法为健康小白菜叶片所提取的特征波段,能够为农作物的病虫害监测、识别提供一种新的方法,同时该方法对地物的高光谱遥感识别、探测等也具有参考作用。
[1]冯洁,李宏宁,杨卫平,等.园艺作物病害的多光谱组合分类[J].光谱学与光谱分析,2010,30(2):426-429.
[2]LUCAS B G,CAMPBELL C L,Lucas L T.Introduction to plant diseases:identification and management[M].US:Springer,1992.
[3]王坤,朱大洲,张东彦,等.成像光谱技术在农作物信息诊断中的研究进展[J].光谱学与光谱分析,2011,31(3):589-594.
[4]刘涛,孙旭东,刘燕德,等.农作物品质的近红外光谱无损检测研究进展[J].食品与机械,2010,26(3):161-166.
[5]刘波,方俊永,刘学,等.基于成像光谱技术的作物杂草识别研究[J].光谱学与光谱分析,2010,30(7):1830-1833.
[6]林文鹏,王长耀,储德平,等.基于光谱特征分析的主要秋季作物类型提取研究[J].农业工程学报,2006,22(9):128-132.
[7]范世福.光谱技术和仪器的新发展[J].光学仪器,2000,22(4):35-40.
[8]ZJAKIC I,PARAC-OSTERMAN D,BATES I.New approach to metamerism measurement on halftone color images[J].Measurement,2011,44(8):1441-1447.
[9]杜丽丽,易维宁,张冬英,等.基于LCTF的多光谱面阵CCD相机的辐射定标[J].光谱学与光谱分析,2011,31(1):272-276.
[10]童庆禧,张兵,郑兰芬.高光谱遥感-原理、技术与应用[M].北京:高等教育出版社,2006:40-44.
[11]李珂,张康伟,罗淼.基于LCTF的大幅面高分辨率多光谱仪光学系统设计[J].红外与激光工程,2013,42(3):675-679.
[12]王庆光.多光谱遥感数据最佳波段选择研究[J].广东水利电力职业技术学院学报,2007,5(2):63-65.
[13]杨金红.高光谱遥感数据最佳波段选择方法研究[D].南京:南京信息工程大学,2005.
[14]田方,樊治飞,焦多礼,等.甘肃药用植物地理分布的聚类分析及分区[J].甘肃科技,2011,27(6):55-57.
[15]沈光,佟斌,孙波.论植物数据库的建立[J].国土与自然资源研究,2007(1):89-90.
[16]赵春晖,陈万海,杨雷.高光谱遥感图像最优波段选择方法的研究进展与分析[J].黑龙江大学自然科学学报,2007,24(5):592-602.
[17]郝学元,季鑫源,陈家胜.用线阵CCD实现光谱测试的虚拟仪器[J].光学仪器,2001,23(3):3-6.
(编辑:刘铁英)
The research of extracting the characteristic wavebands for pakchoi applying for ABS algorithm based on the LCTF imaging techniques
DONG Ming,ZHANG Yang,ZHANG Xiaoxiao,ZHU Jiawei,YANG Weiping
(School of Physics and Electronic Information,Yunnan Normal University,Kunming 650500,China)
Gray information of pakchoi leaves at the wavelength range of 435~720 nm with an interval of 5 nm are captured by using a multi-spectral imaging system which mainly consists of liquid crystal tunable filter (LCTF) and CMOS camera.Then,the average value,standard deviation and correlation coefficient for the gray value of images are calculated,and it has extracted the waveband index of pakchoi leaves by the way of adaptive band selection (ABS) method. Finally,the effective characteristic wavebands of pakchoi leaves are extracted through the sorting of waveband index.The experimental result shows that ABS algorithm using the extracted feature bands can quickly and efficiently obtain spectral information pakchoi leaves,and there are ideal waveband index values at 445 nm,450 nm,455 nm,680 nm,685 nm,690 nm,695 nm and 710 nm respectively,these wavebands have much discrete degree,rich multispectral information.therefore,these wavebands can be used as effective characteristic wavebands identification for pakchoi leaves.
liquid crystal tunable filter(LCTF); ABS algorithm; feature bands; pakchoi leaves; waveband index
2015-09-11
国家自然科学基金(61168003、60968001);云南省自然科学基金(2011FZ079、2009CD047);地方高校国家级大学创新创业训练计划(201310681004)
董明(1990—),男,硕士研究生,主要从事多光谱成像技术方面的研究。E-mail:dm704655044@163.com
杨卫平(1958—),男,教授,主要从事多光谱图像颜色真实再现等方面的研究。E-mail:yangwpkm@126.com
1005-5630( 2016) 03-0243-05
O439
A
10.3969/j.issn.1005-5630.2016.03.010