WO2021185035A1 - Method for detecting abnormal wear spot image on basis of appearance feature - Google Patents

Method for detecting abnormal wear spot image on basis of appearance feature Download PDF

Info

Publication number
WO2021185035A1
WO2021185035A1 PCT/CN2021/077660 CN2021077660W WO2021185035A1 WO 2021185035 A1 WO2021185035 A1 WO 2021185035A1 CN 2021077660 W CN2021077660 W CN 2021077660W WO 2021185035 A1 WO2021185035 A1 WO 2021185035A1
Authority
WO
WIPO (PCT)
Prior art keywords
wear scar
row
length
column
wear
Prior art date
Application number
PCT/CN2021/077660
Other languages
French (fr)
Chinese (zh)
Inventor
肖梅
杜开瑞
张雷
徐婷
Original Assignee
长安大学
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 长安大学 filed Critical 长安大学
Publication of WO2021185035A1 publication Critical patent/WO2021185035A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw 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
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume

Definitions

  • the invention relates to the extended application of a four-ball friction tester to determine the lubricating performance of lubricating oil, in particular to a method for detecting abnormal wear scar images based on shape characteristics.
  • the lubricating oil with good lubricity can protect the machine and prolong the working life. It is often measured by the four-ball wear tester. According to China's petrochemical industry standards (GB-T 12583-1998 and HT 0762-2005), the determination method of lubricant anti-wear performance is: clamp three steel balls with a diameter of 12.7 mm in an oil box, and test After oil immersion, place a steel ball with a diameter on the top of the three balls.
  • the top ball rotates at a certain speed for 60min, and then take out the three steel balls at the bottom , Measure the wear scar diameters of three steel balls under a microscope with a measurement accuracy of 0.01mm. A total of six sets of measurement data are obtained from the three steel balls. The arithmetic mean of the spot diameter is evaluated.
  • the wear scar image of the four-ball friction test is elliptical, but when the experimenter's experience is insufficient or the operation is not standardized during the operation, the wear scar image presents an abnormal shape and cannot be used to determine the friction coefficient of the lubricating oil.
  • testers mainly judge the validity of test data based on experience, and subjective judgment errors are unavoidable. Based on this, we propose a detection method for abnormal wear scar images based on shape features.
  • the purpose of the present invention is to provide a method for detecting abnormal wear scar images based on appearance characteristics, so as to solve the defect of error in the abnormal detection of the existing wear scar images.
  • the present invention provides a method for detecting abnormal wear spot images based on appearance characteristics, including the following steps:
  • Step 1 Preprocess the collected wear scar image to obtain a wear scar area map
  • Step 2 Extract the wear scar area in the wear scar area map
  • Step 3 Calculate the row and column diameters of the wear scar area
  • Step 4 Determine whether the wear scar image is abnormal according to the row and column long diameters of the wear scar level map obtained in step 3. If the wear scar image is abnormal, the algorithm ends; otherwise, go to step 5;
  • Step 5 Perform interval filtering on the row length and column length of the wear spot area to obtain the interval row length and the interval column length;
  • Step 6 Calculate the gradient change value according to the interval row length and interval column length obtained in step 5;
  • Step 7. Determine the adaptive segmentation threshold
  • Step 8 Binarize the gradient change value obtained in step 6 according to the adaptive segmentation threshold obtained in step 7, to obtain gradient binary data;
  • Step 10 Determine whether the shape of the wear scar area is normal according to the gradient binary data obtained in step 9 and the row major diameter and column major diameter of the wear scar area in step 3.
  • step 2 the wear scar area is extracted in the wear scar area map, and the specific method is:
  • step 3 the line length diameter of the wear scar area is calculated according to the following formula:
  • rh is the length of the row
  • D h (i) is the length of the i-th row
  • i 0 is the first row of the wear scar area
  • i 1 is the last row of the wear scar area
  • j 0 is the first row of the wear scar area
  • j 1 is the last row of the wear scar area
  • r l is the length of the column
  • D l (j) is the length of the j-th column
  • step 4 judging whether the wear scar image is abnormal according to the row length diameter and column length diameter of the wear scar level map obtained in step 3, the specific method is that if the row length diameter and the column length diameter satisfy the following formula, it is considered Abnormal wear spot image:
  • is the gap threshold.
  • interval filtering is performed on the row length and column length of the wear spot area in combination with the interval length set by the following formula:
  • s and t are the indexes of row interval and column interval respectively;
  • L h (s) represents the row length of the sth interval;
  • L l (t) is the column length of the t-th interval;
  • q is the set interval length; [] Represents the round-down operator.
  • step 6 the gradient change value is calculated by the interval row length and interval column length obtained in step 5 of the following structure:
  • I h (s) is the row gradient in the s interval
  • I l (t) is the column gradient in the t interval.
  • step 7 the adaptive segmentation threshold is determined, and the specific method is:
  • the order of The corresponding element value is used as the column segmentation threshold T 1 ;
  • step 8 the gradient change obtained in step 6 is binarized according to the adaptive segmentation threshold obtained in step 7 to obtain gradient binary data, and the specific method is calculated using the following formula:
  • B h (s) is the row gradient binary data
  • B l (t) is the column gradient binary data.
  • step 10 judging whether the shape of the wear scar area is normal according to the gradient binary data obtained in step 9 and the row length diameter and column length diameter of the wear scar area in step 3, the specific method is:
  • the integers m and n are arbitrarily selected, and when the integers satisfy the following formula at the same time, it is considered that the shape of the row direction is abnormal:
  • is the length threshold
  • the present invention provides a method for detecting abnormal wear scar images based on shape features. Based on the segmentation of the wear scar area and detection of the direction angle, the row and column lengths are used to characterize the row and column shape characteristics of the wear scar area. ; On the basis of interval filtering, the gradient change values of row length and column length are calculated separately; adaptive threshold is used to binarize the gradient change value, and based on the monotonicity of the gradient, it is determined whether the wear spot image is abnormal.
  • the method of the present invention is based on the abnormal state of the wear spot image based on the shape feature, and has more universal applicability, avoids the error of subjective judgment, and has higher accuracy. Provide a theoretical basis for the availability of test data.
  • Figure 1 is a wear spot image F
  • Figure 2 is a diagram f of the wear scar area
  • Figure 3 is a graph g of the level of wear scars.
  • the present invention provides a method for detecting abnormal wear spot images based on appearance characteristics, including the following steps:
  • i and j represent the row and column of the pixel respectively, and i and j are integers, satisfying: 1 ⁇ i ⁇ M And 1 ⁇ j ⁇ N holds.
  • the wear spot image F is shown in FIG. 1.
  • Step 2 Automatic segmentation of the wear scar area.
  • the wear scar area map f is shown in FIG. 2.
  • Step 3 Detect the direction angle of the wear scar.
  • a four-ball friction test wear scar image wear scar direction angle automatic measurement method ZL201710018314.8
  • Step 4 Rotate the wear scar to the direction of the line axis according to the direction angle of the wear scar. Rotate the wear scar area map f clockwise around the center of the image by w degrees to obtain the wear scar horizontal map, which is represented by the symbol g. At this time, the wear scar direction is consistent with the row axis direction (horizontal direction).
  • the size of the image is changed from M ⁇ N to M w ⁇ N w , and i 0 and i 1 are used to represent the first and last rows of the wear scar area in the wear scar level map, which satisfies the calculation formula (1-2) ;
  • Step 5 Calculate the line length.
  • the wear spot area is elliptical, and the row length increases with the increase of the row number to present an inverted U shape, that is, the row length curve first rises and then drops, and the largest row length is the row length diameter.
  • the row length and row length diameter are shown in formula (5-6).
  • D h (i) represents the length of the i-th row; rh is the length of the row.
  • Step 6 Calculate the column length.
  • the column-wise length calculation method is the same as the column-wise length.
  • the length of each column represents the column-wise length.
  • the column-wise length of the wear scar area increases with the column number in an inverted U shape.
  • the length of the jth column is D l (j) and the column
  • the long diameter r l is shown in formula (7-8):
  • Step 7 Detection of abnormal images based on the long diameter. If there is a large difference between the row length diameter and the column length diameter, and the equation (9) is satisfied, the wear spot image is considered to be abnormal, and then go to step 13; otherwise, go to step 8:
  • is the gap threshold, and the value is 0.167 based on experience.
  • Step 8 Perform interval filtering on the row length and column length of the wear spot area.
  • the wear scar image is composed of wear scars, and the edges of the detected wear scar area are often not smooth.
  • the row length and column length are filtered according to the set interval length to make the small jagged edges smoother. After interval filtering, the interval row length and interval column length of the wear spot area are obtained, and the calculation formula (10-11) is as follows:
  • s and t are the indexes of the row interval and column interval respectively, which are integers;
  • L h (s) represents the row length of the s-th interval;
  • L l (t) is the column length of the t-th interval;
  • q is the set interval length , Usually between 3-20; [] represents the round-down operator.
  • Step 9 Calculate the gradient change value.
  • the gradient change of the data can reflect the monotonicity of the data curve. When the gradient value is positive, the curve increases monotonously, and vice versa, it can be judged whether the shape of the wear scar area exhibits first increase and then decrease characteristics.
  • the gradient of the interval row length and column length after interval filtering is shown in formula (12-13):
  • I h (s) is the row gradient in the s interval
  • I l (t) is the column gradient in the t interval.
  • Step 10 Determine the adaptive segmentation threshold.
  • the method uses an adaptive threshold to perform binary segmentation on the gradient change value.
  • the specific operation is: The elements in are arranged in ascending order, and the order is in the first
  • the element of the bit is selected as the row segmentation threshold, denoted as Th ; in the same way, the column segmentation threshold T 1 can be determined: set After the elements are arranged in ascending order, the order of The corresponding element value is the column segmentation threshold T 1 .
  • Step 11 Gradient binarization. Binarize the gradient data based on the adaptive threshold to obtain gradient binary data.
  • the data with a value of 1 corresponds to a monotonic increasing curve (that is, the rising section of the curve), and the rest are monotonic decreasing curves. If the monotonic increasing curve segment is not monotonously decreasing
  • the curve separation means that the data change conforms to the change characteristics of first increase and then decrease.
  • Step 12 Detection of abnormal wear scar images.
  • Normal image data can be used for subsequent analysis, while abnormal image data cannot be used for subsequent analysis and need to be re-tested.
  • is the length threshold, and its value is 2-10.
  • Step 13 The method ends.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Geometry (AREA)
  • Constituent Portions Of Griding Lathes, Driving, Sensing And Control (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Analysis (AREA)

Abstract

Provided in the present invention is a method for detecting an abnormal wear spot image on the basis of an appearance feature. The method comprises: on the basis of segmenting a wear mark region and measuring a direction angle, representing row-direction and column-direction appearance features of the wear mark region by means of row-direction and column-direction lengths; on the basis of interval filtering, respectively calculating gradient change values of a row length and a column length; and performing binarization on the gradient change values by using an adaptive threshold value, and determining whether a wear spot image is abnormal on the basis of the monotonicity of the gradient. In the present invention, the abnormal state of a wear spot image is determined on the basis of an appearance feature. The invention is more universal, avoids errors of subjective determination, achieves higher accuracy, and provides a theoretical basis for the availability of test data.

Description

一种基于外形特征的异常磨斑图像的检测方法A Detection Method of Abnormal Wear Spot Images Based on Shape Features 技术领域Technical field
本发明涉及一种四球摩擦试验机测定润滑油润滑性能的扩展应用,具体涉及一种基于外形特征的异常磨斑图像的检测方法。The invention relates to the extended application of a four-ball friction tester to determine the lubricating performance of lubricating oil, in particular to a method for detecting abnormal wear scar images based on shape characteristics.
背景技术Background technique
润滑性好的润滑油可以保护机械、延长工作寿命,常常通过四球磨损试验机来测量。根据我国石油化工行业标准(GB-T 12583-1998和H-T 0762-2005),润滑剂抗摩损性能测定方法为:将三个直径为12.7mm的钢球夹紧在一油盒中,并用试油浸没,在三球顶部放置一个直径的钢球,在试油温度达到75℃±2℃后,施加147N或392N作用力,顶球在一定转速下旋转60min,随后取出底部的三个钢球,在测量精度为0.01mm的显微镜下测量三个钢球的磨斑直径,三个钢球共得到六组测量数据,润滑油或润滑脂的抗摩性能通过三个球的六次测量的摩斑直径的算术平均值来评价。理想状态下,四球摩擦试验的磨斑图像呈现椭圆形,但当操作过程中实验人员的经验不足或操作不规范时,磨斑图像呈现异常形状,不能用于测定润滑油的摩擦系数。目前主要由试验员依据经验判定试验数据的有效性,不可避免的产生主观判定的误差。基于此,我们提出了一种基于外形特征的异常磨斑图像的检测方法。The lubricating oil with good lubricity can protect the machine and prolong the working life. It is often measured by the four-ball wear tester. According to China's petrochemical industry standards (GB-T 12583-1998 and HT 0762-2005), the determination method of lubricant anti-wear performance is: clamp three steel balls with a diameter of 12.7 mm in an oil box, and test After oil immersion, place a steel ball with a diameter on the top of the three balls. After the oil test temperature reaches 75℃±2℃, apply 147N or 392N force, the top ball rotates at a certain speed for 60min, and then take out the three steel balls at the bottom , Measure the wear scar diameters of three steel balls under a microscope with a measurement accuracy of 0.01mm. A total of six sets of measurement data are obtained from the three steel balls. The arithmetic mean of the spot diameter is evaluated. In an ideal state, the wear scar image of the four-ball friction test is elliptical, but when the experimenter's experience is insufficient or the operation is not standardized during the operation, the wear scar image presents an abnormal shape and cannot be used to determine the friction coefficient of the lubricating oil. At present, testers mainly judge the validity of test data based on experience, and subjective judgment errors are unavoidable. Based on this, we propose a detection method for abnormal wear scar images based on shape features.
发明内容Summary of the invention
本发明的目的在于提供一种基于外形特征的异常磨斑图像的检测方法,解决现有的磨斑图像的异常检测存在误差的缺陷。The purpose of the present invention is to provide a method for detecting abnormal wear scar images based on appearance characteristics, so as to solve the defect of error in the abnormal detection of the existing wear scar images.
为了达到上述目的,本发明采用的技术方案是:In order to achieve the above objective, the technical solution adopted by the present invention is:
本发明提供的一种基于外形特征的异常磨斑图像的检测方法,包括如下步骤:The present invention provides a method for detecting abnormal wear spot images based on appearance characteristics, including the following steps:
步骤1,对采集到的磨斑图像进行预处理,得到磨痕区域图;Step 1. Preprocess the collected wear scar image to obtain a wear scar area map;
步骤2,在磨痕区域图中提取磨痕区域;Step 2: Extract the wear scar area in the wear scar area map;
步骤3,计算磨痕区域的行长径和列长径;Step 3. Calculate the row and column diameters of the wear scar area;
步骤4,根据步骤3中得到的磨痕水平图的行长径和列长径判断磨斑图像是否异常,若磨斑图像为异常时,算法结束;否则转入步骤5;Step 4. Determine whether the wear scar image is abnormal according to the row and column long diameters of the wear scar level map obtained in step 3. If the wear scar image is abnormal, the algorithm ends; otherwise, go to step 5;
步骤5,对磨斑区域的行长度和列长度进行区间滤波,得到区间行长和区间列长;Step 5: Perform interval filtering on the row length and column length of the wear spot area to obtain the interval row length and the interval column length;
步骤6,根据步骤5中得到的区间行长和区间列长计算梯度变化值;Step 6. Calculate the gradient change value according to the interval row length and interval column length obtained in step 5;
步骤7,确定自适应分割阈值;Step 7. Determine the adaptive segmentation threshold;
步骤8,根据步骤7中得到的自适应分割阈值对步骤6中得到的梯度变化值进行二值化处理,得到梯度二值数据;Step 8. Binarize the gradient change value obtained in step 6 according to the adaptive segmentation threshold obtained in step 7, to obtain gradient binary data;
步骤10,根据步骤9得到的梯度二值数据和步骤3中的磨痕区域的行长径和列长径判断磨痕区域的形状是否正常。Step 10: Determine whether the shape of the wear scar area is normal according to the gradient binary data obtained in step 9 and the row major diameter and column major diameter of the wear scar area in step 3.
优选地,步骤2中,在磨痕区域图中提取磨痕区域,具体方法是:Preferably, in step 2, the wear scar area is extracted in the wear scar area map, and the specific method is:
S1,确定步骤1中得到的磨痕区域图的磨痕方向角;S1, determine the wear scar direction angle of the wear scar area map obtained in step 1;
S2,将磨痕区域图绕图像中心顺时针旋转磨痕方向角度,得到磨痕水平图;S2: Rotate the wear scar area map clockwise around the center of the image to obtain a wear scar level map;
S3,在磨痕水平图中提取磨痕区域,其中,利用i 0和i 1表示磨痕水平图中磨痕区域的首行和尾行;利用j 0和j 1表示磨痕水平图中磨痕区域的首列和尾列: S3, extract the wear scar area in the wear scar level map, where i 0 and i 1 are used to represent the first and last rows of the wear scar area in the wear scar level graph; j 0 and j 1 are used to represent the wear scar in the wear scar level graph The first and last columns of the area:
i 0<i且f(i 0,j)=f(i,j)=1,
Figure PCTCN2021077660-appb-000001
i 1>i且f(i 1,j)=f(i,j)=1,
Figure PCTCN2021077660-appb-000002
j 0<j且f(i,j 0)=f(i,j)=1,
Figure PCTCN2021077660-appb-000003
j 1>j且f(i,j 1)=f(i,j)=1,
Figure PCTCN2021077660-appb-000004
i 0 <i and f(i 0 ,j)=f(i,j)=1,
Figure PCTCN2021077660-appb-000001
i 1 >i and f(i 1 ,j)=f(i,j)=1,
Figure PCTCN2021077660-appb-000002
j 0 <j and f(i,j 0 )=f(i,j)=1,
Figure PCTCN2021077660-appb-000003
j 1 >j and f(i,j 1 )=f(i,j)=1,
Figure PCTCN2021077660-appb-000004
优选地,步骤3中,根据下式计算磨痕区域的行长径:Preferably, in step 3, the line length diameter of the wear scar area is calculated according to the following formula:
Figure PCTCN2021077660-appb-000005
Figure PCTCN2021077660-appb-000005
其中,r h为行长径;D h(i)为第i行的长度,
Figure PCTCN2021077660-appb-000006
i 0为磨痕区域的首行;i 1为磨痕区域的尾行;j 0为磨痕区域的首列;j 1为磨痕区域的尾列;
Among them, rh is the length of the row; D h (i) is the length of the i-th row,
Figure PCTCN2021077660-appb-000006
i 0 is the first row of the wear scar area; i 1 is the last row of the wear scar area ; j 0 is the first row of the wear scar area; j 1 is the last row of the wear scar area;
通过下式计算列长径:Calculate the column long diameter by the following formula:
Figure PCTCN2021077660-appb-000007
Figure PCTCN2021077660-appb-000007
其中,r l为列长径;D l(j)为第j列的长度,
Figure PCTCN2021077660-appb-000008
Among them, r l is the length of the column; D l (j) is the length of the j-th column,
Figure PCTCN2021077660-appb-000008
优选地,步骤4中,根据步骤3中得到的磨痕水平图的行长径和列长径判断磨斑图像是否异常,具体方法是,若行长径和列长径满足下式,则认为磨斑图像异常:Preferably, in step 4, judging whether the wear scar image is abnormal according to the row length diameter and column length diameter of the wear scar level map obtained in step 3, the specific method is that if the row length diameter and the column length diameter satisfy the following formula, it is considered Abnormal wear spot image:
Figure PCTCN2021077660-appb-000009
Figure PCTCN2021077660-appb-000009
其中,β为差距阈值。Among them, β is the gap threshold.
优选地,步骤5中,通过下式结合设定的区间长度对磨斑区域的行长度和列长度进行区间滤波:Preferably, in step 5, interval filtering is performed on the row length and column length of the wear spot area in combination with the interval length set by the following formula:
Figure PCTCN2021077660-appb-000010
Figure PCTCN2021077660-appb-000010
Figure PCTCN2021077660-appb-000011
Figure PCTCN2021077660-appb-000011
其中,s和t分别行区间和列区间的索引;L h(s)表示第s个区间行长;L l(t)为第t个区间列长;q为设定的区间长度;[]表示向下取整运算符。 Among them, s and t are the indexes of row interval and column interval respectively; L h (s) represents the row length of the sth interval; L l (t) is the column length of the t-th interval; q is the set interval length; [] Represents the round-down operator.
优选地,步骤6中,通过下式结构步骤5中得到的区间行长和区间列长计算梯度变化值:Preferably, in step 6, the gradient change value is calculated by the interval row length and interval column length obtained in step 5 of the following structure:
Figure PCTCN2021077660-appb-000012
Figure PCTCN2021077660-appb-000012
Figure PCTCN2021077660-appb-000013
Figure PCTCN2021077660-appb-000013
其中,I h(s)为s区间的行梯度,I l(t)为t区间的列梯度。 Among them, I h (s) is the row gradient in the s interval, and I l (t) is the column gradient in the t interval.
优选地,步骤7中,确定自适应分割阈值,具体方法是:Preferably, in step 7, the adaptive segmentation threshold is determined, and the specific method is:
将集合
Figure PCTCN2021077660-appb-000014
中的元素按升序排列,将位序
Figure PCTCN2021077660-appb-000015
位对应的元素作为行分割阈值T h
Will gather
Figure PCTCN2021077660-appb-000014
The elements in are arranged in ascending order, and the order of
Figure PCTCN2021077660-appb-000015
The element corresponding to the bit is used as the row segmentation threshold Th ;
将集合
Figure PCTCN2021077660-appb-000016
中元素升序排列后,将位序
Figure PCTCN2021077660-appb-000017
对应的元素值作为列分割阈值T 1
Will gather
Figure PCTCN2021077660-appb-000016
After the elements are arranged in ascending order, the order of
Figure PCTCN2021077660-appb-000017
The corresponding element value is used as the column segmentation threshold T 1 ;
其中,
Figure PCTCN2021077660-appb-000018
为四舍五入取整运算。
in,
Figure PCTCN2021077660-appb-000018
It is a rounding operation.
优选地,步骤8中,根据步骤7中得到的自适应分割阈值对步骤6中得到的梯度变化进行二值化,得到梯度二值数据,具体方法利用下式进行计算:Preferably, in step 8, the gradient change obtained in step 6 is binarized according to the adaptive segmentation threshold obtained in step 7 to obtain gradient binary data, and the specific method is calculated using the following formula:
Figure PCTCN2021077660-appb-000019
Figure PCTCN2021077660-appb-000019
Figure PCTCN2021077660-appb-000020
Figure PCTCN2021077660-appb-000020
其中,B h(s)为行梯度二值数据;B l(t)为列梯度二值数据。 Among them, B h (s) is the row gradient binary data; B l (t) is the column gradient binary data.
优选地,步骤10中,根据步骤9得到的梯度二值数据和步骤3中的磨痕区域的行长径和列长径判断磨痕区域的形状是否正常,具体方法是:Preferably, in step 10, judging whether the shape of the wear scar area is normal according to the gradient binary data obtained in step 9 and the row length diameter and column length diameter of the wear scar area in step 3, the specific method is:
分别判断行向形状和列向形状是否均为异常,当行向形状和列向形状均正常时,则认定该磨斑图像为正常,否则,该磨斑图像为异常。Determine whether the row-direction shape and the column-direction shape are both abnormal. When the row-direction shape and the column-direction shape are both normal, the wear spot image is determined to be normal; otherwise, the wear spot image is abnormal.
优选地,任意选取整数m,n,该整数同时满足下式时,则认为行向形状异常:Preferably, the integers m and n are arbitrarily selected, and when the integers satisfy the following formula at the same time, it is considered that the shape of the row direction is abnormal:
Figure PCTCN2021077660-appb-000021
Figure PCTCN2021077660-appb-000021
Figure PCTCN2021077660-appb-000022
Figure PCTCN2021077660-appb-000022
Figure PCTCN2021077660-appb-000023
Figure PCTCN2021077660-appb-000023
其中,α为长度阈值;Among them, α is the length threshold;
任意选取整数k,p,该整数同时满足下式时,则认为行向形状异常:If the integers k and p are selected arbitrarily, when the integers satisfy the following formula at the same time, it is considered that the shape of the row direction is abnormal:
Figure PCTCN2021077660-appb-000024
Figure PCTCN2021077660-appb-000024
Figure PCTCN2021077660-appb-000025
Figure PCTCN2021077660-appb-000025
Figure PCTCN2021077660-appb-000026
Figure PCTCN2021077660-appb-000026
与现有技术相比,本发明的有益效果是:Compared with the prior art, the beneficial effects of the present invention are:
本发明提供的一种基于外形特征的异常磨斑图像的检测方法,在磨痕区域分割和方向角检测的基础上,用行向和列向长度表征磨痕区域的行向和列向外形特征;在区间滤波的基础上,分别计算行长和列长的梯度变化值;采用自适应阈值对梯度变化值进行二值化,并基于梯度的单调性判定磨斑图像是否是异常。和目前的只利用两个方向的直径差值来判定的方法相比,本发明方法基于外形特征对磨斑图像的异常状态,更具有普适性,避免主观判定的误差,精度也更高,为试验数据的可用性提供理论依据。The present invention provides a method for detecting abnormal wear scar images based on shape features. Based on the segmentation of the wear scar area and detection of the direction angle, the row and column lengths are used to characterize the row and column shape characteristics of the wear scar area. ; On the basis of interval filtering, the gradient change values of row length and column length are calculated separately; adaptive threshold is used to binarize the gradient change value, and based on the monotonicity of the gradient, it is determined whether the wear spot image is abnormal. Compared with the current method that only uses the diameter difference in two directions to determine, the method of the present invention is based on the abnormal state of the wear spot image based on the shape feature, and has more universal applicability, avoids the error of subjective judgment, and has higher accuracy. Provide a theoretical basis for the availability of test data.
附图说明Description of the drawings
图1是磨斑图像F;Figure 1 is a wear spot image F;
图2是磨痕区域图f;Figure 2 is a diagram f of the wear scar area;
图3是磨痕水平图g。Figure 3 is a graph g of the level of wear scars.
具体实施方式Detailed ways
下面结合附图,对本发明进一步详细说明。The present invention will be further described in detail below in conjunction with the accompanying drawings.
本发明提供的一种基于外形特征的异常磨斑图像的检测方法,包括如下步骤:The present invention provides a method for detecting abnormal wear spot images based on appearance characteristics, including the following steps:
步骤1:通过扫面电子显微镜采集试验钢球的磨斑图像。具体是:在四球摩擦试验结束后,分别将试验所使用的三个底部钢球取出并置于扫面电镜中,并调节扫面电镜的光照和放大倍数等参数,以便清晰地采集到磨斑图像,所采集的磨斑图像用F表示。同时,得到的所述磨斑图像F的像素大小为M×N,例如磨斑图像的大小为768×1024,即M=768,N=1024。同时,利用(i,j)表示磨斑图像F的任一像素点的坐标,则i和j分别表示该像素点的行和列,且i和j均为整数,满足:1≤i≤M和1≤j≤N成立。Step 1: Acquire the wear scar image of the test steel ball through a scanning electron microscope. Specifically: after the four-ball friction test, the three bottom steel balls used in the test were taken out and placed in the scanning electron microscope, and parameters such as the illumination and magnification of the scanning electron microscope were adjusted to clearly collect the wear spots. Image, the collected wear spot image is denoted by F. At the same time, the pixel size of the obtained wear spot image F is M×N, for example, the size of the wear spot image is 768×1024, that is, M=768, N=1024. At the same time, use (i, j) to represent the coordinates of any pixel of the wear spot image F, then i and j represent the row and column of the pixel respectively, and i and j are integers, satisfying: 1≤i≤M And 1≤j≤N holds.
本实施例中,磨斑图像F如图1所示。In this embodiment, the wear spot image F is shown in FIG. 1.
步骤2:磨痕区域的自动分割。利用发明人在专利中提及“一种基于钢球磨痕梯度的磨斑直径测量方法”(ZL201310752721.3)对磨痕区域实现自动的分割,分割出来的磨痕区域,得到磨痕区域图,用f表示,f(i,j)=1的像素点表示磨痕,f(i,j)=0的像素点表示非磨痕。Step 2: Automatic segmentation of the wear scar area. Using the inventor mentioned in the patent "a method for measuring wear scar diameter based on steel ball wear scar gradient" (ZL201310752721.3), the wear scar area is automatically segmented, and the divided wear scar area is obtained to obtain a wear scar area map. It is denoted by f, the pixel with f(i,j)=1 represents the wear scar, and the pixel with f(i,j)=0 represents the non-wear scar.
本实施例中,磨痕区域图f如图2所示。In this embodiment, the wear scar area map f is shown in FIG. 2.
步骤3:检测磨痕方向角。利用发明人在专利中提及“一种四球摩擦试验磨斑图像磨痕方向角自动测定方法”(ZL201710018314.8)对磨痕方向角进行检测,检测得到的磨痕方向角用w表示。Step 3: Detect the direction angle of the wear scar. Using the inventor mentioned in the patent "a four-ball friction test wear scar image wear scar direction angle automatic measurement method" (ZL201710018314.8) to detect the direction angle of the wear scar, the detected direction angle of the wear scar is represented by w.
本实施例中,磨痕方向角为w=35°。In this embodiment, the direction angle of the wear scar is w=35°.
步骤4:根据磨痕方向角将磨痕旋转至行轴方向。将磨痕区域图f绕图像中心顺时针旋转w度,得到磨痕水平图,用符号g表示,此时磨痕方向和行轴方向(水平方向)一致。经过旋转变换后,图像的大小由M×N变成了M w×N w,用i 0和i 1表示磨痕水平图中磨痕区域的首行和尾行,满足计算式(1-2);用j 0和j 1表示磨痕水平图中磨痕区域的首列和尾列,满足计算 式(3-4);后续的操作仅针对磨痕区域进行,可减少方法的计算量。 Step 4: Rotate the wear scar to the direction of the line axis according to the direction angle of the wear scar. Rotate the wear scar area map f clockwise around the center of the image by w degrees to obtain the wear scar horizontal map, which is represented by the symbol g. At this time, the wear scar direction is consistent with the row axis direction (horizontal direction). After the rotation transformation, the size of the image is changed from M×N to M w ×N w , and i 0 and i 1 are used to represent the first and last rows of the wear scar area in the wear scar level map, which satisfies the calculation formula (1-2) ; Use j 0 and j 1 to denote the first and last columns of the wear scar area in the wear scar level map, which satisfies the calculation formula (3-4); the subsequent operations are only performed on the wear scar area, which can reduce the calculation amount of the method.
i 0<i且f(i 0,j)=f(i,j)=1,
Figure PCTCN2021077660-appb-000027
i 0 <i and f(i 0 ,j)=f(i,j)=1,
Figure PCTCN2021077660-appb-000027
i 1>i且f(i 1,j)=f(i,j)=1,
Figure PCTCN2021077660-appb-000028
i 1 >i and f(i 1 ,j)=f(i,j)=1,
Figure PCTCN2021077660-appb-000028
j 0<j且f(i,j 0)=f(i,j)=1,
Figure PCTCN2021077660-appb-000029
j 0 <j and f(i,j 0 )=f(i,j)=1,
Figure PCTCN2021077660-appb-000029
j 1>j且f(i,j 1)=f(i,j)=1,
Figure PCTCN2021077660-appb-000030
j 1 >j and f(i,j 1 )=f(i,j)=1,
Figure PCTCN2021077660-appb-000030
本实施例中,磨痕水平图g如图3所示,图像尺寸为:M w=1217,N w=1281,i 0=206,i 1=934,j 0=230,j 1=1013。 In this embodiment, the wear scar level map g is shown in FIG. 3, and the image size is: M w =1217, N w =1281, i 0 =206, i 1 =934, j 0 =230, j 1 =1013.
步骤5:计算行向长度。当磨痕处于水平方向时,通过每行磨痕的长度来表征行向长度。理想状态下,磨斑区域为椭圆形,行向长度随行号的增加呈现倒U形,即行向长度曲线先上升后下降,其中最大的行向长度为行长径。行向长度和行长径如式(5-6)所示。Step 5: Calculate the line length. When the wear marks are in the horizontal direction, the length of each line of wear marks is used to characterize the line length. In an ideal state, the wear spot area is elliptical, and the row length increases with the increase of the row number to present an inverted U shape, that is, the row length curve first rises and then drops, and the largest row length is the row length diameter. The row length and row length diameter are shown in formula (5-6).
Figure PCTCN2021077660-appb-000031
Figure PCTCN2021077660-appb-000031
Figure PCTCN2021077660-appb-000032
Figure PCTCN2021077660-appb-000032
其中,D h(i)表示第i行的长度;r h为行长径。 Among them, D h (i) represents the length of the i-th row; rh is the length of the row.
本实施例中,其行长径为r h=784。 In this embodiment, the row length diameter is rh =784.
步骤6:计算列向长度。类似地,列向长度的计算方法同列向长度,每列的长度表示列向长度,磨斑区域列向长度随列号的增加呈倒U形,第j列的长度D l(j)和列长径r l如式(7-8)所示: Step 6: Calculate the column length. Similarly, the column-wise length calculation method is the same as the column-wise length. The length of each column represents the column-wise length. The column-wise length of the wear scar area increases with the column number in an inverted U shape. The length of the jth column is D l (j) and the column The long diameter r l is shown in formula (7-8):
Figure PCTCN2021077660-appb-000033
Figure PCTCN2021077660-appb-000033
Figure PCTCN2021077660-appb-000034
Figure PCTCN2021077660-appb-000034
本实施例中,其列长径为r l=728。 In this embodiment, the column long diameter is r l =728.
步骤7:基于长径的异常图像检测。若行长径和列长径差距较大,满足式(9)时,则认为磨斑图像异常,转入步骤13;否则转入步骤8:Step 7: Detection of abnormal images based on the long diameter. If there is a large difference between the row length diameter and the column length diameter, and the equation (9) is satisfied, the wear spot image is considered to be abnormal, and then go to step 13; otherwise, go to step 8:
Figure PCTCN2021077660-appb-000035
Figure PCTCN2021077660-appb-000035
其中,β为差距阈值,按经验取值为0.167。Among them, β is the gap threshold, and the value is 0.167 based on experience.
本实施例中,
Figure PCTCN2021077660-appb-000036
In this embodiment,
Figure PCTCN2021077660-appb-000036
步骤8:对磨斑区域的行长度和列长度进行区间滤波。磨斑图像由磨痕组成,检出的磨痕区域的边缘常常不光滑,按设定的区间长度对行长度和列长度均进行均值滤波处理,可以使细小的锯齿边缘变得更光滑。区间滤波后,得到磨斑区域的区间行长和区间列长,其计算式(10-11)所示:Step 8: Perform interval filtering on the row length and column length of the wear spot area. The wear scar image is composed of wear scars, and the edges of the detected wear scar area are often not smooth. The row length and column length are filtered according to the set interval length to make the small jagged edges smoother. After interval filtering, the interval row length and interval column length of the wear spot area are obtained, and the calculation formula (10-11) is as follows:
Figure PCTCN2021077660-appb-000037
Figure PCTCN2021077660-appb-000037
Figure PCTCN2021077660-appb-000038
Figure PCTCN2021077660-appb-000038
其中,s和t分别行区间和列区间的索引,为整数;L h(s)表示第s个区间行长;L l(t)为第t个区间列长;q为设定的区间长度,通常为3-20之间;[]表示向下取整运算符。 Among them, s and t are the indexes of the row interval and column interval respectively, which are integers; L h (s) represents the row length of the s-th interval; L l (t) is the column length of the t-th interval; q is the set interval length , Usually between 3-20; [] represents the round-down operator.
本实施例中,区间长度q取值为10,s=0,1,…,72和t=0,1,…,78。In this embodiment, the interval length q is 10, s=0,1,...,72 and t=0,1,...,78.
步骤9:计算梯度变化值。数据的梯度变化可以反映出数据曲线的单调性,当梯度值为正时,曲线单调增,反之为单调减,进而可以判定磨痕区域的形状是否呈现先增后减的特性。对区间滤波后的区间行长和列长进行梯度如式(12-13)所示:Step 9: Calculate the gradient change value. The gradient change of the data can reflect the monotonicity of the data curve. When the gradient value is positive, the curve increases monotonously, and vice versa, it can be judged whether the shape of the wear scar area exhibits first increase and then decrease characteristics. The gradient of the interval row length and column length after interval filtering is shown in formula (12-13):
Figure PCTCN2021077660-appb-000039
Figure PCTCN2021077660-appb-000039
Figure PCTCN2021077660-appb-000040
Figure PCTCN2021077660-appb-000040
其中,I h(s)为s区间的行梯度,I l(t)为t区间的列梯度。 Among them, I h (s) is the row gradient in the s interval, and I l (t) is the column gradient in the t interval.
步骤10:自适应分割阈值的确定。方法利用自适应的阈值对梯度变化值进行二值分割,具体操作为:将集合
Figure PCTCN2021077660-appb-000041
中的元素按升序排列,位序排在第
Figure PCTCN2021077660-appb-000042
位的元素选定为行分割阈值,记为T h;同理可以确定列分割阈值T 1:集合
Figure PCTCN2021077660-appb-000043
中元素升序排列后,位序
Figure PCTCN2021077660-appb-000044
对应的元素值为列分割阈值T 1
Step 10: Determine the adaptive segmentation threshold. The method uses an adaptive threshold to perform binary segmentation on the gradient change value. The specific operation is:
Figure PCTCN2021077660-appb-000041
The elements in are arranged in ascending order, and the order is in the first
Figure PCTCN2021077660-appb-000042
The element of the bit is selected as the row segmentation threshold, denoted as Th ; in the same way, the column segmentation threshold T 1 can be determined: set
Figure PCTCN2021077660-appb-000043
After the elements are arranged in ascending order, the order of
Figure PCTCN2021077660-appb-000044
The corresponding element value is the column segmentation threshold T 1 .
其中,
Figure PCTCN2021077660-appb-000045
为四舍五入取整运算。
in,
Figure PCTCN2021077660-appb-000045
It is a rounding operation.
本实施例中,T h=14.4和T 1=5.7。 In this embodiment, T h = 14.4 and T 1 = 5.7.
步骤11:梯度二值化。基于自适应阈值对梯度数据进行二值化,得到梯度二值数据,值为1的数据对应单调增曲线(即曲线的上升段),其余为单调减曲线,若单调增曲线段未被单调减曲线分隔,则说明数据变化符合先增后减的变化特性。用B h(s)和B l(t)分别表示行、列梯度二值数据,计算式如(14-15)所示: Step 11: Gradient binarization. Binarize the gradient data based on the adaptive threshold to obtain gradient binary data. The data with a value of 1 corresponds to a monotonic increasing curve (that is, the rising section of the curve), and the rest are monotonic decreasing curves. If the monotonic increasing curve segment is not monotonously decreasing The curve separation means that the data change conforms to the change characteristics of first increase and then decrease. Use B h (s) and B l (t) to represent the binary data of row and column gradients respectively, and the calculation formula is as shown in (14-15):
Figure PCTCN2021077660-appb-000046
Figure PCTCN2021077660-appb-000046
Figure PCTCN2021077660-appb-000047
Figure PCTCN2021077660-appb-000047
本实施例中,行、列梯度二值数据如表1所示。In this embodiment, the row and column gradient binary data are shown in Table 1.
表1行、列梯度二值数据B h(s)和B l(t) Table 1 Row and column gradient binary data B h (s) and B l (t)
ss 1-311-31 32-3732-37 38-4738-47 48-5248-52 53-5453-54 55-5655-56 57-5857-58 59-6559-65 6666
B h(s) B h (s) 11 00 11 00 11 00 11 00 11
tt 1-501-50 51-5251-52 53-5453-54 5555 5656 57-7857-78 79-8079-80 81-10281-102 103103
B l(t) B l (t) 11 00 11 00 11 00 11 00 11
步骤12:异常磨斑图像的检测。分别利用梯度二值数据和步骤3中的磨痕区域的行长径和列长径值判断磨痕区域的形状是否正常:当存在任意整数m,n同时满足式(16-18)时,曲线的上升段被下降段分隔,曲线呈现上升、下降、上升和下降,这和标准现状曲线不符,说明行向形状异常,反之行向形状正常;同理,若存在任意整k,p同时满足式(19-21),说明列向形状异常;行向或列向为异常时则为异常磨斑图像,只有当行向和列向均为正常时,则认为是正常磨斑图像。正常的图像数据可用于后续的分析,而异常的图像数据不能用于后续的分析,需要重新试验。Step 12: Detection of abnormal wear scar images. Use the gradient binary data and the row length diameter and column length diameter of the wear scar area in step 3 to judge whether the shape of the wear scar area is normal: when there is any integer m, n satisfies formula (16-18), the curve The ascending segment of is separated by the descending segment, and the curve shows ascending, descending, ascending and descending, which is inconsistent with the standard current curve, indicating that the line shape is abnormal, otherwise the line shape is normal; similarly, if there is any integer k, p satisfies the formula at the same time (19-21), indicating that the column-wise shape is abnormal; when the row or column direction is abnormal, it is an abnormal wear spot image, and only when the row and column directions are both normal, it is regarded as a normal wear spot image. Normal image data can be used for subsequent analysis, while abnormal image data cannot be used for subsequent analysis and need to be re-tested.
Figure PCTCN2021077660-appb-000048
Figure PCTCN2021077660-appb-000048
Figure PCTCN2021077660-appb-000049
Figure PCTCN2021077660-appb-000049
Figure PCTCN2021077660-appb-000050
Figure PCTCN2021077660-appb-000050
Figure PCTCN2021077660-appb-000051
Figure PCTCN2021077660-appb-000051
Figure PCTCN2021077660-appb-000052
Figure PCTCN2021077660-appb-000052
Figure PCTCN2021077660-appb-000053
Figure PCTCN2021077660-appb-000053
其中,α为长度阈值,其取值为2-10。Among them, α is the length threshold, and its value is 2-10.
本实施例中,α取值为4;存在m=32,33;n=38,39,40,41,42,43满足式(16-18),说明行向形状异常;不存在k和l满足式(19-21),说明列向形状正常,故为异常磨斑图像,需要重新试验。In this embodiment, the value of α is 4; there is m=32, 33; n=38, 39, 40, 41, 42, 43 satisfies formula (16-18), indicating that the shape of the row is abnormal; there is no k and l Satisfies the formula (19-21), indicating that the column-direction shape is normal, so it is an abnormal wear spot image and needs to be tested again.
步骤13:方法结束。Step 13: The method ends.

Claims (10)

  1. 一种基于外形特征的异常磨斑图像的检测方法,其特征在于,包括如下步骤:A method for detecting abnormal wear scar images based on shape features is characterized in that it comprises the following steps:
    步骤1,对采集到的磨斑图像进行预处理,得到磨痕区域图;Step 1. Preprocess the collected wear scar image to obtain a wear scar area map;
    步骤2,在磨痕区域图中提取磨痕区域;Step 2: Extract the wear scar area in the wear scar area map;
    步骤3,计算磨痕区域的行长径和列长径;Step 3. Calculate the row and column diameters of the wear scar area;
    步骤4,根据步骤3中得到的磨痕水平图的行长径和列长径判断磨斑图像是否异常,若磨斑图像为异常时,算法结束;否则转入步骤5;Step 4. Determine whether the wear scar image is abnormal according to the row and column long diameters of the wear scar level map obtained in step 3. If the wear scar image is abnormal, the algorithm ends; otherwise, go to step 5;
    步骤5,对磨斑区域的行长度和列长度进行区间滤波,得到区间行长和区间列长;Step 5: Perform interval filtering on the row length and column length of the wear spot area to obtain the interval row length and the interval column length;
    步骤6,根据步骤5中得到的区间行长和区间列长计算梯度变化值;Step 6. Calculate the gradient change value according to the interval row length and interval column length obtained in step 5;
    步骤7,确定自适应分割阈值;Step 7. Determine the adaptive segmentation threshold;
    步骤8,根据步骤7中得到的自适应分割阈值对步骤6中得到的梯度变化值进行二值化处理,得到梯度二值数据;Step 8. Binarize the gradient change value obtained in step 6 according to the adaptive segmentation threshold obtained in step 7, to obtain gradient binary data;
    步骤10,根据步骤9得到的梯度二值数据和步骤3中的磨痕区域的行长径和列长径判断磨痕区域的形状是否正常。Step 10: Determine whether the shape of the wear scar area is normal according to the gradient binary data obtained in step 9 and the row major diameter and column major diameter of the wear scar area in step 3.
  2. 根据权利要求1所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于,步骤2中,在磨痕区域图中提取磨痕区域,具体方法是:The method for detecting abnormal wear scar images based on appearance features according to claim 1, characterized in that, in step 2, the wear scar area is extracted in the wear scar area map, and the specific method is:
    S1,确定步骤1中得到的磨痕区域图的磨痕方向角;S1, determine the wear scar direction angle of the wear scar area map obtained in step 1;
    S2,将磨痕区域图绕图像中心顺时针旋转磨痕方向角度,得到磨痕水平图;S2: Rotate the wear scar area map clockwise around the center of the image to obtain a wear scar level map;
    S3,在磨痕水平图中提取磨痕区域,其中,利用i 0和i 1表示磨痕水平图中磨痕区域的首行和尾行;利用j 0和j 1表示磨痕水平图中磨痕区域的首列和尾列: S3, extract the wear scar area in the wear scar level map, where i 0 and i 1 are used to represent the first and last rows of the wear scar area in the wear scar level graph; j 0 and j 1 are used to represent the wear scar in the wear scar level graph The first and last columns of the area:
    i 0<i且
    Figure PCTCN2021077660-appb-100001
    i 1>i且
    Figure PCTCN2021077660-appb-100002
    j 0<j且
    Figure PCTCN2021077660-appb-100003
    j 1>j且
    Figure PCTCN2021077660-appb-100004
    i 0 <i and
    Figure PCTCN2021077660-appb-100001
    i 1 >i and
    Figure PCTCN2021077660-appb-100002
    j 0 <j and
    Figure PCTCN2021077660-appb-100003
    j 1 >j and
    Figure PCTCN2021077660-appb-100004
  3. 根据权利要求1所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于, 步骤3中,根据下式计算磨痕区域的行长径:The method for detecting abnormal wear scar images based on shape features according to claim 1, characterized in that, in step 3, the line length diameter of the wear scar area is calculated according to the following formula:
    Figure PCTCN2021077660-appb-100005
    Figure PCTCN2021077660-appb-100005
    其中,r h为行长径;D h(i)为第i行的长度,
    Figure PCTCN2021077660-appb-100006
    i 0为磨痕区域的首行;i 1为磨痕区域的尾行;j 0为磨痕区域的首列;j 1为磨痕区域的尾列;
    Among them, rh is the length of the row; D h (i) is the length of the i-th row,
    Figure PCTCN2021077660-appb-100006
    i 0 is the first row of the wear scar area; i 1 is the last row of the wear scar area ; j 0 is the first row of the wear scar area; j 1 is the last row of the wear scar area;
    通过下式计算列长径:Calculate the column long diameter by the following formula:
    Figure PCTCN2021077660-appb-100007
    Figure PCTCN2021077660-appb-100007
    其中,r l为列长径;D l(j)为第j列的长度,
    Figure PCTCN2021077660-appb-100008
    Among them, r l is the length of the column; D l (j) is the length of the j-th column,
    Figure PCTCN2021077660-appb-100008
  4. 根据权利要求1所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于,步骤4中,根据步骤3中得到的磨痕水平图的行长径和列长径判断磨斑图像是否异常,具体方法是,若行长径和列长径满足下式,则认为磨斑图像异常:The method for detecting abnormal wear scar images based on appearance characteristics according to claim 1, characterized in that, in step 4, the wear scar is judged according to the row length diameter and the column length diameter of the wear scar level map obtained in step 3. Whether the image is abnormal, the specific method is, if the row length diameter and column length diameter meet the following formula, the wear spot image is considered abnormal:
    Figure PCTCN2021077660-appb-100009
    Figure PCTCN2021077660-appb-100009
    其中,β为差距阈值。Among them, β is the gap threshold.
  5. 根据权利要求1所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于,步骤5中,通过下式结合设定的区间长度对磨斑区域的行长度和列长度进行区间滤波:The method for detecting abnormal wear scar images based on appearance features according to claim 1, characterized in that, in step 5, the row length and column length of the wear scar area are divided by the following formula combined with the set interval length: Filtering:
    Figure PCTCN2021077660-appb-100010
    Figure PCTCN2021077660-appb-100010
    Figure PCTCN2021077660-appb-100011
    Figure PCTCN2021077660-appb-100011
    其中,s和t分别行区间和列区间的索引;L h(s)表示第s个区间行长;L l(t)为第t个区间列长;q为设定的区间长度;[ ]表示向下取整运算符。 Among them, s and t are the indexes of the row interval and column interval respectively; L h (s) represents the row length of the s-th interval; L l (t) is the column length of the t-th interval; q is the set interval length; [] Represents the round-down operator.
  6. 根据权利要求1所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于, 步骤6中,通过下式结构步骤5中得到的区间行长和区间列长计算梯度变化值:The method for detecting abnormal wear scar images based on appearance characteristics according to claim 1, wherein in step 6, the gradient change value is calculated by the interval row length and interval column length obtained in step 5 of the following structure:
    Figure PCTCN2021077660-appb-100012
    Figure PCTCN2021077660-appb-100012
    Figure PCTCN2021077660-appb-100013
    Figure PCTCN2021077660-appb-100013
    其中,I h(s)为s区间的行梯度,I l(t)为t区间的列梯度。 Among them, I h (s) is the row gradient in the s interval, and I l (t) is the column gradient in the t interval.
  7. 根据权利要求6所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于,步骤7中,确定自适应分割阈值,具体方法是:The method for detecting abnormal wear spot images based on shape features according to claim 6, characterized in that, in step 7, the adaptive segmentation threshold is determined, and the specific method is:
    将集合
    Figure PCTCN2021077660-appb-100014
    中的元素按升序排列,将位序
    Figure PCTCN2021077660-appb-100015
    位对应的元素作为行分割阈值T h
    Will gather
    Figure PCTCN2021077660-appb-100014
    The elements in are arranged in ascending order, and the order of
    Figure PCTCN2021077660-appb-100015
    The element corresponding to the bit is used as the row segmentation threshold Th ;
    将集合
    Figure PCTCN2021077660-appb-100016
    中元素升序排列后,将位序
    Figure PCTCN2021077660-appb-100017
    对应的元素值作为列分割阈值T 1
    Will gather
    Figure PCTCN2021077660-appb-100016
    After the elements are arranged in ascending order, the order of
    Figure PCTCN2021077660-appb-100017
    The corresponding element value is used as the column segmentation threshold T 1 ;
    其中,
    Figure PCTCN2021077660-appb-100018
    为四舍五入取整运算。
    in,
    Figure PCTCN2021077660-appb-100018
    It is a rounding operation.
  8. 根据权利要求1所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于,步骤8中,根据步骤7中得到的自适应分割阈值对步骤6中得到的梯度变化进行二值化,得到梯度二值数据,具体方法利用下式进行计算:The method for detecting abnormal wear scar images based on shape features according to claim 1, characterized in that, in step 8, the gradient change obtained in step 6 is binary-valued according to the adaptive segmentation threshold obtained in step 7 To obtain the gradient binary data, the specific method uses the following formula to calculate:
    Figure PCTCN2021077660-appb-100019
    Figure PCTCN2021077660-appb-100019
    Figure PCTCN2021077660-appb-100020
    Figure PCTCN2021077660-appb-100020
    其中,B h(s)为行梯度二值数据;B l(t)为列梯度二值数据。 Among them, B h (s) is the row gradient binary data; B l (t) is the column gradient binary data.
  9. 根据权利要求1所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于,步骤10中,根据步骤9得到的梯度二值数据和步骤3中的磨痕区域的行长径和列长径判断磨痕区域的形状是否正常,具体方法是:The method for detecting abnormal wear scar images based on shape features according to claim 1, characterized in that, in step 10, according to the gradient binary data obtained in step 9 and the line length diameter of the wear scar area in step 3 Determine whether the shape of the wear scar area is normal by using the long diameter and the column diameter. The specific method is:
    分别判断行向形状和列向形状是否均为异常,当行向形状和列向形状均正常时,则认定该磨斑图像为正常,否则,该磨斑图像为异常。Determine whether the row-direction shape and the column-direction shape are both abnormal. When the row-direction shape and the column-direction shape are both normal, the wear spot image is determined to be normal; otherwise, the wear spot image is abnormal.
  10. 根据权利要求9所述的一种基于外形特征的异常磨斑图像的检测方法,其特征在于,任意选取整数m,n,该整数同时满足下式时,则认为行向形状异常:The method for detecting abnormal wear scar images based on appearance features according to claim 9, characterized in that integers m and n are arbitrarily selected, and when the integers satisfy the following formula at the same time, it is considered that the shape of the row direction is abnormal:
    Figure PCTCN2021077660-appb-100021
    Figure PCTCN2021077660-appb-100021
    Figure PCTCN2021077660-appb-100022
    Figure PCTCN2021077660-appb-100022
    Figure PCTCN2021077660-appb-100023
    Figure PCTCN2021077660-appb-100023
    其中,α为长度阈值;Among them, α is the length threshold;
    任意选取整数k,p,该整数同时满足下式时,则认为行向形状异常:If the integers k and p are selected arbitrarily, when the integers satisfy the following formula at the same time, it is considered that the shape of the row direction is abnormal:
    Figure PCTCN2021077660-appb-100024
    Figure PCTCN2021077660-appb-100024
    Figure PCTCN2021077660-appb-100025
    Figure PCTCN2021077660-appb-100025
    Figure PCTCN2021077660-appb-100026
    Figure PCTCN2021077660-appb-100026
PCT/CN2021/077660 2020-03-20 2021-02-24 Method for detecting abnormal wear spot image on basis of appearance feature WO2021185035A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202010203325.5 2020-03-20
CN202010203325.5A CN111429428B (en) 2020-03-20 2020-03-20 Method for detecting abnormal grinding image based on appearance characteristics

Publications (1)

Publication Number Publication Date
WO2021185035A1 true WO2021185035A1 (en) 2021-09-23

Family

ID=71548346

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2021/077660 WO2021185035A1 (en) 2020-03-20 2021-02-24 Method for detecting abnormal wear spot image on basis of appearance feature

Country Status (2)

Country Link
CN (1) CN111429428B (en)
WO (1) WO2021185035A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117371993A (en) * 2023-12-05 2024-01-09 延安随缘科技发展有限公司 Method and system for improving service cycle of lubricating grease

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111429428B (en) * 2020-03-20 2022-05-13 长安大学 Method for detecting abnormal grinding image based on appearance characteristics

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107358604A (en) * 2017-07-20 2017-11-17 长安大学 A kind of four ball friction tests mill spot image abnormity detection method
CN107403445A (en) * 2017-07-20 2017-11-28 长安大学 A kind of pattern comparative analysis method for lubricant wear resistence
WO2020003762A1 (en) * 2018-06-28 2020-01-02 株式会社Screenホールディングス Shape measuring device and shape measuring method
CN111429428A (en) * 2020-03-20 2020-07-17 长安大学 Abnormal grinding image detection method based on appearance characteristics

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0318827D0 (en) * 2003-08-11 2003-09-10 Bg Intellectual Pty Ltd Dip value in seismic images
US8406859B2 (en) * 2008-08-10 2013-03-26 Board Of Regents, The University Of Texas System Digital light processing hyperspectral imaging apparatus
JP5529660B2 (en) * 2010-07-20 2014-06-25 パナソニック株式会社 Pupil detection device and pupil detection method
CN103712565B (en) * 2013-12-31 2016-05-25 长安大学 A kind of wear scar diameter measuring method based on steel ball polishing scratch gradient
IN201921008342A (en) * 2019-03-04 2019-03-15

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107358604A (en) * 2017-07-20 2017-11-17 长安大学 A kind of four ball friction tests mill spot image abnormity detection method
CN107403445A (en) * 2017-07-20 2017-11-28 长安大学 A kind of pattern comparative analysis method for lubricant wear resistence
WO2020003762A1 (en) * 2018-06-28 2020-01-02 株式会社Screenホールディングス Shape measuring device and shape measuring method
CN111429428A (en) * 2020-03-20 2020-07-17 长安大学 Abnormal grinding image detection method based on appearance characteristics

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117371993A (en) * 2023-12-05 2024-01-09 延安随缘科技发展有限公司 Method and system for improving service cycle of lubricating grease
CN117371993B (en) * 2023-12-05 2024-03-08 延安随缘科技发展有限公司 Method and system for improving service cycle of lubricating grease

Also Published As

Publication number Publication date
CN111429428B (en) 2022-05-13
CN111429428A (en) 2020-07-17

Similar Documents

Publication Publication Date Title
WO2021185035A1 (en) Method for detecting abnormal wear spot image on basis of appearance feature
CN109816031B (en) Transformer state evaluation clustering analysis method based on data imbalance measurement
CN109115800B (en) Method for rapidly detecting burrs of product and accurately measuring length
CN113838054B (en) Mechanical part surface damage detection method based on artificial intelligence
CN105445306A (en) Method for evaluating element segregation degree in steel
CN113340598B (en) Rolling bearing intelligent fault diagnosis method based on regularized sparse model
CN111445438B (en) Double-neighborhood grinding crack angle automatic detection method
CN110738346A (en) batch electric energy meter reliability prediction method based on Weibull distribution
CN115330799B (en) Automatic fault diagnosis method for instrument
CN112131797B (en) Main shaft bearing life prediction and reliability assessment method based on stress analysis
CN109712141A (en) Continuous casting steel billet center segregation fault detection quantifies ranking method
WO2021169996A1 (en) Grinding mark angle automatic detection method based on grinding mark gray level similarity
CN103712565B (en) A kind of wear scar diameter measuring method based on steel ball polishing scratch gradient
CN110749443B (en) Rolling bearing fault diagnosis method and system based on high-order origin moment
CN117218124A (en) Copper-aluminum switching surface processing quality detection method
CN103472070B (en) Method for determining detergency and dispersibility of engine oil
CN106841575A (en) A kind of four ball friction tests mill spot image polishing scratch direction automatic positioning method
CN107358604A (en) A kind of four ball friction tests mill spot image abnormity detection method
CN116644368A (en) Outlier identification method based on improved Grabbs test method
CN114037215A (en) Fault severity evaluation method of variable-step multi-scale complexity fusion index
CN114972353A (en) Quality evaluation method and system for metal high-pressure special pipe fitting
CN107403445B (en) Morphology comparison analysis method for wear resistance of lubricant
CN113588259B (en) Equipment vibration signal scale curve turning point detection method and working condition monitoring device
CN114782303A (en) Method for detecting abnormal grinding deformation degree based on boundary distribution characteristics
CN117491422B (en) Method for detecting high heat-conducting property of aluminum alloy material

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21772174

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21772174

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 21772174

Country of ref document: EP

Kind code of ref document: A1