Method and system for determining geological structure line of area to be explored, computer storage medium

Abstract

Embodiments of the present application relate to the field of image data processing, and particularly relate to a method and system for determining a geological structure line of a to-be-explored area and a computer storage medium. The method and system can determine the dominant direction of the topography of each position and the corresponding response amplitude by extracting the response amplitude of the topography of the to-be-explored area to geological structures in different directions from the surface elevation data, thereby accurately determining the main structure area and the main geological structure line therein, and then extracting the secondary geological structure line on this basis, realizing the extraction of a complete geological structure line network, so that the extraction result can truly reflect the distribution of the geological structure line of the to-be-explored area, and can improve the objectivity and resolution of the extraction of the structure line, and avoid the confusion of the main and secondary structure lines. Compared with the subjectivity and limitations of traditional manual extraction, the method can eliminate topographic noise and highlight the topographic linear features consistent with the structure trend, and objectively identify geological structure lines of different scales and levels.

Description

Technical Field

[0001] The embodiments of this application relate to the field of image data processing, and in particular to a method and system for determining geological structural lines of an area to be explored, as well as a computer storage medium. Background Technology

[0002] The statements herein are provided merely as background information in connection with this application and do not necessarily constitute prior art.

[0003] Geological structural lines include fault zones, linear structures, lithological contact boundaries, and fold axes. In the field of mineral resource exploration, structural lines are often the direct manifestation of ore-controlling structures. In earthquake disaster early warning, identifying active faults can enable geological risk assessment. Therefore, geological structural lines can depict the trajectory of crustal movement, reveal deep mineralization patterns, and evaluate engineering geological stability. Summary of the Invention

[0004] A brief overview of this application is provided below to offer a basic understanding of certain aspects thereof. It should be understood that this overview is not an exhaustive summary of the application. It is not intended to identify key or essential parts of the application, nor is it intended to limit its scope. Its purpose is merely to present certain concepts in a simplified form as a prelude to the more detailed description that follows.

[0005] This application provides a method for determining the geological structural lines of an area to be explored, comprising the following steps: S10: determining the surface elevation data of the area to be explored; S20: determining the response amplitude of the topography at each location in the area to be explored to geological structural lines in different predetermined directions based on the surface elevation data; S30: determining the response amplitude distribution corresponding to the dominant topographic direction at each location based on the response amplitude determined in step S20; S40: determining the main structural region of the area to be explored based on the response amplitude distribution determined in step S30; S50: determining the main geological structural lines of the area to be explored based on the main structural regions; S60: determining the secondary geological structural lines of the area to be explored based on the main geological structural lines determined in step S50; S70: determining the geological structural lines of the area to be explored based on the main geological structural lines and the secondary geological structural lines.

[0006] The method for determining geological structural lines in an area to be explored, provided in the embodiments of this application, extracts the response amplitude of the terrain to geological structures in different directions from surface elevation data. This allows for the determination of the dominant terrain direction and its corresponding response amplitude at each location, thereby accurately identifying the main structural region and its primary geological structural lines. Secondary geological structural lines are then extracted based on this, achieving a complete network of geological structural lines. This ensures that the extraction results accurately reflect the distribution of geological structural lines throughout the entire area to be explored, improving the objectivity and resolution of structural line extraction and avoiding confusion between primary and secondary structural lines. Compared to the subjectivity and limitations of traditional manual extraction, the method provided in the embodiments of this application can eliminate terrain noise and highlight linear features of the terrain consistent with structural trends, enabling the objective identification of geological structural lines at different scales and levels, thus improving the accuracy, efficiency, and repeatability of geological structural line extraction.

[0007] This application also provides a system for determining geological structural lines in an area to be explored, comprising: a sensor configured to determine surface elevation data of the area to be explored; and a processor configured to: determine, based on the surface elevation data, the response amplitude of the terrain at each location in the area to be explored to geological structural lines in different predetermined directions; determine, based on the response amplitude, the distribution of the response amplitude corresponding to the dominant direction of the terrain at each location; determine, based on the response amplitude distribution, the main structural region of the area to be explored; determine, based on the main structural region, the main geological structural line of the area to be explored; determine, based on the main geological structural line, the secondary geological structural line of the area to be explored; and determine, based on the main geological structural line and the secondary geological structural line, the geological structural line of the area to be explored.

[0008] The system for determining geological structural lines in an area to be explored, provided in the embodiments of this application, extracts the response amplitude of the terrain to geological structures in different directions from surface elevation data. This allows for the determination of the dominant topographic direction and its corresponding response amplitude at each location, thereby accurately identifying the main structural region and its primary geological structural lines. Secondary geological structural lines are then extracted based on this, achieving a complete network of geological structural lines. This ensures that the extraction results accurately reflect the distribution of geological structural lines throughout the entire area to be explored, improving the objectivity and resolution of structural line extraction and avoiding confusion between primary and secondary structural lines. Compared to the subjectivity and limitations of traditional manual extraction, the method provided in the embodiments of this application can eliminate topographic noise and highlight linear topographic features consistent with structural trends, enabling the objective identification of geological structural lines at different scales and levels, thus improving the accuracy, efficiency, and repeatability of geological structural line extraction.

[0009] Embodiments of this application also provide a computer storage medium storing a computer program, the computer program including program instructions, which, when executed by a processor, perform the aforementioned method. Attached Figure Description

[0010] To further illustrate the above and other advantages and features of this application, the specific embodiments of this application will be described in more detail below with reference to the accompanying drawings. The drawings, together with the following detailed description, are included in and form a part of this specification. Elements having the same function and structure are indicated by the same reference numerals. It should be understood that these drawings only depict typical examples of this application and should not be considered as limiting the scope of this application.

[0011] Figure 1 This is a schematic diagram of the geological structure lines of the area to be explored, determined using the method provided in the embodiments of this application; Figure 2 This is a schematic diagram of the response amplitude distribution corresponding to the dominant terrain direction at each location in the method provided by the embodiments of this application; Figure 3 This is a schematic diagram of the surface elevation data preprocessing method provided in the embodiments of this application; Figure 4 This is a schematic diagram of the main structural region obtained using the method provided in the embodiments of this application; Figure 5 This is a schematic diagram showing the distribution of pixel lines obtained by processing the main structure region using the method provided in the embodiments of this application. Detailed Implementation

[0012] Exemplary embodiments of this application will be described below with reference to the accompanying drawings. For clarity and brevity, not all features of actual implementations are described in the specification. However, it should be understood that many implementation-specific decisions must be made in the development of any such actual embodiment to achieve the developer's specific goals, such as complying with constraints related to the system and business, and these constraints may vary depending on the implementation. Furthermore, it should be understood that while development work can be very complex and time-consuming, such development work is merely a routine task for those skilled in the art who benefit from the content of this application.

[0013] It should also be noted that, in order to avoid obscuring this application with unnecessary details, only the equipment structure and / or processing steps closely related to the solution according to this application are shown in the accompanying drawings, while other details that are not closely related to this application are omitted.

[0014] The following disclosure provides several different implementations or examples for carrying out this application. To simplify the disclosure of this application, specific examples of components and methods are described below. Of course, these are merely examples and are not intended to limit this application. In the description of the embodiments of this application, "multiple" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0015] Currently, the acquisition of geological structural lines mainly relies on manual visual extraction by experienced experts, which must be combined with the regional geological background and comprehensively judged based on tonal features, texture orientation, and linear discontinuities of topography in the image. However, this method suffers from drawbacks such as extremely low efficiency, high subjectivity and uncertainty in the extraction results, and a lack of objective quantitative evaluation standards. Existing automatic or semi-automatic extraction methods, when faced with complex remote sensing geoscientific scenes, suffer from insurmountable shortcomings such as lack of direction awareness, fragmented extracted structural lines, and output extraction results in formats that are difficult to directly apply.

[0016] To address the aforementioned problems, one embodiment of this application provides a method for determining the geological structural lines of an area to be explored, comprising the following steps: S10: determining the surface elevation data of the area to be explored; S20: determining the response amplitude of the terrain at each location in the area to be explored to geological structural lines in different predetermined directions based on the surface elevation data; S30: determining the response amplitude distribution corresponding to the dominant direction of the terrain at each location based on the response amplitude determined in step S20; S40: determining the main structural region of the area to be explored based on the response amplitude distribution determined in step S30; S50: determining the main geological structural lines of the area to be explored based on the main structural regions; S60: determining the secondary geological structural lines of the area to be explored based on the main geological structural lines determined in step S50; S70: determining the geological structural lines of the area to be explored based on the main geological structural lines and the secondary geological structural lines.

[0017] The method for determining geological structural lines in an area to be explored, provided in the embodiments of this application, extracts the response amplitude of the terrain to geological structures in different directions from surface elevation data. This allows for the determination of the dominant terrain direction and its corresponding response amplitude at each location, thereby accurately identifying the main structural region and its primary geological structural lines. Secondary geological structural lines are then extracted based on this, achieving a complete network of geological structural lines. This ensures that the extraction results accurately reflect the distribution of geological structural lines throughout the entire area to be explored, improving the objectivity and resolution of structural line extraction and avoiding confusion between primary and secondary structural lines. Compared to the subjectivity and limitations of traditional manual extraction, the method provided in the embodiments of this application can eliminate terrain noise and highlight linear features of the terrain consistent with structural trends, enabling the objective identification of geological structural lines at different scales and levels, thus improving the accuracy, efficiency, and repeatability of geological structural line extraction.

[0018] Figure 1 This is a schematic diagram of the geological structural lines of the area to be explored, determined using the method provided in the embodiments of this application, as shown below. Figure 1 As shown in the figure, the green lines represent the construction lines determined by the method provided in the embodiments of this application, and the red lines represent the construction lines determined by the traditional manual visual extraction method. Compared with the traditional visual extraction method, the construction lines determined by the method provided in the embodiments of this application are thinner, more numerous and denser, and their extension direction is basically consistent with the construction lines of the manual visual extraction results. This indicates that the method provided in the embodiments of this application can extract smaller secondary construction lines and obtain more local features in the area to be explored, thus improving the extraction accuracy.

[0019] In some embodiments, step S50 further includes the following steps: S51: determining the pixel lines of the main structure; S52: determining the local connectivity relationships of the pixels in the pixel lines; S53: determining the weights of the local connectivity relationships; S54: determining the endpoints of the main geological structural lines based on the weights of the local connectivity relationships; S55: determining the main geological structural lines based on the endpoints. By quantifying the connectivity relationships between pixels in the pixel lines of the main structure and assigning weights, the determined endpoints can reflect the overall connectivity structure of the pixel lines, thereby reflecting the true start and end positions of the structural lines and effectively avoiding false branches or breaks caused by noise or local discontinuities.

[0020] In some embodiments, in step S53, the weights of local connectivity relationships are determined as follows: ;in, and This indicates two pixels that are adjacent to each other. and These represent the spatial coordinates of two adjacent pixels; The weights represent the local connectivity relationships. By quantizing the spatial distance between pixels as the weights of local connectivity relationships, connections in different directions can be distinguished more accurately, improving the geometric accuracy of the main geological structural lines extraction and the accuracy of continuity determination.

[0021] In some embodiments, during step S52, when determining the local connectivity of a pixel, the connectivity can be determined in an 8-neighborhood geometric space, which refers to the set of pixels in eight adjacent directions surrounding a given pixel. 8-neighborhood connectivity considers the adjacency relationships of a pixel in the top, bottom, left, right, and diagonal directions. Compared to 4-neighborhood, it can more comprehensively capture the changes in the direction of the construction line, especially the tilt direction. This allows the weights of the subsequently determined local connectivity relationships to more accurately reflect the actual distance, thereby improving the continuity and completeness of the construction line extraction and avoiding breaks caused by ignoring diagonal connections.

[0022] In some embodiments, in step S54, the endpoints of the main geological structural line are determined as follows: ;in, and Indicates the endpoints of the main geological structural lines; represents the starting pixel; x represents any pixel; V represents the set of all pixels. This represents the sum obtained by accumulating weights along the local connectivity relationship from the starting pixel to any other pixel; This represents the sum of weights accumulated along the local connectivity from the endpoint to any pixel. This represents the pixel that maximizes the sum of the accumulated weights.

[0023] By accumulating weights along local connectivity starting from the initial pixel to determine the farthest point as one endpoint, and repeating the above operation with this farthest point as the starting pixel to determine the farthest point as the other endpoint, the two endpoints with the greatest geometric distance in the main geological structural line can be automatically determined. The above method utilizes the connectivity weights between all pixels along the entire structural line, effectively eliminating the interference of local short branches or noisy pixels on the determination of the endpoints of the main geological structural line, so that the final obtained main geological structural line has the maximum continuous length and clear start and end positions in space, reflecting the main extension direction and range of the actual geological body.

[0024] Figure 5 This is a schematic diagram showing the distribution of pixel lines obtained by processing the main structure region using the method provided in the embodiments of this application, as shown below. Figure 5 As shown, in step S51, the edge pixels of the main structure can be removed along the extension direction of the main structure to compress the main structure into pixel lines with a width of one pixel, so as to facilitate the subsequent extraction of geological structure lines.

[0025] In some embodiments, in step S20, the response amplitude of the terrain at each location to geological structural lines in different predetermined directions is determined in the following manner: ;in, These are angles from different predetermined directions; Represents the spatial pixel coordinates of each location in the surface elevation data; This indicates the terrain at each location for different predetermined directions. The response amplitude of geological structural lines; and These represent the real and imaginary parts of the response amplitude, respectively. This method allows for the acquisition of phase-independent response amplitudes, enabling stable quantification of the intensity of structural line features in different directions and ensuring the comparability and directional selectivity of response amplitudes in subsequent steps.

[0026] Figure 3 This is a schematic diagram of the surface elevation data preprocessing method provided in the embodiments of this application. In some embodiments, such as... Figure 3 As shown, in step S20, the surface elevation data can be linearly stretched before determining the response amplitude of the terrain at each location to geological structural lines in different predetermined directions, in order to eliminate extreme value interference and enhance local contrast.

[0027] In some embodiments, in step S20, the stretched surface elevation data may be filtered.

[0028] In some embodiments, in step S20, the filter may take the following form: ;in, It is the angle of the predetermined direction. ; Represents the spatial pixel coordinates of each location in the surface elevation data; Represents spatial pixel coordinates through The corresponding coordinates after the angle rotation; This indicates the spatial aspect ratio and controls the elliptical shape of the filter in space; it can take a value of 1.0. It is the standard deviation of the Gaussian envelope function, which controls the effective spatial range of the filter; The spatial frequency of the sinusoidal wave carrier is used to determine the scale of the geological texture that can be extracted, and can be 0.1. It is the phase shift of a sine function, which includes both 0 and 1. Two phases; It is the imaginary unit; Indicates direction as The filtering kernel function can detect texture features in the horizontal, diagonal, and vertical directions by setting the predetermined directional angle to 0°, 45°, and 90° respectively, covering the main extension directions of the terrain, enhancing the perception of anisotropic geological textures, and avoiding the loss of subtle structural features with specific orientations in complex terrain backgrounds.

[0029] In some embodiments, spatial pixel coordinates are determined in the following manner through Corresponding coordinates after angle rotation ,in , .

[0030] In some embodiments, in step S30, the response amplitude distribution corresponding to the dominant terrain direction at each location is determined in the following manner: ;in, Indicates angles in different predetermined directions; Represents the spatial pixel coordinates of each location in the surface elevation data; This indicates the terrain at each location for different predetermined directions. The response amplitude of geological structural lines; Indicates that all directional response amplitude Operations on the maximum value; This represents the set of maximum response amplitudes for each location in all predetermined directions, and represents the response amplitude distribution for each location in the terrain-dominant direction.

[0031] Figure 2 This is a schematic diagram of the response amplitude distribution corresponding to the dominant terrain direction at each location in the method provided by the embodiments of this application, as shown below. Figure 2 As shown, by determining the maximum response amplitude of each location in all predetermined directions as the response amplitude corresponding to the dominant topographic direction of that location, the most significant geological structural direction information at each pixel can be highlighted, while suppressing the responses of other non-dominant directions. The set of maximum response amplitudes helps to identify the main structural area in the subsequent process, further eliminating noise interference in the extraction of geological structural lines.

[0032] In some embodiments, step S40 further includes the following steps: S41: determining a predetermined threshold based on the response amplitude distribution determined in step S30; S42: removing the background texture of the response amplitude distribution based on the predetermined threshold to obtain a binary construction map; S43: determining connected regions with main structures based on the binary construction map determined in step S42; S44: determining discrete regions near the connected regions based on the connected regions; S45: determining the main construction region of the area to be explored based on the connected regions and the discrete regions.

[0033] By converting the response amplitude distribution into a binary map according to a predetermined threshold, it is possible to effectively separate the main structure with strong response from the background noise or non-structural texture with weak response, so as to more clearly identify the connected regions where the main structure exists. Further, the discrete regions near the connected regions are determined and included in the main structure region. This can preserve the continuity of the main structure and avoid missing the surrounding related secondary or isolated structural units, so that the finally determined main structure region is more complete and conforms to the actual distribution of the structure in the area to be explored.

[0034] In some embodiments, in step S43, the connected regions with a main structure can be determined in the following manner: ;in, Indicates the first A connected region with a main structure; Indicates the first The area of ​​each connected region; The union symbol represents merging all objects that satisfy the area condition into a single entity; B represents a set of connected regions with a main structure. This represents the threshold used to determine the connected regions containing the main structure. By retaining connected regions with an area greater than or equal to the set threshold as the connected regions of the main structure, interference from small discrete regions can be avoided in the identification of the main structure, thus improving the reliability of the determination of the main structure region.

[0035] In some embodiments, in step S44, the discrete regions near the connected region can be determined in the following manner: ;in, A binary matrix representing a connected region with a main structure; Represents the spatial coordinates of the binary matrix; It is the standard deviation of the two-dimensional Gaussian convolution kernel; It is a two-dimensional discrete space convolution operator; It is the coordinates after convolution processing. The output value at that location.

[0036] Figure 4 This is a schematic diagram of the main structural region obtained using the method provided in the embodiments of this application, as shown below. Figure 4 As shown, the white area represents the main structural region. By utilizing Gaussian convolution, a smooth spatial diffusion from the connected region to its neighborhood is achieved, which can quantify the proximity between the discrete region and the main structure, thereby enhancing the completeness and spatial continuity of the extracted main structural region.

[0037] In some embodiments, the strength of the connection between discrete regions and connected regions can be adjusted by adjusting the standard deviation of the Gaussian kernel: a larger standard deviation of the Gaussian kernel will cause the area of ​​the discrete region to expand, while a smaller standard deviation of the Gaussian kernel will cause the discrete region to become isolated pixel noise and unable to connect with the connected region.

[0038] In some embodiments, the threshold for determining the existence of a connected region with a principal structure is determined in the following manner: ;in, This represents the threshold used to determine the presence of a connected region with a principal structure. Quantiles representing the area values ​​of multiple connected regions; This represents a minimum area constraint. It avoids the problem of being too large or too small that may be caused by a fixed threshold, while eliminating isolated noise points with small areas and preserving continuous, geologically significant main structures.

[0039] In some embodiments, step S41 may further include the following steps: S411: determining the effective pixels in the response amplitude distribution determined in step S30; S412: determining the cumulative distribution function of the effective pixels; S413: determining the higher quantile of the cumulative distribution function and determining a predetermined threshold. The higher quantile, as a dynamic threshold, can adaptively separate the target region with the strongest response among different images, suppress background noise, preserve the subject structure, and reduce the complexity of subsequent processing.

[0040] In some embodiments, step S60 further includes the following steps: S61: Processing the main geological structural lines to expand their coverage area and obtain a buffer zone; S62: Removing pixel lines outside the buffer zone; S63: Removing the main geological structural lines determined in step S50 from the buffer zone; S64: Repeating steps S20-S50 for the pixel lines within the buffer zone to determine secondary geological structural lines. By expanding the main geological structural lines to form a buffer zone and excluding pixel lines outside the zone, the analysis scope can be focused on the area near the main structural lines, effectively eliminating noise or irrelevant structural information far from the main structures; removing the identified main geological structural lines themselves from the buffer zone ensures that the remaining pixel lines are mainly features that are spatially adjacent to the main structural lines but are not the main structures, achieving hierarchical extraction of primary and secondary structures, avoiding secondary structures being masked or confused by the signals of the main structures in the global analysis, while ensuring that the identification criteria of secondary structures are consistent with those of the main structures, thus improving the overall consistency of geological structural line network extraction.

[0041] In some embodiments, in step S20, invalid values ​​in the stretched surface elevation data can also be removed to ensure that subsequent processing can focus on the effective imaging area.

[0042] In some embodiments, such as Figure 1 As shown, in step S70, after determining the main geological structural lines and secondary geological structural lines, these lines can be smoothed to obtain... Figure 1 Smooth geological structural lines in the middle.

[0043] In some embodiments, in step S70, the main geological structural lines and secondary geological structural lines can be smoothed in the following manner: S71: Determine the spatial trajectory of the main geological structural lines and secondary geological structural lines based on the pixels of the main geological structural lines and secondary geological structural lines; S72: Remove redundant nodes in the spatial trajectory; S73: Determine the smoothed structural lines based on the spatial trajectory processed in step S72; S74: Determine the spatial distribution of the smoothed structural lines in the geographic coordinate system based on the surface elevation data in step S10; S75: Determine the minimum structural line length threshold; S76: Remove small and isolated structural lines from the smoothed structural lines based on the minimum structural line length threshold. By removing redundant nodes and smoothing the main and secondary geological structural lines, the jagged turns and redundant inflection points caused by raster discretization can be eliminated, making the structural line shape closer to the natural curvature of the real geological structural lines. The smoothed structural lines are accurately mapped to the geographic coordinate system based on the surface elevation data, and small isolated line segments are removed according to the minimum length threshold. This avoids the interference of short and noisy line segments on the overall structural distribution and improves the macroscopic accuracy of geological structural line extraction.

[0044] Another embodiment of this application provides a system for determining geological structural lines in an area to be explored, comprising: a sensor configured to determine surface elevation data of the area to be explored; and a processor configured to: determine, based on the surface elevation data, the response amplitude of the terrain at each location in the area to be explored to geological structural lines in different predetermined directions; determine, based on the response amplitude, the response amplitude distribution corresponding to the dominant direction of the terrain at each location; determine, based on the response amplitude distribution, the main structural region of the area to be explored; determine, based on the main structural region, the main geological structural line of the area to be explored; determine, based on the main geological structural line, the secondary geological structural line of the area to be explored; and determine, based on the main geological structural line and the secondary geological structural line, the geological structural line of the area to be explored.

[0045] The system for determining geological structural lines in an area to be explored, provided in the embodiments of this application, extracts the response amplitude of the terrain to geological structures in different directions from surface elevation data. This allows for the determination of the dominant topographic direction and its corresponding response amplitude at each location, thereby accurately identifying the main structural region and its primary geological structural lines. Secondary geological structural lines are then extracted based on this, achieving a complete network of geological structural lines. This ensures that the extraction results accurately reflect the distribution of geological structural lines throughout the entire area to be explored, improving the objectivity and resolution of structural line extraction and avoiding confusion between primary and secondary structural lines. Compared to the subjectivity and limitations of traditional manual extraction, the method provided in the embodiments of this application can eliminate topographic noise and highlight linear topographic features consistent with structural trends, enabling the objective identification of geological structural lines at different scales and levels, thus improving the accuracy, efficiency, and repeatability of geological structural line extraction.

[0046] Embodiments of this application also provide a computer storage medium storing a computer program, the computer program including program instructions, which, when executed by a processor, perform the aforementioned method.

[0047] Regarding the embodiments of this application, it should also be noted that, without conflict, the embodiments of this application and the features in the embodiments can be combined with each other to obtain new embodiments.

[0048] The above are merely specific embodiments of this application, but the scope of protection of this application is not limited thereto. The scope of protection of this application shall be determined by the scope of the claims.

Claims

1. A method for determining geological structural lines in an area to be explored, characterized in that, It includes the following steps: S10: Determine the surface elevation data of the area to be explored; S20: Based on the surface elevation data, determine the response amplitude of the terrain at each location in the area to be explored to geological structural lines in different predetermined directions; S30: Based on the response amplitude determined in step S20, determine the response amplitude distribution corresponding to the dominant terrain direction at each location; S40: Based on the response amplitude distribution determined in step S30, determine the main structural region of the area to be explored; S50: Based on the main structural region, determine the main geological structural line of the area to be explored; S60: Based on the main geological structural line determined in step S50, determine the secondary geological structural line of the area to be explored; S70: Determine the geological structure line of the area to be explored based on the main geological structure line and the secondary geological structure line; Step S60 also includes the following steps: S61: Process the main geological structural line to expand its coverage area and obtain a buffer zone; S62: Remove pixel lines outside the buffer area; S63: Remove the main geological structural line determined in step S50 from the buffer zone; S64: Repeat steps S20-S50 for the pixel lines in the buffer to determine the secondary geological structure lines.

2. The method according to claim 1, characterized in that, Step S50 also includes the following steps: S51: Determine the pixel lines of the main body structure; S52: Determine the local connectivity of the pixels in the pixel line; S53: Determine the weights of the local connectivity relationships; S54: Determine the endpoints of the main geological structural line based on the weights of the local connectivity relationships; S55: Determine the main geological structural line based on the endpoints.

3. The method according to claim 2, characterized in that, In step S53, the weights of the local connectivity relationships are determined as follows: ； in, and This indicates two pixels that are adjacent to each other. and These represent the spatial coordinates of two adjacent pixels; The weight represents the local connectivity relationship.

4. The method according to claim 2, characterized in that, In step S54, the endpoints of the main geological structural line are determined as follows: ； in, and Indicates the endpoints of the main geological structural line; represents the starting pixel; x represents any pixel; V represents the set of all pixels. This represents the sum obtained by accumulating the weights along the local connectivity relationship from the starting pixel to any pixel. This represents the sum obtained by accumulating the weights along the local connectivity relationship from the endpoint to any pixel. This represents the pixel that maximizes the sum obtained by accumulating the weights.

5. The method according to claim 1, characterized in that, In step S20, the response amplitude of the terrain at each location to geological structural lines in different predetermined directions is determined as follows: ； in, These are the angles of the different predetermined directions; This represents the spatial pixel coordinates of each location in the surface elevation data; This indicates the terrain at each location for different predetermined directions. The response amplitude of geological structural lines; and These are the real and imaginary parts of the response amplitude, respectively.

6. The method according to claim 1, characterized in that, In step S30, the response amplitude distribution corresponding to the dominant terrain direction at each location is determined as follows: ； in, Indicates the angles of the different predetermined directions; This represents the spatial pixel coordinates of each location in the surface elevation data; This indicates the terrain at each location for different predetermined directions. The response amplitude of geological structural lines; Indicates that all The response amplitude in the direction Operations on the maximum value; This represents the set of maximum response amplitudes for each location in all predetermined directions, and represents the response amplitude distribution corresponding to the terrain-dominant direction for each location.

7. The method according to claim 1, characterized in that, Step S40 also includes the following steps: S41: Determine a predetermined threshold based on the response amplitude distribution determined in step S30; S42: Remove the background texture of the response amplitude distribution according to the predetermined threshold to obtain a binary constructed map; S43: Based on the binary construction graph determined in step S42, determine the connected regions where the main structure exists; S44: Based on the connected region, determine the discrete region near the connected region; S45: Determine the main structural region of the area to be explored based on the connected region and the discrete region.

8. The method according to claim 7, characterized in that, In step S43, the existence of the connected regions of the main structure is determined in the following manner: ； in, Indicates the first A connected region exists within the main structure; Indicates the first The area of ​​each connected region; B represents the union symbol; B represents the set of connected regions that exist within the main structure. This represents the threshold used to determine the existence of a connected region in the main structure.

9. The method according to claim 7, characterized in that, In step S44, the discrete regions near the connected region are determined in the following manner: ； in, A binary matrix representing the existence of a connected region constructed by the main body; Represents the spatial coordinates of the binarized matrix; It is the standard deviation of the two-dimensional Gaussian convolution kernel; It is a two-dimensional discrete space convolution operator; It is the coordinates after convolution processing. The output value at that location.

10. The method according to claim 8, characterized in that, The threshold used to determine the existence of a connected region of the main structure is determined in the following manner: ； in, This represents the threshold used to determine the existence of a connected region of the main structure; Quantiles representing the area values ​​of multiple connected regions; This represents the minimum area constraint.

11. A system for determining geological structural lines in an area to be explored, applicable to performing the method as described in any one of claims 1-10, characterized in that, It includes: Sensors configured to determine surface elevation data of the area to be surveyed; The processor is configured to: Based on the surface elevation data, determine the response amplitude of the topography at each location in the area to be explored to geological structural lines in different predetermined directions; Based on the response amplitude, determine the response amplitude distribution corresponding to the dominant topographic direction at each location; based on the response amplitude distribution, determine the main structural region of the area to be explored; Based on the main structural region, determine the main geological structural line of the area to be explored; based on the main geological structural line, determine the secondary geological structural line of the area to be explored; based on the main geological structural line and the secondary geological structural line, determine the geological structural line of the area to be explored.

12. A computer storage medium, characterized in that, The computer storage medium stores a computer program, the computer program including program instructions, which, when executed by a processor, perform the method as described in any one of claims 1-10.

Images