CN111650618A - GNSS occultation detection signal processing method - Google Patents
GNSS occultation detection signal processing method Download PDFInfo
- Publication number
- CN111650618A CN111650618A CN202010615777.4A CN202010615777A CN111650618A CN 111650618 A CN111650618 A CN 111650618A CN 202010615777 A CN202010615777 A CN 202010615777A CN 111650618 A CN111650618 A CN 111650618A
- Authority
- CN
- China
- Prior art keywords
- data
- signals
- signal processing
- processing method
- baseband
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a GNSS occultation detection signal processing method, which comprises the following steps: s1: receiving signals of a plurality of antennas; s2: respectively carrying out radio frequency processing on the antenna signals to convert the antenna signals into digital signals; s3: respectively carrying out baseband processing on the digital signals to form baseband data; s4: preprocessing the baseband data and eliminating invalid data; s5: respectively encoding the data to form measurement data; s6: the measurement data is subjected to time sequence processing, and a plurality of measurement data in the same time sequence are combined into one data. The invention separately carries out the processes of radio frequency processing, baseband processing and coding of signals, reduces the mutual coupling of the signals and improves the reliability of signal processing; the preprocessing eliminates invalid data before coding, reduces the data processing burden and improves the processing efficiency; the time sequence processing completes the integration of multiple paths of different signals, and the data are judged and processed in multiple ways, so that the continuity and the measurement precision of the measurement data are improved.
Description
Technical Field
The invention belongs to the technical field of space detection, and particularly relates to a GNSS occultation detection signal processing method.
Background
In recent years, GNSS is rapidly developed in various countries around the world, wherein the aerospace industry has stepped into the commercialization process in various fields, and the application of GNSS remote sensing technology has made a breakthrough progress. The GNSS remote sensing technology is used to invert meteorological parameters such as the earth atmosphere, the ocean, and the soil by using changes in physical quantities such as signal amplitude, phase, etc., of radio waves propagated in the earth atmosphere or reflected by ground objects. The inversion products of GNSS occultation detection are ionospheric and atmospheric parameters. The ionosphere inversion product comprises an ionosphere electron density profile and an ionosphere scintillation index, and the atmosphere inversion product comprises an atmosphere refractive index, density, pressure, temperature and humidity profile, and has the outstanding characteristics of low cost, rich and rapid data acquisition quantity, high measurement precision and the like, and has extremely wide application prospect. GNSS occultation detection is currently implemented using a GNSS occultation receiver mounted on earth's low orbit satellites, a positioning antenna, a pair of forward occultation antennas and a pair of backward occultation antennas. In the universe environment, due to the influence of on-orbit single event upset and other reasons, data can be abnormal, and the accuracy of detection results is influenced.
Disclosure of Invention
Aiming at the technical problems in the prior art, the invention provides the GNSS occultation detection signal processing method, which optimizes the processing process of each antenna signal, performs multiple judgment and processing on the measurement data, and improves the continuity and the measurement precision of the measurement data.
The technical scheme adopted by the invention is as follows: a GNSS occultation detection signal processing method comprises the following steps:
s1: receiving signals of a plurality of antennas;
s2: respectively carrying out radio frequency processing on the antenna signals to convert the antenna signals into digital signals;
s3: respectively carrying out baseband processing on the digital signals to form baseband data;
s4: preprocessing the baseband data and eliminating invalid data;
s5: respectively encoding the data to form measurement data;
s6: the measurement data is subjected to time sequence processing, and a plurality of measurement data in the same time sequence are combined into one data.
The antenna signal in step S1 includes one positioning signal and two occultation signals.
The radio frequency processing and the baseband processing of the two masker signals in steps S2 and S3 are the same.
The preprocessing in step S4 is to determine the quality of the baseband data and reject data that will not meet the requirements.
In step S5, the encoding modes of the baseband data of the two occultants are the same, and all the measurement data transmit the measurement data at the same time sequence to the next step in a ping-pong buffer mode.
In step S6, the frame headers of the data at the same time sequence are respectively determined, and the data that do not meet the conditions are deleted;
the frame count value of the data is then determined and compared: if the frame count values are the same, all the data are synthesized into a new data; if the frame count values are not completely the same, deleting the data of the single frame count value, performing interpolation supplement on the missing data, and synthesizing all the data into a new data.
The number of data in the same time sequence is three; if the frame count values of the three data are the same, the data are not deleted; deleting the data with different frame count values if the frame count values of the two data are the same; and if the frame count values of the three data are different, deleting the three data.
When new data is synthesized, all frames of the data are filled to the same length and are arranged in sequence to be used as the data information part of the frames of the new data.
The other data is filled with 16-ary 5A5A at the end of the data with the longest data as a reference.
And the missing data is subjected to interpolation supplementation according to the corresponding data content of the previous period and a satellite orbit dynamics model.
Compared with the prior art, the invention has the beneficial effects that:
1. the invention separately carries out the radio frequency processing, the baseband processing and the coding process of the signals, reduces the mutual coupling of the signals and improves the reliability of the signal processing;
2. according to the invention, by preprocessing, invalid data is removed before encoding, so that the data processing burden is reduced, and the processing efficiency is improved;
3. the invention completes the integration of a plurality of paths of different signals through time sequence processing, thereby improving the transmission efficiency and the measurement precision;
4. in the time sequence processing process, the reliability of data is increased through frame header judgment and frame counting judgment, data abnormity caused by on-orbit single particle overturn and the like is avoided, meanwhile, the judgment of data quality is realized through comparing frame counting values, a proper data packaging mode is selected according to the judgment result, and the measurement precision is ensured as far as possible under the condition of data loss through an interpolation method adopting a track dynamics model.
Drawings
FIG. 1 is a flow chart of an embodiment of the present invention;
FIG. 2 is a flow chart of a timing process according to an embodiment of the invention.
Detailed Description
In order to make the technical solutions of the present invention better understood, the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
The embodiment of the invention discloses a GNSS occultation detection signal processing method, which comprises the following steps as shown in figures 1-2:
s1: receiving signals of three antennas, wherein the antenna signals comprise a positioning signal and two occultation signals; the signal frequency band is GPS L1/L2/L5 or Beidou B1/B2/B3;
s2: respectively carrying out radio frequency processing on the antenna signals to obtain digital signals, wherein the radio frequency processing modes of the two occultation signals are the same;
s3: respectively carrying out baseband processing on the digital signals to form baseband data, wherein the baseband processing modes of the two occultation signals are the same;
s4: preprocessing the baseband data, judging the quality of the baseband data, and eliminating data which do not meet the requirements;
s5: respectively encoding the data to form measurement data; the coding modes of the baseband data of the two occultants are the same, and all the measurement data respectively adopt a ping-pong cache mode to transmit the measurement data in the same time sequence to the next step;
s6: performing time sequence processing on the measurement data, and combining three measurement data in the same time sequence into one data: respectively judging frame headers of the three data at the same time sequence, and deleting the data which do not meet the conditions;
the frame count values of the three data are then determined and compared: if the frame count values of the three data are the same, the data are not deleted; deleting the data with different frame count values if the frame count values of the two data are the same; if the frame count values of the three data are different, deleting the three data;
if the deleted data exists in the three measured data in the same time sequence, carrying out interpolation supplement on the missing data according to the corresponding data content in the previous period and a satellite orbit dynamics model;
after confirming that all three data exist, filling 16-system 5A5A at the end of the other two data by taking the longest data as a reference for length filling; arranging the three data filled to the same length in sequence to serve as a data information part of a new data frame, adding a new total frame header, and calculating the checksum of the combined data to form a new data frame;
s7: and storing new data and waiting for the instruction to be transmitted.
The present invention has been described in detail with reference to the embodiments, but the description is only illustrative of the present invention and should not be construed as limiting the scope of the present invention. The scope of the invention is defined by the claims. The technical solutions of the present invention or those skilled in the art, based on the teaching of the technical solutions of the present invention, should be considered to be within the scope of the present invention, and all equivalent changes and modifications made within the scope of the present invention or equivalent technical solutions designed to achieve the above technical effects are also within the scope of the present invention.
Claims (10)
1. A GNSS occultation detection signal processing method is characterized in that: the method comprises the following steps:
s1: receiving signals of a plurality of antennas;
s2: respectively carrying out radio frequency processing on the antenna signals to convert the antenna signals into digital signals;
s3: respectively carrying out baseband processing on the digital signals to form baseband data;
s4: preprocessing the baseband data and eliminating invalid data;
s5: respectively encoding the data to form measurement data;
s6: the measurement data is subjected to time sequence processing, and a plurality of measurement data in the same time sequence are combined into one data.
2. The GNSS masker detect signal processing method of claim 1, wherein: the antenna signal in step S1 includes one positioning signal and two occultation signals.
3. The GNSS masker detect signal processing method of claim 2, wherein: the radio frequency processing and the baseband processing of the two masker signals in steps S2 and S3 are the same.
4. The GNSS masker detect signal processing method of claim 1, wherein: the preprocessing in step S4 is to determine the quality of the baseband data and reject data that will not meet the requirements.
5. The GNSS masker detect signal processing method of claim 2, wherein: in step S5, the encoding modes of the baseband data of the two occultants are the same, and all the measurement data transmit the measurement data at the same time sequence to the next step in a ping-pong buffer mode.
6. The GNSS masker detect signal processing method of claim 1, wherein: in step S6, the frame headers of the data at the same time sequence are respectively determined, and the data that do not meet the conditions are deleted;
the frame count value of the data is then determined and compared: if the frame count values are the same, all the data are synthesized into a new data; if the frame count values are not completely the same, deleting the data of the single frame count value, performing interpolation supplement on the missing data, and synthesizing all the data into a new data.
7. The GNSS masker detect signal processing method of claim 6, wherein: the number of data in the same time sequence is three; if the frame count values of the three data are the same, the data are not deleted; deleting the data with different frame count values if the frame count values of the two data are the same; and if the frame count values of the three data are different, deleting the three data.
8. The GNSS masker detection signal processing method of claim 6 or 7, wherein: when new data is synthesized, all frames of the data are filled to the same length and are arranged in sequence to be used as the data information part of the frames of the new data.
9. The GNSS masker detect signal processing method of claim 8, wherein: the other data is filled with 16-ary 5A5A at the end of the data with the longest data as a reference.
10. The GNSS masker detection signal processing method of claim 6 or 7, wherein: and the missing data is subjected to interpolation supplementation according to the corresponding data content of the previous period and a satellite orbit dynamics model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010615777.4A CN111650618B (en) | 2020-06-30 | 2020-06-30 | GNSS occultation detection signal processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010615777.4A CN111650618B (en) | 2020-06-30 | 2020-06-30 | GNSS occultation detection signal processing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111650618A true CN111650618A (en) | 2020-09-11 |
| CN111650618B CN111650618B (en) | 2023-06-20 |
Family
ID=72349852
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010615777.4A Active CN111650618B (en) | 2020-06-30 | 2020-06-30 | GNSS occultation detection signal processing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111650618B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115826007A (en) * | 2023-02-17 | 2023-03-21 | 天津云遥宇航科技有限公司 | Method for forecasting neutral atmosphere occultation event of satellite-based global navigation satellite system |
| CN116338735A (en) * | 2023-05-11 | 2023-06-27 | 天津云遥宇航科技有限公司 | Ionosphere occultation flicker index S4 calculation method based on Butterworth filtering |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101217628A (en) * | 2008-01-08 | 2008-07-09 | 北京航空航天大学 | A diversity combination device and the corresponding method of multi-aerial ground digital television receiver |
| CN102288974A (en) * | 2011-05-10 | 2011-12-21 | 航天恒星科技有限公司 | System for generating satellite navigation occultation signal |
| CN102510493A (en) * | 2011-10-09 | 2012-06-20 | 航天恒星科技有限公司 | Remote sensing satellite image fault-tolerance display method |
| CN105141352A (en) * | 2015-07-24 | 2015-12-09 | 哈尔滨工业大学 | Satellite high-speed data transmission baseband data error statistics and frame sequencing processing system and method |
| CN105607102A (en) * | 2016-01-06 | 2016-05-25 | 中国人民解放军海军航空工程学院青岛校区 | Multi-source-positioning-data-fusion-processing-based flight reappearance system and flight reappearance method thereof |
| CN105891856A (en) * | 2016-04-07 | 2016-08-24 | 中国科学院国家空间科学中心 | GNSS radio occultation detector |
| US20160327686A1 (en) * | 2015-05-08 | 2016-11-10 | European Space Agency | Interferometric radio occultation |
| CN107679152A (en) * | 2017-09-26 | 2018-02-09 | 中国科学院电子学研究所 | Data processing method based on multi-layer information joint index |
| CN107885772A (en) * | 2017-09-29 | 2018-04-06 | 北京空间机电研究所 | A kind of dynamically changeable multiple data format satellite assistance data processing system |
| CN109061682A (en) * | 2018-09-30 | 2018-12-21 | 中国气象局气象探测中心 | LEO occultation system and method suitable for aerostatics |
| CN111092650A (en) * | 2019-12-26 | 2020-05-01 | 北京遥测技术研究所 | X frequency channel measurement and control communication integrated system based on data system fusion |
-
2020
- 2020-06-30 CN CN202010615777.4A patent/CN111650618B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101217628A (en) * | 2008-01-08 | 2008-07-09 | 北京航空航天大学 | A diversity combination device and the corresponding method of multi-aerial ground digital television receiver |
| CN102288974A (en) * | 2011-05-10 | 2011-12-21 | 航天恒星科技有限公司 | System for generating satellite navigation occultation signal |
| CN102510493A (en) * | 2011-10-09 | 2012-06-20 | 航天恒星科技有限公司 | Remote sensing satellite image fault-tolerance display method |
| US20160327686A1 (en) * | 2015-05-08 | 2016-11-10 | European Space Agency | Interferometric radio occultation |
| CN105141352A (en) * | 2015-07-24 | 2015-12-09 | 哈尔滨工业大学 | Satellite high-speed data transmission baseband data error statistics and frame sequencing processing system and method |
| CN105607102A (en) * | 2016-01-06 | 2016-05-25 | 中国人民解放军海军航空工程学院青岛校区 | Multi-source-positioning-data-fusion-processing-based flight reappearance system and flight reappearance method thereof |
| CN105891856A (en) * | 2016-04-07 | 2016-08-24 | 中国科学院国家空间科学中心 | GNSS radio occultation detector |
| CN107679152A (en) * | 2017-09-26 | 2018-02-09 | 中国科学院电子学研究所 | Data processing method based on multi-layer information joint index |
| CN107885772A (en) * | 2017-09-29 | 2018-04-06 | 北京空间机电研究所 | A kind of dynamically changeable multiple data format satellite assistance data processing system |
| CN109061682A (en) * | 2018-09-30 | 2018-12-21 | 中国气象局气象探测中心 | LEO occultation system and method suitable for aerostatics |
| CN111092650A (en) * | 2019-12-26 | 2020-05-01 | 北京遥测技术研究所 | X frequency channel measurement and control communication integrated system based on data system fusion |
Non-Patent Citations (6)
| Title |
|---|
| 于长海等: "多站遥测数据处理中的数据融合方法研究", 《遥测遥控》 * |
| 于长海等: "多站遥测数据处理中的数据融合方法研究", 《遥测遥控》, no. 01, 15 January 2018 (2018-01-15) * |
| 刘夫体;唐海波;张香成;: "数据修补技术在GPS测量数据处理中的应用", 遥测遥控, no. 05 * |
| 吴悠等: "多站点组网的实时遥测数据选定技术", 《测控技术》 * |
| 吴悠等: "多站点组网的实时遥测数据选定技术", 《测控技术》, no. 06, 18 June 2016 (2016-06-18) * |
| 徐洪洲;梁红;韩成柱;: "多站位测量体制下遥测预处理数据提取研究", 信息技术, no. 02 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115826007A (en) * | 2023-02-17 | 2023-03-21 | 天津云遥宇航科技有限公司 | Method for forecasting neutral atmosphere occultation event of satellite-based global navigation satellite system |
| CN116338735A (en) * | 2023-05-11 | 2023-06-27 | 天津云遥宇航科技有限公司 | Ionosphere occultation flicker index S4 calculation method based on Butterworth filtering |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111650618B (en) | 2023-06-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107102346B (en) | Multi-antenna attitude measurement method based on Beidou system | |
| CN108072887B (en) | Single-base-station long-distance offshore real-time dynamic positioning method | |
| CN108828640B (en) | Method and device for weighting satellite navigation positioning observation values | |
| CN111650618A (en) | GNSS occultation detection signal processing method | |
| CN112014860A (en) | Low-orbit satellite space-time reference establishment method based on Beidou PPP-RTK | |
| CN105549056A (en) | Relative positioning device and carrier wave integer ambiguity calculation method thereof | |
| CN103728643B (en) | With the Big Dipper three network RTK blur level single epoch fixing means frequently that wide lane retrains | |
| CN110146906A (en) | Remote time transmission method based on single poor carrier phase observation data | |
| CN113391334B (en) | Space signal ranging error correction coding method based on Beidou global short message | |
| CN107765275A (en) | GPS wide area differential GPS localization method, device, terminal and computer-readable recording medium | |
| WO2022156480A1 (en) | Clock error predicting method and device | |
| KR101419339B1 (en) | An apparatus for processing differential information of differential global navigation satellite system and the method thereof | |
| CN113267790A (en) | Satellite orbit clock error correction number coding method based on Beidou global short message communication | |
| CN113671537B (en) | Three-frequency beacon signal ionosphere channel simulation method | |
| CN107966722A (en) | A kind of GNSS satellite clock solutions method | |
| CN105045086A (en) | Improved high precision RDSS time service method | |
| CN109752737B (en) | Preprocessing method for inter-satellite Ka-band bidirectional measurement pseudo range of navigation satellite | |
| CN103592653A (en) | Ionized layer delay correction method for local area single-frequency satellite navigation user | |
| CN111596328B (en) | Time sequence control method for data of GNSS occultation receiver | |
| CN115980317A (en) | Foundation GNSS-R data soil moisture estimation method based on corrected phase | |
| CN116068670A (en) | Real-time reconstruction method and device for Beidou water vapor field suitable for ground-free network area | |
| CN114061622B (en) | Error calibration method for deep space three-dimensional ranging system | |
| KR101344426B1 (en) | An apparatus for processing differential information of differential global navigation satellite system and the method thereof | |
| CN114280650A (en) | Method and system for acquiring inter-frequency deviation of GNSS satellite based on LEO | |
| CN115561793A (en) | Real-time Beidou phase decimal deviation rapid estimation method based on parallel computation |
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 |