Geological defect area boundary extraction method, device, equipment, medium and product

Abstract

The invention relates to the technical field of tunnel engineering, and discloses a boundary extraction method and device of a geological defect area, equipment, a medium and a product. The boundary extraction method of the geological defect area comprises the steps of obtaining point cloud data of a tunnel where the geological defect area is located, target point cloud data with the over-excavation value larger than a first preset threshold value are selected from the point cloud data of the tunnel where the geological defect area is located; selecting a plurality of boundary starting points according to the dimension direction of the target point cloud data; searching a plurality of defect points; selecting a target defect point from the plurality of defect points; judging whether the target defect point exists in the plurality of boundary starting points or not, and if yes, screening the plurality of boundary starting points to obtain a plurality of target boundary points; and sequentially connecting the plurality of target boundary points to obtain a target geological defect boundary. According to the method, the points in the point cloud data of the tunnel where the geological defect area is located are extracted and judged, so that the accuracy of boundary extraction of the geological defect area is improved.

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 is used in tunnel engineering, if unfavorable geological bodies are encountered, there will be geological defect over-excavation. Geological defect over-excavation refers to the situation where the excavation cannot be carried out according to the pre-designed excavation method and size due to defects in geological conditions (such as faults, fractured zones, weak interlayers and other unfavorable geological conditions), resulting in the excavation contour significantly exceeding the pre-designed excavation contour. When geological defect over-excavation occurs, it is necessary to measure the geological defect area and accurately divide the boundaries in order to handle it and settle the costs.

[0003] In the related art, the method for identifying geological defect areas is to use a measuring instrument to obtain multiple measuring points in the geological defect area, estimate the geological defect boundary based on the multiple measuring points, and delineate the geological defect area. This method has the problems of sparse measuring points and randomness, resulting in the problem that the obtained geological defect boundary is not accurate enough and the delineated geological defect area does not conform to the actual situation. 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 inaccuracy of the geological defect boundary 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, including: 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 with the smallest corresponding numerical value in the first dimensional direction among the multiple boundary starting points as the center. A plurality of defect points whose over-excavation values ​​are 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, 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 azimuth; 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 connected in sequence to obtain a target geological defect boundary.

[0006] The present invention selects 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, screens out target point cloud data with significant over-excavation characteristics, excludes interference information and error data, and also provides a more reliable data basis for more accurately determining the boundary starting point. The present invention selects multiple boundary starting points based on the dimensional direction and the preset range, and can accurately determine 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 takes the first target boundary starting point with the smallest corresponding value in the first dimensional direction among multiple boundary starting points as the center, searches for multiple defect points with an over-excavation value greater than the first preset threshold value within the third preset range, and gradually explores the geological defect area in an orderly manner to ensure that no defect points are missed. At the same time, the setting of the third preset range further limits the search area, improves the search efficiency while ensuring comprehensiveness, avoids blind search in an overly large range, and makes the search of defect points more accurate and efficient. The present invention selects target defect points by clear conditions such as azimuth, and can determine 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, 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 finally 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 being connected in sequence 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 the accuracy of geological defect boundary identification is improved.

[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 the 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 ​​less 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, taking the first target boundary starting point with 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: taking the first target boundary starting point with the smallest corresponding value in the first dimensional direction among multiple boundary starting points as the center, and 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 among 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 corresponding value in the first dimensional 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 among the multiple boundary starting points.

[0012] If the target defect point does not exist in multiple boundary starting points, the present invention uses the target defect point as the first target boundary starting point, returns 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 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 the entire area where geological defects may exist is fully covered by a cycle, thereby improving the completeness and accuracy of target geological defect extraction.

[0013] In a second aspect, the present invention provides a device for extracting the boundary of 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 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 the value at which the excavation contour exceeds 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 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; 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, and 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 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 are communicatively connected to each other, the memory stores computer instructions, and the processor executes the boundary extraction method of the geological defect area of ​​the first aspect or any corresponding embodiment thereof by executing the computer instructions.

[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 boundaries of geological defect areas 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 the related technologies, the drawings required for use in the specific embodiments or the related technical descriptions will be briefly introduced below. 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 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 4 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 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 is a schematic diagram of a target geological defect boundary according to an embodiment of the present invention.

[0025] Figure 8 It 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] Fig. 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] Fig.10It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0028] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the 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 creative work are within the scope of protection of the present invention.

[0029] When drilling and blasting is used in tunnel construction, if unfavorable geological bodies are encountered, geological defect over-excavation will occur. Geological defect over-excavation refers to the inability to carry out construction according to the designed excavation method and size during the excavation process due to defects in geological conditions (such as faults, broken zones, weak interlayers and other unfavorable geological conditions), resulting in over-excavation. In order to prevent geological defect over-excavation, it is necessary to measure the geological defect area and accurately divide the boundaries so as to manage it.

[0030] The related technology mainly uses a total station or other measuring instruments to obtain multiple measuring points in the geological defect area, draws the cross-sectional area of ​​the defect over-excavation area at preset intervals, and manually estimates the geological defect boundary. The related technology has the disadvantages of sparse measuring points and randomness, which leads to inaccurate determination of the over-excavation area boundary, cumbersome calculation process and low efficiency.

[0031] The embodiment of the present invention provides a method for extracting the boundary of a geological defect area, which achieves the effect of improving 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 boundaries 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 a computer device. Figure 1 is a flow chart of a method for extracting a boundary of a geological defect area according to an embodiment of the present invention. Figure 1 As shown, the process includes the following steps: 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 a value when the excavation contour exceeds the preset excavation contour during the tunnel excavation process.

[0034] 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, 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, but 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.

[0035] 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.

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

[0037] Step S102, 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 whose over-excavation values ​​are less than a first preset threshold, and whose corresponding values ​​in the second dimensional direction are the largest.

[0038] 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 a plane rectangular coordinate system. The corresponding value in the second dimensional direction is a 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.

[0039] In some optional embodiments, the first preset range is the length range between two preset coordinate points in the first dimensional direction. Specifically, the first preset range is a range with 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 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. 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.

[0040] In some optional embodiments, the second preset range is the length range between two preset coordinate points in the second dimension direction. 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 direction 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 direction 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.

[0041] 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 can be any one value from 1 to 5, the initial value of i is 0, and can be increased sequentially.

[0042] 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 the 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 ​​less than the first preset threshold, and with the largest corresponding values ​​in the second dimensional direction as multiple boundary starting points.

[0043] The average point spacing is determined by using a sampling method or a point density calculation method. Specifically, the process of determining the average point spacing by using the sampling method is as follows: randomly extract multiple point samples from the point cloud data of the tunnel where the geological defect area is located, calculate the distance between the multiple point samples, and perform statistical analysis on the distance between the multiple point samples to obtain the average point spacing. The process of determining the average point spacing by using the point density calculation method is as follows: statistically calculate the number of point clouds and the area of ​​the region within the point cloud data of the tunnel where the geological defect area is located, determine the point density based on the area and the number of point clouds, and estimate the average point spacing by the point density.

[0044] For example, Figure 3 As shown in FIG. 1 , 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 The multiple dark black solid points shown in the figure are selected from 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 less 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.

[0045] In the 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.

[0046] In some optional embodiments, when the first points are concentrated within the second preset range but there are no points with over-excavation values ​​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 points are concentrated within the second preset range and there are points with over-excavation values ​​less than the first preset threshold.

[0047] 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 whose over-excavation values ​​are greater than the first preset threshold within a third preset range.

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

[0049] In some optional embodiments, taking the first target boundary starting point with the smallest corresponding value in the first dimensional direction among 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 includes: taking the first target boundary starting point with 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, and judging whether there are multiple points whose over-excavation values ​​are greater than the first preset threshold; if there are multiple points whose over-excavation values ​​are greater than the first preset threshold, multiple defect points are obtained.

[0050] 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.

[0051] 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 whose over-excavation values ​​are 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 whose over-excavation values ​​are greater than the first preset threshold, increase the value of α until there are multiple points whose over-excavation values ​​are greater than the first preset threshold.

[0052] Step S104, selecting a target defect point from among multiple defect points; the target defect point is a defect point with the smallest azimuth 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.

[0053] 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 value 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.

[0054] For example, Figure 4 As shown in FIG. 1 , it is a schematic diagram of the target defect point selection process, where the point set m is a set of multiple defect points, with the first target boundary starting point Qn as the center point, and the direction from Qn to the second target boundary starting point Qn-1 as the starting direction. The azimuths 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.

[0055] 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.

[0056] 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.

[0057] 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 the 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 the multiple target boundary points is (Q1, ..., Qn-1, Qn), where Q1 is Qi in the point set composed of multiple boundary starting points.

[0058] 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 numerical 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.

[0059] 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 step of returning 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 searching for multiple defect points whose over-excavation values ​​are greater than the first preset threshold within a third preset range, until the target defect point exists in multiple boundary starting points, can ensure that no possible defect points are missed, and the entire area where geological defects may exist is fully covered by looping, thereby improving the completeness and accuracy of target geological defect extraction.

[0060] Step S106, connecting multiple target boundary points in sequence to obtain a target geological defect boundary.

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

[0062] 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.

[0063] The embodiment of the present invention selects 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, screens out target point cloud data with significant over-excavation characteristics, excludes interference information and error data, and also provides a more reliable data basis for more accurately determining the boundary starting point. The embodiment of the present invention selects multiple boundary starting points based on the dimensional direction and the preset range, and can accurately determine 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 embodiment of the present invention takes the first target boundary starting point with the smallest corresponding value in the first dimensional direction among multiple boundary starting points as the center, searches for multiple defect points with an over-excavation value greater than the first preset threshold value within the third preset range, and gradually explores the geological defect area in an orderly manner to ensure that no defect points are missed. At the same time, the setting of the third preset range further limits the search area, improves the search efficiency while ensuring comprehensiveness, avoids blind search in an overly large range, and makes the search of defect points more accurate and efficient. The embodiment of the present invention selects the target defect point by clear conditions such as azimuth, and can determine the most representative and critical target defect point among multiple defect points. This selection method ensures the uniqueness and certainty of the target defect point. 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, 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 finally 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.

[0064] In this embodiment, a method for extracting the boundary of a geological defect area is provided, which can be used in a computer device. Figure 6 FIG. 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. Figure 6 As shown, the process includes the following steps: Step S601, obtaining point cloud data of the tunnel where the geological defect area is located, and projecting the point cloud data onto a preset section to obtain a projection result.

[0065] 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 by using three-dimensional laser scanning technology.

[0066] 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.

[0067] 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.

[0068] Step S602, calculating the distance between the point cloud data and the projection result to obtain the over-excavation value.

[0069] 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.

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

[0071] In some optional implementations, the plane may be selected as a specific plane perpendicular to the tunnel axis, or a plane having a certain angle relationship with the tunnel direction, etc. Each point in the point cloud data is projected onto the selected plane using a vector projection method.

[0072] Step S604: coloring and rendering the target point cloud model according to the over-excavation value to obtain a coloring and rendering result.

[0073] 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.

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

[0075] 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.

[0076] In the embodiment of the present invention, Figure 7 The figure shows the target geological defect boundary diagram. 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 being processed by the embodiment of the present invention, the three-dimensional point cloud data of the tunnel is transformed into a target point cloud model with target geological defect boundaries. Figure 7In the figure, the irregular solid and dotted lines in the target point cloud model are the target geological defect boundaries.

[0077] In this embodiment, a method for extracting the boundary of a geological defect area is provided, which can be used in a computer device. Figure 8 FIG. 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. Figure 8 As shown, the process includes the following steps: The first method for extracting the target geological defect boundary is: collecting the three-dimensional point cloud of the tunnel, converting the coordinates of the three-dimensional point cloud of the tunnel into the engineering coordinate system, calculating the over-excavation value of the three-dimensional point cloud of the tunnel, mapping the three-dimensional point cloud of the tunnel to a plane, performing visual rendering based on the over-excavation value, and manually or automatically drawing the geological defect boundary according to the rendering result, wherein the geological defect boundary here is the target geological defect boundary.

[0078] 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 effectiveness 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.

[0079] In the 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, and the boundary starting point is extracted according to the average point spacing, including: selecting the 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 the first preset range in the first dimensional direction of the target point cloud data, within the second preset range in the second dimensional direction, and have the largest corresponding values ​​in the second dimensional direction. For details, please refer to Figure 1 Step S102 of the illustrated embodiment will not be described in detail here.

[0080] In some optional implementations, 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 multiple boundary starting points as the center, searching for multiple defect points with over-excavation 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 among the multiple defect points, taking 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. For details, please refer to Figure 1 Step S103 and step S104 of the illustrated embodiment will not be described in detail here.

[0081] 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, and 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 refer to Figure 1 Step S105 of the illustrated embodiment will not be described in detail here.

[0082] In this embodiment, a device for extracting the boundary of a geological defect area is also provided, which is used to implement the above-mentioned embodiments and preferred embodiments, and the descriptions that have been made will not be repeated. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0083] This embodiment provides a device for extracting the boundary of a geological defect area. Fig. 9 As shown, including: 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.

[0084] 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, within a second preset range in the second dimensional direction, and have the largest corresponding values ​​in the second dimensional direction.

[0085] The defect point search module 903 is used to search for multiple defect points whose over-excavation values ​​are 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.

[0086] The target defect point selection module 904 is used to select a target defect point from 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 first target boundary starting point to the second target boundary starting point adjacent to the first target boundary starting point as the starting direction.

[0087] 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.

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

[0089] In some optional implementations, the boundary starting point selection module 902 includes: 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.

[0090] 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.

[0091] The first point set determination 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.

[0092] 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.

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

[0094] In some optional implementations, the defect point finding module 903 includes: A 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 a fourth preset multiple as the radius to determine whether there are multiple points whose over-excavation values ​​are greater than the first preset threshold; if there are multiple points whose over-excavation values ​​are greater than the first preset threshold, multiple defect points are obtained; if there are no multiple points whose over-excavation values ​​are greater than the first preset threshold, the fourth preset multiple is adjusted until there are multiple points whose over-excavation values ​​are greater than the first preset threshold.

[0095] In some optional implementations, the determination module 905 includes: The target boundary point determination unit is used to eliminate the 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.

[0096] In some optional implementations, the device for extracting the boundary of the geological defect area further includes: A cyclic judgment unit is used to take the target defect point as the first target boundary starting point according to the fact that the target defect point does not exist 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.

[0097] 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.

[0098] The boundary extraction device of 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.

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

[0100] See also Fig.10 , Fig.10 is a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention, such as Fig.10 As 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 are connected to each other using different buses for communication, 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. Similarly, multiple computer devices can be connected, and each device provides part of the necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Fig.10 A processor 1010 is taken as an example.

[0101] The processor 1010 may be a central processing unit, a network processor or a combination thereof. The processor 1010 may further include a hardware chip. The hardware chip may be a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic or any combination thereof.

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

[0103] The memory 1020 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function; the data storage area may store data created according to 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 arranged 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.

[0104] 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.

[0105] The computer device also includes a communication interface 1030, which is used for the computer device to communicate with other devices or a communication network.

[0106] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, 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 hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. 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.

[0107] A part 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 existence of the computer program instruction in a computer-readable medium includes, but is not limited to, a source file, an executable file, an installation package file, etc., and 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 accessible to the computer.

[0108] Although the embodiments of the present invention have been described in conjunction with 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, and 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: Acquire point cloud data of the tunnel where the geological defect area is located, and select target point cloud data whose over-excavation value is 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 a value at which the excavation contour exceeds the preset excavation contour during the tunnel excavation process; Select 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 corresponding value in the first dimensional 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; Determine whether the target defect point exists in the multiple boundary starting points. If the target defect point exists in the 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 dimensional direction are less than or equal to the corresponding values ​​of the target defect point in the first dimensional 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 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; 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 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; In the first point set, multiple points that are within the second preset range, whose over-digging values ​​are smaller 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 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 and searching for multiple defect points whose over-excavation values ​​are greater than the first preset threshold within a third preset range comprises: 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 and the average point spacing of the fourth preset multiple as the radius, searching is performed to determine whether there are multiple points whose over-digging values ​​are greater than the first preset threshold; If there are multiple points whose over-digging values ​​are 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 a plurality of points where the over-digging value is greater than the first preset threshold, the fourth preset multiple is adjusted until there are a plurality of 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. The method according to any one of claims 1 to 4, characterized in that The method further comprises: If the target defect point does not exist among the multiple boundary starting points, take the target defect point as the first target boundary starting point, return to the first target boundary starting point with the smallest corresponding value in the first dimensional direction among the multiple boundary starting points, and search for multiple defect points whose over-excavation values ​​are greater than the first preset threshold within a third preset range until the target defect point exists among the multiple boundary starting points.

7. A device for extracting the boundary of a geological defect area, characterized in that: The device comprises: The target point cloud data determination module 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 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; A boundary starting point selection module 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; A defect point search module, configured to search for multiple defect points whose over-excavation values ​​are greater than the first preset threshold value 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; A target defect point selection module is used to select a target defect point 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 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 judgment module, used for judging 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, 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 dimensional direction are less than or equal to the corresponding values ​​of the target defect point in the first dimensional 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.

8. 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 6 by executing the computer instructions.

9. 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 a geological defect area according to any one of claims 1 to 6.

10. 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 6.

Images