可见/近红外高光谱成像技术的西甜瓜糖度检测研究
发布时间:2023-04-07
浏览次数:487
目前大多数研究者利用光谱仪检测果类品质、糖度、酸度等时,仅对单一品种或者某一类水果进行研究,所构建的模型有一定的极限性。针对这种情况,本研究提出了以不同品种的西瓜、甜瓜为研究对象,利用可见光/近红外成像高光谱技术分析不同糖度的西瓜、甜瓜的光谱差异及西瓜、甜瓜糖度在可见光/近红外波段范围的响应。
目前大多数研究者利用高光谱成像仪检测果类品质、糖度、酸度等时,仅对单一品种或者某一类水果进行研究,所构建的模型有一定的极限性。针对这种情况,本研究提出了以不同品种的西瓜、甜瓜为研究对象,利用可见光/近红外成像高光谱技术分析不同糖度的西瓜、甜瓜的光谱差异及西瓜、甜瓜糖度在可见光/近红外波段范围的响应。实验结果表明,不同品种不同糖度的西瓜,其光谱曲线变化趋势相似,在650 -800 nm处,随着糖度升高其光谱反射率呈增加的趋势,在570 - 650 nm 区间反射率呈直线上升;不同品种不同糖度的甜瓜,其光谱曲线变化趋势相似,随着糖度的升高反射率呈上升趋势,在550 -580 nm区间反射率上升较快,形成陡坡;不同品种的西瓜其糖度在可见光/近红外波段相关性较高的波段主要集中在650 - 700 nm之间,不同品种的甜瓜其糖度在可见光/近红外波段相关性较高的波段则主要集中在670 - 690 nm之间,而综合西瓜、甜瓜两种水果,其糖度在可见光/近红外波段相关性较高的集中在630 - 650 nm,其中相关系数最高的是639.3 nm处,其相关系数为0.951;利用独立的建模与检验数据,运用639.3 nm处的反射率构建监测模型,并对该模型进行检验,其建模R2为0.904,检验R2为0.847,斜率Slope为0.905,截距为0.867,相对均方根误差RRMSE为6.78%。研究表明应用成像高光谱技术检测西甜瓜糖度具有可行性,为进一步研究不同水果糖度高精度模型奠定基础。
0 引言
西瓜、甜瓜含糖量较高,营养丰富,深受消费者的喜爱。糖度是影响西甜瓜营养与风味的主要因素,既可表征成熟度,又可体现其品质[1]。传统的品质化学检测方法是一种费时费力的破坏性检测技术,随着成像和光谱技术的快速发展,高光谱成像技术(Hyperspectral imaging)已经广泛应用于农产品品质的快速无损检测[2]。高光谱成像技术融合了传统的成像和光谱技术的优点,可以同时获取被检测物体的空间信息和光谱信息,因此该技术既可以检测物体的外部品质,也可以检测物体的内部品质如水果的糖度、坚实度、酸度等[3-10]。
董一威等[11]以红富士苹果为研究对象,应用自行搭建的CCD近红外光谱系统检测苹果的糖酸度,研究结果表明在630 nm-1030 nm范围内,可实现对苹果糖度、酸度的无损检测研究;金同铭等[12]利用近红外光谱法检测苹果中蔗糖、葡萄糖、果糖三种组分,筛选出914 nm、950 nm、897 nm分别代表蔗糖、葡萄糖和果糖的第一特征波长;赵丽丽等[13]运用近红外光谱仪获取苹果和番茄的光谱,并建立苹果糖度、硬度及番茄的红素的数学模型,所有模型的相关系数均在0.8以上,证明利用红外区域光谱可进行果品类品质的检测;田海清等[14]利用透射光谱仪获取蜜瓜的光谱信息,运用经典最小二乘法、逐步多元线性回归、主成分分析、偏最小二乘法等方法构建模型并检验,结果表明PLS法的建模与预测结果较好;Polder 等[15]利用可见光(波长范围为396-736 nm)高光谱成像系统对西红柿的成熟度进行判别研究;Lu[16]利用高光谱成像系统研究了苹果的坚实度;Masry 等[17]采用波长范围为400-1000 nm 的近红外高光谱成像技术检测了草莓内部含水量和可溶性固形物含量。赵杰文等[18]采集了苹果的高光谱图像(685-900nm),研究发现利用PLS回归模型预测糖度的优波长为704.48-805.26nm,预测的相关系数为0.907;马本学等[1]利用漫反射高光谱成像技术研究了网纹类甜瓜糖度的无损检测方法。
然而,大多数研究者利用光谱仪检测果类品质、糖度、酸度等时,仅针对单一品种或者某一类水果进行研究,很少有研究者利用多类水果或多品种水果进行研究分析其不同糖度、酸度的光谱差异及特征响应波段范围。因此本研究以无籽/有籽西瓜、青色/黄色甜瓜为研究对象,利用可见/近红外高光谱成像技术,探索不同种类西甜瓜不同糖度光谱差异及特征响应波段范围,为探索不同水果糖度高精度检测模型奠定基础。
1 材料与方法
1.1 试验材料
实验材料2015年9月取自大兴庞各庄实验基地。包括不同类型西甜瓜品种(有籽西瓜、无籽西瓜);甜瓜(青皮甜瓜、黄皮甜瓜)。进行瓜瓤的光谱采集和含糖量测定。
1.2 高光谱图像采集
高光谱图像数据采集采用广东赛斯拜克厂家的高光谱成像仪。该系统主要由CCD相机、卤素灯光源、暗箱、计算机组成,如图1。实验仪器参数设置如表1。
图1 高光谱成像仪
在进行高光谱图像采集时,需要设置相机曝光时间,平台移动速度以及物镜之间的距离。这 3 个参数相互影响,图像调节的目的是使采集的图像大小合适,清晰,不变形失真。经过反复尝试,物镜高度设置为 26 cm,曝光时间设置为10ms,平台移动速度设置为 4.3 mm/s。图像采集软件采用广东赛斯拜克厂家提供的高光谱成像系统采集软件完成。图像处理采用软件进行处理。在进行图像处理之前,先要对采集的光谱图像进行图像校正,图像校正公式如下:
(1)
式中,Rref 是校正过的图像,DNraw 是原始图像,DNwhite为白板校正图像,DNdark 是黑板校正图像。
1.3 糖度测定与光谱反射率的获取
西瓜、甜瓜样本光谱数据采集后,为避免光照对糖度造成一定的变化,因此在微光条件下,在西瓜、甜瓜样品上分别作五处标志,并立即用PR-101数字式折射仪(AtagoCo, Ltd. , Tokyo, Japan)对标志部位进行糖度测量,分切取约1 5 mm×10 mm×5 mm中部内缘果肉挤出汁液滴到折光仪镜面上,覆盖住底部,测量其糖度并记录。
1.4 光谱获取
试验得到光谱含有由仪器和试验条件等引起的噪声,对这些噪声的处理有助于减少噪声对光谱分析的影响,突出光谱的有效信息。Savitzky-Golay (SG)平滑算法可以有效消减光谱数据中的随机噪声,消噪效果受平滑点数的影响,本文中选择SG二次多项式7点平滑对光谱数据进行处理[19]。在高光谱图像上找到西瓜、甜瓜的五个标志处,运用ENVI的感兴趣区域提取功能,获取糖度测量区域,提取该区域的平均光谱,将此平均光谱与其对应的糖度建立一一对应的关系。
1.5 建模与检验数据
本研究一共获取了西甜瓜共36处的糖度值并获取其对应区域的平均光谱反射率,其中西瓜共有20个糖度值,甜瓜共有16个糖度值,采用随机抽样的方式获取11个西瓜糖度值、9个甜瓜糖度值及其光谱反射率作为建模数据,其余西瓜、甜瓜糖度值及对应的光谱反射率作为检验数据。西甜瓜建模与检验集的糖度值统计信息如表2所示,本研究采用建模R2、检验R2、相对均方根误差(RRMSE)、Slope和截距等来评价模型性能[20] 。
表2 西甜瓜建模与检验集的糖度值统计信息
数据集 | 西甜瓜类别 | 个数 | 最小值 (%) | 最大值 (%) | 均值 (%) | 标准差 (%) |
建模集 | 西瓜 | 11 | 6.40 | 11.80 | 9.31 | 2.08 |
甜瓜 | 9 | 8.00 | 14.60 | 10.73 | 2.21 | |
西甜瓜 | 20 | 6.40 | 14.60 | 10.06 | 2.16 | |
检验集 | 西瓜 | 9 | 6.60 | 10.20 | 8.96 | 2.14 |
甜瓜 | 7 | 8.60 | 11.82 | 10.04 | 1.88 | |
西甜瓜 | 16 | 6.60 | 11.82 | 9.38 | 1.95 | |
全部 | 西瓜 | 20 | 6.40 | 11.80 | 9.15 | 2.11 |
甜瓜 | 16 | 8.00 | 14.60 | 10.42 | 2.09 | |
西甜瓜 | 36 | 6.40 | 14.60 | 9.76 | 2.10 |
2 结果与分析
2.1 有籽西瓜不同含糖量的光谱分析
图2为有籽西瓜不同含糖量的光谱曲线图,从图中可知,不同含糖量在红光至近红外区间差异明显,随着糖度升高其光谱反射率呈增加的趋势;在550 nm、700 nm、800 nm附近出现小峰值,这由于有籽西瓜的有机分子中含氢基团振动的合频、各级倍频的吸收作用引起的;在570-650 nm 红光区域反射率呈直线上升,这与西瓜含有红色素有关。
图2 有籽西瓜不同含糖量的光谱曲线差异
2.2 无籽西瓜不同含糖量的光谱分析
图3为无籽西瓜不同含糖量的光谱曲线图,从图中可知,在400 -600 nm范围内,高糖度的光谱反射率与低糖度的光谱反射率差异显著,高糖度的光谱反射率大于低糖度的光谱反射率;在650-900 nm处,随着糖度升高其光谱反射率升高;在570-650 nm 红光区间反射率呈直线上升,这与西瓜含有红色素有关;在550 nm、700 nm、800 nm附近出现小峰值。
综合图2和图3两种西瓜品种的光谱反射率曲线可知,西瓜瓜瓤在550 nm、700 nm、800 nm附近有峰值;在570-650 nm 红光区域反射率呈直线上升,因此570-650 nm波段范围与550 nm、700 nm、800 nm三个波段都可以作为西瓜瓜瓤的特征波段。
图3 无籽西瓜不同含糖量的光谱曲线差异
2.3 黄色甜瓜不同含糖量的光谱分析
图4为黄色甜瓜不同含糖量的光谱曲线图,从图中可知,随着含糖量的升高在可见光、近红外区间其光谱反射率呈增加的趋势。在410-600 nm 区间反射率上升速度较快;在600-680 nm区间反射率上升缓慢;在680-720 nm区间缓慢下降;在720-800 nm处呈上升趋势;到了800 nm后,光谱反射率呈下降趋势;与西瓜瓜瓤相似,在700 nm和800 nm附近有2个峰值,与西瓜一样,这也是由于甜瓜的有机分子中含氢基团振动的合频、各级倍频的吸收作用引起的。总体而言,黄色甜瓜糖度不同,其光谱曲线差异显著。
图4 黄色甜瓜不同含糖量的光谱曲线差异
2.4 青色甜瓜不同含糖量的光谱分析
图5为青色甜瓜不同含糖量的光谱曲线图,从图中可知,在430-550 nm范围内,随着含糖量的升高其光谱反射率减小;在550-700 nm区间随着甜瓜糖度的升高,其反射率升高;在550-580 nm区间反射率上升较快,形成陡坡;在680 nm附近,低糖度的甜瓜瓜瓤出现低谷,含糖量越高,谷值越深;与西瓜瓜瓤相似,在700 nm和800 nm附近也有2个峰值;在700-1000 nm区间,糖度与光谱反射率曲线无显著规律变化。
综合图4和图5可知,在400-550 nm范围内,黄色甜瓜的光谱反射率显著高于青色甜瓜,甜瓜在700 nm、800 nm附近有峰值;在500-600 nm 范围内反射率呈直线上升,形成陡坡,因此500-600 nm波段范围与700 nm、800 nm两个波段都可以作为甜瓜瓜瓤的特征波段。
图5青色甜瓜不同含糖量的光谱曲线差异
2.5 糖度相近时不同品种不同种类水果的光谱曲线分析
图6为糖度相近时不同品种不同种类水果的光谱曲线图。从图中可知,西瓜虽然品种不同,其光谱曲线相近,在580-650 nm范围内光谱反射率直线上升,形成陡坡,在550 nm、700 nm、800 nm三个波长附近有三个吸收峰。两种甜瓜在400-600 nm范围内,光谱曲线变化规律相似,在520-580 nm附近光谱反射率曲线直线上升,形成陡坡;但在600-720 nm范围内,两个品种的甜瓜光谱曲线差异显著,青色甜瓜在680 nm附近有吸收谷,而黄瓜甜瓜没有出现吸收谷,在720-1000 nm范围内,两个品种的甜瓜光谱曲线变化趋势相似。比较西瓜、甜瓜两种水果的瓜瓤光谱曲线可知,两种水果都陡坡形成,只是陡坡形成的位置有所差异,在700 nm和800 nm附近均有吸收峰,在720-1000 nm范围内,不同品种的西瓜、甜瓜光谱反射率曲线变化趋势相似。
图6 糖度相近不同西甜瓜的光谱曲线图
2.6 水果糖度与各波段的光谱反射率的相关性分析及模型构建
图7为不同水果的糖度与各波段反射率的相关系数图。从图中可知,四种不同品种的水果,其光谱反射率对糖度的响应差异较大。对两个品种的西瓜而言,无籽西瓜在400-600 nm范围内,其光谱反射率与含糖量呈负相关关系,在600-1000 nm范围内,其光谱反射率与含糖量呈正相关关系,其中600 nm处为零界点;有籽西瓜在400-630 nm范围内,其光谱反射率与含糖量呈负相关关系,在620-700 nm与820-1000 nm范围内,其光谱反射率与含糖量呈正相关关系,在700-820 nm范围内,其光谱反射率与含糖量相关性较小,接近为0;在660 nm附近,有籽和无籽西瓜的光谱反射率与含糖量形成一个峰值,因此,660 nm可作为西瓜瓜瓤对糖度的一个特征响应波段。而对于两个品种的甜瓜,黄色甜瓜的光谱反射率与糖度在400-1000nm范围内呈正相关关系,但在520 nm附近出现一个谷值,该谷值为0,即与糖度相关性为0;在550-1000 nm范围内,黄色甜瓜的光谱反射率与糖度相关性较为稳定,保持在0.8左右;青色甜瓜在400-570 nm范围内,其光谱反射率与含糖量呈负相关关系,在570-1000 nm范围内,其光谱反射率与含糖量呈正相关关系,与黄色甜瓜类似,在520 nm附近出现一个谷值,该谷值为0.8,即相关系数为0.8;青色甜瓜的光谱反射率与含糖量的零界点在570 nm,即相关系数为0。
图 7 西甜瓜的糖度与各波段反射率的相关系数
图8为综合两种西瓜、两种甜瓜的糖度与各波段反射率的相关系数图,从图中可知,西瓜、甜瓜的糖度与各波段的光谱反射率相关性较高的集中在630-650 nm,其中相关系数最高的是639.3 nm处,其相关系数为0.951。
图 8 综合西甜瓜的糖度与其波段反射率的相关系数
图9为两种西瓜、两种甜瓜在639.3 nm处的光谱反射率与其糖度的散点分布图。根据散点图的分布进行线性拟合,其一次线性拟合的模型为:y=30.941x-11.05,决定系数R2为0.904。
图 9 西瓜、甜瓜在639.3 nm处光谱反射率与其糖度的散点图
运用独立的数据对该线性模型y=30.941x-11.05进行检验,如图10所示,其实测值与预测值之间的决定系数R2为0.847,斜率Slope为0.905,截距为0.867,相对均方根误差为6.78%。结果表明由639.3 nm处的光谱反射率与其糖度构建的一次线性模型对于西瓜、甜瓜等糖度的预测具有较好的效果。
图10 639.3 nm处的光谱反射率与糖度的模型检验
3 结论
本研究通过分析不同品种的西瓜、甜瓜的光谱曲线可知,对于不同品种的西瓜而言,其光谱曲线总体变化趋势相似,在550 nm、700 nm、800 nm附近有吸收峰,在570-650 nm 红光区域反射率呈直线上升,形成陡坡,因此570-650 nm波段范围与550 nm、700 nm、800 nm三个波段都可以作为西瓜瓜瓤的特征波段。对于不同品种的甜瓜,在400-550 nm范围内,黄色甜瓜的光谱反射率显著高于青色甜瓜,甜瓜在700 nm、800 nm附近也有吸收峰;在500-600 nm 范围内反射率呈直线上升,形成陡坡,因此500-600 nm波段范围与700 nm、800 nm两个波段都可以作为甜瓜的特征波段。比较相似糖度下的西瓜和甜瓜的光谱曲线可知,两种水果都陡坡形成,只是陡坡形成的位置有所差异,甜瓜的光谱曲线陡坡在500-600 nm附近,西瓜的光谱曲线陡坡在570-650 nm附近;在700 nm和800 nm附近均有吸收峰,在720-1000 nm范围内,不同品种的西瓜、甜瓜光谱反射率曲线变化趋势相似。
分析不同品种的西瓜、甜瓜的光谱反射率与糖度的相关系数可知,不同品种的西瓜、甜瓜,其光谱反射率对糖度的响应差异较大,但在红光范围630-670范围内均有较高的相关性。综合不同品种、不同种类的西甜瓜光谱反射率,分析其对糖度的响应波段范围,研究表明,波段响应最高的是639.3 nm,其相关系数为0.951,利用独立的建模与检验数据,运用639.3 nm处的光谱构建监测模型,并对该模型进行检验,其建模R2为0.904,检验R2为0.847,斜率Slope为0.905,截距为0.867,相对均方根误差为6.78%,取得较好的研究效果,研究结果表明应用成像高光谱技术检测西甜瓜糖度具有可行性,为进一步研究不同水果糖度高精度模型奠定理论基础。
相关产品
-
干旱区基于高光谱的棉花遥感估产研究
棉花产量与冠层光谱植被指数相关关系,建立棉花高光谱估算模型,可以促进高光谱技术在棉花长势监测和估产中应用。本文简单总结了干旱区基于高光谱的棉花遥感估产研究。..
-
近红外光谱分析技术在精准农业中的应用
近红外光谱分析技术在精准农业中的应用非常广泛且深入,主要体现在以下几个方面:..
-
近红外高光谱相机的工作原理
近红外高光谱相机以其高分辨率、高光谱分辨率以及广泛的应用领域,在科研、工业检测、环境监测等多个领域发挥着重要作用。本文简单总结了近红外高光谱相机的工作原理。..
-
高光谱图像处理的主要方法有哪些?
对高光谱图像进行处理,可以获取物体在不同波段的反射光谱信息,提供更加详细和全面的图像数据。那么,什么是高光谱图像处理?高光谱图像处理的主要方法有哪些?下文对此做..