Method for automatically obtaining geological point properties and geological point lithology according to geological regions

Abstract

The application discloses a kind of according to geological area automatically obtaining geological point property and geological point lithology, it is related to geological exploration information technology field, comprising: in GIS based on project number and route number, from geological area and geological point layer, the elements associated with the geological point to be analyzed are extracted, set buffer radius to carry out analysis, generate geological buffer zone, and the intersection of geological area, geological point and buffer zone is obtained by intersection analysis;According to the number of intersection elements, judge whether the topological error and lithology exist in geological map, if there is no topological error, directly determine the lithology of geological point, if there is error, calculate and retain the two area elements with the largest area outside buffer zone, convert them into point elements and inherit all attributes, calculate the point type, radian, degree, eight azimuth of the geological point, and obtain the lithology analysis result;The application solves the efficiency bottleneck, precision defect and quality control dilemma existing in traditional method, and provides a replicable, generalizable automatic solution for geological survey work.

Description

Technical Field

[0001] This invention relates to the field of information technology for geological exploration, and in particular to a method for automatically obtaining the properties and lithology of geological points based on geological regions. Background Technology

[0002] In the field of geological surveys, the determination of the nature of geological points (lithological control points / boundary points) and lithology is the foundation for constructing geological maps. Traditional methods rely on manual operation: geologists need to locate the GPS coordinates collected in the field on the geological map and visually determine the location of the point. If the point is located entirely within a single geological area (such as a granite area), it is marked as a "lithological control point"; if it is close to a geological boundary (such as a fault line), the distance is measured using a scale bar—those with a distance ≤1mm are marked as "boundary points," otherwise they are classified as control points. Finally, the lithological code and description are manually copied from the map legend manual.

[0003] However, this process has systemic flaws: First, there's the efficiency bottleneck; processing a single point takes 3-5 minutes, including 1.2 minutes for coordinate transformation, 1.5 minutes for boundary measurement, and 0.8 minutes for lithology lookup. Second, there's the accuracy issue; paper map distortion leads to field offsets, resulting in high rates of boundary point misjudgment, lithology copying errors, and topological omissions. Third, there's the quality control dilemma; manual records cannot be traced, and early errors are easily amplified during map production. Therefore, the properties and lithology of geological points obtained through manual operation suffer from low reliability and poor consistency, leading to limited accuracy in geological information and hindering geological surveys. Summary of the Invention

[0004] This invention provides a method for automatically obtaining the properties and lithology of geological points based on geological regions, thereby overcoming the aforementioned technical problems.

[0005] To achieve the above objectives, the technical solution of the present invention is as follows:

[0006] A method for automatically obtaining the properties and lithology of geological points based on a geological region includes:

[0007] S1: Obtain topologically verified geological region feature layer data and geological point feature layer data in GIS. The geological region feature layer data includes geological body code, lithology, lithology description and project number field where the geological region is located. The geological point feature layer data includes point number of geological point and route number field of geological point route.

[0008] S2: Based on the project number field and route number field, use the GIS attribute selection tool to extract the geological area elements and geological point elements associated with the geological point to be analyzed from the geological area element layer data and geological point element layer data;

[0009] S3: Set the buffer radius, perform GIS buffer analysis on geological point features within the buffer radius, and generate a geological buffer.

[0010] S4: Analyze the geological region features and geological point features in S2 with the geological buffer zone in S3 using a GIS intersection analysis tool to obtain the intersection of the geological region features and geological point features with the geological buffer zone, thus obtaining the number of intersecting features. Based on the number of intersecting features, determine whether the lithology of the geological point to be analyzed and whether the geological map where the geological point is located has a topological error. If the number of intersecting features is less than or equal to 2, the geological map has no topological error; if the number of intersecting features is greater than 2, the geological map has a topological error. Integrate all attribute fields of the geological region features and geological point features in the intersection to obtain a region feature class file.

[0011] S5: If the judgment result is that there is no topological error on the geological map, the lithology of the geological point to be analyzed is directly judged. If the judgment result is that there is a topological error on the geological map, the area of ​​each zone element located outside the geological buffer zone is calculated, and the areas are sorted from largest to smallest. Only the zone elements with the top two areas are retained. The zone elements are the geological zone elements located outside the geological buffer zone generated after the intersection analysis of the geological zone elements and geological point elements with the geological buffer zone.

[0012] S6: Use the GIS feature to point conversion tool to convert the retained area features into point features, and retain all attribute fields integrated in the area feature class file;

[0013] S7: Based on the conversion results, use GIS to calculate the point type of the geological point to be analyzed, and use the geological route as the baseline to calculate the relative azimuth (radians), relative azimuth (degrees), and relative azimuth (octagonal) between the geological point to be analyzed and the baseline, so as to obtain the corresponding geological point properties and geological point lithology analysis results.

[0014] Furthermore, the radius R of the buffer zone is equal to a preset value multiplied by the denominator of the mapping scale.

[0015] Furthermore, if the determination result indicates a topological error in the geological map, the area of ​​each feature located outside the geological buffer zone is calculated, and the areas are sorted from largest to smallest. Only the features in the top two areas are retained, including:

[0016] S51. Group the area elements with the same geological point number in the intersection results according to the geological point number, and calculate the area of ​​the area elements in each group.

[0017] S52. For groups with ≥3 elements under the same geological point number, sort the area elements in descending order by area.

[0018] S53. Elements with a sorted sequence number ≥ 3 are designated as deletion objects;

[0019] S54. Perform the GIS feature deletion operation to delete the marked area features.

[0020] Furthermore, if the judgment result indicates that there are no topological errors on the geological map, the lithology of the geological points is directly determined, specifically including:

[0021] When the intersection result shows that there is only one element in the geological buffer zone, the geological point is a lithological control point;

[0022] When the intersection result shows that there are two elements in the geological buffer zone, the geological point is the lithological boundary point.

[0023] Furthermore, based on the conversion results, GIS is used to calculate the point type, relative azimuth in radians, relative azimuth in degrees, and relative azimuth in octahedron of the geological points to be analyzed, obtaining the corresponding geological point properties and lithological analysis results, including:

[0024] S71. Construct an SQL query expression based on the geological point route numbers;

[0025] S72. Use the SearchCursor cursor to read the transformation results based on the SQL query expression and save the read data; the read data includes the OID value of the geological points in the point features, the X and Y coordinates of the geometric center point in the point features, and the geological point number field in the point features;

[0026] S73. First traversal of the conversion results: Traverse the conversion results, create a point_dict dictionary and a p_num list, store the geological point numbers in the point_dict dictionary, and store the corresponding OID list in the p_num list to form a number index;

[0027] S74. Perform a second traversal of the conversion results: Use the UpdateCursor cursor to read all data from the point_dict dictionary, traverse each geological point, and determine the number of OIDs stored in the p_num list. When the number of OIDs is 1, the geological point is a geological lithology control point; when the number of OIDs is 2, the geological point is a geological lithology boundary point. The geological lithology boundary point includes the start and end points of a boundary line; thus, the lithology of the geological point is obtained.

[0028] S75. When the number of OIDs is 2, find the end point / start point paired with the current start point / end point by using the number index according to the current point number and the value of the OID in the p_num list. Subtract the X coordinate of the paired end point / start point from the X coordinate of the current start point / end point to obtain dx. Subtract the Y coordinate of the paired end point / start point from the Y coordinate of the current start point / end point to obtain dy. Call the math.atan2(dx, dy) function to obtain the relative radian of the two points. Call the math.degrees method to convert the radian to an angle value. If the angle value is positive, save the angle value directly. If the angle value is negative, add 360 degrees to obtain the final angle value and save it.

[0029] S76. Calculate the angle values ​​using GIS to transform them into eight directions: North, Northeast, East, Southeast, South, Southwest, West, and Northwest, thereby obtaining the properties of the geological points.

[0030] Beneficial effects: This invention provides a method for automatically obtaining the properties and lithology of geological points based on geological regions, which has the following advantages:

[0031] 1. The precision and accuracy of this invention are greatly improved: Through GIS spatial analysis, human error can be eliminated, and the logical consistency of the relationship between geological boundaries and geological bodies can be ensured by using topologically validated data sources;

[0032] 2. Efficiency Improvement: This invention uses GIS to process data, supporting the simultaneous processing of hundreds of geological points, greatly improving the efficiency of regional geological surveys. The processing efficiency is 10-20 times higher than traditional methods, and it is especially suitable for rapid analysis of a large number of geological points, significantly shortening the work cycle.

[0033] 3. Reduced resource and cost usage: The entire process from data extraction and buffer analysis to lithology determination is automated, reducing the need for manual intervention by more than 90% and significantly lowering labor costs; redundant calculations are reduced through algorithm optimization, enabling efficient operation on ordinary computing devices and wider applicability;

[0034] This invention solves the efficiency bottlenecks, accuracy deficiencies, and quality control difficulties of traditional methods, and provides a replicable and scalable automated solution for geological survey work, which has important practical application value and technical demonstration significance. Attached Figure Description

[0035] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0036] Figure 1 This invention provides a flowchart of a method for automatically obtaining the properties and lithology of geological points based on a geological region;

[0037] Figure 2 This is a schematic diagram of the geological buffer zone in the GIS of this invention;

[0038] Figure 3 This is a schematic diagram of the intersection results in the GIS of this invention. Detailed Implementation

[0039] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0040] This embodiment provides a method for automatically obtaining the properties and lithology of geological points based on geological regions, such as... Figure 1 As shown, it includes:

[0041] S1: Obtain topologically verified geological region feature layer data and geological point feature layer data in GIS. The geological region feature layer data includes geological body code, lithology, lithology description and project number field where the geological region is located. The geological point feature layer data includes point number of geological point and route number field of geological point route.

[0042] S2: Based on the project number field and route number field, use the GIS attribute selection tool to extract the geological area elements and geological point elements associated with the geological point to be analyzed from the geological area element layer data and geological point element layer data;

[0043] S3: Set the buffer radius, perform GIS buffer analysis on geological point features within the buffer radius, and generate a geological buffer.

[0044] S4: Analyze the geological region features and geological point features in S2 with the geological buffer zone in S3 using a GIS intersection analysis tool to obtain the intersection of the geological region features and geological point features with the geological buffer zone, thus obtaining the number of intersecting features. Based on the number of intersecting features, determine whether the lithology of the geological point to be analyzed and whether the geological map where the geological point is located has a topological error. If the number of intersecting features is less than or equal to 2, the geological map has no topological error; if the number of intersecting features is greater than 2, the geological map has a topological error. Integrate all attribute fields of the geological region features and geological point features in the intersection to obtain a region feature class file.

[0045] S5: If the judgment result is that there is no topological error on the geological map, the lithology of the geological point to be analyzed is directly judged. If the judgment result is that there is a topological error on the geological map, the area of ​​each zone element located outside the geological buffer zone is calculated, and the areas are sorted from largest to smallest. Only the zone elements with the top two areas are retained. The zone elements are the geological zone elements located outside the geological buffer zone generated after the intersection analysis of the geological zone elements and geological point elements with the geological buffer zone.

[0046] S6: Use the GIS feature to point conversion tool to convert the retained area features into point features, and retain all attribute fields integrated in the area feature class file;

[0047] S7: Based on the conversion results, use GIS to calculate the point type of the geological point to be analyzed, and use the geological route as the baseline to calculate the relative azimuth (radians), relative azimuth (degrees), and relative azimuth (octagonal) between the geological point to be analyzed and the baseline, so as to obtain the corresponding geological point properties and geological point lithology analysis results.

[0048] Specifically, firstly, topologically validated geological region feature layer data and geological point feature layer data are obtained from GIS. The geological region feature layer data includes fields for geological body code, lithology, lithology description, and project number of the geological region. The geological point feature layer data includes field for point number of the geological point and field for route number of the geological point route. Based on the project number field and route number field, GIS attribute selection tools are used to extract geological region features and geological point features associated with the geological point to be analyzed from the geological region feature layer data and geological point feature layer data. The logical consistency of the relationship between geological boundaries and geological bodies is ensured by using a topologically validated data source in GIS.

[0049] Secondly, a buffer radius is set, and GIS buffer analysis is performed on the geological point features within this radius to generate a geological buffer. The geological area features and geological point features in S2 are then analyzed with the geological buffer in S3 using a GIS intersection analysis tool to obtain the intersection of the geological area features and geological point features with the geological buffer, thus determining the number of intersecting features. All attribute fields of the geological area features and geological point features in the intersection are then integrated to obtain a region feature class file. Based on the number of intersecting features, the lithology of the geological point to be analyzed and whether there are topological errors on the geological map where the geological point is located are determined. If the number of intersecting features is less than... If the number of intersecting elements is greater than 2, then the geological map has no topological errors. If the number of intersecting elements is greater than 2, then the geological map has topological errors. All attribute fields of geological area elements and geological point elements in the intersection are integrated to obtain a regional element class file. By generating a geological buffer, the subjective and tedious manual operation of "measuring distances ≤1mm using a scale bar" in traditional methods can be replaced, eliminating subjective judgment standards and measurement errors. Intersection analysis using GIS replaces the subjective process of "visually judging point location," providing objective and quantitative judgment standards and realizing the detection of topological errors, avoiding the possibility of errors being amplified in a chain reaction.

[0050] If the geological map shows no topological errors, the lithology of the geological point to be analyzed is directly determined. If the geological map shows topological errors, the area of ​​each feature outside the geological buffer zone is calculated, and the areas are sorted from largest to smallest, retaining only the top two features. These features are those outside the geological buffer zone, generated after intersection analysis of the geological zone features, geological point features, and the geological buffer zone. By calculating and comparing the areas outside the buffer zone, the two most likely related geological bodies can be selected, improving processing efficiency and avoiding the confusion of dealing with multiple complex geological bodies during manual judgment. This provides a clear and quantifiable target for subsequent orientation calculations.

[0051] Finally, the retained area features are converted into point features using a GIS feature conversion tool, while preserving all attribute fields integrated into the area feature class file. Based on the conversion results, GIS is used to calculate the point type of the geological points to be analyzed, and the relative azimuth (radians), relative azimuth (degrees), and relative octahedron between the geological points and the baseline are calculated using the geological route as a baseline, thus obtaining the geological point properties and lithology. The automated conversion and calculation process using GIS ensures lossless inheritance of attributes and guarantees the standardization and accuracy of the results.

[0052] In a specific embodiment, the geological region feature layer data and geological point feature layer data, which have been topologically verified, are obtained in the GIS. The geological region feature layer data includes fields for geological body code, lithology, lithology description, and the project number of the geological region. The geological point feature layer data includes fields for the point number of the geological point and the route number of the geological point route.

[0053] 1. Prepare a geological body feature layer. The geological body features must be verified by GIS topology tools to be free of overlap / gaps. The attributes contain the geological body code, lithology, lithology description and project number fields, in shp or gdb format.

[0054] 2. Prepare a geological point feature layer. The attributes should contain the point number and route number fields. When a geological point is a boundary point, there can be a certain error with the boundary of the geological area, but the error distance cannot exceed 1mm on the map. The format should be shp or gdb.

[0055] Both geological body feature layers and geological point feature layers are essential data in the geological mapping process. They are verified using GIS topology tools to ensure the logical consistency between geological boundaries and the relationship between geological bodies.

[0056] In a specific embodiment, the scheme for extracting geological area elements and geological point elements associated with the geological point to be analyzed from the geological area element layer data and geological point element layer data using GIS attribute selection tools based on the project number field and route number field is as follows:

[0057] The user provides the project number and route number. The topology-verified geological feature layer is used to extract the geological area portion for this project via GIS based on the project number. Similarly, the topology-verified geological point feature layer is used to extract the geological point portion for this project via GIS based on the route number. Both are stored in a temporary database for later use.

[0058] In a specific embodiment, a buffer radius is set, and GIS buffer analysis is performed on geological point features within the buffer radius to generate the geological buffer.

[0059] like Figure 2 As shown, a buffer zone is generated in the GIS with the geological point as the center, and the buffer zone radius is set to 1mm. The denominator of the mapping scale is set as follows: for example, when the mapping scale is 1:10000, the buffer distance is set to 10m, and the results are stored in a temporary database for later use.

[0060] By generating geological buffer zones, the traditional method of "measuring distances ≤1mm using a scale" can be replaced by subjective and tedious manual operations, eliminating subjective judgment standards and measurement errors.

[0061] In a specific embodiment, the geological region elements and geological point elements in S2 are analyzed with the geological buffer zone in S3 using a GIS intersection analysis tool to obtain the intersection of the geological region elements and geological point elements with the geological buffer zone, thus determining the number of intersecting elements. All attribute fields of the geological region elements and geological point elements in the intersection are then integrated to obtain a region element class file. Based on the number of intersecting elements, it is determined whether the lithology of the geological point to be analyzed and the geological map containing the geological point have topological errors. If the number of intersecting elements is less than or equal to 2, the geological map has no topological errors; if the number of intersecting elements is greater than 2, the geological map has topological errors. The scheme for integrating all attribute fields of the geological region elements and geological point elements in the intersection to obtain the region element class file is as follows:

[0062] An intersection analysis was performed on the overlay of the geological region results in S2 and the geological point buffer zone results in S3, such as... Figure 3 As shown, the intersection of the geological region results in S2 and the geological point buffer in S3 is obtained, and all attribute fields of both are integrated. The final result is a region feature class. If there is only one or two features, the lithology of the geological point to be analyzed can be directly determined. If there are more than two features, the geological point is a lithological boundary point and the geological map has a topological error. The intersection result is stored in a temporary database for later use.

[0063] Intersection analysis using GIS replaces the subjective process of manually "visually judging point ownership," providing objective and quantitative judgment criteria and enabling the detection of topological errors, thus avoiding the possibility of errors being amplified in a chain reaction.

[0064] In a specific embodiment, if the judgment result indicates that the geological map has no errors, the lithology of the geological point to be analyzed is directly determined. If the judgment result indicates that the geological map has a topological error, the area of ​​each zone element located outside the geological buffer zone is calculated, and the areas are sorted from largest to smallest, retaining only the zone elements with the top two areas. The zone elements are the geological zone element scheme located outside the geological buffer zone generated after the intersection analysis of the geological zone elements, geological point elements, and the geological buffer zone.

[0065] When the intersection result shows that there is only one element in the geological buffer zone, the geological point is a lithological control point;

[0066] like Figure 3 As shown, when the intersection result indicates that there are two elements in the geological buffer zone, the geological point is the lithological boundary point;

[0067] like Figure 3As shown, when there are more than two elements in the intersection result, the geological point is a lithological boundary point and the geological map has a topological error. The intersection result is then stored in a temporary database for later use. The presence of more than two elements indicates that the buffer zone simultaneously covers different geological regions. The two geological regions with the largest area outside the buffer zone, i.e., the two elements with the largest area, need to be retained. The specific steps are as follows:

[0068] S51. Group the area elements with the same geological point number in the intersection results according to the geological point number, and calculate the area of ​​the area elements in each group.

[0069] S52. For groups with ≥3 elements under the same geological point number, sort the area elements in descending order by area.

[0070] S53. Elements with a sorted sequence number ≥ 3 are designated as deletion objects;

[0071] S54. Perform the GIS feature deletion operation to delete the marked area features.

[0072] Specifically, the intersection results are loaded into memory using Python. Each geological point number is used as a unit, with a maximum of two features retained for the same geological point number. If more than two features are retained, the two features with the largest areas are chosen. The core algorithm is as follows:

[0073] (1) Construct an SQL expression based on the geological route number provided by the user, for example: where_clause = r"line_id = 'abc'";

[0074] (2) Use the SearchCursor cursor to read all features in the geological point feature layer into memory. The SQL expression is set to the result in (1), the query fields are set to "OID@", "SHAPE@AREA" and the geological point number field, and a list of dictionary objects is used to store the query results. The dictionary object format is {pid:items}, where pid is the geological point id, and items is a tuple set with the content [(oid value 1, area value 1), (oid value 2, area value 2), ...].

[0075] (3) Traverse the list of dictionary objects, select objects with more than 2 items, and then sort them in reverse order according to the value of item[1]. Use list comprehension ([item[0]for item in sorted_items[2:]]) to get the list of OIDs after the first 2 in terms of area and store them in the oids_to_delete set.

[0076] (4) Concatenate the oids_to_delete collection into an IN type SQL clause, with data separated by commas, and finally use the arcpy.management.DeleteRows method to delete redundant polygon features.

[0077] In this scheme, by calculating and comparing the area outside the buffer zone, the two main geological bodies most likely related to the point to be analyzed can be screened out, improving processing efficiency and avoiding the helplessness of facing multiple complex geological bodies when making manual judgments. This provides a clear and quantifiable target for subsequent orientation calculations.

[0078] In a specific embodiment, the solution of using a GIS feature-to-point tool to convert the retained area features into point features while retaining all attribute fields integrated into the area feature class file is as follows:

[0079] In GIS, select the area features to retain, specify the name and save path of the output point features, save the attribute fields contained in all area features, generate point features using the GIS feature-to-point tool, and copy all attribute fields to the corresponding generated points.

[0080] This solution uses GIS to automate the conversion process, ensuring lossless inheritance of attributes and standardization of results.

[0081] In a specific embodiment, based on the conversion results, GIS is used to calculate the point type of the geological point to be analyzed, and the relative azimuth (radians), relative azimuth (degrees), and relative azimuth (octahedron) between the geological point to be analyzed and the baseline are calculated using the geological route as the baseline, to obtain the corresponding geological point properties and lithological analysis results, including:

[0082] S71. Construct an SQL query expression based on the geological point route number; for example: where_clause = r"line_id = 'abc'";

[0083] S72. Use the SearchCursor cursor to read the transformation results based on the SQL query expression and save the read data; the read data includes the OID value of the geological points in the point features, the X and Y coordinates of the geometric center point in the point features, and the geological point number field in the point features;

[0084] Specifically, the SearchCursor cursor is used to read the transformation results and store the data in memory; the SQL expression is the result in S71, the query fields are set to "OID@", "SHAPE@XY" and the geological point number field, and a list of dictionary objects is used to store the query results;

[0085] S73. First traversal of the conversion results: Traverse the conversion results, create a point_dict dictionary and a p_num list, store the geological point numbers in the point_dict dictionary, and store the corresponding OID list in the p_num list to form a number index;

[0086] Specifically, a point number index is established, a point_dict object is created, a p_num list is created during traversal, and then the oid of the geological points is written into the list to form a number index, for example: {"D001":[1001,1002]};

[0087] S74. Perform a second traversal of the conversion results: Use the UpdateCursor cursor to read all data from the point_dict dictionary, traverse each geological point, and determine the number of OIDs stored in the p_num list. When the number of OIDs is 1, the geological point is a geological lithology control point; when the number of OIDs is 2, the geological point is a geological lithology boundary point. The geological lithology boundary point includes the start and end points of a boundary line; thus, the lithology of the geological point is obtained.

[0088] Specifically, the UpdateCursor cursor is used to retrieve all data into memory. The query field is set to all fields of the feature class. An "OID@" field is added at the beginning and an "SHAPE@XY" field is added at the end. When traversing, the number of OIDs in the p_num list is first checked. When it equals 1, it is a geological lithology control point. When it equals 2, it is a geological lithology boundary point.

[0089] S75. When the number of OIDs is 2, find the end point / start point paired with the current start point / end point by using the number index according to the current point number and the value of the OID in the p_num list. Subtract the X coordinate of the paired end point / start point from the X coordinate of the current start point / end point to obtain dx. Subtract the Y coordinate of the paired end point / start point from the Y coordinate of the current start point / end point to obtain dy. Call the math.atan2(dx, dy) function to obtain the relative radian of the two points. Call the math.degrees method to convert the radian to an angle value. If the angle value is positive, save the angle value directly. If the angle value is negative, add 360 degrees to get the final angle value and save it.

[0090] S76. Calculate the angle values ​​using GIS to transform them into eight directions, namely North, Northeast, East, Southeast, South, Southwest, West, and Northwest, to obtain the properties of the geological points; and save all calculation results by writing them into the feature class attributes using the UpdateCursor cursor.

[0091] Finally, the SearchCursor cursor is used again to load the data into memory, constructing a list of object result sets in the format [{"attributes": attributes, "geometry": geometry, "analysis": analysis}], which is then provided to the user.

[0092] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for automatically obtaining the properties and lithology of geological points based on a geological region, characterized in that, include: S1: Obtain topologically verified geological region feature layer data and geological point feature layer data in GIS. The geological region feature layer data includes geological body code, lithology, lithology description and project number field where the geological region is located. The geological point feature layer data includes point number of geological point and route number field of geological point route. S2: Based on the project number field and route number field, use the GIS attribute selection tool to extract the geological area elements and geological point elements associated with the geological point to be analyzed from the geological area element layer data and geological point element layer data; S3: Set the buffer radius, perform GIS buffer analysis on geological point features within the buffer radius, and generate a geological buffer. S4: Analyze the geological region features and geological point features in S2 with the geological buffer zone in S3 using a GIS intersection analysis tool to obtain the intersection of the geological region features and geological point features with the geological buffer zone, thus obtaining the number of intersecting features. Based on the number of intersecting features, determine whether the lithology of the geological point to be analyzed and whether the geological map where the geological point is located has a topological error. If the number of intersecting features is less than or equal to 2, the geological map has no topological error; if the number of intersecting features is greater than 2, the geological map has a topological error. Integrate all attribute fields of the geological region features and geological point features in the intersection to obtain a region feature class file. S5: If the judgment result is that there is no topological error on the geological map, the lithology of the geological point to be analyzed is directly judged. If the judgment result is that there is a topological error on the geological map, the area of ​​each zone element located outside the geological buffer zone is calculated, and the areas are sorted from largest to smallest. Only the zone elements with the top two areas are retained. The zone elements are the geological zone elements located outside the geological buffer zone generated after the intersection analysis of the geological zone elements and geological point elements with the geological buffer zone. S6: Use the GIS feature to point conversion tool to convert the retained area features into point features, and retain all attribute fields integrated in the area feature class file; S7: Based on the conversion results, use GIS to calculate the point type of the geological point to be analyzed, and use the geological route as the baseline to calculate the relative azimuth (radians), relative azimuth (degrees), and relative azimuth (octagonal) between the geological point to be analyzed and the baseline, so as to obtain the corresponding geological point properties and geological point lithology analysis results.

2. The method for automatically obtaining the properties and lithology of geological points based on a geological region according to claim 1, characterized in that, The radius R of the buffer zone is equal to a preset value multiplied by the denominator of the mapping scale.

3. The method for automatically obtaining the properties and lithology of geological points based on a geological region according to claim 1, characterized in that, If the judgment result indicates a topological error in the geological map, then calculate the area of ​​each feature located outside the geological buffer zone, sort the areas from largest to smallest, and retain only the top two feature areas, including: S51. Group the area elements with the same geological point number in the intersection results according to the geological point number, and calculate the area of ​​the area elements in each group. S52. For groups with ≥3 elements under the same geological point number, sort the area elements in descending order by area. S53. Elements with a sorted sequence number ≥ 3 are designated as deletion objects; S54. Perform the GIS feature deletion operation to delete the marked area features.

4. The method for automatically obtaining the properties and lithology of geological points based on a geological region according to claim 1, characterized in that, If the result indicates that there are no topological errors on the geological map, then the lithology of the geological points is directly determined, specifically including: When the intersection result shows that there is only one element in the geological buffer zone, the geological point is a lithological control point; When the intersection result shows that there are two elements in the geological buffer zone, the geological point is the lithological boundary point.

5. The method for automatically obtaining the properties and lithology of geological points based on a geological region according to claim 1, characterized in that, Based on the transformation results, GIS was used to calculate the point type of the geological points to be analyzed. Using the geological route as a baseline, the relative azimuth (radians), relative azimuth (degrees), and relative octahedron between the geological points and the baseline were calculated to obtain the corresponding geological point properties and lithological analysis results, including: S71. Construct an SQL query expression based on the geological point route numbers; S72. Use the SearchCursor cursor to read the transformation results based on the SQL query expression and save the read data; the read data includes the OID value of the geological points in the point features, the X and Y coordinates of the geometric center point in the point features, and the geological point number field in the point features; S73. First traversal of the conversion results: Traverse the conversion results, create a point_dict dictionary and a p_num list, store the geological point numbers in the point_dict dictionary, and store the corresponding OID list in the p_num list to form a number index; S74. Perform a second traversal of the conversion results: Use the UpdateCursor cursor to read all data from the point_dict dictionary, traverse each geological point, and determine the number of OIDs stored in the p_num list. When the number of OIDs is 1, the geological point is a geological lithology control point; when the number of OIDs is 2, the geological point is a geological lithology boundary point. The geological lithology boundary point includes the start and end points of a boundary line; thus, the lithology of the geological point is obtained. S75. When the number of OIDs is 2, find the end point / start point paired with the current start point / end point by using the number index according to the current point number and the value of the OID in the p_num list. Subtract the X coordinate of the paired end point / start point from the X coordinate of the current start point / end point to obtain dx. Subtract the Y coordinate of the paired end point / start point from the Y coordinate of the current start point / end point to obtain dy. Call the math.atan2(dx, dy) function to obtain the relative radian of the two points. Call the math.degrees method to convert the radian to an angle value. If the angle value is positive, save the angle value directly. If the angle value is negative, add 360 degrees to obtain the final angle value and save it. S76. Calculate the angle values ​​using GIS to transform them into eight directions: North, Northeast, East, Southeast, South, Southwest, West, and Northwest, thereby obtaining the properties of the geological points.

Images