A joint correction method and system for GNSS non-tectonic deformation and a storage medium

Abstract

The application discloses a joint correction method and system of GNSS non-tectonic deformation and a storage medium. The method obtains GNSS pseudo-range and carrier phase observation data and precise orbit, clock error and other precise products, constructs three-component displacement time series caused by non-tidal atmospheric load, non-tidal ocean load and land hydrology load and temperature-driven thermal elastic displacement time series, and unifies to a reference frame and a time system consistent with GNSS solution. The non-tectonic displacement vector is interpolated to the observation epoch by time, converted from the local coordinate system ENU to the geocentric rectangular coordinate system, and projected to the line-of-sight direction of the satellite to the reference point of the station to form an equivalent geometric distance correction term; the correction term is applied to the pseudo-range observation and the carrier phase equivalent meter domain observation in the same domain as the observation, and then GNSS daily arc segment precise solution is performed to output the coordinate time sequence, which is preferably used for precise point positioning integer ambiguity fixing technology.

Description

Technical Field

[0001] This invention belongs to the field of satellite geodesy and high-precision GNSS data processing, and involves non-tectonic deformation modeling and correction in GNSS coordinate time series. In particular, it is a processing method and system that jointly introduces non-tidal load displacement and temperature-driven thermoelastic displacement into the solution of GNSS observation equations in the form of observation layers. Background Technology

[0002] Continuous GNSS reference stations, capable of providing millimeter-level three-dimensional deformation monitoring, have been widely used for reference frame maintenance, tectonic deformation monitoring, and geophysical process research. However, the redistribution of surface mass, including atmospheric, oceanic, and terrestrial hydrological data, can cause elastic load displacements; simultaneously, near-surface temperature changes drive thermal expansion and contraction responses in the station's observation pier and its surrounding medium. These non-tectonic deformation signals can reach considerable levels on hourly to interannual scales, and are particularly prone to manifesting as inter-station common-mode errors and spatially correlated noise in dense regional station networks, thus affecting vertical coordinate repeatability, velocity estimation, and its uncertainties.

[0003] In existing engineering applications, non-tidal loads and thermoelastic displacements are more often handled using a "coordinate post-processing subtraction" method: first, coordinate estimation is performed according to a predetermined observation model, and then the displacement of the external model is subtracted at the daily coordinate time series level. This method is simple to implement, but when non-tectonic signals are not explicitly modeled during parameter estimation, they may couple with parameters such as station coordinates, tropospheric delay, clock error, and ambiguity, and be partially absorbed, resulting in residual low-frequency system errors and common-mode components after correction, making it difficult to fundamentally suppress related noise and signal leakage. In addition, although hourly resolution station displacement sequences can be obtained in recent years through layered heat conduction-elastic coupling models for thermoelastic displacement, the consistent introduction and joint correction path in the solution chain still lack reusable implementations. Summary of the Invention

[0004] To overcome the shortcomings of the prior art, this invention provides a joint correction method, system, and storage medium for GNSS non-structural deformation. The purpose of this invention is to provide a joint correction method and system for GNSS non-structural displacement observation layers. By applying geometric distance corrections to epoch-by-epoch pseudorange observations and carrier phase observations at the observation equation level, non-structural deformation is explicitly interpreted before parameter estimation, thereby reducing coupling with the parameters to be estimated and signal "partial absorption," and improving the vertical repeatability and spatial consistency of dense station network calculation results.

[0005] According to one aspect of the present invention, a joint correction method for GNSS non-structural deformation is provided, comprising:

[0006] Acquire GNSS observation data and precision products for the site to be processed; the GNSS observation data includes pseudorange observations and carrier phase observations of at least one satellite system; the precision products include orbit, clock bias, and deviation products for ambiguity fixing;

[0007] Constructing a priori corrections for nontectonic displacements includes: acquiring or calculating three-component displacement time series caused by nontectonic atmospheric load, nontectonic ocean load, and land hydrological load, as well as a thermoelastic displacement time series generated by temperature forcing; summing the three-component displacement time series to obtain a nontectonic load displacement time series; unifying the nontectonic load displacement time series and the thermoelastic displacement time series into a reference frame and time system consistent with GNSS calculation, and constructing a nontectonic displacement vector in the local coordinate system of the station.

[0008] The unconstructed displacement vector is interpolated over time to obtain the displacement increment for each observation epoch, and the displacement increment in the local coordinate system is converted into the displacement increment in the geocentric rectangular coordinate system.

[0009] For each observation epoch and for each satellite, calculate the unit vector of the line of sight from the satellite to the station reference point, and project the displacement increment in the geocentric rectangular coordinate system onto the line of sight direction to form an equivalent geometric distance correction term;

[0010] The equivalent geometric distance correction term is applied to the pseudorange observation and carrier phase observation in the same domain as the observation to obtain the corrected observation values.

[0011] Based on the corrected observations, perform precise GNSS day arc calculations and output coordinate time series.

[0012] As a further technical solution, the unconstructed displacement vector is represented as:

[0013] ,

[0014] in, It is the sum of the three components of displacement caused by non-tidal atmospheric load, non-tidal ocean load, and terrestrial hydrological load. It is a thermoelastic displacement, which includes at least a vertical component.

[0015] As a further technical solution, when the thermoelastic displacement only provides a vertical component, the horizontal component of the non-constructed displacement vector is zero, and the vertical component is the sum of the non-tidal load vertical displacement and the thermoelastic vertical displacement.

[0016] As a further technical solution, time interpolation is performed on the unconstructed displacement vector, including: calculating the displacement increment for each observation epoch by time interpolation of the unconstructed displacement vector with a sampling rate of whole-hour.

[0017] As a further technical solution, the local coordinate system is an East-North-Sky coordinate system, and the conversion of the displacement increment in the local coordinate system to the displacement increment in the geocentric rectangular coordinate system is achieved by the following formula:

[0018] ,

[0019] in, These are the unit vectors of the station in the local East-North-Sky coordinate system. They are respectively in the observed epochs The displacement of the station in the local east, north, and vertical directions.

[0020] As a further technical solution, the equivalent geometric distance correction term is calculated using the following formula:

[0021] ,

[0022] in, This is the unit vector of the line of sight from the satellite to the station reference point. This represents the displacement increment in the geocentric rectangular coordinate system.

[0023] As a further technical solution, the equivalent geometric distance correction term is applied to the pseudorange observation and carrier phase observation in the same domain as the observation, including: directly correcting the pseudorange observation in the meter domain; and multiplying the carrier phase observation by the wavelength to convert it into an equivalent meter domain quantity before performing the same domain correction.

[0024] As a further technical solution, the execution of precise GNSS day arc segment calculation includes using precise single-point positioning ambiguity fixing technology to further complete the integer ambiguity fixing and output the day coordinate time series.

[0025] According to one aspect of the present invention, a joint correction system for GNSS non-structural deformation is provided, comprising:

[0026] Data acquisition module: used to acquire GNSS observation data and precision products of the site to be processed; the GNSS observation data includes pseudorange observations and carrier phase observations of at least one satellite system; the precision products include orbit, clock bias, and deviation products for ambiguity fixing;

[0027] Non-tectonic displacement construction module: used to acquire or calculate the three-component displacement time series caused by non-tidal atmospheric load, non-tidal ocean load and land hydrological load, as well as the thermoelastic displacement time series generated by temperature forcing; sum the three-component displacement time series to obtain the non-tidal load displacement time series; unify the non-tidal load displacement time series and the thermoelastic displacement time series into a reference frame and time system consistent with GNSS calculation, and construct a non-tectonic displacement vector in the local coordinate system of the station;

[0028] Time interpolation and coordinate transformation module: used to perform time interpolation on the unconstructed displacement vector to obtain the displacement increment for each observation epoch, and convert the displacement increment in the local coordinate system into the displacement increment in the geocentric rectangular coordinate system;

[0029] Line-of-sight projection module: used to calculate the line-of-sight unit vector from the satellite to the station reference point for each observation epoch and each satellite, and project the displacement increment in the geocentric rectangular coordinate system onto the line-of-sight direction to form an equivalent geometric distance correction term;

[0030] Observation correction module: used to apply the equivalent geometric distance correction term to the pseudorange observation and carrier phase observation in the same domain as the observation, to obtain the corrected observation values;

[0031] Precision GNSS data processing module: used to perform precise GNSS day arc segment calculation based on the corrected observation values ​​and output coordinate time series.

[0032] According to one aspect of the present invention, a non-transitory computer-readable storage medium is provided, the non-transitory computer-readable storage medium storing computer instructions that cause the computer to perform the steps of the method.

[0033] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0034] (1) Apply geometric correction to the epoch-by-epoch observations at the level of the observation equations so that non-tidal loads and thermoelastic displacements are explicitly explained before parameter estimation, thereby reducing the coupling with parameters such as station coordinates, troposphere and ambiguity and reducing signal leakage.

[0035] (2) By jointly introducing multiple types of environmental displacements and ensuring the self-consistency between the reference frame and the time system, the common mode error of the station network and spatial correlation noise can be suppressed, and the consistency and repeatability of the vertical time series can be improved.

[0036] (3) Input displacement supports hourly resolution and multi-source products. Interpolation and projection are simple to implement, reusable, and easy to batch process, making it suitable for regional dense station networks and global scale operational deployment. Attached Figure Description

[0037] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the accompanying drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0038] Figure 1 A flowchart of a joint correction method for non-structural deformations in GNSS provided in an embodiment of the present invention;

[0039] Figure 2 This is a schematic diagram of the system module composition provided in an embodiment of the present invention;

[0040] Figure 3 A comparative schematic diagram of the GNSS station vertical daily coordinate sequences under three conditions: uncorrected, parameter-layer corrected, and observation-layer combined corrected, provided for embodiments of the present invention, with non-tidal load (NTAL, NTOL, HYDL) and thermoelastic displacement time series;

[0041] Figure 4 This is a statistical comparison diagram of the RMS and WRMS of the vertical residuals of a GNSS station under three processing methods: uncorrected, parameter-level corrected, and observation-level combined corrected according to the present invention. It includes a box plot comparison of the WRMS of the residuals of 15 stations to illustrate the differences in the vertical residual suppression effect of different correction methods. Detailed Implementation

[0042] The terms “comprising” and “having”, and any variations thereof, in the specification, claims, and accompanying drawings of this invention are intended to cover a non-exclusive inclusion, such as a process, method, system, product, or apparatus that includes a series of steps or units, not necessarily limited to those explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0043] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention. In addition, the technical features of the various embodiments or individual embodiments provided by the present invention can be arbitrarily combined to form new technical solutions. Such combinations are not bound by the order of steps and / or structural composition patterns, but must be based on the ability of those skilled in the art to implement them. When the combination of technical solutions is contradictory or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection claimed by the present invention.

[0044] This invention provides a joint correction method for non-structural deformations in GNSS, such as... Figure 1 As shown, it includes the following steps:

[0045] Step 1: Data Acquisition and Standard Model Setup.

[0046] Acquire GNSS observation data and precision products for the site to be processed. The GNSS observation data includes pseudorange observations and carrier phase observations from at least one satellite system; the precision products include orbit, clock bias, and deviation products used for ambiguity fixing, and set up conventional tidal correction models such as solid tides, polar tides, and ocean tidal loads.

[0047] Step 2: Construct the a priori correction for non-constructive displacement.

[0048] Acquire or calculate the three-component displacement time series caused by non-tidal atmospheric load, non-tidal ocean load, and terrestrial hydrological load, as well as the thermoelastic displacement time series generated by temperature forcing; unify the two types of displacements into a reference frame and time system consistent with GNSS calculations, and construct a non-structural displacement vector in the station's local coordinate system (ENU):

[0049] ,

[0050] in, The sum of three types of non-tidal load displacements (providing three ENU components); It is a thermoelastic displacement, which can be a vertical component or a three-component component. These represent the displacements of the station in the local east, north, and vertical directions, respectively, caused by the combined effects of non-tidal loads and thermoelasticity at observation epoch t.

[0051] If thermoelasticity only provides the vertical component, then

[0052] ,

[0053] in, This represents the vertical displacement of the non-tidal load. This represents the thermoelastic displacement in the vertical direction. This represents the thermoelastic displacement vector at time t. This represents the vertical component of the non-constructive total displacement vector at time t.

[0054] Step 3: Time interpolation and coordinate transformation.

[0055] The observation layer scheme introduces geometric corrections for displacement into the observation equations in a "forward" manner. Because... Typically, sampling is done at the hour (e.g., 1 hour), while GNSS observations use high-frequency epochs (e.g., 30 seconds). To introduce corrections into the observation equations, it is necessary to... Time interpolation is performed to obtain each observation epoch. of This embodiment preferably uses linear interpolation, but other interpolation methods such as spline interpolation can also be used.

[0056] To facilitate subsequent projection onto the satellite line of sight, the displacement increments in the site's local coordinate system (ENU) need to be converted to displacement increments in the geocentric rectangular coordinate system (ECEF):

[0057] ,

[0058] in, is the unit vector in the local ENU coordinate system at the station.

[0059] Step 4: Project the line of sight and construct the geometric distance correction term.

[0060] For each observation epoch For each satellite, calculate the unit vector of the line of sight from the satellite to the station reference point. The displacement of the station reference point is projected onto the line of sight to form an equivalent geometric distance correction term. Based on first-order geometric relations, when the station undergoes displacement... When the geometric distance correction is applied, it can be expressed as:

[0061] ,

[0062] The physical meaning of this formula is: when the station undergoes displacement... At that time, the geometric distance between the satellite and the station will shorten or increase accordingly, and the change is... .

[0063] Step 5: Apply corrections in the same domain and perform daily arc segment calculation.

[0064] The geometric distance correction term is applied to the original observations in the same domain as the observations. Specifically, pseudorange observations are directly corrected in the meter domain; carrier phase observations are corrected using an equivalent meter domain quantity. Perform local correction, among which For carrier wavelength, For phase observations.

[0065] After obtaining the corrected observations, precise GNSS daily arc segment calculations are performed under the same observation model, parameter estimation strategy, and stochastic model settings as the uncorrected scheme. The observation model includes GNSS pseudorange observation equations and carrier phase observation equations, and the stochastic model includes observation weight matrices set for pseudorange and carrier phase observations. Parameters estimated during the calculation process include: station coordinates, tropospheric delay, clock bias (or receiver clock bias parameters), and ambiguity. When using Precise Point Positioning Ambiguity Fixing (PPP-AR) technology, the integer characteristics of ambiguity recovery from the deviation product obtained in step 1 are utilized to complete integer ambiguity fixing, outputting a high-precision daily coordinate time series.

[0066] By applying the above strategy of directly correcting at the observation equation level, non-structural deformation is explicitly introduced into the geometric model before parameter estimation, thereby minimizing its coupling and "partial absorption" with parameters such as station coordinates, tropospheric delay, and ambiguity, and reducing the possibility that it is retained in the residuals in the form of common mode.

[0067] Step 6: Output and Evaluation.

[0068] The corrected station coordinate time series and residual statistics are output. To further quantify the improvement effect of the method in this embodiment, the weighted root mean square, spectral energy, band-limited root mean square, common mode error and spatial correlation can be calculated to evaluate the degree of improvement of the accuracy and noise structure of the dense station network solution by the joint correction of the observation layer.

[0069] Figure 3 The paper presents a comparison of the GNSS station vertical daily coordinate series under three conditions: non-tidal load (NTAL, NTOL, HYDL) and thermoelastic displacement time series, as well as the GNSS station vertical daily coordinate series under three conditions: uncorrected, parameter-layer corrected, and observation-layer combined corrected according to the embodiments of the present invention. Figure 4 This paper presents a statistical comparison diagram of the RMS and WRMS of the vertical residuals of a GNSS station under three processing methods: uncorrected, parameter-level corrected, and observation-level combined corrected according to the embodiments of the present invention. The diagram includes a box plot comparison of the WRMS of the residuals of 15 stations to illustrate the differences in the vertical residual suppression effect of different correction methods.

[0070] This invention provides a joint correction system for GNSS non-structural deformations, used to implement the methods described in the foregoing embodiments. For example... Figure 2 As shown, the system includes:

[0071] Data acquisition module: Used to acquire GNSS observation data and precision products from the site to be processed. The GNSS observation data includes pseudorange observations and carrier phase observations from at least one satellite system; the precision products include orbit, clock bias, and deviation products used for ambiguity fixing.

[0072] The non-tidal displacement construction module is used to read in or calculate non-tidal load displacement and thermoelastic displacement and complete reference frame unification. Specifically, it includes: acquiring or calculating the three-component displacement time series caused by non-tidal atmospheric load, non-tidal ocean load, and land hydrological load, as well as the thermoelastic displacement time series generated by temperature forcing; summing the three-component displacement time series to obtain the non-tidal load displacement time series; unifying the non-tidal load displacement time series and the thermoelastic displacement time series into a reference frame and time system consistent with GNSS calculation, and constructing a non-tidal displacement vector in the station's local coordinate system.

[0073] The time interpolation module is used to convert hourly sampled displacements into observed epoch displacements. Specifically, it performs time interpolation on the unconstructed displacement vector to obtain the displacement increment for each observed epoch.

[0074] Coordinate transformation module: used to convert the displacement increment in the local coordinate system into the displacement increment in the geocentric rectangular coordinate system.

[0075] Line-of-sight projection module: For each observation epoch and for each satellite, it calculates the unit vector of the line of sight from the satellite to the station reference point, and projects the displacement increment in the geocentric rectangular coordinate system onto the line of sight direction to form an equivalent geometric distance correction term.

[0076] Observation correction module: used to apply the equivalent geometric distance correction term to the pseudorange observation and carrier phase observation in the same domain as the observation, to obtain the corrected observation values.

[0077] Precision GNSS data processing module: used to perform precise GNSS day arc calculations based on the corrected observations.

[0078] Output and Evaluation Module: Used to output the corrected station coordinate time series and residual statistics; and to further calculate indicators such as weighted root mean square, spectral energy, band-limited root mean square, common mode error and spatial correlation, to evaluate the degree of improvement of the accuracy and noise structure of dense station network solution by the joint correction of the observation layer.

[0079] The specific working principles of each module correspond one-to-one with the method steps in the aforementioned embodiments, and will not be repeated here.

[0080] This invention also provides a computer-readable storage medium storing a computer program thereon. When executed by a processor, the computer program implements the steps of the joint correction method for GNSS non-structural deformation as described in the foregoing embodiments.

[0081] Specifically, when the computer program is executed by the processor, it performs the following steps:

[0082] Acquire GNSS observation data and precision products for the site to be processed; the GNSS observation data includes pseudorange observations and carrier phase observations of at least one satellite system; the precision products include orbit, clock bias, and deviation products for ambiguity fixing;

[0083] Constructing a priori corrections for nontectonic displacements includes: acquiring or calculating three-component displacement time series caused by nontectonic atmospheric load, nontectonic ocean load, and land hydrological load, as well as a thermoelastic displacement time series generated by temperature forcing; summing the three-component displacement time series to obtain a nontectonic load displacement time series; unifying the nontectonic load displacement time series and the thermoelastic displacement time series into a reference frame and time system consistent with GNSS calculation, and constructing a nontectonic displacement vector in the local coordinate system of the station.

[0084] The unconstructed displacement vector is interpolated over time to obtain the displacement increment for each observation epoch, and the displacement increment in the local coordinate system is converted into the displacement increment in the geocentric rectangular coordinate system.

[0085] For each observation epoch and for each satellite, calculate the unit vector of the line of sight from the satellite to the station reference point, and project the displacement increment in the geocentric rectangular coordinate system onto the line of sight direction to form an equivalent geometric distance correction term;

[0086] The equivalent geometric distance correction term is applied to the pseudorange observation and carrier phase observation in the same domain as the observation to obtain the corrected observation values.

[0087] Based on the corrected observations, perform precise GNSS day arc calculations and output coordinate time series.

[0088] Preferably, the precise calculation of the GNSS day arc segment includes using Precise Point Positioning Ambiguity Fixing (PPP-AR) technology, utilizing the deviation product used for ambiguity fixing to complete the integer ambiguity fixing, and outputting the day coordinate time series.

[0089] The computer-readable storage media include, but are not limited to, any non-transitory medium capable of storing program code, such as read-only memory (ROM), random access memory (RAM), magnetic disk, optical disk, flash memory, solid-state drive (SSD).

[0090] In summary, this invention discloses a joint correction method for non-tectonic displacement observation layers for Global Navigation Satellite Systems (GNSS). This method acquires GNSS pseudorange and carrier phase observation data, as well as precise products such as precise orbits and clock errors. It constructs a three-component displacement time series caused by non-tidal atmospheric loads, non-tidal ocean loads, and land hydrological loads, along with a temperature-driven thermoelastic displacement time series, and unifies them to a reference frame and time system consistent with GNSS calculations. The non-tectonic displacement vectors are interpolated to the observation epoch, transformed from the station's local coordinate system (ENU) to a geocentric rectangular coordinate system, and projected onto the line-of-sight direction from the satellite to the station's reference point to form an equivalent geometric distance correction term. This correction term is applied to the pseudorange observations and the carrier phase equivalent meter-domain observations in the same domain as the observations, respectively, and then a precise GNSS daily arc segment calculation is performed to output the coordinate time series. Preferably, this method can be used for precise single-point positioning integer ambiguity fixing technology. This invention can explicitly explain environmental displacement before parameter estimation, reducing its coupling with station coordinates, tropospheric delay, and ambiguity waiting parameters, as well as signal leakage. It also suppresses common-mode errors and spatially correlated noise in dense station networks, improves vertical solution stability and repeatability, and is suitable for regional dense station networks and operational high-precision GNSS processing.

[0091] It should be noted that the joint correction of the observation layer in this invention is not limited to a specific load product or thermoelastic model; as long as a station displacement sequence aligned with the observation time can be provided and can be represented within a unified reference frame, the in-domain correction and PPP-AR solution can be completed according to the steps of this invention. Furthermore, this invention can also be extended to data processing in network solutions or other spatial geodesy techniques to assess or suppress associated noise caused by environmental displacement.

[0092] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the technical solutions of the embodiments of the present invention.

Claims

1. A joint correction method for non-structural deformation in GNSS, characterized in that, include: Acquire GNSS observation data and precision products for the site to be processed; The GNSS observation data includes pseudorange observations and carrier phase observations from at least one satellite system. The precision products include tracks, clock errors, and deviation products used for ambiguity fixing; Constructing a priori corrections for nontectonic displacements includes: acquiring or calculating three-component displacement time series caused by nontectonic atmospheric load, nontectonic ocean load, and land hydrological load, as well as a thermoelastic displacement time series generated by temperature forcing; summing the three-component displacement time series to obtain a nontectonic load displacement time series; unifying the nontectonic load displacement time series and the thermoelastic displacement time series into a reference frame and time system consistent with GNSS calculation, and constructing a nontectonic displacement vector in the local coordinate system of the station. The unconstructed displacement vector is interpolated over time to obtain the displacement increment for each observation epoch, and the displacement increment in the local coordinate system is converted into the displacement increment in the geocentric rectangular coordinate system. For each observation epoch and for each satellite, calculate the unit vector of the line of sight from the satellite to the station reference point, and project the displacement increment in the geocentric rectangular coordinate system onto the line of sight direction to form an equivalent geometric distance correction term; The equivalent geometric distance correction term is applied to the pseudorange observation and carrier phase observation in the same domain as the observation to obtain the corrected observation values. Based on the corrected observations, perform precise GNSS day arc calculations and output coordinate time series.

2. The joint correction method for GNSS non-structural deformation according to claim 1, characterized in that, The unconstructed displacement vector is represented as: ， in, It is the sum of the three components of displacement caused by non-tidal atmospheric load, non-tidal ocean load, and terrestrial hydrological load. It is a thermoelastic displacement, which includes at least a vertical component.

3. The joint correction method for GNSS non-structural deformation according to claim 2, characterized in that, When the thermoelastic displacement provides only the vertical component, the horizontal component of the unconstructed displacement vector is zero, and the vertical component is the sum of the non-tidal load vertical displacement and the thermoelastic vertical displacement.

4. The joint correction method for GNSS non-structural deformation according to claim 1, characterized in that, The unconstructed displacement vector is subjected to time interpolation, which includes: calculating the displacement increment for each observation epoch from the unconstructed displacement vector with a sampling rate of whole-hour time through time interpolation.

5. The joint correction method for GNSS non-structural deformation according to claim 1, characterized in that, The local coordinate system is the East-North-Sky coordinate system. The conversion of displacement increments in the local coordinate system to displacement increments in the geocentric rectangular coordinate system is achieved by the following formula: ， in, These are the unit vectors of the station in the local East-North-Sky coordinate system. They are respectively in the observed epochs The displacement of the station in the local east, north, and vertical directions.

6. The joint correction method for GNSS non-structural deformation according to claim 1, characterized in that, The equivalent geometric distance correction term is calculated using the following formula: ， in, This is the unit vector of the line of sight from the satellite to the station reference point. This represents the displacement increment in the geocentric rectangular coordinate system.

7. The joint correction method for GNSS non-structural deformation according to claim 1, characterized in that, Applying the equivalent geometric distance correction term to pseudorange and carrier phase observations in the same domain as the observation includes: directly correcting pseudorange observations in the meter domain; and multiplying carrier phase observations by wavelength to convert them into equivalent meter domain quantities before applying the same domain correction.

8. The joint correction method for GNSS non-structural deformation according to claim 1, characterized in that, The execution of precise GNSS day arc segment calculation includes using precise single-point positioning ambiguity fixing technology to further complete integer ambiguity fixing and outputting day coordinate time series.

9. A joint correction system for non-structural deformation in GNSS, characterized in that, include: Data acquisition module: used to acquire GNSS observation data and precision products from the site to be processed; The GNSS observation data includes pseudorange observations and carrier phase observations from at least one satellite system; the precision products include orbit, clock bias, and deviation products for ambiguity fixing. Non-tectonic displacement construction module: used to acquire or calculate the three-component displacement time series caused by non-tidal atmospheric load, non-tidal ocean load and land hydrological load, as well as the thermoelastic displacement time series generated by temperature forcing; sum the three-component displacement time series to obtain the non-tidal load displacement time series; unify the non-tidal load displacement time series and the thermoelastic displacement time series into a reference frame and time system consistent with GNSS calculation, and construct a non-tectonic displacement vector in the local coordinate system of the station; Time interpolation and coordinate transformation module: used to perform time interpolation on the unconstructed displacement vector to obtain the displacement increment for each observation epoch, and convert the displacement increment in the local coordinate system into the displacement increment in the geocentric rectangular coordinate system; Line-of-sight projection module: used to calculate the line-of-sight unit vector from the satellite to the station reference point for each observation epoch and each satellite, and project the displacement increment in the geocentric rectangular coordinate system onto the line-of-sight direction to form an equivalent geometric distance correction term; Observation correction module: used to apply the equivalent geometric distance correction term to the pseudorange observation and carrier phase observation in the same domain as the observation, to obtain the corrected observation values; Precision GNSS data processing module: used to perform precise GNSS day arc segment calculation based on the corrected observation values ​​and output coordinate time series.

10. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer instructions that cause the computer to perform the steps of the method according to any one of claims 1 to 8.

Images