CN109407089A - A kind of spaceborne wide swath interference imaging altimeter baseline angle estimation method of bilateral view - Google Patents
A kind of spaceborne wide swath interference imaging altimeter baseline angle estimation method of bilateral view Download PDFInfo
- Publication number
- CN109407089A CN109407089A CN201811123044.8A CN201811123044A CN109407089A CN 109407089 A CN109407089 A CN 109407089A CN 201811123044 A CN201811123044 A CN 201811123044A CN 109407089 A CN109407089 A CN 109407089A
- Authority
- CN
- China
- Prior art keywords
- value
- vector
- baseline angle
- difference
- offing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
Abstract
The invention discloses a kind of spaceborne bilaterals to regard wide swath interference imaging altimeter baseline angle estimation method, which comprises according to the range and step-length of the baseline angle measured value setting baseline angle value of right side view, thus obtains several baseline angle values；The inverting sea level height under each baseline angle value respectively obtains left side offing height vector and right side offing height vector；Difference is done to left side offing height vector and right side offing height vector respectively, left side parallax is obtained and divides view difference vector on the right side of vector sum；The average value that left side parallax divides on the right side of vector sum in the swath for regarding difference vector is calculated separately, two difference curves are drawn；It takes any one difference curves progress x-axis symmetrically to invert and obtains reversion difference curves, reversion difference curves baseline angle value corresponding with the intersection point of another difference curves is calculated as the right side of altimeter and regards baseline angle estimated value, left side regards baseline angle estimated value and subtracts right side view baseline angle estimated value as π.
Description
Technical field
The present invention relates to interference imaging altimeter signal processing technology fields, and in particular to a kind of spaceborne wide swath of bilateral view
Interference imaging altimeter baseline angle estimation method.
Background technique
Interference imaging altimeter is a kind of New Type Radar altimeter developed in recent years, deviates zenith point using low-angle and sees
High-precision sea level altitude, it can be achieved that in tens kilometers of breadth even kilometer ranges up to a hundred is surveyed with short base line interference measurement method
Measurement.It is observed if regarded using left and right bilateral, swath width can be further added by one times.The height observed compared to traditional substar
Degree meter, the efficiency of oceanographic observation and monitoring can be greatly improved in interference imaging altimeter, and meet mesoscale and sub- mesoscale with
The observation requirements of lower marine environment.For interference imaging altimeter, baseline angle error is to influence the master of Absolute height accuracy
One of factor is wanted, is the prerequisite for realizing the measurement of high-precision sea level altitude to accurately measuring or estimating for baseline angle.And
Directly measurement means are unable to satisfy required precision at present, it is necessary to be estimated.Wide swath interference imaging altimeter is regarded for bilateral
For, bilateral observation increases swath width, but also brings exclusive advantage to baseline angle estimation, but at present also
Not just for the baseline angle estimation method of bilateral view interference imaging altimeter.Existing interference imaging altimeter baseline parameter
There are two main classes for estimation method, unilateral depending on that can use with bilateral depending on interference imaging altimeter.One kind is to rely on ground
The external calibration method at control point, this method all have higher requirements to the selection of ground control point, layout and measurement, implement
Often has any problem and cannot often revisit；It is another kind of be carried out based on altimeter observation data itself baseline length and baseline angle,
The Combined estimator of interferometric phase system deviation, this kind of methods are complex, operand is big.
Summary of the invention
It is an object of the invention to overcome above-mentioned technological deficiency, the baseline of wide swath interference imaging altimeter is regarded based on bilateral
Error of tilt proposes a kind of base just for bilateral view interference imaging altimeter the characteristics of the left and right sides is just phase inverse
Line inclination angle estimation method, this method are simple and effective.
To achieve the goals above, the present invention provides a kind of spaceborne bilaterals inclines depending on wide swath interference imaging altimeter baseline
Angular estimation method, which comprises
According to the range and step-length of the baseline angle measured value setting baseline angle value of right side view, several are thus obtained
Baseline angle value；
The inverting sea level height under each baseline angle value respectively obtains left side offing height vector and right side view sea
Face height vector；
Difference is done to left side offing height vector and right side offing height vector respectively, left side is obtained and regards difference vector
Difference vector is regarded with right side；
The average value that left side parallax divides on the right side of vector sum in the swath for regarding difference vector is calculated separately, two difference songs are drawn
Line；
Take any one difference curves to carry out x-axis and symmetrically invert to obtain reversion difference curves, calculate reversion difference curves and
The corresponding baseline angle value of the intersection point of another difference curves regards baseline angle estimated value as the right side of altimeter, altimeter
Left side regards baseline angle estimated value and subtracts right side view baseline angle estimated value as π.
As a kind of improvement of the above method, the method is specifically included:
Step 1) obtains the baseline angle measured value that regards of right side by satellite positioning device as αc, then baseline angle value αR's
Value range is set as [αc-N*Δα,αc+ N* Δ α], i.e. αRWith αcCentered on, by step-length of Δ α take 2N+1 value, N is nature
Number；
Step 2) inverting sea level height under each baseline angle value obtains 2N+1 left side offing height vectorWith 2N+1 right side offing height vectorWherein i=1,2 ... 2N+1, θjFor altimeter visual angle, j=
1 ... M, M are the sum of visual angle value；
Step 3) does difference to left side offing height vector and right side offing height vector respectively, obtains 2N left side
Offing height difference vectorWith 2N right side offing height difference vectorI.e.Wherein i=1,2 ... 2N, j=1 ... M；
Step 4) is with view angle thetajFor variable, to 2N left side offing height difference vectorWith 2N right side view sea
Face height difference vectorIt averages respectively, obtains 2N left side offing height difference valueWith 2N right side view sea
Face height difference value
Then using i as the value of x-axis,WithFor the value of y-axis, two difference curves are drawn respectively；
Step 5) takesDifference curves carry out x-axis and symmetrically invert and obtain reverse curve, calculate reverse curve andDifference
The corresponding baseline angle value of the intersection point of component curve regards baseline angle estimated value as the right side of altimeter, and the left side of altimeter regards base
Line inclination angle estimated value is that π subtracts right side view baseline angle estimated value；
Or:
It takesDifference curves carry out x-axis and symmetrically invert and obtain reverse curve, calculate reverse curve andDifference curves
The corresponding baseline angle value of intersection point regard baseline angle estimated value as the right side of altimeter, the left side of altimeter regards baseline angle
Estimated value is that π subtracts right side view baseline angle estimated value.
Present invention has an advantage that
Method of the invention does not depend on ground control point, only observes data using altimeter itself, method is simple, and can
Independent estimations baseline angle, estimated accuracy are high.
Detailed description of the invention
Fig. 1 is the high schematic illustration of interference imaging altitude measuring；
Fig. 2 Simulation results: the difference curves of left and right side view.
Specific embodiment
The present invention will be further described with reference to the accompanying drawing.
Fig. 1 is the high schematic illustration of interference imaging altitude measuring.The solution of object height h can use following formula subrepresentation:
Wherein α is baseline angle, and ξ is the angle of baseline and primary antenna sight, and θ is altimeter visual angle, and r is altimeter antenna
Range-to-go, H are orbit altitude, R0For the reference ellipsoid radius of target position.
Due to calculate the object height value under different baseline angles, h can be denoted as to the function of α:
H=f (α) (2)
The baseline angle α of left and right side viewLAnd αRThere is following relationship: αL=π-αR, when right side depending on baseline angle there are errors
αe, i.e. αR=αR0+αe, wherein αR0True baseline angle is regarded for right side, then the baseline angle α that left side regardsL=π-αR0-αe=
αL0-αe, wherein αL0True baseline angle is regarded for left side.It follows that baseline angle error is just opposite in the left and right sides
Amount.The object height h of left and right side viewR、hLIt is as follows respectively:
hR=f (αR)=f (αR0+αe) (3)
hL=f (αL)=f (π-αR)=f (αL0-αe) (4)
The first step of baseline angle estimation is setting independent variable value range, we are by αRAs independent variable.Pass through first
The means such as satellite positioning device obtain the baseline angle measured value α that right side regardsc, with αcCentered on α is setRZone of reasonableness and
Step-length.Using Δ α as step-length, with [αc-N*Δα,αc+ N* Δ α] it is that value range obtains 2N+1 αRValue, and calculate it is right therewith
The α answeredLValue.Next to each αRAnd αLCalculate corresponding hRAnd hL, then to hRAnd hLDerivation can obtain respectively
By upper two formula it is found that such as willIt is denoted as g (αe), thenIt can easily be seen that working as αe=0, i.e. baseline
When error of tilt is 0, following formula is set up
And because g (αe) it is monotonic function, so αe=0 is the unique solution of formula (7), at this time αR=αR0。
g(αe) monotonicity prove it is as follows:
To g (αe) again derivation obtain
(H+R is apparent from by the geometrical relationship of Fig. 10) cos θ-r > 0 and H+R0- rcos θ > 0, so g ' (αe) > 0, so g
(αe) it is monotonic function.
In actual treatment, aforementioned derivation operations are realized by seeking difference.By formula (5) and formula (6) it is found that difference
The result is that vector related with view angle theta, this step obtains 2N left side offing height difference vectorWith 2N right side
Offing height difference vectorM is the sum of visual angle value.To each sea
Face height difference vector is with view angle thetajIt averages for variable, obtains 2N left side offing height difference valueWith 2N right side
Offing height difference valueFinal step solve formula (7), concrete methods of realizing be byWithIt is drawn as curve respectively,
Two curves take any one to carry out symmetrical reversion about x-axis and obtain reverse curve, calculate reverse curve and another curve
Intersection point, its corresponding αROr αLValue be baseline angle accurate estimated value.
Fig. 2, which gives, carries out the emulation examination that bilateral regards wide swath interference imaging altimeter baseline angle estimation using this method
Test result.Main simulation parameter such as table 1:
The main simulation parameter of table 1
| Altimeter carrier frequency | 13.58GHz |
| Angular field of view | 1.5 °~8.5 ° |
| Satellite platform height | 800km |
| Baseline angle (right side view) | 90.15° |
Solid line is the difference curves that right side regards in Fig. 2, and dotted line is symmetrical curve of the difference curves of left side view about x-axis, x
Axis is the baseline angle α that right side regardsRVariation, the corresponding α of the intersection point of two curvesRValue be obtained baseline angle estimated value
It is 90.149919 °, compared to 90.15 ° of simulation parameter setting, evaluated error is only 8.1 × 10-5°
It should be noted last that the above examples are only used to illustrate the technical scheme of the present invention and are not limiting.Although ginseng
It is described the invention in detail according to embodiment, those skilled in the art should understand that, to technical side of the invention
Case is modified or replaced equivalently, and without departure from the spirit and scope of technical solution of the present invention, should all be covered in the present invention
Scope of the claims in.
Claims (2)
1. a kind of spaceborne bilateral regards wide swath interference imaging altimeter baseline angle estimation method, which comprises
According to the range and step-length of the baseline angle measured value setting baseline angle value of right side view, several baselines are thus obtained
Inclination value；
The inverting sea level height under each baseline angle value, respectively obtains left side offing height vector and right side offing is high
Spend vector；
Difference is done to left side offing height vector and right side offing height vector respectively, left side parallax is obtained and divides vector sum right
Side view difference vector；
The average value that left side parallax divides on the right side of vector sum in the swath for regarding difference vector is calculated separately, two difference curves are drawn；
Take any one difference curves to carry out x-axis and symmetrically invert to obtain reversion difference curves, calculate invert difference curves with it is another
The corresponding baseline angle value of intersection point of difference curves regards baseline angle estimated value, the left side of altimeter as the right side of altimeter
It regards baseline angle estimated value and subtracts right side view baseline angle estimated value as π.
2. spaceborne bilateral according to claim 1 regards wide swath interference imaging altimeter baseline angle estimation method, special
Sign is that the method specifically includes:
Step 1) obtains the baseline angle measured value that regards of right side by satellite positioning device as αc, then baseline angle value αRValue
Range is set as [αc-N*Δα,αc+ N* Δ α], i.e. αRWith αcCentered on, by step-length of Δ α take 2N+1 value, N is natural number；
Step 2) inverting sea level height under each baseline angle value obtains 2N+1 left side offing height vector
With 2N+1 right side offing height vectorWherein i=1,2 ... 2N+1, θjFor altimeter visual angle, j=1 ... M, M
For the sum of visual angle value；
Step 3) does difference to left side offing height vector and right side offing height vector respectively, obtains 2N left side view sea
Face height difference vectorWith 2N right side offing height difference vectorI.e.Wherein i=1,2 ... 2N, j=1 ... M；
Step 4) is with view angle thetajFor variable, to 2N left side offing height difference vectorIt is high with 2N right side offing
Spend difference vectorIt averages respectively, obtains 2N left side offing height difference valueIt is high with 2N right side offing
Spend difference value
Then using i as the value of x-axis,WithFor the value of y-axis, two difference curves are drawn respectively；
Step 5) takesOrDifference curves carry out x-axis and symmetrically invert and obtain reverse curve, calculate reverse curve with separately
The corresponding baseline angle value of the intersection point of one difference curves for not doing to invert regards baseline angle estimated value as the right side of altimeter,
The left side of altimeter regards baseline angle estimated value and subtracts right side view baseline angle estimated value as π.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811123044.8A CN109407089B (en) | 2018-09-26 | 2018-09-26 | A kind of spaceborne wide swath interference imaging altimeter baseline angle estimation method of bilateral view |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811123044.8A CN109407089B (en) | 2018-09-26 | 2018-09-26 | A kind of spaceborne wide swath interference imaging altimeter baseline angle estimation method of bilateral view |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109407089A true CN109407089A (en) | 2019-03-01 |
| CN109407089B CN109407089B (en) | 2019-09-24 |
Family
ID=65465195
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811123044.8A Active CN109407089B (en) | 2018-09-26 | 2018-09-26 | A kind of spaceborne wide swath interference imaging altimeter baseline angle estimation method of bilateral view |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109407089B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111045099A (en) * | 2019-12-27 | 2020-04-21 | 武汉大学 | Method for inverting ocean gravity field by imaging type altimeter data |
| CN113589282A (en) * | 2021-07-12 | 2021-11-02 | 中国科学院国家空间科学中心 | Method for removing flat ground effect of spaceborne interference imaging altimeter based on image domain transformation |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6362776B1 (en) * | 2000-02-04 | 2002-03-26 | Honeywell International Inc. | Precision radar altimeter with terrain feature coordinate location capability |
| CN103364766A (en) * | 2012-04-01 | 2013-10-23 | 中国科学院电子学研究所 | External calibration method of satellite-borne interferometric synthetic aperture radar (InSAR) system |
| CN107102329A (en) * | 2017-04-17 | 2017-08-29 | 王辉 | Ka wave bands TOPS InSAR systems and its method of work |
| CN108007476A (en) * | 2017-11-20 | 2018-05-08 | 中国科学院空间应用工程与技术中心 | The interference calibrating method and system of a kind of space-based Interferometric Radar Imaging Altimeter |
| CN108061891A (en) * | 2017-12-04 | 2018-05-22 | 上海无线电设备研究所 | A kind of interference SAR base linc vector estimating method at no control point |
-
2018
- 2018-09-26 CN CN201811123044.8A patent/CN109407089B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6362776B1 (en) * | 2000-02-04 | 2002-03-26 | Honeywell International Inc. | Precision radar altimeter with terrain feature coordinate location capability |
| CN103364766A (en) * | 2012-04-01 | 2013-10-23 | 中国科学院电子学研究所 | External calibration method of satellite-borne interferometric synthetic aperture radar (InSAR) system |
| CN107102329A (en) * | 2017-04-17 | 2017-08-29 | 王辉 | Ka wave bands TOPS InSAR systems and its method of work |
| CN108007476A (en) * | 2017-11-20 | 2018-05-08 | 中国科学院空间应用工程与技术中心 | The interference calibrating method and system of a kind of space-based Interferometric Radar Imaging Altimeter |
| CN108061891A (en) * | 2017-12-04 | 2018-05-22 | 上海无线电设备研究所 | A kind of interference SAR base linc vector estimating method at no control point |
Non-Patent Citations (2)
| Title |
|---|
| YUNHUA ZHANG ET AL.: "Experimental demonstration for the attitude measurement capability of interferometric radar altimeter", 《2014 15TH INTERNATIONAL RADAR SYMPOSIUM》 * |
| 徐永生 等: "美国新一代测高卫星SWOT_评述我国宽刈幅干涉卫星的发展借鉴", 《遥感技术与应用》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111045099A (en) * | 2019-12-27 | 2020-04-21 | 武汉大学 | Method for inverting ocean gravity field by imaging type altimeter data |
| CN113589282A (en) * | 2021-07-12 | 2021-11-02 | 中国科学院国家空间科学中心 | Method for removing flat ground effect of spaceborne interference imaging altimeter based on image domain transformation |
| CN113589282B (en) * | 2021-07-12 | 2023-08-08 | 中国科学院国家空间科学中心 | Method for removing ground effect by satellite-borne interference imaging altimeter based on image domain transformation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109407089B (en) | 2019-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106595583B (en) | RTK measurement receiver inclination measurement method | |
| CN106443722B (en) | A kind of antenna phase center variation detection method | |
| CN107167790A (en) | A kind of step scaling method of laser radar two based on Calibration Field | |
| CN108037374B (en) | Array antenna near field calibration method | |
| CN105424058A (en) | Digital camera projection center position precision calibration method based on photogrammetric technology | |
| CN109407089B (en) | A kind of spaceborne wide swath interference imaging altimeter baseline angle estimation method of bilateral view | |
| CN110045340A (en) | The coordinate measurement and calculation method of beacon antenna in Planar Phased Array Antenna | |
| CN104614736A (en) | Calibration method of GPS receiver | |
| CN108738132A (en) | A kind of three base station movement communication positioning methods based on TDOA | |
| CN107991691B (en) | Satellite navigation positioning accuracy verification equipment and method | |
| CN113884977A (en) | Correction method for one-dimensional interferometer direction finding cone effect | |
| CN103335632A (en) | High-precision high-frequency response platform inclined angle measuring device and method for measuring platform inclined angle by adopting same | |
| CN108458710A (en) | Pose measuring method | |
| CN108154535A (en) | Camera Calibration Method Based on Collimator | |
| CN107144815B (en) | Three-dimensional positioning method based on one-dimensional direction finding | |
| CN103727962A (en) | Precision calibration method of wide-field infrared photoelectric theodolite | |
| CN104316027B (en) | Composite level own-side height determination method | |
| CN102519469A (en) | Planetary vehicle positioning method based on computer vision and VLBI combined adjustment | |
| CN115792800A (en) | Grid search-based double-station three-dimensional cross positioning method | |
| CN106646364B (en) | A kind of method for measuring spatial location | |
| CN109283501A (en) | A kind of dimensional turntable baseline alignment methods | |
| Aslakson et al. | Use of shoran in geodetic control | |
| CN116953746B (en) | Method and device for orienting satellite navigation antenna based on single phase distortion | |
| CN117670991A (en) | Method for measuring monocular vision distance of short-distance target between high-orbit satellites | |
| Biagi et al. | Relative antenna calibration for mass-market GNSS receivers: A case study |
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 |