CN109212524A - It is a kind of consider correcting value error correlation satellite altimetry precision determine method - Google Patents

It is a kind of consider correcting value error correlation satellite altimetry precision determine method Download PDF

Info

Publication number
CN109212524A
CN109212524A CN201811242023.8A CN201811242023A CN109212524A CN 109212524 A CN109212524 A CN 109212524A CN 201811242023 A CN201811242023 A CN 201811242023A CN 109212524 A CN109212524 A CN 109212524A
Authority
CN
China
Prior art keywords
correcting value
sea
correcting
value
dry
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201811242023.8A
Other languages
Chinese (zh)
Other versions
CN109212524B (en
Inventor
李洋
韩浪
袁仕耿
王韶波
陈娟
张严
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Aerospace Dongfanghong Satellite Co Ltd
Original Assignee
Aerospace Dongfanghong Satellite Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Aerospace Dongfanghong Satellite Co Ltd filed Critical Aerospace Dongfanghong Satellite Co Ltd
Priority to CN201811242023.8A priority Critical patent/CN109212524B/en
Publication of CN109212524A publication Critical patent/CN109212524A/en
Application granted granted Critical
Publication of CN109212524B publication Critical patent/CN109212524B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/882—Radar or analogous systems specially adapted for specific applications for altimeters
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating

Landscapes

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

Abstract

It is a kind of consider correcting value error correlation satellite altimetry precision determine method, the analytic modell analytical model determined first according to sea high-precision, determine the high observed quantity in sea, correcting value, and obtain the variance of correcting value, then calculating observation amount, the covariance of correcting value, and the covariance of wantonly two correcting values, finally obtained variance, covariance are summed and extracted square root, and then show that the sea after considering error correlation is high.Correlated error between the high observed quantity in the abundant analytical calculation sea of the present invention and correcting value and correcting value, the high precision analysis equation in sea is obtained based on this, compared with the analysis method for ignoring error correlation in the prior art, sea high-precision increases, and has good use value.

Description

It is a kind of consider correcting value error correlation satellite altimetry precision determine method
Technical field
The present invention relates to space technology field, especially a kind of satellite altimetry precision for considering correcting value error correlation is true Determine method.
Background technique
The high meaning in sea is elevation distance of the mean sea level relative to reference ellipsoid.The high-precision high information in sea is The subjects such as geodesy, oceanography bring revolutionary variation.Satellite altimetry utilizes precise orbit determination system and spaceborne Radar altimeter measures satellite to the distance and satellite of reference ellipsoid to the distance on sea respectively, and the two subtracts each other to obtain sea Elevation of the face relative to reference ellipsoid.The measurement data that satellite altimetry obtains has uniform spatial distribution, nearly whole world sea The features such as domain covering, measuring basis are unified, it has also become obtain the most effective measurement method of the high information in sea.
During satellite altimetry measures sea height, is also influenced and caused by factors such as atmosphere, marine environment These deviations need to be corrected by measured deviation.Current main correcting value has: sea situation deviation correcting value, inverse air pressure correcting value, Dry tropospheric correction amount, moist convection layer correcting value, ionosphere corrections amount.Therefore, it needs to consider in the high precision determination in sea above-mentioned The error of correcting value.But the high precision analytical method in sea has the disadvantage that in existing literature
(1) think that the error source in the determination of sea height is independent from each other in existing literature, have ignored correcting value and observation Correlativity that may be present between amount and correcting value.It is stepped up as sea is high-precision, ignored correlated error can It can cause the deviation of precision analysis.
(2) existing literature relates generally to correcting value and the high error propagation relationship in sea, in correcting value determination process Error propagation relationship lack research.For analytical error correlation, need to error source to correcting value and propagation characteristic grind Study carefully.
Summary of the invention
Technical problem solved by the present invention is having overcome the deficiencies of the prior art and provide a kind of consideration correcting value error phase The satellite altimetry precision of closing property determines method, overcomes the shortcomings of to ignore each correcting value error correlation in the prior art, provides one Kind surveys high model based on high-precision satellite, accurate to sea high-precision true in the case where fully considering correcting value error correlation It is fixed.
The technical solution of the invention is as follows: a kind of satellite altimetry precision determination side considering correcting value error correlation Method includes the following steps:
(1) analytic modell analytical model determined according to sea high-precision, determines the high observed quantity in sea, N correcting values；
(2) variance of N correcting values is determined respectively；
(3) judge whether observed quantity, correcting value have identical or associated variable, if so, the covariance of the two is then calculated, it is no Then without operation, all correcting values are traversed；
(4) judge that wantonly two correcting values are with the presence or absence of identical or associated variable in N correcting values, and if it exists, calculate The covariance of the two traverses all correcting values otherwise without operation；
(5) variance, covariance that step (2), (3), (4) obtain are summed and is extracted square root, and then show that consideration misses Sea after poor correlation is high.
The analytic modell analytical model that the described sea high-precision determines is
hSSH=rOrb-rAlt+rSSB+rIB+rDry+rWet+rIon
Wherein, hSSHFor sea height, rOrbFor orbit altitude, rAltFor altimeter distance measurement value, rSSBFor sea situation deviation correcting value, rIBFor inverse air pressure correcting value, rDryTo do tropospheric correction amount, rWetFor moist convection layer correcting value rWet, rIonFor ionosphere corrections Amount.
The N item correcting value includes sea situation deviation correcting value rSSB, inverse air pressure correcting value, dry tropospheric correction amount rDry, Moist convection layer correcting value rWet, ionosphere corrections amount rIon。
The variance of the sea situation deviation correcting value includes the sea deviation correcting value σ of Ku wave band2(rSSB), the sea of C-band Face deviation correcting value σ2(rSSB,C), in which:
σ2(rSSB)=(1.758)2σ2(HSW)+(0.5195)2σ2(UWS)
σ2(rSSB,C)=(2.5897)2σ2(HSW)+(0.8658)2σ2(UWS)
HSWFor sea significant wave height, UWSFor ocean surface wind speed.
The inverse air pressure correcting value rIBVariance is σ2(rIB)=0.99482σ2(PAtm)
Wherein, PAtmFor ocean surface atmospheric pressure.
The dry tropospheric correction amount rDryVariance is σ2(rDry)=0.22772σ2(PAtm)。
The ionosphere corrections amount rIonVariance is
σ2(rIon)=(0.333)2σ2(rAlt)+(0.14712)2σ2(HSW)+(0.17441)2σ2(UWS)。
- 2 σ (r of covariance of the ionosphere corrections amount and altimeter distance measurement valueAlt,rIon) be
-2σ(rAlt,rIonThe σ of)=0.35182(rAlt)
Wherein,It is altitude measuring away from mean value,For ionosphere corrections mean value.
The covariance of the ionosphere corrections amount and sea situation deviation correcting value is
Wherein, HSWFor significant wave height, UWSFor ocean surface wind speed,For significant wave height mean value,For ocean surface wind speed Mean value.
2 σ (the r of covariance of the dry tropospheric correction amount and inverse air pressure correcting valueDry,rIB) be
2σ(rDry,rIBThe σ of)=- 0.45302(PAtm)。
The advantages of the present invention over the prior art are that:
High accuracy analysis method in sea proposed by the present invention, the abundant high observed quantity in analytical calculation sea and correcting value and school Correlated error between positive quantity obtains the high precision analysis equation in sea based on this, in the prior art with ignore correlated error Analysis method compare, sea high-precision increase, have good use value.
Detailed description of the invention
Fig. 1 is the flow chart of the method for the present invention.
Specific embodiment
The present invention proposes a kind of satellite altimetry precision analysis method for considering correcting value error correlation, abundant analytical calculation Correlated error between the high observed quantity in sea and correcting value and correcting value, obtains the high precision analysis side in sea based on this Journey, in the prior art compared with the analysis method for ignoring correlated error, sea high-precision increases, and has and uses valence well Value, includes the following steps:
(1) analytic modell analytical model is determined according to sea high-precision, determines the high observed quantity in sea, N correcting value correcting values；
(2) model is determined using i-th (1≤i≤N) correcting value, application error communication theory analytical analysis i-th (1≤ I≤N) correcting value variance；
(3) step (2) are repeated, i traverses N from 1, exports the variance of N correcting values altogether；
(4) model is determined using i-th (1≤i≤N) correcting value, it is identical or associated whether judgement has with observed quantity Variable, if so, calculating the covariance of the two；If it is not, without operation；
(5) step (4) are repeated, i traverses N from 1, exports all correcting values covariance two-by-two relevant to observed quantity；
(6) using i-th and jth item (1≤i≤N, 1≤j≤N, i ≠ j) correcting value cover half type really, judge above-mentioned two With the presence or absence of identical or associated variable in model, if so, calculating the covariance of the two；If it is not, without operation；
(7) step (6) are repeated, i traverses N from 1, and j traverses N, i ≠ j from 1, exports all correcting values covariance two-by-two；
(8) result in step (3), (5) and (7) is summed and is extracted square root, obtained and consider that the sea of error correlation is high-precision Degree.Detailed explanation and illustration is carried out to the present invention with reference to the accompanying drawing.
1, the high observed quantity in clear sea and correcting value
Corresponding step (1) specifies the high observed quantity in sea and correcting value according to the analytic modell analytical model that sea high-precision determines, Share N correcting values.
The "high" that satellite ocean surveys senior middle school is that sea is high, i.e. elevation of the sea relative to reference ellipsoid.Its measurement process Are as follows: precise orbit determination technology and onboard GPS are utilized, the distance and satellite for measuring satellite to reference ellipsoid respectively arrive The distance on sea, the two subtract each other to obtain sea height.
Measurement process is influenced by sea environment and ocean-atmosphere interaction, is referred to as sea situation deviation (Sea State Bias) and it is inverse atmospheric pressure effect (Inverted Barometer).The electromagnetic wave of radar altimeter transmitting is by atmosphere Stemness gas, moist gas, electron content influence and generate ranging deviation, be referred to as dry troposphere (Wet Troposphere) transmission deviation, moist convection layer (Wet Troposphere) transmission deviation, ionosphere (Ionosphere) transmission Deviation.
The high accurate determination in sea need to correct above-mentioned bias term.The model that sea high-precision determines is established as a result, are as follows:
hSSH=rOrb-rAlt+rSSB+rIB+rDry+rWet+rIon
Wherein, hSSHIt is high (Sea Surface Height) for sea, there it can be seen that height determination in sea was related to Main observed quantity is orbit altitude rOrb, altimeter distance measurement value rAlt.Main correcting value is sea situation deviation correcting value rSSB, inverse air pressure Correcting value rIB, dry tropospheric correction amount rDry, moist convection layer correcting value rWet, ionosphere corrections amount rIon, totally 5 correcting value (N= 5)。
2, correcting value variance analysis
Corresponding step (2) and (3), the amount of being corrected variance analysis.
2.1 sea situation deviation correcting value rSSBVariance
Sea situation deviation correcting value rSSBComputation model are as follows:
rSSB=0.0066HSW+0.0015UWSHSW
rSSB,C=0.0074HSW+0.0025UWSHSW
Wherein rSSB、rSSB,CIt is the sea deviation correcting value of Ku, C-band, H respectivelySWFor sea significant wave height, UWSFor sea Wind speed.
Its variance is derived using above-mentioned sea situation buggy model.Total differential operation is carried out to sea situation buggy model, is had:
drSSB=(0.0066+0.0015UWS)dHSW+0.0015HSWdUWS
drSSB,C=(0.0074+0.0025UWS)dHSW+0.0025HSWdUWS
The measurement error of sea significant wave height and wind speed is independent from each other, then covariance sigma (the H of the twoSW,UWS)=0, The error propagation relationship of sea situation deviation correcting value and significant wave height and wind speed are as follows:
Engineering in practice, the error propagation coefficients C of sea situation deviation and significant wave height, ocean surface wind speed need to be calculatedSWH、CSWH,C、 CWS、CWS,CAverage statistical, have:
Probability density p (the U of ocean surface wind speed and significant wave heightWS) and p (HSW) weibull and log series model is respectively adopted, Have:
Using coefficient a=1.981, the b=0.131 obtained in document to Global Sea-level Wind Speed and Significant Wave Height statistics, μ= 0.326, σ=0.585 obtains E [CSWH]=1.758, E [CSWH,C]=2.5897, E [CWS]=0.5195, E [CWS,C]= 0.8658。
To sum up, the variance of sea situation deviation correcting value are as follows:
σ2(rSSB)=(1.758)2σ2(HSW)+(0.5195)2σ2(UWS)
σ2(rSSB,C)=(2.5897)2σ2(HSW)+(0.8658)2σ2(UWS)
Above formula is that the variance of sea situation deviation correcting value determines model, wherein the error of sea situation deviation correcting value and significant wave height The unit of propagation coefficient 1.758 and 2.5897 is cm/m, the error propagation coefficients 0.5195 of sea situation deviation correcting value and ocean surface wind speed Unit with 0.8658 is cm/m/s.
2.2 against air pressure correcting value rIBVariance
Inverse air pressure correcting value determines model are as follows:
rIB=-0.9948 (PAtm-1013.3)
Wherein PAtmFor ocean surface atmospheric pressure, unit mbar, rIBUnit is cm.
Model, which carries out total differential operation, to be determined to inverse air pressure correcting value:
drIB=0.9948dPAtm
Thus inverse air pressure correcting value variance is obtained are as follows:
σ2(rIB)=0.99482σ2(PAtm)
0.9948 unit of error propagation coefficients of inverse air pressure correcting value and atmospheric pressure is cm/mbar in above formula.
2.3 dry tropospheric correction amount rDryVariance
Dry tropospheric correction amount determines model are as follows:
rDry=0.2277PAtm
Wherein rDryUnit is cm, PAtmFor ocean surface atmospheric pressure, unit mbar.
Model, which carries out total differential operation, to be determined to dry tropospheric correction amount:
drDry=0.2277dPAtm
Thus dry tropospheric correction amount variance is obtained are as follows:
σ2(rDry)=0.22772σ2(PAtm)
Above formula is that the variance of dry tropospheric correction amount determines model, wherein dry tropospheric correction amount and the error of atmospheric pressure pass Broadcasting 0.2277 unit of coefficient is cm/mbar.
2.4 moist convection layer correcting value rWetVariance
Moist convection layer correcting value is generally measured using satellite-borne microwave radiometer, the index of variance and microwave radiometer Bright temperature measurement precision is related, uncorrelated to other correcting values, is denoted as σ2(rWet)。
2.5 ionosphere corrections amount rIonVariance
Ionosphere corrections amount and wave frequency square inversely, carrying out range measurement on both frequencies can be with It is estimated.Wherein, after the distance measure of two frequencies need to first carry out sea situation deviation correction, it is inclined ionosphere could to be used for The calculating of poor correcting value, has:
R=rAlt-rIon-rSSB
R=rAlt,C-rIon,C-rSSB,C
Wherein, r measures true value, r to sea for satelliteAlt、rAlt,C、rIon、rIon,C、rSSB、rSSB,CRespectively Ku and C-band Distance measurement value, ionosphere corrections value and sea situation deviation corrected value, ionosphere corrections amount be inversely proportional with frequency square, have:
K is constant, is had:
Two formulas subtract each other Ku wave band ionosphere corrections amount are as follows:
Sea situation deviation is determined that model substitutes into:
Above formula is ionosphere corrections amount cover half type really, with the distance measurement value of Ku and C-band, significant wave height, ocean surface wind speed phase It closes.Model, which carries out total differential operation, to be determined to ionosphere corrections amount, is had:
The distance measurement value of Ku and C-band, significant wave height, ocean surface wind speed error par wise irrelevance, have:
Wherein there are following relationships for the ranging variance of Ku and C-band:
Have:
When taking Ku frequency fKu=13.75GHz, fCThe error of=5.25GHz, altimeter distance measurement value and ionosphere corrections value passes Broadcast coefficient CrAlt=0.333.
Calculate ionosphere corrections amount and significant wave height, ocean surface wind speed error propagation coefficients CSWHAnd CWSAverage statistical side Method are as follows:
The wherein probability density p (U of ocean surface wind speed and significant wave heightWS) and p (HSW) determined in model with sea situation offset error Method, obtain E [CSWH]=0.14712, E [CWS]=0.17441.
In conclusion ionosphere corrections amount variance are as follows:
σ2(rIon)=(0.333)2σ2(rAlt)+(0.14712)2σ2(HSW)+(0.17441)2σ2(UWS)
Above formula is ionosphere corrections amount variance cover half type really, wherein altitude measuring away from the error of ionosphere corrections amount biography Broadcast 0.333 dimensionless of coefficient, the unit of the error propagation coefficients 0.14712 of significant wave height and ionosphere corrections amount is cm/m, sea The unit of the error propagation coefficients 0.17441 of face wind speed and ionosphere corrections amount is cms/m.
3, the analysis of covariance between correcting value and observed quantity
Corresponding step (4) and (5), the analysis of covariance between correcting value and observed quantity.
By each correcting value really cover half type it is found that only ionosphere corrections amount determines and is related to observed quantity altimeter in model Distance measurement value, other correcting values determine the variable being not directed in model with observed quantity.Therefore, only need to calculate ionosphere corrections amount with The covariance of altimeter distance measurement value.
Determine that model is known by sea high-precision, the covariance of ionosphere corrections amount and altimeter distance measurement value is -2 σ (rAlt, rIon), have:
WhereinRespectively altitude measuring has away from mean value and ionosphere corrections mean value:
NoteHave:
Altitude measuring is away from error delta rAltWith significant wave height measurement error Δ HSWIt is independent from each other, altimeter range error ΔrAltWith significant wave height measured value HSW, ocean surface wind speed measured value UWSIt is independent from each other.Δ r in above formulaAltMean valueThe 2nd in above formula and the 3rd are as follows: E (Δ rAltΔHSW)=E (ΔrAlt)E(ΔHSW)=0
Have:
Same E [Δ rAltThe calculation method of]=0 acquires E [Δ rAlt,C]=0 has using the property of variance:
There are following relationships for the ranging variance of Ku and C-band:
Have:
It is obtained by mean value propertiesHave:
As E (Δ rAltΔrAlt,CThe σ of)=1.6082(rAlt) when, have:
As E (Δ rAltΔrAlt,CThe σ of)=- 1.6082(rAlt) when, have:
Due toMiddle Δ rAlt,CWith Δ rAltA possibility that two kinds of value relationships is identical , therefore two kinds of -2 σ of expression formula of ionosphere corrections amount and the covariance of altimeter distance measurement value1(rAlt,rIon) and -2 σ2(rAlt, rIon) probability be 0.5.The covariance of ionosphere corrections amount and observed quantity altimeter distance measurement value as a result, are as follows:
-2σ(rAlt,rIon)=- 0.5 × 2 σ1(rAlt,rIon)-0.5×2σ2(rAlt,rIon)
=0.3518 σ2(rAlt)
0.3518 dimensionless of above formula coefficient.
4, the analysis of covariance between correcting value
It corresponds to step (6) and (7), the analysis of covariance between correcting value.
By each correcting value really cover half type it is found that sea situation deviation determine model and ionosphere corrections amount determine model all with have It is related to imitate wave height, ocean surface wind speed, inverse air pressure correcting value determine model and dry tropospheric correction amount determine model all with atmospheric pressure phase It closes, therefore the 2 σ (r of covariance of ionosphere corrections amount and sea situation deviation correcting value need to be calculatedIon,rSSB) and inverse air pressure correcting value With the 2 σ (r of covariance of dry tropospheric correction amountIB,rDry)。
The analysis of covariance of 4.1 ionosphere corrections amounts and sea situation deviation correcting value
2 σ (r of ionosphere corrections amount and the covariance of sea situation deviation correcting valueIon,rSSB) are as follows:
Due to Δ rAlt-ΔrAlt,C、ΔHSW、HSW、UWSIt is mutually indepedent two-by-two, E [Δ rAlt-ΔrAlt,C]=0, E [Δ HSW] =0, HSWFor sea significant wave height,For sea significant wave height mean value,For the mean value of ocean surface wind speed, have:
The 1st and the 2nd in above formula is respectively as follows:
The covariance of ionosphere corrections amount and sea situation deviation correcting value are as follows:
The order of magnitude of the 1st and the 2nd term coefficient is 10 in above formula-6With 10-7, illustrate ionosphere corrections amount and sea situation deviation The error correlation very little of correcting value can be neglected in the high error determination process in cm grades of seas.
The analysis of covariance of 4.2 dry tropospheric correction amounts and inverse air pressure correcting value
2 σ (the r of covariance of dry tropospheric correction amount and inverse air pressure correcting valueDry,rIB) are as follows:
WhereinThe measurement mean value of respectively dry tropospheric correction amount and inverse air pressure correcting value, has:
Δ P in above formulaAtmMean valueHave:
:
2σ(rDry,rIBThe σ of)=- 0.45302(PAtm)
Above formula is the covariance of dry tropospheric correction amount and inverse air pressure correcting value, and wherein the unit of coefficient -0.4530 is cm2/mbar2。
5, sea high accuracy analysis
Corresponding step (8), the variance that above-mentioned steps are calculated, covariance summation extraction of square root, obtain consideration correcting value The sea high-precision of error correlation, has:
σ(hSSH)=[σ2(rOrb)+σ2(rAlt)+σ2(rSSB)+σ2(rIB)+σ2(rDry)+σ2(rWet)+σ2(rIon)
+2σ(rIB,rDry)-2σ(rAlt,rIon)]0.5
=[σ2(rOrb)+σ2(rAlt)+(1.758)2σ2(HSW)+(0.5195)2σ2(UWS)
+(0.9948)2σ2(PAtm)+(0.2277)2σ(PAtm)+σ2(rWet)
+(0.333)2σ2(rAlt)+(0.14712)2σ2(HSW)+(0.17441)2σ2(UWS)
-0.4530σ2(PAtm)+0.3518σ2(rAlt)]0.5
=[σ2(rOrb)+(1.2094)2σ2(rAlt)+(1.7641)2σ2(HSW)+(0.5480)2σ2(UWS)
+(0.7671)2σ2(PAtm)+σ2(rWet)]0.5
Wherein 1.2094 dimensionless of error propagation coefficients of sea height and altimeter distance measurement value is passed with the error of significant wave height Broadcasting 1.764 unit of coefficient is cm/m, and the unit with the error propagation coefficients 0.5480 of ocean surface wind speed is cms/m, with atmospheric pressure 0.7671 unit of error propagation coefficients be cm/mbar.As can be seen from the above equation, radial track height error, altitude measuring away from Error, significant wave height error, ocean surface wind speed error, atmospheric pressure error and moist convection layer correction error are that influence sea high-precision is true Fixed independent error source.
By taking TOPEX/Poseidon satellite (abbreviation T/P satellite) and Jason-2 satellite data as an example, sea high-precision is influenced Determining independent error source is shown in Table 1.By the independent error source in table 1, using above-mentioned error analysis method, obtains and be related to sea height The precision of determining 2 observed quantities and 5 correcting values, covariance and sea high-precision, and do not consider error in conventional method The sea high-precision of correlation is compared, and is shown in Table 2.
From Table 2, it can be seen that comparison the method for the present invention and traditional sea high accuracy analysis method, T/P satellite and 6.5% and 11.1% has been respectively increased in the sea high-precision of Jason-2 satellite, shows using after the method for the present invention, not only sea High error separation more accurateization, but also sea high-precision is improved to a certain extent.
The independent error source table that 1 sea height of table determines
2 sea high accuracy analysis method comparison sheet (unit: cm) of table
The content that description in the present invention is not described in detail belongs to the well-known technique of those skilled in the art.

Claims (10)

1. a kind of satellite altimetry precision for considering correcting value error correlation determines method, it is characterised in that include the following steps:
(1) analytic modell analytical model determined according to sea high-precision, determines the high observed quantity in sea, N correcting values；
(2) variance of N correcting values is determined respectively；
(3) judge whether observed quantity, correcting value have identical or associated variable, if so, the covariance both then calculated, otherwise without Operation, traverses all correcting values；
Both (4) judge that wantonly two correcting values are with the presence or absence of identical or associated variable in N correcting values, and if it exists, calculate Covariance traverse all correcting values otherwise without operation；
(5) variance, covariance that step (2), (3), (4) obtain are summed and is extracted square root, and then obtained and consider error phase Sea after closing property is high.
2. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 1 determines method, special Sign is: the analytic modell analytical model that the sea high-precision determines is
hSSH=rOrb-rAlt+rSSB+rIB+rDry+rWet+rIon
Wherein, hSSHFor sea height, rOrbFor orbit altitude, rAltFor altimeter distance measurement value, rSSBFor sea situation deviation correcting value, rIBFor Inverse air pressure correcting value, rDryTo do tropospheric correction amount, rWetFor moist convection layer correcting value rWet, rIonFor ionosphere corrections amount.
3. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 2 determines method, special Sign is: the N item correcting value includes sea situation deviation correcting value rSSB, inverse air pressure correcting value, dry tropospheric correction amount rDry, wet Tropospheric correction amount rWet, ionosphere corrections amount rIon。
4. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 1 determines method, special Sign is: the variance of the sea situation deviation correcting value includes the sea deviation correcting value σ of Ku wave band2(rSSB), the sea of C-band Face deviation correcting value σ2(rSSB,C), in which:
σ2(rSSB)=(1.758)2σ2(HSW)+(0.5195)2σ2(UWS)
σ2(rSSB,C)=(2.5897)2σ2(HSW)+(0.8658)2σ2(UWS)
HSWFor sea significant wave height, UWSFor ocean surface wind speed.
5. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 1 or 2 determines method, It is characterized in that: the inverse air pressure correcting value rIBVariance is σ2(rIB)=0.99482σ2(PAtm)
Wherein, PAtmFor ocean surface atmospheric pressure.
6. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 1 or 2 determines method, It is characterized in that: the dry tropospheric correction amount rDryVariance is σ2(rDry)=0.22772σ2(PAtm)。
7. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 1 or 2 determines method, It is characterized in that: the ionosphere corrections amount rIonVariance is
σ2(rIon)=(0.333)2σ2(rAlt)+(0.14712)2σ2(HSW)+(0.17441)2σ2(UWS)。
8. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 1 or 2 determines method, It is characterized in that: -2 σ (r of covariance of the ionosphere corrections amount and altimeter distance measurement valueAlt,rIon) be
-2σ(rAlt,rIonThe σ of)=0.35182(rAlt)
Wherein,It is altitude measuring away from mean value,For ionosphere corrections mean value.
9. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 1 or 2 determines method, Be characterized in that: the covariance of the ionosphere corrections amount and sea situation deviation correcting value is
Wherein, HSWFor significant wave height, UWSFor ocean surface wind speed,For significant wave height mean value,For the equal of ocean surface wind speed Value.
10. a kind of satellite altimetry precision for considering correcting value error correlation according to claim 1 or 2 determines method, It is characterized by: the 2 σ (r of covariance of the dry tropospheric correction amount and inverse air pressure correcting valueDry,rIB) be
2σ(rDry,rIBThe σ of)=- 0.45302(PAtm)。
CN201811242023.8A 2018-10-24 2018-10-24 Satellite height measurement precision determination method considering correction value error correlation Active CN109212524B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811242023.8A CN109212524B (en) 2018-10-24 2018-10-24 Satellite height measurement precision determination method considering correction value error correlation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811242023.8A CN109212524B (en) 2018-10-24 2018-10-24 Satellite height measurement precision determination method considering correction value error correlation

Publications (2)

Publication Number Publication Date
CN109212524A true CN109212524A (en) 2019-01-15
CN109212524B CN109212524B (en) 2021-07-13

Family

ID=64997007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811242023.8A Active CN109212524B (en) 2018-10-24 2018-10-24 Satellite height measurement precision determination method considering correction value error correlation

Country Status (1)

Country Link
CN (1) CN109212524B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110553609A (en) * 2019-10-21 2019-12-10 中国人民解放军61540部队 Method and system for determining sea depth measurement precision of deep sea pressure instrument

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6380886B1 (en) * 1998-10-15 2002-04-30 The United States Of America As Represented By The Secretary Of The Navy Consistent combination of altimeter data from multiple satellites
US20070194980A1 (en) * 2005-08-23 2007-08-23 Mitchell Douglas A System and Method for Removal of Sea-State Bias in Satellite Altimetry Data
CN103760537A (en) * 2014-01-23 2014-04-30 中国科学院遥感与数字地球研究所 Tide correction method based on satellite altimetry data
CN106990404A (en) * 2017-03-30 2017-07-28 南京信息工程大学 A kind of autoscale algorithm using X-band radar inverting sea wave height of navigating

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6380886B1 (en) * 1998-10-15 2002-04-30 The United States Of America As Represented By The Secretary Of The Navy Consistent combination of altimeter data from multiple satellites
US20070194980A1 (en) * 2005-08-23 2007-08-23 Mitchell Douglas A System and Method for Removal of Sea-State Bias in Satellite Altimetry Data
CN103760537A (en) * 2014-01-23 2014-04-30 中国科学院遥感与数字地球研究所 Tide correction method based on satellite altimetry data
CN106990404A (en) * 2017-03-30 2017-07-28 南京信息工程大学 A kind of autoscale algorithm using X-band radar inverting sea wave height of navigating

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NIGEL T. PENNA ET AL.: "Sea Surface Height Measurement Using a GNSS Wave Glider", 《GEOPHYSICAL RESEARCH LETTERS》 *
李洋等: "卫星海洋测高技术体制发展研究", 《航天器工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110553609A (en) * 2019-10-21 2019-12-10 中国人民解放军61540部队 Method and system for determining sea depth measurement precision of deep sea pressure instrument
CN110553609B (en) * 2019-10-21 2021-07-23 中国人民解放军61540部队 Method and system for determining sea depth measurement precision of deep sea pressure instrument

Also Published As

Publication number Publication date
CN109212524B (en) 2021-07-13

Similar Documents

Publication Publication Date Title
Iijima et al. Automated daily process for global ionospheric total electron content maps and satellite ocean altimeter ionospheric calibration based on Global Positioning System data
Fu et al. Determining the response of sea level to atmospheric pressure forcing using TOPEX/POSEIDON data
Rocken et al. Analysis and validation of GPS/MET data in the neutral atmosphere
Ao et al. Lower troposphere refractivity bias in GPS occultation retrievals
Puysségur et al. Tropospheric phase delay in interferometric synthetic aperture radar estimated from meteorological model and multispectral imagery
Ménard et al. Calibration of the TOPEX/POSEIDON altimeters at Lampedusa: Additional results at Harvest
Ao et al. Rising and setting GPS occultations by use of open‐loop tracking
Bock et al. Validation of precipitable water from ECMWF model analyses with GPS and radiosonde data during the MAP SOP
Labroue et al. First quality assessment of the Cryosat-2 altimetric system over ocean
Liu et al. Quantifying residual ionospheric errors in GNSS radio occultation bending angles based on ensembles of profiles from end-to-end simulations
Garner et al. Total electron content measurements in ionospheric physics
Maus et al. Improved ocean-geoid resolution from retracked ERS-1 satellite altimeter waveforms
Liu et al. A quality control procedure based on bending angle measurement uncertainty for radio occultation data assimilation in the tropical lower troposphere
Bai et al. Evaluation of atmospheric profiles derived from single-and zero-difference excess phase processing of BeiDou radio occultation data from the FY-3C GNOS mission
CN109212524A (en) It is a kind of consider correcting value error correlation satellite altimetry precision determine method
Krasnopolsky et al. A multi-parameter empirical ocean algorithm for SSM/I retrievals
Glover et al. Air–sea gas transfer velocity estimates from the Jason-1 and TOPEX altimeters: Prospects for a long-term global time series
Eastwood et al. Algorithm theoretical basis document for the OSI SAF global reprocessed sea ice concentration product
Tang et al. High-spatial-resolution mapping of precipitable water vapour using SAR interferograms, GPS observations and ERA-Interim reanalysis
Mao et al. Time series insar ionospheric delay estimation, correction, and ground deformation monitoring with reformulating range split-spectrum interferometry
Zhu et al. Independent validation of Jason-2/3 and HY-2B microwave radiometers using Chinese coastal GNSS
Feltz et al. A methodology for the validation of temperature profiles from hyperspectral infrared sounders using GPS radio occultation: Experience with AIRS and COSMIC
Minkwitz et al. Application of SWACI products as ionospheric correction for single-point positioning: a comparative study
CN112528213B (en) Global ionosphere total electron content multilayer analysis method based on low earth orbit satellite
Gendt Comparison of IGS troposphere estimations

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant