CN115082418B - 一种汽车零部件精密识别方法 - Google Patents

一种汽车零部件精密识别方法 Download PDF

Info

Publication number
CN115082418B
CN115082418B CN202210822912.1A CN202210822912A CN115082418B CN 115082418 B CN115082418 B CN 115082418B CN 202210822912 A CN202210822912 A CN 202210822912A CN 115082418 B CN115082418 B CN 115082418B
Authority
CN
China
Prior art keywords
defect
window
run
gray level
frequency
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
Application number
CN202210822912.1A
Other languages
English (en)
Other versions
CN115082418A (zh
Inventor
李振鹏
张延飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Liaocheng Fufeng Auto Parts Co ltd
Original Assignee
Shandong Liaocheng Fufeng Auto Parts Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong Liaocheng Fufeng Auto Parts Co ltd filed Critical Shandong Liaocheng Fufeng Auto Parts Co ltd
Priority to CN202210822912.1A priority Critical patent/CN115082418B/zh
Publication of CN115082418A publication Critical patent/CN115082418A/zh
Application granted granted Critical
Publication of CN115082418B publication Critical patent/CN115082418B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20076Probabilistic image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30164Workpiece; Machine component

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)

Abstract

本发明涉及图像处理领域,具体涉及一种汽车零部件精密识别方法。获取齿轮灰度图像进行窗口分割;对每个窗口进行边缘检测,获取缺陷窗口;计算缺陷窗口中边缘像素点的缺陷程度,获取缺陷像素点;以缺陷像素点为中心,根据邻域上各个方向的缺陷像素点个数计算方向概率;获取缺陷窗口各个方向的灰度游程矩阵,获取每个缺陷窗口的最终灰度游程矩阵;将最终灰度游程矩阵划分为长游程区域和短游程区域;计算游程频数差异度;将最终灰度游程矩阵划分为高灰度阶区域和低灰度阶区域;计算灰度频数差异度;计算每个缺陷窗口的置信度判断该窗口的缺陷类别。本发明结合矩阵中像素的灰度频数差异特征、游程长度差异特征,能够更加准确的识别缺陷类型。

Description

一种汽车零部件精密识别方法
技术领域
本发明涉及图像处理领域,具体涉及一种汽车零部件精密识别方法。
背景技术
齿轮是组成汽车的一个很重要的零部件,但其在生产加工过程中由于锻造力度、温度等原因造成生产缺陷,包括齿轮表面裂纹或划痕等,一旦汽车采用了具有裂纹的残次品齿轮,在重载和应力过大情况下极有可能导致齿轮的破坏,甚至导致整个机械设备的故障和报废,严重时还会危及人身安全,而齿轮表面是划痕缺陷时,可能是仅由摩擦产生的表面缺陷,若误判为裂纹缺陷,则需对其进行更换,对车主造成不必要的经济损失,因此需要齿轮的缺陷类型进行准确判断。
目前对于齿轮表面缺陷主要依据人工检测,人工检测效率低下且不容易准确区分出裂纹与划痕,容易造成漏检和误检,主观性太强,且当前大多是对齿轮使用后的磨损程度进行分析,对于齿轮的表面生产质量分析较少,同时常规的图像处理技术无法对齿轮的缺陷类型做出准确的区分,存在一定的误差。
由此本发明提出一种齿轮表面缺陷的精密识别方法,结合齿轮图像中缺陷区域的特征构建游程矩阵,从而根据游程矩阵中各个区域的像素点特征识别齿轮缺陷是裂纹还是划痕。
发明内容
本发明提供一种汽车零部件精密识别方法,以解决现有的问题,包括:获取齿轮灰度图像进行窗口分割;对每个窗口进行边缘检测,获取缺陷窗口;计算缺陷窗口中边缘像素点的缺陷程度,获取缺陷像素点;以缺陷像素点为中心,根据邻域上各个方向的缺陷像素点个数计算方向概率;获取缺陷窗口各个方向的灰度游程矩阵,获取每个缺陷窗口的最终灰度游程矩阵;将最终灰度游程矩阵划分为长游程区域和短游程区域;计算游程频数差异度;将最终灰度游程矩阵划分为高灰度阶区域和低灰度阶区域;计算灰度频数差异度;计算每个缺陷窗口的置信度判断该窗口的缺陷类别。
本发明首先通过边缘检测对齿轮图像中的像素点进行初步判断,从而对获得的疑似缺陷窗口进行灰度游程矩阵的建立,大大减少了计算量,同时在建立灰度游程矩阵时通过对窗口中缺陷像素点在各个方向的分布对窗口各个方向的游程矩阵进行加权,得到最终的灰度游程矩阵能够更加完整的表征缺陷窗口中像素点的特征,进而结合矩阵中像素的灰度频数差异特征、游程长度差异特征,双重特征紧密结合纹理信息,能够更加准确的识别缺陷类型。
本发明采用如下技术方案,一种汽车零件精密识别方法,包括:
获取齿轮灰度图像,对灰度图像进行窗口分割,得到多个窗口。
对每个窗口进行边缘检测,将检测得到的边缘像素点的个数占比大于阈值的窗口作为缺陷窗口。
计算每个缺陷窗口中的边缘像素点的缺陷程度,将缺陷程度大于阈值的边缘像素点作为缺陷像素点,获取每个缺陷窗口中的所有缺陷像素点。
以每个缺陷像素点为中心,统计每个缺陷像素点邻域上各个方向的缺陷像素点个数,根据每个方向上缺陷像素点的个数计算每个方向的方向概率。
获取每个缺陷窗口各个方向的灰度游程矩阵,根据每个窗口各个方向的灰度游程矩阵以及对应方向的方向概率获取每个缺陷窗口的最终灰度游程矩阵。
将每个缺陷窗口的最终灰度游程矩阵平均划分为长游程区域和短游程区域。
根据最终灰度游程矩阵中长游程区域和短游程区域中元素的频数计算每个缺陷窗口的游程频数差异度。
将每个缺陷窗口的最终灰度游程矩阵平均划分为高灰度阶区域和低灰度阶区域。
根据最终灰度游程矩阵中高灰度阶区域和低灰度阶区域中元素的频数计算每个缺陷窗口的灰度频数差异度。
根据每个缺陷窗口的游程频数差异度和灰度频数差异度计算每个缺陷窗口的置信度,根据每个缺陷窗口的置信度判断该窗口的缺陷类别。
进一步的,一种汽车零部件精密识别方法,计算每个缺陷窗口中的边缘像素点的缺陷程度的方法为:
获取每个边缘像素点的邻域内与该边缘像素点灰度值相同的像素点个数;
对每个边缘像素点的邻域中像素点依次进行遍历,获取该边缘像素点灰度值相同的连续像素点个数,作为每个边缘像素点的像素聚集程度;
根据每个边缘像素点的像素聚集程度以及邻域内与该边缘像素点灰度值相同的像素点个数计算每个边缘像素点的缺陷程度。
进一步的,一种汽车零部件精密识别方法,获取每个缺陷窗口的最终灰度游程矩阵的方法为:
根据每个缺陷窗口中0°、45°、90°、135°方向上的像素点数量获取对应的方向概率;
获取每个缺陷窗口在0°、45°、90°、135°方向上的灰度游程矩阵,根据每个方向上的灰度游程矩阵及对应方向的方向概率进行加权求和,得到每个缺陷窗口的最终灰度游程矩阵。
进一步的,一种汽车零部件精密识别方法,计算每个缺陷窗口的游程频数差异度的方法为:
根据游程长度将最终灰度游程矩阵平均分为长游程区域和短游程区域,分别统计长游程区域和短游程区域中所有元素的频数;
根据长游程区域中所有元素的频数和短游程区域中所有元素的频数计算每个缺陷窗口的游程频数差异度,表达式为:
Figure 100002_DEST_PATH_IMAGE001
其中,
Figure 206369DEST_PATH_IMAGE002
表示每个缺陷窗口的游程频数差异度,
Figure 581986DEST_PATH_IMAGE003
表示短游程的频数,
Figure 487756DEST_PATH_IMAGE004
表示长游程的频数。
进一步的,一种汽车零部件精密识别方法,计算每个缺陷窗口的灰度频数差异度的方法为:
根据灰度阶将最终灰度游程矩阵平均划分为高灰度阶区域和低灰度阶区域,分别统计高灰度阶区域和低灰度阶区域中所有元素的频数;
根据高灰度阶区域和低灰度阶区域中所有元素的频数计算每个缺陷窗口的灰度频数差异度,表达式为:
Figure 632430DEST_PATH_IMAGE005
其中,
Figure DEST_PATH_IMAGE006
表示每个缺陷窗口的灰度频数差异度,
Figure 401672DEST_PATH_IMAGE007
表示游程长度为j,灰度阶为i的像素频数,m表示最终灰度游程矩阵的游程长度,
Figure 592482DEST_PATH_IMAGE008
表示最终灰度游程矩阵的灰度阶长度,
Figure 20052DEST_PATH_IMAGE009
表示最终灰度游程矩阵的第一阶灰度阶。
进一步的,一种汽车零部件精密识别方法,计算每个缺陷窗口的置信度的方法为:
根据每个缺陷窗口的游程频数差异度和灰度频数差异度计算每个缺陷窗口的置信度,表达式为:
Figure 347128DEST_PATH_IMAGE010
其中,
Figure 834741DEST_PATH_IMAGE011
表示每个缺陷窗口的置信度,
Figure 512847DEST_PATH_IMAGE002
表示每个缺陷窗口的游程频数差异度,
Figure 503630DEST_PATH_IMAGE006
表示每个缺陷窗口的灰度频数差异度。
进一步的,一种汽车零部件精密识别方法,根据每个缺陷窗口的置信度判断该窗口的缺陷类别的方法为:
当缺陷窗口的置信度大于阈值时,该缺陷窗口中存在裂纹缺陷;
当缺陷窗口的置信度小于阈值时,该缺陷窗口中存在划痕。
本发明的有益效果是:本发明首先通过边缘检测对齿轮图像中的像素点进行初步判断,从而对获得的疑似缺陷窗口进行灰度游程矩阵的建立,大大减少了计算量,同时在建立灰度游程矩阵时通过对窗口中缺陷像素点在各个方向的分布对窗口各个方向的游程矩阵进行加权,得到最终的灰度游程矩阵能够更加完整的表征缺陷窗口中像素点的特征,进而结合矩阵中像素的灰度频数差异特征、游程长度差异特征,双重特征紧密结合纹理信息,能够更加准确的识别缺陷类型。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的一种汽车零部件精密识别方法结构示意图;
图2为本发明实施例将每个缺陷窗口的最终灰度游程矩阵沿竖直方向划分示意图;
图3为本发明实施例将每个缺陷窗口的最终灰度游程矩阵沿水平方向划分示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,给出了本发明实施例的一种汽车零部件精密识别方法结构示意图,包括:
101.获取齿轮灰度图像,对灰度图像进行窗口分割,得到多个窗口。
对每个窗口进行边缘检测,将检测得到的边缘像素点的个数占比大于阈值的窗口作为缺陷窗口。
本发明所针对的具体场景为:齿轮在生产加工过程中产生的裂纹或划痕缺陷,人眼难以准确高效区分,因此使用图像处理技术分析缺陷图像,做到计算机对缺陷类型进行识别。
齿轮横截面朝上水平摆放,固定光源于齿轮上方,采集过程中使用工业专业采集相机,位于齿轮正上方开始采集,由于采集到的图像受到噪声干扰,因此对采集到的图像进行高斯滤波降噪处理,将采集到的图像进行语义分割去除背景区域,保留齿轮区域,并将齿轮区域图像转换为灰度图像。
由于齿轮横截面可能存在多处缺陷,我们需要对每一处缺陷结合该缺陷的灰度游程矩阵分析缺陷的纹理特征,因此我们需要分割图像,尽可能的将每一处缺陷分割到一个窗口,进而筛选存在缺陷的窗口,根据窗口灰度共生矩阵内反应的纹理信息差异来分析缺陷的裂纹类型。
正常无缺陷区域灰度均匀无异常,而缺陷区域的像素灰度值与正常区域的像素灰度值存在差异,主要表现为裂纹区域的像素灰度值低于正常区域像素灰度值,划痕区域像素灰度值高于正常区域像素灰度值,因此没有必要获取整个齿轮的灰度共生矩阵,只需要获得存在缺陷的窗口区域即可。
对每一窗口都使用sobel算子进行边缘检测,当滑块中检测出较多的边缘时,认为该滑块存在缺陷,反之认为当前滑块无缺陷。
设定当检测出的边缘上的像素数量占窗口像素总数量的大于阈值
Figure DEST_PATH_IMAGE012
时认为需要对缺陷进行后续识别处理,
Figure 639207DEST_PATH_IMAGE013
>0.03,n为边缘检测线上的像素数量,N为窗口的像素数量。
102.计算每个缺陷窗口中的边缘像素点的缺陷程度,将置信度大于阈值的边缘像素点作为缺陷像素点,获取每个缺陷窗口中的所有缺陷像素点。
由于裂纹的是沿着晶粒裂开的,所以统计裂纹灰度游程矩阵时的统计方向应该更加精细的为每个方向划分不同的权重,而划痕区域的灰度共生矩阵也没有必要朝各个方向都进行计算再加权求和。综上所述,我们通过分析裂纹与划痕的实时图像纹理特征获得自适应的灰度共生矩阵统计方向,这样能够使得后续的灰度游程矩阵的频数信息更加真实。
由于本发明灰度级量化是根据常规的分段量化,所以可能存在将非缺陷区域的像素点与缺陷区域像素点的值归并为同一类,且伪缺陷点恰好出现在缺陷点的八邻域内,故我们需要在八邻域内对灰度值相同的点进行筛选,伪缺陷点一般为孤立零散的,周围没有连续的伪缺陷点,真实缺陷点一般邻域一侧或整个邻域大部分区域都是缺陷点。
首先筛选出当前缺陷点A邻域内灰度值与点A相同的点:
Figure 32143DEST_PATH_IMAGE014
,式中H为像素点对应的灰度值,将缺陷点A邻域内不同灰度值的点标记为0,相同灰度值的标记为1。
分析A点邻域内像素的分布特征,即统计邻域同灰度像素数目以及这些像素的分布聚集性,从而计算点A为缺陷点的缺陷程度为:
Figure 197545DEST_PATH_IMAGE015
,式中
Figure 216185DEST_PATH_IMAGE016
为点A邻域内同灰度值的像素数量,HP为像素点
Figure 783433DEST_PATH_IMAGE017
邻域内同灰度值像素的聚集程度,计算公式为:
Figure 347269DEST_PATH_IMAGE018
,HP为像素聚集程度,初始值为0,从像素邻域的第一个像素值为0的点展开开始统计,当下一个像素灰度值与前一个像素灰度值相同且都不为0时,中心像素的聚集程度加1,当出现下一个像素与前一像素灰度值不同时聚集程度保持不变,最终获得当前点为真缺陷点的缺陷程度
Figure 999967DEST_PATH_IMAGE019
,本发明中当像素点的缺陷程度
Figure 573031DEST_PATH_IMAGE020
时,认为该像素点为缺陷像素点。
计算每个缺陷窗口中的边缘像素点的缺陷程度的方法为:
获取每个边缘像素点的邻域内与该边缘像素点灰度值相同的像素点个数;
对每个边缘像素点的邻域中像素点依次进行遍历,获取该边缘像素点灰度值相同的连续像素点个数,作为每个边缘像素点的像素聚集程度;
根据每个边缘像素点的像素聚集程度以及邻域内与该边缘像素点灰度值相同的像素点个数计算每个边缘像素点的缺陷程度。
103.以每个缺陷像素点为中心,统计每个缺陷像素点邻域上各个方向的缺陷像素点个数,根据每个方向上缺陷像素点的个数计算每个方向的方向概率。
以缺陷像素点A为中心,获取缺陷像素点邻域内所有的缺陷像素点,计算这些真缺陷点在点A邻域内的分布角度:
Figure 729206DEST_PATH_IMAGE021
,式中
Figure 729523DEST_PATH_IMAGE022
为邻域同灰度值真缺陷点像素的横纵坐标值,x,y为中心像素点A的坐标。
统计所有真缺陷点邻域内同灰度值像素相对于中心像素点的方向角度:
Figure 603938DEST_PATH_IMAGE023
,对这些方向角度做依次概率分布统计,分别计算出0°方向概率为:
Figure 728496DEST_PATH_IMAGE024
,其中
Figure 4756DEST_PATH_IMAGE025
,式中
Figure 175974DEST_PATH_IMAGE026
为统计0度和180度角的像素数量,i为真缺陷像素点总数量,同理计算45°方向概率为:
Figure 537686DEST_PATH_IMAGE027
,其中
Figure 452552DEST_PATH_IMAGE028
;90°方向概率为:
Figure 583319DEST_PATH_IMAGE029
,其中
Figure 925439DEST_PATH_IMAGE030
;135°方向概率为:
Figure 774446DEST_PATH_IMAGE031
,其中
Figure 476692DEST_PATH_IMAGE032
104.获取每个缺陷窗口各个方向的灰度游程矩阵,根据每个窗口各个方向的灰度游程矩阵以及对应方向的方向概率获取每个缺陷窗口的最终灰度游程矩阵。
获取每个缺陷窗口的最终灰度游程矩阵的方法为:
根据每个缺陷窗口中0°、45°、90°、135°方向上的像素点数量获取对应的方向概率;
获取每个缺陷窗口在0°、45°、90°、135°方向上的灰度游程矩阵,根据每个方向上的灰度游程矩阵及对应方向的方向概率进行加权求和,得到每个缺陷窗口的最终灰度游程矩阵。
根据各个方向的概率作为各方向灰度游程矩阵的权值,最终统计灰度游程矩阵:
Figure 461965DEST_PATH_IMAGE033
,式中P为各方向上的权重,P(i,j)为各方向得到的灰度游程矩阵。
105.将每个缺陷窗口的最终灰度游程矩阵平均划分为长游程区域和短游程区域。
将每个缺陷窗口的最终灰度游程矩阵沿竖直方向划分为2个大小相同的矩形,划分示意图如图2所示。
图中A1区域矩阵的游程短,A2区域矩阵的游程较长,游程的越长说明当前缺陷为划痕的可能性越大,反之游程越短时说明当前缺陷为裂纹的可能性越大。
统计当前窗口灰度游程矩阵各区域像素长游程与短游程的频数差异,根据频数差异初次分析当前缺陷的类型,统计A1区域的短游程频数计算公式为:
Figure 37303DEST_PATH_IMAGE034
,式中
Figure 311290DEST_PATH_IMAGE003
表示短游程的频数,i为第i个灰度阶,j为游程长度,区间为(1,
Figure 630275DEST_PATH_IMAGE035
),P(i,j)为灰度值为i,长度为j的游程频数,
Figure 673318DEST_PATH_IMAGE003
的值越大,短游程频数越大。
统计A2区域的长游程频数计算公式:
Figure 311235DEST_PATH_IMAGE036
Figure 72517DEST_PATH_IMAGE004
表示长游程的频数,i为第i个灰度阶,j为游程长度,区间为(
Figure 195194DEST_PATH_IMAGE037
),P(i,j)为灰度值为i,长度为j的游程频数,
Figure 92743DEST_PATH_IMAGE004
的值越大,长游程频数越大。
根据最终灰度游程矩阵中长游程区域和短游程区域中元素的频数计算每个缺陷窗口的游程频数差异度。
计算每个缺陷窗口的游程频数差异度的方法为:
根据游程长度将最终灰度游程矩阵平均分为长游程区域和短游程区域,分别统计长游程区域和短游程区域中所有元素的频数;
根据长游程区域中所有元素的频数和短游程区域中所有元素的频数计算每个缺陷窗口的游程频数差异度,表达式为:
Figure 744304DEST_PATH_IMAGE001
其中,
Figure 55200DEST_PATH_IMAGE002
表示每个缺陷窗口的游程频数差异度,
Figure 902939DEST_PATH_IMAGE003
表示短游程的频数,
Figure 717311DEST_PATH_IMAGE004
表示长游程的频数,当
Figure 8615DEST_PATH_IMAGE002
>1时说明短游程的频数更高,反映出当前窗口内的缺陷是裂纹的可能性更大。反之,当
Figure 541228DEST_PATH_IMAGE002
<1时说明长游程的频数更高,反映出当前窗口内的缺陷是划痕的可能性更高。
106.将每个缺陷窗口的最终灰度游程矩阵平均划分为高灰度阶区域和低灰度阶区域。
将窗口的灰度游程矩阵沿水平方向划分为两个大小相等的区域,划分示意图如图3所示。
图中A3区域为低灰度区域,A4为高灰度区域。灰度值越高说明当前缺陷为划痕的可能性也越高,反之当灰度值越低是认为当前缺陷为裂纹的可能性越高。
根据最终灰度游程矩阵中高灰度阶区域和低灰度阶区域中元素的频数计算每个缺陷窗口的灰度频数差异度。
计算每个缺陷窗口的灰度频数差异度的方法为:
根据灰度阶将最终灰度游程矩阵平均划分为高灰度阶区域和低灰度阶区域,分别统计高灰度阶区域和低灰度阶区域中所有元素的频数;
根据高灰度阶区域和低灰度阶区域中所有元素的频数计算每个缺陷窗口的灰度频数差异度,表达式为:
Figure 943390DEST_PATH_IMAGE005
其中,
Figure 612269DEST_PATH_IMAGE006
表示每个缺陷窗口的灰度频数差异度,
Figure 74474DEST_PATH_IMAGE007
表示游程长度为j,灰度阶为i的像素频数,m表示最终灰度游程矩阵的游程长度,
Figure 94383DEST_PATH_IMAGE008
表示最终灰度游程矩阵的灰度阶长度,
Figure 782460DEST_PATH_IMAGE009
表示最终灰度游程矩阵的第一阶灰度阶,当
Figure 571424DEST_PATH_IMAGE038
时说明当前缺陷为裂纹的可能性越大,反之
Figure 204531DEST_PATH_IMAGE039
1说明当前缺陷为划痕的可能性越大。
107.根据每个缺陷窗口的游程频数差异度和灰度频数差异度计算每个缺陷窗口的置信度,根据每个缺陷窗口的置信度判断该窗口的缺陷类别。
裂纹与齿痕的区别无法根据单一的元素直接识别出来,需要灰度值和游程长度二者结合,能够很大程度上提高识别准确率。
计算每个缺陷窗口的置信度的方法为:
根据每个缺陷窗口的游程频数差异度和灰度频数差异度计算每个缺陷窗口的置信度,表达式为:
Figure 446157DEST_PATH_IMAGE010
其中,
Figure 190122DEST_PATH_IMAGE011
表示每个缺陷窗口的置信度,
Figure 99172DEST_PATH_IMAGE002
表示每个缺陷窗口的游程频数差异度,
Figure 637601DEST_PATH_IMAGE006
表示每个缺陷窗口的灰度频数差异度,当
Figure 632101DEST_PATH_IMAGE011
的值越大说明矩阵内的低灰度短游程的频数较多,反映缺陷趋向于裂纹的可能性更大,反之当
Figure 163446DEST_PATH_IMAGE011
越小时说明矩阵内高灰度长游程的像素频数较多,反映缺陷趋向于划痕的可能性更大。
像素低灰度短游程区域的特征与裂纹的纹理特征十分匹配,高灰度短游程区域可能包含一小部***纹区域的像素也可能包含一小部分划痕区域的像素,像素低灰度长游程区域同样可能包含一小部***纹区域的像素也可能包含一小部分划痕区域的像素,高灰度长游程区域的特征与划痕的纹理特征十分匹配,因此根据每个窗口的置信度反应该窗口区域中的游程差异和灰度差异的大小,从而判断缺陷窗口中的缺陷类型。
根据每个缺陷窗口的置信度判断该窗口的缺陷类别的方法为:
当缺陷窗口的置信度大于阈值时,该缺陷窗口中存在裂纹缺陷;
当缺陷窗口的置信度小于阈值时,该缺陷窗口中存在划痕。
置信度
Figure 927002DEST_PATH_IMAGE011
趋近于1时认为当前缺陷为裂纹的可能性大,反之
Figure 636332DEST_PATH_IMAGE011
趋向于0时认为当前缺陷为划痕的可能性大,本发明设定经验阈值0.5,当
Figure 118129DEST_PATH_IMAGE040
时认为当前缺陷为裂纹,当
Figure 203897DEST_PATH_IMAGE041
时把当前缺陷当作划痕来处理。
本发明首先通过边缘检测对齿轮图像中的像素点进行初步判断,从而对获得的疑似缺陷窗口进行灰度游程矩阵的建立,大大减少了计算量,同时在建立灰度游程矩阵时通过对窗口中缺陷像素点在各个方向的分布对窗口各个方向的游程矩阵进行加权,得到最终的灰度游程矩阵能够更加完整的表征缺陷窗口中像素点的特征,进而结合矩阵中像素的灰度频数差异特征、游程长度差异特征,双重特征紧密结合纹理信息,能够更加准确的识别缺陷类型。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种汽车零部件精密识别方法,其特征在于,包括:
获取齿轮灰度图像,对灰度图像进行窗口分割,得到多个窗口;
对每个窗口进行边缘检测,将检测得到的边缘像素点的个数占比大于阈值的窗口作为缺陷窗口;设定当检测出的边缘上的像素数量占窗口像素总数量的大于阈值
Figure DEST_PATH_IMAGE001
时认为需要对缺陷进行后续识别处理;
计算每个缺陷窗口中的边缘像素点的缺陷程度,将缺陷程度大于阈值的边缘像素点作为缺陷像素点,获取每个缺陷窗口中的所有缺陷像素点;
以每个缺陷像素点为中心,统计每个缺陷像素点邻域上各个方向的缺陷像素点个数,根据每个方向上缺陷像素点的个数计算每个方向的方向概率,方法如下:
统计所有缺陷点邻域内同灰度值像素相对于中心像素点的方向角度,对这些方向角度做依次概率分布统计,根据方向角度为0°和180°的缺陷像素点数量与所有缺陷像素点数量的比值计算0°方向概率;
同理,根据方向角度为45°和225°的缺陷像素点数量与所有缺陷像素点数量的比值计算45°方向概率;根据方向角度为90°和270°的缺陷像素点数量与所有缺陷像素点数量的比值计算90°方向概率;根据方向角度为135°和315°的缺陷像素点数量与所有缺陷像素点数量的比值计算135°方向概率;
获取每个缺陷窗口各个方向的灰度游程矩阵,根据每个窗口各个方向的灰度游程矩阵以及对应方向的方向概率获取每个缺陷窗口的最终灰度游程矩阵;
将每个缺陷窗口的最终灰度游程矩阵平均划分为长游程区域和短游程区域;
根据最终灰度游程矩阵中长游程区域和短游程区域中元素的频数计算每个缺陷窗口的游程频数差异度;
将每个缺陷窗口的最终灰度游程矩阵平均划分为高灰度阶区域和低灰度阶区域;
根据最终灰度游程矩阵中高灰度阶区域和低灰度阶区域中元素的频数计算每个缺陷窗口的灰度频数差异度;
根据每个缺陷窗口的游程频数差异度和灰度频数差异度计算每个缺陷窗口的置信度,根据每个缺陷窗口的置信度判断该窗口的缺陷类别。
2.根据权利要求1所述的一种汽车零部件精密识别方法,其特征在于,计算每个缺陷窗口中的边缘像素点的缺陷程度的方法为:
获取每个边缘像素点的邻域内与该边缘像素点灰度值相同的像素点个数;
对每个边缘像素点的邻域中像素点依次进行遍历,获取该边缘像素点灰度值相同的连续像素点个数,作为每个边缘像素点的像素聚集程度;
根据每个边缘像素点的像素聚集程度以及邻域内与该边缘像素点灰度值相同的像素点个数计算每个边缘像素点的缺陷程度。
3.根据权利要求1所述的一种汽车零部件精密识别方法,其特征在于,获取每个缺陷窗口的最终灰度游程矩阵的方法为:
根据每个缺陷窗口中0°、45°、90°、135°方向上的像素点数量获取对应的方向概率;
获取每个缺陷窗口在0°、45°、90°、135°方向上的灰度游程矩阵,根据每个方向上的灰度游程矩阵及对应方向的方向概率进行加权求和,得到每个缺陷窗口的最终灰度游程矩阵。
4.根据权利要求1所述的一种汽车零部件精密识别方法,其特征在于,计算每个缺陷窗口的游程频数差异度的方法为:
根据游程长度将最终灰度游程矩阵平均分为长游程区域和短游程区域,分别统计长游程区域和短游程区域中所有元素的频数;
根据长游程区域中所有元素的频数和短游程区域中所有元素的频数计算每个缺陷窗口的游程频数差异度,表达式为:
Figure DEST_PATH_IMAGE003
其中,
Figure 799825DEST_PATH_IMAGE004
表示每个缺陷窗口的游程频数差异度,
Figure DEST_PATH_IMAGE005
表示短游程的频数,
Figure 601559DEST_PATH_IMAGE006
表示长游程的频数。
5.根据权利要求1所述的一种汽车零部件精密识别方法,其特征在于,计算每个缺陷窗口的灰度频数差异度的方法为:
根据灰度阶将最终灰度游程矩阵平均划分为高灰度阶区域和低灰度阶区域,分别统计高灰度阶区域和低灰度阶区域中所有元素的频数;
根据高灰度阶区域和低灰度阶区域中所有元素的频数计算每个缺陷窗口的灰度频数差异度,表达式为:
Figure 166664DEST_PATH_IMAGE008
其中,
Figure DEST_PATH_IMAGE009
表示每个缺陷窗口的灰度频数差异度,
Figure 661230DEST_PATH_IMAGE010
表示游程长度为j,灰度阶为i的像素频数,m表示最终灰度游程矩阵的游程长度,
Figure DEST_PATH_IMAGE011
表示最终灰度游程矩阵的灰度阶长度,
Figure 865946DEST_PATH_IMAGE012
表示最终灰度游程矩阵的第一阶灰度阶。
6.根据权利要求1所述的一种汽车零部件精密识别方法,其特征在于,计算每个缺陷窗口的置信度的方法为:
根据每个缺陷窗口的游程频数差异度和灰度频数差异度计算每个缺陷窗口的置信度,表达式为:
Figure 533688DEST_PATH_IMAGE014
其中,
Figure DEST_PATH_IMAGE015
表示每个缺陷窗口的置信度,
Figure 140250DEST_PATH_IMAGE004
表示每个缺陷窗口的游程频数差异度,
Figure 602455DEST_PATH_IMAGE009
表示每个缺陷窗口的灰度频数差异度。
7.根据权利要求1所述的一种汽车零部件精密识别方法,其特征在于,根据每个缺陷窗口的置信度判断该窗口的缺陷类别的方法为:
当缺陷窗口的置信度大于阈值时,该缺陷窗口中存在裂纹缺陷;
当缺陷窗口的置信度小于阈值时,该缺陷窗口中存在划痕。
CN202210822912.1A 2022-07-14 2022-07-14 一种汽车零部件精密识别方法 Active CN115082418B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210822912.1A CN115082418B (zh) 2022-07-14 2022-07-14 一种汽车零部件精密识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210822912.1A CN115082418B (zh) 2022-07-14 2022-07-14 一种汽车零部件精密识别方法

Publications (2)

Publication Number Publication Date
CN115082418A CN115082418A (zh) 2022-09-20
CN115082418B true CN115082418B (zh) 2022-11-04

Family

ID=83259462

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210822912.1A Active CN115082418B (zh) 2022-07-14 2022-07-14 一种汽车零部件精密识别方法

Country Status (1)

Country Link
CN (1) CN115082418B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115270894B (zh) * 2022-09-27 2022-12-27 山东行创科技有限公司 用于甲醇发电机的异常振动信号提取方法
CN115330773B (zh) * 2022-10-12 2023-03-28 南通红运金属科技有限公司 一种金属研磨麻点缺陷检测方法
CN115359053B (zh) * 2022-10-19 2023-12-19 天津云象科技发展有限公司 一种金属板材缺陷智能检测方法及***
CN115526890B (zh) * 2022-11-25 2023-03-24 深圳市腾泰博科技有限公司 一种电唱机唱头故障因素识别方法
CN115731239B (zh) * 2023-01-06 2023-04-11 山东建筑大学鉴定检测中心有限公司 一种钢材质量的检测方法
CN115953398B (zh) * 2023-03-13 2023-05-26 无锡斯达新能源科技股份有限公司 一种用于带钢表面的缺陷识别方法
CN116229335B (zh) * 2023-05-10 2023-07-21 山东大佳机械有限公司 基于图像数据的畜禽养殖场环境识别方法
CN116309561B (zh) * 2023-05-17 2023-08-04 威海赛宝工业信息技术研究院有限公司 一种基于防漏电绝缘材料的表面缺陷识别方法
CN116468728B (zh) * 2023-06-20 2023-08-29 深圳市锦红兴科技有限公司 一种五金冲压件模具生产质量检测方法
CN116993722B (zh) * 2023-09-26 2023-12-08 山东奥斯登房车有限公司 一种拖挂式房车车身缺陷视觉检测方法
CN117173703B (zh) * 2023-11-02 2024-01-16 温州华嘉电器有限公司 一种隔离开关状态识别方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006331163A (ja) * 2005-05-27 2006-12-07 Matsushita Electric Ind Co Ltd 画像処理装置
CN112998649A (zh) * 2015-01-06 2021-06-22 大卫·伯顿 移动式可穿戴的监控***
CN114049489A (zh) * 2022-01-13 2022-02-15 山东鹰联光电科技股份有限公司 一种光纤电缆护套破裂预警方法
CN114202541A (zh) * 2022-02-17 2022-03-18 湖北中海电线电缆有限公司 基于人工智能的电缆缺陷检测方法
CN114399635A (zh) * 2022-03-25 2022-04-26 珞石(北京)科技有限公司 基于特征定义和深度学习的图像二分类集成学习方法
CN114419025A (zh) * 2022-01-27 2022-04-29 江苏泰和木业有限公司 一种基于图像处理的纤维板质量评估方法
CN114529549A (zh) * 2022-04-25 2022-05-24 南通东德纺织科技有限公司 一种基于机器视觉的布匹缺陷标注方法及***
CN114708325A (zh) * 2022-03-31 2022-07-05 南通华烨塑料工业有限公司 一种基于橡胶喷霜缺陷的橡胶生产问题快速定位方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2793164A1 (en) * 2013-04-19 2014-10-22 Stichting Maastricht Radiation Oncology "Maastro Clinic" Method and system for determining a phenotype of a neoplasm in a human or animal body
CN114240888B (zh) * 2021-12-17 2022-09-20 沭阳安宏木制品厂 一种基于图像处理的家具组件喷漆缺陷修复方法及***
CN114445387A (zh) * 2022-01-29 2022-05-06 泗阳富艺木业股份有限公司 一种基于机器视觉的纤维板质量分类方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006331163A (ja) * 2005-05-27 2006-12-07 Matsushita Electric Ind Co Ltd 画像処理装置
CN112998649A (zh) * 2015-01-06 2021-06-22 大卫·伯顿 移动式可穿戴的监控***
CN114049489A (zh) * 2022-01-13 2022-02-15 山东鹰联光电科技股份有限公司 一种光纤电缆护套破裂预警方法
CN114419025A (zh) * 2022-01-27 2022-04-29 江苏泰和木业有限公司 一种基于图像处理的纤维板质量评估方法
CN114202541A (zh) * 2022-02-17 2022-03-18 湖北中海电线电缆有限公司 基于人工智能的电缆缺陷检测方法
CN114399635A (zh) * 2022-03-25 2022-04-26 珞石(北京)科技有限公司 基于特征定义和深度学习的图像二分类集成学习方法
CN114708325A (zh) * 2022-03-31 2022-07-05 南通华烨塑料工业有限公司 一种基于橡胶喷霜缺陷的橡胶生产问题快速定位方法
CN114529549A (zh) * 2022-04-25 2022-05-24 南通东德纺织科技有限公司 一种基于机器视觉的布匹缺陷标注方法及***

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Real-time textural defect detection based on label run length co-occurrence matrix;Bi Mingde等;《Proceedings 2011 International Conference on Transportation, Mechanical, and Electrical Engineering (TMEE)》;20120514;第271-274页 *
Texture features based microscopic image classification of liver cellular granuloma using artificial neural networks;Fuqian Shi等;《2019 IEEE 8th Joint International Information Technology and Artificial Intelligence Conference (ITAIC)》;20190805;第432-439页 *
基于灰度―游程累加的皮革缺陷自动检测;于彩香等;《皮革与化工》;20081225(第06期);第36-38、43页 *
基于纹理分析的皮革可视缺陷检测方法的研究;于彩香;《中古优秀硕士学位论文全文数据库 信息科技辑》;20100315;第2010年卷(第3期);I138-750 *

Also Published As

Publication number Publication date
CN115082418A (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
CN115082418B (zh) 一种汽车零部件精密识别方法
CN115170576B (zh) 基于机器视觉的铝管表面缺陷检测方法
CN115082467B (zh) 基于计算机视觉的建材焊接表面缺陷检测方法
CN115239704B (zh) 一种木材表面缺陷的精准检测修复方法
CN114549522B (zh) 基于目标检测的纺织品质量检测方法
CN116758061B (zh) 基于计算机视觉的铸件表面缺陷检测方法
CN114529549B (zh) 一种基于机器视觉的布匹缺陷标注方法及***
CN115345885B (zh) 一种金属健身器材外观质量检测方法
CN114842017A (zh) 一种hdmi线缆表面质量检测方法及***
CN115018838B (zh) 用于被氧化钢管材料表面麻点瑕疵的识别方法
CN114170208A (zh) 一种基于人工智能的纸制品缺陷检测方法
CN115953398B (zh) 一种用于带钢表面的缺陷识别方法
CN115100221A (zh) 一种玻璃缺陷分割方法
CN115063430B (zh) 基于图像处理的电气管道裂纹检测方法
CN114943739A (zh) 一种铝管质量检测方法
CN114820625B (zh) 一种汽车顶块缺陷检测方法
CN117011292B (zh) 一种复合板表面质量快速检测方法
CN115797361B (zh) 一种铝模板表面缺陷检测方法
CN114782330B (zh) 基于人工智能的炉排异常检测方法及***
CN117152161B (zh) 一种基于图像识别的刨花板质量检测方法及***
CN114972892B (zh) 基于计算机视觉的汽车刹车片缺陷分类方法
CN109239073A (zh) 一种用于汽车车身的表面缺陷检测方法
CN116883408B (zh) 基于人工智能的积算仪壳体缺陷检测方法
CN117635609B (zh) 一种塑胶产品生产质量视觉检测方法
CN115266732A (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