CN104678371A - Device for measuring sea surface height based on time-delay modification - Google Patents
Device for measuring sea surface height based on time-delay modification Download PDFInfo
- Publication number
- CN104678371A CN104678371A CN201510071278.2A CN201510071278A CN104678371A CN 104678371 A CN104678371 A CN 104678371A CN 201510071278 A CN201510071278 A CN 201510071278A CN 104678371 A CN104678371 A CN 104678371A
- Authority
- CN
- China
- Prior art keywords
- signal
- module
- input
- satellite
- processing module
- 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
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a device for measuring sea surface height based on time-delay modification. The device comprises a navigation receiver, a right-hand antenna, a left-hand antenna, a radio-frequency front end module, an intermediate frequency signal processing module and an inversion module. The navigation receiver is connected with the right-hand antenna, the intermediate frequency signal processing module and the inversion module, and the right-hand antenna and the left-hand antenna are independent with each other and connected with the radio-frequency front end module, the intermediate frequency signal processing module and the inversion module sequentially; a direct radio-frequency voltage signal transmitted from the right-hand antenna and a reflective radio-frequency voltage signal transmitted from the left-hand antenna are adopted as the input of the radio-frequency front end module, a direct digital intermediate frequency signal and the reflective digital intermediate frequency signal transmitted from the radio-frequency front end module are adopted as the input of the intermediate frequency signal processing module, time-delay doppler two-dimensional related power transmitted from the intermediate frequency signal processing module is adopted as the input of the inversion module, the direct radio-frequency voltage signal transmitted from the right-hand antenna is adopted as the input of the navigation receiver, and navigation information transmitted from the navigation receiver is adopted as the input of the intermediate frequency signal processing module and the inversion module.
Description
Technical field
The present invention relates to a kind of sea level height measurement mechanism based on time delay correction.Specifically, utilize in GNSS-R (Global Navigation Satellites System-Reflectmetry) ocean surface wind retrieving result due to time delay error that sea surface roughness causes, revise and survey due to the time delay error that sea surface roughness causes in high process, thus improve altimetry precision.
Background technology
GNSS-R, as a kind of new remote sensing, is subject to the extensive concern of Chinese scholars.At present, GNSS-R gradual perfection in theory, technology and data inversion etc.Ocean remote sensing based on GNSS-R is a major part of GNSS-R remote sensing fields, and its main development direction is Ocean Wind-field and sea level height.Up to now, in succession carry out the experiment repeatedly detecting Ocean Wind-field and sea level height in the world, theory and the algorithm of GNSS-R inverting Ocean Wind-field are shaped substantially.
In inverting Ocean Wind-field, theoretical related power model is comparatively ripe, adopts least square method to mate on the rear edge of theoretical related power curve, can estimate ocean surface wind speed with the rear edge of actual measurement related power curve.Using the wind speed be finally inversed by as auxiliary, wind direction of ocean surface can be estimated by calculating residual sum of squares (RSS).In inverting sea level height, the key of altimetry precision is the mistiming of accurate Calculation direct signal and the reflected signal correlated power peak through mirror-reflection.But due to the impact of sea surface roughness, there is diffuse reflection point around specular reflection point, the peak point position causing reflected signal related power offsets backward.If at this moment use the half-power point method measuring-signal of traditional tracking peak point or waveform to postpone, then can bring larger deviation.
In sea level height is measured, the measuring error caused for solution sea surface roughness lacks a kind of simple directly solution.
Summary of the invention
The object of this invention is to provide a kind of sea level height measurement mechanism based on time delay correction.This device not only can exact inversion ocean surface wind speed, wind direction, can also exact inversion sea level height.During inverting Ocean Wind-field, can obtain and the theoretical related power curve of surveying related power curve and matching.This theoretical related power curve is the function about time delay, therefrom can estimate the time delay of reflected signal peak value point relative to specular reflection point place, i.e. time delay error.The present invention utilizes this time delay error to revise the time delay that senior middle school is surveyed on sea, survey high real time by sea to postpone to deduct this time delay error and obtain the mistiming, mistiming will be multiplied by the light velocity, obtain path difference, and then complete elevation carrection according to the geometric relationship that GNSS satellite, receiver and specular reflection point are formed.
The present invention takes following technical scheme:
A kind of sea level height measurement mechanism based on time delay correction of the present invention, it is made up of a navigation neceiver, dextrorotation antenna, left-handed antenna, RF front-end module, an IF signal processing module and an inverting module.Position annexation between them is: navigation neceiver is connected with inverting module with dextrorotation antenna, IF signal processing module respectively, dextrorotation antenna and left-handed antenna independently of one another, and are connected with inverting module with RF front-end module, IF signal processing module together successively.Signal trend between them is: the direct projection radio frequency voltage signal DRF that dextrorotation antenna exports and the reflected radio voltage signal RRF that left-handed antenna exports is as the input of RF front-end module, the direct projection digital medium-frequency signal DDRF that RF front-end module exports and reflection digital intermediate-freuqncy signal RDIF is as the input of IF signal processing module, the delay-Doppler two-dimensional correlation power DDM (τ that IF signal processing module exports, f) as the input of inverting module, the direct projection radio frequency voltage signal DRF that dextrorotation antenna exports is as the input of navigation neceiver, the navigation information NAV_inf that navigation neceiver exports is as the input of IF signal processing module and inverting module.
Described navigation neceiver is used for the navigation information NAV_inf of output packet containing information such as the height of receiver and speed, elevation of satellite, position angles, as the input of IF signal processing module and inverting module.
Described dextrorotation antenna adopts commercial GNSS dextrorotation antenna for receiving the direct projection electromagnetic signal of Navsat, and this signal is changed into the input of direct projection radio frequency voltage signal DRF as radio-frequency front-end process.
Described left-handed antenna adopts the left-handed antenna of many array element array for receiving the reflected electromagnetic signal of Navsat, and this signal is changed into the input of reflected radio voltage signal RRF as radio-frequency front-end process.
The direct projection radio frequency voltage signal DRF received and reflected radio voltage signal RRF is down-converted to intermediate-freuqncy signal by described RF front-end module, and intermediate-freuqncy signal is converted into direct projection digital medium-frequency signal DDRF and reflection digital intermediate-freuqncy signal RDIF, as the input of IF signal processing module.
Described IF signal processing module is made up of direct projection passage and reflection channel.Relation is between the two: be connected with each other.IF signal processing module is as the input of inverting module.
This direct projection passage by visible star anticipation, catch and follow the tracks of three parts and form.Relation between three is: visible star anticipation, catch and be connected successively with following the tracks of.Visible star anticipation is predicted left-handed antenna satellite within the vision, reduces direct projection passage satellite acquisition scope, exports No. PRN, left-handed antenna satellite within the vision, be designated as (PRN
1..., PRN
n), as the input of catching; Catch (PRN
1..., PRN
n) in satellite search for, judge that satellite is whether true visible, if truly visible, then guestimate carried out to direct signal code phase and Doppler frequency, as the input of tracking; Follow the tracks of and accurate code phase and Doppler's estimation are carried out to the direct signal of acquisition success, be designated as (τ
di, f
di), as the input of reflection channel.
By selecting, star, compensation rate calculate this reflection channel, carrier wave occurs, pseudo-code occurs and relevant treatment five part forms.Relation between five is: select star to calculate with compensation rate and be connected, and compensation rate calculates and occurs to be connected with pseudo-code with carrier wave respectively, and carrier wave occurs and pseudo-code occurs independent of one another, and carrier wave occurs to occur all to be connected with relevant treatment with pseudo-code.Select star to select star to the satellite traced into according to elevation of satellite, select the satellite that elevation of satellite is maximum, its direct signal correspondence code phase place and Doppler's valuation are designated as (τ
d, f
d), as the input that compensation rate calculates; Compensation rate calculates the navigation information NAV_inf exported according to navigation neceiver and calculates compensation rate Δ τ and Δ f, by direct signal correspondence code phase place and Doppler's valuation (τ
d, f
d) obtain the local pseudo-code phase τ of reflection channel
d+ Δ τ and local carrier frequency f
d+ Δ f, as the input that pseudo-code occurs and carrier wave occurs; Pseudo-code occurs to produce local pseudo-code and is used for the stripping of reflected signal pseudo-code; Carrier wave occurs to produce local carrier and is used for the stripping of reflected signal carrier wave; Local carrier and local pseudo-code and reflected signal RDIF are carried out related operation by relevant treatment, export delay-Doppler two-dimensional correlation power DDM (τ, f).
Described inverting module is made up of ocean surface wind retrieving module and sea level height inverting module.Relation is between the two: be connected with each other.
This ocean surface wind retrieving module calculates five parts by the generation of theoretical model storehouse, noise reduction process, normalization, Waveform Matching and residual sum of squares (RSS) and forms.Relation between five is: the generation of theoretical model storehouse, noise reduction process, normalization, Waveform Matching and residual sum of squares (RSS) calculate and be connected successively.It is from DDM (τ that theoretical model storehouse generates, f) information such as the height of receiver, speed, elevation of satellite and position angle are extracted in the navigation information NAV_inf exported with navigation neceiver, information extraction is input in surface scattering signal correction power module, inputs initial wind speed v
0, calculate the surface scattering signal correction power under different wind friction velocity, generative theory related power waveform library; Noise reduction process carries out noise reduction process to actual measurement related power curve D DM (τ, f), and stress release treatment affects; Normalization be to noise reduction process after actual measurement related power and the theoretical related power waveform library of generation be normalized respectively; Waveform Matching be by normalization after actual measurement related power curve and theoretical related power curve adopt least square method to carry out Waveform Matching, export the wind speed v corresponding with surveying related power Curve Matching successful theoretical related power Curves; It is by given initial wind direction φ that residual sum of squares (RSS) calculates
0calculate the theory of many Navsats and the residual sum of squares (RSS) of actual measurement related power waveform, adjustment wind direction, exports the wind direction φ that residual sum of squares (RSS) local minimum is corresponding.
This sea level height inverting module is by noise reduction process, normalization, time-delay calculation, and time delay correction and high computational five part form.Relation between five is: noise reduction process, normalization, time-delay calculation, time delay correction are connected successively with high computational.Noise reduction process carries out noise reduction process to actual measurement related power curve D DM (τ, f), and stress release treatment affects; Normalization be to noise reduction process after actual measurement related power be normalized; Time-delay calculation is by direct signal correlated power peak point code phase and the reflected signal correlated power peak point place code phase computational reflect signal time delay value τ relative to direct signal
e1; Time delay correction is the time delay error τ caused by the theoretical related power curve estimation sea surface roughness that the match is successful
c, time delay value revises the delay inequality that senior middle school is surveyed on sea thus, by revised time delay value τ
eobtain the path delay of direct signal and reflected signal; High computational is from navigation information NAV_inf, extract the information such as position and elevation of satellite of navigation neceiver, can try to achieve the distance of navigation neceiver to the earth's core, thus can calculate sea level height H by navigation neceiver position.
Principle of the present invention and operating mode brief introduction as follows:
When utilizing GNSS-R inverting Ocean Wind-field, the time delay of reflected signal peak value point relative to specular reflection point is obtained from the theoretical related power curve that coupling obtains, revise due to the time delay error that sea surface roughness causes in inverting sea level height with this length of delay, thus improve altimetry precision.Dextrorotation antenna and left-handed antenna receive direct signal and reflected signal respectively and convert it into radio frequency voltage signal, export two-dimensional correlation powertrace through RF front-end module and IF signal processing module again, the navigation information exported by two-dimensional correlation powertrace and navigation neceiver carries out inverting to Ocean Wind-field and sea level height in inverting module.
The invention has the advantages that:
One. this device only needs receiving trap, reduces complicacy and the cost of device.
Two. this device can complete the inverting of Ocean Wind-field and sea level height, strengthens the monitoring capability to ocean.
Three. the time delay error value that this device ocean surface wind retrieving obtains directly revises the time delay error in sea level height measurement, reduces the operand in sea level height measurement.
Four. this device can reduce the error because sea surface roughness causes, and improves sea level height measuring accuracy.
Accompanying drawing explanation
Fig. 1 is the overall framework figure of measurement mechanism of the present invention.
Fig. 2 is the IF signal processing module by signal process block diagram of measurement mechanism of the present invention.
Fig. 3 is Ocean Wind-field and the sea level height inverting block diagram of measurement mechanism of the present invention.
In figure, symbol description is as follows: NAV_inf is the navigation information of the height comprising receiver that exports of navigation neceiver and the information such as speed, elevation of satellite, position angle, DRF is direct projection radio frequency voltage signal, RRF is reflected radio voltage signal, DDRF is direct projection digital medium-frequency signal, RDIF is reflection digital intermediate-freuqncy signal, and DDM (τ, f) is actual measurement delay-Doppler two-dimensional correlation powertrace, H is the sea level height needing to export, v
0for the initial ocean surface wind speed of input, φ
0for the initial wind direction of ocean surface of input, v is the ocean surface wind speed needing to export, and φ is the wind direction of ocean surface needing to export, (PRN
1..., PRN
n) be No. PRN, left-handed antenna satellite within the vision, (τ
di, f
di) be the code phase of the direct signal of satellite that traces into and Doppler's estimated value, (τ
d, f
d) be the code phase of the direct signal of the maximum satellite of elevation of satellite and Doppler's estimated value, τ
d+ Δ τ is the local pseudo-code phase of reflection channel, f
d+ Δ f is the local carrier frequency of reflection channel, τ
cfor the time delay error that sea surface roughness causes, τ
e1for revising the time delay value of front-reflection signal relative to direct signal, τ
efor revising the time delay value of back reflection signal relative to direct signal.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail.
If Fig. 1 is the overall construction drawing of this device.
Apparatus of the present invention are made up of a navigation neceiver, commercial GNSS dextrorotation antenna, the left-handed antenna of array element array more than one, RF front-end module, an IF signal processing module and an inverting module.Position annexation between them is: navigation neceiver is connected with inverting module with dextrorotation antenna, IF signal processing module respectively, dextrorotation antenna and left-handed antenna independently of one another, and are connected with inverting module with RF front-end module, IF signal processing module together successively.Signal trend between them is: the direct projection radio frequency voltage signal DRF that dextrorotation antenna exports and the reflected radio voltage signal RRF that left-handed antenna exports is as the input of RF front-end module, the direct projection digital medium-frequency signal DDRF that RF front-end module exports and reflection digital intermediate-freuqncy signal RDIF is as the input of IF signal processing module, the delay-Doppler two-dimensional correlation power DDM (τ that IF signal processing module exports, f) as the input of inverting module, the direct projection radio frequency voltage signal DRF that dextrorotation antenna exports is as the input of navigation neceiver, the navigation information NAV_inf that navigation neceiver exports is as the input of IF signal processing module and inverting module.Described navigation neceiver is used for the navigation information NAV_inf of output packet containing information such as the height of this receiver and speed, elevation of satellite, position angles, as the input of IF signal processing module and inverting module.This navigation neceiver output packet is containing the navigation information NAV_inf of the information such as the height of receiver and speed, elevation of satellite, position angle.
This signal for receiving the direct projection electromagnetic signal of Navsat, and is changed into the input of direct projection radio frequency voltage signal DRF as radio-frequency front-end process by this dextrorotation antenna.
This signal for receiving the reflected electromagnetic signal of Navsat, and is changed into the input of reflected radio voltage signal RRF as radio-frequency front-end process by this left-handed antenna.
This RF front-end module comprises direct signal radio-frequency front-end and reflected signal radio-frequency front-end.The direct projection radio frequency voltage signal DRF received and reflected radio voltage signal RRF is down-converted to intermediate-freuqncy signal by RF front-end module, and intermediate-freuqncy signal is converted into direct projection digital medium-frequency signal DDRF and reflection digital intermediate-freuqncy signal RDIF, as the input of IF signal processing module.
It is the signal transacting block diagram of the IF signal processing module of this device as Fig. 2.
IF signal processing module comprises direct projection passage and reflection channel.Direct projection passage mainly completes catching of direct signal and follows the tracks of, and time delay tracking obtained, Doppler's estimated value export the reference information as reflected signal process; Reflection channel mainly completes the relevant treatment of reflected signal, obtains delay-Doppler two-dimensional correlation power.
Direct projection passage by visible star anticipation, catch and follow the tracks of three parts and form.
The navigation information NAV_inf that visible star anticipation part receives navigation neceiver output predicts left-handed antenna satellite within the vision, reduces direct projection passage satellite acquisition scope, exports No. PRN, left-handed antenna satellite within the vision, be designated as (PRN
1..., PRN
n);
Catch part to (PRN
1..., PRN
n) in satellite search for, judge that satellite is whether true visible, if truly visible, then guestimate carried out to direct signal code phase and Doppler frequency;
Tracking section carries out accurate code phase to the direct signal of acquisition success and Doppler estimates, is designated as (τ
di, f
di).
By selecting, star, compensation rate calculate reflection channel, carrier wave occurs, pseudo-code occurs and relevant treatment five part forms.
Select star to be select star to the satellite traced into according to elevation of satellite, select the satellite that elevation of satellite is maximum, its direct signal correspondence code phase place and Doppler's valuation are (τ
d, f
d), as the reference information of reflected signal process.
Compensation rate calculates the navigation information NAV_inf exported according to navigation neceiver and calculates compensation rate Δ τ and Δ f, is (τ by direct signal correspondence code phase place and Doppler's valuation
d, f
d) try to achieve the local pseudo-code phase τ of reflection channel
d+ Δ τ and local carrier frequency f
d+ Δ f.
Carrier wave occurs to produce local carrier and is used for the stripping of reflected signal carrier wave.
Pseudo-code occurs to produce local pseudo-code and is used for the stripping of reflected signal pseudo-code.
Local carrier and local pseudo-code and reflected signal RDIF are carried out related operation by relevant treatment, export delay-Doppler two-dimensional correlation power DDM (τ, f).
As Ocean Wind-field and sea level height inverting block diagram that Fig. 3 is this device.
Ocean surface wind retrieving calculates five parts by the generation of theoretical model storehouse, noise reduction process, normalization, Waveform Matching and residual sum of squares (RSS) and forms.This module is mainly estimated ocean surface wind speed, wind direction, can obtain theoretical two-dimensional correlation power in Wind-field Retrieval process by Waveform Matching, from then on can obtain the reflection delay time because sea surface roughness causes in theoretical two-dimensional correlation power.
It is from DDM (τ that theoretical model storehouse generates, f) extract the information such as the height of receiver, speed, elevation of satellite and position angle in the navigation information NAV_inf exported with navigation neceiver, information extraction be input to surface scattering signal correction power module:
In, input initial wind speed v
0, calculate the surface scattering signal correction power under different wind friction velocity, generative theory related power waveform library;
Noise reduction process carries out noise reduction process to actual measurement related power curve D DM (τ, f), and stress release treatment affects;
Normalization be to noise reduction process after actual measurement related power and the theoretical related power waveform library of generation be normalized respectively;
Waveform Matching be by normalization after actual measurement related power curve and theoretical related power curve adopt least square method to carry out Waveform Matching, export the wind speed v corresponding with surveying related power Curve Matching successful theoretical related power Curves;
It is by given initial wind direction φ that residual sum of squares (RSS) calculates
0calculate the theory of many Navsats and the residual sum of squares (RSS) of actual measurement related power waveform, adjustment wind direction, exports the wind direction φ that residual sum of squares (RSS) local minimum is corresponding.
Sea level height inverting is by noise reduction process, normalization, time-delay calculation, and time delay correction and high computational five part form.This module is mainly measured sea level height, by the time delay obtained in WIND FIELDS for revising the time delay in sea level height measurement, thus accurately measures sea level height.
Noise reduction process carries out noise reduction process to actual measurement related power curve D DM (τ, f), and stress release treatment affects, and completes during WIND FIELDS;
Normalization be to noise reduction process after actual measurement related power be normalized, complete during WIND FIELDS;
Time-delay calculation is by direct signal correlated power peak point code phase and the reflected signal correlated power peak point place code phase computational reflect signal time delay value τ relative to direct signal
e1;
Time delay correction is the time delay error τ caused by the theoretical related power curve estimation sea surface roughness that the match is successful
c, time delay value revises the delay inequality that senior middle school is surveyed on sea thus, then revised time delay value is: τ
e=τ
e1-τ
c.By time delay value τ
ethe path delay of trying to achieve direct signal and reflected signal is ρ
e=τ
e× c;
High computational is from navigation information NAV_inf, extract the information such as position and elevation of satellite of navigation neceiver, and trying to achieve receiver to the distance in the earth's core by navigation neceiver position is h, then sea level height H can by formula
try to achieve.
Claims (1)
1. based on a sea level height measurement mechanism for time delay correction, it is characterized in that: it is made up of a navigation neceiver, dextrorotation antenna, left-handed antenna, RF front-end module, an IF signal processing module and an inverting module, navigation neceiver is connected with inverting module with dextrorotation antenna, IF signal processing module respectively, and dextrorotation antenna and left-handed antenna independently of one another, and are connected with inverting module with RF front-end module, IF signal processing module together successively, the direct projection radio frequency voltage signal DRF that dextrorotation antenna exports and the reflected radio voltage signal RRF that left-handed antenna exports is as the input of RF front-end module, the direct projection digital medium-frequency signal DDRF that RF front-end module exports and reflection digital intermediate-freuqncy signal RDIF is as the input of IF signal processing module, the delay-Doppler two-dimensional correlation power DDM (τ that IF signal processing module exports, f) as the input of inverting module, the direct projection radio frequency voltage signal DRF that dextrorotation antenna exports is as the input of navigation neceiver, the navigation information NAV_inf that navigation neceiver exports is as the input of IF signal processing module and inverting module,
Described navigation neceiver is used for the navigation information NAV_inf of output packet containing the height of this receiver and speed, elevation of satellite, azimuth information, as the input of IF signal processing module and inverting module;
Described dextrorotation antenna adopts GNSS dextrorotation antenna for receiving the direct projection electromagnetic signal of Navsat, and this signal is changed into the input of direct projection radio frequency voltage signal DRF as radio-frequency front-end process;
Described left-handed antenna adopts the left-handed antenna of many array element array for receiving the reflected electromagnetic signal of Navsat, and this signal is changed into the input of reflected radio voltage signal RRF as radio-frequency front-end process;
The direct projection radio frequency voltage signal DRF received and reflected radio voltage signal RRF is down-converted to intermediate-freuqncy signal by described RF front-end module, and intermediate-freuqncy signal is converted into direct projection digital medium-frequency signal DDRF and reflection digital intermediate-freuqncy signal RDIF, as the input of IF signal processing module;
Described IF signal processing module is made up of direct projection passage and reflection channel, is connected with each other between the two; IF signal processing module is as the input of inverting module;
This direct projection passage by visible star anticipation, catch and follow the tracks of three parts and form, visible star anticipation, catch and be connected successively with following the tracks of; Visible star anticipation is predicted left-handed antenna satellite within the vision, reduces direct projection passage satellite acquisition scope, exports No. PRN, left-handed antenna satellite within the vision, be designated as (PRN
1..., PRN
n), as the input of catching; Catch (PRN
1..., PRN
n) in satellite search for, judge that satellite is whether true visible, if truly visible, then guestimate carried out to direct signal code phase and Doppler frequency, as the input of tracking; Follow the tracks of and accurate code phase and Doppler's estimation are carried out to the direct signal of acquisition success, be designated as (τ
di, f
di), as the input of reflection channel;
By selecting, star, compensation rate calculate this reflection channel, carrier wave occurs, pseudo-code occurs and relevant treatment five part forms, select star to calculate with compensation rate to be connected, compensation rate calculates and occurs to be connected with pseudo-code with carrier wave respectively, carrier wave occurs and pseudo-code occurs independent of one another, and carrier wave occurs to occur all to be connected with relevant treatment with pseudo-code; Select star to select star to the satellite traced into according to elevation of satellite, select the satellite that elevation of satellite is maximum, its direct signal correspondence code phase place and Doppler's valuation are designated as (τ
d, f
d), as the input that compensation rate calculates; Compensation rate calculates the navigation information NAV_inf exported according to navigation neceiver and calculates compensation rate Δ τ and Δ f, by direct signal correspondence code phase place and Doppler's valuation (τ
d, f
d) obtain the local pseudo-code phase τ of reflection channel
d+ Δ τ and local carrier frequency f
d+ Δ f, as the input that pseudo-code occurs and carrier wave occurs; Pseudo-code occurs to produce local pseudo-code and is used for the stripping of reflected signal pseudo-code; Carrier wave occurs to produce local carrier and is used for the stripping of reflected signal carrier wave; Local carrier and local pseudo-code and reflected signal RDIF are carried out related operation by relevant treatment, export delay-Doppler two-dimensional correlation power DDM (τ, f);
Described inverting module is made up of ocean surface wind retrieving module and sea level height inverting module, is connected with each other between the two;
This ocean surface wind retrieving module calculates five parts by the generation of theoretical model storehouse, noise reduction process, normalization, Waveform Matching and residual sum of squares (RSS) and forms, and the generation of theoretical model storehouse, noise reduction process, normalization, Waveform Matching and residual sum of squares (RSS) calculate and be connected successively; It is from DDM (τ that theoretical model storehouse generates, f) height of receiver, speed, elevation of satellite and azimuth information is extracted in the navigation information NAV_inf exported with navigation neceiver, information extraction is input in surface scattering signal correction power module, inputs initial wind speed v
0, calculate the surface scattering signal correction power under different wind friction velocity, generative theory related power waveform library; Noise reduction process carries out noise reduction process to actual measurement related power curve D DM (τ, f), and stress release treatment affects; Normalization be to noise reduction process after actual measurement related power and the theoretical related power waveform library of generation be normalized respectively; Waveform Matching be by normalization after actual measurement related power curve and theoretical related power curve adopt least square method to carry out Waveform Matching, export the wind speed v corresponding with surveying related power Curve Matching successful theoretical related power Curves; It is by given initial wind direction φ that residual sum of squares (RSS) calculates
0calculate the theory of many Navsats and the residual sum of squares (RSS) of actual measurement related power waveform, adjustment wind direction, exports the wind direction φ that residual sum of squares (RSS) local minimum is corresponding;
This sea level height inverting module is by noise reduction process, normalization, time-delay calculation, and time delay correction and high computational five part form, and noise reduction process, normalization, time-delay calculation, time delay correction are connected successively with high computational; Noise reduction process carries out noise reduction process to actual measurement related power curve D DM (τ, f), and stress release treatment affects; Normalization be to noise reduction process after actual measurement related power be normalized; Time-delay calculation is by direct signal correlated power peak point code phase and the reflected signal correlated power peak point place code phase computational reflect signal time delay value τ relative to direct signal
e1; Time delay correction is the time delay error τ caused by the theoretical related power curve estimation sea surface roughness that the match is successful
c, time delay value revises the delay inequality that senior middle school is surveyed on sea thus, by revised time delay value τ
eobtain the path delay of direct signal and reflected signal; High computational is from navigation information NAV_inf, extract position and the satellite altitude angle information of navigation neceiver, tries to achieve the distance of navigation neceiver to the earth's core by navigation neceiver position, thus calculates sea level height H.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510071278.2A CN104678371B (en) | 2015-02-11 | 2015-02-11 | A kind of sea level height measurement apparatus based on time delay amendment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510071278.2A CN104678371B (en) | 2015-02-11 | 2015-02-11 | A kind of sea level height measurement apparatus based on time delay amendment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104678371A true CN104678371A (en) | 2015-06-03 |
| CN104678371B CN104678371B (en) | 2018-01-09 |
Family
ID=53313699
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510071278.2A Active CN104678371B (en) | 2015-02-11 | 2015-02-11 | A kind of sea level height measurement apparatus based on time delay amendment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104678371B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105301622A (en) * | 2015-11-02 | 2016-02-03 | 北京航大泰科信息技术有限公司 | Wind speed detection device based on navigational satellite reflected signal |
| CN107017707A (en) * | 2015-09-15 | 2017-08-04 | 艾诺格思公司 | Recognize the receiver in wireless charging transmission field |
| CN110824512A (en) * | 2019-11-26 | 2020-02-21 | 中国科学院国家空间科学中心 | Non-uniform chip real-time delay Doppler mapping data generator |
| CN110927221A (en) * | 2018-10-08 | 2020-03-27 | 中国石油化工股份有限公司 | Sea surface oil spill detection device |
| CN112880633A (en) * | 2021-01-12 | 2021-06-01 | 上海海洋大学 | Sea surface height measuring method based on Berger algorithm |
| CN113031015A (en) * | 2021-03-05 | 2021-06-25 | 北京航空航天大学 | Sea ice detection device and method based on GNSS-R carrier phase |
| CN113049062A (en) * | 2021-03-12 | 2021-06-29 | 北京航空航天大学 | Device for measuring lake water level through GNSS direct reflection signal carrier interference |
| CN113671210A (en) * | 2021-08-19 | 2021-11-19 | 北京航空航天大学 | Device for measuring river flow velocity by land-based GNSS direct reflection signal carrier interference |
| CN116539913A (en) * | 2023-05-04 | 2023-08-04 | 极诺星空(北京)科技有限公司 | Method and device for on-board real-time inversion of sea surface wind speed |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101833090A (en) * | 2010-03-12 | 2010-09-15 | 中国科学院遥感应用研究所 | Airborne ocean microwave remote sensing system utilizing signal sources of global satellite positioning system |
| CN101975959A (en) * | 2010-09-03 | 2011-02-16 | 北京航空航天大学 | Device for acquiring inversion significant wave height and relative elevation in delay mapping receiver |
| CN102749067A (en) * | 2012-06-21 | 2012-10-24 | 北京航空航天大学 | Telemetering device for accurately inverting sea level elevation with self-adaptive multistage calculation |
-
2015
- 2015-02-11 CN CN201510071278.2A patent/CN104678371B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101833090A (en) * | 2010-03-12 | 2010-09-15 | 中国科学院遥感应用研究所 | Airborne ocean microwave remote sensing system utilizing signal sources of global satellite positioning system |
| CN101975959A (en) * | 2010-09-03 | 2011-02-16 | 北京航空航天大学 | Device for acquiring inversion significant wave height and relative elevation in delay mapping receiver |
| CN102749067A (en) * | 2012-06-21 | 2012-10-24 | 北京航空航天大学 | Telemetering device for accurately inverting sea level elevation with self-adaptive multistage calculation |
Non-Patent Citations (3)
| Title |
|---|
| 周兆明等: ""渤海GNSS-R机载测风试验及其反演结果"", 《武汉大学学报》 * |
| 李伟强等: ""面向遥感的GNSS反射信号接收处理系统及实验"", 《武汉大学学报》 * |
| 路勇等: ""基于GNSS_R的海面风场监测系统研究"", 《武汉大学学报》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107017707A (en) * | 2015-09-15 | 2017-08-04 | 艾诺格思公司 | Recognize the receiver in wireless charging transmission field |
| CN105301622A (en) * | 2015-11-02 | 2016-02-03 | 北京航大泰科信息技术有限公司 | Wind speed detection device based on navigational satellite reflected signal |
| CN110927221A (en) * | 2018-10-08 | 2020-03-27 | 中国石油化工股份有限公司 | Sea surface oil spill detection device |
| CN110824512A (en) * | 2019-11-26 | 2020-02-21 | 中国科学院国家空间科学中心 | Non-uniform chip real-time delay Doppler mapping data generator |
| CN110824512B (en) * | 2019-11-26 | 2022-01-25 | 中国科学院国家空间科学中心 | Non-uniform chip real-time delay Doppler mapping data generator |
| CN112880633A (en) * | 2021-01-12 | 2021-06-01 | 上海海洋大学 | Sea surface height measuring method based on Berger algorithm |
| CN113031015A (en) * | 2021-03-05 | 2021-06-25 | 北京航空航天大学 | Sea ice detection device and method based on GNSS-R carrier phase |
| CN113049062A (en) * | 2021-03-12 | 2021-06-29 | 北京航空航天大学 | Device for measuring lake water level through GNSS direct reflection signal carrier interference |
| CN113049062B (en) * | 2021-03-12 | 2022-04-15 | 北京航空航天大学 | Device for measuring lake water level through GNSS direct reflection signal carrier interference |
| CN113671210A (en) * | 2021-08-19 | 2021-11-19 | 北京航空航天大学 | Device for measuring river flow velocity by land-based GNSS direct reflection signal carrier interference |
| CN116539913A (en) * | 2023-05-04 | 2023-08-04 | 极诺星空(北京)科技有限公司 | Method and device for on-board real-time inversion of sea surface wind speed |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104678371B (en) | 2018-01-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104678371A (en) | Device for measuring sea surface height based on time-delay modification | |
| CN106646538B (en) | A kind of deformation monitoring GNSS signal multipath correcting method based on single poor filtering | |
| CN103033828B (en) | High-sensitivity compass-assisted time servicing device, time service receiver and time service method | |
| Pirazzi et al. | Preliminary performance analysis with a GPS+ Galileo enabled chipset embedded in a smartphone | |
| CN103842846B (en) | A kind of for low-power consumption with the system of low cost GNSS receiver, method and computer program | |
| CN104677463B (en) | A kind of GNSS R multi-sources sea height measuring device based on statistic | |
| US7030814B2 (en) | System and method to estimate the location of a receiver in a multi-path environment | |
| CN109709591A (en) | A kind of GNSS high-precision locating method towards intelligent terminal | |
| CN102033236B (en) | Position and speed combined estimation method for satellite navigation | |
| CN102749067A (en) | Telemetering device for accurately inverting sea level elevation with self-adaptive multistage calculation | |
| CN103605140A (en) | Occultation signal rapid recovery method | |
| CN109444935A (en) | A kind of Doppler's detection and reparation for cycle slips method of low sampling rate | |
| CN105425258A (en) | Highly-dynamical weak signal GPS capturing method assisted by inertial navigation system | |
| CN105182380A (en) | Hardware receiver and method for achieving GNSS-R phase difference extraction | |
| US20140266910A1 (en) | Techniques to improve the performance of a fixed, timing-based radio positioning network using external assistance information | |
| Bahrami et al. | Instantaneous Doppler-aided RTK positioning with single frequency receivers | |
| CN105044744A (en) | Non-cooperative navigation satellite identification method based on Doppler frequency shift characteristic | |
| CN110346816B (en) | Real-time differential positioning method of ionosphere modeling constraint under short-to-long baseline scene | |
| US10598757B2 (en) | Systems and methods for improving the performance of a timing-based radio positioning network using estimated range biases | |
| Narula et al. | Accelerated collective detection technique for weak GNSS signal environment | |
| CN111007490B (en) | Sky wave over-the-horizon radar coordinate registration method based on buoy geographic information | |
| CN103940424A (en) | Detection and estimation method for integer ambiguity of XNAV (X-ray pulsar-based navigation) signals | |
| Li et al. | Implementation and performance of clock steering in a software GPS L1 single frequency receiver | |
| Xiang et al. | A new inter-satellite ranging method based on pseudo-range and dual-frequency carrier phase | |
| EP3748401A1 (en) | System and method for position determination of a stationary gnss receiver using a distributed time signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211029 Address after: 272001 fifth floor, building A3, industry university research base, high tech Zone, Jining City, Shandong Province Patentee after: SHANDONG HANGXIANG ELECTRONIC SCIENCE & TECHNOLOGY Co.,Ltd. Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road Patentee before: BEIHANG University |
|
| TR01 | Transfer of patent right |