CN112380310A - GNSS high-precision anti-aliasing calculation result smoothing method - Google Patents
GNSS high-precision anti-aliasing calculation result smoothing method Download PDFInfo
- Publication number
- CN112380310A CN112380310A CN202011348445.0A CN202011348445A CN112380310A CN 112380310 A CN112380310 A CN 112380310A CN 202011348445 A CN202011348445 A CN 202011348445A CN 112380310 A CN112380310 A CN 112380310A
- Authority
- CN
- China
- Prior art keywords
- smoothing
- calculation result
- result
- period
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 83
- 238000009499 grossing Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000001914 filtration Methods 0.000 claims description 11
- 238000013480 data collection Methods 0.000 claims description 2
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 238000013215 result calculation Methods 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 abstract description 22
- 230000008859 change Effects 0.000 abstract description 13
- 238000012544 monitoring process Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/215—Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
- Emergency Alarm Devices (AREA)
Abstract
The invention discloses a GNSS high-precision anti-sawtooth calculation result smoothing method which comprises the steps of determining time duration, collecting smoothing time period calculation result data, calculating a smoothing time period calculation result arithmetic mean, eliminating data with the maximum distance from the mean, and calculating final result values of all components by using a trend smoothing formula to output a smoothing result. By carrying out smoothing processing on the high-precision calculation result, the problem of peak protrusion of the calculation result is solved, the problem of saw-toothed shape is solved, and meanwhile, the displacement trend of the high-precision calculation result is kept, so that the displacement change value is closer to the real displacement value when alarm information is generated by using a threshold model, and the false alarm and the missing alarm are effectively reduced.
Description
Technical Field
The invention relates to the technical field of geological monitoring, in particular to a GNSS high-precision anti-aliasing calculation result smoothing method.
Background
In geological disaster and major engineering project deformation monitoring projects, GNSS high-precision resolving software gives displacement change values of a monitoring station relative to a reference station at short time intervals. The relative displacement variation value error influence factors are as follows: A. station distance: increasing with increasing distance, B, station environment: for example, the shielding condition around the measuring station, the multipath influence of the measuring station and the like are combined, most of positioning errors are about 2-5 mm, and part of conditions can only reach the centimeter-level precision. When the monitored object has no displacement change, the displacement change value is displayed in a visual chart with the horizontal axis of time and the vertical axis of the displacement change value, and the displacement change value curve is jagged and occasionally has peak protrusions.
The current calculation software generally directly outputs the calculation result or takes the average value of the latest N calculation results as the calculation result of the current time period, so that the wave crest and the wave trough of the saw-toothed form are reduced to a certain extent, but the existence of the saw-toothed form cannot be completely eliminated. The existence of the jagged shape and the peak protrusions brings certain false alarm and missed alarm for alarm based on threshold value, which is not beneficial to the realization of alarm strategy and the group detection and group defense work of monitoring units.
Therefore, how to solve the problem of peak protrusion and the problem of zigzag morphology of the calculation result, and meanwhile, the method retains the displacement trend of the high-precision calculation result and reduces the false alarm and the missing report, and is a technical problem which needs to be solved urgently.
The above is only for the purpose of assisting understanding of the technical aspects of the present invention, and does not represent an admission that the above is prior art.
Disclosure of Invention
The invention mainly aims to provide a GNSS high-precision anti-sawtooth calculation result smoothing method, and aims to solve the technical problems that sawtooth shapes cannot be completely eliminated, and certain false alarm and false alarm are caused in the prior art.
In order to achieve the above object, a first aspect of the present invention provides a method for smoothing a GNSS high-precision anti-aliasing solution result, where the method for smoothing the high-precision anti-aliasing solution result includes the following steps:
determining time period duration comprising a single epoch calculation time period, a calculation result smooth output time period and a data smooth time period;
collecting smooth time interval calculation result data;
calculating the arithmetic mean of the calculation results in the smoothing time period, and calculating the distance between each calculation result and the arithmetic mean of the results;
sorting according to the distance between the calculation result and the arithmetic mean value of the result, and rejecting 5% of data with the maximum distance from the mean value;
removing special cases from the output time interval of the calculation result, and calculating the average value of each component of the filtered result according to the output time interval of the calculation result;
and (5) solving a final result value of each component by using a trend smoothing formula, finishing smoothing and outputting a result.
Preferably, the time duration of the determined time period is determined according to a rule that 1 smooth output time period comprises 10 or more single-epoch calculation time periods, and 1 data smooth time period comprises 60-100 smooth output time periods.
Preferably, the high-precision anti-aliasing calculation result smoothing method uses the mean value position vector in the smoothing time period as a reference, calculates the distances between all calculation results and the mean value, sequences the obtained distances, and eliminates the distances from large to small according to a percentage threshold so as to remove unstable values.
Preferably, the method for smoothing the high-precision anti-aliasing calculation result generates a special case during filtering, namely, for the calculation value in the last result output period, when the distance percentage threshold rejection rate in the output period reaches more than 30%, the rejection is not carried out.
Preferably, the high-precision anti-aliasing calculation result smoothing method calculates the average value of the filtered position vector data in X, Y, Z three directions according to the calculation result output time periods, and each calculation result output time period obtains an average position vector data; and then performing trend smoothing on the group of position vector data, wherein a smoothing formula is as follows:
wherein n is the number of output intervals of the calculation result in the smoothing period, x0、x1……xn-1The x-direction average of the interval is output for each solution from new to old, respectively.
In a second aspect of the present invention, a GNSS high-precision antialiasing calculation result smoothing system is provided, where the high-precision antialiasing calculation result smoothing system includes:
a data collection unit: collecting smooth time interval calculation result data;
a smoothing period calculation result calculation unit: calculating the arithmetic mean of the calculation results in the smoothing time period, calculating the distance between each calculation result and the arithmetic mean of the results, and sequencing the distances between the calculation results and the arithmetic mean of the results;
a data filtering unit: eliminating 5% of data with the maximum distance from the mean value, and performing special case elimination processing on the current resolving result output time period;
a trend smoothing unit: and (4) solving the average value of each component according to the output time period of the calculation result after filtering, solving the final result value of each component by using a trend smoothing formula, and outputting a smooth result.
According to the invention, the problem of peak protrusion of the calculation result is solved and the problem of saw-toothed shape is solved by smoothing the high-precision calculation result, and the displacement trend of the high-precision calculation result is kept, so that the displacement change value is closer to the real displacement value when the alarm information is generated by using the threshold model, and the false alarm and the missing alarm are effectively reduced.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the structures shown in the drawings without creative efforts.
FIG. 1 is a schematic diagram illustrating a step principle of a GNSS high-precision anti-aliasing solution result smoothing method provided by the present invention;
FIG. 2 is a schematic diagram of a time interval relationship proposed by the present invention;
FIG. 3 is a schematic representation of the present invention before it is smoothly practiced;
FIG. 4 is a schematic diagram showing the present invention after being smoothly implemented;
FIG. 5 is a monitoring diagram of the landslide process after smoothing according to the present invention.
The implementation, functional features and advantages of the objects of the present invention will be further explained with reference to the accompanying drawings.
Detailed Description
It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention provides an embodiment, and referring to fig. 1, fig. 1 is a schematic diagram illustrating a step principle of a GNSS high-precision anti-aliasing solution result smoothing method provided by the invention.
As shown in fig. 1, in this embodiment, a method for smoothing a high-precision antialiasing calculation result of a GNSS includes the following steps:
determining time period duration comprising a single epoch calculation time period, a calculation result smooth output time period and a data smooth time period;
collecting smooth time interval calculation result data;
calculating the arithmetic mean of the calculation results in the smoothing time period, and calculating the distance between each calculation result and the arithmetic mean of the results;
sorting according to the distance between the calculation result and the arithmetic mean value of the result, and rejecting 5% of data with the maximum distance from the mean value;
removing special cases from the output time interval of the calculation result, and calculating the average value of each component of the filtered result according to the output time interval of the calculation result;
and (5) solving a final result value of each component by using a trend smoothing formula, finishing smoothing and outputting a result.
Further, three time periods need to be determined in advance before smoothing of the high-precision calculation result, as shown in fig. 2, fig. 2 is a time period relationship diagram provided by the present invention, and includes a single epoch calculation time period, a calculation result smooth output time period, and a data smooth time period; the three time interval durations are configured according to the specific situation of the project, generally 1 smooth output time interval includes more than 10 single epoch solution time intervals, and the data smooth time interval is 60 to 100 smooth output time intervals.
Furthermore, the main technical key points of the smoothing of the high-precision calculation result are filtering and trend smoothing, and through data processing in the two processes, the purposes of outputting the high-precision calculation result without saw teeth or with few saw teeth and retaining the displacement change trend are achieved.
Further, the high-precision anti-aliasing calculation result smoothing method uses the mean value position vector in the smoothing time period as a reference, calculates the distances between all calculation results and the mean value, sequences the obtained distances, and eliminates the distances from large to small according to the percentage threshold so as to remove unstable values.
Furthermore, for the items with light multipath effect, not serious shielding and short baseline, the filtering percentage can be below 5 percent, even filtering is not needed; on the contrary, the filtration percentage is 5 to 10 percent according to specific items.
Further, there is a special case during filtering, that is, for a solution value in a last result output period, if the rejection rate reaches more than 30% according to the distance percentage threshold in the period, the solution value is preferably not rejected, which may be true displacement, and needs to be preserved to ensure the sensitivity of deformation monitoring.
Further, the high-precision anti-aliasing calculation result smoothing method calculates the average value of the filtered position vector data in X, Y, Z three directions according to the calculation result output time periods, and each calculation result output time period obtains an average position vector data; and then performing trend smoothing on the group of position vector data, wherein a smoothing formula is as follows:
wherein n is the number of output intervals of the calculation result in the smoothing period, x0、x1……xn-1The x-direction average of the interval is output for each solution from new to old, respectively. In the smoothing process, the latest calculation value occupies a larger weight, and the current displacement change trend is more easily reflected, so that the displacement change of the monitored object can be more sensitively reflected.
Specifically, the smoothing method of the present invention is used in a certain geological disaster slip monitoring project, fig. 3 and 4 are respectively a comparison graph before and after smoothing, and fig. 5 is a monitoring graph of a slip process. The smooth process is explained to keep the displacement change trend, the displacement change can be quickly reflected, meanwhile, the chart is smoother, and the client is not disturbed by the false report of the jagged data.
The method and the device perform smoothing processing on the high-precision calculation result, solve the problem of peak protrusion of the calculation result, solve the problem of saw-toothed shape, and simultaneously keep the displacement trend of the high-precision calculation result, so that the displacement change value is closer to the real displacement value when the alarm information is generated by using a threshold model, and the false alarm and the missing alarm are effectively reduced.
The foregoing is illustrative of the preferred embodiments of this invention, and it is to be understood that the invention is not limited to the precise form disclosed herein and that various other combinations, modifications, and environments may be resorted to, falling within the scope of the concept as disclosed herein, either as described above or as apparent to those skilled in the relevant art. And that modifications and variations may be effected by those skilled in the art without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (6)
1. A GNSS high-precision anti-aliasing calculation result smoothing method is characterized by comprising the following steps:
determining time period duration comprising a single epoch calculation time period, a calculation result smooth output time period and a data smooth time period;
collecting smooth time interval calculation result data;
calculating the arithmetic mean of the calculation results in the smoothing time period, and calculating the distance between each calculation result and the arithmetic mean of the results;
sorting according to the distance between the calculation result and the arithmetic mean value of the result, and eliminating the data with the maximum distance from the mean value;
removing special cases from the output time interval of the calculation result, and calculating the average value of each component of the filtered result according to the output time interval of the calculation result;
and (5) solving a final result value of each component by using a trend smoothing formula, finishing smoothing and outputting a result.
2. The method of claim 1, wherein the duration of the determined period comprises 10 or more single epoch solution periods according to 1 smooth output period, and 1 data smoothing period comprises a rule of 60-100 smooth output periods, and the single epoch solution period, the solution result smooth output period, and the data smoothing period are determined.
3. The method of claim 1, wherein the method of smoothing the high-precision anti-aliasing solution results of the GNSS uses a mean position vector in a smoothing period as a reference, calculates distances between all solution results and a mean, sorts the obtained distances, and eliminates the distances from large to small according to a percentage threshold to remove unstable values.
4. The method of claim 1, wherein the method of smoothing the high-accuracy anti-aliasing solution results of GNSS is characterized in that a special case is generated during filtering, that is, for a solution value in a last result output period, when a distance percentage threshold rejection rate in the output period reaches more than 30%, the solution value is not rejected.
5. The method of claim 1, wherein the method of smoothing the high-accuracy anti-aliasing solution results averages filtered position vector data in X, Y, Z three directions according to solution result output time periods, and each solution result output time period obtains an averaged position vector data; and then performing trend smoothing on the group of position vector data, wherein a smoothing formula is as follows:
wherein n is the number of output intervals of the calculation result in the smoothing period, x0、x1……xn-1The x-direction average of the interval is output for each solution from new to old, respectively.
6. A GNSS high-precision anti-aliasing solution smoothing system, comprising:
a data collection unit: collecting smooth time interval calculation result data;
a smoothing period calculation result calculation unit: calculating the arithmetic mean of the calculation results in the smoothing time period, calculating the distance between each calculation result and the arithmetic mean of the results, and sequencing the distances between the calculation results and the arithmetic mean of the results;
a data filtering unit: eliminating 5% of data with the maximum distance from the mean value, and performing special case elimination processing on the current resolving result output time period;
a trend smoothing unit: and (4) solving the average value of each component according to the output time period of the calculation result after filtering, solving the final result value of each component by using a trend smoothing formula, and outputting a smooth result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011348445.0A CN112380310B (en) | 2020-11-26 | 2020-11-26 | GNSS high-precision anti-aliasing resolving result smoothing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011348445.0A CN112380310B (en) | 2020-11-26 | 2020-11-26 | GNSS high-precision anti-aliasing resolving result smoothing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112380310A true CN112380310A (en) | 2021-02-19 |
| CN112380310B CN112380310B (en) | 2023-12-01 |
Family
ID=74588696
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011348445.0A Active CN112380310B (en) | 2020-11-26 | 2020-11-26 | GNSS high-precision anti-aliasing resolving result smoothing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112380310B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101794345A (en) * | 2009-12-30 | 2010-08-04 | 北京世纪高通科技有限公司 | Data processing method and device |
| CN103792558A (en) * | 2014-01-10 | 2014-05-14 | 中国人民解放军63921部队 | GNSS carrier phase smoothness pseudo-range processing method |
| CN105425254A (en) * | 2015-12-28 | 2016-03-23 | 西北工业大学 | Dynamic GNSS measurement data anti-outlier bidirectional smoothing filtering method |
| JP2016065742A (en) * | 2014-09-24 | 2016-04-28 | 株式会社トプコン | Vibration measurement method, vibration detector, and vibration measurement program |
| CN106443735A (en) * | 2016-11-01 | 2017-02-22 | 国网新疆电力公司信息通信公司 | Beidou high-precision measurement method based on comprehensive analysis method |
| CN107798871A (en) * | 2017-10-27 | 2018-03-13 | 云南大学 | A kind of freeway toll station traffic flow forecasting method and system |
-
2020
- 2020-11-26 CN CN202011348445.0A patent/CN112380310B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101794345A (en) * | 2009-12-30 | 2010-08-04 | 北京世纪高通科技有限公司 | Data processing method and device |
| CN103792558A (en) * | 2014-01-10 | 2014-05-14 | 中国人民解放军63921部队 | GNSS carrier phase smoothness pseudo-range processing method |
| JP2016065742A (en) * | 2014-09-24 | 2016-04-28 | 株式会社トプコン | Vibration measurement method, vibration detector, and vibration measurement program |
| CN105425254A (en) * | 2015-12-28 | 2016-03-23 | 西北工业大学 | Dynamic GNSS measurement data anti-outlier bidirectional smoothing filtering method |
| CN106443735A (en) * | 2016-11-01 | 2017-02-22 | 国网新疆电力公司信息通信公司 | Beidou high-precision measurement method based on comprehensive analysis method |
| CN107798871A (en) * | 2017-10-27 | 2018-03-13 | 云南大学 | A kind of freeway toll station traffic flow forecasting method and system |
Non-Patent Citations (1)
| Title |
|---|
| 高宝: "基于指数平滑模型的学习能力变化趋势评价方法研究", 《中国林业教育》, vol. 36, no. 5, pages 1 - 4 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112380310B (en) | 2023-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105808759B (en) | Data preprocessing method and system for track mining | |
| CN110286656B (en) | False alarm filtering method and device for tolerance of error data | |
| CN110648480B (en) | Single variable alarm system and method based on change rate | |
| CN112067015B (en) | Step counting method and device based on convolutional neural network and readable storage medium | |
| CN117439827B (en) | Network flow big data analysis method | |
| CN113237432B (en) | Laser welding penetration extraction method and device | |
| CN109087490B (en) | A kind of interference alarm recognition methods based on alarm duration features | |
| CN110889597A (en) | Method and device for detecting abnormal business timing sequence indexes | |
| CN112016206A (en) | Method and system for judging instability state of tower, computer equipment and application | |
| CN116880356A (en) | Method and device for monitoring machining state of numerical control machine tool | |
| CN112332807A (en) | Weak envelope signal detection method and system | |
| CN112380310A (en) | GNSS high-precision anti-aliasing calculation result smoothing method | |
| CN112781556A (en) | Well lid transaction monitoring method and device based on multi-data fusion filtering | |
| CN103020907B (en) | Based on the DSPI striped filtering system of two-dimensional ensemble empirical mode decomposition | |
| CN114722866B (en) | Management and control system is tracked to equipment based on industry internet | |
| CN114879033A (en) | Motor fault detection method, system and storage medium | |
| CN116123459A (en) | Judgment method for pressure fluctuation of heat supply pipe network | |
| CN113465826B (en) | Gas leakage detection method and device | |
| CN115424075A (en) | Method and system for monitoring pipeline state | |
| CN110631650B (en) | Data cleaning method and device based on time series data self-increment characteristics | |
| CN110134069B (en) | Self-diagnosis method and system of numerical control machine tool | |
| CN114021421B (en) | Wireless sensor availability evaluation method and device based on hidden Markov model | |
| CN112061080B (en) | Vehicle abnormal motion detection method, device, equipment and medium | |
| He et al. | On the estimation error in zero-inflated Poisson model for process control | |
| CN113450000B (en) | Index data anomaly detection method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |