CN102519990B - Fault on-line detection method of steel wire rope conveyer belt based on texture regularity analysis - Google Patents

Fault on-line detection method of steel wire rope conveyer belt based on texture regularity analysis Download PDF

Info

Publication number
CN102519990B
CN102519990B CN201110396031.XA CN201110396031A CN102519990B CN 102519990 B CN102519990 B CN 102519990B CN 201110396031 A CN201110396031 A CN 201110396031A CN 102519990 B CN102519990 B CN 102519990B
Authority
CN
China
Prior art keywords
image
lrb
value
fault
mrow
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
CN201110396031.XA
Other languages
Chinese (zh)
Other versions
CN102519990A (en
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.)
TIANJIN HENG YI MECHANICAL AND ELECTRONIC TECHNOLOGY Co Ltd
Tianjin Polytechnic University
Taiyuan University of Technology
Original Assignee
TIANJIN HENG YI MECHANICAL AND ELECTRONIC TECHNOLOGY Co Ltd
Tianjin Polytechnic University
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 TIANJIN HENG YI MECHANICAL AND ELECTRONIC TECHNOLOGY Co Ltd, Tianjin Polytechnic University filed Critical TIANJIN HENG YI MECHANICAL AND ELECTRONIC TECHNOLOGY Co Ltd
Priority to CN201110396031.XA priority Critical patent/CN102519990B/en
Publication of CN102519990A publication Critical patent/CN102519990A/en
Application granted granted Critical
Publication of CN102519990B publication Critical patent/CN102519990B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The invention relates to a fault on-line detection method of a steel wire rope conveyer belt based on texture regularity analysis, belonging to the field of fault detection. The problems that in the current fault detection of the steel wire rope conveyer belt based on X-ray, the real-time performance is worse, the accuracy rate is not high, the fault detection is mainly done basing on manual analysis, and the like are solved. According to the characteristics that the gray scale of an X-ray image of the steel wire rope conveyer belt is in regular variation and the texture of the X-ray image of the steel wire rope conveyer belt is in regular strike-shaped distribution, after the input image is subjected to normalization and standardization processes, parameter values of LRB (Light Regular Band) and DRB (Dark Regular Band) of each pixel can be progressively scanned and calculated, the parameter values of the LRB and the DRB of each pixel can be subtracted from the parameter values of the LRB and the DRB of the whole image, then the difference value can be compared with a threshold value obtained in a training stage, thus whether the fault exists or not can be judged, and a fault area can be detected from the fault image. According to the fault on-line detection method of the steel wire rope conveyer belt based on the texture regularity analysis, disclosed by the invention, the fault detection of the steel wire rope conveyer belt can be online, automatically and accurately realized.

Description

On-line detection method for faults of steel wire rope core conveying belt based on texture regularity analysis
Technical Field
The invention belongs to the field of fault detection, and particularly relates to a method for obtaining a threshold value and detecting a fault on line based on texture regularity analysis in fault detection of a steel wire rope core conveying belt.
Background
The steel wire rope core conveying belt is used in a large number in coal mine production, faults occur due to various reasons in use, if the faults cannot be detected and processed in time, safety accidents such as belt breakage or longitudinal tearing can be caused, and safety production is seriously influenced[1-3]. Therefore, each coal mine sets up a corresponding steel wire rope core conveying belt detection system.
At present, researchers have proposed effective methods such as eddy current detection, magnetic induction detection, shock wave detection, X-ray detection and the like besides manual detection for detection of steel wire rope core conveying belts, and have obtained certain application[4-9]. The X-ray detection method can acquire and display the internal steel wire rope core image of the steel wire rope core conveying belt in real time, can detect faults by using an image processing technology, and is simple to operate and high in detection speed, so that the detection device based on the method is valued by researchers and popular with conveying belt users. The inventor develops a nondestructive testing system of a steel wire rope core conveying belt based on X-ray in earlier stage, and researches an automatic detection method of the tensile failure of a joint of the steel wire rope core conveying belt[9,10]. In the use of the detection system, the user urgently requires that the faults except the joint stretching fault can be automatically detected on line, such as the faults of the conveyer belt breakage, the disconnection of a single steel wire rope core, the intersection of the ends of the steel wire rope cores and the like.
The X-ray image of the steel wire rope core conveying belt has strong periodic texture characteristics, and the texture is distributed in a regular stripe shape. In the aspect of regularity analysis of texture detection, although there are several effective methods, most of them are very sensitive to noise, the accuracy of detection results is not high, and in addition, the computation amount is often large, and therefore, the method is not suitable for on-line rapid detection. The invention can realize the fault detection of the steel wire rope core conveying belt on line, automatically and accurately.
Reference documents:
[1] mine steel rope core belt conveyor real-time working condition monitoring and fault diagnosis technology [ J ] coal science report, 2005, 30 (2): 245-250.
[2] Analysis of the cause of longitudinal tearing of a mining conveyor belt and prevention thereof [ J ] coal mine machinery, 2007, 28 (10): 182-183.
[3] Longitudinal tear monitoring method research [ J ] of university of mineral china, 2002, 31 (1): 49-52.
[4]Harrison A.A new technique for measuring loss of adhesion in conveyor belt splices[J].Australian J.Coal Mining Technology and Research,1983,(4):27-34.
[5]Sukhorukov V.Steel-cord conveyor belt NDT[J].International Journal of Materials&ProductTechnology,2006,27(3/4):236-246.
[6]Alport M,Basson J F,Padayachee T.Digital magnetic imaging of steel cord belts[J].BulkSolids Handling,2008,28(3):182-185.
[7] Design and implementation of a longitudinal anti-tear system of a heavy conveyor belt [ J ] mining machinery, 2006, 34 (6): 93-95.
[8]He Zhi-qiang,Wang Bai-yan,Gao Yu-lin.Research of algorithm for judgment of joint elongationof conveyer Belt[J].Journal of China University of Mining&Technology,1999,9(1):51-54.
[9]Li Xian-Guo,Miao Chang-Yun,Zhang Yan.An algorithm for selecting and stitching theconveyer belt joint images based on X-ray[A].2010International Conference on IntelligentComputation Technology and Automation[C].IEEE Computer Society,2010.vol.1,474-477.
[10] The method for automatically detecting the tensile failure of the joint of the steel wire rope core conveying belt by X-ray imaging comprises the following steps: china, 201010146881.X [ P ].2010-4-15.
Disclosure of Invention
The invention aims to provide an X-ray imaging on-line detection method for faults of a steel wire rope core conveying belt, which solves the problems that the current X-ray-based fault detection of the steel wire rope core conveying belt is poor in real-time performance, low in accuracy rate, mainly dependent on manual analysis and the like.
The technical scheme adopted by the invention is divided into two stages, namely a training stage and a detection stage.
The specific operation steps in the training phase are as follows:
step 1, reading a training image without fault information.
And 2, normalizing the input image, and enabling the pixel values of the image to be in the range of (0, 1) through linear transformation.
Step 3, normalizing the normalized image, and subtracting the pixel mean value of the normalized image
Figure BSA00000627313000021
All the pixel values are at
Figure BSA00000627313000022
Within the range.
Step 4, line scanning, calculating LRB (Light regulated Band, bright ruled Band, which measures bright defects by using regularity of texture image) and DRB (Dark regulated Band, which measures Dark defects by using regularity of texture image) parameter values of each pixel point of each line, wherein when the LRB parameter values are calculated on a certain line of the image, the definition formula is as follows
When calculating DRB parameter value, the formula is
Wherein the moving average
Figure BSA00000627313000025
And standard deviation of
Figure BSA00000627313000026
Is expressed as
Figure BSA00000627313000027
Figure BSA00000627313000028
Where n is an integer representing the dimension of the repeated elements in the input image, i.e. the texture period, xijIs the gray value of pixel points in the ith row and the jth column after the image X with the size of p multiplied by q is normalized, and r is more than or equal to 01≤rn≤q-1,rn=r1+n-1,r1∈[0,q-n],rn∈[n-1,q-1]And is and
Figure BSA00000627313000029
Figure BSA000006273130000211
Figure BSA000006273130000212
and 5, scanning lines and calculating the mean value of LRB parameter values of each line. For example, for a training image of size p × q, the mean of the LRB parameter values for row i is
Figure BSA00000627313000031
Wherein
Figure BSA00000627313000032
The LRB parameter values of the (i, j) pixel points obtained by calculation in the row direction are used. And calculating the mean value of the DRB values of each row in the same way.
Step 6, line scanning, calculating two thresholds of LRB parameterAnd
Figure BSA00000627313000034
and two thresholds for DRB parameters
Figure BSA00000627313000035
And
Figure BSA00000627313000036
during calculation, the maximum deviation amount of the mean value of the LRB parameter values and the DRB parameter values corresponding to the positive deviation and the negative deviation of the LRB parameter values and the DRB parameter values of the training images is used as the reference of a threshold value, and an adjusting parameter is added for self-adaptive adjustment and used as the threshold value. For example, for a training image of size p × q, the LRB parameter values for row i are biased forward from the maximum value of the mean of the LRB parameter values for that row
l max + i = max ( L r i , j - u L , r i ) = max ( L r i , j ) - u L , r i
The forward deviation of the LRB parameter values obtained by all the training rows is as the maximum value
l max + = max ( l max + i )
In practice, the maximum value of the forward deviation of the LRB parameter value used is
<math> <mrow> <msubsup> <mi>l</mi> <mrow> <mi>max</mi> <mo>+</mo> </mrow> <mi>use</mi> </msubsup> <mo>=</mo> <msub> <mi>l</mi> <mrow> <mi>max</mi> <mo>+</mo> </mrow> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </math>
Wherein alpha is an adjusting parameter, alpha is more than or equal to 0 and is defaulted to be 1, and the maximum value of the negative deviation of the LRB parameter value is
l max - i = min ( L r i , j - u L , r i ) = min ( L r i , j ) - u L , r i
l max - = min ( l max - i )
<math> <mrow> <msubsup> <mi>l</mi> <mrow> <mi>max</mi> <mo>-</mo> </mrow> <mi>use</mi> </msubsup> <mo>=</mo> <msub> <mi>l</mi> <mrow> <mi>max</mi> <mo>-</mo> </mrow> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </math>
For obtaining DRB parameter values in the same wayTwo thresholds, i.e. maximum of the DRB parameter value in the forward direction
Figure BSA000006273130000313
And maximum value of negative deviation of DRB parameter value
Figure BSA000006273130000314
In the training process, the mean value of the LRB and DRB parameter values is independently calculated for each row, and the maximum value of the mean value of the corresponding row with positive deviation and negative deviation is calculated respectively, so that the method is equivalent to training with all rows respectively, which is equivalent to increasing the number of training samples and reducing the accidental influence of training sample selection.
The specific operation steps in the detection stage are as follows:
step 1, reading an image to be detected.
And 2, normalizing the image read in the step 1, and enabling each pixel value of the image to be in a (0, 1) range through linear transformation.
Step 3, normalizing the normalized image obtained in the step 2, and subtracting the pixel mean value of the normalized image
Figure BSA00000627313000041
All the pixel values are at
Figure BSA00000627313000042
Within the range.
In particular, the pixel mean subtracted here
Figure BSA00000627313000043
Is the normalized pixel mean value of the image to be detected
Figure BSA00000627313000044
So as to adapt to different changes of the image to be detected and keep the image to be detected and the training image independent in operation.
And 4, scanning the normalized image obtained in the step 3, and calculating LRB and DRB parameter values of each pixel point in each line.
And 5, calculating the mean value of the LRB and DRB parameter values of each row according to the result obtained in the step 4. Because the difference between the mean values of the LRB parameter values and the DRB parameter values in the X-ray image of the non-joint area of the steel wire rope core conveying belt is not large, and the size of the fault area is much smaller than that of the whole image, in order to reduce the influence of the LRB or DRB parameter values of the fault area on the line where the fault is positioned on the mean value of the line, the mean value of the mean values of the LRB and DRB parameter values of all lines of the detected image is used as the mean value for actually judging all lines, and for the image to be detected with the size of p multiplied by q, the mean value for actually judging and using the LRB is used as the mean
Figure BSA00000627313000045
And mean value used for DRB actual discrimination
Figure BSA00000627313000046
Are respectively calculated as
<math> <mrow> <msubsup> <mi>u</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>use</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>u</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>i</mi> </msubsup> </mrow> <mi>p</mi> </mfrac> </mrow> </math>
<math> <mrow> <msubsup> <mi>u</mi> <mrow> <mi>D</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>use</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>u</mi> <mrow> <mi>D</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>i</mi> </msubsup> </mrow> <mi>p</mi> </mfrac> </mrow> </math>
Wherein,
Figure BSA00000627313000049
is the mean of the LRB parameter values of row i;
Figure BSA000006273130000410
is the mean of the DRB parameter values of row i.
Step 6, using the result of step 5 and two thresholds of LRB parameter value obtained in training phase And two thresholds for DRB parameter values
Figure BSA000006273130000413
Line scanning, performing threshold segmentation on the original image, if the (i, j) pixel point <math> <mrow> <msubsup> <mi>L</mi> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>u</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>use</mi> </msubsup> <mo>&NotElement;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>l</mi> <mrow> <mi>max</mi> <mo>-</mo> </mrow> <mi>use</mi> </msubsup> <mo>,</mo> </mtd> <mtd> <msubsup> <mi>l</mi> <mrow> <mi>max</mi> <mo>+</mo> </mrow> <mi>use</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Or <math> <mrow> <msubsup> <mi>D</mi> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>u</mi> <mrow> <mi>D</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>use</mi> </msubsup> <mo>&NotElement;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>d</mi> <mrow> <mi>max</mi> <mo>-</mo> </mrow> <mi>use</mi> </msubsup> <mo>,</mo> </mtd> <mtd> <msubsup> <mi>d</mi> <mrow> <mi>max</mi> <mo>+</mo> </mrow> <mi>use</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> </math> Making the gray value of the pixel point be white, otherwise setting the gray value to be black, whereinIn lineThe LRB parameter values of the (i, j) pixel points obtained by calculation in the direction,
Figure BSA000006273130000417
and calculating the DRB parameter value of the (i, j) pixel point in the row direction.
And 7, counting white pixel points according to the result of thresholding division in the step 6, and when the number of the white pixel points exceeds the set number, considering the image as an image containing fault information.
And 8, reporting the detection result according to the judgment of the step 7.
The invention has the advantages that the practical application shows that the invention can automatically and accurately realize the on-line detection of the faults of the steel wire rope core conveyer belt breakage, the single steel wire rope core breakage, the steel wire rope core end crossing and the like, and has good use value in industrial application. The technology of the invention has the following advantages:
(1) the method has the advantages of high automation degree, simple use and few parameters needing to be set, and the detection result is not very sensitive to the value of n (n is the texture period of the input image).
(2) In the training process, each row and each column are processed independently, which is equivalent to increase the number of training samples, reduces the accidental influence of the selection of the training samples, and has high accuracy of fault detection.
(3) The method can detect fine faults, has high image contrast, can well highlight the fault area, and the detected fault area is closer to the appearance of the fault, thereby being beneficial to further segmentation and fault identification of the fault area.
(4) The method has small calculation amount and is suitable for the on-line detection of the fault of the steel wire rope core conveying belt.
(5) The texture regularity analysis-based method can also be applied to the online detection of other articles with regular texture characteristics, such as fabrics, wallpaper, ceramics and the like.
Drawings
FIG. 1 is an operation flow chart of a fault on-line detection training phase of a steel wire rope core conveyer belt;
FIG. 2 is an operation flow chart of a fault on-line detection stage of the steel wire rope core conveying belt;
FIG. 3 is a fault image of a wire rope core conveyor belt;
fig. 4 is a detection result image of the failure image of the wire rope core conveyor belt shown in fig. 3.
Detailed Description
The following detailed description of the embodiments of the invention refers to the accompanying drawings and examples.
The method provided by the invention is divided into two stages when being implemented specifically, namely, firstly, a threshold value required for discrimination is acquired through a training stage, and an operation flow chart of the on-line detection training stage of the fault of the steel wire rope core conveying belt is shown as an attached figure 1; then, the detection of the image to be detected is realized in the detection stage, and the operation flow chart of the on-line detection stage of the fault of the steel wire rope core conveying belt is shown in the attached figure 2.
The steel wire rope cores in the steel wire rope core conveying belt are longitudinally arranged according to a certain distance, and the X-ray image has strong texture characteristics, namely the gray scale of the image is regularly changed, and the texture is regularly distributed in a stripe shape (vertical stripe shape). In order to reduce the complexity of implementation and improve the execution efficiency, the fault on-line detection of the steel wire rope core conveying belt only needs to be carried out in the row direction in the actual operation.
The method is used for detecting three typical faults of common conveyor belt breakage of the steel wire rope core conveyor belt, single steel wire rope core breakage, steel wire rope core end intersection and the like. The detection database comprises 2577X-ray images of four steel wire rope core conveyor belts with large differences of four coal mines, wherein the X-ray images comprise 8 conveyor belt fracture images (wherein a certain conveyor belt fracture image is shown in a figure 3 (a)), 33 single steel wire rope core fracture images (wherein a certain single steel wire rope core fracture image is shown in a figure 3 (b)), and 22 steel wire rope core end crossing images (wherein a certain steel wire rope core end crossing image is shown in a figure 3 (c)). In order to enhance the display effect, the image shown in fig. 3 is a partial image cut from the original complete image and containing the fault region with a size of 160 × 160 pixels. And when each conveying belt is detected, taking the first image without fault information as a training image. n is 12, alpha is a default value of 1, and the statistical threshold of the white pixel point is 5.
The operation flow of the training phase is as follows:
step 1, reading a first training image without fault information in a database of a current conveying belt.
And 2, normalizing the input image, and enabling the pixel values of the image to be in the range of (0, 1) through linear transformation.
Step 3, normalizing the normalized image, and subtracting the pixel mean value of the normalized image
Figure BSA00000627313000051
All the pixel values are at
Figure BSA00000627313000052
Within the range.
Step 4, line scanning, calculating LRB (Light regulated Band, bright ruled Band, which measures bright defects by using regularity of texture image) and DRB (Dark regulated Band, which measures Dark defects by using regularity of texture image) parameter values of each pixel point of each line, wherein when the LRB parameter values are calculated on a certain line of the image, the definition formula is as follows
When calculating DRB parameter value, the formula is
Figure BSA00000627313000062
Wherein the moving averageAnd standard deviation of
Figure BSA00000627313000064
Is expressed as
Figure BSA00000627313000065
Figure BSA00000627313000066
Where n is an integer representing the dimension of the repeated elements in the input image, i.e. the texture period, xijIs the gray value of pixel points in the ith row and the jth column after the image X with the size of p multiplied by q is normalized, and r is more than or equal to 01≤rn≤q-1,rn=r1+n-1,r1∈[0,q-n],rn∈[n-1,q-1]And is and
Figure BSA00000627313000067
Figure BSA00000627313000068
Figure BSA00000627313000069
Figure BSA000006273130000610
and 5, scanning lines and calculating the mean value of LRB parameter values of each line. For example, for a training image of size p × q, the mean of the LRB parameter values for row i is
Figure BSA000006273130000611
Wherein
Figure BSA000006273130000612
The LRB parameter values of the (i, j) pixel points obtained by calculation in the row direction are used. And calculating the mean value of the DRB values of each row in the same way.
Step 6, line scanning, calculating two thresholds of LRB parameter
Figure BSA000006273130000613
And
Figure BSA000006273130000614
and two thresholds for DRB parameters
Figure BSA000006273130000615
And
Figure BSA000006273130000616
for a training image of size p × q, the LRB parameter values for row i deviate positively from the maximum value of the mean of the LRB parameter values for that row
l max + i = max ( L r i , j - u L , r i ) = max ( L r i , j ) - u L , r i
The forward deviation of the LRB parameter values obtained by all the training rows is as the maximum value
l max + = max ( l max + i )
In practice, the maximum value of the forward deviation of the LRB parameter value used is
<math> <mrow> <msubsup> <mi>l</mi> <mrow> <mi>max</mi> <mo>+</mo> </mrow> <mi>use</mi> </msubsup> <mo>=</mo> <msub> <mi>l</mi> <mrow> <mi>max</mi> <mo>+</mo> </mrow> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </math>
Wherein, alpha is an adjusting parameter, alpha is 1, and the maximum value of the negative deviation of the LRB parameter value is
l max - i = min ( L r i , j - u L , r i ) = min ( L r i , j ) - u L , r i
l max - = min ( l max - i )
<math> <mrow> <msubsup> <mi>l</mi> <mrow> <mi>max</mi> <mo>-</mo> </mrow> <mi>use</mi> </msubsup> <mo>=</mo> <msub> <mi>l</mi> <mrow> <mi>max</mi> <mo>-</mo> </mrow> </msub> <mo>&times;</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </math>
Two threshold values of DRB parameter value, namely the maximum value of the positive deviation of DRB parameter value
Figure BSA00000627313000075
And maximum value of negative deviation of DRB parameter value
Figure BSA00000627313000076
Four thresholds are established by the training phase
Figure BSA00000627313000077
And then, the detection stage can be entered, and the fault on-line detection of the steel wire rope core conveying belt is realized.
The operation flow of the detection stage is as follows:
step 1, reading an image to be detected, wherein the image to be detected may be a fault image or a non-fault image.
And 2, normalizing the image read in the step 1, and enabling each pixel value of the image to be in a (0, 1) range through linear transformation.
Step 3, normalizing the normalized image obtained in the step 2, and subtracting the pixel mean value of the normalized image
Figure BSA00000627313000078
All the pixel values are at
Figure BSA00000627313000079
Within the range.
And 4, scanning the normalized image obtained in the step 3, and calculating LRB and DRB parameter values of each pixel point in each line.
And 5, calculating the mean value of the LRB and DRB parameter values of each row according to the result obtained in the step 4. For the image to be detected with the size of p multiplied by q, the average value used by LRB actual discrimination
Figure BSA000006273130000710
And mean value used for DRB actual discrimination
Figure BSA000006273130000711
Are respectively calculated as
<math> <mrow> <msubsup> <mi>u</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>use</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>u</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>i</mi> </msubsup> </mrow> <mi>p</mi> </mfrac> </mrow> </math>
<math> <mrow> <msubsup> <mi>u</mi> <mrow> <mi>D</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>use</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>u</mi> <mrow> <mi>D</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>i</mi> </msubsup> </mrow> <mi>p</mi> </mfrac> </mrow> </math>
Wherein,
Figure BSA000006273130000714
is the mean of the LRB parameter values of row i;
Figure BSA000006273130000715
is the mean of the DRB parameter values of row i.
Step 6, using the result of step 5 and two thresholds of LRB parameter value obtained in training phase
Figure BSA000006273130000716
And two thresholds for DRB parameter valuesLine scanning, performing threshold segmentation on the original image, if the (i, j) pixel point <math> <mrow> <msubsup> <mi>L</mi> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>u</mi> <mrow> <mi>L</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>use</mi> </msubsup> <mo>&NotElement;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>l</mi> <mrow> <mi>max</mi> <mo>-</mo> </mrow> <mi>use</mi> </msubsup> <mo>,</mo> </mtd> <mtd> <msubsup> <mi>l</mi> <mrow> <mi>max</mi> <mo>+</mo> </mrow> <mi>use</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> Or <math> <mrow> <msubsup> <mi>D</mi> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>u</mi> <mrow> <mi>D</mi> <mo>,</mo> <mi>r</mi> </mrow> <mi>use</mi> </msubsup> <mo>&NotElement;</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>d</mi> <mrow> <mi>max</mi> <mo>-</mo> </mrow> <mi>use</mi> </msubsup> <mo>,</mo> </mtd> <mtd> <msubsup> <mi>d</mi> <mrow> <mi>max</mi> <mo>+</mo> </mrow> <mi>use</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> </math> Making the gray value of the pixel point be white, otherwise setting the gray value to be black, wherein
Figure BSA00000627313000083
For the LRB parameter values of the (i, j) pixel points calculated in the row direction,
Figure BSA00000627313000084
and calculating the DRB parameter value of the (i, j) pixel point in the row direction.
And 7, counting white pixel points according to the result of thresholding and dividing in the step 6, and when the number of the white pixel points exceeds a set threshold value, determining that the image is an image containing fault information.
And 8, reporting the detection result according to the judgment of the step 7.
Fig. 4 is a detection result image of the fault image of the steel wire rope core conveyer belt shown in fig. 3 (in order to enhance the display effect, fig. 4 is a local image cut from the result image and corresponding to the fault area shown in fig. 3 and having a size of 160 × 160 pixels), from which it can be seen that the method of the present invention can detect a fine fault (as shown in fig. 3(b), that is, a single steel wire rope core disconnection fault caused by corrosion or manufacturing process, etc., because the disconnection degree is low, the gray value has no large change, and the position thereof is close to the edge of the image, it is generally difficult for human eyes to identify it, but the method of the present invention can still correctly identify it), the image contrast is high, and the fault area can be well highlighted, and the detected fault area is closer to the appearance of the fault, thereby being beneficial to further segmentation and fault identification of the fault area.
For four conveyer belt images of the detection database, the three types of faults contained in the images are successfully detected, and the accuracy rate is 100%.
The detection stage program provided by the invention is written by C # language, and runs on a PC with Intel Core 2Duo CPU T58702.0GHz main frequency and 1GB memory, the detection time of an X-ray image of a 600X 960 steel wire Core conveying belt is 630ms, which is less than the acquisition time of 750ms, so the method is suitable for online detection. The method can also be applied to the online detection of other articles with regular texture characteristics, such as fabrics, wallpaper, ceramics and the like.
Practical application shows that the invention can realize the on-line, automatic and accurate detection of faults of the steel wire rope core conveyer belt breakage, single steel wire rope core breakage, steel wire rope core end crossing and the like, and has good use value in industrial application.

Claims (1)

1. A fault on-line detection method of a steel wire rope core conveying belt based on texture regularity analysis is characterized by comprising the following steps:
step 1, reading an image to be detected;
step 2, normalizing the image read in the step 1, and enabling each pixel value of the image to be in a (0, 1) range through linear transformation;
step 3, normalizing the normalized image obtained in the step 2, and subtracting the normalized pixel of the image to be detected from the normalized imageMean value
Figure FSB00001175213900000117
All the pixel values are at
Figure FSB00001175213900000118
Within the range;
step 4, scanning the normalized image obtained in the step 3, and calculating LRB and DRB parameter values of each pixel point of each line, wherein LRB is Light regulated Band initial letter abbreviation, meaning is bright Regular Band, bright defect is measured by utilizing regularity of texture image, DRB is DarkRegularBand initial letter abbreviation, meaning is dark Regular Band, dark defect is measured by utilizing regularity of texture image, and when LRB parameter values are calculated on a certain line of the image, a defined formula is
When calculating DRB parameter value, the formula is
Figure FSB0000117521390000012
Wherein the moving average
Figure FSB0000117521390000013
And standard deviation of
Figure FSB0000117521390000014
Is expressed as
Figure FSB0000117521390000015
Figure FSB0000117521390000016
Where n is an integer representing the dimension of the repeated elements in the input image, i.e. the texture period, xijIs the gray value of pixel points in the ith row and the jth column after the image X with the size of p multiplied by q is normalized, and r is more than or equal to 01≤rn≤q-1,rn=r1+n-1,r1∈[0,q-n],rn∈[n-1,q-1]And is and
Figure FSB0000117521390000017
Figure FSB0000117521390000018
Figure FSB0000117521390000019
Figure FSB00001175213900000110
step 5, calculating the mean value of the LRB and DRB parameter values of each line according to the result obtained in the step 4, and for the image to be detected with the size of p multiplied by q, actually judging the mean value used by the LRBAnd mean value used for DRB actual discrimination
Figure FSB00001175213900000112
Are respectively calculated as
Figure FSB00001175213900000113
Figure FSB00001175213900000114
Wherein,
Figure FSB00001175213900000115
is the mean value of the LRB parameter values of the ith row, and the calculation formula is
In the above formula
Figure FSB0000117521390000021
For the LRB parameter values of the (i, j) pixel points calculated in the row direction,is the mean value of the DRB parameter values of the ith row, and the calculation formula is
Figure FSB0000117521390000023
In the above formula
Figure FSB0000117521390000024
Calculating the DRB parameter value of the (i, j) pixel point in the row direction;
step 6, utilizing the result of step 5 and the maximum value of the negative deviation of the LRB parameter value obtained in the training stage
Figure FSB0000117521390000025
Maximum value of forward deviation of LRB parameter value
Figure FSB0000117521390000026
Maximum value of negative deviation of DRB parameter valueAnd maximum value of the forward deviation of the DRB parameter value
Figure FSB0000117521390000028
Performing line scanning, performing threshold segmentation on the original image, and determining the (i, j) pixel point
Figure FSB00001175213900000210
Or
Figure FSB00001175213900000211
Making the gray value of the pixel point be white, otherwise setting the gray value to be black, wherein
Figure FSB00001175213900000212
Figure FSB00001175213900000213
The method comprises selecting a fault-free image as a training image, and performing the same operations as in steps 2 to 4, wherein for the training image with the size of p × q, the maximum value of the forward deviation of the LRB parameter value of the ith row from the mean value of the LRB parameter values of the row is
The forward deviation of the LRB parameter values obtained by all the training rows is as the maximum value
Figure FSB00001175213900000215
In practice, the maximum value of the forward deviation of the LRB parameter value used is
Figure FSB00001175213900000216
Wherein alpha is an adjusting parameter, alpha is more than or equal to 0 and is defaulted to be 1, and the maximum value of the negative deviation of the LRB parameter value is
Figure FSB00001175213900000217
Obtaining the maximum value of the negative deviation of the DRB parameter value in the same way
Figure FSB00001175213900000218
And maximum value of the forward deviation of the DRB parameter value
Figure FSB00001175213900000219
Step 7, carrying out white pixel point statistics according to the result of thresholding division in the step 6, and when the number of the white pixel points exceeds the set number, considering the image as an image containing fault information;
and 8, reporting the detection result according to the judgment of the step 7.
CN201110396031.XA 2011-12-05 2011-12-05 Fault on-line detection method of steel wire rope conveyer belt based on texture regularity analysis Active CN102519990B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110396031.XA CN102519990B (en) 2011-12-05 2011-12-05 Fault on-line detection method of steel wire rope conveyer belt based on texture regularity analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110396031.XA CN102519990B (en) 2011-12-05 2011-12-05 Fault on-line detection method of steel wire rope conveyer belt based on texture regularity analysis

Publications (2)

Publication Number Publication Date
CN102519990A CN102519990A (en) 2012-06-27
CN102519990B true CN102519990B (en) 2014-03-26

Family

ID=46290980

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110396031.XA Active CN102519990B (en) 2011-12-05 2011-12-05 Fault on-line detection method of steel wire rope conveyer belt based on texture regularity analysis

Country Status (1)

Country Link
CN (1) CN102519990B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102910436A (en) * 2012-10-26 2013-02-06 昆山特力伯传动科技有限公司 Automatic detection method for conveying belt
CN104016099B (en) * 2014-05-23 2016-04-13 爱德森(厦门)电子有限公司 The damage of steel cord conveyor belt longitudinal tear is in line vortex monitoring method
CN104535356B (en) * 2015-01-19 2017-08-29 中南大学 A kind of roller steel wire rope rope fault monitoring method and system based on machine vision
CN106706238B (en) * 2016-11-22 2019-01-29 山西大学 Steel cable core conveying belt joint overlap joint label and recognition methods
CN107423744A (en) * 2017-03-23 2017-12-01 北京环境特性研究所 The Seam tracking and damage positioning method of steel rope core conveying belt
CN108169332A (en) * 2017-12-14 2018-06-15 太原理工大学 Steel cord detection method and device based on thermal acoustic imaging
JP7134793B2 (en) 2018-08-29 2022-09-12 オーチス エレベータ カンパニー Elevator rope elongation measuring device and elevator rope elongation measuring method
CN111901559B (en) * 2020-07-08 2021-06-01 南京幸庄科技创新产业园管理有限公司 System and method for recognizing breaking state of protective rope body

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1910448A (en) * 2004-06-11 2007-02-07 株式会社石田 X-ray inspection apparatus
CN101986143A (en) * 2010-03-17 2011-03-16 燕山大学 Machine vision belt tear detection and protective device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7564943B2 (en) * 2004-03-01 2009-07-21 Spectramet, Llc Method and apparatus for sorting materials according to relative composition

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1910448A (en) * 2004-06-11 2007-02-07 株式会社石田 X-ray inspection apparatus
CN101986143A (en) * 2010-03-17 2011-03-16 燕山大学 Machine vision belt tear detection and protective device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Henry Y. T. Ngan ET AL.Regularity Analysis for Patterned Texture Inspection.《IEEE TRANSACTIONS ON AUTOMATION SCIENCE AND ENGINEERING》.2009,第6卷(第1期),131-144.
Regularity Analysis for Patterned Texture Inspection;Henry Y. T. Ngan ET AL;《IEEE TRANSACTIONS ON AUTOMATION SCIENCE AND ENGINEERING》;20090131;第6卷(第1期);131-144 *
荣锋等.强力输送带X光无损检测仪的研制.《光学精密工程》.2011,第19卷(第10期),2393-2401. *

Also Published As

Publication number Publication date
CN102519990A (en) 2012-06-27

Similar Documents

Publication Publication Date Title
CN102519990B (en) Fault on-line detection method of steel wire rope conveyer belt based on texture regularity analysis
CN115082418B (en) Precise identification method for automobile parts
CN109879005A (en) Device for detecting belt tearing and method
CN110111303A (en) A kind of large-scale carrier strip tearing intelligent fault detection method based on dynamic image
CN101995412B (en) Robust glass scratch defect detection method and device thereof
CN104535356A (en) Method and system for monitoring rope arrangement faults of drum steel wire rope on basis of machine vision
CN113657221B (en) Power plant equipment state monitoring method based on intelligent sensing technology
CN116839682A (en) Cable processing and manufacturing real-time monitoring system based on Internet of things
CN102590247A (en) Steel wire rope defect detection method based on X-ray image processing technology
CN108844612B (en) Transformer internal fault identification method based on mathematical statistics probability model
CN103593668B (en) Automatic crack identification method in metal plate stamping connector mechanical property test
CN106815831B (en) Method and device for detecting abrasion of steel wire rope
CN115876258B (en) Livestock and poultry breeding environment abnormity monitoring and alarming system based on multi-source data
CN112173636B (en) Method for detecting faults of belt conveyor carrier roller by inspection robot
CN101710081A (en) Method for detecting surface defects of mechanical parts based on image texture and fractal dimension
CN106442122B (en) Steel material drop hammer test fracture toughness section method for detecting percentage based on image segmentation and identification
CN117350967B (en) Optical detection-based wafer appearance defect detection system and method
CN115861294A (en) Computer vision-based concrete production abnormity detection method and device
CN110935646A (en) Full-automatic crab grading system based on image recognition
CN113720914B (en) Ultrasonic flaw detection system and ultrasonic flaw detection method
CN101509870B (en) Carbon black content grade analysis method in rubber goods
CN105083916A (en) Method for online detection of longitudinal tear fault of conveyor belt
CN109459489B (en) Elevator crack detection method based on magnetic memory detection principle
CN105021630A (en) Automatic detection method for surface damages of conveyor belt
CN110017777A (en) A kind of wire rope lay pitch Real-time and Dynamic Detection appraisal procedure

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
CB03 Change of inventor or designer information

Inventor after: Miao Changyun

Inventor after: Li Xianguo

Inventor after: Qiao Tiezhu

Inventor after: Shen Qinglei

Inventor after: Zhong Weiliang

Inventor after: Yu Jingyi

Inventor before: Miao Changyun

Inventor before: Li Xianguo

Inventor before: Shen Qinglei

Inventor before: Zhong Weiliang

Inventor before: Yu Jingyi

COR Change of bibliographic data
TR01 Transfer of patent right

Effective date of registration: 20160622

Address after: 300160 Tianjin City Hedong District Forest Road No. 63

Patentee after: Tianjin Polytechnic University

Patentee after: Tianjin Heng Yi Mechanical and Electronic Technology Co., Ltd.

Patentee after: Taiyuan University of Technology

Address before: 300160 Tianjin City Hedong District Forest Road No. 63

Patentee before: Tianjin Polytechnic University

Patentee before: Tianjin Heng Yi Mechanical and Electronic Technology Co., Ltd.