基于高光谱漫透射成像整体检测番茄可溶性固形物含量.pdf
《基于高光谱漫透射成像整体检测番茄可溶性固形物含量.pdf》由会员分享,可在线阅读,更多相关《基于高光谱漫透射成像整体检测番茄可溶性固形物含量.pdf(6页珍藏版)》请在咨信网上搜索。
1、第 29 卷 第 23 期 农 业 工 程 学 报 Vol.29 No.23 2013年 12月 Transactions of the Chinese Society of Agricultural Engineering Dec.2013 247 基于高光谱漫透射成像整体检测番茄可溶性固形物含量 张若宇1,2,饶秀勤1,3,高迎旺1,3,胡 栋1,3,应义斌1,3(1.浙江大学生物系统工程与食品科学学院,杭州 310058;2.石河子大学机械电气工程学院,石河子 832000;3.农业部设施农业装备与信息化重点实验室,杭州 310058)摘 要:为了实现番茄可溶性固形物含量(soluble
2、 solids content,SSC)的有效检测,提出高光谱漫透射成像检测方法,对比该成像方式下不同姿态(果脐端面姿态 BS、赤道圆周 3 姿态 C1、C2、C3 以及组合姿态 C1C2C3)的检测效果。首先对采集的不同姿态光谱图像,通过剪裁消除图像边缘噪声。针对圆周赤道面姿态 C1、C2 和 C3,进行了拼接处理,获得组合姿态图像 C1C2C3。其后对以上 5 种姿态图像进行单波段背景分割,获取目标区域,并统计不同姿态下番茄漫透射平均光谱。最后利用漫透射光谱结合偏最小二乘回归(partial least squares,PLS)方法,对番茄 SSC 分别在 450720、720990、45
3、0990 nm 3 个波段进行定量分析。结果表明,组合姿态 C1C2C3在 3 个波段区域上整体检测效果优于单个姿态的检测效果,其模型验证集均方根误差(root mean squared error of prediction,RMSEP)分别为 0.299%、0.133%、0.151%;相关系数 rp分别为 0.42,0.89,0.90。说明利用高光谱漫透射成像,获取组合姿态光谱图像,可以有效检测番茄 SSC。关键词:图像技术,光谱分析,果实,高光谱漫透射成像,成像姿态,番茄,可溶性固形物 doi:10.3969/j.issn.1002-6819.2013.23.034 中图分类号:TP75
4、1;S641.2 文献标志码:A 文章编号:1002-6819(2013)-23-0247-07 张若宇,饶秀勤,高迎旺,等.基于高光谱漫透射成像整体检测番茄可溶性固形物含量J.农业工程学报,2013,29(23):247252.Zhang Ruoyu,Rao Xiuqin,Gao Yingwang,et al.Inspection of soluble solid content for tomatoes in different positions based on hyperspectral diffuse transmittance imagingJ.Transactions of t
5、he Chinese Society of Agricultural Engineering(Transactions of the CSAE),2013,29(23):247252.(in Chinese with English abstract)0 引 言 番茄(Lycopersicum esculentum Mill.)因其营养丰富,又具药用价值1,深受人们喜爱,已成为世界第二大蔬菜消费品2。可溶性固形物含量是番茄制品品质评价的重要指标之一3,决定番茄制品的等级和价格。传统的光电和视觉技术能够实现番茄的外部颜色和缺陷识别,但难以有效检测可溶性固形物等内部品质4-5。近红外光谱分析技术作
6、为一种无损检测技术已成功用于多种果蔬内部品质的检测,包括苹果6、脐橙7、土豆8、西瓜9和番茄10-12等。并且上述研究表明,漫透射检测方式较反射和全透射检测方式更具优势。该检测方式不仅可以获取果蔬 收稿日期:2013-08-13 修订日期:2013-10-20 基金项目:国家自然基金项目(30825027),国家科技支撑项目(2011BAD20B12)作者简介:张若宇(1980),男,安徽临泉人,博士生,副教授,中国农业工程学会高级会员(E041200545S),主要从事农产品无损检测研究。杭州 浙江大学生物系统工程与食品科学学院,310058。Email: 通信作者:应义斌(1964),男,
7、浙江宁海人,教授,博士生导师,主要从事农产品/食品品质与安全快速检测技术和智能装备方面的研究。杭州 浙江大学生物系统工程与食品科学学院,310058。Email: 深层信息,还可以避免其形状大小及果核的影响。虽然近红外光谱能较好地表征果蔬内部品质信息,但是这种信息是一维的。由于果蔬形状、大小各异,内部成分和组织特征空间分布不同,而光纤探头检测的部位很小,所以光谱表达的信息不全面,预测模型缺乏稳健性13。高光谱成像技术融合了图像和光谱技术,可同时获取检测对象的空间和光谱信息,具有全波段连续光谱信息及识别目标对象不同化学组分的能力,已被逐渐引入果蔬品质检测领域14:如对土豆含水率15,香蕉总糖、含
8、水率、硬度16,雪花梨含糖量、含水率17等果蔬内部品质的检测。但前期研究发现,番茄表皮光亮,易形成表面亮斑,若反射成像,采集的光谱图像亮斑现象明显,严重影响后续番茄内部品质的定量光谱和图像分析。高光谱全透射成像可以获取果蔬深层次内部信息,但易受果核的影响,积分时间长,噪声大,此成像方法很少用于检测果蔬的内部组分,多用于检测果蔬的内部缺陷,如樱桃象鼻虫害18、马铃薯黑心病19、腌制黄瓜的内部缺陷等20。为了发挥高光谱透射方式成像优势,同时规避果核等因素影响,本文提出新的高光谱成像方式即农业工程学报 2013 年 248 高光谱漫透射成像。搭建了高光谱漫透射成像系统,并在此系统上对番茄可溶性固形物
9、进行检测,以验证该高光谱成像方式对番茄可溶性固形物含量检测的可行性。与此同时,通过对比不同姿态下番茄可溶性固形物含量的检测结果,进而证实番茄可溶性固形物空间分布差异,以寻求较为合理的检测姿态。1 材料与方法 1.1 试验样本 试验所用番茄购买自杭州蔬菜批发市场,产地为山东苍山。选取 110 个外观色泽基本一致的成熟番茄,用白色纱布将其表面擦拭干净,依次做好标记。番茄质量范围为 80164 g。1.2 高光谱漫透射成像系统 自行搭建的高光谱漫透射成像系统组成如图 1所示。硬件系统主要包括:行扫描成像光谱仪(ImSpector V10E-QE,Spectral Imaging Ltd.,Oulu,
10、Finland)、CCD 相机(C8484-05G,Hamamatsu Photonics,Japan)、50W 卤素灯 2 盏、自制激光传感器控制的限高板、电动位移平台(PSA200-11-X,Zolix Instruments Co.,Ltd.,北京)、电动升降平台(GZ02DS20,Guangzheng Instruments Co.,Ltd.,北京)。成像系统由 PC 机统一控制,由专业高光谱图像采集软件 SpectralCube(Spectral Imaging Ltd.,Finland)进行信息采集。成像光谱仪光谱范围为 4001 000 nm,曝光时间设定为 5 ms,水平位移速
11、度为 4.5 mm/s。采用 22 binning 压缩采集方式,得到拥有 477 波段,672120 像素大小的高光谱图。图 1 高光谱漫透射成像系统 Fig.1 Hyperspectral diffuse transmittance imaging system 1.3 光谱图像的采集 1.3.1 番茄光谱成像的姿态 获取不同姿态下番茄漫透射高光谱图像,并对比各姿态的检测效果。图 2a 所示将番茄果蒂端竖直向下放置,获取果脐端面图像,该姿态简称 BS;图 2b 所示将番茄果蒂端水平放置,获取番茄的圆周赤道面,包含 3 个姿态(C1、C2、C3),每 2个相邻位置圆周赤道面间隔 120。另外
12、在后续分析中,将单姿态 C1、C2、C3 进行组合形成组合姿态 C1C2C3。a.姿态 BS a.Position BS b.姿态 C1、C2 和 C3 b.Position C1,C2 and C3 图 2 番茄高光谱漫透射成像姿态 Fig.2 Positions of tomato in hyperspectral diffuse transmittance imaging 1.3.2 光谱图像黑白校正 由于各波段下光源强度分布不均匀及传感器中暗电流的存在,导致在光照强度分布较弱的波段下获得的图像含有较大噪声,因此必须进行标定21。本文采用高光谱图像采集软件 SpectralCube 对采
13、集的光谱图像进行黑白校正。首先获取聚四氟乙烯薄板的全白漫透射图像 IW,然后盖上镜头,获取全黑图像 ID,最后根据公式(1)进行黑白校正,获取校正后高光谱漫透射图像 IT。4095()()DTWDIIIII=(1)式中,IT为校正后高光谱漫透射图像,I为采集的原始图像,IW为全白图像,ID为全黑图像。4 095是系数22,因为高光谱系统数据输出为12位,数值(digital number,DN)范围为04095。1.4 可溶性固形物含量测定 应用手持式糖度计(PR-101,Atago Co.,Ltd.,Japan)测量番茄可溶性固形物含量。该糖度计的测量范围是045%,具有温度自动补偿功能。将
14、番茄不可食用部分切除后余下部分榨汁,滤纸过滤,滴入糖度计镜面并读取数值。最终样本数量为109个。番茄样本的可溶性固形物质量分数范围为3.20%5.50%,标准偏差为0.42%。2 结果与讨论 2.1 光谱图像的剪裁 为了消除图像边缘噪声,方便后续单波段背景第 23 期 张若宇等:基于高光谱漫透射成像整体检测番茄可溶性固形物含量 249 分割,同时为了减少数据运算量,将黑白校正后的光谱图像剪裁成200100像素大小。2.2 不同姿态光谱图像的拼接 为了考察组合姿态C1C2C3图像对可溶性固形物含量的检测效果,对C1、C2和C3姿态图像进行拼接(图3所示),然后整体进行背景分割,提取所有感兴趣区域
15、像素点,计算组合姿态图像的平均光谱。图 3 姿态 C1、C2 和 C3 图像的拼接 Fig.3 Images mosaicing of tomato in position C1,C2 and C3 2.3 单波段背景分割及感兴趣区域提取与统计 采用较少的波段进行背景分割,可以大大减小分割难度,提高分割效率。由番茄原始平均光谱可以看出,第250个波段(714.37 nm)是番茄漫透射光谱出现最大波峰处。并且发现番茄漫透射成像区域(regions of interest,ROIs)与背景在该波段的DN值差异较大便于分割。因此,选取第250个波段图像用来进行背景分割和感兴趣区域提取。番茄漫透射图像
16、采集时,限高板漏光孔边缘偶尔密封不严,引起非漫透射光进入相机,导致饱和区出现。为了避免饱和像素点被统计,设定DN值上限为4090,一旦某像素DN值大于4090意味着该像素逼近饱和,在计算平均光谱时将其忽略。因此DN值大于4090判别为饱和区,如图4a所示。另外背景在该波段DN值远小于500,最终选取500DN4090作为感兴趣区域提取的阈值范围。图4中白色区域为提取的目标区域,黑色区域为分割后的背景区域。a.单个姿态 a.Single position b.组合姿态 C1C2 C3 b.Integrated position C1C2 C3 图 4 背景分割和感兴趣区域提取 Fig.4 Bac
17、kground segmentation and extraction of regions of interest 统计感兴趣区域所有像素点,然后计算该区域在4001 000 nm内的平均光谱曲线,如图5a所示。2.4 光谱预处理 2.4.1 光谱归一化 本文采用的归一化方法是寻找每个样本透射率曲线的最大值,并计算各波段透射率与该最大值的比值,如式(2)所示。max()iNTTT=(2)式中,TN为归一化透射率,Ti为各波长透射率,max(T)为样本透射率曲线的最大值。这种方法可以得到各波长上透射率数据之间的相对关系,一定程度上可消除打光和采样环境等因素带来的影响。图5b为归一化后的番茄漫透
18、射光谱曲线。a.原始光谱 a.Original spectra b.归一化光谱 b.Normalized spectra 图 5 番茄漫透射光谱曲线(4001 000 nm)Fig.5 Diffuse transmittance spectra(4001 000 nm)of tomato 2.4.2 光谱校正 标准正态变量(standard normal variate,SNV)处理是光谱数据预处理常用方法之一。这种变换是将每个光谱减去平均值后再除以光谱的标准差,使得处理后的数据满足正态归一化特性,可以消除光谱中光程变化、光源功率变化和光强衰减等引起的噪声23。试验中番茄尺寸的差异、光程和测试
19、环境的轻微变化会引起光谱的噪声和基线的偏移,为消农业工程学报 2013 年 250 除这些因素的影响,将归一化后的光谱曲线再次采用 标 准 正 态 变 量 变 换 和 二 次 基 线 漂 移 去 除(quadratic removed baseline)进行校正处理。2.5 异常样本剔除 为了减少异常样本对检测效果的影响,建立校正模型前需剔除异常样本。同时为避免在判别过程中把非异常样本错误地当作异常样本剔除,本文采用Chauvenet检验方法和光谱杠杆值分布综合对异常样本进行判定并剔除24。首先将预处理后获得的样本光谱根据Chauvenet检验方法进行异常光谱判别。在95%的置信水平下,自动计
20、算出各样本光谱的马氏距离,超过临界线的样本视为异常样本,暂不做处理。与此同时,利用软件查看光谱杠杆值分布。一般当样本具有较高的光谱杠杆值时,视为异常样本。光谱杠杆阈值等于2nf/ns,(ns为样本数,nf为因子数)25。当2种方法都判定为异常样本时方做剔除处理。最后将余下的样本留作建模分析,其中75%样本作为校正集,25%样本作为模型验证集。2.6 可溶性固形物偏最小二乘回归模型建立 图6是不同姿态下获取的番茄漫透射光谱与其可溶性固形物含量在全波段(4001 000 nm)的相关系数分布图。本文选择450720、720990和 450990 nm 3个波段范围作为最重要的光谱区域进行偏最小二乘
21、回归(PLS)对比分析。图 6 不同姿态番茄漫透射光谱与其 可溶性固形物含量相关系数 Fig.6 Correlation coefficients between tomato spectra and soluble solids content over entire wavelength range 采用PLS分别对测得的可溶性固形物含量与预处理后不同姿态下番茄漫透射光谱进行回归分析。模型的准确性用校正集相关系数(rc)、校正集均方根误差(RMSEC)、验证集相关系数(rp)和验证集均方根误差(RMSEP)来描述。结果如表1所示。表 1 PLS 模型预测可溶性固形物结果 Table 1 P
22、erformance of PLS models for predicting soluble solids content 校正集 Calibration set 验证集 Validation set 姿态 Position 异常样本 Outlier samples 波段范围 Wavelength/nm 因子数 Number of factors均方根误差RMSEC/%相关系数rc 均方根误差 RMSEP/%相关系数rp 450720 4 0.340 0.60 0.379 0.24 720990 9 0.213 0.87 0.243 0.76 BS 2 450990 8 0.302 0.71
23、 0.342 0.44 450720 5 0.326 0.66 0.395 0.23 720990 11 0.125 0.96 0.197 0.88 C1 2 450990 17 0.114 0.96 0.224 0.82 450720 4 0.377 0.55 0.342 0.07 720990 14 0.056 0.99 0.267 0.57 C2 3 450990 4 0.367 0.58 0.330 0.13 450720 3 0.343 0.60 0.377 0.45 720990 8 0.226 0.85 0.360 0.47 C3 3 450990 17 0.119 0.96 0
24、.301 0.67 450720 9 0.314 0.72 0.299 0.42 720990 10 0.192 0.90 0.133 0.89 C1C2C3 1 450990 17 0.157 0.94 0.151 0.90 由表1可以看出,对于校正集,在450720波段区域,组合姿态C1C2C3最优,PLS模型的RMSEC=0.314%;rc=0.72。在720990 nm、450990 nm 2个波段区域,分别是单姿态C2和C1最优。而对于验证集,在450720、720990、450990 nm 3个波段区域上,组合姿态C1C2C3最优,RMSEP=0.299%、0.133%、0.151
25、%;rp=0.42、0.89、0.90。虽然单姿态C1、C2和C3在某个波段区域上,校正集建模效果优于组合姿态C1C2C3,但他们的验证集与校正集相关系数和均方根误差差别较大,说明其模型稳定性不好,可能是由于该姿态下软件自动划分的验证集没能很好覆盖校正集参数分布区域。从整体看组合姿态C1C2C3检测效果明显优于其他4个独立姿态的检测效果,这与图6所示番茄的漫透射光谱与SSC的相关性趋势基本吻合。对于上述姿态在波段范围450720 nm的检测第 23 期 张若宇等:基于高光谱漫透射成像整体检测番茄可溶性固形物含量 251 效果均比较差,但在波段区域720990 nm,除了C3外,其他姿态的检测效
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 光谱 透射 成像 整体 检测 番茄 可溶性 固形物 含量
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【二***】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【二***】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。