CN102346013A - Tunnel lining crack width measuring method and device - Google Patents

Tunnel lining crack width measuring method and device Download PDF

Info

Publication number
CN102346013A
CN102346013A CN2010102405432A CN201010240543A CN102346013A CN 102346013 A CN102346013 A CN 102346013A CN 2010102405432 A CN2010102405432 A CN 2010102405432A CN 201010240543 A CN201010240543 A CN 201010240543A CN 102346013 A CN102346013 A CN 102346013A
Authority
CN
China
Prior art keywords
crack
image
pixel
zone
tunnel lining
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.)
Pending
Application number
CN2010102405432A
Other languages
Chinese (zh)
Inventor
朱合华
刘学增
叶康
罗仁立
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji 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 Tongji University filed Critical Tongji University
Priority to CN2010102405432A priority Critical patent/CN102346013A/en
Publication of CN102346013A publication Critical patent/CN102346013A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention relates to a tunnel lining crack width measuring method, which comprises the follow steps of: (1) picking up crack images through a digital camera and measuring pickup distance through a laser distance measuring instrument; (2) selecting a region required to be measured on the picked images and converting selected region images into gray level images; (3) determining the threshold value of the gray level images selected in the second step, cutting binary image cutting, and obtaining a target region; (4) extracting a white region only comprising crack in the target region; (5) extracting the edge of the crack by a sub pixel method, and obtaining sub pixel level edge images; (6) circulating the pixel width of the crack by a minimum distance method; (7) circulating the actual width of the obtained crack according to the calibration proportion of the predetermined actual pixel size and the image pickup distance; and (8) storing the obtained actual width data into a crack data base. Compared with the prior art, the tunnel lining crack width measuring method has the advantages of low cost, high speed, convenience, high precision and the like.

Description

A kind of measuring method of Tunnel Lining Cracks width and device
Technical field
The present invention relates to a kind of Tunnel Lining Cracks measuring technique, especially relate to a kind of measuring method and device of Tunnel Lining Cracks width.
Background technology
The quantification of defectives such as crack is the target that Non-Destructive Testing is pursued always.Progress along with detection technique; Traditional human eye that passes through is visual or use the method for the mensuration fracture width that simple instrument such as reading glass estimate; Because its personnel's subjectivity is bigger, and precision and efficient are lower, with being replaced by new method gradually.The fast development of Along with computer technology and correlation theory constantly perfect, digital image processing techniques receive extensive attention and have obtained great pioneering achievement in many applications.And be deep into already in the field of civil engineering such as fracture width measurement, deformation monitoring, the identification of rock mass rubble, and bringing into play noncontact, convenient relatively, directly perceived and accurate advantage based on the photographic measurement technology of Digital Image Processing.
2004, people's independent developments such as the Huang Zhanhua of University Of Tianjin, Li Meng one cover crack identification and analysis software.In the same year, the Zou Yiqun of BJ University of Aeronautics & Astronautics, Hou Guicang, Yang Feng have proposed a kind of surface crack detection method based on Digital Image Processing.The same year, Zhang Juan, Sha Aimin, Gao Huaigang, grandson towards cloud analysis based on the principle of work of the pavement crack of Digital Image Processing identification with evaluation system.2005, the integrated approach that bridge has proposed to judge concrete cracks was talked by the Liu Qing of Wuhan University of Technology unit.2006, the Yin Lan of Southeast China University, He Xiaoyuan utilized based on the digital image processing techniques on the flash spotting basis concrete surface crack width characteristics are measured and analyzed.
Yet, more than the method studied mainly be the crack enlarged image of gathering to through contact scanning or shooting at close range.And in the tunnel, for impalpable high-order crack pattern picture, hand-held contact gatherer process is loaded down with trivial details.And the image of wide-long shot has a strong impact on the image imaging quality to noise, light sensitive, makes the follow-up image pre-service grain granite that becomes.The area of crack in the wide-long shot image occupies the ratio of ratio in the picture that contact is gathered and wants much little; The background image that remaining large tracts of land is complicated and changeable; Add the defective of some edge detection methods; Make that the edge extracting method in roomy crack is inapplicable in the recent photograph, need to seek a kind of suitable new method.In addition, measure differently with the contact of fixed distance, the randomness that the position appears in the crack causes the randomness of shooting distance, needs to seek a kind of new scaling method and replace traditional spacing to demarcate.
Therefore, a kind of image detecting method and the detection system that in the tunnel, can measure the FRACTURE CHARACTERISTICS value convenient, quantitatively, quickly and accurately of research become one of pressing for of Tunnel Engineering structure Non-Destructive Testing field.
Summary of the invention
The object of the invention is exactly measuring method and the device that a kind of low cost, quick, convenient, high-precision Tunnel Lining Cracks width are provided for the defective that overcomes above-mentioned prior art existence.
The object of the invention can be realized through following technical scheme:
A kind of measuring method of Tunnel Lining Cracks width is characterized in that, may further comprise the steps:
(1) gathers the image in crack through digital camera, and pass through laser range finder and measure the collection distance;
(2) zone of on the image that collects, select to need measuring, and transfer the area image of choosing to gray level image;
(3) confirm the threshold values of the gray level image that step (2) is chosen, and carry out binary image and cut apart, obtain the target area;
(4) in the target area, extract the white portion that only comprises the crack;
(5) utilization sub-pix method is extracted the edge in crack, obtains sub-pix rank outline map;
(6) pixel wide of employing minimum distance method calculating fracture;
(7), calculate the developed width in crack according to the demarcation ratio of predetermined pixel physical size and photograph distance;
(8) the developed width data that obtain are deposited in the database of crack.
Digital camera then carries out the segmentation collection to large fracture in the described step (1).
The concrete steps of described step (3) are following:
(31) obtain the average gray value of entire image according to following formula:
T ave=∑Pixels/N
Wherein ∑ Pixels is the gray scale summation of each pixel in the image, and N is a total pixel number; With T AveBe made as initial threshold T k, this moment k=1;
(32) according to T 1Image segmentation is become two parts of target and background, with gray scale less than T 1The zone be called the target area, greater than T 1The zone be called the background area, and then obtain two the zone average gray be respectively T LowAnd T High
(33) obtain new threshold value according to following formula:
T k+1=(T low+T high)/2
If T K+1Be not equal to T k, then with T K+1Replace T k, return step (32), make k=k+1 simultaneously, up to T K+1=T k, execution in step (34);
(34) with T K+1Be used as final threshold value, carry out binary image and cut apart, comprise the approximate pocket in crack area and noise spot, gray-scale value and crack in the target area, picture inversion.
The concrete steps of described step (4) are following:
(40) find out white portion maximum in the target area, remainder is treated to black, extract the zone that comprises the crack;
(41) operation is corroded in the zone that extracts of square unit matrix and the step (40) through a n*n, detects the zone that whether exists with its coupling, if yes, execution in step (42), if not, execution in step
(43);
(42) average to carrying out gray scale corresponding to all points in the zone that is extracted in the original image: T j=∑ Pixels Target/ N Target, again according to T jThe image two-value is divided into two parts of target and background; After carrying out picture inversion; Corrode operation once more, judge whether it exists the zone of mating with square unit matrix, and circulate; Till the zone that in image, does not have to mate with square unit matrix, execution in step (43);
(43) the utilization morphological method is carried out open and close computing, cavity filling and is rejected burrs on edges image, seeks out maximum white portion, and all the other zones are treated to black, is only comprised the white portion in crack.
The concrete steps of described step (5) are following:
(51) { M20} carries out convolution algorithm to each pixel of image for M00, M11 through 7*7 Zernike template; To obtain corresponding image Zernike square { Z00, Z11, Z20}; Calculate parameter (φ, h, the l of each pixel through the Zernike square of image; K), to judge whether this pixel is marginal point
The M00 template:
0 0.0287 0.0686 0.0807 0.0686 0.0287 0 0.0287 0.0815 0.0816 0.0816 0.0816 0.0815 0.0287 0.0686 0.0816 0.0816 0.0816 0.0816 0.0816 0.0686 0.0807 0.0816 0.0816 0.0816 0.0816 0.0816 0.0807 0.0686 0.0816 0.0816 0.0816 0.0816 0.0816 0.0868 0.0287 0.0815 0.0816 0.0816 0.0816 0.0815 0.0287 0 0.0287 0.0686 0.0807 0.0686 0.0287 0
M11 real number template:
0 0.0150 0.0190 0 - 0.0190 - 0.0150 0 0.0220 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0220 0.0570 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0570 0.0700 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0700 0.0570 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0570 0.0220 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0220 0 0.0150 0.0190 0 - 0.0190 - 0.0150 0
M11 imaginary number template:
0 - 0.0220 - 0.0570 - 0.0700 - 0.0570 - 0.0220 0 - 0.0150 - 0.0470 - 0.0470 - 0.0470 - 0.0470 - 0.0470 - 0.0150 - 0.0190 - 0.0230 - 0.0230 - 0.0230 - 0.0230 - 0.0230 - 0.0190 0 0 0 0 0 0 0 0.0190 0.0230 0.0230 0.0230 0.0230 0.0230 0.0190 0.0150 0.0470 0.0470 0.0470 0.0470 0.0470 0.0150 0 0.0220 0.0570 0.0700 0.0570 0.0220 0
The M20 template:
0 0.0230 0.0390 0.0410 0.0390 0.0230 0 0.0230 0.0270 - 0.0130 - 0.0260 - 0.0130 0.0270 0.0230 0.0390 - 0.0130 - 0.0530 - 0.0660 - 0.0530 - 0.0130 0.0390 0.0410 - 0.0260 - 0.0660 - 0.0810 - 0.0660 - 0.0260 0.0410 0.0390 - 0.0130 - 0.0530 - 0.0660 - 0.0530 - 0.0130 0.0390 0.0230 0.0270 - 0.0130 - 0.0260 - 0.0130 0.0270 0.0230 0 0.0230 0.0390 0.0410 0.0390 0.0230 0
(52) according to following formula
φ=arctan[Im(Z11)/Re(Z11)]
Calculate angle φ, wherein Im (Z11) and Re (Z11) are respectively imaginary part and the real part of Z11;
(53) according to computes Z ' 11:
Z′11=Re(Z11)cosφ+Im(Z11)sinφ
According to formula l=Z20/Z ' 11, obtain l then;
(54) according to computes step height k:
k=3Z′11/2(1-l 2) 3/2
(55) according to computes background gray scale h:
h = [ Z 00 - kπ 2 + k · arcsin l + kl ( 1 - l 2 ) 1 / 2 ] π ;
(56) obtained the edge parameters of each pixel after, if the parameter of pixel satisfies k>=k t∩ l≤l t, then this pixel is a marginal point, utilizes following formula:
x s=x+l·cos(φ)
y s=y+l·sin(φ)
Calculate the sub-pixel edge point coordinate and obtain final sub-pix rank outline map, wherein k t, l tBe decision threshold.
Described k tBe 0.3, described l tBe gradation of image peaked 1/10th.
Described step (7) concrete steps are following:
(71) try to achieve under the fixed lens focal length through test, when shooting distance was L, the corresponding physical size of each pixel was a in the image;
(72) with L as horizontal ordinate, a is an ordinate, makes calibration curve.
A kind of measurement mechanism of Tunnel Lining Cracks width; It is characterized in that; Comprise digital camera, laser range finder, connecting screw rod, metering computer and crack database; Described digital camera is located at the top of laser range finder through connecting screw rod, and is connected with metering computer through data line, and described crack database is connected with metering computer.
Described digital camera adopts sony α 350 single anti-digital cameras, and the camera lens of this digital camera is fixed as 300mm length.
Compared with prior art; The present invention is different from traditional closely hand-held contact fracture width measuring method; But the technology that digital photographing combines with Digital Image Processing applied to first the measurement of Tunnel Lining Cracks width; Utilize ordinary digital camera remote digital photographic means to collect the crack picture, to some harmful effects aspect Flame Image Process of remote captured picture, propose a kind of Tunnel Lining Cracks measuring method and device again based on Digital photographic and image processing techniques; For miscellaneous tunnel slot surveying work provides a kind of low cost, handled means fast, easily, and obtained satisfactory accuracy.
Description of drawings
Fig. 1 is a process flow diagram of the present invention;
Fig. 2 is a hardware configuration synoptic diagram of the present invention;
Fig. 3 is the synoptic diagram of the present invention through minimum distance method calculating fracture pixel wide;
The desirable step edge model synoptic diagram that Fig. 4 extracts for sub-pixel edge of the present invention.
Embodiment
Below in conjunction with accompanying drawing and specific embodiment the present invention is elaborated.
Embodiment
As shown in Figure 1, a kind of measuring method of Tunnel Lining Cracks width may further comprise the steps:
Step 101 is gathered the image in crack through digital camera, and measures the collection distance through laser range finder;
Step 102, selection needs the zone of measurement on the image that collects, and transfers the area image of choosing to gray level image;
Step 103 is confirmed the threshold values of the gray level image that step 102 is chosen, and carries out binary image and cut apart, and obtains the target area;
Step 104 extracts the white portion that only comprises the crack in the target area;
Step 105, utilization sub-pix method is extracted the edge in crack, obtains sub-pix rank outline map;
Step 106, the pixel wide of employing minimum distance method calculating fracture;
Step 107 according to the demarcation ratio of predetermined pixel physical size and photograph distance, calculates the developed width in crack;
Step 108 deposits the developed width data that obtain in the database of crack in.
The concrete steps of described step 103 are following:
(31) obtain the average gray value of entire image according to following formula:
T ave=∑Pixels/N
Wherein ∑ Pixels is the gray scale summation of each pixel in the image, and N is a total pixel number; With T AveBe made as initial threshold T k, this moment k=1;
(32) according to T 1Image segmentation is become two parts of target and background, with gray scale less than T 1The zone be called the target area, greater than T 1The zone be called the background area, and then obtain two the zone average gray be respectively T LowAnd T High
(33) obtain new threshold value according to following formula:
T k+1=(t low+T high)/2
If T K+1Be not equal to T k, then with T K+1Replace T k, return step (32), make k=k+1 simultaneously, up to T K+1=T k, execution in step (34);
(34) with T K+1Be used as final threshold value, carry out binary image and cut apart, comprise the approximate pocket in crack area and noise spot, gray-scale value and crack in the target area, picture inversion.
The concrete steps of described step 104 are following:
(40) find out white portion maximum in the target area, remainder is treated to black, extract the zone that comprises the crack;
(41) operation is corroded in the zone that extracts of square unit matrix and the step (40) through a n*n, detects the zone that whether exists with its coupling, if yes, execution in step (42), if not, execution in step
(43);
(42) average to carrying out gray scale corresponding to all points in the zone that is extracted in the original image: T j=∑ Pixels Target/ N Target, again according to T jThe image two-value is divided into two parts of target and background; After carrying out picture inversion; Corrode operation once more, judge whether it exists the zone of mating with square unit matrix, and circulate; Till the zone that in image, does not have to mate with square unit matrix, execution in step (43);
(43) the utilization morphological method is carried out open and close computing, cavity filling and is rejected burrs on edges image, seeks out maximum white portion, and all the other zones are treated to black, is only comprised the white portion in crack.
The concrete steps of described step 105 are following:
(51) { M20} carries out convolution algorithm to each pixel of image for M00, M11 through 7*7 Zernike template; To obtain corresponding image Zernike square { Z00, Z11, Z20}; Calculate parameter (φ, h, the l of each pixel through the Zernike square of image; K), to judge whether this pixel is marginal point
The M00 template:
0 0.0287 0.0686 0.0807 0.0686 0.0287 0 0.0287 0.0815 0.0816 0.0816 0.0816 0.0815 0.0287 0.0686 0.0816 0.0816 0.0816 0.0816 0.0816 0.0686 0.0807 0.0816 0.0816 0.0816 0.0816 0.0816 0.0807 0.0686 0.0816 0.0816 0.0816 0.0816 0.0816 0.0868 0.0287 0.0815 0.0816 0.0816 0.0816 0.0815 0.0287 0 0.0287 0.0686 0.0807 0.0686 0.0287 0
M11 real number template:
0 0.0150 0.0190 0 - 0.0190 - 0.0150 0 0.0220 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0220 0.0570 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0570 0.0700 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0700 0.0570 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0570 0.0220 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0220 0 0.0150 0.0190 0 - 0.0190 - 0.0150 0
M11 imaginary number template:
0 - 0.0220 - 0.0570 - 0.0700 - 0.0570 - 0.0220 0 - 0.0150 - 0.0470 - 0.0470 - 0.0470 - 0.0470 - 0.0470 - 0.0150 - 0.0190 - 0.0230 - 0.0230 - 0.0230 - 0.0230 - 0.0230 - 0.0190 0 0 0 0 0 0 0 0.0190 0.0230 0.0230 0.0230 0.0230 0.0230 0.0190 0.0150 0.0470 0.0470 0.0470 0.0470 0.0470 0.0150 0 0.0220 0.0570 0.0700 0.0570 0.0220 0
The M20 template:
0 0.0230 0.0390 0.0410 0.0390 0.0230 0 0.0230 0.0270 - 0.0130 - 0.0260 - 0.0130 0.0270 0.0230 0.0390 - 0.0130 - 0.0530 - 0.0660 - 0.0530 - 0.0130 0.0390 0.0410 - 0.0260 - 0.0660 - 0.0810 - 0.0660 - 0.0260 0.0410 0.0390 - 0.0130 - 0.0530 - 0.0660 - 0.0530 - 0.0130 0.0390 0.0230 0.0270 - 0.0130 - 0.0260 - 0.0130 0.0270 0.0230 0 0.0230 0.0390 0.0410 0.0390 0.0230 0
(52) according to following formula
φ=arctan[Im(Z11)/Re(Z11)]
Calculate angle φ, wherein Im (Z11) and Re (Z11) are respectively imaginary part and the real part of Z11;
(53) according to computes Z ' 11:
Z′11=Re(Z11)cosφ+Im(Z11)sinφ
According to formula l=Z20/Z ' 11, obtain l then;
(54) according to computes step height k:
k=3Z′11/2(1-l 2) 3/2
(55) according to computes background gray scale h:
h = [ Z 00 - kπ 2 + k · arcsin l + kl ( 1 - l 2 ) 1 / 2 ] π ;
(56) obtained the edge parameters of each pixel after, if the parameter of pixel satisfies k>=k t∩ l≤l t, then this pixel is a marginal point, utilizes following formula:
x s=x+l·cos(φ)
y s=y+l·sin(φ)
Calculate the sub-pixel edge point coordinate and obtain final sub-pix rank outline map, wherein k t, l tBe decision threshold.
Described k tBe 0.3, described l tBe gradation of image peaked 1/10th.
Like Fig. 4, (l k) can represent that more intuitively k is a step height through desirable step edge model to the parameter of pixel for φ, h; H is the background gray scale; L is the vertical range of disc centre to the edge; φ is edge and y axle angulation.
Crack of the present invention pixel width gauge is calculated and is adopted minimum distance method, is specially the last lower limb of distinguishing the crack, chooses the each point of coboundary respectively, adopts minimum distance method to calculate the width in target crack.
As shown in Figure 3; Be synoptic diagram through minimum distance method calculating fracture pixel wide; According to the coordinate of marginal point about the vertical direction, begin from the first point of coboundary earlier, the coordinate of using marginal point utilizes the range formula of the point-to-point transmission in the higher mathematics to calculate with each coordinate points of lower limb respectively; Calculate minimum value put the distance of lower limb as this, can be expressed as:
w i = min ( ( x i - x k ) 2 + ( y i - y k ) 2 )
Wherein, k=0,1,2,3......
The i that formulate is got coboundary puts the minimum value and value that lower limb is had a few.
Calculate every of the coboundary distance value to lower limb successively, it is average or ask maximum that these are calculated distance value, is the mean value and the maximal value of fracture width, is expressed as:
Figure BSA00000210453900083
w max=max(w i)。
Described step 107 concrete steps are following:
(71) try to achieve under the fixed lens focal length through test, when shooting distance was L, the corresponding physical size of each pixel was a in the image;
(72) with L as horizontal ordinate, a is an ordinate, makes calibration curve.
The measurement mechanism of Tunnel Lining Cracks width of the present invention; Comprise digital camera 1, laser range finder 2, connecting screw rod 3, metering computer 4 and crack database 5; Described digital camera 1 is located at the top of laser range finder 2 through connecting screw rod 3; And be connected with metering computer 4 through a data line 6, described crack database 5 is connected with metering computer 4.
Described digital camera 1 adopts sony α 350 single anti-digital cameras, and its valid pixel is 1,420 ten thousand, adopts the 70-300mm camera lens, and the camera lens of this digital camera 1 is fixed as 300mm length.
Under the situation of conditions permit, can utilize utility appliance such as illuminating lamp and camera trivets to take, prevent because of light or artificial shake etc. former thereby cause the picture quality of clapping not high, influence the subsequent image processing effect.

Claims (8)

1. the measuring method of a Tunnel Lining Cracks width is characterized in that, may further comprise the steps:
(1) gathers the image in crack through digital camera, and pass through laser range finder and measure the collection distance;
(2) zone of on the image that collects, select to need measuring, and transfer the area image of choosing to gray level image;
(3) confirm the threshold values of the gray level image that step (2) is chosen, and carry out binary image and cut apart, obtain the target area;
(4) in the target area, extract the white portion that only comprises the crack;
(5) utilization sub-pix method is extracted the edge in crack, obtains sub-pix rank outline map;
(6) pixel wide of employing minimum distance method calculating fracture;
(7), calculate the developed width in crack according to the demarcation ratio of predetermined pixel physical size and photograph distance;
(8) the developed width data that obtain are deposited in the database of crack.
2. the measuring method of a kind of Tunnel Lining Cracks width according to claim 1 is characterized in that, the concrete steps of described step (3) are following:
(31) obtain the average gray value of entire image according to following formula:
T ave=∑Pixels/N
Wherein ∑ Pixels is the gray scale summation of each pixel in the image, and N is a total pixel number; With T AveBe made as initial threshold T k, this moment k=l;
(32) according to T 1Image segmentation is become two parts of target and background, with gray scale less than T 1The zone be called the target area, greater than T 1The zone be called the background area, and then obtain two the zone average gray be respectively T LowAnd T High
(33) obtain new threshold value according to following formula:
T k+1=(T low+T high)/2
If T K+1Be not equal to T k, then with T K+1Replace T k, return step (32), make k=k+1 simultaneously, up to T K+1=T k, execution in step (34);
(34) with T K+1Be used as final threshold value, carry out binary image and cut apart, comprise the approximate pocket in crack area and noise spot, gray-scale value and crack in the target area, picture inversion.
3. the measuring method of a kind of Tunnel Lining Cracks width according to claim 1 is characterized in that, the concrete steps of described step (4) are following:
(40) find out white portion maximum in the target area, remainder is treated to black, extract the zone that comprises the crack;
(41) operation is corroded in the zone that extracts of square unit matrix and the step (40) through a n*n, detects the zone that whether exists with its coupling, if yes, execution in step (42), if not, execution in step
(43);
(42) average to carrying out gray scale corresponding to all points in the zone that is extracted in the original image: T j=∑ Pixels Target/ N TargetAccording to Tj the image two-value is divided into two parts of target and background again, carry out picture inversion after, corrode operation once more; Judge whether it exists the zone of mating with square unit matrix; And circulate, till the zone that in image, does not have to mate with square unit matrix, execution in step (43);
(43) the utilization morphological method is carried out open and close computing, cavity filling and is rejected burrs on edges image, seeks out maximum white portion, and all the other zones are treated to black, is only comprised the white portion in crack.
4. the measuring method of a kind of Tunnel Lining Cracks width according to claim 1 is characterized in that, the concrete steps of described step (5) are following:
(51) { M20} carries out convolution algorithm to each pixel of image for M00, M11 through 7*7 Zernike template; To obtain corresponding image Zernike square { Z00, Z11, Z20}; Calculate parameter (φ, h, the l of each pixel through the Zernike square of image; K), to judge whether this pixel is marginal point
The M00 template:
0 0.0287 0.0686 0.0807 0.0686 0.0287 0 0.0287 0.0815 0.0816 0.0816 0.0816 0.0815 0.0287 0.0686 0.0816 0.0816 0.0816 0.0816 0.0816 0.0686 0.0807 0.0816 0.0816 0.0816 0.0816 0.0816 0.0807 0.0686 0.0816 0.0816 0.0816 0.0816 0.0816 0.0868 0.0287 0.0815 0.0816 0.0816 0.0816 0.0815 0.0287 0 0.0287 0.0686 0.0807 0.0686 0.0287 0
M11 real number template:
0 0.0150 0.0190 0 - 0.0190 - 0.0150 0 0.0220 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0220 0.0570 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0570 0.0700 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0700 0.0570 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0570 0.0220 0.0470 0.0230 0 - 0.0230 - 0.0470 - 0.0220 0 0.0150 0.0190 0 - 0.0190 - 0.0150 0
M11 imaginary number template:
0 - 0.0220 - 0.0570 - 0.0700 - 0.0570 - 0.0220 0 - 0.0150 - 0.0470 - 0.0470 - 0.0470 - 0.0470 - 0.0470 - 0.0150 - 0.0190 - 0.0230 - 0.0230 - 0.0230 - 0.0230 - 0.0230 - 0.0190 0 0 0 0 0 0 0 0.0190 0.0230 0.0230 0.0230 0.0230 0.0230 0.0190 0.0150 0.0470 0.0470 0.0470 0.0470 0.0470 0.0150 0 0.0220 0.0570 0.0700 0.0570 0.0220 0
The M20 template:
0 0.0230 0.0390 0.0410 0.0390 0.0230 0 0.0230 0.0270 - 0.0130 - 0.0260 - 0.0130 0.0270 0.0230 0.0390 - 0.0130 - 0.0530 - 0.0660 - 0.0530 - 0.0130 0.0390 0.0410 - 0.0260 - 0.0660 - 0.0810 - 0.0660 - 0.0260 0.0410 0.0390 - 0.0130 - 0.0530 - 0.0660 - 0.0530 - 0.0130 0.0390 0.0230 0.0270 - 0.0130 - 0.0260 - 0.0130 0.0270 0.0230 0 0.0230 0.0390 0.0410 0.0390 0.0230 0
(52) according to following formula
φ=arctan[Im(Z11)/Re(Z11)]
Calculate angle φ, wherein Im (Z11) and Re (Z11) are respectively imaginary part and the real part of Z11;
(53) according to computes Z ' 11:
Z′11=Re(Z11)cosφ+Im(Z11)sinφ
According to formula l=Z20/Z ' 11, obtain l then;
(54) according to computes step height k:
k=3Z′11/2(1-l 2) 3/2
(55) according to computes background gray scale h:
h = [ Z 00 - kπ 2 + k · arcsin l + kl ( 1 - l 2 ) 1 / 2 ] π ;
(56) obtained the edge parameters of each pixel after, if the parameter of pixel satisfies k>=k t∩ l≤l t, then this pixel is a marginal point, utilizes following formula:
x s=x+l·cos(φ)
y s=y+l·sin(φ)
Calculate the sub-pixel edge point coordinate and obtain final sub-pix rank outline map, wherein k t, l tBe decision threshold.
5. the measuring method of a kind of Tunnel Lining Cracks width according to claim 4 is characterized in that, described k tBe 0.3, described l tBe gradation of image peaked 1/10th.
6. the measuring method of a kind of Tunnel Lining Cracks width according to claim 1 is characterized in that, described step (7) concrete steps are following:
(71) try to achieve under the fixed lens focal length through test, when shooting distance was L, the corresponding physical size of each pixel was a in the image;
(72) with L as horizontal ordinate, a is an ordinate, makes calibration curve.
7. the measurement mechanism of a Tunnel Lining Cracks width; It is characterized in that; Comprise digital camera, laser range finder, connecting screw rod, metering computer and crack database; Described digital camera is located at the top of laser range finder through connecting screw rod, and is connected with metering computer through data line, and described crack database is connected with metering computer.
8. the measurement mechanism of a kind of Tunnel Lining Cracks width according to claim 1 is characterized in that, described digital camera adopts sony α 350 single anti-digital cameras, and the camera lens of this digital camera is fixed as 300mm length.
CN2010102405432A 2010-07-29 2010-07-29 Tunnel lining crack width measuring method and device Pending CN102346013A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102405432A CN102346013A (en) 2010-07-29 2010-07-29 Tunnel lining crack width measuring method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102405432A CN102346013A (en) 2010-07-29 2010-07-29 Tunnel lining crack width measuring method and device

Publications (1)

Publication Number Publication Date
CN102346013A true CN102346013A (en) 2012-02-08

Family

ID=45544907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102405432A Pending CN102346013A (en) 2010-07-29 2010-07-29 Tunnel lining crack width measuring method and device

Country Status (1)

Country Link
CN (1) CN102346013A (en)

Cited By (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102679914A (en) * 2012-05-14 2012-09-19 同济大学 Method and device for measuring percolating water area of shield tunnel lining segment
CN102768022A (en) * 2012-07-20 2012-11-07 长安大学 Tunnel surrounding rock deformation detection method adopting digital camera technique
CN103268409A (en) * 2013-05-14 2013-08-28 长安大学 Method for measuring fissure widths of tunnel lining rocks
CN103486971A (en) * 2013-08-14 2014-01-01 北京交通大学 Subway tunnel crack width detecting and correcting algorithm
CN103528527A (en) * 2013-10-15 2014-01-22 北京交通大学长三角研究院 Area selection-based automatic crack size measurement method
CN103673858A (en) * 2013-12-12 2014-03-26 同济大学 Dynamic measurement device for shield tunnel lining segment joint wrong lifting
CN103776839A (en) * 2014-02-10 2014-05-07 宁波工程学院 Novel surface crack detection algorithm
CN104165892A (en) * 2014-07-25 2014-11-26 中国科学院东北地理与农业生态研究所 Soil crack parameter online measuring system and method for realizing soil crack parameter extraction by system
CN105606030A (en) * 2015-12-30 2016-05-25 浙江建设职业技术学院 Building surface crack detection device based on embedded mode
CN105741285A (en) * 2016-01-28 2016-07-06 北京工业大学 Image processing device used for tunnel diseases
CN105765374A (en) * 2014-07-25 2016-07-13 西日本高速道路工程四国株式会社 Tunnel lining surface examination system and vehicle used in tunnel lining surface examination system
CN106504246A (en) * 2016-11-08 2017-03-15 太原科技大学 The image processing method of tunnel slot detection
CN107747907A (en) * 2017-09-28 2018-03-02 中国十七冶集团有限公司 A kind of detection method in underground pipe gallery crack
CN107894384A (en) * 2017-11-10 2018-04-10 石家庄铁道大学 Hydrous fluids crack rock tunnel-liner distribution of water pressure Experimental mimic system
CN108596869A (en) * 2018-01-08 2018-09-28 东北大学 A kind of Tunnel Lining Cracks rapid detection method
CN109341542A (en) * 2018-10-29 2019-02-15 福建省智能养护工程有限公司 Method and its monitoring device based on digital signal sequences length identification fracture width
CN109374644A (en) * 2018-11-09 2019-02-22 中国矿业大学 Tunnel-liner defect intelligent measurement imitative experimental appliance based on image recognition
CN109444150A (en) * 2018-11-06 2019-03-08 西北工业大学 Contactless crackle measurement method and its device
CN109579709A (en) * 2018-12-25 2019-04-05 陕西圆周率文教科技有限公司 A kind of Immovable Cultural Relics Crack Monitoring device and method
CN109632103A (en) * 2018-11-22 2019-04-16 西安理工大学 High vacant building Temperature Distribution and surface crack remote supervision system and monitoring method
CN109631771A (en) * 2018-12-29 2019-04-16 西安联创兴科测控科技有限公司 Crack dynamic change sensor and crack measurement method based on the processing of more target images
CN110082354A (en) * 2019-03-29 2019-08-02 同济大学 A kind of potential peeling removing infrared detection device of tunnel-liner shallow-layer
CN110207592A (en) * 2019-04-15 2019-09-06 深圳高速工程检测有限公司 Building cracks measurement method, device, computer equipment and storage medium
CN110298802A (en) * 2019-06-28 2019-10-01 广东诚泰交通科技发展有限公司 Road surface crack detection method and device
CN110533713A (en) * 2019-07-10 2019-12-03 湖南交工智能技术有限公司 Bridge Crack width high-precision measuring method and measuring device
CN110660056A (en) * 2019-10-10 2020-01-07 昆山市建设工程质量检测中心 Building crack width measurement algorithm and method based on image processing
CN110738646A (en) * 2019-10-12 2020-01-31 北京城市排水集团有限责任公司 automatic identification method for pipeline breakage
CN110926342A (en) * 2019-11-27 2020-03-27 北京工业大学 Crack width measuring method and device
CN111102932A (en) * 2019-12-13 2020-05-05 中国二十冶集团有限公司 Automatic inspection method and system for foundation pit safety
CN111739043A (en) * 2020-04-13 2020-10-02 北京京东叁佰陆拾度电子商务有限公司 Parking space drawing method, device, equipment and storage medium
CN112560587A (en) * 2020-11-27 2021-03-26 贵州中建建筑科研设计院有限公司 Dynamic early warning method and system for convolutional neural network slope crack change
CN113076959A (en) * 2021-04-27 2021-07-06 中国矿业大学 Concrete structure surface microcrack feature extraction method based on convolutional neural network
CN113252700A (en) * 2021-07-01 2021-08-13 湖南大学 Structural crack detection method, equipment and system
CN114353670A (en) * 2022-01-07 2022-04-15 深圳市北斗云信息技术有限公司 Method for measuring cracks through visual multi-parameter intelligent crack meter
CN116739994A (en) * 2023-05-18 2023-09-12 山东新宝龙工业科技有限公司 Digital conveyer belt graphic data monitoring system with early warning function

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006166212A (en) * 2004-12-09 2006-06-22 Seiko Epson Corp Image processing method, system, apparatus and program, and digital camera
JP2008139285A (en) * 2006-11-29 2008-06-19 Korea Research Inst Of Standards & Science Construct using picture processing technique, and crack width measuring method of product
CN101706259A (en) * 2009-11-25 2010-05-12 天津大学 Concrete crack width test method based on wavefront coding technology and hand-held tester
CN101957178A (en) * 2009-07-17 2011-01-26 上海同岩土木工程科技有限公司 Method and device for measuring tunnel lining cracks

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006166212A (en) * 2004-12-09 2006-06-22 Seiko Epson Corp Image processing method, system, apparatus and program, and digital camera
JP2008139285A (en) * 2006-11-29 2008-06-19 Korea Research Inst Of Standards & Science Construct using picture processing technique, and crack width measuring method of product
CN101957178A (en) * 2009-07-17 2011-01-26 上海同岩土木工程科技有限公司 Method and device for measuring tunnel lining cracks
CN101706259A (en) * 2009-11-25 2010-05-12 天津大学 Concrete crack width test method based on wavefront coding technology and hand-held tester

Cited By (53)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102679914B (en) * 2012-05-14 2014-08-13 同济大学 Method and device for measuring percolating water area of shield tunnel lining segment
CN102679914A (en) * 2012-05-14 2012-09-19 同济大学 Method and device for measuring percolating water area of shield tunnel lining segment
CN102768022A (en) * 2012-07-20 2012-11-07 长安大学 Tunnel surrounding rock deformation detection method adopting digital camera technique
CN102768022B (en) * 2012-07-20 2015-04-15 长安大学 Tunnel surrounding rock deformation detection method adopting digital camera technique
CN103268409A (en) * 2013-05-14 2013-08-28 长安大学 Method for measuring fissure widths of tunnel lining rocks
CN103486971B (en) * 2013-08-14 2016-07-13 北京交通大学 A kind of subway tunnel fracture width detection and correcting algorithm
CN103486971A (en) * 2013-08-14 2014-01-01 北京交通大学 Subway tunnel crack width detecting and correcting algorithm
CN103528527A (en) * 2013-10-15 2014-01-22 北京交通大学长三角研究院 Area selection-based automatic crack size measurement method
CN103673858A (en) * 2013-12-12 2014-03-26 同济大学 Dynamic measurement device for shield tunnel lining segment joint wrong lifting
CN103776839A (en) * 2014-02-10 2014-05-07 宁波工程学院 Novel surface crack detection algorithm
CN103776839B (en) * 2014-02-10 2016-05-04 湖州师范学院 A kind of Surface Crack Inspection Algorithm
CN105765374A (en) * 2014-07-25 2016-07-13 西日本高速道路工程四国株式会社 Tunnel lining surface examination system and vehicle used in tunnel lining surface examination system
CN105765374B (en) * 2014-07-25 2020-04-10 西日本高速道路工程四国株式会社 Tunnel lining surface examination system and vehicle used for tunnel lining surface examination system
CN104165892B (en) * 2014-07-25 2016-10-26 中国科学院东北地理与农业生态研究所 Soil Crack Parameters on-line measurement system and use this system to realize soil Crack Parameters extracting method
CN104165892A (en) * 2014-07-25 2014-11-26 中国科学院东北地理与农业生态研究所 Soil crack parameter online measuring system and method for realizing soil crack parameter extraction by system
CN105606030B (en) * 2015-12-30 2018-08-24 浙江建设职业技术学院 One kind being based on Embedded building surface crack detecting device
CN105606030A (en) * 2015-12-30 2016-05-25 浙江建设职业技术学院 Building surface crack detection device based on embedded mode
CN105741285A (en) * 2016-01-28 2016-07-06 北京工业大学 Image processing device used for tunnel diseases
CN105741285B (en) * 2016-01-28 2018-09-04 北京工业大学 Image processing apparatus for tunnel defect
CN106504246A (en) * 2016-11-08 2017-03-15 太原科技大学 The image processing method of tunnel slot detection
CN106504246B (en) * 2016-11-08 2019-04-30 太原科技大学 The image processing method of tunnel slot detection
CN107747907A (en) * 2017-09-28 2018-03-02 中国十七冶集团有限公司 A kind of detection method in underground pipe gallery crack
CN107894384A (en) * 2017-11-10 2018-04-10 石家庄铁道大学 Hydrous fluids crack rock tunnel-liner distribution of water pressure Experimental mimic system
CN108596869A (en) * 2018-01-08 2018-09-28 东北大学 A kind of Tunnel Lining Cracks rapid detection method
CN108596869B (en) * 2018-01-08 2020-06-16 东北大学 Rapid detection method for tunnel lining cracks
CN109341542A (en) * 2018-10-29 2019-02-15 福建省智能养护工程有限公司 Method and its monitoring device based on digital signal sequences length identification fracture width
CN109444150A (en) * 2018-11-06 2019-03-08 西北工业大学 Contactless crackle measurement method and its device
CN109374644A (en) * 2018-11-09 2019-02-22 中国矿业大学 Tunnel-liner defect intelligent measurement imitative experimental appliance based on image recognition
CN109374644B (en) * 2018-11-09 2021-08-03 中国矿业大学 Tunnel lining defect intelligent detection simulation experiment device based on image recognition
CN109632103A (en) * 2018-11-22 2019-04-16 西安理工大学 High vacant building Temperature Distribution and surface crack remote supervision system and monitoring method
CN109632103B (en) * 2018-11-22 2020-02-14 西安理工大学 High-altitude building temperature distribution and surface crack remote monitoring system and monitoring method
CN109579709A (en) * 2018-12-25 2019-04-05 陕西圆周率文教科技有限公司 A kind of Immovable Cultural Relics Crack Monitoring device and method
CN109631771A (en) * 2018-12-29 2019-04-16 西安联创兴科测控科技有限公司 Crack dynamic change sensor and crack measurement method based on the processing of more target images
CN109631771B (en) * 2018-12-29 2024-01-16 西安联创兴科测控科技有限公司 Crack dynamic change sensor based on multi-target image processing and crack measuring method
CN110082354A (en) * 2019-03-29 2019-08-02 同济大学 A kind of potential peeling removing infrared detection device of tunnel-liner shallow-layer
CN110207592B (en) * 2019-04-15 2021-11-09 深圳高速工程检测有限公司 Building crack measuring method and device, computer equipment and storage medium
CN110207592A (en) * 2019-04-15 2019-09-06 深圳高速工程检测有限公司 Building cracks measurement method, device, computer equipment and storage medium
CN110298802A (en) * 2019-06-28 2019-10-01 广东诚泰交通科技发展有限公司 Road surface crack detection method and device
CN110298802B (en) * 2019-06-28 2021-07-02 广东诚泰交通科技发展有限公司 Pavement crack detection method and device
CN110533713A (en) * 2019-07-10 2019-12-03 湖南交工智能技术有限公司 Bridge Crack width high-precision measuring method and measuring device
CN110660056B (en) * 2019-10-10 2020-08-18 昆山市建设工程质量检测中心 Building crack width measuring method based on image processing
CN110660056A (en) * 2019-10-10 2020-01-07 昆山市建设工程质量检测中心 Building crack width measurement algorithm and method based on image processing
CN110738646A (en) * 2019-10-12 2020-01-31 北京城市排水集团有限责任公司 automatic identification method for pipeline breakage
CN110926342A (en) * 2019-11-27 2020-03-27 北京工业大学 Crack width measuring method and device
CN110926342B (en) * 2019-11-27 2021-07-23 北京工业大学 Crack width measuring method and device
CN111102932A (en) * 2019-12-13 2020-05-05 中国二十冶集团有限公司 Automatic inspection method and system for foundation pit safety
CN111739043A (en) * 2020-04-13 2020-10-02 北京京东叁佰陆拾度电子商务有限公司 Parking space drawing method, device, equipment and storage medium
CN111739043B (en) * 2020-04-13 2023-08-08 北京京东叁佰陆拾度电子商务有限公司 Parking space drawing method, device, equipment and storage medium
CN112560587A (en) * 2020-11-27 2021-03-26 贵州中建建筑科研设计院有限公司 Dynamic early warning method and system for convolutional neural network slope crack change
CN113076959A (en) * 2021-04-27 2021-07-06 中国矿业大学 Concrete structure surface microcrack feature extraction method based on convolutional neural network
CN113252700A (en) * 2021-07-01 2021-08-13 湖南大学 Structural crack detection method, equipment and system
CN114353670A (en) * 2022-01-07 2022-04-15 深圳市北斗云信息技术有限公司 Method for measuring cracks through visual multi-parameter intelligent crack meter
CN116739994A (en) * 2023-05-18 2023-09-12 山东新宝龙工业科技有限公司 Digital conveyer belt graphic data monitoring system with early warning function

Similar Documents

Publication Publication Date Title
CN101957178B (en) Method and device for measuring tunnel lining cracks
CN102346013A (en) Tunnel lining crack width measuring method and device
CN111855664B (en) Adjustable three-dimensional tunnel defect detection system
CN106969706A (en) Workpiece sensing and three-dimension measuring system and detection method based on binocular stereo vision
CN102297660B (en) Measuring method of shield tunnel lining segment seam open width and apparatus thereof
CN103928363B (en) A kind of silicon chip crystal orientation detection method and detection device
CN113252700B (en) Structural crack detection method, equipment and system
CN104173054A (en) Measuring method and measuring device for height of human body based on binocular vision technique
CN104034733A (en) Service life prediction method based on binocular vision monitoring and surface crack image recognition
CN109559348A (en) A kind of contactless deformation measurement method of bridge based on tracing characteristic points
CN102680480A (en) Intelligent detecting method for cracks of concrete structures
CN104700395A (en) Method and system for detecting appearance crack of structure
CN106017325A (en) Improved non-contact optical measurement method for complex surface and irregular object volume
CN111612846A (en) Concrete crack width measuring method based on U-net CNN image recognition and pixel calibration
CN102359966A (en) Positioning system for micro-cracks on concrete surfaces
CN106677037B (en) Portable asphalt pavement disease detection method and device based on machine vision
CN102635056A (en) Measuring method for construction depth of asphalt road surface
EP3023912A1 (en) Crack data collection apparatus and server apparatus to collect crack data
CN104501735A (en) Method for observing three-dimensional deformation of side slope by utilizing circular marking points
CN106841214A (en) A kind of non-contact wind power blade dust storm erosion degree detection method
CN107238375A (en) Detect the one camera photogrammetric survey method of bridge rubber bearing displacement
CN204405556U (en) A kind of glue into concrete beam cracks pick-up unit based on Digital Image Processing
CN112254675B (en) Space occupancy rate acquisition and judgment equipment and method containing moving object
CN113884011A (en) Non-contact concrete surface crack measuring equipment and method
CN111473776A (en) Landslide crack monitoring method based on single-image close-range photogrammetry

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120208