Drawing method, device and equipment for dynamic evaluation of GNSS station data integrity and medium

Abstract

The invention provides a drawing method, device and equipment for dynamic evaluation of GNSS station data integrity and a medium. The drawing method comprises the following steps: grouping a current data set according to a near space; judging whether current abnormal data exists in the current data set of each group or not; if yes, acquiring data of the GNSS site with the current abnormal data in a preset time before the current moment; judging whether the data is abnormal data; if yes, discarding the current abnormal data in the current data set, and determining that the to-be-judged sensor data set is a data set obtained after the GNSS site data set of the current abnormal data is excluded from the current data set; if not, determining that the to-be-judged sensor data set is the current data set; judging whether a GNSS site corresponding to the to-be-judged sensor data set is greater than a preset proportion of a total GNSS site or not; if yes, grid interpolation, space contour surface drawing and tile slicing are carried out. According to the method, abnormal data can be eliminated, and the accuracy of GNSS site meteorological data drawing is improved.

Description

Technical Field

[0001] The present invention relates to the field of meteorological monitoring technology, and in particular to a drawing method, device, electronic device and computer-readable storage medium for dynamic evaluation of GNSS site data integrity. Background Art

[0002] Multiple stations are installed within a certain geographical space (generally, the interval between stations is 5km to 15km). Each station contains GNSS (Global Navigation Satellite System) observation equipment and meteorological sensor equipment. The data collected by the GNSS observation equipment can be used to calculate the atmospheric precipitation content, and the meteorological sensor equipment can directly collect data such as precipitation, humidity, temperature, wind speed, and wind direction around the current station. Since the equipment may fail after long-term operation, resulting in abnormal GNSS observation data and meteorological collection data, local abnormalities will occur when drawing meteorological maps of sensors such as atmospheric precipitation, precipitation, humidity, temperature, wind speed and direction. Therefore, how to eliminate abnormal data and then draw meteorological maps has become an urgent problem to be solved. Summary of the invention

[0003] The present invention aims to provide a mapping method, device, electronic device and computer-readable storage medium for dynamic evaluation of GNSS site data integrity that overcome the above problems or at least partially solve the above problems.

[0004] To achieve the above object, the technical solution of the present invention is specifically implemented as follows: A first aspect of the present invention provides a mapping method for dynamic evaluation of GNSS site data integrity, comprising: Obtain the current data set collected by the meteorological environment sensors of N GNSS stations within the preset spatial range at the current moment; Taking each of the GNSS stations as the center of a circle, grouping the GNSS stations within a preset radius to obtain N groups; Processing the current data set collected by the meteorological environment sensor of each group to determine whether there is current abnormal data in the current data set of each group; If the current abnormal data exists, obtaining data of the GNSS site having the current abnormal data at a preset time before the current moment; Determine whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data; If the data of the GNSS site with the current abnormal data at a preset time before the current moment is abnormal data, the current abnormal data in the current data set is discarded, and the sensor data set to be judged is determined, wherein the sensor data set to be judged is the data set after the current data set excludes the GNSS site data set of the current abnormal data; if the data of the GNSS site with the current abnormal data at a preset time before the current moment is not abnormal data, the sensor data set to be judged is determined, wherein the sensor data set to be judged is the current data set; Determine whether the GNSS site corresponding to the sensor data set to be determined is greater than a preset ratio of the total GNSS sites; If it is greater than the preset ratio, the sensor data set to be judged is subjected to grid interpolation, and spatial isosurface drawing is performed to perform tile slicing.

[0005] Optionally, the method further comprises: If the current abnormal data does not exist, the sensor data set to be determined is determined, wherein the sensor data set to be determined is the current data set.

[0006] Optionally, the processing of the current data set collected by the meteorological environment sensor of each group to determine whether there is current abnormal data in the current data set of each group includes: Using the Z-Score function to process the current data set collected by the meteorological environment sensor of each group; Determine whether there is a difference in the current data set collected by the GNSS station that meets the preset rules; If there is a current data set collected by a GNSS site that is different from the preset rule, it is determined that the current data set collected by the GNSS site that is consistent with the preset rule has current abnormal data.

[0007] Optionally, the determining whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data includes: A sliding T test method is used to determine whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data.

[0008] A second aspect of the present invention provides a mapping device for dynamic evaluation of GNSS site data integrity, comprising: An acquisition module is used to acquire the current data set collected by the meteorological environment sensors of N GNSS stations within a preset spatial range at the current moment; A grouping module, used to group the GNSS sites within a preset radius with each of the GNSS sites as the center of the circle to obtain N groups; A first judgment module is used to process the current data set collected by the meteorological environment sensor of each group to determine whether there is current abnormal data in the current data set of each group; The acquisition module is further configured to, when the judgment module determines that the current abnormal data exists, acquire data of the GNSS site having the current abnormal data at a preset time before the current moment; A second judgment module is used to judge whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data; A processing module, configured to determine, in the second judgment module, that the data of the GNSS site with the current abnormal data at a preset time before the current moment is abnormal data, discard the current abnormal data in the current data set, and determine a sensor data set to be determined, wherein the sensor data set to be determined is a data set after the current data set excludes the GNSS site data set with the current abnormal data; in the second judgment module, determine that the data of the GNSS site with the current abnormal data at a preset time before the current moment is not abnormal data, and determine the sensor data set to be determined, wherein the sensor data set to be determined is the current data set; A third judgment module is used to judge whether the GNSS sites corresponding to the sensor data set to be judged are greater than a preset ratio of the total GNSS sites; A drawing module is used to perform grid interpolation on the sensor data set to be judged, draw a spatial isosurface, and perform tile slicing when the third judgment module judges that the ratio is greater than the preset ratio.

[0009] Optionally, the processing module is further configured to notify the processing module to determine the sensor data set to be determined if the first determination module determines that there is no current abnormal data, wherein the sensor data set to be determined is the current data set.

[0010] Optionally, the first judgment module processes the current data set collected by the meteorological environment sensor of each group in the following manner to judge whether there is current abnormal data in the current data set of each group: Using the Z-Score function to process the current data set collected by the meteorological environment sensor of each group; Determine whether there is a difference in the current data set collected by the GNSS station that meets the preset rules; If there is a current data set collected by a GNSS site that is different from the preset rule, it is determined that the current data set collected by the GNSS site that is consistent with the preset rule has current abnormal data.

[0011] Optionally, the second judgment module judges whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data in the following manner: A sliding T test method is used to determine whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data.

[0012] A third aspect of the present invention provides an electronic device, comprising: a processor, a memory; The memory is used to store computer programs; The processor is used to execute the above-mentioned drawing method for dynamic evaluation of GNSS site data integrity by calling the computer program.

[0013] A fourth aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned drawing method for dynamic evaluation of GNSS site data integrity.

[0014] It can be seen that, through the drawing method, device, electronic device and computer-readable storage medium for dynamic evaluation of GNSS site data integrity provided by the present invention, abnormal data can be eliminated and the accuracy of GNSS site meteorological data drawing can be improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0016] Figure 1 A flowchart of a drawing method for dynamic evaluation of GNSS site data integrity provided by an embodiment of the present invention; Figure 2 A specific flow chart of a drawing method for dynamic evaluation of GNSS site data integrity provided by an embodiment of the present invention; Figure 3 A schematic diagram of a specific example of a drawing method for dynamic evaluation of GNSS site data integrity provided by an embodiment of the present invention; Figure 4 A schematic diagram of the structure of a drawing device for dynamic evaluation of GNSS site data integrity provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0017] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0018] Figure 1 A flowchart of a drawing method for dynamic evaluation of GNSS site data integrity provided by an embodiment of the present invention is shown. Figure 2 A specific flow chart of a drawing method for dynamic evaluation of GNSS site data integrity provided by an embodiment of the present invention is shown; see Figure 1 and Figure 2 The embodiment of the present invention provides a drawing method for dynamic evaluation of GNSS site data integrity, including: S1, obtaining the current data set collected by the meteorological environment sensors of N GNSS stations within a preset spatial range at the current moment.

[0019] Specifically, the present invention obtains a data set of all meteorological environment sensors within a certain spatial range at a certain time, and the data set may include, for example, temperature, humidity, atmospheric water content, gas propagation, etc. Wherein, N≥2.

[0020] S2, taking each GNSS station as the center of a circle, grouping the GNSS stations within a preset radius to obtain N groups.

[0021] Specifically, the present invention groups the sensors in spatial proximity based on the spatial dimension, for example, grouping the sensors with a radius of 10 km as the center of each site, wherein the collected values ​​of the gas phase environment sensors in the adjacent spaces are not much different.

[0022] S3, processing the current data set collected by each group of meteorological environment sensors, and determining whether there is current abnormal data in the current data set of each group.

[0023] As an optional implementation manner of an embodiment of the present invention, the current data set collected by each group of meteorological environment sensors is processed, and it is determined whether there is current abnormal data in the current data set of each group, including: using the Z-Score function to process the current data set collected by each group of meteorological environment sensors; determining whether there is a difference in the current data set collected by the GNSS station that meets the preset rules; if there is a difference in the current data set collected by the GNSS station that meets the preset rules, determining that the current data set collected by the GNSS station that meets the preset rules has current abnormal data.

[0024] Specifically, the present invention can use the Z-Score function to process each set of data sets to find out whether each set has abnormal data. The present invention finds out the data with differences between the sensors of the site and the surrounding sites through group calculation, and determines whether there are site data with large differences. Large differences can be specified according to actual needs. If there are sites with large differences, it can be determined that there are abnormal data.

[0025] S4, if there is current abnormal data, obtain data of the GNSS station with the current abnormal data at a preset time before the current moment.

[0026] Specifically, if the present invention finds abnormal data, the time series data of the site of the abnormal data in the previous period of time is obtained.

[0027] As an optional implementation of an embodiment of the present invention, the drawing method for dynamic evaluation of GNSS site data integrity provided by the embodiment of the present invention also includes: if there is no current abnormal data, determining a sensor data set to be judged, wherein the sensor data set to be judged is the current data set.

[0028] Specifically, if no abnormal data is found, the present invention performs a subsequent grid interpolation and drawing process.

[0029] S5, determining whether the data of the GNSS station with current abnormal data at a preset time before the current moment is abnormal data.

[0030] As an optional implementation party of an embodiment of the present invention, determining whether the data of a GNSS station with current abnormal data at a preset time before the current moment is abnormal data includes: determining whether the data of a GNSS station with current abnormal data at a preset time before the current moment is abnormal data through a sliding T test method.

[0031] S6, if the data of a GNSS site with current abnormal data at a preset time before the current moment is abnormal data, the current abnormal data in the current data set is discarded, and the sensor data set to be judged is determined, wherein the sensor data set to be judged is the data set after the current data set excludes the GNSS site data set of the current abnormal data; if the data of a GNSS site with current abnormal data at a preset time before the current moment is not abnormal data, then the sensor data set to be judged is determined, wherein the sensor data set to be judged is the current data set.

[0032] Specifically, the present invention verifies whether the data generates abnormal data by a sliding T test method. If it is abnormal data, the current abnormal data is discarded and does not participate in drawing. The purpose of excluding abnormal data is to prevent the interpolation process from using abnormal data, which affects the interpolation data in the local area. If it is not abnormal data, it can participate in drawing. If there is abnormal data, the grid interpolation of the present invention can make up for the current abnormal site data through the surrounding normal data. If there is no abnormal data, the data set of all meteorological environment sensors at a certain moment can be grid interpolated.

[0033] S7, determining whether the GNSS sites corresponding to the sensor data set to be determined are greater than a preset proportion of the total GNSS sites.

[0034] Specifically, before executing grid interpolation, the present invention also determines whether the site corresponding to the data set of the site sensor excluding abnormal data is a preset ratio of the total site. If it is greater than the preset ratio, the grid interpolation and drawing operation is performed, otherwise an error prompt message is thrown. Among them, the preset ratio in the present invention can be set to 1 / 3, of course, it can also be set to other ratios, which is not specifically limited in the present invention.

[0035] S8, if it is greater than the preset ratio, grid interpolation is performed on the sensor data set to be judged, and spatial isosurface drawing is performed to perform tile slicing.

[0036] Specifically, the present invention performs spatial isosurface drawing after grid interpolation, and finally performs tile slicing.

[0037] Below, with Figure 3 As an example, the specific implementation of the present invention is further described, but the present invention is not limited to this: After collecting the data set, it is determined that the current data of 5 points is quite different from that of other points, and there is a sudden change with the data of the previous moment of 5 points. In this case, it is necessary to draw a circle with 5 points as the center and a radius of 10km (larger than the station spacing, such as 1 or 2 times to ensure that other stations can be obtained) to obtain the current data of other stations (2, 3, 4, 8, 9 points) in the circle, and then determine whether the data of 5 points is available. If it is not available, it is discarded and it is considered that there is no data at 5 points at the current moment. When drawing the meteorological map, interpolate the point data around the missing points to avoid the influence of wrong values. The interpolation can include but is not limited to IDW (inverse distance weighted interpolation) and Kriging (Kriging interpolation).

[0038] It can be seen that the drawing method for dynamic evaluation of GNSS site data integrity provided by the embodiment of the present invention can eliminate abnormal data and improve the accuracy of GNSS site meteorological data drawing.

[0039] Figure 4The structure diagram of the drawing device for dynamic evaluation of GNSS site data integrity provided by an embodiment of the present invention is shown. The drawing device for dynamic evaluation of GNSS site data integrity applies the above method. The following is only a brief description of the structure of the drawing device for dynamic evaluation of GNSS site data integrity. For other matters not covered, please refer to the relevant description in the drawing method for dynamic evaluation of GNSS site data integrity. Figure 4 , a drawing device for dynamic evaluation of GNSS site data integrity provided by an embodiment of the present invention comprises: An acquisition module is used to acquire the current data set collected by the meteorological environment sensors of N GNSS stations within a preset spatial range at the current moment; A grouping module, used to group the GNSS stations within a preset radius with each GNSS station as the center of the circle to obtain N groups; The first judgment module is used to process the current data set collected by the meteorological environment sensor of each group, and judge whether there is current abnormal data in the current data set of each group; The acquisition module is further used to obtain data of the GNSS station with the current abnormal data at a preset time before the current moment when the judgment module determines that there is current abnormal data; The second judgment module is used to judge whether the data of the GNSS station with current abnormal data at a preset time before the current moment is abnormal data; A processing module, used for determining in the second judgment module that the data of the GNSS station with current abnormal data at a preset time before the current moment is abnormal data, discarding the current abnormal data in the current data set, and determining the sensor data set to be determined, wherein the sensor data set to be determined is the data set after the current data set excludes the GNSS station data set with the current abnormal data; determining in the second judgment module that the data of the GNSS station with current abnormal data at a preset time before the current moment is not abnormal data, and determining the sensor data set to be determined, wherein the sensor data set to be determined is the current data set; A third judgment module is used to judge whether the GNSS sites corresponding to the sensor data set to be judged are greater than a preset ratio of the total GNSS sites; The drawing module is used to perform grid interpolation on the sensor data set to be judged, draw a spatial isosurface, and perform tile slicing when the third judgment module judges that the ratio is greater than a preset ratio.

[0040] As an optional implementation of the embodiment of the present invention, the processing module is further used to notify the processing module to determine the sensor data set to be determined when the first judgment module determines that there is no current abnormal data, wherein the sensor data set to be determined is the current data set.

[0041] As an optional implementation of the embodiment of the present invention, the first judgment module processes the current data set collected by each group of meteorological environment sensors in the following manner to determine whether there is current abnormal data in the current data set of each group: Use the Z-Score function to process the current data set collected by each group of meteorological environment sensors; Determine whether there is a difference in the current data set collected by the GNSS station that meets the preset rules; If there is a difference with the current data set collected by the GNSS site that meets the preset rule, it is determined that the current data set collected by the GNSS site that meets the preset rule has current abnormal data.

[0042] As an optional implementation of the embodiment of the present invention, the second judgment module judges whether the data of the GNSS station with current abnormal data at a preset time before the current moment is abnormal data in the following manner: The sliding T test method is used to determine whether the data of the GNSS station with current abnormal data at a preset time before the current moment is abnormal data.

[0043] It can be seen that, through the GNSS site data integrity dynamic evaluation drawing device provided by the embodiment of the present invention, abnormal data can be eliminated and the accuracy of GNSS site meteorological data drawing can be improved.

[0044] The present invention also provides an electronic device, comprising: a processor, a memory; Memory for storing computer programs; The processor is used to execute the above-mentioned drawing method for dynamic evaluation of GNSS site data integrity by calling the computer program.

[0045] It can be seen that, through the electronic device provided by the embodiment of the present invention, abnormal data can be eliminated and the accuracy of GNSS site meteorological data mapping can be improved.

[0046] The present invention also provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the above-mentioned drawing method for dynamic evaluation of GNSS site data integrity is implemented.

[0047] It can be seen that, through the computer-readable storage medium provided by the embodiment of the present invention, abnormal data can be eliminated and the accuracy of GNSS site meteorological data mapping can be improved.

[0048] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.

Claims

1. A mapping method for dynamic evaluation of GNSS station data integrity, characterized in that: include: Obtain the current data set collected by the meteorological environment sensors of N GNSS stations within the preset spatial range at the current moment; Taking each of the GNSS stations as the center of a circle, grouping the GNSS stations within a preset radius to obtain N groups; Processing the current data set collected by the meteorological environment sensor of each group to determine whether there is current abnormal data in the current data set of each group; If the current abnormal data exists, obtaining data of the GNSS site having the current abnormal data at a preset time before the current moment; Determine whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data; If the data of the GNSS site with the current abnormal data at a preset time before the current moment is abnormal data, the current abnormal data in the current data set is discarded, and the sensor data set to be judged is determined, wherein the sensor data set to be judged is the data set after the current data set excludes the GNSS site data set of the current abnormal data; if the data of the GNSS site with the current abnormal data at a preset time before the current moment is not abnormal data, the sensor data set to be judged is determined, wherein the sensor data set to be judged is the current data set; Determine whether the GNSS site corresponding to the sensor data set to be determined is greater than a preset ratio of the total GNSS sites; If it is greater than the preset ratio, the sensor data set to be judged is subjected to grid interpolation, and spatial isosurface drawing is performed to perform tile slicing.

2. The method according to claim 1, characterized in that Also includes: If the current abnormal data does not exist, the sensor data set to be determined is determined, wherein the sensor data set to be determined is the current data set.

3. The method according to claim 2, characterized in that The processing of the current data set collected by the meteorological environment sensor of each group to determine whether there is current abnormal data in the current data set of each group includes: Using the Z-Score function to process the current data set collected by the meteorological environment sensor of each group; Determine whether there is a difference in the current data set collected by the GNSS station that meets the preset rules; If there is a current data set collected by a GNSS site that is different from the preset rule, it is determined that the current data set collected by the GNSS site that is consistent with the preset rule has current abnormal data.

4. The method according to claim 3, characterized in that The determining whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data comprises: A sliding T test method is used to determine whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data.

5. A mapping device for dynamic evaluation of GNSS station data integrity, characterized in that: include: An acquisition module is used to acquire the current data set collected by the meteorological environment sensors of N GNSS stations within a preset spatial range at the current moment; A grouping module, used to group the GNSS sites within a preset radius with each of the GNSS sites as the center of the circle to obtain N groups; A first judgment module is used to process the current data set collected by the meteorological environment sensor of each group to determine whether there is current abnormal data in the current data set of each group; The acquisition module is further configured to, when the judgment module determines that the current abnormal data exists, acquire data of the GNSS site having the current abnormal data at a preset time before the current moment; A second judgment module is used to judge whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data; A processing module, configured to determine, in the second judgment module, that the data of the GNSS site with the current abnormal data at a preset time before the current moment is abnormal data, discard the current abnormal data in the current data set, and determine a sensor data set to be determined, wherein the sensor data set to be determined is a data set after the current data set excludes the GNSS site data set with the current abnormal data; in the second judgment module, determine that the data of the GNSS site with the current abnormal data at a preset time before the current moment is not abnormal data, and determine the sensor data set to be determined, wherein the sensor data set to be determined is the current data set; A third judgment module is used to judge whether the GNSS sites corresponding to the sensor data set to be judged are greater than a preset ratio of the total GNSS sites; A drawing module is used to perform grid interpolation on the sensor data set to be judged, draw a spatial isosurface, and perform tile slicing when the third judgment module judges that the ratio is greater than the preset ratio.

6. The device according to claim 5, characterized in that The processing module is further configured to notify the processing module to determine the sensor data set to be determined if the first determination module determines that there is no current abnormal data, wherein the sensor data set to be determined is the current data set.

7. The device according to claim 6, characterized in that The first judgment module processes the current data set collected by the meteorological environment sensor of each group in the following manner to judge whether there is current abnormal data in the current data set of each group: Using the Z-Score function to process the current data set collected by the meteorological environment sensor of each group; Determine whether there is a difference in the current data set collected by the GNSS station that meets the preset rules; If there is a current data set collected by a GNSS site that is different from the preset rule, it is determined that the current data set collected by the GNSS site that is consistent with the preset rule has current abnormal data.

8. The device according to claim 7, characterized in that The second judgment module judges whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data in the following manner: The sliding T test method is used to determine whether the data of the GNSS station having the current abnormal data at a preset time before the current moment is abnormal data.

9. An electronic device, characterized in that: include: Processor, memory; The memory is used to store computer programs; The processor is configured to execute the drawing method for dynamic evaluation of GNSS site data integrity according to any one of claims 1 to 4 by calling the computer program.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the drawing method for dynamic evaluation of GNSS site data integrity as described in any one of claims 1 to 4 is implemented.

Images