Geographic space coordinate system identification method based on spatial range verification and visual confirmation

Abstract

The invention relates to the technical field of processing of various coordinate systems, and discloses a geographic space coordinate system identification method based on spatial range verification and visual confirmation. The method comprises the following steps: judging the type of an initial coordinate system according to information of coordinate points to be identified, further analyzing according to the type of the initial coordinate system to obtain a coordinate system identification library, and performing coordinate conversion and range verification on the initial coordinate system in the coordinate system identification library to obtain a target data coordinate set; and further loading a basic map based on the target data coordinate set and performing visualization processing to identify a final coordinate system name, so that a candidate coordinate system is screened through range checking by combining a user-guided sequence process, range verification and visualization confirmation, and then a result is presented in a visualization mode for a user to intuitively compare and judge, thereby improving the user experience. And the identification accuracy and reliability are improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of coordinate system processing, and in particular to a geographic coordinate system identification method based on spatial range verification and visual confirmation. BACKGROUND

[0002] In geographic spatial data processing, accurate identification of coordinate reference system (CRS) is the core prerequisite for data integration, analysis and visualization. Coordinate reference system includes geographic coordinate system (GCS) and projection coordinate system (PCS), which establishes a unified reference framework for spatial data. However, diversified data sources are often accompanied by problems of missing, incorrect or incomplete metadata, making it difficult for users to determine the true coordinate system of the data. Such identification deviation will cause serious consequences: for example, misuse of coordinate system in GIS analysis will cause point displacement, resulting in invalid spatial analysis; in the automatic driving scene, it may trigger navigation accidents due to coordinate errors, threatening safety; in the unmanned aerial vehicle operation, it is easy to cause track deviation due to coordinate system confusion, reducing task reliability.

[0003] Currently, the manual inspection or professional software (such as ArcGIS) identification scheme relied on by the industry has significant limitations: first, the operation is tedious and time-consuming, requiring users to be proficient in various coordinate system parameters and definitions; second, the high threshold of professional ability excludes non-technical users; third, the quality defects of metadata itself will directly lead to identification failure; fourth, even if the identification is completed, the cross-coordinate system data conversion still needs to accurately match the source system and target system parameters, the process is complex and the fault tolerance is low.

[0004] Therefore, there is an urgent need for a method that can efficiently and reliably achieve accurate identification of coordinate system. SUMMARY

[0005] The present application provides a geographic coordinate system identification method based on spatial range verification and visual confirmation, which can efficiently and reliably achieve accurate identification of coordinate system.

[0006] The present application provides a geographic coordinate system identification method based on spatial range verification and visual confirmation, which includes: S1: obtaining coordinate point information to be identified and background information of an initial coordinate system corresponding to the coordinate point information to be identified; S2: determining the type of the initial coordinate system according to the coordinate point information to be identified, wherein the type of the initial coordinate system includes a geographic coordinate system and a projection coordinate system; S3: analyzing the initial coordinate system according to the type of the initial coordinate system to obtain a coordinate system identification library; S4: performing coordinate conversion and range verification on the initial coordinate system in the coordinate system identification library to obtain a target data coordinate set; S5: loading a basic map based on the target data coordinate set with a specified spatial range as the center, drawing the coordinate points converted in coordinate on the basic map, and performing visual processing S6: identifying the final coordinate system name by comparing the coordinate points converted in coordinate of different initial coordinate systems through the visual result.

[0007] The present application has the following beneficial effects: The geographic coordinate system identification method based on spatial range verification and visual confirmation provided in the present application judges the type of the initial coordinate system according to the information of the coordinate points to be identified, further analyzes the coordinate system identification library according to the type of the initial coordinate system, performs coordinate conversion and range verification on the initial coordinate system in the coordinate system identification library to obtain a target data coordinate set, further loads a basic map based on the target data coordinate set and performs visual processing to identify the final coordinate system name. In this way, through the combination of the sequential process guided by the user, range verification and visual confirmation, the candidate coordinate system is first screened through range checking, and then the result is presented in a visual manner for the user to intuitively compare and judge, thereby improving the accuracy and reliability of the identification.

[0008] The present application will be further described below with reference to the drawings. BRIEF DESCRIPTION OF DRAWINGS

[0009] The accompanying drawings, which form a part of this application, are included to provide a further understanding of the application, and are incorporated in and constitute a part of this application. The embodiments illustrated in the drawings are provided to explain the present application and should not be considered limiting of the present application. In the drawings: Figure 1 is one of the flowcharts of the geographic coordinate system identification method based on spatial range verification and visual confirmation according to a preferred embodiment of the present application; Figure 2 is another flowchart of the geographic coordinate system identification method based on spatial range verification and visual confirmation according to a preferred embodiment of the present application; Figure 3 is a commonly used geographic coordinate system and part of the projection coordinate system in A country according to a preferred embodiment of the present application. DETAILED DESCRIPTION

[0010] The technical solutions of the present application will be described in detail below. Obviously, the described embodiments are only a part of the embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0011] Unless otherwise defined, technical terms or scientific terms used in the present application shall have the ordinary meanings to those skilled in the art to which the present application pertains. The terms "first", "second", and similar terms used in the present application do not denote any order, quantity, or importance, but are used to distinguish different components. Similarly, the terms "one" or "a" or similar terms do not denote a quantity restriction, but mean at least one.

[0012] Currently, in the field of geospatial data processing, due to the lack of metadata, record errors or non-standard management, the identification of coordinate system has long been a serious challenge. Professionals need to rely on experience to try repeatedly in professional software, which is tedious and inefficient; non-professionals are almost unable to complete independently. The existing technology relies too much on the completeness of metadata, and lacks intuitive result verification mechanism.

[0013] It should be understood that the present application covers a wide range of applications in the fields of geographic information systems, remote sensing technology, surveying and mapping, autonomous driving, and unmanned vehicle navigation, and can process geospatial data of different scales and types worldwide. Especially in the absence of explicit coordinate system information, the present application can effectively identify the coordinate system through coordinate feature analysis and range constraints.

[0014] Please refer to Figures 1-2 The present application provides a geospatial coordinate system identification method based on spatial range verification and visualization confirmation, comprising: S1: obtaining coordinate point information to be identified and background information of an initial coordinate system corresponding to the coordinate point information to be identified.

[0015] In this step, the coordinate point information to be identified can include a group of coordinate points or multiple groups of coordinate points These coordinate points can correspond to longitude / latitude in the geographic coordinate system (GCS) or east / north offset in the projected coordinate system (PCS). Meanwhile, there can be cases of order reversal, confusion, and offset encryption.

[0016] S2: determining the type of the initial coordinate system according to the coordinate point information to be identified, wherein the type of the initial coordinate system includes a geographic coordinate system and a projected coordinate system.

[0017] S3: analyzing the initial coordinate system according to the type of the initial coordinate system to obtain a coordinate system identification library.

[0018] S4: performing coordinate transformation and range verification on the initial coordinate system in the coordinate system identification library to obtain a target data coordinate set.

[0019] S5: Loading a basic map based on the target data coordinate set with a specified spatial range as the center, drawing the coordinate points after coordinate conversion on the basic map, and performing visual processing.

[0020] S6: Identifying the final coordinate system name by comparing the coordinate points after conversion of different initial coordinate systems through the visual results.

[0021] In this step, the final coordinate system name refers to the specific name of the coordinate system to which the coordinate points belong, for example, "CGCS2000_3_degree_Gk_Zone_36, Beijing_1954_3_degree_GK_Zone_36, Xian_1980_3_degree_GK_Zone_36". This is only an example and is not limited. In actual situations, geospatial data, due to its diverse sources, will inevitably involve a variety of coordinate reference systems established in different historical periods, based on different measurement benchmarks and applicable purposes, in various parts of the world. Therefore, specifying a clear coordinate system name for coordinate points is a fundamental prerequisite to ensure data readability, convertibility and interoperability.

[0022] The above-mentioned geospatial coordinate system identification method based on spatial range verification and visual confirmation determines the type of initial coordinate system according to the information of the coordinate points to be identified, further analyzes the coordinate system identification library according to the type of the initial coordinate system, performs coordinate conversion and range verification on the initial coordinate system in the coordinate system identification library to obtain a target data coordinate set, and further loads a basic map based on the target data coordinate set and performs visual processing to identify the final coordinate system name. In this way, by combining the sequential process guided by the user, range verification and visual confirmation, the candidate coordinate systems are first screened through range checking, and then the results are presented in a visual manner for the user to intuitively compare and judge, thereby improving the accuracy and reliability of the identification.

[0023] Optionally, the S1 comprises: S11: Obtaining information of coordinate points to be identified ; S12: Obtaining a spatial range in which the information of the coordinate points to be identified approximately exists , ; S13: Determining a target area based on the information of the coordinate points to be identified

[0024] In this optional embodiment, the spatial range in which the information of the coordinate points to be identified approximately exists refers to the range in latitude and longitude. The spatial coordinate points to be verified refer to coordinate points whose coordinate systems are temporarily unknown, which are one of the main input data and can be a single coordinate point or multiple groups of coordinate points The spatial range of coordinate points mainly refers to a rough estimate of the range of coordinate points. This range has tolerance; it can cover both small and large areas, but latitude and longitude in WGS84 coordinates are required. The selected target area mainly refers to a large region or country, which is used to determine the possible coordinate system to be used. This is fundamental information that is often considered in coordinate recognition. Custom coordinate systems can also be defined (e.g., the proj4 string), and these custom coordinate systems can be directly included in the candidate list for verification.

[0025] In addition, you can enter a specified coordinate type, mainly to confirm the input. This step, specifying whether to use a geographic coordinate system (GCS) or a projected coordinate system (PCS), reduces the need for subsequent data type determination. For example... Figure 3 The image shows the commonly used geographic coordinate system and some projected coordinate systems in Country A. Specifying the coordinate type refers to the direct qualitative analysis of the input coordinates. Since geographic coordinates and projected coordinates differ significantly in format and data size, their properties can be directly compared, but this requires relevant knowledge. Therefore, you can choose to define your own coordinates or use the default system definition. If no input is provided, the analysis will proceed in subsequent steps.

[0026] Optionally, S2 includes: The initial coordinate system type is inferred based on the numerical range of the coordinate point information to be identified. Coordinates of all points If the set range conditions are met, the initial coordinate system type is determined to be a geographic coordinate system, and the coordinate format is converted to degrees; if the set range conditions are not met, it is determined to be a projected coordinate system. The range conditions are as follows: ; or: .

[0027] In this optional implementation, the determination of coordinate type mainly considers that geographic coordinates are in latitude and longitude form, and latitude and longitude have clear data boundaries. Since projected data is generally large in volume, determining whether it falls within a geographic coordinate range in this step can quickly distinguish data types, unify units, and facilitate subsequent data processing.

[0028] Optionally, S3 includes: A commonly used coordinate series table is determined based on the target area, and the commonly used coordinate series table is incorporated into a custom coordinate system to form a coordinate system identification library; The process of defining a custom coordinate system is as follows: If the type is a geographic coordinate system, firstly consider whether the geographic coordinate system is an encrypted coordinate system, and then count the number of bits of the integer part of the input coordinates According to the number of bits of the integer part, select whether to integrate the numbered projection coordinate system into the coordinate system identification library; If the type is a projection coordinate system, count the number of bits of the integer part of the input coordinates According to the number of bits of the integer part, select whether to integrate the numbered projection coordinate system into the coordinate system identification library; According to the number of bits of the integer part, select whether to integrate the numbered projection coordinate system into the coordinate system identification library, including: If the number of bits of the integer part of one of x and y is 7 and the number of bits of the integer part of the other is 8, then select the numbered projection coordinate system to integrate into the coordinate system identification library; If the number of bits of the integer part of one of x and y is 6 and the number of bits of the integer part of the other is 7, then select the unnumbered projection coordinate system to integrate into the coordinate system identification library.

[0029] In this embodiment, first determine the type of the coordinate system, specifically whether it is a geographic coordinate system or a projection coordinate system, further determine whether the geographic coordinate system is an encrypted coordinate system, and then perform the step of selecting whether to integrate the numbered projection coordinate system into the coordinate system identification library according to the number of bits of the integer part. In this way, by identifying the coordinate type and whether it is encrypted in the early stage, some unnecessary data processing processes can be reduced, the calculation amount can be reduced, and the identification process can be accelerated.

[0030] According to the target area, determine a list of commonly used coordinate systems, for example, A country mainly uses WGS84 (EPSG: 4326), CGCS2000 (EPSG: 4490), Gauss-Kruger projection (EPSG: 4526), etc., and B country mainly uses WGS84 (EPSG: 4326), NAD83 (EPSG: 4269), UTM zone (such as EPSG: 32614), etc. Since the definitions of various geographic coordinate systems and projection coordinate systems are closely related to the location of the region, the commonly used geographic or projection coordinate systems in the corresponding region are basically of a fixed type.

[0031] Specifically, unencrypted coordinate systems such as WGS84, CGCS2000; encrypted coordinate systems such as commonly used GCJ-02 coordinate system, BD-09 coordinate system. Select the projection coordinate system with number and incorporate into the identification library, such as Xian_1980_GK_Zone_20, Xian_1980_3_Degree_GK_Zone_33; Select the projection coordinate system without number and incorporate into the identification library, such as Xian_1980_GK_CM_105E, Xian_1980_3_Degree_GK_CM_108E.

[0032] In this step, the information of the target area provided in the previous step is used, and the possible coordinate system is screened according to the preliminary analysis of the coordinates, and the self-defined input coordinate system is combined into the coordinate system data identification library, and the subsequent verification will also be carried out around the coordinate system data identification library.

[0033] Optionally, the S4 comprises: S41: increase the coordinate and replace the coordinate point position: At the same time, the verification coordinate set is incorporated: ; S42: for unencrypted coordinate system conversion, use the standard coordinate conversion method for conversion, project the input coordinates from the candidate coordinate system to WGS84, and the conversion formula is as follows: ; Wherein, is a standard coordinate conversion function, is a candidate coordinate system, is an input coordinate, is the corresponding converted WGS84 latitude and longitude; For encrypted coordinate system conversion, use the coord-convert library to project the encrypted coordinates to WGS84, and the conversion formula is as follows: ; S43: check whether the converted coordinates fall within the specified WGS84 space range, and the verification formula is as follows: ; Wherein, represents logical and operation, represents the data set of points in the range after the projection of the coordinate point corresponding to the coordinate system; Take as the target data coordinate set.

[0034] In this embodiment, the input coordinates may be reversed or confused, so the coordinates are increased and the coordinate point positions are replaced. The coordinate data is exchanged to increase the fault tolerance of the data. The proj4 is used to realize the coordinate system transformation of the coordinate point. It is worth explaining that the proj4 is a widely used coordinate conversion library in geographic information system (GIS), which can realize the conversion between different coordinate reference systems (CRS), cover map projection and geodetic coordinate transformation, support more than 100 kinds of map projection methods, and provide key support for various applications involving geographic spatial data processing. The coord-convert library is a special python library for converting encrypted coordinate system to WGS84, which can realize the conversion of encrypted coordinates with an actual measurement error of less than 2m, and is suitable for most GIS application scenarios. Then, based on the numerical characteristics and spatial range of the WGS84 geographic coordinates, the converted data is screened and identified to obtain the coordinate system information that meets the range correction.

[0035] Optionally, the S5 comprises: Loading a basic map based on the target data coordinate set and taking the specified spatial range as the center; For each candidate coordinate system, the converted coordinate points are plotted on the basic map, different point sets corresponding to different coordinate systems are distinguished by different marking methods, and position labels are added to the basic map to specify the spatial range boundary.

[0036] In this step, the specified spatial range refers to the specified latitude and longitude range; the basic map can be loaded in the form of a Web map service provided by OpenStreetMap, or a local remote sensing image uploaded. For each retained candidate coordinate system, the converted coordinate points are plotted on the map, and different point sets corresponding to different coordinate systems are distinguished by different colors or symbols. Adding position labels to the map to specify the spatial range boundary can strengthen the range reference.

[0037] In this step, the converted WGS84 data points are loaded using online maps or local maps, and the coordinate system information before conversion is attached, which can effectively display the positions of the coordinate points and intuitively understand the position information of the data points through data visualization, improving the data readability.

[0038] Optionally, the S6 comprises: Based on the visualization result, the most suitable coordinate system point in the data coordinates is selected according to the judgment basis, the final coordinate system name is identified according to the most suitable coordinate system point, and the corresponding parameter code is determined.

[0039] The judgment basis includes: The matching degree of the converted point with the spatial distribution of the known geographical features; the alignment accuracy of the converted point with other existing data layers; whether the actual position of the converted point conforms to geographical logic.

[0040] In this optional implementation, the known geographical features are, for example, roads and rivers; the existing data layers are, for example, vector layers and remote sensing images; and the geographical logic is, for example, that points in urban areas should not be located in the ocean.

[0041] Finally, by clicking to determine the final corresponding coordinate system, the recognized coordinate system name and its corresponding EPSG code (such as “WGS84 (EPSG: 4326)”) are obtained, and the correct coordinate points and coordinate system are finally exported in TXT format or JSON format. This is only an example and is not limited.

[0042] In summary, the present application combines user-guided sequential processes with system automation processing, integrates multiple-link operations such as coordinate input, type judgment, candidate screening, conversion verification, and visual presentation, constructs a “input - analysis - screening - confirmation - output” whole-process coordinate reference system identification system, accurately covers the identification needs of multiple coordinate systems in the global range through standard conversion technology and range verification, and especially fills the data processing gap in the scenario lacking clear coordinate system information. In specific implementation, the input coordinates are analyzed for numerical range to determine the type, the proj4 library and the coord-convert library are used to realize coordinate conversion and range checking, and the coordinate distribution of different candidate coordinate systems is simultaneously presented through map visualization. The user completes the confirmation based on geographical feature matching degree, etc., providing reliable technical support for accurate identification of coordinate reference systems.

[0043] Next, an example of an automated identification and visual selection platform based on python is described, and the coordinate data of the verified coordinate system is used for verification.

[0044] The known coordinate data in the example is (3284973.21, 36441039.05). According to data collection, the data point is the endpoint of a river section of a certain river in a certain city, and the specific position and explicit coordinate system are unknown (the coordinate system is Beijing_1954_3_degree_GK_Zone_36 projection coordinates after manual traversal and query). Professional analysis of the coordinate data can determine that the data point is a projection coordinate, and there is a data bit inversion, i.e., the data format is .

[0045] Combined with part of the information of the known coordinate point, the input is set as the coordinate point (3284973.21, 36441039.05), and it is known that the range of the entire city in the WGS84 coordinate system is:

[0046] According to the actual position of a river, the final selection range is , . The known coordinates are projection coordinates, the region belongs to country A, and the corresponding known conditions are input. The number of coordinate bits is counted, and the integer parts of the coordinates are 7 and 8 bits, respectively, so a numbered projection coordinate system is selected and incorporated into the recognition library. Through the definition of the coordinate points, projection, projection transformation, range judgment and other operations are performed.

[0047] This example does not add a custom coordinate system, and an online map is selected for data display. Through calculation, it is found that there are three coordinate systems that belong to the interval range, which are CGCS2000_3_degree_Gk_Zone_36, Beijing_1954_3_degree_GK_Zone_36 and Xian_1980_3_degree_GK_Zone_36, which meet the screening conditions. Through the visualization map, the positions of the three data points are displayed, and it is found that the coordinate points of the Beijing_1954_3_degree_GK_Zone_36 projection are closer to the real demand (the river section endpoints are located on both sides of the river). Finally, the coordinate system of the determined metadata is clicked. Subsequently, the relevant coordinate points and coordinate systems can be exported.

[0048] As can be seen, the present application performs data analysis on coordinate points of unknown coordinate systems and range determination, and finally obtains three coordinate systems closest to the coordinate points. By marking the coordinate points under the transformation of the three coordinate systems on the visual online map, the final coordinate system of the coordinate points is determined, and the practical value of the method in similar regions is verified.

[0049] In summary, the present application performs multi-dimensional data analysis on coordinate points of unknown coordinate systems, including coordinate value range judgment, data type feature extraction and bit inversion identification, and combines spatial range constraints and geographical region limitations to construct a hierarchical screening mechanism, accurately locates and selects the three most likely candidate coordinate systems. On this basis, online map services are called to perform visual mapping, and the coordinate points transformed by different coordinate systems are intuitively displayed on the geographical background layer to form a multi-dimensional spatial comparison. Based on the geographical feature matching degree, the logic rationality of geographical objects and professional knowledge, the candidate coordinate systems are comprehensively evaluated and decided. This process not only realizes the intelligent recognition of the coordinate reference system, but also significantly improves the accuracy and reliability of the coordinate system determination through the verification mechanism of human-computer collaboration. The example shows that this method can maintain the recognition accuracy in complex scenarios lacking clear metadata, effectively solving the problems of long time consumption and easy errors in traditional methods relying on manual experience, and verifying the practical value and promotion potential of the method in similar regional geographical spatial data processing.

[0050] The above merely provides the specific implementation of the present disclosure, but the protection scope of the present disclosure is not limited thereto, any person skilled in the art can easily think of the changes or replacements within the technical range disclosed by the present disclosure, which should be covered in the protection scope of the present disclosure. Therefore, the protection scope of the present disclosure should be subject to the protection scope of the claims.

Claims

1. A geospatial coordinate system identification method based on spatial extent verification and visual confirmation, characterized in that, Comprise: S1: obtain the to be identified coordinate point information and the background situation of the initial coordinate system corresponding to the to be identified coordinate point information; S2: determine the type of the initial coordinate system according to the to be identified coordinate point information, wherein the type of the initial coordinate system comprises a geographic coordinate system and a projection coordinate system; S3: analyze the initial coordinate system according to the type of the initial coordinate system to obtain a coordinate system identification library; S4: for the initial coordinate system in the coordinate system identification library, coordinate conversion and range verification are carried out to obtain a target data coordinate set; S5: based on the target data coordinate set, a basic map is loaded with a specified spatial range as the center, the converted coordinate points are drawn on the basic map, and visual processing is performed; S6: the final coordinate system name is identified by comparing the converted coordinate points of different initial coordinate systems through the visual results.

2. The spatial range based validation and visualization confirmation enabled geospatial coordinate system identification method of claim 1, wherein, The S1 comprises: S11: Obtain coordinate point information to be identified ; S12: Obtain a space range in which the to-be-identified coordinate point information is approximately located. , ; S13: determine the target area based on the to be identified coordinate point information.

3. The spatial range based verification and visualization confirmation enabled geospatial coordinate system identification method of claim 1, wherein, The S2 comprises: The initial coordinate system type is inferred based on the numerical range of the coordinate point information to be identified. Coordinates of all points If the set range conditions are met, the initial coordinate system type is determined to be a geographic coordinate system, and the coordinate format is converted to degrees; if the set range conditions are not met, it is determined to be a projected coordinate system. Wherein, the range condition is as follows: ； Or: 。 4. The spatial range based verification and visualization confirmation enabled geospatial coordinate system identification method of claim 2, wherein, The S3 comprises: Determine a commonly used coordinate system list according to the target area, and incorporate the commonly used coordinate system list into a self-defined coordinate system to form a coordinate system identification library; Wherein, the process of self-defining the coordinate system is as follows: If the type is geographic coordinate system, firstly consider whether the geographic coordinate system is encrypted coordinate system, and then input the coordinates Statistical bit number of integer part, select whether to integrate the numbered projection coordinate system into the coordinate system identification library according to the statistical results of the bit number of the integer part; If the type is a projection coordinate system, the input coordinates are counted for the number of bits of the integer part, and a projection coordinate system with or without a number is selected into the coordinate system identification library according to the counting result of the number of bits of the integer part. Wherein, according to the bit number statistical result of the integer part, it is determined whether to incorporate the numbered projection coordinate system into the coordinate system identification library, comprising: If the bit number of the integer part of one of x and y is 7, and the bit number of the integer part of the other is 8, then the numbered projection coordinate system is selected to be incorporated into the coordinate system identification library; If the bit number of the integer part of one of x and y is 6, and the bit number of the integer part of the other is 7, then the unnumbered projection coordinate system is selected to be incorporated into the coordinate system identification library.

5. The spatial range based verification and visualization confirmation enabled geospatial coordinate system identification method of claim 1, wherein, The S4 comprises: S41: increase the coordinates and replace the coordinate point position: incorporate the verification coordinate set simultaneously: ； S42: for the unencrypted coordinate system conversion, use the standard coordinate conversion method to convert, project the input coordinates from the candidate coordinate system to WGS84, and the conversion formula is as follows: ； wherein, is a standard coordinate conversion function, is a candidate coordinate system, is an input coordinate, is the corresponding converted WGS84 longitude and latitude; For encrypted coordinate system conversion, use the coord-convert library to project the encrypted coordinates to WGS84, and the conversion formula is as follows: ； S43: check whether the converted coordinates fall within the specified WGS84 spatial range, and the verification formula is as follows: ； wherein, represents a logical AND operation, represents a data set of points in the range after projection of the coordinate points corresponding to the coordinate system; To as the target data coordinate set.

6. The spatial range based validation and visualization confirmation enabled geospatial coordinate system identification method of claim 5, wherein, The S5 comprises: Based on the target data coordinate set, a basic map is loaded with a specified spatial range as the center; For each candidate coordinate system, the converted coordinate points are drawn on the basic map, different point sets corresponding to different coordinate systems are distinguished by different marking methods, and position annotations are added on the basic map to specify the spatial range boundary.

7. The spatial range based validation and visualization confirmation enabled geospatial coordinate system identification method of claim 1, wherein, The S6 comprises: Based on the visual results, the point of the coordinate system that best fits the data coordinates is selected according to the judgment basis, the final coordinate system name is identified according to the point of the coordinate system that best fits, and the corresponding parameter code is determined.

8. The spatial range based verification and visual confirmation enabled geospatial coordinate system identification method as claimed in claim 7, wherein, The judgment basis comprises: The matching degree of the converted points and the spatial distribution of known geographic features; the alignment accuracy of the converted points and other existing data layers; whether the actual position of the converted points conforms to the geographic logic.

Images