CN102201118A  Method for positioning boundaries of belt layers at tyre crown part of Xray image of tyre  Google Patents
Method for positioning boundaries of belt layers at tyre crown part of Xray image of tyre Download PDFInfo
 Publication number
 CN102201118A CN102201118A CN2010106160468A CN201010616046A CN102201118A CN 102201118 A CN102201118 A CN 102201118A CN 2010106160468 A CN2010106160468 A CN 2010106160468A CN 201010616046 A CN201010616046 A CN 201010616046A CN 102201118 A CN102201118 A CN 102201118A
 Authority
 CN
 China
 Prior art keywords
 belt
 matrix
 image
 information
 gray
 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.)
 Granted
Links
 238000001228 spectrum Methods 0.000 claims abstract description 50
 238000010183 spectrum analysis Methods 0.000 claims abstract description 3
 239000011159 matrix material Substances 0.000 claims description 67
 230000000875 corresponding Effects 0.000 claims description 22
 238000001914 filtration Methods 0.000 claims description 22
 239000000284 extract Substances 0.000 claims description 9
 101700045921 fft2 Proteins 0.000 claims description 9
 229910000831 Steel Inorganic materials 0.000 claims description 7
 239000010959 steel Substances 0.000 claims description 7
 235000020127 ayran Nutrition 0.000 claims description 6
 238000004364 calculation method Methods 0.000 claims description 3
 238000006243 chemical reaction Methods 0.000 claims description 3
 238000006073 displacement reaction Methods 0.000 claims 1
 230000002950 deficient Effects 0.000 description 6
 238000010586 diagram Methods 0.000 description 4
 238000000034 method Methods 0.000 description 4
 239000000463 material Substances 0.000 description 3
 230000000694 effects Effects 0.000 description 2
 238000004458 analytical method Methods 0.000 description 1
 230000019771 cognition Effects 0.000 description 1
 239000000470 constituent Substances 0.000 description 1
 238000005516 engineering process Methods 0.000 description 1
 238000003384 imaging method Methods 0.000 description 1
 238000000465 moulding Methods 0.000 description 1
 238000000926 separation method Methods 0.000 description 1
 230000003595 spectral Effects 0.000 description 1
 238000005987 sulfurization reaction Methods 0.000 description 1
Abstract
The invention relates to a method for positioning the boundaries of belt layers at a tyre crown part of an Xray image of a tyre. The method comprises the steps of: positioning belt layers of all the belt layers at the tyre crown part of the Xray image of a tyre, carrying out frequencyspectrum analysis on the tyre crown part of the Xray image of the tyre with overlapped multiple belt layer textures, respectively processing frequencyspectrum information by using a Gaussian wedgeshaped filter with the angle between 0 degree and 179 degrees, restoring the extracted frequencyspectrum information into images, carrying out grey statistics to obtain gray statistical information vectors, combining all the gray statistical information vectors and converting the combined gray statistical information vectors into a gray information image, observing and selecting the gray information vectors which can be used for measuring and calculating the boundary information of all the belt layers, and calculating the boundary position coordinates of all the belt layers according to the selected gray information vectors.
Description
[technical field]
The invention belongs to the digital image processing techniques field, be specifically related to be applied to many texture images that crown position in the tire xray image forms because of many belts texture mutual superposition are carried out each belt boundary alignment.
[background technology]
Tire is that member is carried out in the main action of motor vehicles, and stablizing and whether meeting the Safety Design standard of tyre performance will directly determine to use the personal security of motor vehicles.Tire is the rotatable body of a kind of circular ring type of tubular section, and the rubber preformed material, the compounded rubber preformed material that are had allsteel cord by multilayer process through applying, moulding, sulfuration typing.The tire of described structure, the constituent material at its crown position inevitably exist problems such as allsteel cord density unevenness, wire of steel wire tire cord fracture, allsteel cord disappearance.Perspective imaging principle according to Xray can obtain inside tires allsteel cord state and distribution situation, by to the state of allsteel cord and the analysis of distribution situation, obtains the information of defective, thereby realizes tire is declared level.
Existing tire Xray machine test unit, in process of the test, tire is mounted, is clamped on four rotatable reference columns.In the time of the reference column rotation, drive the tire rotation, at this moment, the Xray generator tube emission Xray line by being positioned at the tire inboard penetrates tire, is received by the U type receiver that is positioned at the tire outside.
Existing tire Xray machine test unit after obtaining the xray image of tire, is a foundation with the xray image of tire, declares level.Declare the observation of the main people of dependence of level at present, the tire xray image defective that is presented on the computer screen is carried out artificial cognition,, confirm the grade of tire by artificial intuitive judgment by naked eyes.At the crown position, because a plurality of belt mutual superposition are arranged, belt offset or scarce layer, the defective that can have a strong impact on tire xray image crown position is declared level.The belt that overlaps can be covered the distribution situation of each belt each other mutually, therefore at the crown position of tire xray image, with the naked eye observe the boundary position that is difficult to tell belt intuitively, the unascertainable words of the boundary position of belt, automatically declaring level at the defective of belt just is difficult to carry out, therefore need each belt be positioned before the computing machine automatic defect declares level carrying out, followup like this declare level work and could shoot the arrow at the target.
[summary of the invention]
The present invention seeks to overcome the prior art above shortcomings, a kind of tire xray image crown part position belt boundary alignment method that is applied to is provided.
The method that the present invention adopts Digital Image Processing is carried out boundary alignment to the belt at crown position, mark off the effective coverage of each belt, at followup declaring in the level work, help in concrete zone, to declare targetedly level operation and defect recognition, defective is divided by diverse location, different ownership, simultaneously help also to judge whether belt is offset, whether have the layer of lacking phenomenon to take place.
Spectrum analysis is carried out at the tire xray image crown position that the present invention is directed to many belts texture mutual superposition, Gauss's wedge shape wave filter is accepted or rejected spectrum information when adopting angle to be 0 ° to 179 ° respectively, the spectrum information that extracts is reduced to image and carries out grayscale statistical, obtain the grayscale statistical information vector, all grayscale statistical information vectors are made up and are converted into the halftone information image, by observing the halftone information vector of selecting to can be used for calculating each belt boundary information, according to the halftone information vector that is selected, calculate the boundary position coordinate of each belt, thereby each belt Boundary Recognition of mutual superposition is come out, realize the boundary alignment of each belt, implementation step is:
The first step is obtained a complete tire xray image I
_{X}, at this xray image I
_{X}Inside also is simultaneously that its crown position intercepting size is 2
^{N}* 2
^{N}Squareshaped image I
_{O}The value of N requires: N is the integer greater than 1, has 2 simultaneously
^{N}Less than tire xray image I
_{X}Width and 2
^{N}Less than xray image I
_{X}Height, and make I
_{O}Width can cover whole crown position just, image I simultaneously
_{O}Four limits should be parallel to tire xray image I respectively
_{X}Four limits;
Second step is to first step truncated picture I
_{O}Carry out 2 dimension Fast Fourier Transform (FFT)s, obtain spectrum information matrix I
_{F}, with spectrum information matrix I
_{F}Carry out shifting function, and then to obtain with frequency (0,0) be 2 dimension spectrum information matrix I of central point
_{FS}, I
_{F}, I
_{FS}All and image I
_{O}With size, I
_{F}, I
_{FS}With image I
_{O}Meet following relation:
I
_{F}＝fft2(I
_{O})
I
_{FS}＝fftshift(I
_{F})
Wherein fft2 represents to carry out quick 2 dimension Fourier transforms, and in the result images after the fft2 conversion, the pairing value of information of frequency (0,0) is corresponding to spectrum information matrix I
_{F}Header element position, the upper left corner;
Fftshift represents spectrum information matrix I
_{F}Carrying out shifting function, is initial point with the matrix center promptly, does level and vertical coordinate axis, with the whole exchange of 1,3 quadrants, with the whole matrix that obtains that exchanges of 2,4 quadrants; After the exchange, the pairing value of information of frequency (0,0) is corresponding to the center of matrix;
In the 3rd step, utilize following Gauss's wedge shape wave filter formula, with spectrum information matrix I
_{FS}Carry out the filtering of Gauss's wedge shape,
G
_{2}(x，y)＝G
_{1}(x，y)
G(x，y)＝G
_{1}(x，y)+G
_{2}(x，y)
Wherein, θ
_{C}Be the angle center, its span is 0 °≤θ
_{c}＜180 °
σ is a standard deviation, has determined the width of Gaussian curve on the circumferential direction of certain radius
θ is spectrum information matrix I
_{FS}In certain first vegetarian refreshments be the summit with its frequency central point (0,0), be 0 ° angle value with horizontal dextrad direction, as 0 °≤θ
_{c}In the time of＜180 °, the θ span is90 °≤θ≤270 °, and the computing method of θ are:
Work as I
_{FS}In the coordinate of certain element (x when y) being positioned at 1,4 quadrants, has
Work as I
_{FS}In the coordinate of certain pixel (x when y) being positioned at 2,3 quadrants, has
G is 2 dimension Gauss wedge shape wave filters, and this is one is the filtering template of former point symmetry, i.e. G about frequency central point (0,0)
_{1}And G
_{2}Be about frequency central point (0,0) symmetry;
At first make θ
_{C}=0 °, and according to the design parameter of each belt, the σ parameter of given Gauss's wedge shape wave filter is utilized Gauss's wedge shape filtering formula, generates the Gauss's wedge shape wave filter G about (0,0) former point symmetry;
The 4th step is with spectrum information matrix I
_{FS}G carries out dot product with Gauss's wedge shape wave filter, obtains spectrum information matrix I
_{FL}, this matrix has at utmost kept direction, the spectrum information of expectation texture, has abandoned direction, the spectrum information of other textures, and the dot product formula is as follows;
I
_{FL}＝I
_{FS}·G
The 5th step is with spectrum information matrix I
_{FL}Carry out shifting function, and carry out inverse fast Fourier transform, thereby with I
_{FL}Be reduced to image I
_{R}, I
_{FL}With I
_{R}Satisfy following relation:
I
_{R}＝ifft2(fftshift(I
_{FL}))
Wherein, fftshift represents spectrum information matrix I
_{FL}Carrying out shifting function, is initial point with the matrix center promptly, does level and vertical coordinate axis, with the whole exchange of 1,3 quadrants, the whole exchange of 2,4 quadrants is obtained new matrix; Before the exchange, the pairing value of information of frequency (0,0) is corresponding to spectrum information matrix I
_{FL}The center; After the exchange, the pairing value of information of frequency (0,0) is corresponding to the header element position, the upper left corner of matrix of consequence,
Ifft2 represents to carry out quick 2 dimension inverse Fourier transforms,
Through the image I after the filtering reduction
_{R}, farthest kept the texture of certain direction, and the texture information of other directions is weakened;
In the 6th step, calculate I
_{R}Pixel grey scale average on the vertical direction of each horizontal ordinate position obtains one 1 vector of tieing up
Computing formula is as follows:
Wherein, I
_{R}(i, m) presentation video I
_{R}In certain pixel (i represents abscissa value for i, grey scale pixel value m), and m represents ordinate value;
In the 7th step, make θ again
_{C}Be respectively 1 °, 2 °, 3 ° ... 179 °, repeat the 3rd and went on foot for the 6th step, obtain
V
_{θ}(i) θ＝0，1，…，179；i＝1，2，…2
^{N}
And with following formula vector set V
_{θ}(i) θ=0,1 ..., 179; I=1,2 ... 2
^{N}Be rewritten as 2 dimension matrixes
V(θ，i) θ＝0，1，…，179；i＝1，2，…2
^{N}
Being that matrix V is one 2 dimension matrix, highly is 180, and width is 2
^{N}, this 2 dimension matrix is shown with image mode, can observe each belt boundary position distribution situation;
In the 8th step, the observation by to 2 dimension image V can pick out and the horizontal strip region of the some one to one high grayscale values of each belt, and hypothesis has M belt here, and the horizontal strip region of then high grayscale value is made as Zone
_{i}I=1,2, M, the number of the horizontal strip region of high grayscale value is consistent with the belt number, simultaneously the horizontal strip region of each high grayscale value is corresponding with some belts, and the width of the horizontal strip region of high grayscale value equals the width of the belt corresponding with it, and the left and right sides boundary coordinate of the horizontal strip region of high grayscale value equals the belt left and right sides border corresponding with it at I
_{O}In coordinate, will with the horizontal strip region Zone of high grayscale value
_{i}I=1,2 ..., the abscissa value β of M
_{i}I=1,2 ..., M extracts, can extract with each belt one to one
I=1,2 ..., the M vector is at M
I=1,2 ..., the M vector can calculate M belt boundary information;
In the 9th step, calculate
I=1,2 ..., the average of M
I=1,2 ..., M is with average
I=1,2 ..., M multiplication by constants coefficient lambda
_{i}I=1,2 ..., M obtains threshold value T
_{i}I=1,2 ..., M, promptly satisfy following relational expression:
At each vector
I=1,2 ..., M, (being x=0) begins search from the left side, when the first fit condition
The time, then Ci Shi x is designated as the left margin F on the border of asking
_{Li}I=1,2 ..., M;
In like manner, from vector
I=1,2 ..., the right side of M (is x=2
^{N}) beginning search left, when the first fit condition
The time, then Ci Shi x is designated as the right margin E on the border of asking
_{Ri}I=1,2 ..., M;
At gray level image I
_{O}In, respectively with I
_{O}(E
_{Li}, y), I
_{O}(E
_{Ri}, y) i=1,2 ..., M; Y=1,2 ..., 2
^{N}The pixel value grayscale value is set to (255,255,255), i.e. boundary line of drawing each belt with white lines;
So far, finish the belt boundary alignment.
Advantage of the present invention and good effect
The present invention can pick out the border of each belt accurately, mark off the effective coverage of each belt, at followup declaring in the level work, help in concrete zone, to declare targetedly level operation and defect recognition, defective is divided by diverse location, different ownership, simultaneously help also to judge whether belt is offset, whether have the layer of lacking phenomenon to take place; What the method helped the tire xray image accurately declares level, improves accuracy rate and efficient that tire is declared level.
[description of drawings]
Fig. 1 is a tire xray image crown position intercepting synoptic diagram;
Fig. 2 is Gauss's wedge shape filtering template synoptic diagram;
Fig. 3 is the process flow diagram of described belt texture boundary alignment method;
Fig. 4 is the average statistical graph that obtains after filtering under certain angle direction;
Fig. 5 is after adopting 0 °～179 ° Gauss's wedge shape wave filters that the spectrum information of tire xray image crown part bit image is carried out filtering respectively, the grayscale statistical information matrix that shows with image format that obtains;
Fig. 6 is boundary alignment final effect figure.
[embodiment]
As shown in Figure 1, tyre crown position many belts of the xray image texture separation method that is applied to of the present invention, its effective object is the subimage of the suitable size in tire xray image crown position.The size of subimage is 2
^{9}* 2
^{9}, subimage has covered the left hand edge and the right hand edge at crown position; Four limits of subimage are parallel to four limits of tire xray image simultaneously;
As shown in Figure 2, this is to adopt Gauss's wedge shape filtering formula, Gauss's wedge shape filtering template synoptic diagram of generation.This Gauss's wedge shape filtering template is about central point as can be known.
As shown in Figure 3, the concrete steps of tire xray image crown part position belt boundary alignment method are as follows:
The first step is obtained a complete tire xray image I
_{X}, at this xray image I
_{X}Inside also is simultaneously that its crown position intercepting size is 2
^{N}* 2
^{N}Squareshaped image I
_{O}The value of N requires: N is the integer greater than 1, has 2 simultaneously
^{N}Less than tire xray image I
_{X}Width and 2
^{N}Less than xray image I
_{X}Height, and make I
_{O}Width can cover whole crown position just, image I simultaneously
_{O}Four limits should be parallel to tire xray image I respectively
_{X}Four limits;
Second step is to first step truncated picture I
_{O}Carry out 2 dimension Fast Fourier Transform (FFT)s, obtain spectrum information matrix I
_{F}, with spectrum information matrix I
_{F}Carry out shifting function, and then to obtain with frequency (0,0) be 2 dimension spectrum information matrix I of central point
_{FS}, I
_{F}, I
_{FS}All and image I
_{O}With size, I
_{F}, I
_{FS}With image I
_{O}Meet following relation:
I
_{F}＝fft2(I
_{O})
I
_{FS}＝fftshift(I
_{F})
Wherein fft2 represents to carry out quick 2 dimension Fourier transforms, and in the result images after the fft2 conversion, the pairing value of information of frequency (0,0) is corresponding to spectrum information matrix I
_{F}Header element position, the upper left corner.
Fftshift represents I
_{F}Carrying out shifting function, is initial point with the matrix center promptly, does level and vertical coordinate axis, with the whole exchange of 1,3 quadrants, with the whole matrix that obtains that exchanges of 2,4 quadrants.After the exchange, the pairing value of information of frequency (0,0) is corresponding to the center of matrix.
In the 3rd step, utilize following Gauss's wedge shape wave filter formula, with spectrum information matrix I
_{FS}Carry out the filtering of Gauss's wedge shape,
G
_{2}(x，y)＝G
_{1}(x，y)
G(x，y)＝G
_{1}(x，y)+G
_{2}(x，y)
Wherein, θ
_{C}Be the angle center, its span is 0 °≤θ
_{c}＜180 °
σ is a standard deviation, has determined the width of Gaussian curve on the circumferential direction of certain radius
θ is I
_{FS}In certain first vegetarian refreshments be the summit with its frequency central point (0,0), be 0 ° angle value with horizontal dextrad direction, as 0 °≤θ
_{c}In the time of＜180 °, the θ span is90 °≤θ≤270 °, and the computing method of θ are:
Work as I
_{FS}In the coordinate of certain element (x when y) being positioned at 1,4 quadrants, has
Work as I
_{FS}In the coordinate of certain pixel (x when y) being positioned at 2,3 quadrants, has
G is 2 dimension Gauss wedge shape wave filters, and this is one is the filtering template of former point symmetry, i.e. G about frequency central point (0,0)
_{1}And G
_{2}Be about frequency central point (0,0) symmetry;
At first make θ
_{C}=0 °, and according to the design parameter of each belt, the σ parameter of given Gauss's wedge shape wave filter is utilized Gauss's wedge shape filtering formula, generates the Gauss's wedge shape wave filter G about (0,0) former point symmetry;
The 4th step is with spectral matrix I
_{FS}G carries out dot product with Gauss's wedge shape wave filter, obtains spectrum information matrix I after the filtering
_{FL}, this matrix has at utmost kept direction, the spectrum information of expectation texture, has abandoned direction, the spectrum information of other textures, and the dot product formula is as follows;
I
_{FL}＝I
_{FS}·G
The 5th step is with filtered spectrum information matrix I
_{FL}Carry out shifting function, and carry out inverse fast Fourier transform, thereby with I
_{FL}Be reduced to image I
_{R}, I
_{FL}With I
_{R}Satisfy following relation:
I
_{R}＝ifft2(fftshift(I
_{FL}))
Wherein, fftshift represents I
_{FL}Carrying out shifting function, is initial point with the matrix center promptly, does level and vertical coordinate axis, with the whole exchange of 1,3 quadrants, the whole exchange of 2,4 quadrants is obtained new matrix.Before the exchange, the pairing value of information of frequency (0,0) is corresponding to I
_{FL}The center; After the exchange, the pairing value of information of frequency (0,0) is corresponding to the header element position, the upper left corner of matrix of consequence,
Ifft2 represents to carry out quick 2 dimension inverse Fourier transforms,
Through the image I after the filtering reduction
_{R}, farthest kept the texture of certain direction, and the texture information of other directions is weakened;
In the 6th step, calculate I
_{R}Pixel grey scale average on the vertical direction of each horizontal ordinate position obtains one 1 vector of tieing up
Computing formula is as follows:
Wherein, I
_{R}(i, m) presentation video I
_{R}In certain pixel (i represents abscissa value for i, grey scale pixel value m), and m represents ordinate value;
In the 7th step, make θ again
_{C}Be respectively 1 °, 2 °, 3 ° ... 179 °, repeat the 3rd and went on foot for the 6th step, obtain
V
_{θ}(i) θ＝0，1，…，179；i＝1，2，…2
^{N}
And with following formula vector set V
_{θ}(i) θ=0,1 ..., 179; I=1,2 ... 2
^{N}Be rewritten as 2 dimension matrixes
V(θ，i) θ＝0，1，…，179；i＝1，2，…2
^{N}
Being that matrix V is one 2 dimension matrix, highly is 180, and width is 2
^{N}, this 2 dimension matrix is shown with image mode, can observe each belt boundary position distribution situation;
In the 8th step, the observation by to 2 dimension image V can pick out and the horizontal strip region of the some one to one high grayscale values of each belt, and hypothesis has M belt here, and the horizontal strip region of then high grayscale value is made as Zone
_{i}
I=1,2, M, the number of the horizontal strip region of high grayscale value is consistent with the belt number, simultaneously the horizontal strip region of each high grayscale value is corresponding with some belts, and the width of the horizontal strip region of high grayscale value equals the width of the belt corresponding with it, and the left and right sides boundary coordinate of the horizontal strip region of high grayscale value equals the belt left and right sides border corresponding with it at I
_{O}In coordinate, will with the horizontal strip region Zone of high grayscale value
_{i}I=1,2 ..., the abscissa value β of M
_{i}
I=1,2 ..., M extracts, can extract with each belt one to one
I=1,2 ..., the M vector is at M
I=1,2 ..., the M vector can calculate M belt boundary information;
In the 9th step, calculate
I=1,2 ..., the average of M
I=1,2 ..., M is with average
I=1,2 ..., M multiplication by constants coefficient lambda
_{i}I=1,2 ..., M obtains threshold value T
_{i}I=1,2 ..., M, promptly satisfy following relational expression:
At each vector
I=1,2 ..., M, (being x=0) begins search from the left side, when the first fit condition
The time, then Ci Shi x is designated as the left margin E on the border of asking
_{Li}I=1,2 ..., M;
In like manner, from vector
I=1,2 ..., the right side of M (is x=2
^{N}) beginning search left, when the first fit condition
The time, then Ci Shi x is designated as the right margin E on the border of asking
_{Ri}I=1,2 ..., M;
At gray level image I
_{O}In, respectively with I
_{O}(E
_{Li}, y), I
_{O}(E
_{Ri}, y) i=1,2 ..., M; Y=1,2 ..., 2
^{N}The pixel value grayscale value is set to (255,255,255), i.e. boundary line of drawing each belt with white lines;
So far, finish the belt boundary alignment.
As shown in Figure 4, this is after adopting certain Gauss's wedge shape wave filter that filtering is carried out at tire xray image crown position, the grayscale statistical curve of the filtering result images that obtains; Horizontal ordinate has been represented the width of image, and ordinate is represented the pixel grey scale average at certain horizontal ordinate place of filtering result images.
As shown in Figure 5, this is the grayscale statistical matrix of being made up of 180 grayscale statisticals vector that shows with image format, 1,2, No. 3 high grayscale value zone correspondence wherein 3 different belts.As can be seen, the boundary position in 3 high grey value profile zones is that significantly the border, the left and right sides of extracting these 3 high grayscale value zones can obtain the border, the left and right sides of 3 belts.
As shown in Figure 6, the 3 pairs of vertical straight lines of white that on the crown part bit image of tire xray image, draw, the corresponding border, the left and right sides of 3 belts.
Claims (2)
1. one kind is applied to tire xray image crown part position belt boundary alignment method, it is characterized in that: spectrum analysis is carried out at the tire xray image crown position at many belts texture mutual superposition, Gauss's wedge shape wave filter is accepted or rejected spectrum information when adopting angle to be 0 ° to 179 ° respectively, the spectrum information that extracts is reduced to image and carries out grayscale statistical, obtain the grayscale statistical information vector, all grayscale statistical information vectors are made up and are converted into the halftone information image, by observing the halftone information vector of selecting to can be used for calculating each belt boundary information, according to the halftone information vector that is selected, calculate the boundary position coordinate of each belt, thereby each belt Boundary Recognition of mutual superposition is come out, realize the boundary alignment of each belt, implementation step is:
The first step is obtained a complete tire xray image I
_{X}, at this xray image I
_{X}Inside also is simultaneously that its crown position intercepting size is 2
^{N}* 2
^{N}Squareshaped image I
_{O}The value of N requires: N is the integer greater than 1, has 2 simultaneously
^{N}Less than tire xray image I
_{X}Width and 2
^{N}Less than xray image I
_{X}Height, and make I
_{O}Width can cover whole crown position just, image I simultaneously
_{O}Four limits should be parallel to tire xray image I respectively
_{X}Four limits;
Second step is to first step truncated picture I
_{O}Carry out 2 dimension Fast Fourier Transform (FFT)s, obtain spectrum information matrix I
_{F}, with spectrum information matrix I
_{F}Carry out shifting function, and then to obtain with frequency (0,0) be 2 dimension spectrum information matrix I of central point
_{FS}, I
_{F}, I
_{FS}All and image I
_{O}With size, I
_{F}, I
_{FS}With image I
_{O}Meet following relation:
I
_{F}＝fft2(I
_{O})
I
_{FS}＝fftshift(I
_{F})
Wherein fft2 represents to carry out quick 2 dimension Fourier transforms, and in the result images after the fft2 conversion, the pairing value of information of frequency (0,0) is corresponding to spectrum information matrix I
_{F}Header element position, the upper left corner;
Fftshift represents spectrum information matrix I
_{F}Carrying out shifting function, is initial point with the matrix center promptly, does level and vertical coordinate axis, with the whole exchange of 1,3 quadrants, with the whole matrix that obtains that exchanges of 2,4 quadrants; After the exchange, the pairing value of information of frequency (0,0) is corresponding to the center of matrix;
In the 3rd step, utilize following Gauss's wedge shape wave filter formula, with spectrum information matrix I
_{FS}Carry out the filtering of Gauss's wedge shape,
G
_{2}(x，y)＝G
_{1}(x，y)
G(x，y)＝G
_{1}(x，y)+G
_{2}(x，y)
Wherein, θ
_{C}Be the angle center, its span is 0 °≤θ
_{c}＜180 °
σ is a standard deviation, has determined the width of Gaussian curve on the circumferential direction of certain radius
θ is spectrum information matrix I
_{FS}In certain first vegetarian refreshments be the summit with its frequency central point (0,0), be 0 ° angle value with horizontal dextrad direction, as 0 °≤θ
_{c}In the time of＜180 °, the θ span is90 °≤θ≤270 °, and the computing method of θ are:
Work as I
_{FS}In the coordinate of certain element (x when y) being positioned at 1,4 quadrants, has
Work as I
_{FS}In the coordinate of certain pixel (x when y) being positioned at 2,3 quadrants, has
G is 2 dimension Gauss wedge shape wave filters, and this is one is the filtering template of former point symmetry, i.e. G about frequency central point (0,0)
_{1}And G
_{2}Be about frequency central point (0,0) symmetry;
At first make θ
_{C}=0 °, and according to the design parameter of each belt, the σ parameter of given Gauss's wedge shape wave filter is utilized Gauss's wedge shape filtering formula, generates the Gauss's wedge shape wave filter G about (0,0) former point symmetry;
The 4th step is with spectrum information matrix I
_{FS}G carries out dot product with Gauss's wedge shape wave filter, obtains spectrum information matrix I
_{FL}, this matrix has at utmost kept direction, the spectrum information of expectation texture, has abandoned direction, the spectrum information of other textures, and the dot product formula is as follows;
I
_{FL}＝I
_{FS}·G
The 5th step is with spectrum information matrix I
_{FL}Carry out shifting function, and carry out inverse fast Fourier transform, thereby with I
_{FL}Be reduced to image I
_{R}, I
_{FL}With I
_{R}Satisfy following relation:
I
_{R}＝ifft2(fftshift(I
_{FL}))
Wherein, fftshift represents spectrum information matrix I
_{FL}Carrying out shifting function, is initial point with the matrix center promptly, does level and vertical coordinate axis, with the whole exchange of 1,3 quadrants, the whole exchange of 2,4 quadrants is obtained new matrix; Before the exchange, the pairing value of information of frequency (0,0) is corresponding to spectrum information matrix I
_{FL}The center; After the exchange, the pairing value of information of frequency (0,0) is corresponding to the header element position, the upper left corner of matrix of consequence,
Ifft2 represents to carry out quick 2 dimension inverse Fourier transforms,
Through the image I after the filtering reduction
_{R}, farthest kept the texture of certain direction, and the texture information of other directions is weakened;
In the 6th step, calculate I
_{R}Pixel grey scale average on the vertical direction of each horizontal ordinate position obtains one 1 vector of tieing up
Computing formula is as follows:
Wherein, I
_{R}(i, m) presentation video I
_{R}In certain pixel (i represents abscissa value for i, grey scale pixel value m), and m represents ordinate value;
In the 7th step, make θ again
_{C}Be respectively 1 °, 2 °, 3 ° ... 179 °, repeat the 3rd and went on foot for the 6th step, obtain
V
_{θ}(i) θ＝0，1，…，179；i＝1，2，…2
^{N}
And with following formula vector set V
_{θ}(i) θ=0,1 ..., 179; I=1,2 ... 2
^{N}Be rewritten as 2 dimension matrixes
V(θ，i) θ＝0，1，…，179；i＝1，2，…2
^{N}
Being that matrix V is one 2 dimension matrix, highly is 180, and width is 2
^{N}, this 2 dimension matrix is shown with image mode, can observe each belt boundary position distribution situation;
In the 8th step, the observation by to 2 dimension image V can pick out and the horizontal strip region of the some one to one high grayscale values of each belt, and hypothesis has M belt here, and the horizontal strip region of then high grayscale value is made as Zone
_{i}I=1,2, M, the number of the horizontal strip region of high grayscale value is consistent with the belt number, simultaneously the horizontal strip region of each high grayscale value is corresponding with some belts, and the width of the horizontal strip region of high grayscale value equals the width of the belt corresponding with it, and the left and right sides boundary coordinate of the horizontal strip region of high grayscale value equals the belt left and right sides border corresponding with it at I
_{O}In coordinate, will with the horizontal strip region Zone of high grayscale value
_{i}I=1,2 ..., the abscissa value β of M
_{i}I=1,2 ..., M extracts, can extract with each belt one to one
I=1,2 ..., the M vector is at M
I=1,2 ..., the M vector can calculate M belt boundary information;
In the 9th step, calculate
I=1,2 ..., the average of M
I=1,2 ..., M is with average
I=1,2 ..., M multiplication by constants coefficient lambda
_{i}I=1,2 ..., M obtains threshold value T
_{i}I=1,2 ..., M, promptly satisfy following relational expression:
At each vector
I=1,2 ..., M, (being x=0) begins search from the left side, when the first fit condition
The time, then Ci Shi x is designated as the left margin E on the border of asking
_{Li}I=1,2 ..., M;
In like manner, from vector
I=1,2 ..., the right side of M (is x=2
^{N}) beginning search left, when the first fit condition
The time, then Ci Shi x is designated as the right margin E on the border of asking
_{Ri}I=1,2 ..., M;
At gray level image I
_{O}In, respectively with I
_{O}(E
_{Li}, y), I
_{O}(E
_{Ri}, y) i=1,2 ..., M; Y=1,2 ..., 2
^{N}The pixel value grayscale value is set to (255,255,255), i.e. boundary line of drawing each belt with white lines;
So far, finish the belt boundary alignment.
2. method according to claim 1 is characterized in that: be applied to separability that tire xray image crown part position belt boundary alignment method utilized spectrum information to realize the belt boundary alignment, separability is meant:
Crown position at the tire xray image, a plurality of belt texture mutual superposition, what therefore the crown position of the xray image of tire demonstrated is the texture of mutual each the belt allsteel cord that overlaps, the border of each belt is covered by each belt texture of mutual superposition, tire xray image crown part bit image I
_{O}In each belt include different texture informations, contain frequency spectrum, direction etc., to tire xray image crown part bit image I
_{O}After carrying out Fast Fourier Transform (FFT) and displacement, the spectrum information matrix I that obtains
_{FS}In, the texture information correspondence that different belt comprised spectrum information matrix I
_{FS}In different local locations, the spectrum information that is in different local locations is separable, extracts spectrum information matrix I by Gauss's wedge shape wave filter
_{FS}In the information of certain regional area, just can obtain the texture information of expectation.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN 201010616046 CN102201118B (en)  20101230  20101230  Method for positioning boundaries of belt layers at tyre crown part of Xray image of tyre 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN 201010616046 CN102201118B (en)  20101230  20101230  Method for positioning boundaries of belt layers at tyre crown part of Xray image of tyre 
Publications (2)
Publication Number  Publication Date 

CN102201118A true CN102201118A (en)  20110928 
CN102201118B CN102201118B (en)  20130109 
Family
ID=44661768
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN 201010616046 Active CN102201118B (en)  20101230  20101230  Method for positioning boundaries of belt layers at tyre crown part of Xray image of tyre 
Country Status (1)
Country  Link 

CN (1)  CN102201118B (en) 
Cited By (4)
Publication number  Priority date  Publication date  Assignee  Title 

CN102592292A (en) *  20111212  20120718  河南理工大学  Symmetric figure center positioning method based on inner integral operation 
CN105793697A (en) *  20131203  20160720  横滨橡胶株式会社  Tire inspection method and device therefor 
CN107833209A (en) *  20171027  20180323  浙江大华技术股份有限公司  A kind of xray image detection method, device, electronic equipment and storage medium 
CN109829899A (en) *  20190118  20190531  创新奇智(广州)科技有限公司  A kind of background restrainable algorithms for steel coil endface defects detection 
Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN101173907A (en) *  20061102  20080507  张芝泉  X optical detection machine for project tire 
JP2008309646A (en) *  20070614  20081225  Yokohama Rubber Co Ltd:The  Method and device for inspecting tire 
KR20090023974A (en) *  20070903  20090306  금호타이어 주식회사  Xray inspecting divice for tire 
CN101603931A (en) *  20080612  20091216  青岛高校软控股份有限公司  Tyre Xray detection device and method thereof 
CN101685073A (en) *  20080926  20100331  青岛高校软控股份有限公司  Xray apparatus test device and method for truck tire 

2010
 20101230 CN CN 201010616046 patent/CN102201118B/en active Active
Patent Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN101173907A (en) *  20061102  20080507  张芝泉  X optical detection machine for project tire 
JP2008309646A (en) *  20070614  20081225  Yokohama Rubber Co Ltd:The  Method and device for inspecting tire 
KR20090023974A (en) *  20070903  20090306  금호타이어 주식회사  Xray inspecting divice for tire 
CN101603931A (en) *  20080612  20091216  青岛高校软控股份有限公司  Tyre Xray detection device and method thereof 
CN101685073A (en) *  20080926  20100331  青岛高校软控股份有限公司  Xray apparatus test device and method for truck tire 
NonPatent Citations (1)
Title 

黄战华等: "基于统计特征的轮胎纹理缺陷在线检测", 《光学技术》 * 
Cited By (8)
Publication number  Priority date  Publication date  Assignee  Title 

CN102592292A (en) *  20111212  20120718  河南理工大学  Symmetric figure center positioning method based on inner integral operation 
CN102592292B (en) *  20111212  20140716  河南理工大学  Symmetric figure center positioning method based on inner integral operation 
CN105793697A (en) *  20131203  20160720  横滨橡胶株式会社  Tire inspection method and device therefor 
CN105793697B (en) *  20131203  20181019  横滨橡胶株式会社  Tire checkup method and its device 
CN107833209A (en) *  20171027  20180323  浙江大华技术股份有限公司  A kind of xray image detection method, device, electronic equipment and storage medium 
CN107833209B (en) *  20171027  20200526  浙江大华技术股份有限公司  Xray image detection method and device, electronic equipment and storage medium 
CN109829899A (en) *  20190118  20190531  创新奇智(广州)科技有限公司  A kind of background restrainable algorithms for steel coil endface defects detection 
CN109829899B (en) *  20190118  20200807  创新奇智(广州)科技有限公司  Background suppression algorithm for steel coil end surface defect detection 
Also Published As
Publication number  Publication date 

CN102201118B (en)  20130109 
Similar Documents
Publication  Publication Date  Title 

CN102298779B (en)  Image registering method for panoramic assisted parking system  
CN103886308B (en)  A kind of pedestrian detection method of use converging channels feature and soft cascade grader  
US8774502B2 (en)  Method for image/video segmentation using texture feature  
CN102682455B (en)  Front vehicle detection method based on monocular vision  
CN106934795B (en)  A kind of automatic testing method and prediction technique of glue into concrete beam cracks  
CN104732227B (en)  A kind of Location Method of Vehicle License Plate based on definition and luminance evaluation  
CN102201118B (en)  Method for positioning boundaries of belt layers at tyre crown part of Xray image of tyre  
CN106033610B (en)  Profile testing method based on the modulation of nonclassical receptive field spatial summation  
CN103886612B (en)  Automatic water level extraction method and system based on reservoir monitoring camera  
CN104966285B (en)  A kind of detection method of salient region  
CN105744256A (en)  Threedimensional image quality objective evaluation method based on graphbased visual saliency  
KR20120019425A (en)  Image processing device, method, and program  
CN105976392B (en)  Vehicle tyre detection method and device based on maximum output probability  
CN106778668A (en)  A kind of method for detecting lane lines of the robust of joint RANSAC and CNN  
CN104036480A (en)  Surf algorithm based quick mismatching point eliminating method  
CN106485224A (en)  A kind of seatbelt wearing recognition methodss and device  
CN105930848A (en)  SARSIFT featurebased SAR image target recognition method  
CN101739549A (en)  Face detection method and system  
CN104966054A (en)  Weak and small object detection method in visible image of unmanned plane  
CN104517269A (en)  Terahertz image strip processing method and system  
CN107767387A (en)  Profile testing method based on the global modulation of changeable reception field yardstick  
CN103456006A (en)  Oil depot target positioning method of SAR image  
CN104463896B (en)  Image corner point detection method and system based on kernel similar region distribution characteristics  
CN102081797B (en)  Multiple belted layer texture separation method applied to tyre Xray image crown position  
CN102567990B (en)  Stereo image objective quality estimation method 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
SE01  Entry into force of request for substantive examination  
C10  Entry into substantive examination  
GR01  Patent grant  
C14  Grant of patent or utility model 