本文介绍了一种基于计算机视觉系统(CVS)的实验方案,用于测定多相食品体系的熔融行为。
本文介绍了一种基于计算机视觉系统(CVS)的实验方案,用于测定多相食品体系的熔融行为。
融化特性是冰淇淋最重要的质量指标之一。通常采用重量法进行评估,并以开始融化时间及融化速率来表示。然而,冰淇淋在融化过程中的外观同样重要,因为形状保持能力与产品结构的优良程度密切相关。本文提出的方案展示了一种计算机视觉系统(CVS),可用于补充现有的重力法,以计算两个与形状保持和融化速率相关的新融化指标。在总共90分钟的时间内,每15分钟拍摄一次冰淇淋融化的图像。随后,利用专门开发的图像处理方法对数字图像进行分析,以计算冰淇淋的面积、高度和宽度。在每次融化时间点测得的高度与宽度之比,称为 0 时刻的比值(Rt/ R0),是冰淇淋形状保持性的指标,而不同融化时间下的面积则相对于时间为0时的面积(At/A0与融化速率相关。该计算机视觉系统能够获得高灵敏度和高可靠性的结果,不仅可应用于冰淇淋,还可用于其他多种食品基质,如打发的牛奶奶油或蛋清。
冰淇淋是一种多相体系,其中液相、固相和气相紧密关联。连续的液相包裹着气泡和冰晶,其中含有部分结晶的脂肪、胶体蛋白、盐类、糖类(可能结晶)以及稳定剂。冰淇淋的成分因当地市场需求和可能的法规要求而异。尽管加工工艺会影响最终冰淇淋的特性,但每种成分在决定产品品质方面均起着重要作用1。融化行为是冰淇淋最重要的质量指标之一,需考虑在食用过程中以及口腔内发生的各种现象。当热量传入冰淇淋时,冰晶融化,水分扩散并与血清相混合,可能通过剩余结构渗出2。融化过快的产品不仅不利于舒适食用,也难以保证较高的抗热冲击能力;然而,融化过慢也可能表明配方存在某些缺陷1。已知冰淇淋的微观结构决定了其融化特性3,但迄今为止发表的研究结果存在矛盾,表明目前对微观结构因素如何影响融化过程的认识仍然有限4。因此,有必要开展更多研究以阐明融化机制,这对于新配方的设计也至关重要3。
融化行为通常通过重量法进行评估,并以开始融化时间及融化速率来表示5。将一定量的冰淇淋置于恒温箱内的金属丝网上,记录融化后产物的重量变化。根据重量-时间曲线,可区分出三个阶段:滞后期,即热量渗透阶段;快速融化期,此时稀薄的乳清相以最大速率流经冰淇淋结构;以及稳定期,此时大部分产品已滴落完毕2。
通过重量法可以区分慢融化和快融化产品;然而,冰淇淋在融化过程中的外观同样重要,因为形状保持能力与产品的良好组成和结构密切相关6。因此,基于计算机视觉系统(CVS)的方法可通过研究产品在融化过程中的外观,补充现有的重力测定方法。CVS能够以人眼无法观察到的精确细节获取多种食品属性3(例如尺寸、重量、形状、质地和颜色)。这类系统通常由数字相机和图像处理软件组成7。事实上,基于CVS的方案主要包括两个步骤:1)图像采集和2)图像处理。可采用不同层次的图像处理技术7,从最简单的处理方法到更复杂的深度学习方法,用于人工智能的开发8,9。近年来,食品领域对CVS的关注显著增加,已开发出大量应用于食品安全检测、食品加工监控、异物检测等方面的系统。这些系统快速、高效且具有非破坏性,因此是保障消费者获得高质量安全食品的有效工具10。
在冰淇淋领域,已有研究建议采用图像分析方法,通过光学显微镜观察冰的重结晶过程11。最近,研究人员还通过处理X射线计算机断层扫描图像,分析包括冰淇淋在内的软质多孔材料的三维微观结构3。然而,利用简单的数字电荷耦合器件(CCD)图像在获取便捷性以及呈现消费者所感知的冰淇淋外观方面具有多项优势。一些作者展示了冰淇淋在融化过程中的图像12,但据我们所知,首次从图像中提取数值指标的研究是由 Moriano 和 Alamprese 报道的13。
因此,本方案基于 Moriano 和 Alamprese 的研究工作提出13,展示了可用于支持现有重力法研究冰淇淋融化行为的简单CVS。所提出方法的框图如图所示 图1该系统的使用可计算出两个与形状保持和融化速率相关的熔融指数。特别是,本文首次详细描述了冰淇淋融化过程中进行数字图像采集的实验装置与操作步骤,以及图像处理流程。此外,还报告了使用不同甜味剂(即蔗糖、蔗糖麦芽糖和赤藓糖醇)制成的冰淇淋的实验结果,以展示该方法的潜力。

图1:所提出方法的框图。 所提出的用于研究冰淇淋融化行为的计算机视觉系统与重量分析方法的总体步骤概要。请点击此处查看此图的放大版本。
熔化实验的实验装置与操作步骤

图2:熔化测试装置示意图。 该图展示了如何在恒温箱中设置熔化试验:将一个带刻度的玻璃量筒放在电子天平上,用于收集并称量融化的冰淇淋。冰淇淋样品置于悬于量筒上方漏斗中的金属丝网筛上。 请点击此处查看该图的放大版本。
2. 用于计算熔融指数的图像处理
(1)
(2)
(3)补充图1:图像空间校准。(A)在图像分析软件中进入测量 > 校准 > 空间窗口。选择新建,然后勾选图像以打开比例窗口。此处标明了用于将像素转换为实际单位(例如毫米)的参考长度。(B)将绿色标尺条与对应于该指定长度的参考部分精确重叠,然后点击确定。请点击此处下载该文件。
补充图2:裁剪感兴趣区域(AOI)并将其转换为灰度图像。(A)将感兴趣区域(AOI)转换为灰度图像,以及(B)得到的图像。请点击此处下载该文件。
补充图3:待测参数的选择。 在选择测量窗口中,可选择需要测量的参数;对于冰淇淋融化评估,必须选择面积、框宽度和框高度。请点击此处下载该文件。
补充图4:冰淇淋样品的分割。 在“分割”窗口中,可选择直方图范围,以精确覆盖冰淇淋形状的区域。请点击此处下载该文件。
补充图5:对象筛选与计数功能。 红色线条标出了被识别出的明亮物体。通过应用计数功能并打开“视图,测量数据”窗口,所选参数的结果将被显示(A)。为了仅筛选出冰淇淋形状的物体,可在“选择测量”窗口中设定最小和最大面积范围,从而仅对单一物体的参数进行计数(B)。请点击此处下载该文件。
补充图6:形状保持指数(R)。 用于计算形状保持指数(R)的盒体高度(Ht,红色虚线)和盒体宽度(Wt,黑色实线)如图所示。 请点击此处下载该文件。

图 3:形状和面积保留曲线。 冰淇淋的(A)形状保留曲线和(B)面积保留曲线示例,其中随时间绘制了 Rt/R0 和 At/A0 的平均值;误差条表示由重复分析得到的标准偏差值。 请点击此处查看该图的放大版本。
3. 重力数据的阐述
(4)
(5)
图 4:重量法曲线。 通过重量法测得的冰淇淋融化曲线示例。原始曲线以红色显示;线性部分所选数据序列以绿色显示;计算得到的回归线以黑色圆点表示。回归线的方程及决定系数(R2)也一并标出。 请点击此处查看该图的放大版本。
注意:为了获得可用于统计分析的可靠结果,每种样品的熔点测定实验及图像处理全过程至少需重复三次。
作为所提出的计算机视觉系统(CVS)输出的示例,展示了三种不同冰淇淋配方的融化分析结果,并与重量法获得的数据进行了比较。特别是研究了使用不同甜味剂(即蔗糖、蔗糖麦芽糖和赤藓糖醇)制成的冰淇淋的融化行为。
表1 和 图5A 展示了三种冰淇淋样品在融化过程中形状保持指数(Rt/R0)的结果。在前45分钟内,所有样品表现出相似的形状保持能力;此后,蔗糖样品的形状丧失速度较其他产品更快,根据单因素方差分析(ANOVA)及随后的最小显著差异(LSD)检验(p < 0.05),其 Rt/R0 值显著更低。在整个试验期间,仅有含赤藓糖醇的样品保持了完整形状,因而获得了可测量的 Rt/R0 值。
| 蔗糖 | 蔗糖麦芽糖 | 赤藓糖醇 | ||||
| 时间(分钟) | At/A0 | 标准差 | At/A0 | 标准差 | At/A0 | 标准差 |
| 0 | 1.00 | - | 1.00 | - | 1.00 | - |
| 15 | 0.91a | 0.01 | 0.93a | 0.02 | 0.92a | 0.04 |
| 30 | 0.76a | 0.04 | 0.82a | 0.02 | 0.84a | 0.04 |
| 45 | 0.56a | 0.04 | 0.66ab | 0.02 | 0.74b | 0.03 |
| 60 | 0.32a | 0.01 | 0.46b | 0.02 | 0.62c | 0.01 |
| 75 | 0.13a | 0.01 | 0.24b | 0.02 | 0.52c | 0.02 |
| 90 | n.d. | - | n.d. | - | 0.40 | 0.01 |
| n.d.,不可检出 | ||||||
表1:使用不同甜味剂制备的三种冰淇淋样品的形状保持指数的平均值和标准差。 同一行中不同的小写字母表示平均值之间存在显著差异(p < 0.05)。

图 5:使用不同甜味剂制备的三种冰淇淋样品的形状保持率和面积保持率曲线。(A)形状保持率曲线和(B)面积保持率曲线。粉线代表作为参照的蔗糖样品;黄线代表蔗糖麦芽糖样样品,蓝线代表赤藓糖醇样品。请点击此处查看该图的放大版本。
添加赤藓糖醇的样品也表现出更低的熔化速率,这一点由A结果可见t/A0 熔融指数(表2 和 图5B)。从45分钟开始,赤藓醇样品表现出显著差异(p < 0.05)最高的面积保留率,其次是含有蔗麦芽糖和蔗糖的冰淇淋。
| 蔗糖 | 蔗糖麦芽糖 | 赤藓糖醇 | ||||
| 时间(分钟) | Rt/R0 | 标准差 | Rt/R0 | 标准差 | Rt/R0 | 标准差 |
| 0 | 1.00 | - | 1.00 | - | 1.00 | - |
| 15 | 0.98a | 0.04 | 0.97a | 0.01 | 0.93a | 0.02 |
| 30 | 0.94a | 0.05 | 0.94a | 0.02 | 0.89a | 0.02 |
| 45 | 0.80a | 0.03 | 0.89a | 0.02 | 0.84a | 0.02 |
| 60 | 0.64a | 0.04 | 0.78b | 0.01 | 0.76b | 0.03 |
| 75 | 0.54a | 0.03 | 0.65a | 0.01 | 0.60a | 0.07 |
| 90 | n.d. | - | n.d. | - | 0.45 | 0.06 |
| n.d.,不可检出 | ||||||
表2:使用不同甜味剂制成的三种冰淇淋样品的面积保留指数的平均值和标准差。 同一行中不同的小写字母表示平均值之间存在显著差异(p < 0.05)。
结果还通过重量分析法进行了验证,相关数据以熔化开始时间(tS)和熔化速率(克/分钟)表示,见表3。然而,在此情况下,应用方差分析(ANOVA)未发现显著差异(p > 0.05)。
| 蔗糖 | 低聚麦芽糖醇 | 赤藓糖醇 | ||||
| 均值 | 标准差 | 均值 | 标准差 | 均值 | 标准差 | |
| ts (min) | 15.3a | 1.4 | 16.4a | 1.4 | 17.2a | 2.6 |
| 融化速率 (g/min) | 2.67a | 0.26 | 2.29a | 0.09 | 1.96a | 0.41 |
表3:使用不同甜味剂制备的三种冰淇淋样品通过重量法测定的融化指数的平均值和标准差。 ts,融化起始时间。同行中不同的小写字母表示平均值之间存在显著差异(p < 0.05)。
图6展示了在三种冰淇淋样品的融化试验(步骤1.3)过程中获得的部分已处理图像(步骤2.1);可以看出,所有样品均具有良好的形状保持性,而赤藓糖醇制成的冰淇淋融化速率最低,在90分钟后仍未完全融化。

图6:三种冰淇淋样品的融化过程。 在使用不同甜味剂制备的三种冰淇淋样品融化过程中,每15分钟拍摄一次图像。三种冰淇淋样品在融化过程中的处理后图像。请点击此处查看该图的放大版本。
该提议的计算机视觉系统(CVS)可用于计算冰淇淋样品在融化过程中的形态和面积保持指数,同时可视化其融化过程。该系统可与传统的重量法相结合,用于评估冰淇淋的融化特性5,从而获得与冰淇淋外观相关的结果。这一点非常重要,因为消费者对产品质量的评价也基于产品的视觉外观,而融化过程中保持形状的能力与冰淇淋的成分和结构密切相关1,3。该系统在食品工业中的应用可提升质量控制水平,并促进新配方的研发。
使用不同甜味剂的三种冰淇淋样品的融化实验结果表明,所提出的CVS方法对不同配方具有敏感性(表1 和 表2),其敏感性甚至高于重量法,后者未能凸显样品间的显著差异(表3)。我们先前的研究也已证明形状保持指数对冰淇淋配方的敏感性13,14,15,16。由于该方法具有高敏感性,甚至可替代重量法测定程序。
除了冰淇淋样品的制备外,所建立的CVS方法中最关键的步骤与高质量的图像采集和分割密切相关7。将Rt和At数据相对于时间0时计算的初始值进行归一化处理,可避免因不同重复实验之间初始设置差异而导致的误导性结果。然而,在图像采集过程中可能出现视差误差,尤其是在单次试验中,若相机在不同采集时间点发生位置偏移,会导致所采集图像不一致,从而在融化过程中高估样品面积。因此,必须高度重视图像采集程序的设置,确保放置样品的金属丝网筛与相机镜头完全对齐,并精确定义和固定相机在恒温箱前方的位置。同样,图像分割效果不佳可能导致对冰淇淋面积、以及高度(H)和宽度(W)的高估或低估,进而导致形状保持指数失去意义。在此情况下,主要误差来源可能是位于丝网上冰淇淋样品边缘可能出现的尖状突起,或在融化过程中冰淇淋表面产生的气泡(图7),这些干扰物必须在测量包围盒面积及盒的W和H之前从亮区轮廓中去除。可通过使用图像分析软件中的分割对象(Split Object)功能来解决此问题。此外,高对比度的背景有助于实现更佳的图像分割效果7。背景颜色应根据冰淇淋样品的颜色进行选择,且较高的色彩饱和度有利于获得更优的对比度。

图7:冰淇淋形状分割中可能的误差来源。 红色箭头指示冰淇淋形状上可能出现的尖刺和气泡,这些会导致面积、高度和宽度的错误计算。 请点击此处查看该图的放大版本。
如果图像采集和处理得当,计算机视觉法(CVS)的灵敏度和精确度甚至优于重量法。事实上,在所示的示例中,重量法分析的变异系数介于4%至20%之间,导致所有平均数据之间无显著差异(ANOVA,p > 0.05)。相反,CVS指标的变异系数范围为0.1%至13%,揭示了三种冰淇淋样品在融化过程中存在显著差异(表1 和 表2)。赤藓糖醇导致融化速率降低,可能与冰淇淋冷冻过程中非晶态糖晶体的形成有关,且该甜味剂相较于蔗糖和蔗糖麦芽糖,溶解度较低,从而使最终产品中冻结水的比例更高13。本研究中所提出的CVS方法的其他成功应用,可参见我们此前的研究13,14,15,16。
当图像处理流程实现良好标准化后,可利用图像分析软件的宏记录功能设置自动图像处理,从而显著加快数据处理速度。当然,操作人员仍需对数据进行仔细审核,以避免伪影和错误数据的产生。此外,还可通过分析形状和面积保留曲线来计算动力学参数,以替代传统的重量法。
该方案也可应用于其他对融化速率和形状保持指数有要求的食品,例如慕斯和泡沫类食品(如打发的牛奶奶油和蛋清)。最重要的要点是产品需具备自支撑性,而融化试验过程中的温度则可根据具体产品进行调整。
作者无任何利益冲突需要披露。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 柜体 | Cavallo s.r.l. | FTX700 | 熔点测试放置位置 |
| 数码相机 | Sony Group Corp | DSC-S650 | |
| 电子天平 | Gibertini Elettronica | EU-C 4002 LCD | |
| ImagePro Plus 7.0 | Media Cybernetics, Inc | N/A | 图像分析处理软件 |
| Microsoft Excel | Microsoft | N/A | 数据与图形处理 |
| Scalecom | Gibertini Elettronica | N/A | 电子天平数据采集软件 |
| 三脚架 | Manfratto | #055 |
申请许可以重复使用本 JoVE 文章的文本或图表
申请许可