CN115616630A - Method for reducing tide distribution error in research of sea tide load displacement by utilizing GNSS - Google Patents

Method for reducing tide distribution error in research of sea tide load displacement by utilizing GNSS Download PDF

Info

Publication number
CN115616630A
CN115616630A CN202211260841.7A CN202211260841A CN115616630A CN 115616630 A CN115616630 A CN 115616630A CN 202211260841 A CN202211260841 A CN 202211260841A CN 115616630 A CN115616630 A CN 115616630A
Authority
CN
China
Prior art keywords
tide
partial
sea
point
amplitude
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.)
Pending
Application number
CN202211260841.7A
Other languages
Chinese (zh)
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.)
First Institute of Oceanography MNR
Original Assignee
First Institute of Oceanography MNR
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 First Institute of Oceanography MNR filed Critical First Institute of Oceanography MNR
Priority to CN202211260841.7A priority Critical patent/CN115616630A/en
Publication of CN115616630A publication Critical patent/CN115616630A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/40—Correcting position, velocity or attitude
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention belongs to the field of ocean monitoring, and provides a method for reducing tide distribution errors in the research of sea tide load displacement by utilizing GNSS, which comprises the following steps: determining N 2 To M is aligned with 2 、S 2 To M 2 、K 2 To M 2 、K 2 To S 2 、P 1 To K 1 、O 1 To K 1 、Q 1 To K 1 And Q 1 To O 1 The tide pair of (1); establishing a relative relation of each tide separating pair, and acquiring an amplitude ratio and a lag angle difference of each grid point; calculating a balance tide harmonic constant to obtain the balance tide amplitude of each partial tide after ground tide correction; calculating a GPS data harmonic constant, and acquiring the amplitude of each tide obtained by the GPS data harmonic analysis; and correcting the estimation errors of the K1 and the K2 tide according to the results of the steps and by using the relative admittance. The invention establishes 8 relative relations of tide pairs, and can effectively obtain the relative admittance relation of each tide family. KnotThe ratio of the harmonic constant and the corresponding balance tide harmonic constant is analyzed by the jointly adopted GPS data, so that K can be effectively corrected 1 And K 2 Error in tide division and increase of K 1 And K 2 Accuracy of tide separation.

Description

Method for reducing tide distribution error in research of sea tide load displacement by utilizing GNSS
Technical Field
The invention belongs to the technical field of ocean monitoring, and particularly relates to a method for reducing tide distribution errors in the research of sea tide load displacement by utilizing GNSS.
Background
With the rapid development of the GNSS technology and the data processing method, the displacement deformation of the observation station caused by the sea tide load can be accurately obtained through the GNSS, and then the sea tide load displacement harmonic constant is obtained, so that a new method is provided for researching the sea tide load effect. Generally, when a GNSS is used to research sea tide load displacement, GPS observation data is mainly used, and a period of a main tide should be estimated in advance when the GPS observation data is subjected to harmonic analysis, whereas a period of a GPS satellite orbit (about 11.967 hours) is very similar to a period of a K2 tide (about 11.9672 hours), and a repetition period of a GPS constellation (about 23.9319 hours) is very similar to a period of a K1 tide (about 23.9345 hours), which is easy to generate a coupling phenomenon, and when a sea tide load displacement parameter is extracted by using a GNSS technology, an error with a similar period is easily absorbed by a K1/K2 tide corresponding to the same period, resulting in a large K1/K2 tide error estimated by the harmonic analysis. Therefore, how to solve the problem that the estimation error of the K1/K2 tide division is overlarge is a difficult point and a hot point problem of researching sea tide load displacement by GNSS.
Disclosure of Invention
The invention utilizes S 2 To M 2 ，N 2 To M 2 ，K 2 To S 2 ，O 1 To K 1 ，P 1 To K 1 And Q 1 To O is 1 Correcting GPS observation to obtain K by relative admittance of the tide 1 And K 2 And the tide division estimation error is too large.
The invention provides a method for reducing tide distribution errors in the research of sea tide load displacement by utilizing GNSS, which is characterized by comprising the following steps:
step 1, determining a tide separating pair; the determined tide pair is determined as N 2 To M 2 、S 2 To M 2 、K 2 To M 2 、K 2 To S 2 、P 1 To K 1 、O 1 To K 1 、Q 1 To K 1 And Q 1 To O 1 The tide dividing pair is adopted;
step 2, establishing the relative relation of each tide pair, and acquiring the amplitude ratio of each grid point
Figure BDA0003890906760000021
And difference in retardation
Figure BDA0003890906760000022
Step 3, calculating the equilibrium tide harmonic constant to obtain the equilibrium tide amplitude of each tide after the correction of the ground tide
Figure BDA0003890906760000023
Step 4, calculating a GPS data harmonic constant, and acquiring the amplitude H of each tide obtained by the GPS data harmonic analysis Gpsi (ii) a According to the tide theory, assuming that ζ is the instantaneous displacement of the sea tide load to one direction in the observation station, the harmonic constant of the main tide division can be expressed as,
Figure BDA0003890906760000024
in the formula, S 0 Is the average position of the displacement; h and g are the harmonic constants of the partial tides; f. of i And u i Correcting angle, V, for the cross point factor and cross point of partial tide 0 The phase of the ith partial tide of the equilibrium tide at the moment t = 0; omega i Is the tide division angular rate; where f, ω t and v 0 + u is time dependent, location independent; and H and g are location dependent, time independent; therefore, the above formula can be written as
Figure BDA0003890906760000025
Wherein the content of the first and second substances,
Figure BDA0003890906760000026
Figure BDA0003890906760000027
according to the principle of least square method, solving X i And Y i Then, the harmonic constant H is calculated according to the following formula Gps And g Gps ，
Figure BDA0003890906760000028
Step 5, correcting K by using relative admittance according to the results of the step 2, the step 3 and the step 4 1 And K 2 And (4) tidal score estimation error.
In a possible design, the step 1 specifically comprises the following steps:
for a certain tide i, admittance M i Is composed of
Figure BDA0003890906760000031
Wherein the content of the first and second substances,
Figure BDA0003890906760000032
H i and g i Harmonic constants, amplitudes and retardations, C, for the partial tides i For the tidal power coefficient, the relationship between the main partial tide n and the secondary partial tide m in the same tide family can be expressed as
Figure BDA0003890906760000033
Wherein
Figure BDA0003890906760000034
a m/n ＝g m -g n (4)
H in the formula (3) m /H n I.e. the amplitude ratio, g, in conventional tidal science m -g n Retardation difference, A * m/n Is the amplitude relationship of the primary partial tide n relative to the secondary partial tide m, a m/n The lag angle relation of the main partial tide n relative to the secondary partial tide m is shown; if the amplitude ratio of the actual partial tide is the same as the ratio of the induced tide force, the current tide is divided into a first tide and a second tide
Figure BDA0003890906760000035
Equal to 1; in the real sea
Figure BDA0003890906760000036
It will not be exactly equal to 1, but it is close to 1 for most of the sea; if the two lag angles of partial tide are equal, then a m/n Is equal to zero; in the actual ocean a m/n Generally not equal to zero, and is greater than zero or less than zero for most sea areas;
determining N 2 To M is aligned with 2 、S 2 To M is aligned with 2 、K 2 To M is aligned with 2 、K 2 To S 2 、P 1 To K 1 、O 1 To K 1 、Q 1 To K 1 And Q 1 To O is 1 Dividing the tide pair;
Figure BDA0003890906760000037
g m ＝a m/n +g n (6)
let the amplitude ratio be
Figure BDA0003890906760000041
The difference in retardation angle is respectively
Figure BDA0003890906760000042
Setting the amplitude ratio as H and the lag angle difference as g for two partial tides, and recording A = Hcosg and B = Hsing;
the specific process of the step 2 comprises the following steps:
it is known that there are N observation points, each represented by i (i =1,2.., N), each having a longitude and latitude of
Figure BDA0003890906760000043
Wherein
Figure BDA0003890906760000044
Negative values for south picks; the number of the lambda-beams is increased,
Figure BDA0003890906760000045
degree is taken as a unit;
Figure BDA0003890906760000046
and
Figure BDA0003890906760000047
the amplitude ratio and the lag angle difference of observation points of each tide checking station are obtained;
for observation station i, calculating the distance between the observation station i and other stations j
Figure BDA0003890906760000048
Weighted average of i points
Figure BDA0003890906760000049
Figure BDA00038909067600000410
Wherein the content of the first and second substances,
Figure BDA00038909067600000411
represents summing all observation points, but excluding the i point; r i Is the corrected Kelvin wave length, which is determined according to the water depth near the point i, and the method is as follows:
setting the water depth value of each grid point in and around the known sea area, namely the longitude and latitude of the point m
Figure BDA00038909067600000412
Figure BDA00038909067600000413
Depth of water D m (ii) a For the point i, calculating the distance between the grid point m and the tide checking station i as
Figure BDA00038909067600000414
Select all r im ＜150×10 3 m, the water depth value of which is averaged and is recorded as
Figure BDA00038909067600000415
Calculating corresponding Kelvin wave velocity
Figure BDA0003890906760000051
Kelvin wave length
L i ＝u i T (13)
Wherein u is i At a Kelvin wave velocity and T at a tidal cycle, i.e.
T＝360°/ω (14)
ω is the tidal angular rate;
finally, selecting
R i ＝κL i (κ＝1/4) (15)
Calculate for all points i =1,2
Figure BDA0003890906760000052
And
Figure BDA0003890906760000053
and calculate
Figure BDA0003890906760000054
And root mean square deviation
Figure BDA0003890906760000055
If it is
δ i ≥mσ (18)
The data of the ith point is considered to be abnormal and needs to be omitted;
to verify the accuracy of the removed outliers, the method comprises
Figure BDA0003890906760000056
And
Figure BDA0003890906760000057
inversely calculating amplitude ratio of each observation point
Figure BDA0003890906760000058
And difference of retardation angle
Figure BDA0003890906760000059
Figure BDA00038909067600000510
Figure BDA0003890906760000061
Then by
Figure BDA0003890906760000062
Figure BDA0003890906760000063
Calculating the mean square error of the amplitude ratio and the retardation
Figure BDA0003890906760000064
Figure BDA0003890906760000065
According to the method recorded in the formula (18), unreasonable observed values are abandoned, grid points are set according to a preset resolution, the serial number of each grid point is recorded as k (k =1, 2.) and the corresponding longitude and latitude are recorded as
Figure BDA0003890906760000066
Calculate the distance from point k to each observation point i as
Figure BDA0003890906760000067
Weighted average of k points
Figure BDA0003890906760000068
Figure BDA0003890906760000069
Wherein R is k Calculation method and R in the above formula (15) i The calculation method of (2) is the same,
by
Figure BDA00038909067600000610
And
Figure BDA00038909067600000611
inversely calculating the amplitude ratio of each grid point
Figure BDA00038909067600000612
And difference in retardation
Figure BDA00038909067600000613
Figure BDA00038909067600000614
Figure BDA0003890906760000071
In a possible design, the specific process of step 3 is:
for any observation point
Figure BDA0003890906760000072
Its corresponding equilibrium tidal height ζ EQ Can be expressed as
Figure BDA0003890906760000073
In the formula, f i And u i Correcting angle, V, for the cross point factor and cross point of the partial tide 0 The phase of the ith partial tide of the equilibrium tide at the moment t = 0; omega i Is the tidal angular rate; p is a family number, p =1 for full solar tide, p =2 for half solar tide; λ is longitude; s = -8h is a Beijing standard time zone number;
Figure BDA0003890906760000074
the balanced tide amplitude of each partial tide after the correction of the ground tide,
Figure BDA0003890906760000075
Figure BDA0003890906760000076
Figure BDA0003890906760000077
Figure BDA0003890906760000078
in the formula (I), the compound is shown in the specification,
Figure BDA0003890906760000079
is the latitude.
In a possible design, the specific process of step 5 is:
order to
Figure BDA00038909067600000710
In the formula, H Gpsi For the amplitudes of the partial tides derived from the GPS data harmonic analysis,
Figure BDA00038909067600000711
for each partial tide to balance tide amplitude psi i For modulus, divide each half-day tide by N 2 、M 2 、S 2 And K 2 Is curve fitted to the psi value of (A) and then K is fitted 2 Correcting abnormal psi value of the tide to a normal fitting curve, thereby completing K 2 Correction of moisture separation error, and then increase K 2 Accuracy of the damping and reconciliation constants; in the same way, divide each full day into tides K 1 、P 1 、O 1 And Q 1 Is curve fitted to the psi value of (A) and then K is fitted 1 Correcting abnormal psi value of the tide to a normal fitting curve, thereby completing K 1 Correction of tide-dividing error, and then raising K 1 Accuracy of the damping and harmonic constants.
The second aspect of the present invention also provides an apparatus for reducing tidal separation error when researching sea tide load displacement by using GNSS, the apparatus comprising at least one processor and at least one memory, the processor and the memory being coupled; the memory having stored therein a computer program or instructions; the processor, when executing the computer program or the instructions, may implement the method for reducing tide distribution error when researching sea tide load displacement by using GNSS as described in the first aspect.
The third aspect of the present invention further provides a computer-readable storage medium, in which a computer program or instructions are stored, and when the computer program or instructions are executed by a processor, the method for reducing tide distribution error when researching sea tide load displacement by using GNSS as described in the first aspect may be implemented.
Compared with the prior art, the invention provides a method for reducing tide distribution errors in the research of sea tide load displacement by utilizing GNSS, and has the following beneficial effects:
1. for GPS orbit period and K 2 The tide-separating period is consistent, and the repetition period of the GPS constellation is equal to K 1 The periods of partial tide are consistent to generate coupling effect, resulting in K 2 Divide tide and K 1 And (4) the difficulty of overlarge tide-dividing estimation error. The invention provides a K obtained by correcting GPS observation based on relative admittance 1 And K 2 The method for dividing the tide error can effectively improve the K observed by the GPS 1 And K 2 Error in tide division and increase of K 1 And K 2 Accuracy of tide separation.
2. The displacement deformation of the observation station caused by sea tide load is accurately obtained through GNSS, and further the sea tide load displacement harmonic constant is obtained.
3. In order to improve the accuracy of researching the sea tide load effect, the invention adopts a method for establishing the relative relationship of 8 tide divisions, and can effectively obtain the relative admittance relationship of each tide family. The ratio of the harmonic constant and the corresponding balance tide harmonic constant is analyzed by combining the GPS data adopted by the invention, so that the K can be effectively corrected 1 And K 2 Error in tide division and increase of K 1 And K 2 Accuracy of tide separation.
Drawings
FIG. 1 illustrates the correction of the present inventionK observed by GPS 1 And K 2 And a flow chart of the tide separation error method.
FIG. 2 shows the correction K of the present invention 2 Schematic representation of tidal separation error.
FIG. 3 is a diagram of K obtained by correcting GPS observations in accordance with the present invention 1 And K 2 The simple and easy sketch of equipment structure of dividing the tide error.
Detailed Description
The invention is further illustrated by the following examples.
Example 1:
correction of GPS observations based on relative admittance 1 And K 2 Calculation of tide-dividing error includes calculation of relative admittance and correction of GPS observed K 1 And K 2 The general flow of the calculation steps of the tide separation error is shown in figure 1.
Determining tide separating pair
For a certain tide i, admittance M i Is composed of
Figure BDA0003890906760000091
Wherein
Figure BDA0003890906760000092
H i And g i The harmonic constant, amplitude and lag angle of the tide are taken as the parameters; c i Is the tidal force index. In the same tide family, the relationship between the primary tide n and the secondary tide m can be expressed as
Figure BDA0003890906760000093
Wherein
Figure BDA0003890906760000094
a m/n ＝g m -g n (4) H in the formula (3) m /H n Is to transmitAmplitude ratio in statics, g m -g n Is the retardation angle difference. If the amplitude ratio of the actual partial tide is the same as the ratio of the induced tide force, the current tide is divided into a first tide and a second tide
Figure BDA0003890906760000101
Equal to 1. In the real sea
Figure BDA0003890906760000102
It will not be exactly equal to 1, but it is close to 1 for most of the sea. If the two tide-dividing lags are equal, then a m/n Equal to zero. In the real ocean a m/n It is generally not equal to zero and is greater than zero or less than zero for most of the sea.
Determining N 2 To M is aligned with 2 ，S 2 To M is aligned with 2 ，K 2 To M is aligned with 2 ，K 2 To S 2 ，P 1 To K 1 ，O 1 To K 1 ，Q 1 To K 1 And Q 1 To O is 1 And (5) dividing the tide pair.
Figure BDA0003890906760000103
g m ＝a m/n +g n (6)
Let the amplitude ratio be
Figure BDA0003890906760000104
A difference in retardation angle of
Figure BDA0003890906760000105
For two partial tides, let the amplitude ratio be H, the lag angle difference be g, note A = Hcosg, B = Hsing.
(II) establishing the relative relationship between the partial tide pairs
It is known that there are N observation points, each represented by i (i =1,2.., N), each having a longitude and latitude of
Figure BDA0003890906760000106
Wherein
Figure BDA0003890906760000107
Negative values are indicated for southern picks. The number of the lambda-beams is increased,
Figure BDA0003890906760000108
in degrees. H i And g i And the amplitude ratio and the lag angle difference of observation points of the tide gauging stations are obtained.
For observation station i, calculate its distance to each other station j
Figure BDA0003890906760000109
Weighted average of i points
Figure BDA00038909067600001010
Figure BDA0003890906760000111
Wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003890906760000112
the representation sums all observation points, but does not include the i point. R i The method is determined according to the water depth near the point i and comprises the following steps:
setting the water depth value of each grid point in and around the known sea area, namely the longitude and latitude of the point m
Figure BDA0003890906760000113
Figure BDA0003890906760000114
Water depth D m . For the point i, calculating the distance between the grid point m and the tide gauge point i as
Figure BDA0003890906760000115
Select all r im ＜150×10 3 m points, averaging the water depth values at these points, and recording as
Figure BDA0003890906760000116
Calculating corresponding Kelvin wave velocity
Figure BDA0003890906760000117
Kelvin wave length
L i ＝u i T (13)
Wherein T is the tidal cycle, i.e.
T＝360°/ω (14)
Where ω is the tidal angular rate.
Finally, selecting
R i ＝κL i (κ＝1/4) (15)
For all i =1,2
Figure BDA0003890906760000118
And
Figure BDA0003890906760000119
and calculate
Figure BDA00038909067600001110
And root mean square deviation
Figure BDA00038909067600001111
If it is
δ i ≥mσ (18)
The ith data is considered likely to be anomalous, considering whether to drop it, where m may take 2.
For verificationAccuracy after removal of the outlier, is determined by
Figure BDA0003890906760000121
And
Figure BDA0003890906760000122
inversely calculating amplitude ratio of each observation point
Figure BDA0003890906760000123
And difference of retardation angle
Figure BDA0003890906760000124
Figure BDA0003890906760000125
Figure BDA0003890906760000126
Then is provided with
Figure BDA0003890906760000127
Figure BDA0003890906760000128
Calculating the mean square error of the amplitude ratio and the delay angle difference
Figure BDA0003890906760000129
Figure BDA00038909067600001210
According to the method, unreasonable observation values are abandoned, grid points are set according to a preset resolution, and the serial number of each grid point is recorded as k (k =1, 2)Corresponding to a longitude and latitude of
Figure BDA00038909067600001211
Calculate the distance from point k to each observation point i as
Figure BDA00038909067600001212
Weighted average of k points
Figure BDA0003890906760000131
Figure BDA0003890906760000132
R k The calculation method is similar to the aforementioned method.
By
Figure BDA0003890906760000133
And
Figure BDA0003890906760000134
inverse calculation of amplitude ratio of each grid point
Figure BDA0003890906760000135
And difference of retardation angle
Figure BDA0003890906760000136
Figure BDA0003890906760000137
Figure BDA0003890906760000138
(III) calculating the equilibrium tidal Condition constant
For any observation point
Figure BDA0003890906760000139
Its corresponding equilibrium tidal height ζ EQ Can be represented as, for example,
Figure BDA00038909067600001310
in the formula, f i And u i Correcting angle, V, for the cross point factor and cross point of the partial tide 0 The phase of the ith partial tide of the equilibrium tide at the moment t = 0; omega i Is the tide division angular rate; p is a family number, p =1 for full tides and p =2 for half tides; λ is longitude; s = -8h is a standard time zone number of Beijing;
Figure BDA00038909067600001311
for the balanced tide amplitude after the correction of the ground tide for each partial tide,
Figure BDA00038909067600001312
Figure BDA00038909067600001313
Figure BDA00038909067600001314
Figure BDA00038909067600001315
in the formula (I), the compound is shown in the specification,
Figure BDA00038909067600001316
the latitude is.
(IV) calculating GPS data harmonic constants
Analyzing GPS observation data, according to the tide theory, assuming that zeta is the instantaneous displacement of sea tide load to one of east, north and radial directions of the observation station, the harmonic constant of the main tide can be expressed as,
Figure BDA0003890906760000141
in the formula, S 0 Is the average position of the displacement; h and g are the harmonic constants of the partial tides; f. of i And u i Correcting angle, V, for the cross point factor and cross point of partial tide 0 The phase of the ith partial tide of the balance tide at the time of t =0 is obtained; omega i Is the tide division angular rate; where f, ω t and v 0 + u is time dependent, location independent; and H and g are location dependent, time independent; the above formula can be written
Figure BDA0003890906760000142
Wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003890906760000143
Figure BDA0003890906760000144
according to the principle of least square method, solving X i And Y i Then, the harmonic constant H is calculated according to the following formula Gps And g Gps ，
Figure BDA0003890906760000145
In the actual calculation process, the period of the GPS satellite orbit (about 11.967 hours), the repetition period of the GPS constellation (about 23.9319 hours) and the period of partial tide (K) 2 The period is about 11.9672 hours, K 1 Period of about 23.9345 hours) are very similar and coupling phenomena tend to occur. Errors in the GPS satellite orbit are very easily absorbed by these tides, resulting in K 1 /K 2 The estimation result of the tide distribution load displacement has larger error.
(V) correcting K1 and K2 tide estimation errors by using relative admittance
For the same tide family, the main partial tides and the secondary partial tides are both present, the main partial tide proportion is large, the secondary partial tide proportion is small, but the tide characteristics of the partial tides in the same tide family are basically consistent.
Order to
Figure BDA0003890906760000151
In the formula, H Gpsi The amplitudes of the partial tides are obtained through GPS data harmonic analysis;
Figure BDA0003890906760000152
for each partial tide to balance tide amplitude psi i Is the modulus.
For the same tide family, the relative admittance characteristics of each tide family have similar characteristics, S for the semilunar tide family 2 To M 2 ，N 2 To M 2 ，K 2 To S 2 Relative admittance characteristics are similar, so M 2 、N 2 、S 2 And K 2 The R values corresponding to the partial tides are similar, the psi value of each semidiurnal partial tide is subjected to curve fitting, and K is obtained 2 The abnormal psi value of the tide is corrected to the normal fitting curve, thus solving the problem of K 2 The estimation result of the tide load displacement has a large error, as shown in fig. 2. For the same reason, for the whole climax group, O 1 To K 1 ，P 1 To K 1 And Q 1 To O 1 The relative admittance characteristics of the partial tides are similar, so K 1 、P 1 、O 1 And Q 1 The R values corresponding to the partial tides are similar, the psi value of each semidiurnal partial tide is subjected to curve fitting, and K is obtained 1 The abnormal psi value of the tide is corrected to the normal fitting curve, thus solving the problem of K 1 The estimation result of the tide distribution load displacement has larger error.
Example 2:
as shown in fig. 3, the present invention also provides a device for reducing tide separation error when researching sea tide load displacement by using GNSS, the device includes at least one processor and at least one memory, and also includes a communication interface and an internal bus; the memory having stored therein a computer program or instructions; when the processor executes the computer program or the instructions, the method for reducing the tide distribution error in the research of sea tide load displacement by using GNSS as described in embodiment 1 can be implemented. The internal bus may be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, an Extended ISA (EISA) bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, the buses in the figures of the present application are not limited to only one bus or one type of bus. The memory may include a high-speed RAM memory, and may further include a non-volatile storage NVM, such as at least one magnetic disk memory, and may also be a usb disk, a removable hard disk, a read-only memory, a magnetic disk or an optical disk.
The device may be provided as a terminal, server, or other form of device.
Fig. 3 is a block diagram of an apparatus shown for exemplary purposes. The device may include one or more of the following components: processing components, memory, power components, multimedia components, audio components, interfaces for input/output (I/O), sensor components, and communication components. The processing components typically control overall operation of the electronic device, such as operations associated with display, telephone calls, data communications, camera operations, and recording operations. The processing components may include one or more processors to execute instructions to perform all or a portion of the steps of the methods described above. Further, the processing component can include one or more modules that facilitate interaction between the processing component and other components. For example, the processing component may include a multimedia module to facilitate interaction between the multimedia component and the processing component.
The memory is configured to store various types of data to support operations at the electronic device. Examples of such data include instructions for any application or method operating on the electronic device, contact data, phonebook data, messages, pictures, videos, and the like. The memory may be implemented by any type or combination of volatile and non-volatile memory devices, such as Static Random Access Memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic or optical disks.
The power supply component provides power to various components of the electronic device. The power components may include a power management system, one or more power supplies, and other components associated with generating, managing, and distributing power for an electronic device. The multimedia component comprises a screen providing an output interface between said electronic device and a user. In some embodiments, the screen may include a Liquid Crystal Display (LCD) and a Touch Panel (TP). If the screen includes a touch panel, the screen may be implemented as a touch screen to receive an input signal from a user. The touch panel includes one or more touch sensors to sense touch, slide, and gestures on the touch panel. The touch sensor may not only sense the boundary of a touch or slide action, but also detect the duration and pressure associated with the touch or slide operation. In some embodiments, the multimedia component includes a front facing camera and/or a rear facing camera. The front camera and/or the rear camera may receive external multimedia data when the electronic device is in an operating mode, such as a shooting mode or a video mode. Each front camera and rear camera may be a fixed optical lens system or have a focal length and optical zoom capability.
The audio component is configured to output and/or input an audio signal. For example, the audio assembly includes a Microphone (MIC) configured to receive an external audio signal when the electronic device is in an operational mode, such as a call mode, a recording mode, and a voice recognition mode. The received audio signal may further be stored in a memory or transmitted via a communication component. In some embodiments, the audio assembly further comprises a speaker for outputting audio signals. The I/O interface provides an interface between the processing component and a peripheral interface module, which may be a keyboard, click wheel, button, etc. These buttons may include, but are not limited to: a home button, a volume button, a start button, and a lock button.
The sensor assembly includes one or more sensors for providing various aspects of status assessment for the electronic device. For example, the sensor assembly may detect an open/closed state of the electronic device, the relative positioning of the components, such as a display and keypad of the electronic device, the sensor assembly may also detect a change in the position of the electronic device or a component of the electronic device, the presence or absence of user contact with the electronic device, orientation or acceleration/deceleration of the electronic device, and a change in the temperature of the electronic device. The sensor assembly may include a proximity sensor configured to detect the presence of a nearby object in the absence of any physical contact. The sensor assembly may also include a light sensor, such as a CMOS or CCD image sensor, for use in imaging applications. In some embodiments, the sensor assembly may further include an acceleration sensor, a gyroscope sensor, a magnetic sensor, a pressure sensor, or a temperature sensor.
The communication component is configured to facilitate wired or wireless communication between the electronic device and other devices. The electronic device may access a wireless network based on a communication standard, such as WiFi,2G or 3G, or a combination thereof. In an exemplary embodiment, the communication component receives a broadcast signal or broadcast related information from an external broadcast management system via a broadcast channel. In one exemplary embodiment, the communication component further includes a Near Field Communication (NFC) module to facilitate short-range communications. For example, the NFC module may be implemented based on Radio Frequency Identification (RFID) technology, infrared data association (IrDA) technology, ultra Wideband (UWB) technology, bluetooth (BT) technology, and other technologies.
In an exemplary embodiment, the electronic device may be implemented by one or more Application Specific Integrated Circuits (ASICs), digital Signal Processors (DSPs), digital Signal Processing Devices (DSPDs), programmable Logic Devices (PLDs), field Programmable Gate Arrays (FPGAs), controllers, micro-controllers, microprocessors or other electronic components for performing the above-described methods.
Example 3:
the invention further provides a non-volatile computer-readable storage medium, in which a computer program or instructions are stored, and when the computer program or instructions are executed by a processor, the method for reducing tide distribution error when researching sea tide load displacement by using GNSS as described in embodiment 1 can be implemented.
In particular, a system, apparatus or device may be provided which is provided with a readable storage medium on which software program code implementing the functionality of any of the embodiments described above is stored and which causes a computer or processor of the system, apparatus or device to read out and execute instructions stored in the readable storage medium. In this case, the program code itself read from the readable medium can realize the functions of any of the above-described embodiments, and thus the machine-readable code and the readable storage medium storing the machine-readable code form part of the present invention.
The storage medium may be implemented by any type or combination of volatile or non-volatile memory devices such as Static Random Access Memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic or optical disks (e.g., CD-ROM, CD-R, CD-RW, DVD-20ROM, DVD-RAM, DVD-RW), magnetic tape, or the like. A storage media may be any available media that can be accessed by a general purpose or special purpose computer.
It should be understood that the Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), etc. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of a method disclosed in connection with the present invention may be embodied directly in a hardware processor, or in a combination of hardware and software modules.
It should be understood that a storage medium is coupled to the processor such that the processor can read information from, and write information to, the storage medium. Of course, the storage medium may also be integral to the processor. The processor and the storage medium may reside in an Application Specific Integrated Circuits (ASIC). Of course, the processor and the storage medium may reside as discrete components in a terminal or server.
The computer-readable program instructions described herein may be downloaded from a computer-readable storage medium to a respective computing/processing device, or to an external computer or external storage device over a network, such as the internet, a local area network, a wide area network, and/or a wireless network. The network may include copper transmission cables, fiber optic transmission, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. The network adapter card or network interface in each computing/processing device receives the computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in a computer-readable storage medium in the respective computing/processing device.
The computer program instructions for carrying out operations of the present disclosure may be assembler instructions, instruction Set Architecture (ISA) instructions, machine-related instructions, microcode, firmware instructions, state setting data, or source or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C + + or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The computer readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet service provider). In some embodiments, the electronic circuitry that can execute the computer-readable program instructions implements aspects of the present disclosure by utilizing the state information of the computer-readable program instructions to personalize the electronic circuitry, such as a programmable logic circuit, a Field Programmable Gate Array (FPGA), or a Programmable Logic Array (PLA).
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made to the present application by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Although the present invention has been described with reference to the specific embodiments, it should be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.

Claims (6)

1. A method for reducing tide separation errors in research on sea tide load displacement by utilizing GNSS is characterized by comprising the following steps:
step 1, determining a tide separating pair; the determined tide pair is determined as N 2 To M is aligned with 2 、S 2 To M is aligned with 2 、K 2 To M is aligned with 2 、K 2 To S 2 、P 1 To K 1 、O 1 To K 1 、Q 1 To K 1 And Q 1 To O 1 The tide dividing pair is adopted;
step 2, establishing the relative relation of each tide pair, and acquiring the amplitude ratio of each grid point
Figure FDA0003890906750000011
And difference in retardation
Figure FDA0003890906750000012
Step 3, calculating the equilibrium tide harmonic constant to obtain the equilibrium tide amplitude of each tide after the correction of the ground tide
Figure FDA0003890906750000013
Step 4, calculating a GPS data harmonic constant, and acquiring the amplitude H of each tide obtained by the harmonic analysis of the GPS data Gpsi (ii) a According to the tide theory, assuming that ζ is the instantaneous displacement of the sea tide load to one direction in the observation station, the harmonic constant of the main tide division can be expressed as,
Figure FDA0003890906750000014
in the formula, S 0 Is the average position of the displacement; h and g are the harmonic constants of the partial tides; f. of i And u i Correcting angle, V, for the cross point factor and cross point of the partial tide 0 The phase of the ith partial tide of the equilibrium tide at the moment t = 0; omega i Is the tide division angular rate; where f, ω t and v 0 + u is time dependent, location independent; and H and g are location dependent, time independent; the above formula can be written
Figure FDA0003890906750000015
Wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003890906750000016
Figure FDA0003890906750000017
according to the principle of least square method, solving X i And Y i Then, the harmonic constant H is calculated according to the following formula Gps And g Gps ，
Figure FDA0003890906750000021
Step 5, correcting K by using relative admittance according to the results of the step 2, the step 3 and the step 4 1 And K 2 And (4) a tide estimation error.
2. The method for reducing tide distribution error in the research of sea tide load displacement by using GNSS as claimed in claim 1, wherein the specific process of step 1 is as follows:
for a certain partial tide i, admittance M i Is composed of
Figure FDA0003890906750000022
Wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003890906750000023
H i and g i Harmonic constants, amplitudes and retardations, C, for the partial tides i For the tidal power coefficient, the relationship between the main partial tide n and the secondary partial tide m in the same tide family can be expressed as
Figure FDA0003890906750000024
Wherein
Figure FDA0003890906750000025
a m/n ＝g m -g n (4)
H in the formula (3) m /H n I.e. the amplitude ratio, g, in conventional tidal science m -g n For a difference in retardation, A * m/n With main partial tide n relative to secondary partial tide mAmplitude relation, a m/n The lag angle relation of the main partial tide n relative to the secondary partial tide m is shown; if the amplitude ratio of the actual partial tide is the same as the ratio of the induced tide force, the current tide is divided into a first tide and a second tide
Figure FDA0003890906750000026
Equal to 1; in the real ocean
Figure FDA0003890906750000027
It will not be exactly equal to 1, but it is close to 1 for most of the sea; if the two lag angles of partial tide are equal, then a m/n Is equal to zero; in the real ocean a m/n Generally not equal to zero, and is greater than zero or less than zero for most sea areas;
determining N 2 To M 2 、S 2 To M 2 、K 2 To M 2 、K 2 To S 2 、P 1 To K 1 、O 1 To K 1 、Q 1 To K 1 And Q 1 To O is 1 Dividing the tide pair;
Figure FDA0003890906750000031
g m ＝a m/n +g n (6)
Figure FDA0003890906750000032
for two partial tides, setting the amplitude ratio as H, the lag angle difference as g, and recording A = Hcosg and B = Hsing;
the specific process of the step 2 comprises the following steps:
it is known that there are N observation points, each represented by i (i =1,2.., N), each having a longitude and latitude of
Figure FDA0003890906750000033
Wherein
Figure FDA0003890906750000034
Negative values for south picks; the number of the lambda-beams is increased,
Figure FDA0003890906750000035
degree is taken as a unit;
Figure FDA0003890906750000036
and
Figure FDA0003890906750000037
the amplitude ratio and the lag angle difference of observation points of each tide testing station are obtained;
for observation station i, calculating the distance between the observation station i and other stations j
Figure FDA0003890906750000038
Weighted average of i points
Figure FDA0003890906750000039
Figure FDA00038909067500000310
Wherein the content of the first and second substances,
Figure FDA00038909067500000311
represents summing all observation points, but excluding the i point; r i Is the corrected Kelvin wave length, which is determined according to the water depth near the point i, and the method is as follows:
setting the water depth value of each grid point in and around the known sea area, namely the longitude and latitude of the point m
Figure FDA00038909067500000312
M =1,2, M, water depth D m (ii) a For the point i, calculating the distance between the grid point m and the tide checking station i as
Figure FDA0003890906750000041
Select all r im ＜150×10 3 m points, averaging the water depth values, and recording
Figure FDA00038909067500000411
Calculating corresponding Kelvin wave velocity
Figure FDA0003890906750000042
Kelvin wave length
L i ＝u i T (13)
Wherein u is i At a Kelvin wave velocity and T at a tidal cycle, i.e.
T＝360°/ω (14)
ω is the tidal angular rate;
finally, selecting
R i ＝κL i (κ＝1/4) (15)
For all i =1,2
Figure FDA0003890906750000043
And
Figure FDA0003890906750000044
and calculate
Figure FDA0003890906750000045
And root mean square deviation
Figure FDA0003890906750000046
If it is
δ i ≥mσ (18)
The data of the ith point is considered to be abnormal and needs to be omitted;
to verify the accuracy of the removed outliers, the method comprises
Figure FDA0003890906750000047
And
Figure FDA0003890906750000048
inversely calculating amplitude ratio of each observation point
Figure FDA0003890906750000049
And difference in retardation
Figure FDA00038909067500000410
Figure FDA0003890906750000051
Figure FDA0003890906750000052
Then by
Figure FDA0003890906750000053
Figure FDA0003890906750000054
Calculating the mean square error of the amplitude ratio and the retardation
Figure FDA0003890906750000055
Figure FDA0003890906750000056
Discarding unreasonable observed values according to the method described in formula (18), setting each grid point according to a predetermined resolution, and recording the serial number of each grid point as k (k =1,2.. Multidot.
Figure FDA0003890906750000057
Calculate the distance from point k to each observation point i as
Figure FDA0003890906750000058
Weighted average of k points
Figure FDA0003890906750000059
Figure FDA00038909067500000510
Wherein R is k Calculation method and R in the above formula (15) i The calculation method of (a) is the same,
by
Figure FDA0003890906750000061
And
Figure FDA0003890906750000062
inversely calculating the amplitude ratio of each grid point
Figure FDA0003890906750000063
And difference in retardation
Figure FDA0003890906750000064
Figure FDA0003890906750000065
Figure FDA0003890906750000066
3. The method for reducing tidal separation error in the research of sea tide load displacement using GNSS as claimed in claim 1, wherein the specific process of step 3 is:
for any observation point
Figure FDA0003890906750000067
Its corresponding equilibrium tidal height ζ EQ Can be expressed as
Figure FDA0003890906750000068
In the formula (f) i And u i Correcting angle, V, for the cross point factor and cross point of partial tide 0 The phase of the ith partial tide of the equilibrium tide at the moment t = 0; omega i Is the tide division angular rate; p is a family number, p =1 for full tides and p =2 for half tides; λ is longitude; s = -8h is a standard time zone number of Beijing;
Figure FDA0003890906750000069
for the balanced tide amplitude after the correction of the ground tide for each partial tide,
Figure FDA00038909067500000610
Figure FDA00038909067500000611
Figure FDA00038909067500000612
Figure FDA00038909067500000613
in the formula (I), the compound is shown in the specification,
Figure FDA00038909067500000614
the latitude is.
4. The method for reducing tidal separation error in the research of sea tide load displacement using GNSS as claimed in claim 1, wherein the specific process of step 5 is:
order to
Figure FDA00038909067500000615
In the formula, H Gpsi For the amplitude of each tide obtained from the GPS data harmonic analysis,
Figure FDA0003890906750000071
for each partial tide to balance tide amplitude psi i For modulus, divide each half-day tide by N 2 、M 2 、S 2 And K 2 Is curve fitted to the psi value of (A) and then K is fitted 2 Correcting the abnormal psi value of the tide division to a normal fitting curve, thereby completing K 2 Correction of moisture separation error, and then increase K 2 Accuracy of the damping and reconciliation constants; in the same way, divide each full day into tide K 1 、P 1 、O 1 And Q 1 Is curve fitted to the psi value of (A) and then K is fitted 1 Correcting the abnormal psi value of the tide division to a normal fitting curve, thereby completing K 1 Correction of tide-dividing error, and then raising K 1 Accuracy of the tide blending constant.
5. The utility model provides an utilize equipment that GNSS reduces and divides tide error when research sea tide load displacement which characterized in that: the apparatus comprises at least one processor and at least one memory, the processor and memory coupled; the memory having stored therein a computer program or instructions; the processor, when executing the computer program or instructions, may implement the method for reducing tidal separation error when studying sea tide load displacement according to any of claims 1 to 4.
6. A computer-readable storage medium, in which a computer program or instructions are stored, which, when executed by a processor, implement the method for reducing tidal separation error when studying sea tide load displacement using GNSS according to any of claims 1 to 4.
CN202211260841.7A 2022-08-10 2022-08-10 Method for reducing tide distribution error in research of sea tide load displacement by utilizing GNSS Pending CN115616630A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211260841.7A CN115616630A (en) 2022-08-10 2022-08-10 Method for reducing tide distribution error in research of sea tide load displacement by utilizing GNSS

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202211260841.7A CN115616630A (en) 2022-08-10 2022-08-10 Method for reducing tide distribution error in research of sea tide load displacement by utilizing GNSS
CN202210953013.5A CN115032670B (en) 2022-08-10 2022-08-10 Method and equipment for correcting tide distribution errors obtained by GPS observation

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN202210953013.5A Division CN115032670B (en) 2022-08-10 2022-08-10 Method and equipment for correcting tide distribution errors obtained by GPS observation

Publications (1)

Publication Number Publication Date
CN115616630A true CN115616630A (en) 2023-01-17

Family

ID=83130549

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202211260841.7A Pending CN115616630A (en) 2022-08-10 2022-08-10 Method for reducing tide distribution error in research of sea tide load displacement by utilizing GNSS
CN202210953013.5A Active CN115032670B (en) 2022-08-10 2022-08-10 Method and equipment for correcting tide distribution errors obtained by GPS observation

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202210953013.5A Active CN115032670B (en) 2022-08-10 2022-08-10 Method and equipment for correcting tide distribution errors obtained by GPS observation

Country Status (1)

Country Link
CN (2) CN115616630A (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117572467A (en) * 2023-10-18 2024-02-20 澳门大学 Multi-sea-condition tide level observation method and device, electronic equipment and storage medium

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20170017655A (en) * 2015-08-06 2017-02-15 대한민국(해양수산부 국립해양조사원장) Methods for production of tidal stream atlases using the Digital Tidal Stream Atlas S/W
CN106649991A (en) * 2016-11-15 2017-05-10 国家海洋局第海洋研究所 Calculation method for calculating secondary partial tide experience cotidal chart based on relative admittance
CN109614711B (en) * 2018-12-12 2022-10-21 国家海洋局第一海洋研究所 Method for improving calculation precision of multi-tide harmonic constant based on tidal current numerical simulation
CN113761808B (en) * 2021-11-08 2022-02-11 长沙理工大学 Surface tide displacement acquisition method and application method and system
CN114444330A (en) * 2022-04-02 2022-05-06 自然资源部第一海洋研究所 GNSS sea tide load displacement prediction method considering average effect
CN114662348B (en) * 2022-05-25 2022-09-06 山东省科学院海洋仪器仪表研究所 Method for establishing time-varying sea tide load displacement model based on GNSS technology

Also Published As

Publication number Publication date
CN115032670A (en) 2022-09-09
CN115032670B (en) 2022-11-18

Similar Documents

Publication Publication Date Title
CN111127563A (en) Combined calibration method and device, electronic equipment and storage medium
RU2644076C1 (en) Method of transmission numerical values, terminal and cloud server
CN114814872A (en) Pose determination method and device, electronic equipment and storage medium
US20130027427A1 (en) Associating Digital Images with Waypoints
US20210109228A1 (en) Identifying gnss navigation data as potentially manipulated or as trustworthy at least partially based on an estimated deviation of a second estimate of a satellite state from a first estimate of the satellite state
CN107272038B (en) High-precision positioning method and device
CN115616630A (en) Method for reducing tide distribution error in research of sea tide load displacement by utilizing GNSS
CN105635570B (en) Shooting preview method and system
CN106383355A (en) GNSS time difference monitoring method and system based on high precision time
CN111522045A (en) GPS module cold start method and device of terminal, terminal and storage medium
CN112488783A (en) Image acquisition method and device and electronic equipment
WO2022057304A1 (en) Data correction method and apparatus, terminal device and storage medium
Cortés et al. Performance evaluation of adaptive tracking techniques with direct-state Kalman filter
US20160308984A1 (en) Location accurate mobile events and social content
CN109813300B (en) Positioning method and terminal equipment
EP2177928A1 (en) Method for providing indoor positioning information to a navigation application based on the conversion of the positioning information to a standard format
CN110865170B (en) Method and device for determining sand content and storage medium
Li et al. Unified astronomical positioning and orientation model based on robust estimation
CN116222483A (en) Valve opening detection and correction method and device
CN116088011A (en) Quick calibration method and device for GNSS-R sea surface measurement height and electronic equipment
JP2012211840A (en) Position information correction system
CN116125495A (en) Ionosphere correction determination method, ionosphere correction determination device, ionosphere correction determination storage medium, and ionosphere correction determination program product
CN115830175A (en) Method and device for generating bending angle profile, electronic equipment and readable storage medium
CN115495038A (en) Audio adjusting method and device, electronic equipment and readable storage medium
CN110749912B (en) Cycle slip detection method, device, equipment and storage medium

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