Geological map automatic generation method based on geological survey route observation points

Abstract

The invention provides a geological map automatic generation method based on geological survey route observation points, and relates to the technical field of geological map generation, and the method comprises the steps: collecting observation point coordinate data of each geological survey route in a research area; based on the coordinate data of the observation points, the geological exploration route is divided into a main route and an auxiliary route; constructing a topological relation through an attribute priority matching rule and a suboptimal distance matching rule; carrying out pinout processing on a missing region in the topological relation; drawing a preliminary geological plane graph according to the topological relation and a sharp-out processing result; determining a new route perpendicular to each geological survey route for generating a second geological plan; superposing the preliminary geological plan and the second geological plan to obtain a comprehensive geological map; and performing image morphological processing on the comprehensive geological map to obtain a high-quality geological planar graph. According to the technical scheme, the whole process from observation point data to the high-quality geological map is automatically generated.

Description

Technical Field

[0001] The present application relates to the technical field of geological map generation, and in particular to a method for automatically generating a geological map based on geological survey route observation points. Background Art

[0002] Geological maps are basic tools for studying regional geological structures and are widely used in mineral resource exploration, engineering geological research, disaster prevention and control, and other fields. The traditional method of drawing geological maps usually relies on manual operations by geologists, generating geological plane maps or profiles by empirically connecting geological observation points. This method is time-consuming and inefficient, and the results drawn by different people vary greatly, lacking consistency and reliability, especially in areas with sparse data or complex geological conditions (such as faults and folds).

[0003] Existing technologies have partially realized the automated generation of geological maps, mainly based on 3D modeling methods of borehole data or the generation of profiles. However, there is no complete and efficient solution for the automatic generation of 2D planar geological maps of observation points along geological survey routes. In addition, the ability to handle sparse data and complex geological conditions needs to be improved. Therefore, it is of great practical significance to develop a method for automatically generating geological maps, especially an efficient method that combines the topological relationship of the main / secondary routes and the pinch-out processing technology. Summary of the invention

[0004] The purpose of the present invention is to provide a method for automatically generating a geological map based on geological survey route observation points in order to solve the problems of low efficiency, discontinuous geological interfaces and poor accuracy in sparse data areas in existing geological map generation methods.

[0005] The above-mentioned purpose of the present application is achieved through the following technical solutions:

[0006] S1: Collect the coordinate data of observation points of each geological survey route in the study area and divide the study area into various sub-areas;

[0007] S2: Based on the observation point coordinate data, the geological survey route is divided into a main route and a secondary route;

[0008] S3: Construct the topological relationship between the main route and the secondary route through the attribute priority matching rule and the suboptimal distance matching rule;

[0009] S4: performing a pinch-out process on the missing area in the topological relationship to obtain a pinch-out process result;

[0010] S5: Draw a preliminary geological plane map based on the topological relationship and pinch-out processing results;

[0011] S6: determining a new route perpendicular to each geological survey route; generating a second geological plane map through the new route in combination with the processing of steps S1 to S5; superimposing the preliminary geological plane map and the second geological plane map to obtain a comprehensive geological map;

[0012] S7: Perform image morphological processing on the comprehensive geological map to obtain a high-quality geological plan map.

[0013] Optionally, step S1 includes:

[0014] Lay out geological survey routes according to the scope of the study area;

[0015] Collect observation point coordinate data of each observation point on the geological survey route;

[0016] The coordinates of the observation point coordinate data are expressed as Z(x,y,z), where x and y are the plane coordinates of the geological observation point, x represents longitude and y represents latitude, and z is depth information, which represents the vertical position of the observation point;

[0017] The area and observation point density calculation formula of each sub-area in the study area is:

[0018]

[0019] Where D i is the observation point density of sub-region i; N i is the number of observation points in the sub-region; A i is the area of ​​the sub-region.

[0020] Optionally, step S2 includes:

[0021] In the geological survey route, according to the observation point coordinate data, the number of observation points, spatial distribution characteristics and coverage of the geological survey route are determined, and the spatial position weight W of each geological survey route is calculated:

[0022] W=α·Density norm +β·Centrality norm +γ·Coverage norm

[0023] Where, Density norm Represents the normalized observation point density, Centrality norm Indicates the normalized value of the distance from the center of the route to the geometric center of the study area; Coverage norm Indicates the normalized ratio of the adjacent area covered by the route; α, β, and γ are adjustable weight coefficients, α reflects the density of data points covered by the route, β measures the proximity of the route to the core of the region, and γ describes the range of the adjacent geological area covered by the route;

[0024]

[0025] Where Density is the density of observation points; min is the minimum density of observation points; Density max The maximum density of observation points; d is the distance between the route center and the geometric center of the area; d max is the maximum distance;

[0026] Compare the spatial position weights W of two adjacent geological survey routes, and classify the geological survey route with a smaller W value as the main route, and the geological survey route with a larger W value as the secondary route.

[0027] Optionally, step S3 includes:

[0028] Set the geological attributes of each observation point; geological attributes include: rock mass characteristics, fold structure characteristics and stratigraphic characteristics;

[0029] By using the attribute priority matching rule and combining the geological attributes of each observation point, the attribute matching score A(P m ,P s ),as follows

[0030]

[0031] Where P m Indicates the main line; P s Indicates a secondary route;

[0032] The Euclidean distance D(P m ,P s ):

[0033]

[0034] In the formula, x m ,y m , z m are the coordinates of the observation points on the main route, representing longitude, latitude and elevation respectively; y s ,y s , z s are the coordinates of the observation point on the secondary route;

[0035] Comprehensive score S(P m ,P s ) is calculated as:

[0036] S(P m ,P s )=α×A(Pm ,P s )-β×D(P m ,P s )

[0037] In the formula, α and β are weight coefficients;

[0038] According to the comprehensive score S(P m ,P s ), connect the observation points of the main route and the secondary route, and build the topological relationship between the main route and the secondary route.

[0039] Optionally, step S4 includes:

[0040] The specific steps of pinch-out processing include: adjusting the position of the missing connection point in the topological relationship through the pinch-out direction and the pinch-out distance;

[0041] Pinching distance D 1 The calculation formula is:

[0042] D 1 =n·L

[0043] Where n is the scale factor, L is the average distance between observation points, and the scale factor n ranges from [1,3];

[0044] Determine the pinch-out direction based on the geological properties of the observation point.

[0045] Optionally, step S6 includes:

[0046] Calculate the spatial overlap M between the preliminary geological plane map and the second geological plane map:

[0047]

[0048] In the formula, S overlap The overlapping area of ​​two geological maps, S 1 +S 2 is the total area of ​​the preliminary geological plan and the second geological plan;

[0049] Through the overlap degree M, determine the area with significant difference in the process of superimposing the preliminary geological plane map and the second geological plane map;

[0050] The areas with significant differences were re-pinch-out to obtain the final comprehensive geological map.

[0051] Optionally, step S7 includes:

[0052] Image morphological processing includes: dilation operation, erosion operation and feature extraction;

[0053] The expansion operation formula is as follows:

[0054]

[0055] In the formula, f(x,y) and f(x ′ ,y ′ ) represent the pixel values ​​of the original image and the image after the dilation operation, respectively, (x ′ ,y ′ ) is the local position on the structural element K, f(x ′ ,y ′ ) is the value of the pixel in the new image after the dilation operation;

[0056] K represents the neighborhood range of the structural element in the pixel (x, y), and the structural element is selected as a matrix;

[0057] The corrosion operation formula is as follows:

[0058]

[0059] Where f(x,y) is the grayscale value of the original image at pixel (x,y); f ′ (x, y) is the gray value of the image at that location after corrosion;

[0060] Through the edge extraction algorithm, the key geological features of the comprehensive geological map are extracted to obtain a high-quality geological plane map. The formula is as follows:

[0061]

[0062] Where G is the gradient amplitude, which is used to represent the magnitude of the gradient of each pixel; G x ,G y They represent the gradient components of the comprehensive geological map in the x and y directions, respectively;

[0063]

[0064] In the formula, Indicates the grayscale change rate of the comprehensive geological map in the horizontal direction;

[0065]

[0066] In the formula, Represents the grayscale change rate of the comprehensive geological map in the vertical direction.

[0067] An electronic device includes a processor, a memory, a user interface and a network interface, wherein the memory is used to store instructions, the user interface and the network interface are used to communicate with other devices, and the processor is used to execute the instructions stored in the memory so that the electronic device executes a method for automatically generating a geological map based on observation points of a geological survey route.

[0068] A computer-readable storage medium stores instructions, and when the instructions are executed, a method for automatically generating a geological map based on geological survey route observation points is executed.

[0069] The beneficial effects of the technical solution provided by this application are:

[0070] Based on the observation point coordinate data, the main route and the secondary route are first divided, and the topological relationship of the main / secondary route is established through attribute priority matching and suboptimal distance matching rules; for areas that cannot be directly connected, the pinch-out processing is used to supplement the missing data and generate a preliminary geological plan map; by constructing a vertical route to verify and optimize the preliminary geological map, the generated result is further optimized in combination with image morphology technology, and finally a high-quality geological plan map is generated. By dividing the main / secondary routes, establishing a topological relationship, combining pinch-out processing, vertical route verification and image morphology optimization technology, the whole process from observation point data to high-quality geological maps is automatically generated. The present invention solves the problem of discontinuity of geological interfaces in data-sparse areas through dynamic topological relationship construction and pinch-out processing technology, and improves the reliability and accuracy of geological maps through vertical verification and image morphology optimization. The present invention realizes the automation of the whole process of geological map generation, is suitable for the fields of mineral resource exploration, geological disaster prevention and control, and has important application value and practical significance. BRIEF DESCRIPTION OF THE DRAWINGS

[0071] The present application will be further described below with reference to the accompanying drawings and embodiments, in which:

[0072] Figure 1 is a step diagram in an embodiment of the present application;

[0073] Figure 2 It is a regional division diagram of two survey routes in the embodiment of the present application;

[0074] Figure 3 is the topological relationship in the embodiment of the present application;

[0075] Figure 4 It is a pinch-out processing diagram in an embodiment of the present application;

[0076] Figure 5 It is a connection comparison diagram of two geological maps in the embodiment of the present application;

[0077] Figure 6 It is a connection comparison diagram of two geological maps in the embodiment of the present application;

[0078] Figure 7 It is a schematic diagram of the structure of an electronic device in an embodiment of the present application. DETAILED DESCRIPTION

[0079] In order to have a clearer understanding of the technical features, purposes and effects of the present application, the specific implementation methods of the present application are now described in detail with reference to the accompanying drawings.

[0080] An embodiment of the present application provides a method for automatically generating a geological map based on observation points along a geological survey route.

[0081] Please refer to Figure 1 , Figure 1 This is a step diagram of a method for automatically generating a geological map based on geological survey route observation points in an embodiment of the present application, including:

[0082] S1: Collect the coordinate data of observation points of each geological survey route in the study area and divide the study area into various sub-areas;

[0083] S2: Based on the observation point coordinate data, the geological survey route is divided into a main route and a secondary route;

[0084] S3: Construct the topological relationship between the main route and the secondary route through the attribute priority matching rule and the suboptimal distance matching rule;

[0085] S4: performing a pinch-out process on the missing area in the topological relationship to obtain a pinch-out process result;

[0086] As an example, according to the topological relationship of the main / secondary routes and the pinch-out processing results, the observation points are connected to generate a geological curve. Based on the geological curve, a preliminary geological plane map is drawn. The preliminary geological map shows the spatial distribution trend of the strata, provides a general framework of regional geological characteristics, and lays the foundation for subsequent verification and optimization.

[0087] S5: Draw a preliminary geological plane map based on the topological relationship and pinch-out processing results;

[0088] As an example, according to the topological relationship of the main / secondary routes and the pinch-out processing results, the observation points are connected to generate a geological curve. Based on the geological curve, a preliminary geological plane map is drawn to show the spatial distribution trend of the strata in the study area.

[0089] S6: determining a new route perpendicular to each geological survey route; generating a second geological plane map through the new route in combination with the processing of steps S1 to S5; superimposing the preliminary geological plane map and the second geological plane map to obtain a comprehensive geological map;

[0090] S7: Perform image morphological processing on the comprehensive geological map to obtain a high-quality geological plan map.

[0091] As an example, the comprehensive geological map is subjected to image morphological processing to eliminate noise and smooth the interface, extract key features and output a final high-quality geological plane map, such as Figure 6 shown.

[0092] Step S1 includes:

[0093] Lay out geological survey routes according to the scope of the study area;

[0094] Collect observation point coordinate data of each observation point on the geological survey route;

[0095] The coordinates of the observation point coordinate data are expressed as Z(x,y,z), where x and y are the plane coordinates of the geological observation point, x represents longitude and y represents latitude, and z is depth information, which represents the vertical position of the observation point;

[0096] As an embodiment, through a unified coordinate system and unit, these coordinates provide spatial basic data and geological attribute information such as lithology, fold structure, etc. for the generation of geological maps.

[0097] The area and observation point density calculation formula of each sub-area in the study area is:

[0098]

[0099] Where D i is the observation point density of sub-region i; N i is the number of observation points in the sub-region; A i is the area of ​​the sub-region.

[0100] As an example, the observation point density reflects the distribution of observation points in each sub-region. A higher density value indicates good data coverage, whereas a lower density value indicates that the observation point density of the region needs to be increased. The area of ​​the sub-region is calculated by the boundary polygon of the sub-region. The observation point density D of the sub-region i It is an important indicator for measuring the uniformity of the distribution of observation points. By dividing the area and laying out the routes, we ensure that the observation points cover the study area and provide basic data for subsequent analysis, such as Figure 2 shown.

[0101] Step S2 includes:

[0102] In the geological survey route, according to the observation point coordinate data, the number of observation points, spatial distribution characteristics and coverage of the geological survey route are determined, and the spatial position weight W of each geological survey route is calculated:

[0103] W=α·Density norm +β·Centrality norm +γ·Coverage norm

[0104] Where, Density normRepresents the normalized observation point density, Centrality norm Indicates the normalized value of the distance from the center of the route to the geometric center of the study area; Coverage norm Indicates the normalized ratio of the adjacent area covered by the route; α, β, and γ are adjustable weight coefficients, α reflects the density of data points covered by the route, β measures the proximity of the route to the core of the region, and γ describes the range of the adjacent geological area covered by the route;

[0105]

[0106] Where Density is the density of observation points; min is the minimum density of observation points; Density max The maximum density of observation points; d is the distance between the route center and the geometric center of the area; d max is the maximum distance;

[0107] Compare the spatial position weights W of two adjacent geological survey routes, and classify the geological survey route with a smaller W value as the main route, and the geological survey route with a larger W value as the secondary route.

[0108] As an embodiment, the main route is distributed in the main stratigraphic belt of the target area, and the number of observation points is small and distributed continuously; the sub-route covers the areas on both sides of the main route to provide supplementary information for establishing the topological relationship.

[0109] Step S3 includes:

[0110] Set the geological attributes of each observation point; geological attributes include: rock mass characteristics, fold structure characteristics and stratigraphic characteristics;

[0111] By using the attribute priority matching rule and combining the geological attributes of each observation point, the attribute matching score A(P m ,P s ),as follows

[0112]

[0113] Where P m Indicates the main line; P s Indicates a secondary route;

[0114] The Euclidean distance D(P m ,P s ):

[0115]

[0116] In the formula, x m,y m , z m are the coordinates of the observation points on the main route, representing longitude, latitude and elevation respectively; y s ,y s , z s are the coordinates of the observation point on the secondary route;

[0117] Comprehensive score S(P m ,P s ) is calculated as:

[0118] S(P m ,P S )=α×A(P m ,P S )-β×D(P m ,P S )

[0119] In the formula, α and β are weight coefficients;

[0120] According to the comprehensive score S(P m ,P s ), connect the observation points of the main route and the secondary route, and build the topological relationship between the main route and the secondary route.

[0121] As an example, when setting the attribute priority matching and the secondary distance matching rules, it is achieved by setting α to 1 and β to 0.1, and then matching the comprehensive score S(P m ,P S ) to connect the highest secondary route points and build the topological relationship between the main / secondary routes, such as Figure 3 shown.

[0122] Step S4 includes:

[0123] The specific steps of pinch-out processing include: adjusting the position of the missing connection point in the topological relationship through the pinch-out direction and the pinch-out distance;

[0124] Pinching distance D 1 The calculation formula is:

[0125] D 1 =n·L

[0126] Where n is the scale factor, L is the average distance between observation points, and the scale factor n ranges from [1,3];

[0127] Determine the pinch-out direction based on the geological properties of the observation point.

[0128] As an example, the pinch-out process, such as Figure 4As shown. The pinch-out direction depends on the fold structure or lithology change trend. In the continuous pinch-out section, the direction is kept consistent to ensure the logical coherence of the geological interface. For areas where the main / secondary routes cannot be directly connected, pinch-out processing technology is used to supplement them. According to the pinch-out direction, it is divided into forward pinch-out and backward pinch-out. Forward pinch-out is used to fill in the missing area in front of the main route observation point, and backward pinch-out is used to fill in the missing area behind the main route observation point.

[0129] Step S6 includes:

[0130] Calculate the spatial overlap M between the preliminary geological plane map and the second geological plane map:

[0131]

[0132] In the formula, S overlap The overlapping area of ​​two geological maps, S 1 +S 2 is the total area of ​​the preliminary geological plan and the second geological plan;

[0133] Through the overlap degree M, determine the area with significant difference in the process of superimposing the preliminary geological plane map and the second geological plane map;

[0134] The areas with significant differences were re-pinch-out to obtain the final comprehensive geological map.

[0135] As an example, Figure 5 As shown in the figure, the overlap M can quantify the consistency of the spatial distribution of the two geological maps. The higher the value, the smaller the difference between the two maps. The difference area is dynamically optimized. In the area with significant difference, the topological rules of the main / secondary route can be adjusted, or the pinch-out direction can be readjusted to ensure that the optimized comprehensive geological map is closer to the actual geological conditions.

[0136] As an embodiment, a new survey route perpendicular to the main route is constructed in the study area to capture geological changes in areas not covered by the main route and provide more comprehensive data. By collecting observation point data on the vertical route, the connection relationship of the difference area is optimized by adjusting the weight parameters to ensure the consistency of the two geological maps and generate a second geological map.

[0137] Step S7 includes:

[0138] Image morphological processing includes: dilation operation, erosion operation and feature extraction;

[0139] The expansion operation formula is as follows:

[0140]

[0141] In the formula, f(x,y) and f(x′ ,y ′ ) represent the pixel values ​​of the original image and the image after the dilation operation, respectively, (x ′ ,y ′ ) is the local position on the structural element K, f(x ′ ,y ′ ) is the value of the pixel in the new image after the dilation operation;

[0142] K represents the neighborhood range of the structural element in the pixel (x, y), and the structural element is selected as a matrix;

[0143] As an example, the dilation operation makes the highlight area larger, and the broken noise points in the geological map are connected or filled.

[0144] The corrosion operation formula is as follows:

[0145]

[0146] Where f(x,y) is the grayscale value of the original image at pixel (x,y); f ′ (x, y) is the gray value of the image at that location after corrosion;

[0147] As an example, the erosion operation reduces highlight areas, removes smaller noise points, and filters out unnecessary details in the geological map.

[0148] Through the edge extraction algorithm, the key geological features of the comprehensive geological map are extracted to obtain a high-quality geological plane map. The formula is as follows:

[0149]

[0150] Where G is the gradient amplitude, which is used to represent the magnitude of the gradient of each pixel; G x ,G y They represent the gradient components of the comprehensive geological map in the x and y directions, respectively;

[0151]

[0152] In the formula, Indicates the grayscale change rate of the comprehensive geological map in the horizontal direction;

[0153]

[0154] In the formula, Represents the grayscale change rate of the comprehensive geological map in the vertical direction.

[0155] As an example, key geological features include faults and boundaries of pinch-out horizons.

[0156] The present application also discloses an electronic device. Figure 7 , Figure 7 The electronic device 500 may include: at least one processor 501 , at least one network interface 504 , a user interface 503 , a memory 505 , and at least one communication bus 502 .

[0157] The communication bus 502 is used to realize the connection and communication between these components.

[0158] The user interface 503 may include a display screen, and the optional user interface 503 may also include a standard wired interface or a wireless interface.

[0159] The network interface 504 may optionally include a standard wired interface or a wireless interface (such as a WI-FI interface).

[0160] The present application also discloses a computer-readable storage medium storing a plurality of instructions suitable for loading by a processor to execute the above-mentioned method for automatically generating a geological map based on geological survey route observation points.

[0161] The above are only exemplary embodiments of the present disclosure and cannot be used to limit the scope of the present disclosure. That is, any equivalent changes and modifications made according to the teachings of the present disclosure are still within the scope of the present disclosure.

[0162] This application is intended to cover any variation, use or adaptation of the present disclosure, which follows the general principles of the present disclosure and includes common knowledge or customary technical means in the art not described in the present disclosure. The description and examples are to be regarded as exemplary only, and the scope and spirit of the present disclosure are defined by the claims.

Claims

1. A method for automatically generating a geological map based on observation points of a geological survey route, characterized in that: The method comprises the following steps: S1: Collect the coordinate data of observation points of each geological survey route in the study area and divide the study area into various sub-areas; S2: Based on the observation point coordinate data, the geological survey route is divided into a main route and a secondary route; S3: Construct the topological relationship between the main route and the secondary route through the attribute priority matching rule and the suboptimal distance matching rule; S4: performing a pinch-out process on the missing area in the topological relationship to obtain a pinch-out process result; S5: Draw a preliminary geological plane map based on the topological relationship and pinch-out processing results; S6: determining a new route perpendicular to each geological survey route; generating a second geological plane map through the new route in combination with the processing of steps S1 to S5; superimposing the preliminary geological plane map and the second geological plane map to obtain a comprehensive geological map; S7: Perform image morphological processing on the comprehensive geological map to obtain a high-quality geological plan map.

2. A method for automatically generating a geological map based on geological survey route observation points as claimed in claim 1, characterized in that: Step S1 includes: Lay out geological survey routes according to the scope of the study area; Collect observation point coordinate data of each observation point on the geological survey route; The coordinates of the observation point coordinate data are expressed as Z(x,y,z), where x and y are the plane coordinates of the geological observation point, x represents longitude and y represents latitude, and z is depth information, which represents the vertical position of the observation point; The area and observation point density calculation formula of each sub-area in the study area is: Where D i is the observation point density of sub-region i; N i is the number of observation points in the sub-region; A i is the area of ​​the sub-region.

3. A method for automatically generating a geological map based on geological survey route observation points as claimed in claim 1, characterized in that: Step S2 includes: In the geological survey route, according to the observation point coordinate data, the number of observation points, spatial distribution characteristics and coverage of the geological survey route are determined, and the spatial position weight W of each geological survey route is calculated: W=α·Density norm +β·Centrality norm +γ·Coverage norm Where, Density norm Represents the normalized observation point density, Centrality norm Indicates the normalized value of the distance from the center of the route to the geometric center of the study area; Coverage norm Indicates the normalized ratio of the adjacent area covered by the route; α, β, and γ are adjustable weight coefficients, α reflects the density of data points covered by the route, β measures the proximity of the route to the core of the region, and γ describes the range of the adjacent geological area covered by the route; Where Density is the density of observation points; min is the minimum density of observation points; Density max The maximum density of observation points; d is the distance between the route center and the geometric center of the area; d max is the maximum distance; Compare the spatial position weights W of two adjacent geological survey routes, and classify the geological survey route with a smaller W value as the main route, and the geological survey route with a larger W value as the secondary route.

4. A method for automatically generating a geological map based on geological survey route observation points as claimed in claim 1, characterized in that: Step S3 includes: Set the geological attributes of each observation point; geological attributes include: rock mass characteristics, fold structure characteristics and stratigraphic characteristics; By using the attribute priority matching rule and combining the geological attributes of each observation point, the attribute matching score A(P m ,P s ),as follows Where P m Indicates the main line; P s Indicates a secondary route; The Euclidean distance D(P m ,P s ): In the formula, x m ,y m , z m are the coordinates of the observation points on the main route, representing the longitude, latitude and elevation respectively; y s ,y s , z s are the coordinates of the observation point on the secondary route; Comprehensive score S(P m ,P s ) is calculated as: S(P m ,P s )=α×A(P m ,P s )-β×D(P m ,P s ) In the formula, α and β are weight coefficients; According to the comprehensive score S(P m ,P s ), connect the observation points of the main route and the secondary route, and build the topological relationship between the main route and the secondary route.

5. A method for automatically generating a geological map based on geological survey route observation points as claimed in claim 1, characterized in that: Step S4 includes: The specific steps of pinch-out processing include: adjusting the position of the missing connection point in the topological relationship through the pinch-out direction and the pinch-out distance; The calculation formula of the pinch-out distance D1 is: D1=n·L Where n is the scale factor, L is the average distance between observation points, and the scale factor n ranges from [1,3]; Determine the pinch-out direction based on the geological properties of the observation point.

6. A method for automatically generating a geological map based on geological survey route observation points as claimed in claim 1, characterized in that: Step S6 includes: Calculate the spatial overlap M between the preliminary geological plane map and the second geological plane map: In the formula, S overlap The overlapping area of ​​the two geological maps, S1+S2 is the total area of ​​the preliminary geological plan map and the second geological plan map; Through the overlap degree M, determine the area with significant difference in the process of superimposing the preliminary geological plane map and the second geological plane map; The areas with significant differences were re-pinch-out to obtain the final comprehensive geological map.

7. A method for automatically generating a geological map based on geological survey route observation points as claimed in claim 1, characterized in that: Step S7 includes: Image morphological processing includes: dilation operation, erosion operation and feature extraction; The expansion operation formula is as follows: In the formula, f(x,y) and f(x ′ ,y ′ ) represent the pixel values ​​of the original image and the image after the dilation operation, respectively, (x ′ ,y ′ ) is the local position on the structural element K, f(x ′ ,y ′ ) is the value of the pixel in the new image after the dilation operation; K represents the neighborhood range of the structural element in the pixel (x, y), and the structural element is selected as a matrix; The corrosion operation formula is as follows: Where f(x,y) is the grayscale value of the original image at pixel (x,y); f ′ (x, y) is the gray value of the image at that location after corrosion; Through the edge extraction algorithm, the key geological features of the comprehensive geological map are extracted to obtain a high-quality geological plane map. The formula is as follows: Where G is the gradient amplitude, which is used to represent the magnitude of the gradient of each pixel; G x ,G y They represent the gradient components of the comprehensive geological map in the x and y directions, respectively; In the formula, Indicates the grayscale change rate of the comprehensive geological map in the horizontal direction; In the formula, Represents the grayscale change rate of the comprehensive geological map in the vertical direction.

8. An electronic device, characterized in that: It includes a processor, a memory, a user interface and a network interface, the memory is used to store instructions, the user interface and the network interface are used to communicate with other devices, and the processor is used to execute the instructions stored in the memory so that the electronic device executes the method as described in any one of claims 1-7.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores instructions, and when the instructions are executed by a computer, the method according to any one of claims 1 to 7 is executed.

Images