基于高光谱成像技术的鲜枣裂纹的识别研究
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第34卷,第2期 光谱学与光谱分析Vol畅34,No畅2,pp532‐5372014年2月 SpectroscopyandSpectralAnalysisFebruary,2014
基于高光谱成像技术的鲜枣裂纹的识别研究
余克强1,赵艳茹1,李晓丽1,张淑娟2,何 勇1倡
1畅浙江大学生物系统工程与食品科学学院,浙江杭州 3100582畅山西农业大学工学院,山西太谷 030801
摘 要 裂纹是衡量鲜枣品质的重要指标之一,果皮裂纹加速鲜枣的腐烂,导致鲜枣货架期的缩短,严重降低鲜枣的经济价值。
采用高光谱成像技术在380~1030nm波段范围内对鲜枣裂纹的位置及大小信息特征进行快速识别。
选用偏最小二乘回归(PLSR)、连续投影法(SPA)和全波段图像主成分分析(PCA),得到鲜枣裂纹相关的敏感波段。
然后利用选取的鲜枣裂纹的敏感波段对建模集的132个样本建立最小二乘支持向量机(LS‐SVM)判别模型,并对预测集的44个样本进行判别。
对PLSR‐LS‐SVM,SPA‐LS‐SVM和PCA‐LS‐SVM判别模型采用ROC曲线进行评判,得出PLSR‐LS‐SVM模型对鲜枣裂纹定性判别的结果(area=1,std=0)最佳。
选取PLSR回归系数挑选出的5条鲜枣裂纹敏感波段(467,544,639,673和682nm)对应的单波段图像进行主成分分析,其中将主成分PC4的图像结合图像处理技术,最终识别出鲜枣裂纹的位置、大小信息。
结果表明,采用高光谱成像技术结合光谱图像处理可以实现鲜枣裂纹定性判别和定量识别的研究,为进一步开发相关仪器的研究提供理论方法和依据。
关键词 高光谱成像技术;鲜枣裂纹;定性判别;定量识别
中图分类号:S123;TP391畅4 文献标识码:A DOI:10畅3964/j畅issn畅1000‐0593(2014)02‐0532‐06
收稿日期:2013‐05‐06,修订日期:2013‐07‐18
基金项目:国家“十二五”科技支撑计划课题项目(2011BAD21B04)和国家自然科学基金项目(31071332)资助
作者简介:余克强,1986年生,浙江大学生物系统工程与食品科学学院博士生 e‐mail:yuke406336022@gmail畅com
倡通讯联系人 e‐mail:yhe@zju畅edu畅cn
引 言
鲜枣在成熟期由于降雨,果皮水溶性果胶的增加使得果皮生长应力降低,当果肉迅速膨胀,导致果皮出现裂纹[1]。
鲜枣的果皮开裂加速了鲜枣腐烂变质,降低果实的商用价值。
目前主要采用手工剔除裂枣,不仅费时费力,而且对人工挑选的要求较高。
因此,快速、准确地检测鲜枣裂纹对于提高水果品质具有重要的意义。
高光谱成像(HSI)集光谱和数字图像信息于一体。
目前,该技术广泛在应用于农产品检测领域,不仅可以定量检测农产品内部品质(成熟度[2]、坚实度[3]、SSC[4]等),还能对农产品的外部特征进行定性分析。
近些年,研究人员运用高光谱成像技术进行了大量的研究,Wang等运用高光谱成像技术检测鲜枣的外部缺陷,利用少数敏感波段建立缺陷判别模型[5];Cho等用高光谱荧光成像技术对樱桃西红柿的表皮损伤进行研究,取得较好的检测效果[6]。
本研究采用高光谱成像技术对鲜枣的裂纹进行定性和定量检测。
在获取鲜枣高光谱数据后,提取相关感兴趣区的平均光谱,并采用PLSR,SPA和图像PCA分析得到鲜枣裂纹敏感波段,建立鲜枣裂纹判别的LS‐SVM模型以及用ROC曲线对判别模型进行评判,最终选用最优模型下敏感波段对应的高光谱图像进行分析,以获取鲜枣裂纹的位置、大小信息特征。
1 实验部分
1畅1 鲜枣样本
实验所用的鲜枣(梨枣)采摘于山西太谷果园。
为避免实验中鲜枣表面距高光谱系统镜头的不同而造成的误差,本实验尽量选择大小和形状相近,裂纹特征比较明显的鲜枣作为样本,挑选176个鲜枣(85个裂纹鲜枣和91个正常鲜枣)进行表面去污处理后,进行高光谱数据获取。
1畅2 高光谱成像系统和图像的获取
高光谱成像系统(图1)主要部件包括相机、光谱仪(380~1030nm)、镜头、卤素灯线光源、步进电机控制的传送带、计算机、暗箱等。
为获得清晰、不失真的高光谱图像,需要对传送带的速度、镜头的曝光时间等参数进行调节,本研
究中传送带的移动速度为3畅5mm・s-1,曝光时间为0畅13s,样本与镜头的距离为330mm。
Fig畅1 Schematicdiagramofmaincomponents
ofthehyperspectralimagingsystem
1畅3 高光谱图像的校正
由于镜头中存在的暗电流,光源的强度在各波段分布不均,使得高光谱图像在光源强度较弱的波段含有较大噪音,因此需要对高光谱原始图像进行黑白校正[7]。
在相同的图像采集环境下,扫描标准白色校正板(反射率接近99畅9%)得到标定图像Rw,关闭光源(反射率接近0%)进行图像采集得到标定图像Rd,然后按照式(1)进行校正得到校正后图像(Ic)。
Ic=Rr-Rd
Rw-Rd×DN(1)式中,Rr为原始的高光谱图像,DN为光照强度的最高值(本系统为4095)。
所有高光谱图像数据的获取都是基于Spec‐tralCube软件平台,数据处理和分析是运用ENVI4畅6(ITTvisualinformationsolutions,Boulder,CO,USA),“TheUn‐scramblerX10畅1”(CAMOPROCESSAS,Oslo,Norway)和MatlabR2009a(TheMathWorks,Inc畅,Natick,MA,USA)软件。
1畅4 鲜枣裂纹敏感波段的选择方法
鲜枣的高光谱图像是由380~1030nm内512条波段中图像线性组合而成,提取少量的、带有大量图像信息的敏感波段可以简化鲜枣裂纹的判别模型。
本研究选取敏感波段的方法有:(a)PLSR回归系数,挑选系数局部极值处的波段作为敏感波段;(b)SPA变量选取方法,它是利用向量的投影分析,寻找含有最低限度的冗余信息组,并使变量之间的共线性达到最小,同时能大大减少建模所用变量的个数,提高建模的速度和速率[7];(c)图像PCA分析,该方法是多波段的一个线性变换,变换数据到一个新的坐标系统,以使数据的差异达到最大。
由于每个PC图像都是由原始数据中的各个波段在的图像经过线性组合而成,通过选取线性组合中的权重系数可以对最佳特征图像进行选择[8]。
1畅5 鲜枣裂纹判别模型的建立和评价
本研究选用的LS‐SVM建模方法,它是对经典SVM算法的一种改进,为求方程线性解而对经典SVM内复杂的二次优化问题改进,不仅提升了计算速度,而且简化了问题的复杂性,并能利用RBF核函数对输入变量的非线性和线性信息,实现鲜枣裂纹的判别。
LS‐SVM算法的实现过程见参考文献[9]。
受试者工作特性曲线(ROC)是反映敏感性和特异性连续变量的综合指标[10]。
在ROC曲线中,“area”代表位于ROC曲线下的面积,“std”表示剩余残差的标准偏差,“1‐specifici‐
ty”表示每一个阈值的1‐specificity,“sensitivity”表示每个阈值的敏感性。
假使“area”值为1,则这个分类器对于进行分类判别是非常有效地,假如“area”值为小于或等于0畅5,则这个分类器是失效的[11]。
一个好的分类器对应的ROC曲线应该尽量靠近单位方形的左上角[10]。
1畅6 图像识别鲜枣裂纹
由于利用光谱分析只能对鲜枣裂纹进行定性判别,为了获取鲜枣裂纹定量信息,本文采用图像处理技术对鲜枣裂纹的位置、大小等信息进行识别。
2 结果与讨论
2畅1 鲜枣正常区与裂纹区的光谱特征
选取鲜枣裂纹和正常枣对应位置处感兴趣区(5×5像素点),并求得平均光谱(如图2)。
鲜枣裂纹的平均光谱在680nm处的吸收峰,其主要是由鲜枣果肉细胞中叶绿素对光谱吸收引起的[12,13];960nm左右的吸收峰是由鲜枣所含水分对光谱吸收引起的
[14]。
Fig畅2 Meanspectraoffreshjujubewith
soundandcracksurface
2畅2 特征波段的提取
利用PLSR回归系数(如图3)中局部极值来选取特征波段,可得到5条特征波段(467,544,639,673和682nm);采用SPA方法得到3条特征波段(579,978和1023nm)。
另外,将全波段高光谱数据进行PCA分析,前7个主成分的累积贡献率达99畅4%。
图4是前7个主成分的单波段灰度图像,可以看出前6个主成分图像保留了原来图像主要信息,
335
第2期 光谱学与光谱分析
PC7图像开始为噪声。
而PC6图像中裂纹信息较清晰,比较适合裂纹的识别,故选用该图像的权重系数(图5)来挑选敏感波段,可以得出有关鲜枣裂纹的敏感波段(414,443,555,683和894
nm)。
Fig畅3
RegressioncoefficientsofPLSRmodel
Fig畅4 Single‐bandgrayimageof
first
sevenprincipalcomponents
Fig畅5 WeightedcoefficientsofPC6image
根据以上选取的鲜枣裂纹敏感波段,发现其大部分集中在叶绿素的敏感波段处[15],这为正常表皮(深红色)和裂纹(浅绿色)鲜枣的有效判别模型提供得力基础。
2畅3 LS‐SVM判别模型的建立和ROC曲线的判别
按照大约3∶1的比例,将所有鲜枣样本分为校正集(132个)和预测集(44个),将得到的敏感波段采用LS‐SVM建立鲜枣裂纹的判别模型,统计结果见表1。
从表1中可以得到,对于PLSR‐LS‐SVM和SPA‐LS‐SVM鲜枣裂纹判别模型的错判个数为0,PCA‐LS‐SVM模型的错判数为13。
利用ROC曲线分别对PLSR‐LS‐SVM,SPA‐LS‐SVM和PCA‐LS‐SVM模型进行评价以确定最佳鲜枣裂纹判别模型。
由图6评价模型的ROC曲线可以得到PLSR‐LS‐SVM鲜枣裂纹的判别模型(area=1,std=0)效果最好。
因此选用PLSR回归系数得到的敏感波段对图像进行鲜枣裂纹位置、大小信息的识别。
Table1 StatisticsresultsofLS‐SVMdiscriminantmodelsofcracksfreshjujube
Models
PLSR‐LS‐SVM
SPA‐LS‐SVM
PCA‐LS‐SVM
Discriminantnumber
444444Misjudgmentnumber
0
0
13Fig畅6 ROCcurveofthreediscriminantmodels
(a):ROCcurveofPLSR‐LS‐SVM;(b):ROCcurveofSPA‐LS‐SVM;(c):ROCcurveofPCA‐LS‐SVM
435光谱学与光谱分析 第34卷
3 基于特征波段下鲜枣裂纹的识别过程
3畅1 鲜枣裂纹的识别方法
用光谱模型对鲜枣的裂纹判别只能达到定性的判别效果,采用图像处理的方法可以进一步识别出鲜枣裂纹的位置,大小信息特征。
由图7为PLSR提取敏感波段下的单波段灰度图像。
除了467nm波段下的图像以外,其余波段对应的图像可以看到裂纹,但由于裂纹比较细长不易识别,因此对这5条波段
对应的图像进一步PCA分析。
前5个主成分的图像如图8所示,其累积贡献率为99畅1%。
PC4图像下的鲜枣裂纹最明显,选用PC4图像来识别鲜枣裂纹的位置、大小等信息特征。
图9为鲜枣裂纹识别的图像过程,将PC4图像(1),经掩膜得到图(3),然后阈值分割(图3,图4)过程可以得出(4)图,但发现裂纹处有白色噪声,采用连通区域面积法消除白色噪声后得到了(5),计算裂纹在视觉平面内面积Sc(crack)=
727。
Fig畅7 Single
‐bandgrayimagebasedonsensitivebandsofPLSRregressioncoefficients
Fig畅8 Single
‐bandgrayimageoffirstfiveprincipalcomponents
Fig畅9 Processofimageidentificationofcrackofjujube
对鲜枣掩膜后的图像采用“sobel”算子进行边缘检测,可得到如图10中的鲜枣轮廓图(1),也可得鲜枣在视觉平面内的面积Sw(whole)=46533,计算裂纹占整个鲜枣的面积比为λ=Sc/Sw=1畅56%。
将鲜枣裂纹位置图(2)与鲜枣实物图(3)对比,
识别出裂纹位置与实物一致。
Fig畅10 Identificationofcrack’slocationandestimationarea3畅2 鲜枣裂纹识别的验证
为了验证鲜枣裂纹识别方法的可靠性,用同样方法对8个裂纹鲜枣进行识别,识别的结果和面积统计分析如图11和表2所示。
根据结果(图11和表2)可以看出大部分的裂纹都能被准确检测,但是对于细小的裂纹(g)和反光较强(h)的样本,鲜枣裂纹位置和大小有一定误差。
由于鲜枣呈椭球状,对于实际的相关面积测量和验证需要进行进一步研究,本文初步实现了对鲜枣裂纹的位置识别,以及鲜枣裂纹大小的估算,为以后的鲜枣裂纹定量研究奠定了基础。
4 结 论
采用高光谱成像技术对鲜枣裂纹进行定性判别和定量检测。
运用PLSR回归系数、SPA和图像PCA来提取鲜枣裂纹敏感波段,采用LS‐SVM建立鲜枣裂纹判别模型以及采用ROC曲线评价判别模型,表明PLSR‐LS‐SVM模型对鲜枣裂纹判别效果最佳。
运用该模型敏感波段下的图像进行PCA
5
35第2期 光谱学与光谱分析
Fig畅11 CrackimageinformationofeightfreshjujubesTable2 Statisticsresultsofcrackinformationofeightfreshjujubes
SamplesSc(crack)Sw(whole)PercentageSamplesSc(crack)Sw(whole)Percentage(a1)271295680畅92(e1)343257461畅33(b1)
420308341畅36(f1)296336520畅88(c1)342376110畅91(g1)393253371畅55(b1)596
30186
1畅97
(h1)
640
22860
2畅80
分析,对最能表现鲜枣裂纹特征的PC4图像进行裂纹信息的识别,定量地检测出裂纹在鲜枣的位置、大小信息特征,并将鲜枣裂纹识别模型及算法应用于另外8个样本,同样也得
到了较好的识别效果,这为鲜枣裂纹的检测仪器开发提供理论基础和依据。
References
[1] LIUGuang‐sheng,LIUYuan‐yuan,MAShu‐ping,etal(刘光生,刘源源,马树平,等).JournalofShanxiAgricultureSciences(山西农
业科学),2012,40(2):150.
[2] RajkumarP,WangN,ElmasryG,etal.JournalofFoodEngineering,2012,108:194.[3] GabrielA,LeivaV,LuRF,etal.JournalofFoodEngineering,2013,115:91.
[4] MendozaF,LuRF,AfrianaDW,etal.PostharvestBiologyandTechnology,2011,62:149.[5] WangJ,NakanoK,OhashiS,etal.BiosystemsEngineering,2011,108:345.
635光谱学与光谱分析 第34卷
[6] ChoBK,KimMS,etal.PostharvestBiologyandTechnology,2013,76:40.[7] WuD,ShiH,WangSJ,etal.AnalyticChimicaActa,2012,726:57.
[8] LIJiang‐bo,RAOXiu‐qin,YINGYin‐bin,etal(李江波,饶秀勤,应义斌,等).TransactionsoftheChineseSocietyofAgriculturalEngineering(农业工程学报),2010,26(8):222.
[9] WuD,SunDW.Talanta,2013,111:39.
[10] ZOUHong‐xia,QINFeng,CHENGZe‐kai,etal(邹洪侠,秦 锋,程泽凯,等).ComputerTechnologyandDevelopment(计算机技术
与发展),2009,19(6):109.
[11] LiXL,NiePC,QiuZJ,etal.ExpertSystemswithApplications,2011,38:11149.
[12] HUANGWen‐qian,CHENLi‐ping,LIJiang‐bo,etal(黄文倩,陈立平,李江波,等).TransactionsoftheChineseSocietyofAgricultur‐alEngineering(农业工程学报),2013,29(1):272.[13] AbbottJA,LuR,UpchurchBL,etal.HorticulturalReview,1997,(20):1.
[14] ElMasryG,WangN,VigneaultC,etal.LWT‐FoodScienceandTechnology,2008,41:337.
[15] HeHJ,WuD,SunDW.InnovativeFoodScienceandEmergingTechnologies,2013,http://dx.doi.org/10.1016/j.ifset.2013.02.
009.
StudyonIdentificationtheCrackFeatureofFreshJujubeUsingHyperspectralImaging
YUKe‐qiang1,ZHAOYan‐ru1,LIXiao‐li1,ZHANGShu‐juan2,HEYong1倡
1.CollegeofBiosystemsEngineeringandFoodScience,ZhejiangUniversity,Hangzhou 310058,China2.CollegeofEngineering,ShanxiAgriculturalUniversity,Taigu 030801,China
Abstract Crackisoneofthemostimportantindicatorstoevaluatethequalityoffreshjujube.Cracknotonlyacceleratesthede‐cayoffreshjujube,butalsodiminishestheshelflifeandreducestheeconomicvalueseverely.Inthisstudy,thepotentialofhy‐perspectralimagingcoveredtherangeof380~1030nmwasevaluatedfordiscriminationcrackfeature(locationandarea)offreshjujube.Regressioncoefficientsofpartialleastsquaresregression(PLSR),successiveprojectionanalysis(SPA)andprinci‐palcomponentanalysis(PCA)basedfull‐bandsimagewereadoptedtoextractsensitivebandsofcrackoffreshjujube.Thenleast‐squaressupportvectormachine(LS‐SVM)discriminantmodelsusingtheselectedsensitivebandsforcalibrationset(132samples)wereestablishedforidentificationthepredictionset(44samples).ROCcurvewasusedtojudgethediscriminantmod‐elsofPLSR‐LS‐SVM,SPA‐LS‐SVMandPCA‐LS‐SVMwhichareestablishedbysensitivebandsofcrackoffreshjujube.TheresultsdemonstratedthatPLSR‐LS‐SVMmodelhadanoptimaleffect(area=1,std=0)todiscriminatecrackfeatureoffreshjujube.Next,imagescorrespondingtofivesensitivebands(467,544,639,673and682nm)selectedbyPLSRwereexecutedtoPCA.Finally,theimageofPC4wasemployedtoidentifythelocationandareaofcrackfeaturethroughimagingprocessing.Theresultsrevealedthathyperspectralimagingtechniquecombinedwithimageprocessingcouldachievethequalitativediscriminationandquantitativeidentificationofcrackfeatureoffreshjujube,whichprovidedatheoreticalreferenceandbasisfordevelopinstru‐mentofdiscriminationofcrackofjujubeinfurtherwork.
Keywords Hyperspectralimaging;Cracksoffreshjujube;Qualitativediscrimination;Quantitativeidentification
(ReceivedMay6,2013;acceptedJul.18,2013)
倡Correspondingauthor
7
35第2期 光谱学与光谱分析
基于高光谱成像技术的鲜枣裂纹的识别研究
作者:余克强, 赵艳茹, 李晓丽, 张淑娟, 何勇, YU Ke-qiang, ZHAO Yan-ru, LI Xiao-li , ZHANG Shu-juan, HE Yong
作者单位:余克强,赵艳茹,李晓丽,何勇,YU Ke-qiang,ZHAO Yan-ru,LI Xiao-li,HE Yong(浙江大学生物系统工程与食品科学学院,浙江 杭州,310058), 张淑娟,ZHANG Shu-juan(山西农业大学工学
院,山西 太谷,030801)
刊名:
光谱学与光谱分析
英文刊名:Spectroscopy and Spectral Analysis
年,卷(期):2014(2)
本文链接:/Periodical_gpxygpfx201402055.aspx。