Boundary extraction method, device, equipment, medium and product for geological defect area

Abstract

The present invention relates to the field of tunnel engineering technology, and discloses a method, device, equipment, medium, and product for extracting the boundary of a geological defect area. The method comprises: obtaining point cloud data of a tunnel where the geological defect area is located, selecting target point cloud data with an overexcavation value greater than a first preset threshold value from the point cloud data of the tunnel where the geological defect area is located; selecting multiple boundary starting points based on the dimensional direction of the target point cloud data; searching for multiple defect points; selecting a target defect point from the multiple defect points; determining whether the target defect point exists in the multiple boundary starting points, and if so, screening the multiple boundary starting points to obtain multiple target boundary points; and sequentially connecting the multiple target boundary points to obtain a target geological defect boundary. The present invention improves the accuracy of boundary extraction of the geological defect area by extracting and judging points in the point cloud data of the tunnel where the geological defect area is located.

Description

Technical Field

[0001] The present invention relates to the technical field of tunnel engineering, and in particular to a method, device, equipment, medium and product for extracting the boundary of a geological defect area. Background Art

[0002] When drilling and blasting tunnels are used, encountering unfavorable geological conditions can lead to overexcavation due to geological defects. This occurs when geological defects (such as faults, fracture zones, weak interlayers, and other unfavorable geological conditions) prevent the excavation from proceeding according to the pre-designed excavation method and dimensions, resulting in the excavation profile significantly exceeding the pre-designed excavation profile. When overexcavation due to geological defects occurs, the geological defect area must be measured and its boundaries accurately demarcated for resolution and cost settlement.

[0003] In related technologies, the method for identifying geological defect areas involves using a surveying instrument to obtain multiple measurement points within the geological defect area, estimating the geological defect boundary based on these multiple measurement points, and then delineating the geological defect area. This method suffers from problems such as sparse and random measurement points, resulting in inaccurate geological defect boundaries and the delineated geological defect area being inconsistent with actual conditions. Summary of the Invention

[0004] In view of this, the present invention provides a method, device, equipment, medium and product for extracting the boundary of a geological defect area to solve the problem of inaccurate geological defect boundaries obtained by the boundary extraction method of the geological defect area in the related art.

[0005] In a first aspect, the present invention provides a method for extracting the boundary of a geological defect area, comprising: acquiring point cloud data of a tunnel where the geological defect area is located, and selecting target point cloud data whose over-excavation value is greater than a first preset threshold value in the point cloud data of the tunnel where the geological defect area is located; the over-excavation value is a value at which the excavation contour exceeds the preset excavation contour during tunnel excavation; selecting multiple boundary starting points according to the dimensional direction of the target point cloud data; the multiple boundary starting points are multiple points that are within a first preset range in the first dimensional direction of the target point cloud data, within a second preset range in the second dimensional direction, and have the largest corresponding numerical value in the second dimensional direction; searching within a third preset range with the first target boundary starting point having the smallest corresponding numerical value in the first dimensional direction among the multiple boundary starting points as the center. A plurality of defect points having over-excavation values ​​greater than a first preset threshold value; a target defect point is selected from the plurality of defect points; the target defect point is a defect point with the smallest azimuth in the counterclockwise direction, centered on the first target boundary starting point and starting from the first target boundary starting point to the second target boundary starting point adjacent to the first target boundary starting point; it is determined whether the target defect point exists in the plurality of boundary starting points, and if the target defect point exists in the plurality of boundary starting points, the plurality of boundary starting points are screened to obtain a plurality of target boundary points; the plurality of target boundary points are a plurality of points in the plurality of boundary starting points whose corresponding values ​​in the first dimensional direction are less than or equal to the corresponding values ​​of the target defect point in the first dimensional direction; the plurality of target boundary points are sequentially connected to obtain a target geological defect boundary.

[0006] The present invention selects target point cloud data with overexcavation values ​​greater than a first preset threshold from the point cloud data of a tunnel located in a geological defect area. This screens out target point cloud data with significant overexcavation characteristics, eliminates interference and erroneous data, and provides a more reliable data foundation for more accurately determining boundary starting points. The present invention selects multiple boundary starting points based on dimensional directions and preset ranges, enabling precise determination of multiple boundary starting points from complex target point cloud data, ensuring that the selection of multiple boundary starting points conforms to the characteristics of the geological defect area. The present invention searches for multiple defect points with overexcavation values ​​greater than the first preset threshold within a third preset range, centering on the first target boundary starting point with the smallest corresponding value in the first dimensional direction among the multiple boundary starting points. This method gradually explores the geological defect area in an orderly manner, ensuring that no defect points are missed. Furthermore, the setting of the third preset range further limits the search area, ensuring comprehensiveness while improving search efficiency and avoiding blind searches within an overly large range, making defect point discovery more accurate and efficient. The present invention selects target defect points by using clear conditions such as azimuth angles, enabling identification of the most representative and critical target defect points among multiple defect points. This selection method ensures the uniqueness and certainty of the target defect points. The present invention determines whether the target defect point exists in multiple boundary starting points. If the target defect point exists in multiple boundary starting points, the multiple boundary starting points are screened to obtain multiple target boundary points. By re-screening the boundary starting points, it is ensured that the points ultimately used to connect to form the boundary are points that have a reasonable relationship with the target defect point in the dimensional direction, thereby improving the accuracy and rationality of the geological defect boundary. The present invention selects multiple target boundary points through a series of rigorous steps, and the target geological defect boundary formed after sequential connection is complete and accurate, clearly defining the scope of the geological defect area, and providing intuitive and reliable geological information for engineering personnel. Compared with related technologies, the target geological defect boundary obtained by the present invention is more in line with the actual situation, and improves the accuracy of geological defect boundary identification.

[0007] In an optional embodiment, multiple boundary starting points are selected according to the dimensional direction of the target point cloud data, including: obtaining the average point spacing of the point cloud data of the tunnel where the geological defect area is located; taking the difference between the coordinates corresponding to the target point cloud data in the first dimensional direction and the average point spacing of the first preset multiple as the first boundary, and taking the sum of the coordinates corresponding to the target point cloud data in the first dimensional direction and the average point spacing of the first preset multiple as the second boundary, to obtain a first preset range; in the point cloud data of the tunnel where the geological defect area is located, selecting a first point set within the first preset range; taking the sum of the coordinates corresponding to the target point cloud data in the second dimensional direction and the average point spacing of the second preset multiple as the third boundary, and taking the sum of the coordinates corresponding to the target point cloud data in the second dimensional direction and the average point spacing of the third preset multiple as the fourth boundary, to obtain a second preset range; the second preset multiple is less than the third preset multiple; in the first point set, selecting multiple points within the second preset range, with over-excavation values ​​greater than the first preset threshold, and with the largest corresponding values ​​in the second dimensional direction, as multiple boundary starting points.

[0008] In an optional embodiment, with the first target boundary starting point having the smallest corresponding value in the first dimensional direction among multiple boundary starting points as the center, multiple defect points with over-excavation values ​​greater than the first preset threshold are searched within a third preset range, including: with the first target boundary starting point having the smallest corresponding value in the first dimensional direction among multiple boundary starting points as the center, searching with an average point spacing of a fourth preset multiple as a radius to determine whether there are multiple points with over-excavation values ​​greater than the first preset threshold; if there are multiple points with over-excavation values ​​greater than the first preset threshold, multiple defect points are obtained.

[0009] In some optional implementations, if there are not multiple points where the over-digging value is greater than the first preset threshold, the fourth preset multiple is adjusted until there are multiple points where the over-digging value is greater than the first preset threshold.

[0010] In an optional embodiment, multiple boundary starting points are screened to obtain multiple target boundary points, including: eliminating boundary starting points whose corresponding values ​​in the first dimension direction are greater than the corresponding values ​​of the target defect point in the first dimension direction, and forming multiple target boundary points from the remaining multiple boundary starting points.

[0011] In an optional embodiment, the boundary extraction method of the geological defect area also includes: if the target defect point does not exist in multiple boundary starting points, taking the target defect point as the first target boundary starting point, returning to the first target boundary starting point with the smallest corresponding value in the first dimension direction among the multiple boundary starting points as the center, and searching for multiple defect points with over-excavation values ​​greater than the first preset threshold within a third preset range until the target defect point exists in multiple boundary starting points.

[0012] According to the present invention, if the target defect point does not exist in multiple boundary starting points, the target defect point is used as the first target boundary starting point, and the process returns to the first target boundary starting point with the smallest corresponding value in the first dimensional direction among the multiple boundary starting points, and searches for multiple defect points with over-excavation values ​​greater than the first preset threshold within a third preset range until the target defect point exists in multiple boundary starting points. This ensures that no possible defect point is missed, and comprehensively covers the entire area where geological defects may exist through a cycle, thereby improving the integrity and accuracy of target geological defect extraction.

[0013] In the second aspect, the present invention provides a boundary extraction device for a geological defect area, including: a target point cloud data determination module, used to obtain the point cloud data of the tunnel where the geological defect area is located, and select target point cloud data with an over-excavation value greater than a first preset threshold value in the point cloud data of the tunnel where the geological defect area is located; the over-excavation value is the value of the excavation contour exceeding the preset excavation contour during tunnel excavation; a boundary starting point selection module, used to select multiple boundary starting points according to the dimensional direction of the target point cloud data; the multiple boundary starting points are within a first preset range in the first dimensional direction of the target point cloud data, within a second preset range in the second dimensional direction, and have the largest corresponding numerical value in the second dimensional direction; a defect point search module, used to search within a third preset range with the first target boundary starting point with the smallest corresponding numerical value in the first dimensional direction among the multiple boundary starting points as the center. Find multiple defect points whose over-excavation values ​​are greater than a first preset threshold; a target defect point selection module is used to select a target defect point from multiple defect points; the target defect point is a defect point with the smallest azimuth in the counterclockwise direction, centered on the first target boundary starting point and starting from the first target boundary starting point to the second target boundary starting point adjacent to the first target boundary starting point; a judgment module is used to judge whether the target defect point exists in multiple boundary starting points. If the target defect point exists in multiple boundary starting points, the multiple boundary starting points are screened to obtain multiple target boundary points; the multiple target boundary points are multiple points in the multiple boundary starting points whose corresponding values ​​in the first dimensional direction are less than or equal to the corresponding values ​​of the target defect point in the first dimensional direction; a target geological defect boundary module is used to connect the multiple target boundary points in sequence to obtain a target geological defect boundary.

[0014] In a third aspect, the present invention provides a computer device comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the computer instructions to thereby execute the boundary extraction method for the geological defect area of ​​the first aspect or any corresponding embodiment thereof.

[0015] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the method for extracting the boundary of a geological defect area according to the first aspect or any corresponding embodiment thereof.

[0016] In a fifth aspect, the present invention provides a computer program product comprising computer instructions for causing a computer to execute the method for extracting the boundary of a geological defect area according to the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in related technologies, the following briefly introduces the drawings required for use in the specific embodiments or related technical descriptions. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0018] Figure 1 4 is a flow chart of a method for extracting the boundary of a geological defect area according to an embodiment of the present invention.

[0019] Figure 2 It is a schematic diagram of over-excavation and under-excavation values ​​of a tunnel project according to an embodiment of the present invention.

[0020] Figure 3 2 is a schematic diagram of a process for selecting multiple boundary starting points according to an embodiment of the present invention.

[0021] Figure 4 2 is a schematic diagram of a target defect point selection process according to an embodiment of the present invention.

[0022] Figure 5 2 is a schematic diagram of a process for determining a target geological defect boundary according to an embodiment of the present invention.

[0023] Figure 6 4 is a flow chart of another method for extracting the boundary of a geological defect area according to an embodiment of the present invention.

[0024] Figure 7 Schematic diagram of the target geological defect boundary according to an embodiment of the present invention.

[0025] Figure 8 4 is a flow chart of another method for extracting the boundary of a geological defect area according to an embodiment of the present invention.

[0026] Figure 9 4 is a structural block diagram of a device for extracting a boundary of a geological defect area according to an embodiment of the present invention.

[0027] Figure 10Schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0028] To make the purpose, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without making creative efforts shall fall within the scope of protection of the present invention.

[0029] During tunnel construction, when drilling and blasting is used, unfavorable geological conditions can lead to overexcavation. This occurs when the excavation process is unable to follow the designed excavation method and dimensions due to geological defects (such as faults, fracture zones, and weak interlayers). To prevent overexcavation, it is necessary to measure the geological defect area and accurately demarcate its boundaries for management.

[0030] Related technologies primarily use total stations or other surveying instruments to acquire multiple measurement points within a geological defect area, plot the cross-sectional area of ​​the defect overexcavation at preset intervals, and manually estimate the boundaries of the geological defect. These technologies suffer from sparse and random measurement points, resulting in inaccurate overexcavation boundaries and a cumbersome and inefficient calculation process.

[0031] An embodiment of the present invention provides a method for extracting the boundary of a geological defect area, which improves the accuracy of the boundary extraction of the geological defect area by extracting and judging points in the point cloud data of a tunnel where the geological defect area is located.

[0032] According to an embodiment of the present invention, an embodiment of a method for extracting the boundary of a geological defect area is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0033] In this embodiment, a method for extracting the boundary of a geological defect area is provided, which can be used in computer equipment. Figure 1 FIG. 1 is a flow chart of a method for extracting the boundary of a geological defect area according to an embodiment of the present invention. Figure 1 As shown, the process includes the following steps:

[0034] Step S101, obtaining point cloud data of the tunnel where the geological defect area is located, and selecting target point cloud data with an over-excavation value greater than a first preset threshold value from the point cloud data of the tunnel where the geological defect area is located; the over-excavation value is the value of the excavation contour exceeding the preset excavation contour during the tunnel excavation process.

[0035] Among them, the overbreak value is the value of the actual excavation contour exceeding the preset excavation contour corresponding to the design section during the tunnel excavation process; the underbreak value is the value of the actual excavation contour not reaching the preset excavation contour corresponding to the preset section during the tunnel excavation process; Figure 2 As shown in FIG, it is a schematic diagram of the over-excavation and under-excavation values ​​of a tunnel project. The difference between the measured point cloud and the designed section is the over-excavation and under-excavation value. In the embodiment of the present invention, in the case of under-excavation, the measured point cloud can be made to conform to the designed section by excavating again. However, in the case of over-excavation, the stability of the tunnel may be destroyed and the construction cost may be increased. Therefore, the embodiment of the present invention mainly uses the over-excavation value to study the geological defect area.

[0036] In some optional embodiments, the process of obtaining the over-excavation value is: obtaining point cloud data of the tunnel where the geological defect area is located, projecting the point cloud data onto a preset section, obtaining a projection result, calculating the distance between the point cloud data and the projection result, and obtaining the over-excavation value.

[0037] The first preset value is a preset geological defect threshold. For example, the first preset value may be 15 centimeters.

[0038] In step S102, a plurality of boundary starting points are selected according to the dimensional direction of the target point cloud data; the plurality of boundary starting points are points within a first preset range in the first dimensional direction of the target point cloud data, within a second preset range in the second dimensional direction, and having an over-excavation value greater than a first preset threshold value, and having the largest corresponding value in the second dimensional direction.

[0039] The dimensional direction includes a first dimensional direction and a second dimensional direction perpendicular to the first dimensional direction. The first dimensional direction may be a horizontal direction of a plane rectangular coordinate system, and the second dimensional direction may be a vertical direction of the plane rectangular coordinate system. The corresponding value in the second dimensional direction is the coordinate value in the second dimensional direction. For example, when the second dimensional direction is the y-axis direction, the corresponding value in the second dimensional direction is the y value.

[0040] In some optional embodiments, the first preset range is a length range between two preset coordinate points in the first dimension. Specifically, the first preset range is a range with a first boundary being the difference between the coordinate corresponding to the target point cloud data in the first dimension and the average point spacing of the first preset multiple as a first boundary, and a second boundary being the sum of the coordinate corresponding to the target point cloud data in the first dimension and the average point spacing of the first preset multiple as a second boundary. The average point spacing is the average value of the point spacing of the point cloud data of the tunnel where the geological defect area is located.

[0041] In some optional embodiments, the second preset range is a length range between two preset coordinate points in the second dimension. Specifically, the second preset range is a range with the sum of the coordinates corresponding to the target point cloud data in the second dimension and the average point spacing of the second preset multiple as the third boundary, and with the sum of the coordinates corresponding to the target point cloud data in the second dimension and the average point spacing of the third preset multiple as the fourth boundary. The second preset multiple is greater than the third preset multiple.

[0042] Exemplarily, when the coordinates of the target point cloud data are (xs, ys), the first dimension direction is the x-axis direction, the second dimension direction is the y-axis direction, the average point spacing is D, the first preset multiple is n, the second preset multiple is i, and the third preset multiple is i+1, the first preset range is [xs-nD, xs+nD], and the second preset range is (ys+iD, ys+(i+1)D), where n is any one of 1 to 5, the initial value of i is 0, and it can be increased sequentially.

[0043] In some optional embodiments, multiple boundary starting points are selected according to the dimensional direction of the target point cloud data, including: obtaining the average point spacing of the point cloud data of the tunnel where the geological defect area is located; taking the difference between the coordinates corresponding to the target point cloud data in the first dimensional direction and the average point spacing of the first preset multiple as the first boundary, and taking the sum of the coordinates corresponding to the target point cloud data in the first dimensional direction and the average point spacing of the first preset multiple as the second boundary, to obtain a first preset range; in the point cloud data of the tunnel where the geological defect area is located, selecting a first point set within the first preset range; taking the sum of the coordinates corresponding to the target point cloud data in the second dimensional direction and the average point spacing of the second preset multiple as the third boundary, and taking the sum of the coordinates corresponding to the target point cloud data in the second dimensional direction and the average point spacing of the third preset multiple as the fourth boundary, to obtain a second preset range; the second preset multiple is less than the third preset multiple; in the first point set, selecting multiple points within the second preset range, with over-excavation values ​​greater than the first preset threshold, and with the largest corresponding values ​​in the second dimensional direction, as multiple boundary starting points.

[0044] The average point spacing is determined using a sampling method or a point density calculation method. Specifically, the sampling method involves randomly extracting multiple point samples from the point cloud data of the tunnel where the geological defect area is located, calculating the distances between the multiple point samples, and performing statistical analysis on the distances between the multiple point samples to obtain the average point spacing. The point density calculation method involves statistically analyzing the number of point clouds and the area within the point cloud data of the tunnel where the geological defect area is located, determining the point density based on the area and number of point clouds, and estimating the average point spacing using the point density.

[0045] For example, Figure 3 As shown in the figure, it is a schematic diagram of the process of selecting multiple boundary starting points. First, any target point cloud data S with an over-excavation value greater than 15 cm is selected from the point cloud data of the tunnel where the geological defect area is located. The coordinates of S are (xs, ys). For example, the point cloud data of the tunnel where the geological defect area is located is Figure 3 As shown in the figure, in the point cloud data of the tunnel where the geological defect area is located, the first point set in [xs-nD, xs+nD] is selected, that is, the value of the first point set x is: xs-nD≤x≤xs+nD; in the first point set, the over-excavation value in (ys+iD, ys+(i+1)D) is greater than 15 cm, and the two largest points in the y-axis direction are selected as the boundary starting points, that is, Figure 3 The boundary starting points Q1 and Q0 in the , where the value corresponding to Q0 in the x-axis direction is greater than the value corresponding to Q1 in the x-axis direction.

[0046] In an embodiment of the present invention, there are multiple boundary starting points, and the multiple boundary starting points are sorted according to the size of the corresponding values ​​in the x-axis direction. The larger the corresponding value in the x-axis direction, the higher the sorting.

[0047] In some optional embodiments, when the first point concentration is within the second preset range but there is no point with an over-excavation value less than the first preset threshold, the second preset multiple and the third preset multiple are adjusted. Specifically, the second preset multiple and the third preset multiple are increased until the first point concentration is within the second preset range and there is a point with an over-excavation value less than the first preset threshold.

[0048] Step S103 , taking the first target boundary starting point with the smallest corresponding value in the first dimensional direction among the multiple boundary starting points as the center, searching for multiple defect points with over-excavation values ​​greater than a first preset threshold within a third preset range.

[0049] The corresponding value in the first dimension is the corresponding coordinate value in the first dimension. The third preset range is the range within a circle having a center and a radius. Specifically, the third preset range is the range with the first target boundary point as the center and the average point spacing of the fourth preset multiple as the radius.

[0050] In some optional embodiments, with the first target boundary starting point having the smallest corresponding value in the first dimensional direction among multiple boundary starting points as the center, multiple defect points with over-excavation values ​​greater than the first preset threshold are searched within a third preset range, including: with the first target boundary starting point having the smallest corresponding value in the first dimensional direction among multiple boundary starting points as the center, searching with an average point spacing of a fourth preset multiple as a radius to determine whether there are multiple points with over-excavation values ​​greater than the first preset threshold; if there are multiple points with over-excavation values ​​greater than the first preset threshold, multiple defect points are obtained.

[0051] In some optional implementations, if there are not multiple points where the over-digging value is greater than the first preset threshold, the fourth preset multiple is adjusted until there are multiple points where the over-digging value is greater than the first preset threshold.

[0052] Exemplarily, among the boundary starting points Q1 and Q0, the value corresponding to Q1 in the x-axis direction is the smallest. Then, with Q1 as the center and αD as the radius, search to see whether there are multiple points with over-excavation values ​​greater than the first preset threshold, wherein the value of α is any one in the range of 1 to 50. When there are no multiple points with over-excavation values ​​greater than the first preset threshold, increase the value of α until there are multiple points with over-excavation values ​​greater than the first preset threshold.

[0053] Step S104, selecting a target defect point from among multiple defect points; the target defect point is the defect point with the smallest azimuth angle in the counterclockwise direction, with the first target boundary starting point as the center and the orientation from the first target boundary starting point to the second target boundary starting point adjacent to the first target boundary starting point as the starting orientation.

[0054] Among them, since multiple boundary starting points are sorted according to the size of the corresponding values ​​in the x-axis direction, the larger the corresponding values ​​in the x-axis direction, the higher the sorting, so the first target boundary starting point is the boundary starting point with the last sorting, and the second target boundary starting point is the boundary starting point with the second to last sorting.

[0055] For example, Figure 4 The figure shows a schematic diagram of the target defect point selection process. The point set m is a set of multiple defect points. The starting point Qn of the first target boundary is used as the center point, and the direction from Qn to the starting point Qn-1 of the second target boundary is used as the starting direction. The azimuth angles of multiple defect points are calculated in the counterclockwise direction, and the defect point with the smallest azimuth angle is determined as the target defect point P.

[0056] Step S105, determine whether the target defect point exists in multiple boundary starting points. If the target defect point exists in multiple boundary starting points, screen the multiple boundary starting points to obtain multiple target boundary points; the multiple target boundary points are multiple points in the multiple boundary starting points whose corresponding values ​​in the first dimension direction are less than or equal to the corresponding values ​​of the target defect point in the first dimension direction.

[0057] In some optional embodiments, multiple boundary starting points are screened to obtain multiple target boundary points, including: eliminating boundary starting points whose corresponding values ​​in the first dimension direction are greater than the corresponding values ​​of the target defect point in the first dimension direction, and forming multiple target boundary points from the remaining multiple boundary starting points.

[0058] Exemplarily, if the target defect point exists in multiple boundary starting points, and the corresponding symbol in the multiple boundary starting points is Qi, the point set composed of multiple boundary starting points is (Q0, Q1, ..., Qi, ..., Qn-1, Qn), then the boundary starting point before Qi is eliminated to obtain multiple target boundary points, then the point set composed of multiple target boundary points is (Q1, ..., Qn-1, Qn), where Q1 is Qi in the point set composed of multiple boundary starting points.

[0059] In some optional embodiments, the boundary extraction method of the geological defect area also includes: if the target defect point does not exist in multiple boundary starting points, taking the target defect point as the first target boundary starting point, returning to the first target boundary starting point with the smallest corresponding value in the first dimension direction among the multiple boundary starting points as the center, and searching for multiple defect points with over-excavation values ​​greater than the first preset threshold within a third preset range until the target defect point exists in multiple boundary starting points.

[0060] In an embodiment of the present invention, when the target defect point does not exist in multiple boundary starting points, the target defect point is used as the first target boundary starting point, and the process returns to the first target boundary starting point with the smallest corresponding value in the first dimensional direction among the multiple boundary starting points, and searches for multiple defect points with over-excavation values ​​greater than the first preset threshold within a third preset range until the target defect point exists in multiple boundary starting points. This ensures that no possible defect point is missed, and by cyclically covering the entire area where geological defects may exist, the completeness and accuracy of target geological defect extraction are improved.

[0061] Step S106: Connect multiple target boundary points in sequence to obtain a target geological defect boundary.

[0062] The target geological defect boundary is a boundary of a geological defect area with high accuracy obtained by selecting and connecting target boundary points.

[0063] For example, Figure 5As shown, it is a schematic diagram of the process of determining the target geological defect boundary. Multiple target boundary points (the multiple target boundary points may include Q0, Q1, Q2 and Qn, etc.) are connected in sequence to obtain the target geological defect boundary.

[0064] The present invention selects target point cloud data with overexcavation values ​​greater than a first preset threshold from the point cloud data of a tunnel located in a geological defect area. This screens out target point cloud data with significant overexcavation characteristics, eliminates interference and erroneous data, and provides a more reliable data foundation for more accurately determining boundary starting points. The present invention selects multiple boundary starting points based on dimensional directions and preset ranges, enabling precise determination of multiple boundary starting points from complex target point cloud data and ensuring that the selection of multiple boundary starting points conforms to the characteristics of the geological defect area. The present invention searches for multiple defect points with overexcavation values ​​greater than the first preset threshold within a third preset range, centering on the first target boundary starting point with the smallest corresponding value in the first dimensional direction among the multiple boundary starting points. This method explores the geological defect area in an orderly manner, ensuring that no defect points are missed. Furthermore, the third preset range further limits the search area, ensuring comprehensiveness while improving search efficiency and avoiding blind searches within an overly large range, making defect point discovery more accurate and efficient. The present invention selects target defect points based on clear conditions such as azimuth, enabling identification of the most representative and critical target defect points among multiple defect points. This selection method ensures the uniqueness and certainty of the target defect points. The embodiment of the present invention determines whether the target defect point exists in multiple boundary starting points. If the target defect point exists in multiple boundary starting points, the multiple boundary starting points are screened to obtain multiple target boundary points. By re-screening the boundary starting points, it is ensured that the points ultimately used to connect to form the boundary are points that have a reasonable relationship with the target defect point in the dimensional direction, thereby improving the accuracy and rationality of the geological defect boundary. The embodiment of the present invention selects multiple target boundary points through a series of rigorous steps. The target geological defect boundary formed after sequential connection is complete and accurate, clearly defines the scope of the geological defect area, and provides intuitive and reliable geological information for engineering personnel. Compared with related technologies, the target geological defect boundary obtained by the present invention is more in line with the actual situation and improves the accuracy of geological defect boundary identification.

[0065] In this embodiment, a method for extracting the boundary of a geological defect area is provided, which can be used in computer equipment. Figure 6 FIG. 1 is a flow chart of another method for extracting the boundary of a geological defect area according to an embodiment of the present invention. Figure 6 As shown, the process includes the following steps:

[0066] Step S601: Acquire point cloud data of the tunnel where the geological defect area is located, project the point cloud data onto a preset section, and obtain a projection result.

[0067] Among them, the point cloud data of the tunnel where the geological defect area is located is the three-dimensional point cloud data of the tunnel obtained using three-dimensional laser scanning technology.

[0068] In some optional implementations, after acquiring the point cloud data of the tunnel, the point cloud data of the tunnel is converted into a unified engineering coordinate system through target coordinates, and noise reduction processing is performed on the point cloud data.

[0069] In some optional implementations, projecting the point cloud data onto a preset cross-section to obtain a projection result includes: projecting each point in the point cloud data onto the preset cross-section through vector operations to obtain a projection result.

[0070] Step S602: Calculate the distance between the point cloud data and the projection result to obtain the over-excavation value.

[0071] The Euclidean distance between the coordinates of the point cloud data and the coordinates corresponding to the projection result is calculated to obtain the over-excavation value.

[0072] Step S603: performing plane projection on the point cloud data of the tunnel to obtain a target point cloud model.

[0073] In some optional embodiments, the plane can be selected as a specific plane perpendicular to the tunnel axis, or a plane with a certain angle relationship with the tunnel direction, etc., and each point in the point cloud data is projected onto the selected plane using a vector projection method.

[0074] Step S604 : performing color rendering on the target point cloud model according to the over-excavation value to obtain a color rendering result.

[0075] In some optional embodiments, during the coloring rendering of the target point cloud model of the tunnel, the depth of the rendered color can be adjusted according to the size of the over-excavation value. For example, the larger the over-excavation value, the darker the color, and the smaller the over-excavation value, the lighter the color.

[0076] Step S605 : selecting different color areas according to the coloring rendering result to obtain the target geological defect boundary.

[0077] In the embodiment of the present invention, a boundary detection algorithm is used to select different color regions to obtain the target geological defect boundary.

[0078] In the embodiment of the present invention, Figure 7 The figure shows the target geological defect boundary. The point cloud data of the tunnel is projected on a plane to obtain the target point cloud model. Different color areas are selected to obtain the target geological defect boundary. Figure 7 After processing in the embodiment of the present invention, the three-dimensional point cloud data of the tunnel is transformed into a target point cloud model with the target geological defect boundary. Figure 7 In the figure, the irregular solid and dotted lines in the target point cloud model are the target geological defect boundaries.

[0079] In this embodiment, a method for extracting the boundary of a geological defect area is provided, which can be used in computer equipment. Figure 8 FIG. 1 is a flow chart of another method for extracting the boundary of a geological defect area according to an embodiment of the present invention. Figure 8 As shown, the process includes the following steps:

[0080] The first method for extracting the target geological defect boundary is: collecting the tunnel 3D point cloud, converting the coordinates of the tunnel 3D point cloud into the engineering coordinate system, calculating the over-excavation and under-excavation values ​​of the tunnel 3D point cloud, mapping the tunnel 3D point cloud to a plane, performing visual rendering based on the over-excavation and under-excavation values, and manually or automatically drawing the geological defect boundary based on the rendering results, where the geological defect boundary here is the target geological defect boundary.

[0081] The second method for extracting the target geological defect boundary is to automatically calculate the geological defect boundary. The specific process is: based on the geological defect boundary, the plane average point spacing is calculated, the boundary starting point is extracted according to the average point spacing, the geological defect boundary points are cyclically extracted, and the validity of the geological defect area is judged to complete the geological defect boundary extraction, where the geological defect boundary here is the target geological defect boundary.

[0082] In an embodiment of the present invention, the three-dimensional point cloud of the tunnel is the point cloud data of the tunnel where the geological defect area is located. The boundary starting point is extracted according to the average point spacing, including: selecting target point cloud data with an over-excavation value greater than a first preset threshold in the point cloud data of the tunnel where the geological defect area is located, and selecting multiple boundary starting points according to the dimensional direction of the target point cloud data; the multiple boundary starting points are multiple points that are within a first preset range in the first dimensional direction of the target point cloud data, within a second preset range in the second dimensional direction, and have the largest corresponding value in the second dimensional direction. For details, please see Figure 1 Step S102 of the illustrated embodiment will not be described in detail here.

[0083] In some optional embodiments, the cyclic extraction of geological defect boundary points includes: taking the first target boundary starting point with the smallest corresponding value in the first dimension direction among the multiple boundary starting points as the center, searching for multiple defect points with overexcavation values ​​greater than the first preset threshold within a third preset range, and selecting a target defect point from the multiple defect points; the target defect point is the defect point with the smallest azimuth angle in the counterclockwise direction from the first target boundary starting point to the second target boundary starting point adjacent to the first target boundary starting point as the starting direction among the multiple defect points. For details, please refer to Figure 1 Step S103 and step S104 of the illustrated embodiment will not be described in detail here.

[0084] In some optional implementations, the validity judgment of the geological defect area includes: judging whether the target defect point exists in multiple boundary starting points; if the target defect point exists in multiple boundary starting points, screening the multiple boundary starting points to obtain multiple target boundary points; the multiple target boundary points are multiple points in the multiple boundary starting points whose corresponding values ​​in the first dimension direction are less than or equal to the corresponding values ​​of the target defect point in the first dimension direction. For details, please see Figure 1 Step S105 of the illustrated embodiment will not be described in detail here.

[0085] This embodiment also provides a device for extracting the boundaries of geological defect areas. This device is used to implement the above-mentioned embodiments and preferred embodiments, and details already described will not be repeated. As used below, the term "module" may refer to a combination of software and / or hardware that implements a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation using hardware, or a combination of software and hardware, is also possible and contemplated.

[0086] This embodiment provides a device for extracting the boundary of a geological defect area. Figure 9 Shown, including:

[0087] The target point cloud data determination module 901 is used to obtain the point cloud data of the tunnel where the geological defect area is located, and select the target point cloud data whose over-excavation value is greater than the first preset threshold value from the point cloud data of the tunnel where the geological defect area is located; the over-excavation value is the value of the excavation contour exceeding the preset excavation contour during the tunnel excavation process.

[0088] The boundary starting point selection module 902 is used to select multiple boundary starting points according to the dimensional direction of the target point cloud data; the multiple boundary starting points are multiple points that are within a first preset range in the first dimensional direction of the target point cloud data and within a second preset range in the second dimensional direction, and have the largest corresponding values ​​in the second dimensional direction.

[0089] The defect point search module 903 is configured to search for multiple defect points having over-excavation values ​​greater than a first preset threshold within a third preset range, with the first target boundary starting point having the smallest corresponding value in the first dimensional direction among the multiple boundary starting points as the center.

[0090] The target defect point selection module 904 is used to select a target defect point from multiple defect points; the target defect point is the defect point with the smallest azimuth angle in the counterclockwise direction, with the first target boundary starting point as the center and the first target boundary starting point to the second target boundary starting point adjacent to the first target boundary starting point as the starting direction.

[0091] The judgment module 905 is used to judge whether the target defect point exists in multiple boundary starting points. If the target defect point exists in multiple boundary starting points, the multiple boundary starting points are screened to obtain multiple target boundary points; the multiple target boundary points are multiple points in the multiple boundary starting points whose corresponding values ​​in the first dimension direction are less than or equal to the corresponding values ​​of the target defect point in the first dimension direction.

[0092] The target geological defect boundary module 906 is used to connect multiple target boundary points in sequence to obtain the target geological defect boundary.

[0093] In some optional implementations, the boundary starting point selection module 902 includes:

[0094] The average point spacing determination unit is used to obtain the average point spacing of the point cloud data of the tunnel where the geological defect area is located.

[0095] The first preset range determination unit is used to obtain the first preset range by taking the difference between the coordinates corresponding to the target point cloud data in the first dimensional direction and the average point spacing of the first preset multiple as the first boundary, and taking the sum of the coordinates corresponding to the target point cloud data in the first dimensional direction and the average point spacing of the first preset multiple as the second boundary.

[0096] The first point set determining unit is used to select a first point set within a first preset range from the point cloud data of the tunnel where the geological defect area is located.

[0097] The second preset range determination unit is used to obtain the second preset range by taking the sum of the coordinates corresponding to the target point cloud data in the second dimensional direction and the average point spacing of the second preset multiple as the third boundary, and taking the sum of the coordinates corresponding to the target point cloud data in the second dimensional direction and the average point spacing of the third preset multiple as the fourth boundary; the second preset multiple is smaller than the third preset multiple.

[0098] The boundary starting point determination unit is used to select, from the first point set, multiple points within a second preset range, with over-excavation values ​​greater than a first preset threshold, and with the largest corresponding values ​​in the second dimensional direction as multiple boundary starting points.

[0099] In some optional implementations, the defect point search module 903 includes:

[0100] The defect point search unit is used to search with the first target boundary starting point having the smallest corresponding value in the first dimensional direction among multiple boundary starting points as the center and the average point spacing of the fourth preset multiple as the radius to determine whether there are multiple points with over-excavation values ​​greater than the first preset threshold; if there are multiple points with over-excavation values ​​greater than the first preset threshold, multiple defect points are obtained; if there are no multiple points with over-excavation values ​​greater than the first preset threshold, the fourth preset multiple is adjusted until there are multiple points with over-excavation values ​​greater than the first preset threshold.

[0101] In some optional implementations, the determination module 905 includes:

[0102] The target boundary point determination unit is used to eliminate boundary starting points whose corresponding values ​​in the first dimensional direction are greater than the corresponding values ​​of the target defect point in the first dimensional direction, and form multiple target boundary points from the remaining multiple boundary starting points.

[0103] In some optional embodiments, the device for extracting the boundary of the geological defect area further includes:

[0104] A cyclic judgment unit is used to take the target defect point as the first target boundary starting point based on the target defect point not existing in multiple boundary starting points, return to the first target boundary starting point with the smallest corresponding value in the first dimensional direction among the multiple boundary starting points as the center, and search for multiple defect points with over-excavation values ​​greater than the first preset threshold within a third preset range until the target defect point exists in multiple boundary starting points.

[0105] The further functional description of each of the above modules and units is the same as that of the above corresponding embodiments and will not be repeated here.

[0106] The boundary extraction device for the geological defect area in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that executes one or more software or fixed programs, and / or other devices that can provide the above functions.

[0107] The embodiment of the present invention also provides a computer device having the above Figure 9 The boundary extraction device of the geological defect area is shown.

[0108] See also Figure 10 , Figure 10 is a structural diagram of a computer device provided by an optional embodiment of the present invention, such as Figure 10As shown, the computer device includes: one or more processors 1010, memory 1020, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components utilize different buses to communicate with each other and can be installed on a common mainboard or installed in other ways as needed. The processor can process the instructions executed in the computer device, including instructions stored in the memory or on the memory to display the graphical information of the GUI on an external input / output device (such as, a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Equally, multiple computer devices can be connected, and each device provides the necessary operation of part (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 10 A processor 1010 is taken as an example.

[0109] Processor 1010 may be a central processing unit, a network processor, or a combination thereof. Processor 1010 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a combination thereof. The programmable logic device may be a complex programmable logic device (CPLD), a field programmable gate array (FPGA), a general purpose array logic (GAL), or any combination thereof.

[0110] The memory 1020 stores instructions that can be executed by at least one processor 1010, so that the at least one processor 1010 executes the method shown in the above embodiment.

[0111] The memory 1020 may include a program storage area and a data storage area, wherein the program storage area may store an operating system and application programs required for at least one function; the data storage area may store data created based on the use of the computer device, etc. In addition, the memory 1020 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 1020 may optionally include a memory remotely located relative to the processor 1010, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0112] The memory 1020 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid-state drive; the memory 1020 may also include a combination of the above types of memory.

[0113] The computer device further includes a communication interface 1030 for the computer device to communicate with other devices or a communication network.

[0114] The embodiment of the present invention also provides a computer-readable storage medium. The above-mentioned method according to the embodiment of the present invention can be implemented in hardware, firmware, or implemented as a computer code that can be recorded in a storage medium, or implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and downloaded through a network and will be stored in a local storage medium, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state drive, etc.; further, the storage medium can also include a combination of the above-mentioned types of memory. It can be understood that a computer, a processor, a microprocessor controller or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by a computer, a processor or hardware, the method shown in the above embodiment is implemented.

[0115] A portion of the present invention may be applied as a computer program product, such as a computer program instruction, which, when executed by a computer, can call or provide the method and / or technical solution according to the present invention through the operation of the computer. Those skilled in the art should understand that the form in which the computer program instruction exists in a computer-readable medium includes, but is not limited to, a source file, an executable file, an installation package file, etc. Accordingly, the way in which the computer program instruction is executed by the computer includes, but is not limited to: the computer directly executes the instruction, or the computer compiles the instruction and then executes the corresponding compiled program, or the computer reads and executes the instruction, or the computer reads and installs the instruction and then executes the corresponding installed program. Here, the computer-readable medium may be any available computer-readable storage medium or communication medium that can be accessed by the computer.

[0116] Although the embodiments of the present invention have been described with reference to the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention. Such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A method for extracting the boundary of a geological defect area, characterized in that: The method comprises: Acquiring point cloud data of a tunnel where a geological defect area is located, and selecting target point cloud data having an overexcavation value greater than a first preset threshold value from the point cloud data of the tunnel where the geological defect area is located; the overexcavation value is a value by which an excavation contour exceeds a preset excavation contour during tunnel excavation; Selecting multiple boundary starting points according to the dimensional directions of the target point cloud data; the multiple boundary starting points are multiple points that are within a first preset range in the first dimensional direction of the target point cloud data and within a second preset range in the second dimensional direction, and have the largest corresponding values ​​in the second dimensional direction; Taking the first target boundary starting point having the smallest value corresponding to the first dimension direction among the multiple boundary starting points as the center, searching for multiple defect points whose over-excavation values ​​are greater than the first preset threshold within a third preset range; A target defect point is selected from the multiple defect points; the target defect point is a defect point with the smallest azimuth angle in the counterclockwise direction, with the first target boundary starting point as the center and the direction from the first target boundary starting point to the second target boundary starting point adjacent to the first target boundary starting point as the starting direction; determining whether the target defect point exists in the multiple boundary starting points; if the target defect point does not exist in the multiple boundary starting points, taking the target defect point as the first target boundary starting point, returning to the first target boundary starting point with the smallest value corresponding to the first dimension direction among the multiple boundary starting points as the center, and searching for multiple defect points with over-excavation values ​​greater than the first preset threshold within the third preset range until the target defect point exists in the multiple boundary starting points; If the target defect point exists in the multiple boundary starting points, the multiple boundary starting points are screened to obtain multiple target boundary points; the multiple target boundary points are multiple points in the multiple boundary starting points whose corresponding values ​​in the first dimension direction are less than or equal to the corresponding values ​​of the target defect point in the first dimension direction; The multiple target boundary points are connected in sequence to obtain a target geological defect boundary.

2. The method according to claim 1, characterized in that The step of selecting a plurality of boundary starting points according to the dimensional direction of the target point cloud data includes: Obtaining the average point spacing of the point cloud data of the tunnel where the geological defect area is located; The first preset range is obtained by taking the difference between the coordinates corresponding to the target point cloud data in the first dimensional direction and the average point spacing of the first preset multiple as a first boundary, and taking the sum of the coordinates corresponding to the target point cloud data in the first dimensional direction and the average point spacing of the first preset multiple as a second boundary; Selecting a first point set within the first preset range from the point cloud data of the tunnel where the geological defect area is located; The second preset range is obtained by taking the sum of the coordinates corresponding to the target point cloud data in the second dimensional direction and the average point spacing of the second preset multiple as a third boundary, and taking the sum of the coordinates corresponding to the target point cloud data in the second dimensional direction and the average point spacing of the third preset multiple as a fourth boundary; the second preset multiple is smaller than the third preset multiple; In the first point set, multiple points that are within the second preset range, whose over-excavation values ​​are greater than the first preset threshold, and whose corresponding values ​​in the second dimensional direction are the largest are selected as the multiple boundary starting points.

3. The method according to claim 2, characterized in that The step of searching for multiple defect points having over-excavation values ​​greater than the first preset threshold within a third preset range with the first target boundary starting point having the smallest corresponding value in the first dimensional direction among the multiple boundary starting points as the center includes: Taking the first target boundary starting point with the smallest corresponding value in the first dimension among the multiple boundary starting points as the center and the average point spacing that is a fourth preset multiple as the radius, a search is performed to determine whether there are multiple points whose over-excavation values ​​are greater than the first preset threshold; If there are multiple points where the over-excavation value is greater than the first preset threshold, the multiple defect points are obtained.

4. The method according to claim 3, characterized in that The method further comprises: If there are not multiple points where the over-digging value is greater than the first preset threshold, the fourth preset multiple is adjusted until there are multiple points where the over-digging value is greater than the first preset threshold.

5. The method according to any one of claims 1 to 4, characterized in that The screening of the plurality of boundary starting points to obtain a plurality of target boundary points includes: The boundary starting points whose corresponding values ​​in the first dimension direction are greater than the corresponding values ​​of the target defect point in the first dimension direction are eliminated, and the remaining multiple boundary starting points form the multiple target boundary points.

6. A device for extracting the boundary of a geological defect area, characterized in that: The device comprises: a target point cloud data determination module, configured to obtain point cloud data of a tunnel in which a geological defect region is located, and select target point cloud data having an overexcavation value greater than a first preset threshold value from the point cloud data of the tunnel in which the geological defect region is located; the overexcavation value is a value at which the excavation contour exceeds a preset excavation contour during tunnel excavation; a boundary starting point selection module, configured to select a plurality of boundary starting points according to the dimensional directions of the target point cloud data; the plurality of boundary starting points being a plurality of points that are within a first preset range in the first dimensional direction of the target point cloud data and within a second preset range in the second dimensional direction, and have the largest corresponding value in the second dimensional direction; A defect point search module is configured to search for a plurality of defect points whose over-excavation values ​​are greater than the first preset threshold within a third preset range, with a first target boundary starting point having the smallest corresponding value in the first dimensional direction among the plurality of boundary starting points as the center; a target defect point selection module, configured to select a target defect point from the plurality of defect points; the target defect point being a defect point with the smallest azimuth angle in the counterclockwise direction, with the first target boundary starting point as the center and the direction from the first target boundary starting point to the second target boundary starting point adjacent to the first target boundary starting point as the starting direction; a loop judgment unit configured to, based on the target defect point not existing in the multiple boundary starting points, use the target defect point as a first target boundary starting point, return to the step of searching for multiple defect points having over-excavation values ​​greater than a first preset threshold within a third preset range with the first target boundary starting point having the smallest corresponding value in the first dimension direction among the multiple boundary starting points as the center, until the target defect point exists in the multiple boundary starting points; a determination module, configured to determine whether the target defect point exists in the multiple boundary starting points, and if the target defect point exists in the multiple boundary starting points, filter the multiple boundary starting points to obtain multiple target boundary points; the multiple target boundary points are multiple points in the multiple boundary starting points whose corresponding values ​​in the first dimension direction are less than or equal to the corresponding values ​​of the target defect point in the first dimension direction; The target geological defect boundary module is used to connect the multiple target boundary points in sequence to obtain the target geological defect boundary.

7. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the boundary extraction method of the geological defect area according to any one of claims 1 to 5 by executing the computer instructions.

8. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the boundary extraction method of the geological defect area according to any one of claims 1 to 5.

9. A computer program product, characterized in that The method comprises computer instructions for causing a computer to execute the method for extracting the boundary of a geological defect area according to any one of claims 1 to 5.

Images