综述:数据科学时代下的食品微观结构:从图像到价值

时间:2026年5月18日
来源:Food Structure

编辑推荐:

里娅·G·罗德里格斯(Ria G. Rodriguez)| 伊娜·弗尔默(Ina Folmer)| 詹斯·萨尔布林克(Jens Saalbrink)| 米克·施朗根(Miek Schlangen)| 黛安娜·索托-阿吉拉尔(Diana Soto-Aguilar)| 马蒂亚斯·P·

广告
   X   

里娅·G·罗德里格斯(Ria G. Rodriguez)| 伊娜·弗尔默(Ina Folmer)| 詹斯·萨尔布林克(Jens Saalbrink)| 米克·施朗根(Miek Schlangen)| 黛安娜·索托-阿吉拉尔(Diana Soto-Aguilar)| 马蒂亚斯·P·克劳森(Mathias P. Clausen)| 何塞·C·博尼利亚(José C. Bonilla)
哥本哈根大学食品科学系

**摘要**
量化食品中的微观结构是理解食品材料在不同长度尺度上如何分层组织的关键。本文综述了在食品微观结构定量分析方面所做的努力,并旨在推动该领域向定量研究方向发展。目前,对于食品微观结构的分类仍缺乏共识,这可能导致多种解释。此外,食品的微观结构多种多样,并且通常不像其他学科那样清晰易懂。数据科学方法(包括机器学习和深度学习)并不新鲜,但它们刚刚开始被用于食品微观结构的分割和量化。总体而言,向更定量科学的转变将推动食品科学朝着合理设计食品的方向发展,包括针对食品质地和其他结构相关属性的定制设计。

**1. 引言**
食品结构是分层组织的,因此需要从分子组成到宏观特性等多个长度尺度进行研究。目前,对食品进行微观尺度研究的尝试往往缺乏定量分析。能够观察到微观尺度使得可以研究食品大分子的物理排列及其在基质中的物理化学相互作用。此外,研究表明食品的微观结构会影响其机械性能和感官特性,并且与产品的稳定性和营养素生物利用率有关(Hiolle等人,2020年)。将食品材料的微观特性与其宏观物理特性联系起来的重要性已被广泛认可,自20世纪80年代以来就开始研究食品的结构-功能关系(Davis & Gordon,1982年)。因此,量化和理解食品微观结构有助于推进食品科学的进步,实现食品产品的合理设计(Aguilera,2022年)。

尽管这一认识已存在数十年,但阻碍人们更深入理解食品在微观尺度上如何组织的因素是获取此类数据的可能性有限。最近在光学方面的进展使得能够在多个维度上观察和捕捉食品微观结构,包括超分辨率(长度)、体积、快速动态(时间)等多方面(Bonilla等人,2025年)。然而,如何将微观图像转换为描述食品材料结构特征的数据值仍然是挑战。现代数据科学的进步使得能够处理和分析大量数据,从而可以用于获取食品微观结构数据。本文重点讨论了从微观尺度食品结构图像中提取和处理数据的方法;然而,这些分析方法同样适用于其他长度尺度下生成 的图像。

在处理通过成像技术获得的数据时,必须区分图像数据和微观结构数据:图像数据指的是成像仪器生成的像素;而微观结构数据指的是图像中可见的物理微观结构。要从图像数据获取微观结构数据,通常需要经过以下步骤:图像采集、预处理、分割(包括对象合并)和量化。

**2. 图像预处理:清洗图像数据**
采集图像后,需要将感兴趣的结构与其他部分区分开来,这一过程称为分割。但在许多情况下,显微图像会受到多种因素的影响,从而妨碍分割。例如噪声、漂白、对比度低以及其他光学限制和像差。成熟的图像预处理技术(如强度变换和滤波器)可以减少这些系统引起的强度差异。本节展示了如何将图像预处理技术应用于食品微观结构的研究。

**2.1. 强度变换**
强度变换是一种预处理方法,通过修改原始像素值来提高对比度并突出特征。对比度拉伸或直方图归一化以及直方图均衡化都能增加像素之间的强度差异,从而使高亮部分变得更亮,阴影部分变得更暗,从而增强细节的可见性(Spring等人,1995年)。
对比度拉伸(也称为直方图归一化)将图像直方图扩展到所需的范围内(通常对于8位图像为0-255),同时保持其形状不变——从而增强前景和背景之间的对比度。对比度拉伸已应用于蛋清泡沫、麦芽糊精聚集体和乳制品凝胶的显微图像(Bonilla等人,2022年;García-Armenta等人,2016年;Glover等人,2019年)。通过“直方图均衡化”也可以提高对比度,例如在低脂酸奶的显微图像中应用该方法(Ciron等人,2012年)。在这种情况下,通过最大化动态范围的利用率来增强图像对比度(González等人,2009年;McReynolds和Blythe,2005年)。直方图均衡化生成的直方图比对比度拉伸更离散(图1)。在Ciron等人(2012年)的酸奶图像中,强度变换使得微流化牛奶样品中的蛋白质网络更加密集和相互连接。另一方面,随着微流化压力的增加,脂质相的脂肪颗粒大小减小。作者还发现,叠加图像中的黄色加深表明脂肪(绿色通道)和蛋白质(红色通道)的共定位以及它们之间相互作用的增强。

**2.2. 直方图匹配**
直方图匹配是利用参考图像的直方图形状来匹配后续图像的像素分布的方法。首先计算参考图像(CDF0)和第n个图像(CDFn)的累积分布函数(CDFn),然后根据CDF0更新第n个图像的像素值,从而使两张图像具有统一的对比度(Miura,2020年)。这在长时间曝光实验中观察到的荧光强度变化时非常有用。直方图匹配使得可以量化豌豆蛋白质网络,尽管由于加热和长时间激光照射导致蛋白质网络的对比度初始存在差异(Luecha等人,2024年)。这确保了在整个视频帧中获得的骨架能够代表豌豆蛋白质网络的微观结构。他们观察到了豌豆蛋白质浓缩物凝胶的动态形成过程,并能够量化网络分支和连接点随温度的变化。然而,豌豆蛋白质分离物的微观结构在加热过程中没有变化。

虽然强度变换通过处理前景信号来增强对比度,但也可以用相同的方法处理背景信号。背景校正通常用于解决采集过程中由照明模式引起的图像差异。例如“暗角现象”,即图像中心较亮而边缘较暗(Leong,2003年)。
平场校正可以修正各种乳制品凝胶和悬浮液的TEM显微图像中的阴影(Silva等人,2015年)。可以通过将感兴趣的图像除以经过背景校正的空白图像来进行校正。对于明场显微镜,可以使用没有样本的图像或失焦图像来进行平场校正。而对于荧光显微镜,使用与待检测对象相同荧光染料的载玻片作为平场图像(Sanderson等人,2018年;Zhang等人,2012年)。在所述的乳制品样品中,暗色颗粒被识别为均匀分布的蛋白质基质,而较亮的相位表示水相。也可以不通过获取平场图像来进行校正,这称为伪平场校正。使用同一原始图像,通过大核均值滤波器得到伪平场图像,从而大致掩盖图像细节,同时保留照明梯度(Schindelin等人,2012年)。Keshanidokht等人(2023年)利用伪平场校正假设CoARS(相干反斯托克斯散射)显微图像中的暗区不是空气孔洞,而是局部甘油。背景照明的归一化突显了高甘油浓度对暗区的影响。

另一种背景校正方法是使用“滚动球”算法来减去背景。通过从由滚动球确定的较大邻域内像素的平均值减去背景强度的变化来进行校正。使用半径大于感兴趣对象的球形结构元素,可以实现背景轮廓的归一化(Sternberg,1983年)。Zhu等人(2025年)设计了一种工作流程,利用滚动球方法监测原位木糖的蒸发结晶过程,同时校正图像的不均匀照明。通过边缘检测和强度阈值处理,成功生成了二值图像。由此可以量化结晶初期和末期晶体的面积分数、等效直径和长宽比。

**2.3. 模态运算**
形态学运算基于集合论,使用称为“结构元素”的小滤波器来增强图像中的微小特征(见第3.4.2节)。白色顶帽滤波器通过从原始图像中减去形态学开运算后的图像来突出比结构元素更小的明亮部分;黑色顶帽滤波器则相反,通过从闭合图像中减去原始图像来强调小的暗部分(González等人,2009年)。这些变换特别适用于校正不均匀的照明。例如,白色顶帽滤波器用于去除显微图像中的伪影并增强蛋清凝胶的蛋白质结构(Somaratne等人,2020年)。200×200像素的黑色顶帽滤波器还应用于牛奶悬浮液和凝胶的TEM显微图像,以增强小的暗结构并在Otsu阈值处理前抑制较宽的阴影变化(Silva等人,2015年)。该操作也用于校正苹果切片明场图像中的不均匀照明(Guillemin等人,2008年),有助于细胞的分割,并生成用于从相应免疫荧光图像中分离细胞壁强度信号的二值掩模。

**2.4. 滤波器**
另一种图像预处理方法是使用滤波器。滤波器是对像素数据进行的数学变换,用于抑制噪声、提取边缘、去除周期性元素或增强感兴趣的区域。

**2.4.1. 线性滤波器**
线性滤波器具有两个基本特性:叠加性和均匀性,即输入像素值组合后的响应等于各个输出像素值之和,并且卷积对像素值的数值缩放不敏感。高斯模糊滤波器通过对中心像素赋予更高系数来实现局部和减弱的模糊效果;这样可以减弱尖锐噪声,同时保持较大结构的完整性。它被广泛用于食品微观结构的研究,以去除不需要的噪声并促进2D和3D应用中的高效图像分割。高斯滤波器已应用于旋转模制饼干X射线微计算机断层扫描(X-ray μCT)图像的重建(Molina等人,2021年),观察到低糖样品的孔径更小、壁层更薄,这可能是导致其比对照组更粗糙和更硬的原因。这种方法也被应用于法国炸薯条的X射线μCT图像处理中,有助于可视化其外皮和孔隙结构(Gouyo, Rondet, Mestres, Hofleitner, & Bohuon, 2021)。同样,该滤波器也被用于涂有米粉或小麦面粉的油炸肉类模拟物的X射线μCT成像,以展示其外皮微观结构(Bhuiyan & Ngadi, 2024)。高斯模糊是使用共聚焦激光扫描显微镜(CLSM)捕获的麸质网络图像进行AngioTool基础分割的关键步骤。蛋白质网络的总体形状基于通过滤波器处理后的模糊图像(Bernklau, Lucas, Jekle, & Becker, 2016)。除了网络结构外,高斯滤波器还应用于油水乳液中的脂质滴和充气糖凝胶中的气泡等球形物体(Herremans et al., 2013, Teng et al., 2024)。另一种常用的线性滤波器是均值滤波器,它用每个像素及其邻居的平均值来替换该像素(Schindelin et al., 2012, Van Der Walt et al., 2014)。均值滤波被用于浓缩番茄悬浮液的光学显微照片和油炸火鸡块MRI图像的噪声去除(Bayod and Tornberg, 2011, Ghaitaranpour et al., 2024)。

2.2.2. 非线性滤波器
与线性滤波器不同,非线性滤波器不具备叠加性和均匀性特性,且通常是基于统计原理工作的。例如,中值滤波器(在Fiji和Python图像处理包中也称为去斑滤波器)通过用邻域的中值替换中心像素的强度值来有效处理噪声,尤其是椒盐噪声(Gonzalez et al., 2009, Schindelin et al., 2012, Van Der Walt et al., 2014)。通过中值滤波去噪已被证明可以增强几种食品微观结构研究中的结构可见性(Ando and Nei, 2024, Dreher et al., 2020, Herremans et al., 2013, Keshanidokht et al., 2022; J. D. Torres, Dueik, Contardo, Carré, & Bouchon, 2024)。维纳滤波器则通过估计含噪声估计图像与目标真实图像之间的最小均方误差来工作。空间变维纳滤波器通常通过采用自适应平滑效应来衰减高频成分,该效应取决于局部方差(Lim, 1990)。两项关于凝胶乳制品系统的研究利用维纳滤波和自相关分析提取了与孔间距、孔隙比例和典型蛋白质域长度相关的参数(Brüls et al., 2023, Glover et al., 2019)。通过应用维纳滤波器,异质区域(如蛋白质域边缘)的噪声会经历最小化的平滑效果,而均匀区域(如孔隙和空洞)的噪声则会减少。这被认为可以减少由于信噪比造成的虚假边缘。Kuwahara滤波器以在平滑过程中保留边缘和角落而闻名。使用这种非线性滤波器时,将对称正方形邻域中心的像素值替换为最均匀子区域的平均灰度值,即标准差最小的子区域(Papari, Petkov, & Campisi, 2007)。图2展示了高斯滤波器和Kuwahara滤波器之间的比较,以突出后者的边缘保留特性。非线性Kuwahara滤波器被用于蛋清泡沫微观结构的定量分析,以及其他处理步骤中,以增强气泡的轮廓(Bonilla et al., 2022a)。使用明场显微镜拍摄的蛋泡显微照片对比度较低(图2)。使用Kuwahara滤波器可以清理图像,同时保留气泡和液体基质之间的区别。

下载:下载高分辨率图像(486KB)
下载:下载全尺寸图像

图2. 分别使用Fiji对添加了高斯噪声的蛋白质网络图像(1a)实现了高斯平滑滤波器(1b)和Kuwahara滤波器(1c)。Kuwahara滤波器被用于蛋清泡沫显微照片(2b)上,以增强分割前的边缘。改编自(Bonilla et al., 2022),采用知识共享许可(CC BY 4.0)。 (3)用于过滤蛋白质凝胶CLSM图像的Haar小波及其结果分解的表示。改编自(I. C. Torres, Amigo Rubio, & Ipsen, 2012),版权(2011),经Elsevier许可。

非局部均值滤波器(NLM)被用于在不破坏原始结构完整性的情况下提高信噪比,特别是对于薄切片图像。NLM通过取具有相似特征(无论是强度还是颜色值)的像素的平均值并替换目标像素来进行去噪。这些相似的像素不一定是邻居;因此,会选择一个较大的非局部工作邻域来搜索相似像素。与Kuwahara滤波器类似,这种滤波算法旨在去噪时不丢失显著的边缘。相比之下,NLM不会产生使Kuwahara滤波后的图像具有独特绘画效果的块状伪影(图2)。与传统的高斯平滑滤波器相比,NLM可以保留固有的图像纹理,这一点在超声图像上得到了验证(Buades, Coll, & Morel, 2005)。这些滤波器广泛用于食品微观结构的X射线微计算机断层扫描图像处理中(Ando and Nei, 2024, Derossi et al., 2019; Riley et al., 2023; Riley, Verma, Verboven, Nicolai, & Delcour, 2024)。最后,非线性锐化滤波器常用于增强图像中的细节辨识,特别是在细节明显的情况下强调细节和边缘(Ciron et al., 2012, Derossi et al., 2019, Molina et al., 2021)。在食品微观结构的定量分析中,锐化滤波器很少单独使用(Germain & Aguilera, 2012)。它们通常作为机器学习或深度学习管道中的特征提取器,用于分割和分类。

2.2.3. 频率域中的滤波器
迄今为止讨论的线性和非线性滤波器都是针对空间域的滤波器,但也可以在频率域中对成像数据进行滤波。在这种情况下,图像数据被转换为频率分量,然后在数据转换回空间域之前对频率表示进行滤波。频率域滤波的优点是计算效率高,尤其是使用快速傅里叶变换(FFT)时,以及能够避免空间滤波可能产生的振铃伪影。频率域滤波可以处理图像的平滑和锐化。频率域滤波还广泛用于去除周期性元素,如莫尔纹图案和重复线条(Gonzalez et al., 2009, Russ and Neal, 2018)。低通滤波器(理想滤波器、巴特沃斯滤波器或高斯滤波器)用于去除高频噪声,而高通滤波器通过增强高频成分来突出边缘。带通滤波器可以通过结合这两种方法来隔离特定大小的特征。在乳清蛋白分离物-海藻酸盐(WPI/NaAlg)凝胶的图像分析中,使用了低通巴特沃斯滤波器来在分割前衰减微粒边界(Leon, Medina, Park, & Aguilera, 2016)。作者得出结论,碳酸钙的浓度影响了WPI/NaAlg微粒的大小。较高的海藻酸盐浓度产生了更小且更明确的微粒,这也影响了凝胶的机械性能,因为相同样品的压缩应力值高于海藻酸盐浓度较低的样品。频率域中的另一种滤波器是小波变换。与傅里叶变换相比,小波变换的优点是保留了时间信息。小波变换使用一系列小波(Haar、Daubechies和墨西哥帽)来近似给定函数,类似于傅里叶变换使用一系列正弦波。一系列高通和低通滤波器将输入图像分为四个子图像(最终近似、垂直细节、水平细节和对角细节)进行第一次小波变换(Gonzalez et al., 2009)。在 Different yogurt microstructures 上实施了一步Haar小波变换,以在对其进行分形分析之前压缩图像(I. C. Torres et al., 2012)。结果分解后的图像在保持重要信息的同时减小了图像的大小(图2)。使用从酸奶微结构的分形维度派生的局部指数可以区分全脂、低脂和含有微粒的酸奶样品。

2.3. 反卷积
光学成像系统会受到由衍射、干涉和其他光学系统固有偏差引起的模糊效应的影响。反卷积是一种数学上逆转光学系统引入的卷积效应的过程。反卷积算法(如Richardson-Lucy或Wiener滤波)利用描述光学系统的点扩散函数(PSF)。PSF可以根据已知的光学参数理论估计,或者使用荧光珠实验测量(Scientific Volume Imaging B.V., n.d.)。在食品微观结构成像中——折射率不匹配和多相系统很常见——准确的PSF估计对于恢复样品的“真实”结构至关重要。反卷积已被用于处理激射发光耗尽显微镜下的蛋清凝胶图像(Bonilla et al., 2022)。借助反卷积,作者能够清晰地可视化约50-100纳米的蛋白质纳米聚集物,这些聚集物构成了蛋清凝胶的结构。反卷积是使用3D CLSM验证人乳脂肪球体的重要步骤,与激光衍射测量相比(Verveld, De Wolf, Legtenberg, Knop, & Bosschaart, 2024)。对于每个感兴趣的球体,使用Matlab函数(deconvlucy)在物体的中心Z位置对图像进行反卷积。反卷积图像及其原始图像用于获得体积 milk 尺寸分布。3D成像的结果提供了比传统激光衍射实验更准确的尺寸分布,用于单分散 milk 标准。盲反卷积技术能够使用偏振光显微镜检索模型脂肪系统的体积晶体网络的分形尺度特性(Litwinenko, 2004)。

3. 图像分割:像素分类
图像分割是从处理后的图像中提取有用信息的关键步骤。它将连续图像转换为离散区域,通常产生一个二值图像,将各个部分编码为对象(前景)或背景;因此,这个过程也被称为二值化。在大多数工作流程中,分割通常作为量化食品材料微观结构特性的入口。主要目标是通过将感兴趣的结构从背景中分离出来,从而简化进一步的定量分析,减少图像的复杂性——本质上是数据集的维度。

3.1. 基于强度的分割
3.1.1. 强度阈值
阈值分割是最简单且最广泛使用的分割方法。它使用图像的一个基本属性——强度直方图——作为分割的指导。例如,一种基本的迭代方法从一个初始猜测值(Tinitial)开始,即强度阈值,然后通过下方和上方像素组的统计参数进行调整,直到达到收敛。然后可以在整个图像上实施单一阈值进行二值化——其计算方法已经通过不同的方法进行了优化。在Otsu方法的情况下,最大化背景和前景两类像素之间的方差。当直方图不是明显的双峰时,基于熵的方法(例如最大熵、Renyi熵)或三角形算法(为偏态直方图设计)可能更有效。对于没有明显山谷或突出峰值分离的直方图,Otsu方法比基本阈值算法表现得更好(Gonzalez et al., 2009)。Otsu阈值分割算法是分析食品微观结构最广泛使用的二值化技术之一。结合Otsu阈值分割的法国炸薯条的X射线μCT图像显示,较高的外皮孔隙率转化为样品更高的酥脆度(Gouyo et al., 2021)。此外,与热空气炸薯条相比,油炸薯条的较小孔径和非均匀孔隙大小分布有助于提高其酥脆度。相反,研究表明,涂有米粉或小麦面粉的油炸肉类模拟物(大豆蛋白分离物或基于麸质的)具有均匀的孔隙分布,产生了更好的性能(Bhuiyan & Ngadi, 2024)。米饭裹面肉类的外壳微观结构抑制了质点的重新分布,从而提高了这些样本的硬度、脆性和酥脆度。通过对奶油蛋奶冻、海绵蛋糕、布丁和饼干样品中的脂滴进行粒度分布分析,发现颗粒分割对脂质消化的动力学和程度有显著影响(Hiolle等人,2020年)。这项研究表明,尽管这些复杂的食品成分相同,但其微观结构在宏量和微量营养素的生物利用度方面起着重要作用。同样,也研究了气泡的大小分布及其与充气糖凝胶和蛋清泡沫的质地和流变特性的关系(Bonilla等人,2022年;Herremans等人,2013年)。Otsu阈值算法也被用于蛋白质网络的分割,正如Glover、Bisgaard等人(2019年),Li等人(2023年)和Silva等人(2015年)所展示的那样,他们都在不同的处理条件下检测了基于乳制品的凝胶。Silva及其同事能够利用Otsu阈值生成的二值图像来估算面积分数、孔隙率和单位面积的边界长度。该算法还被用于将CLSM图像二值化,以研究它们的微观结构与胃蛋白酶扩散性的关系(Somaratne等人,2020年)。他们发现,在pH=5时,与pH=9的凝胶相比,更开放(曲率更高)和异质性更强(颗粒大小和颗粒间聚集距离分布更广)的蛋清凝胶使得胃蛋白酶在基质中的移动速度更快。然而,不同的系统和成像模式可能需要不同的二值化算法。Floury等人(2018年)使用三种不同的阈值算法(Otsu、Triangle或Yen)来跟踪脱脂牛奶凝胶的消化过程。通过时间序列同步辐射深紫外荧光显微镜获得的每个视频帧中,都生成了感兴趣区域及其对应的二值图像作为指导。选择分割算法是基于对颗粒轮廓的视觉检查。Otsu算法对于酸凝胶(pH=2)在盐酸消化过程中的分割效果最清晰,而Triangle算法则更适合于同样处理条件下的凝乳酶凝胶。最终,Yen分割算法被用于用胃蛋白酶消化的酸凝胶和凝乳酶凝胶。因此,他们得出结论,胃蛋白酶对酸凝胶的水解速度大约是凝乳酶凝胶的10倍。这是由于凝乳酶凝胶发生了重大的结构变化,导致颗粒变得更加密集。阈值算法被用于通过偏振光显微镜观察到的不同脂肪微晶的分割,如图3所示。这些算法通过峰值信号噪声比、结构相似性指数、晶体尺寸误差等指标进行了评估。从这些性能指标来看,默认算法、IsoData和Huang算法在所有测试的脂肪微晶分割案例中表现最佳(Xiong、Qi和Zhang,2024年)。阈值设置不当会导致对晶体像素的高估或低估,从而影响其分形维数的计算。在其他矩阵中,如蛋白质网络中,如果使用的分割技术不足,也可能出现类似情况,如图3所示。

**图3:** (1) 使用FIJI内置阈值算法处理的脂肪晶体偏振光显微图的二值图像。(改编自Xiong等人,2024年,版权所有©2024年,经AOCS许可。)(2) 通过Triangle、Otsu和Ilastik分割的乳清蛋白凝胶的显微图,测量了总长度。(3) 使用基于迟滞的过滤器对麸质蛋白进行的图像分割(作者未发表的工作)。(4) 使用K-means(k=3)对FLIM显微图进行分割,将图像分为蛋白质(橙色)、蛋白质-空隙界面(绿色)和空隙空间(蓝色)。改编自Mishra等人,2025年,采用知识共享许可(CC BY 4.0)。(5) 在乳液中定量氧化和未氧化脂质的工作流程,使用了StarDist分割技术。改编自Brüls-Gill等人,2024年,采用知识共享许可(CC BY 4.0)。

由于基于强度的分割通常只能区分两类——“前景”和“背景”,因此在复杂矩阵中二值图像无法区分不同成分(例如蛋白质与脂质)。只有使用特定技术(如荧光成像)才能区分成分,其中不同的染料标记不同的成分。在这种情况下,为每种染料生成一个图像通道,并可以为每个通道使用一个阈值。在其他情况下,可以通过多级阈值分割或更复杂的分割技术来实现单个图像中的多成分分割(Liao、Chen和Chung,2001年)。

3.1.2 基于形态的分割
基于边缘的分割利用了边缘对应于强度快速变化的事实。例如,可以使用Canny边缘检测算法来提取边缘图,该算法本质上是一个基于迟滞的过滤器。迟滞分割需要高阈值和低阈值作为指导。迟滞分割基于像素的连续性,因此Tlow < 强度 < Thigh的像素也被归类到对象类中(Condurache和Mertins,2012年)。然而,仅凭边缘无法形成完整的对象,因此通常会使用膨胀和填充空洞等形态学操作来构建封闭对象。图3展示了一个基于麸质的面团切片上的迟滞分割示例。

基于区域的分割是另一种基于形态的分割方法。它从感兴趣对象内部的一个或多个手动指定的种子点开始。然后算法通过迭代包含满足特定相似性标准的相邻像素来“生长”这些区域,例如属于特定的强度或RGB值范围。区域生长特别适用于对象和背景对比度低的图像(Gonzalez等人,2009年)。典型的区域生长分割大多被应用于其他食品相关识别系统,用于饮食评估(Dalakleidi、Papadelli、Kapolos和Papadimitriou,2022年)。为了分析热处理肉乳液的微观结构,使用马尔可夫随机场建模和图割能量最小化结合的方法对多模态X射线µCT获得的图像进行了分割(Einarsdóttir等人,2014年)。为了进一步分离接触的对象,开发了一种区域生长标记算法。该方法消除了距离背景一定欧几里得距离以下的体素,形成了一个新的二值体积。在这个新空间中,从距离背景最远的体素开始,向对象-背景边界扩展一个填充算法。

3.2 机器学习分割
除了经典算法之外,机器学习(ML)已成为图像分割的关键工具。传统的机器学习方法(如k-means聚类、支持向量机和随机森林)基于从图像像素中提取的特征(例如强度、纹理或颜色)进行分割。这些基于规则的算法结合了注释数据来学习不同图像区域之间的分类边界。这些方法的优势在于它们的相对简单性和可解释性,但可能需要手动选择特征和调整参数。

3.2.1 k-means聚类
这种无监督方法根据k-means(也称为质心)与待分类的未标记数据点之间的欧几里得距离最小化来将像素分成簇。当k设置为2时,该方法可以有效地区分前景和背景(Russ和Neal,2018年)。由于荧光寿命值的显著差异,Mishra等人(2025年)成功使用k=3聚类算法分割了酪蛋白网络、界面和空隙空间(图3)。使用k-means聚类使他们能够分离蛋白质-空隙界面,并研究果胶添加对系统中分子相互作用的影响,这是使用简单阈值难以实现的。

3.2.2 随机森林
随机森林是一种非线性分类算法,它使用决策树集合对像素进行分类。决策树从一个包含全部像素的根节点开始,然后将其分割成更小的子组(分支),每次分割都使子样本更加均匀,直到达到某个停止分割的标准。熵和信息增益是确定哪些属性应用于分割节点的常用方法(Berg等人,2019年)。随机森林分割已广泛用于食品微观结构分析,特别是通过Ilastik(Berg等人,2019年)或ImageJ/FIJI中的Trainable WEKA分割插件(Arganda-Carreras等人,2017年)等交互式软件。

Ilastik作为一种重要的生物成像分析工具,也在食品微观结构应用中崭露头角。在训练过程中,用户可以从一组通用特征(平滑的像素强度、边缘特征和纹理)中选择,然后通过过滤器提取这些特征。除了这些特征外,分类器还考虑用户提供的正确类别分配。在像素分类的工作流程中,可以用不同颜色标注感兴趣的区域。Ilastik默认使用100棵树的随机森林进行语义分割。例如,它可以与分水岭算法结合使用,以获得图像的实例分割。Ilastik已成功用于分割球形对象,如细胞、孔隙、脂肪滴和空气细胞(Keshanidokht等人,2022年;Lourenco等人,2022年;Luecha等人,2024年;Ma等人,2020年;Martínez-Sanz等人,2020年;Riley等人,2024年)。Luecha等人(2024年)仅需要标注视频中的6帧(分别对应于30°C、40°C、50°C、60°C、70°C和80°C的凝胶化关键温度),即可快速完成豌豆蛋白网络的分割。之后,所有其他帧都使用训练好的模型自动分割。

Waikato环境用于知识分析(WEKA)是一个支持实现各种机器学习算法(如随机森林、K最近邻、朴素贝叶斯和支持向量机)的数据挖掘工具。为了扩展其在图像分析中的能力,开发了适用于ImageJ的Trainable Weka Segmentation(TWS)插件。该工作流程应用不同的过滤器提取特征,然后将它们转换为向量。使用标记的训练数据,所选算法学习图像中先前未标记像素的隐藏模式(Arganda-Carreras等人,2017年)。TWS被用于研究不同蘑菇物种的微观结构和质量(Hu等人,2022年;Xue等人,2024年)。

基于图割的分割被用于通过衍射X射线计算机断层扫描(PXCT)获得的3D乳液混合物微观结构的分割。作者对数据进行了建模,以获得图割所需的马尔可夫随机场参数,从而找到最佳的分割方案,将图像划分为空气、微液滴、水和纤维素相。研究发现,在棕榈仁油-水乳液中,由脂质聚集体组成的扩展脂质网络占据了大部分脂质,未连接的脂滴数量较少。文章还报告了PXCT分析中所有组分的定量电子密度,这有助于识别样品中的已知相以及检测可能的异常(Nielsen等人,2016年)。

3.3 深度学习在图像分割中的应用
近年来,深度学习彻底改变了图像分割,尤其是在处理复杂的异质图像(如食品微观结构)方面。深度学习具有多种优势,例如高精度、能够从原始数据中自动学习特征以及处理大规模数据集的鲁棒性。卷积神经网络(CNN)是最常用的深度学习模型之一,用于计算机视觉任务。AlexNet、YOLOv5和VGGNet是已建立的深度学习架构,已在食品相关研究中得到应用。这些模型已被有效用于根据微粒的形态特征对微粒和乳液滴进行分类(Z. Huang等人,2022年;Yu、Lu、Han和Xue,2021年)。AlexNet、VGG-16和VGG-19模型在扫描电子显微镜(SEM)图像中对喷雾干燥微藻油的微胶囊进行分类的准确性进行了比较(Yu等人,2021年)。通过迁移学习,他们发现VGG-19模型的性能最好,并被用作分类聚合与非聚合、完整与破损微粒的基准。深度学习模型在这些应用中的表现很大程度上取决于训练数据的可用性和质量。从零开始构建模型通常需要一个大型数据集。然而,迁移学习提供了一个实用的解决方案,允许在较小的数据集上对预训练模型进行微调,从而减少了对大量数据收集的需求。U-net架构已被用于食品微观结构分析,主要通过两种现有算法:StarDist和Cellpose。在StarDist中,物体被表示为星形凸多边形,可以预测物体存在的概率以及从中心点到物体的径向距离。与边界框相比,这种几何形状更加灵活,更适合近似圆形物体(Schmidt, Weigert, Broaddus, & Myers, 2018)。StarDist最初是为细胞分割开发的,但也有效用于荧光成像中类似圆形的物体(如油滴)的分割。这项技术已被几篇论文改进,用于研究不同乳液系统中的脂质滴(Brüls-Gill et al., 2024; Yang, Verhoeff, Merkx, Van Duynhoven, & Hohlbein, 2024)。在Cellpose中,这种分割方法受到了分水岭算法的启发,模拟了一个从中心向外辐射的扩散过程。然后训练U-Net神经网络,不仅预测像素分类(前景与背景),还定义X和Y空间梯度(流场)。在新图像的分割过程中,使用预测的梯度矢量场来引导像素朝向固定边界作为终点进行重建。Cellpose最初被开发为一种通用的细胞分割算法,后来改进为包含人类参与的过程以微调分割性能(Pachitariu & Stringer, 2022),最近还加入了图像修复算法以提高在质量不佳的图像中的性能(Stringer & Pachitariu, 2025)。由于Cellpose是在细胞数据上训练的,因此该模型在球形物体的分割上也表现最佳。这已在一些与食品微观结构相关的论文中得到验证(Bonilla et al., 2022, Chen et al., 2025)。上述的Cellpose、StarDist和Ilastik最近与一个名为MIDAS的改良版本进行了比较,后者可以在多个直径输入下运行(Saalbrink et al., 2025)。尽管这些方法在形状上与细胞相似,但显然,一些具有强烈尺寸分布的乳液和泡沫在使用传统机器学习算法进行分割时表现较差。这可能是因为细胞训练集的尺寸变化有限。与仅基于少量图像在Ilastik中训练的强度阈值分割模型相比,MIDAS显示出更好的结果。这项研究清楚地表明,分割技术的选择及其准确性对微观结构测量有重大影响。关于基于U-Net的乳液和泡沫分割与分析的应用将在相关章节中进一步讨论。另一种用于图像分割的CNN深度学习方法是基于掩模区域的CNN,它可以预测物体边界框、类别标签和分割掩模。这种两阶段过程——首先提出感兴趣区域,然后对其进行精细化——已被证明对于分割重叠物体非常有效(He, Gkioxari, Dollár, & Girshick, 2017)。为了分析高速气体泡流中的重叠气泡,开发了一个名为“Bubble Boundary R-CNN”的多任务模型。该模型可以分割气泡、重建其形状并提取用于进一步分析的相关特征。它包含了一个保边掩模头部,该头部使用拉普拉斯变换的掩模来增强对细微边界细节的捕捉(Kang et al., 2025)。由于这种方法能够准确勾勒出即使是小且部分遮挡的气泡,因此对于分析双乳液特别有价值。其他深度学习方法,如视觉Transformer(ViTs),现在也被用于基于显微镜的分割任务(Archit et al., 2025)。虽然ViTs可能需要更大的数据集和更多的计算资源,但它们具有更高的分割精度和捕捉全局上下文信息的潜力。例如,结合CNN和注意力机制或Transformer架构的方法在医学图像分割中显示出了有希望的结果,并且可能对食品科学应用也有益。3.4. 对象整合虽然已经处理了将感兴趣的结构区分为主观和背景的问题,但在量化之前,这些对象可能仍需要进一步处理。例如,在对象重叠或接触的情况下需要这些操作。这些操作通常被称为形态学操作,因为它们处理分割对象的形状。形态学运算符需要一个结构元素或核来指定操作的过程,这最好通过集合论来描述(Gonzalez et al., 2009)。3.4.1. 分水岭分水岭是一种常用于物体紧密排列或接触的图像的形态学运算符,例如泡沫中的气泡或乳液中的油滴;传统的阈值分割可能无法分离单个结构。为了构建对象边界,分水岭方法将图像视为一地形表面,其中像素强度表示海拔高度(图4)。水“流入”集水区,而流域之间的边界或“坝”则划定了对象(Soille & Vincent, 1990)。尽管有效,但在嘈杂的图像中分水岭算法容易发生过分割,因此通常需要预处理步骤,如高斯平滑或使用标记(也称为种子)来控制分割的范围。下载:下载高分辨率图像(521KB)下载:下载全尺寸图像图4。(1)分水岭算法的图形表示。根据Creative Commons许可(CC BY 3.0)改编自(Fisher, 2014)。(2)根据3D分水岭可视化挤压的螺旋藻-淀粉泡沫图像中的单个孔隙。改编自(Martínez-Sanz et al., 2020),版权(2020),经Elsevier许可。(3)通过Otsu阈值分割的充气糖凝胶气泡,并通过分水岭分离成单个对象。改编自(Herremans et al., 2013),版权(2013),经Elsevier许可。(4)酪蛋白网络的处理步骤包括对悬挂端和环进行修剪,这是一种形态学操作。改编自(Brüls et al., 2024),版权(2023),经Elsevier许可。分水岭算法有助于量化单个对象及其属性,例如计算脂肪滴的数量和构建它们的大小分布,这对于乳液微观结构研究可能很有兴趣(Dreher et al., 2020, Keshanidokht et al., 2022)。虽然分水岭算法广泛用于细胞、气泡(图4)和滴状物等球形对象,但它也可以用于分离长形物体,如木糖晶体(Germain and Aguilera, 2012, Herremans et al., 2013, Rahman et al., 2018, Zhu et al., 2025)。通过分水岭划定对象轮廓的方法可以与其他二值化技术一起使用,无论是简单的阈值分割还是基于机器学习的方法。分水岭算法也可以扩展到3D。这在螺旋藻-淀粉泡沫的一些3D微观结构描述符(孔隙率、比表面积、细胞壁厚度、连接密度等)的估计中有所描述(Martínez-Sanz et al., 2020)。比较了两种从3D重建中估计孔径大小的方法。第一种方法仅使用Pore3D软件库,该库结合了包括分水岭分割在内的图像处理步骤。第二种方法使用ImageJ中的MorphoLibJ包中的“Distance Transform Watershed 3D”算法进行分割,然后使用Pore3D进行量化(图4)。两种方法都支持这样的结论:螺旋藻含量的增加导致孔径整体变小。3.4.2. 粒度测量形态学开运算和闭运算结合了基本侵蚀(从图像边界移除像素)和膨胀(向图像边界添加像素)操作。开运算是为了平滑轮廓、去除细小突起并打破物体内的狭窄连接。另一方面,闭运算是为了填充小孔和间隙,并连接附近的对象。粒度测量是一种使用逐渐增大的结构元素进行连续形态学开运算来估计图像中粒子大小分布的操作。这个过程类似于物理筛分,即将粒子通过逐渐增加孔径尺寸的筛子。通过测量连续开运算之间的面积差异,可以识别出主导的粒子大小。还可以通过将结构元素大小与打开的像素总面积相对应来绘制累积分布函数(CDFs)(Silva et al., 2015, Soille, 2004)。各种研究使用这种方法来确定特征尺寸,例如冻干胡萝卜孔隙的中点(Voda et al., 2012)或乳制品系统中蛋白质聚集体的加权平均值(Silva et al., 2015)。Lassoued, Babin, Della Valle, Devaux和Réguerre(2007)通过在2D和3D面包结构分析中比较粒度测量进一步证明了其有效性。3.4.3. 形态学细化(骨架化)和修剪细化或骨架化在基本形态学中涉及一系列侵蚀和开运算,将对象简化为细长的连接结构(Gonzalez et al., 2009)。然而,单独的形态学运算无法保证连通性,因此将在第4.3节中讨论其他骨架化算法的创新。与此相关的是,修剪作为骨架化的一个整合步骤,用于去除由笔画不规则性引起的“寄生组件”。这是通过使用检测端点的结构元素进行细化,然后通过条件膨胀来重建主要结构来实现的。如图4所示,骨架的端点连接分支(悬挂端)通过“修剪端”函数被移除,以确定酸奶凝胶网络的连通性(Brüls et al., 2024)。4. 图像量化:测量食品微观结构迄今为止讨论的所有图像处理步骤的最终目标都是量化在显微镜下研究的食品或其他软物质的微观尺度的结构元素。在结构元素被分割和识别之后,最后的图像分析任务集中在量化它们、测量它们的空间分布及其形态学特征上。根据这些量化结果,我们将食品微观尺度的结构元素分为三类:1. 滴状物和聚集体,2. 网络,3. 纤维。4.1. 滴状物、细胞和聚集体滴状物和聚集体的分析基于它们是具有高度中心对称性的单个元素这一基础。食品微观结构中的常见例子包括乳液中的胶体滴、泡沫中的气泡或球形蛋白质聚集体。乳液和泡沫是滴状物的分散体,其微观大小、形状和空间排列决定了霜、酱料及相关食品产品的功能特性(Dickinson, 2010)。根据胶体滴的排列方式,它们会形成球形或多边形阵列,并表现出时间依赖的结构变化,包括聚结、奥斯特瓦尔德熟化或絮凝(Langevin, 2023)。因此,显微镜成像提供了宝贵的见解,因为这些事件可以从滴状物的原位结构观察到并量化,而广泛使用的技术,如光散射,则需要大量稀释(Balcaen, De Neve, Dewettinck, & Van Der Meeren, 2021),尤其是对于高度密集的系统,它会忽略天然的胶体排列和时间依赖的变化。虽然第3.3节讨论的分割机器学习算法并非专门为分割滴状物和聚集体而开发,但它们与密集细胞的相似性使它们可以用于此类情况。StarDist结合一个内部的Python脚本被用来研究脂质氧化的定位以及引入生育酚和过渡金属的影响(Brüls-Gill et al., 2024)。使用BODIPY 665/676(氧化脂质:λex = 561 nm 对比非氧化脂质:λex = 640 nm)的荧光强度和蛋白质的自荧光通过CLSM监测了14天内的氧化情况(图5)。由此得出结论,乳液水相中的蛋白质-脂质聚集体比油滴相中的脂质氧化得更快。下载:下载高分辨率图像(699KB)下载:下载全尺寸图像图5。(1)乳液的原始图像,以及根据每个区域的光子计数量化的氧化滴积累和蛋白质及脂质氧化的滴大小依赖性。根据知识共享许可(CC BY 4.0)改编自(Yang等人,2024):(2)在0%和90%压缩下半硬质奶酪的孔径对数正态分布的比较,以及相应的共聚焦显微镜图像,显示连续的蛋白质基质(红色)和渗透的脂肪滴(绿色)。图像是根据知识共享许可(CC BY 4.0)改编自(De Boer等人,2025)。(3)新鲜制备的乳液和经过摩擦学处理后的乳液中内部水滴和外部油滴大小分布的变化。根据知识共享许可(CC BY 4.0)改编自(Tenorio-Garcia等人,2026)。还有两篇论文展示了使用StarDist对乳清蛋白分离稳定性大豆油-水乳液和蛋黄酱中的脂肪滴进行分割的方法(Yang等人,2020;Yang等人,2024)。这两项研究得出结论,脂肪滴的大小调节了脂肪和蛋白质的氧化速率,如图55所示。同样,StarDist也被用来分割硬质奶酪中的微观空穴,以将宏观变形与蛋白质、脂肪和空穴的微观排列联系起来(De Boer等人,2025)。通过共聚焦显微镜图像获得了这些半硬质奶酪在0%和90%压缩下的孔径分布(图5)。另一种用于细胞分割的算法是Cellpose,之前在第3.3节中提到,它已被用于量化蛋清泡沫中的气泡(Bonilla等人,2022)。Cellpose经过改进,以克服基于细胞的分割算法在处理大小分布差异较大的滴液时存在的一些内在限制。Microscopic Droplet AnalysiS(MIDAS)工作流程被开发出来,可以在多个定义的直径设置下迭代执行基于Cellpose的分割(Saalbrink等人,2025)。除了提供更准确的滴液分割外,MIDAS还通过实现一个名为‘SimpleITK’的Python库引入了量化步骤。具体来说,该库为每个分割对象提供了以下信息:x-y质心位置、面积、Ferret直径、扁平度、圆度、长圆形、球形周长、球形半径、x方向的椭球体直径、y方向的椭球体直径。在非稀释系统中量化滴液可以为研究它们的排列和形状提供可能性,而不仅仅是提供简单的颗粒大小分布。作者强调了通过轮廓或质心数据的邻近分布可以研究滴液界面间的聚集或排斥效应。

双乳液是一种更复杂的系统,其特征是嵌套的滴液,MIDAS也可以轻松对其进行量化,因为它对非均匀大小滴液的分割采用了全面的方法。研究表明,该工作流程可以从双Pickering双乳液的显微镜图像中分别量化内部和外部滴液的大小分布(Saalbrink等人,2025)。在相关研究中,Feret直径和大小分布与含有亚微米乳清蛋白微凝胶和可可脂油凝胶基晶体的双Pickering水-油-水双乳液的摩擦学性质有关(Tenorio-Garcia等人,2026)。油水乳液中的较小滴液可以填充舌乳头的间隙,从而减少摩擦并有效提高润滑效率。另一方面,对于双乳液来说,内部滴液大小分布的差异对润滑性能影响不大。润滑机制显著受到油滴和油凝胶内部晶体网络的存在的影响。记录了摩擦学测试前后内部水滴和外部油滴(红色)分布的变化(图5)。

一个专门为显微泡沫图像开发的开源工作流程显示,针对食品特定图像分析工具箱的开发趋势正在增加。一个基于Python的脚本pyFoam被用来分析豌豆白蛋白泡沫的图像,并结合小角散射和导电性测量(Li, Lamolinairie, Chiappisi, & Corredig, 2025)。通过这种方法,可以直接比较不同pH值下泡沫的大小和气泡数量的变化与通过导电性测量的排水情况。同样,还可以比较小角散射和图像分析之间的比表面积、膜厚度和平坦边界。

网络是指由蛋白质(如面团中的麸质或牛奶凝胶中的酪蛋白微胶束)、多糖(如水凝胶中的果胶)、脂肪晶体(如巧克力中的)或这些成分的组合形成的三维互连矩阵(Bonilla等人,2020;Cheng等人,2025;S. Huang等人,2025;Mishra等人,2025;Xia, Drositi, Czaja, Via, & Ahrné, 2024)。在网络的情况下,感兴趣的特定形态描述符是它们的互连性和网络网格的不同度量,因为通过表征和量化这些网络特征,可以分析它们如何影响食品的质地、稳定性和性能。这最终将导致产品配方和加工的优化。骨架提取,也称为骨架化,是指生成对应于结构中心的细线的过程。最常见的方法是保留拓扑的并行(细化)算法,它通过欧几里得距离变换(EDT)逐步从二值结构中移除像素,直到原始对象中心只剩下一条像素线(Lee & Kashyap, 1994)。另一种方法称为最小成本路径,包含在TEASAR算法中(Sato等人,2000),它通过欧几里得距离变换(EDT)识别峰值体素,然后提取骨架路径。对于量化,算法首先将骨架视为一个图,从中识别图中的节点或点。通过体素或像素扫描方法,分析两种类型的节点:末端节点(只有一个邻居的体素)和交汇节点(有两个以上邻居的体素)。然后,识别连接两个节点的路径作为分支或边缘。最后,计算网络参数,包括分支数量及其长度、交汇节点、端点、扭曲度(分支偏离直线的程度)、聚类性和中心性度量等,所有这些都指示了网络的连通性。使用图论和图像分析的算法开始量化。这里列出了一些最具代表性的算法。

使用Skeletonize 3D Fiji插件可以通过细化算法对二值图像进行2D或3D骨架化(Lee, Kashyap, & Chu, 1994)。然后,使用另一个Fiji插件Analyze Skeleton可以量化参数,例如分支长度和数量、交汇节点、端点和扭曲度。具体来说,Analyze Skeleton利用体素连通性和拓扑进行分析。这已被用来量化酪蛋白网络的互连性。如图所示,使用牛奶凝胶的共聚焦图像来测量每个骨架的总长度和分支数量,然后利用分支数量与总分支数量的比例来量化酪蛋白网络的链接密度(Brüls等人,2024;Foroutanparsa等人,2023)。同样的方法也被用来研究纤维添加对酪蛋白凝胶化的影响(Wang, Chen, Zhong, Xu, & Ren, 2025)。同样,如图6所示,Analyze Skeleton被用来测量加热过程中豌豆蛋白凝胶的分支和交汇节点的动态增加(Luecha等人,2024)。此外,Skan是一个Python库,它使用稀疏矩阵,允许在2D或3D中进行更广泛的网络量化。除了常见的分支描述符外,它还存储有关分支类型、节点坐标和节点-分支连通性的信息。在图6中,使用Skan量化了乳清蛋白水凝胶的3D图像,其中分支根据其长度进行颜色编码(Nunez-Iglesias, Blanch, Looker, Dixon, & Tilley, 2018)。

下载:下载高分辨率图像(1MB)
下载:下载全尺寸图像

图6. (1)通过链接密度评估酸奶中酪蛋白网络的连通性。改编自(Brüls等人,2024),版权(2023),经Elsevier许可。
(2)豌豆蛋白凝胶化过程中骨架化图像的分支和交汇节点数量的相对变化。根据知识共享许可(CC BY 4.0)改编自(Luecha等人,2024)。
(3)利用AngioTool量化了体积水浓度对微观结构(例如蛋白质面积百分比、交汇节点密度)的影响。改编自(Brüls等人,2024),版权(2023),经Elsevier许可。
(4)使用Skan对乳清蛋白水凝胶进行3D分割和网络分析(分支长度和节点量化)(作者未发表的工作)。

Angiotool是一个用于分析显微图像中血管网络的Fiji插件(Zudaire, Gambardella, Kurcz, & Vermeren, 2011)。它可以用于灰度或彩色图像,其中结构已标记。通过多尺度Hessian基增强滤波器检测管状结构,然后应用骨架化细化算法。量化的参数包括血管(分支)数量、面积、长度和指数、交汇节点(密度)和端点、血管(网络)密度以及通过快速盒计数算法得到的lacunarity。后者提供了关于网络内部异质性的信息。在食品矩阵中,该工具已被用于量化不同麸质网络的强度和凝聚力(Bernklau等人,2016;Bonilla等人,2020;Lucas等人,2018;Lucas等人,2019;Ogawa和Matsumura,2021),如图6所示。在各种参数中,分支率描述了网络连通性,端点率描述了网络的弱化,而lacunarity则归因于网络间隙的数量和大小,这描述了网络的不规则性(Lucas等人,2019)。

其他基于分割和细化骨架化的插件在量化分支、交汇节点、端点和扭曲度之前包括Neurite tracer(Pool等人,2008)、NeurphologyJ(Ho等人,2011)和simple neurite tracer(Longair等人,2011;Tavares等人,2017)。这些插件专门用于神经元和血管的研究,但它们对食品网络的类似详细量化在食品科学中也有用,无论是2D还是3D网络。

纤维是指细长的丝状结构,它们可以单独存在,也可以作为互连矩阵的一部分,其中每个单独的纤维与其他纤维分开。例如,在基于植物的挤出物或其他食品如马苏里拉奶酪中可以找到纤维(Patel, Barnés-Calle, Rinnan, Engelsen, & Van Den Berg, 2025)。在材料科学中,根据它们在不同长度尺度上的位置, filament、fibril和fiber之间有所区分。Filament代表最小尺度,主要由排列成丝状的蛋白质分子组成。例如,肌动蛋白丝(厚15纳米)和actin丝(厚7纳米)。Fibril稍大一些,通常约为1微米厚,可以由多个filament聚集而成。纤维是最大的结构,有时肉眼可见。它们被描述为具有方向性的线状或绳状结构,例如在结构化的肉制品中找到。在图像分析中,图像的长度尺度并不重要,因此可以按照相同的步骤进行量化。出于同样的原因,这些术语经常可以互换使用,因为类似的图像分析技术通常可以量化所有三个结构层次。为了简便起见,本综述中使用“fibril”一词来指代filament、fibril和fiber。结构中纤维的排列和相互作用可以显著影响食品的属性,如质地、持水能力和感官感知。因此,控制产生这些纤维的过程以及随后对这些纤维的量化在食品科学中是相关的。在fibril中需要量化的一个重要元素是它们的方向性,因为这可以直接关联到宏观结构的机械响应。讨论的一些方法已成功应用于食品微观结构,对于fibril(filament、fibril和fiber)分析具有不同程度的相关性。

BoneJ是一个Fiji插件,主要用于量化食品产品中的气泡而不是纤维(Domander, Felder, & Doube, 2021)。Guillermic等人(2021年)使用BoneJ测量了通过X射线断层扫描成像的低水分豆类基零食中墙壁和细胞的厚度。同样,Metilli等人(2021年)也使用BoneJ量化了基于可可脂和高油酸葵花籽油的油脂泡沫结构中的气泡,具体包括计数和测量气泡壁,并确定每个气泡的轴线。尽管BoneJ可能不完全适合直接测量食品中的纤维,但其在孔隙率和连通性测量中的应用仍然具有相关性。考虑到文献中记录的气泡与肉类替代品纤维结构之间的关系(Schlangen, Van Der Doef, Van Der Goot, Clausen, & Kodger, 2025),这一点特别有趣。DiameterJ已被用于食品中,以测量通过静电纺丝产生的κ-卡拉胶/酪蛋白纳米纤维的直径(Sharma, Ziegler, & Harte, 2025)。这些基于酪蛋白的纳米结构垫的平均孔隙率和纤维直径是通过SEM图像量化的(图7)。DiameterJ主要应用于纳米纤维,这限制了其在食品科学中的应用范围(Hotaling, Bharti, Kriel, & Simon, 2015)。DiameterJ的局限性在于它难以处理重叠或不规则的纤维,并且需要输入图像具有高对比度。下载:下载高分辨率图像(823KB);下载全尺寸图像。图7:(1) 通过SEM图像量化了静电纺丝κ-卡拉胶/酪蛋白纳米纤维的平均纤维直径。改编自(Sharma et al., 2025),版权属于Elsevier(2024年),经许可使用。(2) 用颜色表示方向的马苏里拉奶酪的µCT扫描。根据Creative Commons许可证(CC BY 4.0)改编自(Pieta, Dahl, Frisvad, Bigdeli, & Christensen, 2025)。(3) 在挤出大豆蛋白浓缩物表皮附近的共聚焦显微图上叠加的矢量图,以及相应的重量有序参数图。根据Creative Commons许可证(CC BY 4.0)改编自(Gobes et al., 2025)。OrientationJ已被用于食品科学中分析主要为肉类替代品的纤维结构。该插件用于量化由豌豆蛋白和燕麦蛋白浓缩物制成的高水分挤出肉类替代品的微观结构(Gaber et al., 2023)。使用蛋白质和脂肪染色的共聚焦激光扫描显微镜(CLSM)图像,作者利用OrientationJ的相干性参数量化了纤维的各向异性,得到的相干性值为0.2-0.5,较高的值对应于更多的蛋白质排列。类似地,OrientationJ插件也被用于分析通过剪切细胞技术产生的带蛋白质染色的肉类替代品的CLSM图像中的方向性(Schlangen et al., 2025),他们发现相干性值范围为0.08至0.56,并发现较高的相干性值与增强的宏观纤维评分之间存在强相关性。他们成功地利用相干性值区分了不同肉类替代品的微观结构排列。此外,OrientationJ还成功应用于培养肉类的菠菜支架(Jones, Rebello, & Gaudette, 2021)。通过细胞核和细胞骨架的方向来测量细胞排列,并使用OrientationJ可视化结构的微丝方向并获取角度分布。同样,该插件的相干性值也成功展示了胶原蛋白结构的差异(Clemons et al., 2018)。OrientationJ插件的优点是易于使用、可访问且快速。OrientationJ中的高斯窗口大小参数允许在特定长度尺度上进行分析,尽管这引入了一定的主观性。一种结构张量尺度空间方法用于量化2D和3D CT扫描中马苏里拉奶酪的各向异性(Pieta et al., 2025)。通过带有颜色到方向映射的µCT扫描可视化显示,马苏里拉奶酪具有不同方向的紧密排列特征(图7)。这种方法似乎适用于其他表现出高各向异性的食品结构,如肉类替代品。然而,这种结构张量方法的一个显著缺点是在处理大量图像时计算量较大。开发了一种旋转傅里叶变换(RFT)方法,用于量化高水分植物蛋白挤出物的各向异性结构(Gobes et al., 2025)。该方法基于使用旋转窗口在傅里叶空间提取角度和各向异性信息。他们将RFT方法应用于CLSM蛋白质染色图像和高场MRI测量的大豆蛋白浓缩物挤出物,显示了挤出物的位置依赖性各向异性(图7)。此外,他们还表明pH值会影响挤出物的微观结构,中性pH值制备的挤出物具有最高的各向异性。RFT方法的优点是可以检测多个重叠的角度特征。5. 一般讨论5.1. 图像预处理目前没有通用的图像处理流程,具体方法取决于成像技术和成像实验的目标(例如,在基于断层扫描的成像中使用非局部均值滤波器来保留纹理)。这通常需要通过试错或依赖以往的研究来选择滤波器和背景校正协议。对底层数学运算及其对后续图像分析步骤的影响有基本了解,可以帮助避免意外的图像处理。虽然没有绝对正确的图像预处理方法,但一致性对于可靠的分割非常重要。更精细的分割算法可能受益于输入图像的变异性,从而不需要过多的图像预处理。5.2. 分割每种分割技术都有其优缺点。传统方法(如阈值处理)计算成本较低且相对容易实现,但对于照明不均匀、有噪声或对象重叠的图像可能处理起来比较困难。虽然机器学习方法在这些问题上有所改进,但仍需要仔细的特征工程和参数调整。深度学习模型提供了最先进的准确性和鲁棒性,但代价是计算需求高且需要大量的、标注良好的数据集。迁移学习是一种常见的策略,可以绕过小型数据集的局限性,使得模型无需大量数据收集即可进行微调。在这种方法中,一个在大型多样数据集(通常来自相关领域)上预训练的模型被适配到更具体的任务。通过使用预训练的模型作为基础,只需要重新训练网络的最后几层或某些部分即可。分割中的另一个挑战是自动化与可解释性之间的平衡。尽管深度学习模型通常性能优越,但有时它们可能被视为“黑箱”,其决策过程不立即透明,因此难以理解。当定量严谨性和可重复性对科学研究至关重要时,这种缺乏可解释性的问题会成为障碍。5.3. 图像量化食品微观结构的量化经常借鉴了其他科学领域(如生物学、医学科学、材料科学和地质科学)的方法。对于大多数样本,这些改编的方法仅在食品结构的形态(液滴、聚集物、网络和纤维)清晰可辨时才适用。通常,在应用适当的图像分析方法之前,需要先对这些可见结构进行初步分类。与其他研究领域相比,食品基质由于通常具有许多组分而具有特定的挑战性,这使得某些分析方法的应用变得复杂。例如,当食品系统中存在纤维时,这些纤维通常不是均匀的。同样,对于液滴、聚集物和网络,它们的分布在样本中的分布可能不够均匀。因此,在应用为均匀系统(例如单细胞菌落)开发的现有方法时,这些方法可能会失败(例如,将分割伪影误认为结构元素)。仅使用来自其他领域的现有工具可能会限制通过图像分析获得的知识,因为这将局限于与分析领域相关的预设值。因此,也需要专门为食品基质设计和优化的技术。许多用于获取不同形态描述符(如形状、大小、分形维数等)的图像分析库也存在,这些库可以用于开发分析方法,例如在MIDAS中用于液滴分析(Saalbrink et al., 2024),其中可以使用‘Simple-ITK’ Python库量化许多其他形状描述符。Brüls-Gill等人(2024)还开发了一个内部Python脚本,用于计算液滴的交叉面积和平均强度。5.4. 体积成像(3D)本综述中提到的大多数分析方法仅限于分析二维(2D)光学平面或投影,然而食品微观结构是三维(3D)的。有可用的3D成像技术,如X射线断层扫描(Guo et al., 2017, Nielsen et al., 2016),用于研究乳液和冰淇淋;中子断层扫描(Martínez-Sanz et al., 2020)用于研究挤出泡沫;或共聚焦成像技术(Ogawa and Matsumura, 2021, Xia et al., 2022)。对于某些技术,由于其特性、施加的剪切或温度敏感性,可使用的样本类型有限。此外,光学技术用于分析和量化3D体积的应用在存在性和数据输出方面都有限,且数据处理通常比2D光学平面的成熟方法要求更高。尽管如此,已经有一些3D图像分析在食品微观结构中的应用。在乳品科学中,Imaris软件(Oxford Instruments, n.d.)已被用于分析共聚焦图像堆栈中的奶酪孔隙度(Ong, Li, Ong, and Gras, 2022; Gras et al., 2024),后来Masiá等人(2024)也用于植物基奶酪的分析。Imaris是一种集成的多模态分析技术,可以将分割和分析结合在一个工作流程中,通常用于细胞生物学、神经科学和癌症研究领域,但随着2D分析的成功适应,这些方法也可能适用于食品科学。其他新兴的开源工具,如Python中的napari(Chiu et al., 2022)和Ilastik中的3D分割(Berg et al., 2019),也可能使3D图像分析更加容易使用。此外,由于食品结构是三维的,因此非常需要验证2D图像分析的发现和相关性是否也适用于3D分析。5.5. 无需分割的量化尽管分割作为常见图像分析工作流程中的有价值步骤越来越受重视,但之前已经讨论过,这一过程可能繁琐,本质上是通过降低数据集的维度来工作的。在此过程中,既有可能引入错误,也可能会排除有价值的数据。为了规避这一风险,也可以直接对原始图像数据进行图像分析,从而包含每个像素的强度值中的信息。例如,分形维数分析是一种广泛用于研究食品微观结构-功能关系的工具,可以量化图像的变异性和异质性(García-Armenta & Gutiérrez-López, 2022)。对于乳胶凝胶,除了分形图像分析(Torres et al., 2012)外,互相关分析是另一种无需分割的工作流程,可以量化蛋白质网络中的相对空间分布(Glover et al., 2019)。然而,在进一步研究这些方法适用于食品基质的适用性及其解释之前,还需要更多关于它们的研究。5.6. 未来展望:微观结构数据在数据驱动的功能预测中的作用本文回顾的图像分析方法实现了食品微观结构的获取、处理和解释。这种新的食品微观结构数据可以与来自其他长度尺度(如X射线散射、流变学等)的结构数据结合,以在更大规模的数据驱动研究中实现多长度尺度的相关性。特别是,使用从机器学习到大型语言模型的各种数据科学方法的新研究可以探索、分析和解释由不同技术生成的复杂数据集中的未见模式。通过这种方式,食品微观结构的量化将显著促进人工智能在技术功能预测中的应用,从而推动合理设计的食品结构的发展。除了结构分析之外,数据科学方法还可以发现食品结构数据与化学、感官和营养科学数据之间的新相关性。这种方法有可能促进数字孪生或计算机模拟框架的发展,从而更高效地开发食品(Hoffman和Cozzolino,2026;Mengucci等人,2022)。6. 结论:先进显微镜的普及以及数据处理能力的提升,尤其是机器学习和深度学习应用的增加,正在推动食品微观结构研究向定量分析方向发展。通常遵循一种通用的图像处理和分析流程:预处理、分割和量化。在预处理步骤中,根据分析目的的不同,可能需要使用不同的强度变换或滤波器来处理复杂的食品基质。同样,也可以采用不同的分割方法,从基于强度的技术到更复杂的机器学习和深度学习算法。虽然后者能够提供更定量的数据,但也应考虑到可能存在的陷阱。这包括简单的基于强度的阈值处理方法与更高级的机器学习/深度学习技术之间的权衡,即它们在准确性方面的差异以及所需的重计算能力。就量化而言,食品微观结构可以被分类为明确的形态(如液滴、细胞、聚集体、网络和纤维)。目前对大多数这类结构的研究方法都是借鉴了其他材料系统的分析技术。因此,现有工作流程和工具的功能可能无法直接应用于量化本质上不同的食品基质。然而,向定量分析的转变正在推动专门为食品微观结构开发的新方法和流程的发展。这将使食品微观结构研究领域能够处理更复杂的数据集,包括3D体积图像和延时序列。这些新的食品微观结构数据维度,结合新的数据科学方法,可以帮助我们更深入地了解食品结构,并获得可用于食品结构合理设计的微观结构原理。

未引用的参考文献(Arganda-Carreras等人,2010;Mishra等人,2025a;Mishra等人,2025b;Oxford Instruments,无日期;Riley等人,2023;Scientific Volume Imaging,无日期;Valente等人,2017;Yang等人,2020;Young和Morrison,2018)

CRediT作者贡献声明:
Ria G. Rodriguez:撰写——初稿编制、数据可视化、研究、概念化。
José C. Bonilla:撰写——审阅与编辑、监督、项目管理、资金争取、概念化。
Mathias P. Clausen:撰写——审阅与编辑、概念化。
Jens Saalbrink:撰写——初稿编制、研究、概念化。
Ina Folmer:撰写——初稿编制、研究、概念化。
Diana Soto-Aguilar:撰写——初稿编制、研究、概念化。
Miek Schlangen:撰写——审阅与编辑、初稿编制、研究、概念化。

生物通微信公众号
微信
新浪微博


生物通 版权所有