[go: up one dir, main page]

CN114236581B - Beidou slope monitoring data post-processing method - Google Patents

Beidou slope monitoring data post-processing method Download PDF

Info

Publication number
CN114236581B
CN114236581B CN202210188780.1A CN202210188780A CN114236581B CN 114236581 B CN114236581 B CN 114236581B CN 202210188780 A CN202210188780 A CN 202210188780A CN 114236581 B CN114236581 B CN 114236581B
Authority
CN
China
Prior art keywords
result
beidou
post
data
historical
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.)
Active
Application number
CN202210188780.1A
Other languages
Chinese (zh)
Other versions
CN114236581A (en
Inventor
梁晓东
雷孟飞
汤金毅
龙兴
黎凯
张涛
周俊华
刘琴
刘�文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan Lianzhi Technology Co Ltd
Original Assignee
Hunan Lianzhi Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan Lianzhi Technology Co Ltd filed Critical Hunan Lianzhi Technology Co Ltd
Priority to CN202210188780.1A priority Critical patent/CN114236581B/en
Publication of CN114236581A publication Critical patent/CN114236581A/en
Application granted granted Critical
Publication of CN114236581B publication Critical patent/CN114236581B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/35Constructional details or hardware or software details of the signal processing chain
    • G01S19/37Hardware or software details of the signal processing chain
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Remote Sensing (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Signal Processing (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Algebra (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention solves the problems of low precision and high delay of the Beidou data post-processing method in the prior art. The invention provides a Beidou slope monitoring data post-processing method, which realizes the resolving of Beidou original monitoring data, carries out gross error judgment and periodic error correction on the current resolving result according to the historical monitoring result of a monitoring point, and effectively corrects the periodic error in the monitoring result; the method also comprises a data smoothing algorithm, wherein the data smoothing algorithm is a novel self-adaptive smoothing algorithm, on one hand, high-frequency errors can be removed, on the other hand, the displacement of the current result can be reserved, and high-precision low-delay post-processing calculation is realized.

Description

Beidou slope monitoring data post-processing method
Technical Field
The invention relates to the field of big data, in particular to a Beidou side slope monitoring data post-processing method.
Background
Along with the construction of the Beidou system, the Beidou system is greatly popularized in side slope monitoring by the characteristics of high precision, all weather, simplicity and convenience in operation and the like.
Most errors in the Beidou data resolving process are eliminated, however, in slope monitoring projects, signal shielding is serious, satellite constellation distribution is uneven, long period errors and high frequency errors are obvious and cannot be effectively eliminated, the periodic errors and the high frequency errors bring great difficulty to slope monitoring, technicians are difficult to extract useful signals from all the errors to judge the state of a slope, the deformation of the slope is slow deformation under the common condition, sudden landslide may occur under extreme weather such as rainstorm, typhoon or earthquake, and therefore slope monitoring needs to be provided with both slow deformation high-precision monitoring and sudden deformation sudden monitoring.
At present, the Beidou calculation in slope monitoring mainly adopts real-time calculation and post-processing calculation, wherein the real-time calculation adopts a dynamic rtk algorithm and can obtain a monitoring result according to the single epoch data calculation of a receiver; the post-processing calculation carries out unified calculation on data of the receiver for a period of time to obtain the base length from the monitoring point to a reference station, then the adjustment algorithm is utilized to carry out adjustment processing on the data of the monitoring point to obtain a high-precision displacement result of each monitoring point, and the real-time dynamic rtk in the prior art has low calculation precision and stability and can not meet the requirements of slope monitoring and early warning; the post-processing calculation precision is high, but the traditional post-processing calculation delay is high, and sudden slope sliding cannot be quickly reflected.
Disclosure of Invention
In order to solve the above problems, the invention provides a post-processing method of a slope monitoring result aiming at the characteristics of a slope, which comprises the following steps:
step S1: acquiring original data, and acquiring observation data and ephemeris data of a monitoring point and a reference station;
step S2: resolving original data, specifically resolving the original data to obtain a time period result;
step S3: screening data, specifically, eliminating gross error results from the time interval results to obtain processed time interval results;
step S4: averaging the processed time period results to obtain an average result;
step S5: periodic error correction, specifically, performing periodic error correction on the average result to obtain a corrected periodic result, and taking the corrected periodic result as a historical result;
step S6: a self-adaptive smoothing algorithm, specifically a dynamic smoothing algorithm or a stable smoothing algorithm is configured in configuration information;
step S7: and (4) performing post-processing, specifically, calculating by using a dynamic smoothing algorithm or a stable smoothing algorithm to obtain final data.
Note that the monitoring point post-processing time interval is set to h hours.
Preferably, in step S2, the raw data solution employs a conventional baseline solution algorithm.
Preferably, in step S3, the coarse error result is eliminated by using a triple error method for the segment result.
Preferably, the latest D days of the period result of the point is queried in step S5, and there are D × 24/h results. The accuracy of error correction is higher along with the increase of the number of days of the time interval result, but the accuracy improvement is limited after the number of days is more than 7 days, and (D is more than or equal to 2 and less than or equal to 7) is taken in practical application in consideration of the calculation complexity, and if (D is more than or equal to 0 and less than 2), periodic error correction is not carried out.
Specifically, the periodic error correction method in step S5 is as follows:
subtracting the average value of the time period result from the time period result, segmenting the time period result according to days, dividing the time period result into D segments on the assumption that D days are totally used as the time period result, respectively averaging the results of the corresponding time periods in the D segment historical result, using the average value as the periodic error correction value of the time period to totally obtain 24/h error values, searching the corresponding error values according to the time of the average value result, and subtracting the periodic error values from the current average value result to obtain the corrected time period result.
Further, the configuration information configures dynamic smoothing calculation or stable smoothing calculation according to the application scenario.
In particular, dynamic smoothing is used in monitoring scenarios where the displacement of the monitored object needs to be reflected quickly, such as: scenes such as monitoring of a sliding slope, monitoring of bridge deviation rectification and the like are generated; the stable mode is suitable for long-term monitoring of a more stable slope.
Preferably, the dynamic smoothing calculation uses a kalman filter algorithm.
Preferably, stable smooth calculation introduces the historical result stability index into the final post-processing result calculation, and the stability of the monitoring result is improved.
Specifically, the stable smoothing algorithm needs to calculate a historical result stationarity index, the historical result stationarity index needs to be calculated by a slope ratio, a median error ratio and a counter ratio, wherein the slope ratio represents stationarity of a historical monitoring result, the median error ratio represents stability degree of a current result, the counter ratio represents reliability of the current abnormal result, and the larger the counter ratio is, the larger the reliability of the current result is.
Further, the historical result stationarity index calculation method comprises the following steps:
Figure GDA0003587251600000031
wherein theta is a stability index of a historical result; k is the slope of the historical result; k M Is a slope threshold, K M 10/(24/h); p is the current time period result; m is a historical result reference value; STD is the error in the historical result; f is the weight of the error ratio; n is the current of the counterA value; n is a radical of max Is the maximum value of the counter; and c is the counter weight.
In the above formula
Figure GDA0003587251600000032
When c is equal to 0, the compound is,
Figure GDA0003587251600000033
when c is 1; when | P-M |>When 16 STD, f is 0, when | P-M ≦ 16 STD, f is 1; when theta is more than or equal to 1, theta is equal to 1;
wherein the historical result reference value M is the average value of the latest 24-hour data in the historical results of D days.
And calculating the slope K of the historical result by using a least square method, wherein the slope K is used for judging whether the historical monitoring result is stable or not.
And the error STD in the historical result is the error in the historical result of D days, and is used for judging whether the result of the current time period is mutated or not.
Wherein, 3 STD < | X-M | in the above formula represents that P is an abnormal result; the value taking method for dividing the abnormal result into the current value n of the counter under three conditions according to the difference value of P and M is as follows:
Figure GDA0003587251600000034
further, the stable smoothing calculation post-processing result calculation formula is as follows:
P'=P*θ+M*(1-θ);
where P' is the final data.
The technical scheme of the invention has the following beneficial effects: the method realizes the resolving of the original Beidou monitoring data, performs gross error judgment and periodic error correction on the current resolving result according to the historical monitoring result of the monitoring point, and effectively corrects the periodic error in the monitoring result; the method also comprises a data smoothing algorithm which is a novel self-adaptive smoothing algorithm, so that on one hand, high-frequency errors can be removed, on the other hand, the displacement of the current result can be kept, and high-precision low-delay post-processing calculation is realized.
In addition to the objects, features and advantages described above, other objects, features and advantages of the present invention are also provided. The invention will now be described in further detail with reference to the accompanying drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the invention and, together with the description, serve to explain the invention and not to limit the invention. In the drawings:
FIG. 1 is a flow chart of a method for post-processing slope monitoring results;
FIG. 2 is a graph of a periodic error correction value model in the X direction;
FIG. 3 is a graph comparing the error in the X-direction post-processing results with the original results.
Detailed Description
Embodiments of the invention will be described in detail below with reference to the drawings, but the invention can be implemented in many different ways, which are defined and covered by the claims.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs, and the terms used herein in the specification of the present invention are for the purpose of describing particular embodiments only and are not intended to limit the present invention.
The invention provides a Beidou side slope monitoring data post-processing method, which realizes high-precision low-delay post-processing of Beidou original monitoring data, and comprises the following specific implementation modes:
as shown in fig. 1, a Beidou side slope monitoring data post-processing method includes the following steps:
step S1: acquiring original data, and acquiring observation data and ephemeris data of a monitoring point and a reference station;
step S2: resolving original data, specifically resolving the original data to obtain a time interval result;
step S3: screening data, specifically, rejecting gross error results from the time interval results to obtain processed time interval results;
step S4: calculating an average result, specifically calculating an average value of the processed time period result to obtain an average result;
step S5: performing periodic error correction on the average result, specifically performing periodic error correction on the average result to obtain a corrected periodic result as a historical result;
step S6: and the adaptive smoothing algorithm is specifically a dynamic smoothing algorithm or a stable smoothing algorithm configured in the configuration information.
Step S7: and (4) performing post-processing, specifically, calculating by using a dynamic smoothing algorithm or a stable smoothing algorithm to obtain final data.
Note that the monitoring point post-processing time interval is set to h hours.
Preferably, in step S2, the raw data is solved using a conventional baseline solution algorithm.
Preferably, in step S3, the gross error result is removed using a triple medium error method on the segment result.
Preferably, the latest D days of the period result of the point is queried in step S5, and there are a total of D × 24/h results. The accuracy of error correction is higher along with the increase of the number of days of the time interval result, but the accuracy improvement is limited after the number of days is more than 7 days, and (D is more than or equal to 2 and less than or equal to 7) is taken in practical application in consideration of the calculation complexity, and if (D is more than or equal to 0 and less than 2), periodic error correction is not carried out.
Specifically, the periodic error correction method in step S5 is as follows:
subtracting the average value of the period result from the period result, segmenting according to the number of days, dividing the period result into D segments on the assumption that the period result has D days, respectively averaging the results of the corresponding period in the D segment historical results, taking the average value as the periodic error correction value of the period to obtain 24/h error values in total, searching the corresponding error value according to the time of the average value result, and subtracting the periodic error value from the current average value result to obtain the corrected period result.
Further, the configuration information configures dynamic smoothing calculation or stable smoothing calculation according to the application scenario.
In particular, dynamic smoothing is used in monitoring scenarios where the displacement of the monitored object needs to be reflected quickly, such as: scenes such as monitoring of a sliding slope, monitoring of bridge deviation rectification and the like are generated; the stabilization mode is suitable for long-term monitoring of more stable slopes.
Preferably, the dynamic smoothing algorithm is a kalman filter algorithm.
Preferably, the stability smoothing algorithm introduces the stability index of the historical result into a post-processing method to calculate final data, and the stability of the monitoring result is improved.
It should be noted that the stable smoothing algorithm needs to calculate a historical result stationarity index, the historical result stationarity index needs to be calculated by a slope ratio, a median error ratio and a counter ratio, wherein the slope ratio represents stationarity of a historical monitoring result, the median error ratio represents stability degree of a current result, the counter ratio represents reliability of a current abnormal result, and the larger the counter ratio is, the larger the reliability of the current result is.
Further, the historical result stationarity index calculation method comprises the following steps:
Figure GDA0003587251600000051
wherein K is the slope of the historical result; k M Is a slope threshold, K M 10/(24/h); p is the current time interval result; n is a radical of hydrogen max Is the maximum value of the counter; n is the current value of the counter; f is the weight of the error ratio; c is the counter weight; theta is a stability index of a historical result; m is a historical result reference value; STD is the error in the historical result.
In the above formula, when
Figure GDA0003587251600000061
When c is 0; when the temperature is higher than the set temperature
Figure GDA0003587251600000062
When c is 1; when | P-M |>When 16 STD, f is 0, when | P-M ≦ 16 STD, f is 1; theta is more than or equal to 1,θ=1;
Wherein the historical result reference value M is the average value of the latest 24-hour data in the historical results of D days.
Wherein the historical result slope is calculated using a least squares method.
Wherein, STD is the error in the D calendar history result.
The value method of the current value n of the counter is as follows:
Figure GDA0003587251600000063
in the above formula, 3 × STD < | P-M | indicates that P is an abnormal result, and the abnormal result is divided into three cases according to the difference between X and M, each case having a different influence on n.
Further, the final data calculation formula of the post-processing method is as follows:
P'=P*θ+M*(1-θ);
the above scheme was applied to the following specific experimental cases:
the data from 2021-09-15 to 2021-09-24 were processed with a post-processing time interval of 0.5 hours.
1. Acquiring observation data and ephemeris data of a monitoring point and a reference station for 30 minutes as original data;
2. resolving the Beidou data by using a conventional baseline resolving algorithm to obtain a 30-minute time period result set of the monitoring points;
3. removing gross error results in a 30-minute time period result set by using a triple error method;
4. averaging the result set in the 30-minute time period to obtain an average result;
5. and carrying out periodic error correction on the time interval result after the gross error is eliminated, wherein the specific correction method comprises the following steps:
(1) the historical results of the last 5 days of the point are queried, and the total number of results is 240.
(2) Subtracting the average value of the time interval result from each time interval result, segmenting the time interval result according to the number of days, dividing the time interval result into 5 segments in total, respectively averaging the time interval results of the corresponding time points in the 5 segments of time interval results, taking the average value as the periodic error correction value of the time point, and obtaining 48 error values in total, wherein the model curve of the periodic error correction value in the X direction is shown in fig. 2.
(3) The current time period results in (120.1,115.3,134.8), the result time is 2021-10-1000:30: 00; the horizontal axis of the error model value corresponding to the time is 2, and the periodic error model value is found to be (5.2, -6.9, 9.7); if the error model value is subtracted from the current result, the corrected result is P ═ 114.9, 122.2, 125.1.
6. And (3) calculating a stationarity index:
calculate the average of the last 24 hours history: m ═ (116.7, 125.5, 121.6) as a reference;
calculate median error for 5 days of historical results: STD ═ (1.3, 1.5, 2.7);
calculating the slope of the curve of the 5-day history result by using a least square method, wherein the slopes in three directions are K (0.009,0.015 and 0.012), and the slope threshold is 0.1, N max Taking 3, | X-M | (1.8,3.3,3.5), | X-M | ≦ 3 ≦ STD, so f ═ 1, c ═ 0;
Figure GDA0003587251600000071
7. and (3) calculating a post-processing result:
P'=P*θ+M*(1-θ)=(116.4,124.6,122.3)。
as shown in fig. 3, a graph comparing the error between the post-processing result in the X direction and the original result, where the solid line is the post-processing result of this embodiment and the dotted line is the error in the original result, it can be seen from the graph that the periodicity of the error in the original result is obvious, the processed result is stable, there is no obvious periodicity, and both the precision and the stability are improved.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (8)

1. A Beidou side slope monitoring data post-processing method is characterized by comprising the following steps:
step S1: acquiring original data, and acquiring observation data and ephemeris data of a monitoring point and a reference station;
step S2: resolving original data, specifically, resolving the original data to obtain a time period result;
step S3: screening data, specifically, removing gross error results from the time interval results to obtain processed time interval results;
step S4: averaging the processed time period results to obtain an average result;
step S5: periodic error correction, specifically, periodic error correction is carried out on the average result to obtain a corrected periodic result, and the corrected periodic result is used as a historical result;
step S6: an adaptive smoothing algorithm, specifically a stable smoothing algorithm is configured in configuration information; the stability smoothing algorithm introduces stability indexes into post-processing result calculation, and the historical result stability index calculation formula is as follows:
Figure FDA0003594005050000011
wherein theta is a stability index of a historical result; k is the slope of the historical result; k M Is a slope threshold; p is the current time period result; m is a historical result reference value; STD is the error in the historical result; f is the weight of the error ratio; n is the current value of the counter; n is a radical of hydrogen max Is the maximum value of the counter; c is the counter weight;
the value method of the current value n of the counter is as follows:
Figure FDA0003594005050000012
the counter counts abnormal values in the historical result;
step S7: and (4) performing post-processing, specifically, calculating by using a stable smoothing algorithm to obtain final data.
2. The Beidou side slope monitoring data post-processing method according to claim 1, characterized in that in the step S2, a conventional baseline solution algorithm is adopted for original data solution.
3. The Beidou slope monitoring data post-processing method according to claim 1, wherein in step S3, a triple error method is adopted for time interval result coarse difference elimination results.
4. The Beidou slope monitoring data post-processing method according to claim 1, wherein in the step S5, the periodic error correction method is as follows:
and subtracting the average value of the time period result from the time period result, segmenting according to the number of days, respectively averaging the results of the corresponding time periods in the time period result of the current day to be used as the periodic error correction value of the time period, subtracting the periodic error value from the average value of the current time period result to obtain a corrected time period result, and using the corrected time period result as a historical result.
5. The Beidou slope monitoring data post-processing method according to claim 1, wherein in the step S6, the configuration information is further configured with a dynamic smoothing algorithm according to an application scenario;
the dynamic smoothing algorithm is used for monitoring scenes needing to rapidly reflect the displacement of the monitored object;
the stability smoothing algorithm is used for long-term monitoring of relatively stable monitored objects.
6. The Beidou side slope monitoring data post-processing method according to claim 5, characterized in that a Kalman filtering algorithm is adopted for the dynamic smoothing algorithm to output results.
7. The Beidou side slope monitoring data post-processing method according to claim 6, characterized in that the historical result slope K is calculated by a least square method.
8. The Beidou side slope monitoring data post-processing method according to claim 7, wherein in the step S7, a calculation formula of post-processing is as follows:
P′=P*θ+M*(1-θ);
where P' is the final data.
CN202210188780.1A 2022-03-01 2022-03-01 Beidou slope monitoring data post-processing method Active CN114236581B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210188780.1A CN114236581B (en) 2022-03-01 2022-03-01 Beidou slope monitoring data post-processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210188780.1A CN114236581B (en) 2022-03-01 2022-03-01 Beidou slope monitoring data post-processing method

Publications (2)

Publication Number Publication Date
CN114236581A CN114236581A (en) 2022-03-25
CN114236581B true CN114236581B (en) 2022-07-26

Family

ID=80748293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210188780.1A Active CN114236581B (en) 2022-03-01 2022-03-01 Beidou slope monitoring data post-processing method

Country Status (1)

Country Link
CN (1) CN114236581B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115508873B (en) * 2022-09-30 2025-07-18 湖南五凌电力科技有限公司 Beidou continuous positioning method and system suitable for landslide geological disaster monitoring
CN118282477B (en) * 2024-03-29 2025-06-13 交通运输部水运科学研究所 A Beidou data communication method and system based on Hongmeng operating system
CN119357547B (en) * 2024-10-12 2025-06-10 长江勘测规划设计研究有限责任公司 Rapid processing method for slope monitoring data rough difference

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110554412A (en) * 2019-09-20 2019-12-10 鞍钢集团矿业有限公司 GNSS original data quality correction method for slope monitoring
CN111222095A (en) * 2019-12-19 2020-06-02 国网电力科学研究院有限公司 Gross error discrimination method, device and system in dam deformation monitoring
CN111323795A (en) * 2020-03-20 2020-06-23 湖南联智科技股份有限公司 Multi-path error weakening method in Beidou deformation monitoring
CN111735381A (en) * 2020-07-21 2020-10-02 湖南联智科技股份有限公司 Beidou monitoring result error elimination method
CN111795639A (en) * 2020-05-29 2020-10-20 湖南联智科技股份有限公司 Infrastructure structure deformation monitoring method based on Beidou high-precision positioning

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8138970B2 (en) * 2003-03-20 2012-03-20 Hemisphere Gps Llc GNSS-based tracking of fixed or slow-moving structures
US7734245B2 (en) * 2006-01-13 2010-06-08 Sai Ravela Statistical-deterministic approach to natural disaster prediction

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110554412A (en) * 2019-09-20 2019-12-10 鞍钢集团矿业有限公司 GNSS original data quality correction method for slope monitoring
CN111222095A (en) * 2019-12-19 2020-06-02 国网电力科学研究院有限公司 Gross error discrimination method, device and system in dam deformation monitoring
CN111323795A (en) * 2020-03-20 2020-06-23 湖南联智科技股份有限公司 Multi-path error weakening method in Beidou deformation monitoring
CN111795639A (en) * 2020-05-29 2020-10-20 湖南联智科技股份有限公司 Infrastructure structure deformation monitoring method based on Beidou high-precision positioning
CN111735381A (en) * 2020-07-21 2020-10-02 湖南联智科技股份有限公司 Beidou monitoring result error elimination method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A Workflow Based on SNAP–StaMPS Open-Source Tools and GNSS Data for PSI-Based Ground Deformation Using Dual-Orbit Sentinel-1 Data: Accuracy Assessment with Error Propagation Analysis;Francesco Mancini 等;《Remote Sensing》;20210218;第13卷(第4期);1-22 *
Extracting small deformation beyond individual station precision from dense Global Navigation Satellite System (GNSS) networks in France and western Europe;Christine Masson 等;《Solid Earth》;20191108;第10卷(第6期);1905-1920 *
受特定变化趋势限制的传感器数据处理方法研究;周俊华 等;《物联网技术》;20201231(第12期);90-92 *
自适应卡尔曼滤波在BDS变形监测数据处理中的应用;雷孟飞 等;《导航定位学报》;20191231;第7卷(第4期);75-79 *

Also Published As

Publication number Publication date
CN114236581A (en) 2022-03-25

Similar Documents

Publication Publication Date Title
CN114236581B (en) Beidou slope monitoring data post-processing method
WO2021184541A1 (en) Method for weakening multipath errors in beidou deformation monitoring
CN110061716B (en) An Improved Kalman Filtering Method Based on Least Squares and Multiple Fading Factors
CN107831516B (en) A method to obtain real-time high-precision displacement of dam by integrating GNSS and ground monitoring network
CN111735381B (en) Beidou monitoring result error elimination method
CN117540220B (en) Near-zero carbon park source network load matching method and system
CN113109837B (en) Data processing method of lidar system
CN110969342B (en) Method for balancing small-reservoir-capacity hydroelectric generation safety and flood control risk management and control
CN117272812A (en) A method for constructing a low-latitude small-area ionospheric model
CN116704007B (en) Sea ice area change statistical method, system, electronic equipment and medium
CN113391328A (en) Method for autonomously monitoring satellite ephemeris integrity by navigation satellite system
CN119276363A (en) A method for improving the initial pointing accuracy of laser communication terminal
CN114998213A (en) A Calculation Method of Flood Detention Capacity of Lake-reservoir Group Based on Runoff Weight Model
CN112799110A (en) Doppler-considered Beidou corrected pseudorange single-point positioning method, system and equipment
CN117554679B (en) Current monitoring method and device, electronic equipment and medium
CN118706865A (en) A GNSS-IR real-time soil moisture inversion method considering inter-system bias correction
CN117390379B (en) On-line signal measuring device and confidence measuring device for signal characteristics
CN108983301B (en) A kind of K index shows reporting method
CN112629703B (en) Method for processing distributed optical fiber temperature measurement data and improving precision
CN112799031B (en) Clutter suppression method for millimeter wave ground-like radar data
CN116680882A (en) Ocean depth benchmark construction method and device based on tide splitting recombination and tide checking constraint
CN112415559B (en) High-order fault-tolerant satellite orbit determination method based on polynomial expansion technology
CN110716215B (en) Method for improving Doppler precision
CN114004248B (en) A method for denoising against outliers based on Kalman filtering and least squares fitting
CN116227190B (en) Coastal tide level extremum prediction method based on extremum mixed distribution

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