A method for quickly correcting a vector file

Abstract

The present application relates to a kind of vector file fast correction method, belong to the development and use technical field of geographic information technology.The present application first extracts the picture frame data of original engineering file, then selects two intersection and can read out its projection plane rectangular coordinate kilometer network as calculation line, mark the selected calculation line on engineering file to form a new engineering file, adjust the direction of calculation line, obtain the projection plane rectangular coordinate of intersection, according to the drawing coordinates of calculation line, calculate the angle of rotation, according to the drawing coordinates of intersection of rotated calculation line, scale and the plane rectangular coordinate of intersection before rotation, calculate the distance of translation, finally, the original engineering file is corrected according to the calculated rotation angle and translation distance.The present application has clear idea, simple, fast and accurate operation, makes non-professional personnel more easily get started, is convenient for promotion, improves the correction efficiency of vector file, has definite practical application prospect.

Description

A fast correction method for vector files Technical Field

[0001] This invention relates to a method for rapid correction of vector files, belonging to the technical field of geographic information technology development and application. Background Technology

[0002] MAPGIS is a geographic information system with completely independent intellectual property rights in my country. It has wide applications, with a large amount of existing data in MAPGIS format in fields such as geology, land, and minerals. With the advent of the "big data" era, almost all big data is closely related to spatial data. Due to different data sources and non-standard data, much spatial data, especially existing data, is not in standard locations and needs to be corrected before it can be overlaid and analyzed with other data. File correction has become the most time-consuming and labor-intensive task in the process of cleaning and integrating existing data. Therefore, it is necessary and beneficial to explore a simple, effective, and easy-to-operate correction method.

[0003] Currently, the main method for calibrating vector files is the corresponding control point calibration method. This method involves selecting corresponding points from both the file to be calibrated and the standard location file to generate a control point file, which is then used to calibrate the vector file. While this method successfully calibrates the file, it requires skilled technicians with a strong theoretical foundation, placing high demands on the personnel. Furthermore, for vector files with tilted kilometer grids or incomplete kilometer grid annotations, repeated comparisons between the tilted and upright files are necessary, and corresponding points need to be labeled separately. This process is cumbersome and inefficient. Summary of the Invention

[0004] The purpose of this invention is to provide a fast vector file correction method to solve the problems of cumbersome process and low efficiency in the current correction process.

[0005] To solve the above-mentioned technical problems, this invention provides a method for rapid calibration of vector files, which includes the following steps:

[0006] 1) Extract the drawing frame data from the project file;

[0007] 2) Select two kilometer grids with intersecting points from the drawing frame data and whose projected rectangular coordinates can be read as calculation lines, and combine them with the project file to form a new project file;

[0008] 3) Adjust the calculation line according to the set direction, obtain the coordinates of the calculation line on the drawing, and calculate the required rotation angle based on the coordinates of the calculation line on the drawing;

[0009] 4) Obtain the graphic coordinates of the intersection points of the calculation lines after rotation, and calculate the required translation distance based on the graphic coordinates of the intersection points after rotation, the scale of the project file, and the Cartesian coordinates of the intersection points of the calculation lines before rotation;

[0010] 5) Based on the calculated rotation angle and translation distance, perform rotation and translation operations on the project file to complete the correction.

[0011] This invention first extracts the frame data from the engineering file, then selects two kilometer grids with intersecting points whose projected rectangular coordinates can be read as calculation lines. These selected calculation lines are marked on the engineering file to form a new engineering file. Based on the coordinates of the selected calculation lines, the required rotation angle and translation distance of the engineering file are calculated, allowing the engineering file to be corrected through rotation and translation without a cumbersome process. This invention has a clear concept, is simple to operate, fast and accurate, making it easier for non-professionals to learn and promote. It improves the correction efficiency of vector files and has clear prospects for practical application.

[0012] Furthermore, the vector file fast correction method also includes a step of verifying the correction result of step 5): generating a standard map frame of the same scale, comparing the map frame of the corrected file with the standard map frame to verify the result. If the accuracy requirements are met, the correction is successful; if the accuracy requirements are not met, batch correction is performed using the control point files of the corrected map frame and the standard map frame.

[0013] Control point files can be automatically collected and generated. Using control point files, vector files that do not meet accuracy requirements can be quickly and accurately corrected.

[0014] Furthermore, the direction set in step 3) is from west to east.

[0015] Furthermore, the formula for calculating the rotation angle in step 3) is as follows:

[0016] α0=arctn[(Y2-Y1) / (X2-X1)]; α=-α0

[0017] α is the angle to be rotated, α0 is the angle, X1 and Y1 are the map coordinates of the starting point of the calculation line, and X2 and Y2 are the map coordinates of the ending point of the calculation line.

[0018] Furthermore, the formula for calculating the translation distance in step 4) is as follows:

[0019] ΔX = 1000 * 0 / MX S ΔY = 1000 * 0 / MY S

[0020] ΔX, ΔY are the translation distances, 1 / M is the scale, (X0, Y0) are the Cartesian coordinates of the intersection points before rotation, (X...S Y S ) represents the coordinates of the intersection point of the calculated lines after rotation.

[0021] The calculation formula used in this invention can quickly and accurately calculate the required rotation angle and translation distance without a tedious calculation process, thus improving the efficiency of document correction. Attached Figure Description

[0022] Figure 1 is a flowchart illustrating the implementation process of the vector file fast correction method of the present invention;

[0023] Figure 2 is a technical flowchart of the vector file fast correction method of the present invention;

[0024] Figure 3 is a diagram of the original project file in an embodiment of the present invention;

[0025] Figure 4 is a diagram of the extracted frame data file in an embodiment of the present invention;

[0026] Figure 5-a is a diagram of two intersecting kilometer grids in an embodiment of the present invention;

[0027] Figure 5-b is a calculation line diagram selected in an embodiment of the present invention;

[0028] Figure 6 is a direction diagram of the calculation line in an embodiment of the present invention;

[0029] Figure 7-a is a diagram showing the determination of the scale and the projection coordinates of the intersection point in an embodiment of the present invention;

[0030] Figure 7-b is a diagram showing the coordinates of the intersection points after rotation in an embodiment of the present invention;

[0031] Figure 8-a is a diagram showing the rotation parameter settings in an embodiment of the present invention;

[0032] Figure 8-b is a diagram showing the rotation result in an embodiment of the present invention;

[0033] Figure 8-c is a diagram showing the translation parameter settings in an embodiment of the present invention;

[0034] Figure 8-d is a translation result diagram in an embodiment of the present invention;

[0035] Figure 9 is a comparison diagram of the correction frame and the standard frame in an embodiment of the present invention;

[0036] Figure 10 is a diagram of the correction results in an embodiment of the present invention. Detailed Implementation

[0037] The specific embodiments of the present invention will be further described below with reference to the accompanying drawings.

[0038] Examples of a fast correction method for vector files

[0039] This invention first extracts the frame data from the original project file; then, it selects two kilometer grids with intersecting points whose projected rectangular coordinates can be read, cuts them, and extracts the calculation lines. The selected calculation lines are marked on the project file to form a new project file. The direction of the calculation lines is adjusted, and the required rotation angle is calculated based on the map coordinates of the calculation lines. The required translation distance is calculated based on the map coordinates of the intersection points of the rotated calculation lines, the scale, and the original rectangular coordinates of the intersection points. The original project file is then corrected according to the calculated rotation angle and translation distance. Finally, the corrected file is checked to ensure it meets the accuracy requirements. The implementation process of this method is shown in Figure 1, and the technical flow is shown in Figure 2. The following explanation uses specific examples.

[0040] 1. Extract the drawing frame data from the original project file.

[0041] Open the original project file, extract the title block data, and save it as a new file. If a north arrow is present, it can be saved as well. The original project file selected in this embodiment is shown in Figure 3, and the extracted title block data file is shown in Figure 4.

[0042] 2. Select the calculation line and create a new project file, which together with the existing project file forms a new project file.

[0043] Two kilometer grids with intersecting points and whose projected Cartesian coordinates can be read are selected from the extracted map frame data file as calculation lines. In this embodiment, lines meeting the above conditions are selected from the map frame data file. These lines are then cut into four lines at their intersection points. The line along a predetermined direction on one side of the intersection point is selected as the calculation line, and the selected calculation lines are marked on the project file to form a new project file. The two intersecting kilometer grids selected according to the above method in this embodiment are shown in Figure 5-a. The east-west line to the right of the intersection point is selected as the calculation line, and the result is shown in Figure 5-b, where TKN is the east-west line. Alternatively, the line to the left of the intersection point can also be selected.

[0044] 3. Adjust the direction of the calculation line to the set direction, obtain the graph coordinates of the calculation line, and calculate the required rotation angle based on the graph coordinates of the calculation line.

[0045] Adjust the direction of the calculation line to the set direction, obtain the coordinates of the start and end points of the calculation line from the plaintext file, and use this to calculate the rotation angle by which the calculation line should be rotated to the set direction.

[0046] Formula for calculating rotation angle:

[0047] α0=arctn[(Y2-Y1) / (X2-X1)]; α=-α0

[0048] α is the angle to be rotated, α0 is the angle, X1 and Y1 are the map coordinates of the starting point of the calculation line, and X2 and Y2 are the map coordinates of the ending point of the calculation line. The map coordinates are read from a plaintext file exported by MAPGIS software. In this embodiment, the direction of the calculation line is from west to east and does not need to be adjusted, as shown in Figure 6. The map coordinates of the calculation line are: starting point: 7719.188, 747630.568; ending point: 7884.704, 747733.995; after calculation, α0 = 32.00°, and the rotation angle is: α = -α0 = -32°.

[0049] 4. Calculate the distance to be translated based on the map coordinates of the intersection points of the rotated calculation lines, the scale, and the Cartesian coordinates of the intersection points before rotation.

[0050] Calculation of translation distance:

[0051] ΔX = 1000 * X0 / MX S ΔY = 1000 * Y0 / MY S

[0052] ΔX, ΔY are the translation distances, 1 / M is the scale, (X0, Y0) are the Cartesian coordinates of the intersection points before rotation, (X... S Y S The coordinates of the intersection points of the calculated lines after rotation are shown on the map. The map coordinates are read from a plaintext file exported by MAPGIS software.

[0053] In this embodiment, the scale is 1:5000, M = 5000, and the unit is mm; the rectangular coordinates of the intersection point are 38439500 and 3743000, respectively, in meters; the MapGIS projection plane does not include the projection zone number, and the coordinates after removing the zone number are 439500 and 3743000; the scale and map coordinates of the intersection point are shown in Figure 7-a. The rotated map coordinates of the intersection point are read as 402730.083 and 629936.133, as shown in Figure 7-b. The translation distance is calculated as follows:

[0054] ΔX=1000*439500 / 5000-402730.083=-314830.083

[0055] ΔY=1000*3743000 / 5000-629936.133=118663.867

[0056] 5. Correct the original project file according to the calculated rotation angle and translation distance.

[0057] First, rotate the project file by the rotation angle, with the rotation parameters set as shown in Figure 8-a, and the rotation result is shown in Figure 8-b; then, translate the rotated file, with the translation parameters set as shown in Figure 8-c, and the translation result is shown in Figure 8-d.

[0058] 6. Verification of calibration results.

[0059] A standard drawing frame with the same scale as the original project file is generated. The drawing frame of the corrected file is compared with the standard drawing frame to verify the result. If the accuracy requirement is met (i.e., the drawing frame of the corrected file completely overlaps with the standard drawing frame), the correction is successful. If the accuracy requirement is not met, batch correction is performed using the control point files of the corrected drawing frame and the standard drawing frame. In this embodiment, the comparison between the corrected drawing frame and the standard drawing frame is shown in Figure 9. The corrected file is completely aligned with the generated 1:5000 standard drawing frame. The final correction result is shown in Figure 10.

Claims

1. A method for fast correction of vector files, characterized in that, The correction method includes the following steps: 1) Extracting the drawing frame data from the project file; 2) Selecting two kilometer grids with intersecting points and whose projected rectangular coordinates can be read from the drawing frame data as lines that meet the conditions, and selecting the line along the set direction on one side of the intersection point as the calculation line, and combining it with the project file to form a new project file; 3) Adjusting the calculation line according to the set direction, obtaining the drawing coordinates of the calculation line, and calculating the required rotation angle based on the drawing coordinates of the calculation line. The formula for calculating the rotation angle is: For the angle that needs to be rotated, For the angle, ( ， ) represents the map coordinates of the starting point of the calculation line. ， 1) Obtain the coordinates of the endpoint of the calculation line on the drawing; 4) Obtain the coordinates of the intersection point of the calculation line after rotation, and calculate the distance to be translated based on the coordinates of the intersection point after rotation, the scale of the project file, and the Cartesian coordinates of the intersection point of the calculation line before rotation; 5) Perform rotation and translation operations on the project file to complete the correction based on the calculated rotation angle and translation distance.

2. The vector file fast correction method according to claim 1, characterized in that, The method also includes a step of verifying the correction results of step 5): generating a standard drawing frame of the same scale, comparing the drawing frame of the corrected file with the standard drawing frame to verify the results. If the accuracy requirements are met, the correction is successful. If the accuracy requirements are not met, batch correction is performed using the control point files of the corrected drawing frame and the standard drawing frame.

3. The vector file fast correction method according to claim 1 or 2, characterized in that, The direction set in step 3) is from west to east.

4. The vector file fast correction method according to claim 1, characterized in that, The map coordinates in step 3) are obtained by reading the plaintext file exported from MAPGIS software.

5. The vector file fast correction method according to claim 1 or 2, characterized in that, The formula for calculating the translation distance in step 4) is as follows: , The translation distance is... For the scale, ( ， ( ) are the Cartesian coordinates of the points of intersection before rotation, ( ， () represents the coordinates of the intersection point on the graph after rotation.

Images