CN116912235B - 基于计算机视觉的脱水蒜状态智能检测方法 - Google Patents
基于计算机视觉的脱水蒜状态智能检测方法 Download PDFInfo
- Publication number
- CN116912235B CN116912235B CN202311152406.7A CN202311152406A CN116912235B CN 116912235 B CN116912235 B CN 116912235B CN 202311152406 A CN202311152406 A CN 202311152406A CN 116912235 B CN116912235 B CN 116912235B
- Authority
- CN
- China
- Prior art keywords
- garlic
- abnormal
- dehydrated
- edge
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 235000004611 garlic Nutrition 0.000 title claims abstract description 336
- 238000001514 detection method Methods 0.000 title claims abstract description 39
- 244000245420 ail Species 0.000 title 1
- 240000002234 Allium sativum Species 0.000 claims abstract description 335
- 238000000034 method Methods 0.000 claims abstract description 36
- 230000008569 process Effects 0.000 claims abstract description 27
- 238000011156 evaluation Methods 0.000 claims abstract description 22
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims abstract description 4
- 230000002159 abnormal effect Effects 0.000 claims description 90
- 238000003708 edge detection Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000004891 communication Methods 0.000 claims description 6
- 230000018044 dehydration Effects 0.000 abstract description 21
- 238000006297 dehydration reaction Methods 0.000 abstract description 21
- 238000012545 processing Methods 0.000 abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 17
- 230000000694 effects Effects 0.000 description 12
- 230000008859 change Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000003746 surface roughness Effects 0.000 description 3
- 238000010606 normalization Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000746 body region Anatomy 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 239000002932 luster Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20076—Probabilistic image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30128—Food products
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
Abstract
本发明涉及图像处理技术领域,具体涉及基于计算机视觉的脱水蒜状态智能检测方法,包括:对脱水蒜片进行抽样检测并采集图像,对图像进行预处理得到待处理图像;对脱水过程中的待处理图像进行分析,然后利用函数量化脱水蒜片的颜色参数和形状参数,利用颜色参数和形状参数得到脱水蒜的评估标准;根据脱水蒜的评估标准对脱水蒜片的脱水状态进行分析。本发明通过根据脱水蒜片的脱水状态中的颜色特征和轮廓边缘特征构建模型得到脱水蒜的评估标准,实现了基于计算机视觉的智能化脱水蒜状态检测方法。
Description
技术领域
本发明涉及图像处理技术领域,具体涉及基于计算机视觉的脱水蒜状态智能检测方法。
背景技术
在对脱水蒜片的检测过程中,因为蒜片的数量比较多,人工检测的工作量较大且人工检测过程中对脱水状态的评估会产生一定的主观性,而且同时对大量的脱水蒜片一起进行测量会大大的浪费人力和物力,并且也会存在一定的误差。在利用边缘检测算法检测脱水蒜状态的过程中,检测得到的蒜片边缘往往不能很好的表现蒜片区域的含水量特征,且对于脱水蒜边缘点的处理不能很好地体现蒜片整体轮廓对脱水蒜状态的影响。
发明内容
本发明提供基于计算机视觉的脱水蒜状态智能检测方法,以解决现有的问题。
本发明的基于计算机视觉的脱水蒜状态智能检测方法采用如下技术方案:
本发明一个实施例提供了基于计算机视觉的脱水蒜状态智能检测方法,该方法包括以下步骤:
对脱水蒜片进行抽样拍摄得到每一片脱水蒜片的采集图像,对每一片脱水蒜片的采集图像进行预处理得到待处理图像;
利用待处理图像的灰度直方图得到蒜片整体最大灰度值和蒜片整体最小灰度值,根据蒜片整体最大灰度值和蒜片整体最小灰度值得到蒜片表面的平均亮度和蒜片表面的逆光泽度;根据蒜片表面的平均亮度和蒜片表面的逆光泽度得到蒜片表面的颜色特征;根据边缘检测得到蒜片边缘曲线和蒜片像素点,利用蒜片边缘曲线和蒜片像素点得到蒜片表面的形状参数;利用蒜片表面的颜色特征和蒜片表面的形状参数计算得到脱水蒜片的评估标准;
根据脱水蒜片的评估标准对脱水蒜状态进行检测。
优选的,所述利用待处理图像的灰度直方图得到蒜片整体最大灰度值和蒜片整体最小灰度值,包括的具体步骤如下:
统计待处理图像的灰度直方图的最大值点,并从最大值点开始向灰度直方图右侧进行灰度值的遍历,计算相邻两个灰度值的频率分布的高度差,直至相邻两个灰度值的频率分布的高度差为0,停止遍历,将停止遍历的灰度值记为蒜片整体最大灰度值,同样地从最大值点开始向灰度直方图左侧进行灰度值的遍历得到蒜片整体最小灰度值。
优选的,所述根据蒜片整体最大灰度值和蒜片整体最小灰度值得到蒜片表面的平均亮度和蒜片表面的逆光泽度,包括的具体步骤如下:
获取待处理图像中灰度值在蒜片整体最小灰度值和蒜片整体最大灰度值之间的像素点的数目,记为区间像素点数;统计灰度值大于等于蒜片整体最小灰度值且小于等于蒜片整体最大灰度值的所有像素点的灰度值之和,将灰度值之和除以区间像素点数得到蒜片表面的平均亮度;
获取待处理图像中灰度值大于蒜片整体最大灰度值的像素点的数目,记为超区间像素点数;计算蒜片表面的逆光泽度,具体的计算公式如下:
其中,表示蒜片表面的逆光泽度,/>表示超区间像素点数,/>表示蒜片整体最大灰度值,/>表示灰度值为/>的像素点个数,且有/>,255表示像素点能取到的最大的灰度值。
优选的,所述根据蒜片表面的平均亮度和蒜片表面的逆光泽度得到蒜片表面的颜色特征的具体计算公式如下:
其中,表示蒜片表面的颜色特征,/>表示蒜片表面的平均亮度,/>表示蒜片表面的逆光泽度,/>表示以自然常数为底的指数函数。
优选的,所述根据边缘检测得到蒜片边缘曲线和蒜片像素点,利用蒜片边缘曲线和蒜片像素点得到蒜片表面的形状参数,包括的具体步骤如下:
获取蒜片边缘曲线和待处理图像上所有的异常像素点,计算每两个异常像素点之间的欧式距离,从每个异常像素点开始沿着顺时针方向和与其欧式距离最小的异常像素点进行直线连接,计算直线在待处理图像中的斜率,记为每个异常像素点的异常斜率;计算每个异常像素点的异常斜率差值,统计所有异常斜率差值的数目,将所有异常斜率差值中数值相等的异常斜率差值归为一类异常斜率差值,计算每一类异常斜率差值在所有异常像素点的异常斜率差值中出现的概率,记为每一类异常斜率差值的概率;获取每个蒜片像素点的梯度幅值;计算蒜片表面的形状参数,具体的计算公式如下:
其中,表示脱水蒜片的形状参数,/>表示蒜片像素点的数目,/>表示第/>个蒜片像素点的梯度幅值,且有/>,/>表示第/>类异常斜率差值的概率,/>表示异常斜率差值的数目,且有/>,/>和/>分别表示以自然常数为底的指数函数和以自然常数为底数的对数函数。
优选的,所述蒜片边缘曲线的具体获取过程如下:
利用边缘检测算法对待处理图像进行边缘检测得到若干边缘像素点,对所有的边缘像素点进行闭合曲线拟合得到蒜片边缘曲线。
优选的,所述待处理图像上所有的异常像素点的具体获取过程如下:
计算每个边缘像素点距离蒜片边缘曲线的欧式距离,计算所有欧式距离的方差,记为边缘整体方差;将每个边缘像素点记为检测像素点,计算除去每个检测像素点的其余所有边缘像素点距离蒜片边缘曲线的欧式距离,记为去点欧式距离,计算所有去点欧式距离的方差并记为每个边缘像素点的去点方差;比较每个边缘像素点的去点方差和边缘整体方差的大小,若每个边缘像素点的去点方差小于边缘整体方差,则将每个边缘像素点记为异常像素点,获取待处理图像上所有的异常像素点。
优选的,所述每个异常像素点的异常斜率差值的具体获取过程如下:
将所有异常斜率按照从任意一个异常像素点开始的顺时针方向进行排列得到异常斜率联通序列,序列中最后一个异常斜率的后一个元素为序列中第一个异常斜率,计算异常斜率联通序列中每个异常斜率和与其相邻的后一个异常斜率的差值,记为每个异常像素点的异常斜率差值。
优选的,所述获取每个蒜片像素点的梯度幅值,包括的具体步骤如下:
统计待处理图像上的蒜片边缘曲线内部的所有像素点,将蒜片边缘曲线内部的所有像素点记为蒜片像素点,通过边缘检测算子检测图像中所有蒜片像素点的梯度得到每个蒜片像素点的梯度幅值。
优选的,所述利用蒜片表面的颜色特征和蒜片表面的形状参数计算得到脱水蒜片的评估标准的具体计算公式如下:
其中,表示脱水蒜片的评估标准,/>表示蒜片表面的颜色特征,/>表示脱水蒜片的形状参数。
本发明的技术方案的有益效果是:脱水效果较好的蒜片的颜色偏暗,形状比较褶皱,并且蒜片表面粗糙程度较大,这些参数都是由于脱水蒜片内的含水量的不同导致发生的变化,主要体现在蒜片表面的逆光泽度和蒜片表面的粗糙程度以及蒜片边缘的连续性,本发明通过对这些参数进行建模,得到脱水蒜片的评估标准。本发明通过拟合颜色指标和形状指标使得脱水蒜状态的评估标准解决了传统边缘检测技术无法表征的物理特征量的技术问题,比传统的利用边缘检测对蒜片的脱水状态进行检测得到的结果更加准确。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明基于计算机视觉的脱水蒜状态智能检测方法的步骤流程图。
具体实施方式
为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的基于计算机视觉的脱水蒜状态智能检测方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
下面结合附图具体的说明本发明所提供的基于计算机视觉的脱水蒜状态智能检测方法的具体方案。
请参阅图1,其示出了本发明一个实施例提供的基于计算机视觉的脱水蒜状态智能检测方法的步骤流程图,该方法包括以下步骤:
S001:对脱水蒜片进行抽样拍摄得到每一片脱水蒜片的采集图像,对每一片脱水蒜片的采集图像进行预处理得到待处理图像。
需要说明的是,由于图像中可能会受到其他因素的影响,为了减少后续处理的计算量,对图像进行分割,用现有算法对图像进行阈值分割,目的就是为了待测蒜片与背景分割开,以便于后续的计算。由于在检测过程中,蒜片的数量很多,甚至蒜片都是堆放在一起的,这样的话对检测工作带来一定的挑战,这对整个脱水过程的检测会造成一定的误差,因为有些蒜片可能都拍摄不到,那么对这些蒜片的状态没办法进行检测。因此为了减小检测过程中的误差,为了保证检测结果的精确度,本实施例对每次脱水过程中的蒜片进行固定样本数量检测,对每一片蒜片进行检测,也会减小检测过程中的误差。
具体的,对脱水过程中的蒜片进行固定样本数量检测,其中样本数量为,对每一片脱水蒜片进行拍摄得到每一片脱水蒜片的采集图像,其中拍摄时放置LED灯来对拍摄环境进行一个补光处理;对采集图像进行灰度化处理,得到灰度图像;在灰度图像上,利用算法对灰度图像进行阈值分割,得到脱水蒜片前景图像,将脱水蒜片前景图像记为待处理图像。
至此,得到了待处理图像。
S002:在待处理图像上,根据整体灰度值分布特性进行建模得到蒜片表面的颜色特征。
需要说明的是,不同脱水状态的蒜片会有不同的效果,比如在形态和颜色观感上的差别。有些蒜片的表面比较平滑,说明蒜片内部还含有水分,而脱水效果较好的蒜片表面看起来是褶皱、粗糙的。在灰度图像中,脱水蒜片中不同的含水量会对应不同的灰度值,一般来说,水分含量越高,灰度值就越低;水分含量越低,灰度值就越高。由于预处理之后的图像是一个二值化图像,并且每个蒜片之间都是互相不接触的,所以蒜片的轮廓在图像中都会表现的很清楚;每个蒜片内部的区域的也只会出现两种情况,一种是灰度比较均匀的分布的情况,意味着完全脱水或者没有脱水;另一种情况就是灰度分布出现变化的情况,这代表着蒜片脱水不完全,蒜片内部还存在着部分的水。因为在灰度图像中,脱水蒜片中含水和不含水的部分的区别较大。脱水蒜片的脱水状态中蒜片中的含水量从多到少进行变化,含水的蒜片通常呈现出看起来是一种偏亮的状态,因为光线的反射造成图像中会出现比较高亮的部分;而含水较少的蒜片整体比较偏暗淡,没有那种反光很好的效果,因此其光泽度会下降,故使用逆光泽度和平均亮度对蒜片表面的颜色特征进行衡量。
需要进一步说明的是,在蒜片的灰度直方图中,因为蒜片的在灰度图像中的表现较为单一,没有太多的灰度级,它在直方图中会聚集在一个比较大的灰度区间内,而由于蒜片内水分含量较高且在拍摄过程中光线反射形成了高亮点,直方图中还存在一些零零散散的点,其余的一些零散的像素点就是蒜片的局部亮点,因为不同的含水量在拍摄的过程中会导致蒜片局部出现比较高亮的像素点,在脱水过程中随着含水量的降低,这些局部高亮点就会减少。
具体的,获得待处理图像的灰度直方图,统计待处理图像的灰度直方图的最大值点,并从最大值点开始向灰度直方图右侧进行灰度值的遍历,计算相邻两个灰度值的频率分布的高度差,直至相邻两个灰度值的频率分布的高度差为0,停止遍历,将停止遍历的灰度值记为蒜片整体最大灰度值,同样地从最大值点开始向灰度直方图左侧进行灰度值的遍历得到蒜片整体最小灰度值;获取待处理图像中灰度值在蒜片整体最小灰度值和蒜片整体最大灰度值之间的像素点的数目,记为区间像素点数;计算蒜片表面的平均亮度,具体的计算公式如下:
其中,表示蒜片表面的平均亮度,/>表示区间像素点数,/>和/>分别表示蒜片整体最小灰度值和蒜片整体最大灰度值,/>表示灰度值为/>的像素点个数,且有。/>的值越高表示主体区域的平均灰度越大,蒜片内的含水量越低。
进一步,获取待处理图像中灰度值大于蒜片整体最大灰度值的像素点的数目,记为超区间像素点数;计算蒜片表面的逆光泽度,具体的计算公式如下:
其中,表示蒜片表面的逆光泽度,/>表示超区间像素点数,/>表示蒜片整体最大灰度值,/>表示灰度值为/>的像素点个数,且有/>,255表示像素点能取到的最大的灰度值。/>实际上表征了除蒜片区域之外的局部高亮点的像素点的平均灰度,/>表示将主体灰度之外的灰度值进行归一化,归一化的结果/>越接近于0,则代表高亮点的亮度越大,也可以反映出蒜片含水量越高。
进一步,根据蒜片表面的平均亮度和蒜片表面的逆光泽度,得到蒜片表面的颜色特征的公式如下:
其中,表示蒜片表面的颜色特征,/>表示蒜片表面的平均亮度,/>表示蒜片表面的逆光泽度,/>表示以自然常数为底的指数函数。/>表示将/>反比例归一化至之间,/>值越大则/>在/>之间越小,表征蒜片内的含水量越低,脱水状态越好,当/>越大时,代表此时的蒜片内的含水量越低,脱水状态越好,两者相加得到脱水蒜片表面的颜色特征/>,当/>越大时,代表此时的蒜片内的含水量越低,脱水状态越好。
至此,得到了蒜片表面的颜色特征。
S003:在待处理图像上获取蒜片区域,在蒜片区域内根据蒜片的灰度值梯度特征以及蒜片区域边缘的混乱程度进行建模得到蒜片表面的形状参数,利用蒜片表面的颜色特征和蒜片表面的形状参数计算得到脱水蒜片的评估标准。
需要说明的是,对于脱水蒜片来说,大小均匀,形状完整的脱水蒜片受到更多人的欢迎,在蒜片脱水过程中,含水量越多,蒜片的边缘轮廓可能会更平滑,蒜片表面的也比较平整;脱水效果也好,蒜片的体积也会变小,轮廓也会变的不平滑,当然蒜片表面也会出现一些褶皱,这些情况表示蒜片中的含水量较少或者几乎没有水,因此,体积小、轮廓不平整、表面出现褶皱都可以是衡量脱水蒜片状态的标准。本实施例将脱水蒜在脱水过程中的形状上的变化叫做形状参数,主要包括蒜片表面的粗糙程度和蒜片边缘的平整情况,因为蒜片在脱水的过程中,脱水质量越好,那么蒜片表面就会越粗糙,并且蒜片的边缘会变的不平整。我们就认为脱水蒜片在脱水过程中,蒜片的表面粗糙度越大,并且蒜片边缘越不平整,这个时候的脱水状态较好。
具体的,利用边缘检测算法对待处理图像进行边缘检测得到若干边缘像素点,对所有的边缘像素点利用最小二乘法进行闭合曲线拟合得到蒜片边缘曲线,计算每个边缘像素点距离拟合曲线的欧式距离,计算所有欧式距离的方差,记为边缘整体方差;将每个边缘像素点记为检测像素点,计算除去每个检测像素点的其余所有边缘像素点距离拟合曲线的欧式距离,记为去点欧式距离,计算所有去点欧式距离的方差并记为每个边缘像素点的去点方差;比较每个边缘像素点的去点方差和边缘整体方差的大小,若每个边缘像素点的去点方差小于边缘整体方差,则将每个边缘像素点记为异常像素点;获取待处理图像上所有的异常像素点,计算所有每两个异常像素点之间的欧式距离,从某个异常像素点开始沿着顺时针方向和与其欧式距离最小的异常像素点进行直线连接,计算直线在待处理图像中的斜率,记为每个异常像素点的异常斜率;将所有异常斜率按照从任意一个异常像素点开始的顺时针方向进行排列得到异常斜率联通序列,序列中最后一个异常斜率的后一个元素为序列中第一个异常斜率,计算异常斜率联通序列中每个异常斜率和与其相邻的后一个异常斜率的差值,记为每个异常像素点的异常斜率差值;统计所有异常斜率差值的数目,将所有异常斜率差值中数值相等的异常斜率差值归为一类异常斜率差值,计算每一类异常斜率差值在所有异常像素点的异常斜率差值中出现的概率,记为每一类异常斜率差值的概率;统计待处理图像上的蒜片边缘曲线内部的所有像素点,将蒜片边缘曲线内部的所有像素点记为蒜片像素点,通过/>算子检测图像中所有蒜片像素点的梯度得到每个蒜片像素点的梯度幅值;计算蒜片表面的形状参数,具体的计算公式如下:
其中,表示脱水蒜片的形状参数,/>表示蒜片像素点的数目,/>表示第/>个蒜片像素点的梯度幅值,且有/>,/>表示第/>类异常斜率差值的概率,/>表示异常斜率差值的数目,且有/>,/>和/>分别表示以自然常数为底的指数函数和以自然常数为底数的对数函数。/>表示蒜片图像内的像素点的灰度梯度的密度,密度越大就表示蒜片表面的纹理越多,即蒜片表面的粗糙度越大且脱水蒜片的脱水状态越好越好,而/>表示脱水蒜片边缘的异常像素点的斜率差值的信息熵,信息熵越大,就代表斜率差值的混乱程度越大,即边缘处变得凹凸不平,这些原因都是由于脱水蒜片在脱水过程中导致的,故形状参数/>越大,就表示脱水蒜片的脱水状态越好。
需要进一步说明的是,由于脱水蒜片在脱水的过程中,随着含水量的减少,蒜片表面会出现一些褶皱的情况,并且蒜片的边缘也会变的凹凸不平;而且,脱水过程中蒜片表面的光泽也会变得较暗,这是因为在拍摄过程中由于含水量的降低,导致光线的反射效果较差,才会使得蒜片表面变得比较暗淡,也就是脱水效果达到较好的一个标准,因此通过对蒜片的观感和形状方面的特征关系来对蒜片的脱水效果进行一个量化处理。针对脱水蒜片的脱水状态从颜色和形状两个方面对脱水状态进行评估,脱水蒜片的颜色是随着脱水蒜片中的含水量变化进行变化,含水量多的蒜片在图像中会显得比较高亮,也就是蒜片表面的光泽度较大,此时蒜片的形状参数较小,即蒜片的表面比较光滑并且蒜片的边缘比较平整;相反的含水量小的蒜片,整体形状会变得比较褶皱,并且表面的粗糙度较大,也就是蒜片内的像素点的梯度整体变化较大,而它的颜色也相对变得较暗,即蒜片表面的光泽度较差,此时证明算片内含水量少,说明脱水效果更佳。
进一步,计算脱水蒜片的评估标准,具体的计算公式如下:
其中,表示脱水蒜片的评估标准,/>表示蒜片表面的颜色特征,/>表示脱水蒜片的形状参数。颜色特征/>越大时表示蒜片内的含水量越小,随着含水量的减少,蒜片的形状参数将会变得越大,形状参数/>越大表示脱水蒜片的脱水状态越好,因此用颜色特征和形状参数的欧式范数作为脱水蒜片脱水状态的评估标准,/>越大表示脱水蒜片的脱水效果越好。
至此,得到了脱水蒜片的脱水状态的评估标准。
S004:根据脱水蒜片的评估标准对脱水蒜状态进行判断并标记得到检测报告,完成基于计算机视觉的脱水蒜状态的智能检测。
预设一个脱水阈值,其中本实施例以/>为例进行叙述,本实施例不进行具体限定,其中脱水阈值/>可根据具体实施情况而定。
具体的,根据每一幅待处理图像计算每一片脱水蒜片的评估标准,将每一片脱水蒜片的评估标准与脱水阈值进行比较,若某一片脱水蒜片的评估标准大于脱水阈值,则将该脱水蒜片的脱水状态标记为不良,否则将脱水蒜片的脱水状态标记为良好,将所有个样本的脱水蒜片的脱水状态打印成检测报告。
至此,完成了基于计算机视觉的脱水蒜状态的智能检测。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.基于计算机视觉的脱水蒜状态智能检测方法,其特征在于,该方法包括以下步骤:
对脱水蒜片进行抽样拍摄得到每一片脱水蒜片的采集图像,对每一片脱水蒜片的采集图像进行预处理得到待处理图像;
利用待处理图像的灰度直方图得到蒜片整体最大灰度值和蒜片整体最小灰度值,根据蒜片整体最大灰度值和蒜片整体最小灰度值得到蒜片表面的平均亮度和蒜片表面的逆光泽度;根据蒜片表面的平均亮度和蒜片表面的逆光泽度得到蒜片表面的颜色特征;根据边缘检测得到蒜片边缘曲线和蒜片像素点,利用蒜片边缘曲线和蒜片像素点得到蒜片表面的形状参数;利用蒜片表面的颜色特征和蒜片表面的形状参数计算得到脱水蒜片的评估标准;
根据脱水蒜片的评估标准对脱水蒜状态进行检测;
所述利用待处理图像的灰度直方图得到蒜片整体最大灰度值和蒜片整体最小灰度值,包括的具体步骤如下:
统计待处理图像的灰度直方图的最大值点,并从最大值点开始向灰度直方图右侧进行灰度值的遍历,计算相邻两个灰度值的频率分布的高度差,直至相邻两个灰度值的频率分布的高度差为0,停止遍历,将停止遍历的灰度值记为蒜片整体最大灰度值,同样地从最大值点开始向灰度直方图左侧进行灰度值的遍历得到蒜片整体最小灰度值;
所述根据蒜片整体最大灰度值和蒜片整体最小灰度值得到蒜片表面的平均亮度和蒜片表面的逆光泽度,包括的具体步骤如下:
获取待处理图像中灰度值在蒜片整体最小灰度值和蒜片整体最大灰度值之间的像素点的数目,记为区间像素点数;统计灰度值大于等于蒜片整体最小灰度值且小于等于蒜片整体最大灰度值的所有像素点的灰度值之和,将灰度值之和除以区间像素点数得到蒜片表面的平均亮度;
获取待处理图像中灰度值大于蒜片整体最大灰度值的像素点的数目,记为超区间像素点数;计算蒜片表面的逆光泽度,具体的计算公式如下:
其中,表示蒜片表面的逆光泽度,/>表示超区间像素点数,/>表示蒜片整体最大灰度值,/>表示灰度值为/>的像素点个数,且有/>,255表示像素点能取到的最大的灰度值。
2.根据权利要求1所述基于计算机视觉的脱水蒜状态智能检测方法,其特征在于,所述根据蒜片表面的平均亮度和蒜片表面的逆光泽度得到蒜片表面的颜色特征的具体计算公式如下:
其中,表示蒜片表面的颜色特征,/>表示蒜片表面的平均亮度,/>表示蒜片表面的逆光泽度,/>表示以自然常数为底的指数函数。
3.根据权利要求1所述基于计算机视觉的脱水蒜状态智能检测方法,其特征在于,所述根据边缘检测得到蒜片边缘曲线和蒜片像素点,利用蒜片边缘曲线和蒜片像素点得到蒜片表面的形状参数,包括的具体步骤如下:
获取蒜片边缘曲线和待处理图像上所有的异常像素点,计算每两个异常像素点之间的欧式距离,从每个异常像素点开始沿着顺时针方向和与其欧式距离最小的异常像素点进行直线连接,计算直线在待处理图像中的斜率,记为每个异常像素点的异常斜率;计算每个异常像素点的异常斜率差值,统计所有异常斜率差值的数目,将所有异常斜率差值中数值相等的异常斜率差值归为一类异常斜率差值,计算每一类异常斜率差值在所有异常像素点的异常斜率差值中出现的概率,记为每一类异常斜率差值的概率;获取每个蒜片像素点的梯度幅值;计算蒜片表面的形状参数,具体的计算公式如下:
其中,表示脱水蒜片的形状参数,/>表示蒜片像素点的数目,/>表示第/>个蒜片像素点的梯度幅值,且有/>,/>表示第/>类异常斜率差值的概率,/>表示异常斜率差值的数目,且有/>,/>和/>分别表示以自然常数为底的指数函数和以自然常数为底数的对数函数。
4.根据权利要求3所述基于计算机视觉的脱水蒜状态智能检测方法,其特征在于,所述蒜片边缘曲线的具体获取过程如下:
利用边缘检测算法对待处理图像进行边缘检测得到若干边缘像素点,对所有的边缘像素点进行闭合曲线拟合得到蒜片边缘曲线。
5.根据权利要求3所述基于计算机视觉的脱水蒜状态智能检测方法,其特征在于,所述待处理图像上所有的异常像素点的具体获取过程如下:
计算每个边缘像素点距离蒜片边缘曲线的欧式距离,计算所有欧式距离的方差,记为边缘整体方差;将每个边缘像素点记为检测像素点,计算除去每个检测像素点的其余所有边缘像素点距离蒜片边缘曲线的欧式距离,记为去点欧式距离,计算所有去点欧式距离的方差并记为每个边缘像素点的去点方差;比较每个边缘像素点的去点方差和边缘整体方差的大小,若每个边缘像素点的去点方差小于边缘整体方差,则将每个边缘像素点记为异常像素点,获取待处理图像上所有的异常像素点。
6.根据权利要求3所述基于计算机视觉的脱水蒜状态智能检测方法,其特征在于,所述每个异常像素点的异常斜率差值的具体获取过程如下:
将所有异常斜率按照从任意一个异常像素点开始的顺时针方向进行排列得到异常斜率联通序列,序列中最后一个异常斜率的后一个元素为序列中第一个异常斜率,计算异常斜率联通序列中每个异常斜率和与其相邻的后一个异常斜率的差值,记为每个异常像素点的异常斜率差值。
7.根据权利要求3所述基于计算机视觉的脱水蒜状态智能检测方法,其特征在于,所述获取每个蒜片像素点的梯度幅值,包括的具体步骤如下:
统计待处理图像上的蒜片边缘曲线内部的所有像素点,将蒜片边缘曲线内部的所有像素点记为蒜片像素点,通过边缘检测算子检测图像中所有蒜片像素点的梯度得到每个蒜片像素点的梯度幅值。
8.根据权利要求1所述基于计算机视觉的脱水蒜状态智能检测方法,其特征在于,所述利用蒜片表面的颜色特征和蒜片表面的形状参数计算得到脱水蒜片的评估标准的具体计算公式如下:
其中,表示脱水蒜片的评估标准,/>表示蒜片表面的颜色特征,/>表示脱水蒜片的形状参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311152406.7A CN116912235B (zh) | 2023-09-08 | 2023-09-08 | 基于计算机视觉的脱水蒜状态智能检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311152406.7A CN116912235B (zh) | 2023-09-08 | 2023-09-08 | 基于计算机视觉的脱水蒜状态智能检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116912235A CN116912235A (zh) | 2023-10-20 |
CN116912235B true CN116912235B (zh) | 2023-11-24 |
Family
ID=88356783
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311152406.7A Active CN116912235B (zh) | 2023-09-08 | 2023-09-08 | 基于计算机视觉的脱水蒜状态智能检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116912235B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118247264B (zh) * | 2024-05-20 | 2024-09-17 | 华域生物科技(天津)有限公司 | 一种受精卵胚胎解冻形态的图像检测评估方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492668A (zh) * | 2018-10-10 | 2019-03-19 | 华中科技大学 | 基于多通道卷积神经网络的mri不同期多模图像表征方法 |
CN113128249A (zh) * | 2019-12-26 | 2021-07-16 | 合肥美亚光电技术股份有限公司 | 大蒜的疤识别方法、系统及分选设备 |
CN114972187A (zh) * | 2022-04-20 | 2022-08-30 | 北京与子成科技有限公司 | 基于人工智能的裂纹缺陷评估方法 |
CN115082470A (zh) * | 2022-08-22 | 2022-09-20 | 启东市固德防水布有限公司 | 一种基于图像识别的防水布性能检测方法及系统 |
CN217609437U (zh) * | 2022-08-11 | 2022-10-21 | 山东凯益食品有限公司 | 一种脱水蒜制品异物剔除设备 |
CN115880282A (zh) * | 2023-01-16 | 2023-03-31 | 四川远方云天食品科技有限公司 | 一种用于火锅底料均匀装料的监测方法及其监测装置 |
CN116645363A (zh) * | 2023-07-17 | 2023-08-25 | 山东富鹏生物科技有限公司 | 一种基于视觉的淀粉生产质量实时检测方法 |
CN116645429A (zh) * | 2023-07-25 | 2023-08-25 | 山东中胜涂料有限公司 | 一种视觉辅助的涂料生产样品色彩分析检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109671046B (zh) * | 2017-10-12 | 2021-08-17 | 精诚工坊电子集成技术(北京)有限公司 | 利用皮肤图像分析皮肤水分的方法及装置 |
-
2023
- 2023-09-08 CN CN202311152406.7A patent/CN116912235B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492668A (zh) * | 2018-10-10 | 2019-03-19 | 华中科技大学 | 基于多通道卷积神经网络的mri不同期多模图像表征方法 |
CN113128249A (zh) * | 2019-12-26 | 2021-07-16 | 合肥美亚光电技术股份有限公司 | 大蒜的疤识别方法、系统及分选设备 |
CN114972187A (zh) * | 2022-04-20 | 2022-08-30 | 北京与子成科技有限公司 | 基于人工智能的裂纹缺陷评估方法 |
CN217609437U (zh) * | 2022-08-11 | 2022-10-21 | 山东凯益食品有限公司 | 一种脱水蒜制品异物剔除设备 |
CN115082470A (zh) * | 2022-08-22 | 2022-09-20 | 启东市固德防水布有限公司 | 一种基于图像识别的防水布性能检测方法及系统 |
CN115880282A (zh) * | 2023-01-16 | 2023-03-31 | 四川远方云天食品科技有限公司 | 一种用于火锅底料均匀装料的监测方法及其监测装置 |
CN116645363A (zh) * | 2023-07-17 | 2023-08-25 | 山东富鹏生物科技有限公司 | 一种基于视觉的淀粉生产质量实时检测方法 |
CN116645429A (zh) * | 2023-07-25 | 2023-08-25 | 山东中胜涂料有限公司 | 一种视觉辅助的涂料生产样品色彩分析检测方法 |
Non-Patent Citations (3)
Title |
---|
Applications of Computer Vision for Assessing Quality of Agri-food Products: A Review of Recent Research Advances;Ji Ma 等;《Critical Reviews in Food Science and Nutrition》;第56卷;113-127 * |
基于计算机视觉的脱水蒜片检测与分级研究;汪源;《中国优秀硕士学位论文全文数据库 信息科技辑》;第2012年卷(第6期);I138-1811 * |
大蒜干燥技术研究进展;郝启栋 等;《中国调味品》;第45期(第1期);193-195 * |
Also Published As
Publication number | Publication date |
---|---|
CN116912235A (zh) | 2023-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115375676B (zh) | 一种基于图像识别的不锈钢制品质量检测方法 | |
CN114723704B (zh) | 一种基于图像处理的纺织品质量评估方法 | |
CN116912235B (zh) | 基于计算机视觉的脱水蒜状态智能检测方法 | |
CN116168026A (zh) | 基于计算机视觉的水质检测方法及系统 | |
CN103034838B (zh) | 一种基于图像特征的特种车辆仪表类型识别与标定方法 | |
CN114170208A (zh) | 一种基于人工智能的纸制品缺陷检测方法 | |
CN116758077B (zh) | 一种冲浪板表面平整度在线检测方法及系统 | |
CN114494210A (zh) | 基于图像处理的塑料膜生产缺陷检测方法及系统 | |
CN117011292B (zh) | 一种复合板表面质量快速检测方法 | |
CN116385450B (zh) | 基于图像处理的ps片材抗耐磨性检测方法 | |
CN117422712B (zh) | 基于图像滤波处理的塑料母粒视觉检测方法及系统 | |
CN111160451A (zh) | 一种柔性材料检测方法及其储存介质 | |
CN116246174B (zh) | 基于图像处理的甘薯种类识别方法 | |
CN117541582B (zh) | 一种高频变频器igbt绝缘质量检测方法 | |
CN116664565A (zh) | 一种光伏太阳能电池片的隐裂检测方法及系统 | |
CN116645363B (zh) | 一种基于视觉的淀粉生产质量实时检测方法 | |
CN115359043A (zh) | 一种pcb板面异物智能检测方法 | |
CN116309599B (zh) | 基于污水前置处理的水质视觉监测方法 | |
CN118115498B (zh) | 用于钛棒的光泽度快速检测方法及系统 | |
CN118279336B (zh) | 一种用于航空发动机叶片尺寸的智能化测量方法 | |
CN105787912B (zh) | 一种基于分类的阶跃型边缘亚像素定位方法 | |
CN109961432A (zh) | 一种滤布破损的检测方法和系统 | |
CN115511907B (zh) | 一种用于led屏幕的划痕检测方法 | |
CN118134938A (zh) | 一种耐火砖表面缺陷视觉检测方法及系统 | |
CN111738931A (zh) | 光伏阵列无人机航拍图像的阴影去除算法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |