THREE-DIMENSIONAL (3D) TERRAIN RECONSTRUCTION METHOD FOR SCOURED AREA AROUND BRIDGE PIER FOUNDATION BASED ON MECHANICAL SCANNED IMAGING SONAR

Abstract

A three-dimensional (3D) terrain reconstruction method for a scoured area around bridge pier foundation based on a mechanical scanned imaging sonar includes scanning an overall terrain of a scoured area around bridge pier foundation with a sonar from different azimuths to acquire n sonar images of a foundation scouring terrain; intercepting multiple analysis sections from each of acquired sonar images at a same distance; extracting key parameters of upper and lower edges on a terrain imaging strip in each of the analysis sections in the image, and transforming extracted parameters to a 3D space, a fan-shaped beam surface of the sonar being represented with a fan-shaped arc; recognizing a scour terrain profile in the analysis section; recognizing terrain profiles one by one, and respectively extracting spatially scattered 3D coordinate data; and performing interpolation and fitting on the spatially scattered data, thus implementing 3D reconstruction for the foundation scouring terrain.

Claims

1. A three-dimensional (3D) terrain reconstruction method for a scoured area around bridge pier foundation based on a mechanical scanned imaging sonar, comprising the following steps: step S1: scanning an overall terrain of a scoured area around bridge pier foundation with a mechanical scanning sonar device from different azimuths to acquire n sonar images of a scoured terrain in the different azimuths; step S2: intercepting multiple analysis sections from each of acquired sonar images at a same distance; step S3: extracting key parameters of upper and lower edges on a terrain imaging strip in each of the analysis sections in the image, and transforming extracted parameters to a 3D space according to an acoustic imaging principle, a fan-shaped beam surface of the sonar being represented with a fan-shaped arc; step S4: recognizing, according to a relative positional relationship of fan-shaped arcs for a group of measuring points in a same analysis section, a scour terrain profile in the analysis section; step S5: recognizing terrain profiles in different azimuths and different analysis sections one by one, and respectively extracting 3D coordinates of key points on each of profile curves to obtain required spatially scattered 3D coordinate data; and step S6: performing interpolation and fitting on the spatially scattered data with a Kriging interpolation method, thus implementing 3D reconstruction for the foundation scouring terrain.

2. The 3D terrain reconstruction method for a scoured area around bridge pier foundation based on a mechanical scanned imaging sonar according to claim 1, wherein the extracting key parameters of upper and lower edges on a terrain imaging strip in each of the analysis sections in the sonar image, and transforming extracted parameters to a 3D space according to an acoustic imaging principle in step S3 specifically comprises: step 3a: labeling, in an analysis section p.sub.j,m.sup.i at a distance d.sub.j.sup.i away from a sonar center, upper and lower edge points M, N on an imaging strip, ranges r.sub.jk.sup.i, R.sub.jk.sup.i from an original sonar center point P.sub.jk.sup.i to the upper and lower edge points M, N, and a large opening angle θ of the fan-shaped beam surface; step 3b: drawing a circular arc with the upper edge point M as a center and the r.sub.jk.sup.i as a radius, an intersection with the analysis section p.sub.j,m.sup.i being an equivalent sonar center point Q.sub.jk.sup.i, and the equivalent sonar center point Q.sub.jk.sup.i being T.sub.jk.sup.i higher than a plane where the original sonar center is located; drawing a fan-shaped arc with the equivalent sonar center point Q.sub.jk.sup.i that is T.sub.jk.sup.i higher than the plane where the sonar is located as a center, the r.sub.jk.sup.i as a radius, and the large opening angle θ of the fan-shaped beam surface to obtain a corresponding fan-shaped arc of the upper edge point M in the 3D space; drawing a circular arc with the lower edge point N as a center and the as a radius, an intersection with the analysis section being an equivalent sonar center point that is higher than the plane where the original sonar center is located; and drawing a fan-shaped arc with the as a center, the as a radius and the angle θ to obtain a corresponding fan-shaped arc of the lower edge point N in the 3D space, where the superscript i represents a pier No., the first subscript j represents an azimuth in which a group of measuring points are located relative to a bridge pier, the second subscript k represents a specific measuring point in the group of measuring points, and k is 0, 1 and 2 representing left, middle and right measuring points in the group of measuring points; and step 3c: applying the same parameter transformation method to another two measuring points in the group of measuring points, to obtain a relative positional relationship in the 3D space for fan-shaped arcs corresponding to images for the group of measuring points in the same analysis section.

3. The 3D terrain reconstruction method for a scoured area around bridge pier foundation based on a mechanical scanned imaging sonar according to claim 1, wherein the recognizing, according to a relative positional relationship of fan-shaped arcs corresponding to a group of measuring points and obtained by parameter transformation in a same analysis section, a terrain profile in the scoured area in step S4 specifically comprises: step 4a: determining whether a terrain type belongs to a depressed terrain, a raised terrain, a continuous terrain, or an uncertain terrain according to a relative positional relationship for a group of fan-shaped arcs, wherein in case of the uncertain terrain, the terrain type is further determined in combination with a terrain profile curve determined in an adjacent azimuth; step 4b: screening key points for forming a terrain profile curve according to following conditions: 1) the terrain type belongs to the depressed terrain, the raised terrain or the continuous terrain; 2) a fitted terrain profile curve is located between two fan-shaped arcs, and an intersected portion between adjacent fan-shaped surface is an area through which the terrain profile curve passes inevitably; and 3) upper and lower arcs restored by the measuring points each at least comprise a terrain profile point; and step 4c: determining key points of the terrain profile curve on each of the fan-shaped arcs while satisfying the above conditions, and connecting the key points through a smooth curve to obtain the terrain profile curve of a beam coverage area in the analysis section.

4. The 3D terrain reconstruction method for a scoured area around bridge pier foundation based on a mechanical scanned imaging sonar according to claim 3, wherein in step 4a, when distances from edge points on upper fan-shaped arcs of beams of auxiliary measuring points on two sides to the plane where the sonar is located are compared, an equation L.sub.02>l.sub.00>l.sub.02 is used if a greater distance is less than a distance from an edge point on a lower fan-shaped arc of an auxiliary measuring point on the other side; and concerning a connecting line a.sub.00a.sub.02 for mid-points on the upper fan-shaped arcs of the beams on the two sides, a mid-point a.sub.01 on an upper fan-shaped arc of a middle beam is located below the connecting line a.sub.00a.sub.02, indicating that the terrain type is the depressed terrain.

5. The 3D terrain reconstruction method for a scoured area around bridge pier foundation based on a mechanical scanned imaging sonar according to claim 1, wherein in step S6, the interpolation and fitting are performed on the massive spatially scattered data through the Kriging interpolation method, thus implementing 3D terrain reconstruction for the bridge pier foundation scoured area, and obtaining a foundation scouring contour map and a 3D visual map; and therefore, terrain information of the bridge pier foundation scoured area is intuitively viewed in the 3D space.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0026] The present disclosure will be further described in detail below with reference to the accompanying drawings and specific implementations:

[0027] FIG. 1 is a schematic flowchart of the present disclosure;

[0028] FIG. 2 is a schematic view for intercepting an analysis section of a sonar image of a scoured terrain around bridge pier foundation according to a method of the present disclosure;

[0029] FIG. 3 is a schematic view for transforming extracted parameters of a sonar image to a 3D space according to a method of the present disclosure; and

[0030] FIG. 4 is a schematic view for recognizing a depressed terrain profile according to a relative positional relationship of fan-shaped arcs for a group of measuring points according to a method of the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0031] As shown in figures, a 3D terrain reconstruction method for a scoured area around bridge pier foundation based on a mechanical scanned imaging sonar includes the following steps:

[0032] Step S1: Scan an overall terrain of a scoured area around bridge pier foundation with a mechanical scanning sonar from different azimuths to acquire n sonar images of the scoured terrain in the different azimuths.

[0033] Step S2: Intercept multiple analysis sections from each of acquired sonar images at a same distance.

[0034] Step S3: Extract key parameters of upper and lower edges on a terrain imaging strip in each of the analysis sections in the image, and transform extracted parameters to a 3D space according to an acoustic imaging principle, a fan-shaped beam surface of the sonar being represented with a fan-shaped arc.

[0035] Step S4: Recognize, according to a relative positional relationship of fan-shaped arcs for a group of measuring points in a same analysis section, a scouring terrain profile in the analysis section.

[0036] Step S5: Recognize terrain profiles in different azimuths and different analysis sections one by one, and respectively extract 3D coordinates of key points on each of profile curves to obtain required spatially scattered 3D coordinate data.

[0037] Step S6: Perform interpolation and fitting on the spatially scattered data with a Kriging interpolation method, thus implementing 3D reconstruction for the foundation scouring terrain.

[0038] The step of extracting key parameters of upper and lower edges on a terrain imaging strip in each of the analysis sections in the sonar image, and transforming extracted parameters to a 3D space according to an acoustic imaging principle in Step S3 specifically includes:

[0039] Step 3a: Label, in an analysis section p.sub.j,m.sup.i at a distance d.sub.j.sup.i away from a sonar center, upper and lower edge points M, N on an imaging strip, ranges r.sub.jk.sup.i, R.sub.jk.sup.i from an original sonar center point P.sub.jk.sup.i to the upper and lower edge points M, N, and a large opening angle θ of the fan-shaped beam surface, as shown in a of FIG. 3(b).

[0040] Step 3b: Draw a circular arc with the upper edge point M as a center and the r.sub.jk.sup.i as a radius, an intersection with the analysis section p.sub.j,m.sup.i being an equivalent sonar center point Q.sub.jk.sup.i, and the equivalent sonar center point Q.sub.jk.sup.i being T.sub.jk.sup.i higher than a plane where the original sonar center is located; draw a fan-shaped arc with the equivalent sonar center point Q.sub.jk.sup.i that is T.sub.jk.sup.i higher than the plane where the sonar is located as a center, the r.sub.jk.sup.i as a radius, and the large opening angle θ of the fan-shaped beam surface to obtain a corresponding fan-shaped arc of the upper edge point M in the 3D space; draw a circular arc with the lower edge point N as a center and the R.sub.jk.sup.i as a radius, an intersection with the analysis section p.sub.j,m.sup.i being an equivalent sonar center point K.sub.jk.sup.i that is t.sub.jk.sup.i higher than the plane where the original sonar center is located; and draw a fan-shaped arc with the as a center, the as a radius and the angle θ to obtain a corresponding fan-shaped arc of the lower edge point N in the 3D space, where the superscript i represents a pier No., the first subscript j represents an azimuth in which a group of measuring points are located relative to a bridge pier, the second subscript k represents a specific measuring point in the group of measuring points, and k is 0, 1 and 2 representing left, middle and right measuring points in the group of measuring points.

[0041] Step 3b is described with b and c in FIG. 3 as an example: Ranges from the original sonar center point to the upper and lower edge points M, N of the analysis section are and . A circular arc is drawn with the upper edge point M as a center and the as a radius, an intersection with the analysis section being an equivalent sonar center point , and the equivalent sonar center point being higher than a plane where the original sonar center is located. A fan-shaped arc is drawn with the equivalent sonar center point that is higher than the plane where the sonar is located as a center, the as a radius, and the large opening angle θ of the fan-shaped beam surface to obtain a corresponding fan-shaped arc of the upper edge point M in FIG. 3 (c) in the 3D space. Likewise, a circular arc is drawn with the lower edge point N as a center and the as a radius, an intersection with the analysis section being an equivalent sonar center point that is higher than the plane where the original sonar center is located. A fan-shaped arc with the as a center, the as a radius and the angle θ to obtain a corresponding fan-shaped arc of the lower edge point N in the 3D space. The superscript i represents a pile No., the first subscript j represents an azimuth in which a group of measuring points are located relative to a bridge pier, the second subscript k represents a specific measuring point in the group of measuring points, and k is 0, 1 and 2 representing left, middle and right measuring points in the group of measuring points.

[0042] Step 3c: Apply the same parameter transformation method to another two measuring points in the group of measuring points, to obtain a relative positional relationship in the 3D space for fan-shaped arcs corresponding to images for the group of measuring points in the same analysis section.

[0043] The step of recognizing, according to a relative positional relationship of fan-shaped arcs corresponding to a group of measuring points and obtained by parameter transformation in a same analysis section, a terrain profile in the scoured area in Step S4 specifically includes:

[0044] Step 4a: Determine whether a terrain type belongs to a depressed terrain, a raised terrain, a continuous terrain, or an uncertain terrain according to a relative positional relationship for a group of fan-shaped arcs, where in case of the uncertain terrain, the terrain type is further determined in combination with a terrain profile curve determined in an adjacent azimuth.

[0045] Step 4b: Screen key points for forming a terrain profile curve according to following conditions: 1) the terrain type belongs to the depressed terrain, the raised terrain or the continuous terrain; 2) a fitted terrain profile curve is located between two fan-shaped arcs, and an intersected portion between adjacent fan-shaped surfaces is an area through which the terrain profile curve passes inevitably; and 3) upper and lower arcs restored by the measuring points each at least include a terrain profile point.

[0046] Step 4c: Determine key points of a terrain profile curve on each of the fan-shaped arcs while satisfying the above conditions, and connecting the key points through a smooth curve to obtain the terrain profile curve of a beam coverage area in the analysis section.

[0047] As shown in FIG. 3, in Step 4a, when distances from edge points on upper fan-shaped arcs of beams of auxiliary measuring points on two sides to the plane where the sonar is located are compared, an equation L.sub.02>l.sub.00>l.sub.02 is used if a greater distance is less than a distance from an edge point of a lower fan-shaped arc of an auxiliary measuring point on the other side; and concerning a connecting line a.sub.00a.sub.02 for mid-points on the upper fan-shaped arcs of the beams on the two sides, a mid-point a.sub.01 on an upper fan-shaped arc of a middle beam is located below the connecting line a.sub.00a.sub.02, indicating that the terrain type is the depressed terrain.

[0048] As shown in FIG. 4, when the terrain type is determined as the depressed terrain, the fitted terrain profile curve is located between two fan-shaped arcs, the intersected portions between adjacent fan-shaped arc surfaces (namely areas A.sub.01C.sub.02D.sub.00c.sub.00 and c.sub.02C.sub.02D.sub.01B.sub.01 in FIG. 4) are areas through which the terrain profile curve passes inevitably, and upper and lower fan-shaped arcs restored by the measuring points each at least include a terrain profile point.

[0049] In Step S6, the interpolation and fitting are performed on the massive spatially scattered data through the Kriging interpolation method, thus implementing 3D terrain reconstruction for the scoured area around bridge pier foundation, and obtaining a scour contour map and a 3D visual map of the scour; and therefore, terrain information of the scoured area around bridge pier foundation can be intuitively viewed in the 3D space.

[0050] The above are merely preferred specific implementations of the present disclosure. These specific implementations are different implementations based on overall concepts, rather than limitations to the protection scope of the present disclosure. The implementations are illustrated only for explaining functions and use methods of the present disclosure rather than limitations to the present disclosure. Various changes or replacements made by a person of ordinary skill in the art without departing from the substantial scope of the present disclosure fall into the protection scope of the present disclosure.