Method for Automatically Identifying Ring Joint of Shield Tunnel Based on Lining Structure

Abstract

The present disclosure provides a method for automatically identifying a ring joint of a shield tunnel based on a lining structure. The method includes the following steps: S1: acquiring a three-dimensional (3D) point cloud of a shield tunnel through a mobile scanning system: S2: generating an orthographic projection image of an inner wall of the tunnel: S3: identifving a feature of a bolt hole; S4: extracting a longitudinal joint of the shield tunnel: S5: generating a prior structural template ring; and S6: extracting a transverse joint of the shield tunnel. The present disclosure has the following advantages. Starting from the features of the lining structure of the shield tunnel, the present disclosure selects a bolt hole with a strong structural feature, takes the structural feature of the bolt hole as an identification feature, and indirectly extracts joint information of straight and staggered joints tunnel.

Claims

1. A method for automatically identifying a ring joint of a shield tunnel based on a lining structure, comprising the following steps: S1: acquiring a three-dimensional (3D) point cloud of a shield tunnel through a mobile scanning system; S2: generating an orthographic projection image of an inner wall of the tunnel; S3: identifying a feature of a bolt hole; S4: extracting a longitudinal joint of the shield tunnel; S5: generating a prior structural template ring; and S6: extracting a transverse joint of the shield tunnel.

2. The method for automatically identifying a ring joint of a shield tunnel based on a lining structure according to claim 1, wherein in step S1, the point cloud of the shield tunnel is acquired by a mobile scanning system, and the point cloud comprises 3D coordinate position information and intensity information.

3. The method for automatically identifying a ring joint of a shield tunnel based on a lining structure according to claim 1, wherein in step S2, a cylindrical projection model is used to perform orthographic projection of the shield tunnel to generate an orthographic image of the inner wall of the tunnel.

4. The method for automatically identifying a ring joint of a shield tunnel based on a lining structure according to claim 1, wherein step S3 comprises the following sub-steps: S31: selecting a bolt hole region in a tunnel image; sampling along a tunnel mileage; calculating a distance from a corresponding point to the center of a cross-section fitting ellipse; selecting a maximum distance from a same cross-section sampling point set as a current cross-section sampling distance to compose a sampling point set G, so as to eliminate an impact of the obstruction of an auxiliary facility; S32: taking a design width t_w and a depth t_d of a bolt hole as thresholds to perform a clustering segmentation algorithm on points in the sampling point set G, and identifying all t clusters to form a cluster centroid point set J; and S33: calculating a mean k of all identified cluster centroids in a sliding window with a width of 6; taking points with a distance less than k in the window as tunnel wall points J_2 and points with a distance greater than k as bolt hole points J_1; composing all J_1 into a bolt hole point set H, wherein the point set J comprises two types of points, namely bolt hole cluster points J_1 and tunnel wall points J_2.

5. The method for automatically identifying a ring joint of a shield tunnel based on a lining structure according to claim 1, wherein, step S4 comprises: composing all theoretical joint mileage positions into l; traversing in the bolt hole point set H to select a point H.sub.i; traversing in l to find a point l.sub.i closest to H.sub.i; putting H.sub.i into a point set H.sub.left if H.sub.i<l.sub.i; putting H.sub.i into a point set right if H.sub.i>l.sub.i, and taking a closest pair of points p.sub.l and p.sub.r from H.sub.leftt and H.sub.right, to obtain a current longitudinal joint position p.sub.h of the shield tunnel: $p_h=\frac{p_l+p_r}{2}.$

6. The method for automatically identifying a ring joint of a shield tunnel based on a lining structure according to claim 1, wherein step S5 comprises the following sub-steps: S51: manually selecting joint positions of 1 to 3 rings as prior position information according to an actual situation of the tunnel, and comprising joint positions of an i-th ring into a point set O.sub.i; S52: extracting a bolt hole point set H.sub.i in the i-th ring through the algorithm in step S3: traversing bolt holes in H.sub.i, and finding a joint closest to a current bolt hole in the point set O.sub.i: storing a current joint-bolt hole positional relationship index h-oi into a positional relationship index set HO.sub.i; and S53: taking a union of the HO.sub.i of all rings to obtain an overall prior structural template ring set HO.

7. The method for automatically identifying a ring joint of a shield tunnel based on a lining structure according to claim 1, wherein step S6 comprises: rotating a to-be-identified ring i by a rotation angle θ of 0° to 360°: extracting a bolt hole point set H.sub.i of the to-be-identified ring i by the algorithm in step S3: traversing the bolt hole set in the structural template ring HO: finding, by matching, bolt holes with a smallest azimuth angle difference in the structural template ring HO corresponding to each bolt hole in H.sub.i: calculating an average angle difference δ under a current rotation angle θ: adding the δ under all values of the θ to an average angle difference set ω; and selecting a smallest average angle difference in the set ω, and obtaining a corresponding rotation angle θ.sub.min; rotating the to-be-identified ring i by θ.sub.min; traversing in the prior structural template ring HO to find a template bolt hole that is closest to each bolt hole in H.sub.i, and directly obtaining a corresponding transverse joint position: taking a mean as a final transverse joint position p_z if there is a repeated joint position.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0026] FIG. 1 is a flowchart of a method according to the present disclosure.

[0027] FIG. 2 shows a structural feature according to an embodiment of the present disclosure.

[0028] FIG. 3 shows identification of a longitudinal joint according to an embodiment of the present disclosure.

[0029] FIG. 4 shows generation of a template ring according to an embodiment of the present disclosure.

[0030] FIG. 5 shows a transverse joint matching curve according to an embodiment of the present disclosure.

[0031] FIG. 6 shows segmentation of cross-section point segments according to an embodiment of the present disclosure.

DETAILED DESCRIPTION

[0032] The present disclosure is described in further detail with reference to the preferred embodiments and accompanying drawings (FIGS. 1 to 6) of the present disclosure.

[0033] FIG. 1 is a flowchart of a method according to the present disclosure. The method includes the following steps:

[0034] S1: Acquire a three-dimensional (3D) point cloud of a shield tunnel through a mobile scanning system.

[0035] S2: Generate an orthographic projection image of an inner wall of the tunnel.

[0036] S3: Identify a feature of a bolt hole.

[0037] S4: Extract a longitudinal joint of the shield tunnel.

[0038] S5: Generate a prior structural template ring.

[0039] S6: Extract a transverse joint of the shield tunnel.

[0040] Further, in step S2, a cylindrical projection model is used to perform orthographic projection of the shield tunnel to generate an orthographic image of the inner wall of the tunnel, which is used for manual prior selection of a joint and verification of a joint identification result.

[0041] FIG. 2 shows a feature of a bolt hole that needs to be identified and transverse and longitudinal joints that need to be extracted in the present disclosure.

[0042] Step S3 is described in detail with reference to parameters shown in FIG. 3, and includes the following sub-steps:

[0043] S31: Select a bolt hole region in a tunnel image; sample along a tunnel mileage; calculate a distance from a corresponding point to the center of a cross-section fitting ellipse; select a maximum distance from a same cross-section sampling point set as a current cross-section sampling distance to compose a sampling point set G, which is a curve composed of original points in FIG. 3, so as to eliminate an impact of the obstruction of an auxiliary facility.

[0044] S32: Take a design width t_w and a depth t_d of a bolt hole as thresholds to perform a clustering segmentation algorithm on points in the sampling point set G, and identify all t clusters to form a cluster centroid point set J.

[0045] S33: Calculate a mean k of all identified cluster centroids in a sliding window with a width of 6; take points with a distance less than k in the window as tunnel wall points J_2 and points with a distance greater than k as bolt hole points J_1; compose all J_1 into a bolt hole point set H. As shown in FIG. 3, the point set J includes two types of points, namely bolt hole cluster points J_1 and tunnel wall points J_2.

[0046] Step S4 is described in detail with reference to parameters shown in FIG. 3. It includes: compose all theoretical joint mileage positions into l; traverse in the bolt hole point set H to select a point H.sub.i; traverse in 1 to find a point l.sub.i closest to H.sub.i; put H.sub.i into a point set H.sub.left if H.sub.i<l.sub.i; put H.sub.i into a point set H.sub.right if H.sub.i>l.sub.i, and take a closest pair of points p.sub.l and p.sub.r from H.sub.left and H.sub.right (as shown in FIG. 3), to obtain a current longitudinal joint position p.sub.h of the shield tunnel:

[00002]$p_h=\frac{p_l+p_r}{2}$

(indicated in the dashed box in FIG. 3)

[0047] Step S5 is described in detail with reference to parameters shown in FIG. 4. The reference numerals 1 to 6 in FIG. 4 are respectively as follows: 1. Bolt hole, 2. Transverse joint, 3. Prior joint selection, 4. Joint-bolt hole correspondence, 5. Prior joint union, 6. Prior structural template ring. This step specifically includes the following sub-steps:

[0048] S51: Select joint positions of two rings as prior position information (as shown in “3. Prior joint selection” in FIG. 4), and compose joint positions of an i-th ring into a point set O.sub.i.

[0049] S52: Extract a bolt hole (“1. Bolt hole” in FIG. 4) in the i-th ring through the algorithm in step S3 to compose a point set H.sub.i; traverse bolt holes in H.sub.i, and find a joint (“2. Transverse joint” in FIG. 4) closest to a current bolt hole in the point set O.sub.i; store a current joint-bolt hole positional relationship index h-oi (“4. Joint-bolt hole correspondence” in FIG. 4) into a positional relationship index set HO.sub.i.

[0050] S53: Take a union of the HO.sub.i of all rings (“5. Prior joint union” shown in FIG. 4) to obtain an overall prior structural template ring set HO (“6. Prior structural template ring” shown in FIG. 4).

[0051] Step S6 is described with reference to FIG. 5. In this embodiment, this step includes: rotate a to-be-identified ring i by a rotation angle θ of 0° to 360°; extract a bolt hole point set H.sub.i of the to-be-identified ring i by the algorithm in step S3; traverse the bolt hole set in the structural template ring HO; find, by matching, bolt holes with a smallest azimuth angle difference in the structural template ring HO corresponding to each bolt hole in H.sub.i; calculate an average angle difference δ under a current rotation angle θ; draw a curve of the average angle difference with the rotation angle (as shown in FIG. 5); add the δ under all values of the θ to an average angle difference set φ; select a smallest average angle difference in the set φ, and obtain a corresponding rotation angle θ.sub.min (indicated by “Min” in FIG. 5); finally, rotate the to-be-identified ring i by θ.sub.min; traverse in the prior structural template ring HO to find a template bolt hole that is closest to each bolt hole in H.sub.i, and directly obtain a corresponding transverse joint position; take a mean as a final transverse joint position p_z if there is a repeated joint position.

[0052] Under the premise of accurate identification of the positions of the transverse and longitudinal joints, the point cloud of the tunnel ring is divided into multiple segment point clouds (as shown in FIG. 6), and an origin of a cross-section coordinate system is the center of an ellipse for fitting a cross-section point cloud. Each segment can be regarded as an arc-shaped rigid body. In order to eliminate noise and misalignment calculations, the point cloud of each segment is fitted with a least squares circle against gross error, to obtain the center (a,b), radius r and fitting accuracy a of the fitting circle of each segment.

[0053] Misalignment refers to an uneven deformation between adjacent segments of the lining ring. The intra-ring misalignment and the inter-ring misalignment are calculated by calculating the difference in a space distance between the points on the segments i and j on both sides of the joint. Therefore, it is necessary to calculate the coordinates of the point p on the arc of the segment when the azimuth angle θ is given in the cross-section coordinate system:

x.sub.p=r*sin(θ)+a

y.sub.p=r*cos(θ)+b

[0054] The misalignment between two adjacent segments is a distance between two points:

δ.sub.ij=∥p.sub.i−p.sub.j∥

[0055] The calculation accuracy of the misalignment can be obtained by the circle fitting accuracy of the corresponding segment:

σ.sub.ij=√{square root over (σ.sub.i.sup.2+σ.sub.j.sup.2)}

[0056] The intra-ring misalignment is calculated by adjacent segments in the ring, and the inter-ring misalignment is calculated by adjacent segments between the rings. Because there are multiple scan lines in a single ring and the fitting accuracy of different scan lines is different, the calculation accuracy of the misalignment is also different. Therefore, the scan lines in the same ring can be sampled multiple times to calculate the amount and accuracy of the misalignment respectively, and the misalignment value with the highest accuracy is selected as the final misalignment calculation result.