CN112087255A - Ground inversion method for GNSS occultation data - Google Patents
Ground inversion method for GNSS occultation data Download PDFInfo
- Publication number
- CN112087255A CN112087255A CN202010933835.8A CN202010933835A CN112087255A CN 112087255 A CN112087255 A CN 112087255A CN 202010933835 A CN202010933835 A CN 202010933835A CN 112087255 A CN112087255 A CN 112087255A
- Authority
- CN
- China
- Prior art keywords
- data
- frame
- occultation
- frame header
- file
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18513—Transmission in a satellite or space-based system
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L69/00—Network arrangements, protocols or services independent of the application payload and not provided for in the other groups of this subclass
- H04L69/22—Parsing or analysis of headers
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- Detection And Prevention Of Errors In Transmission (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
The invention provides a ground inversion method of GNSS occultation data, which comprises the following steps: s1, receiving the data transmitted by the ground station, and decrypting and restoring the occultation data; s2, performing sticky packaging processing on the occultation data according to the frame header, and storing each frame of data in a classified manner; s3, packaging the data, and storing the data belonging to the same occultation event into the same RINEX format file; s4, repairing missing values and abnormal values in the masker data file in the RINEX format; s5 inversion calculation is carried out through the ephemeris file and the repaired occultation data file in the RINEX format; s6 stores the process data and the result data of the inversion calculation. The invention realizes that data abnormity does not occur in the process of processing broken packets and sticky packets, improves the inversion success rate in a plurality of modes, and meets the real-time processing requirement of mass occultation data caused by multi-channel, multi-system and multi-type occultation events.
Description
Technical Field
The invention belongs to the technical field of GNSS signal processing, and particularly relates to a ground inversion method of GNSS occultation data.
Background
The current GNSS systems mainly include beidou in china, GPS in the united states, GLONASS in russia, and GALILEO in the european union, which are four major global satellite navigation systems. There are hundreds of satellites in orbit in four satellite navigation systems. The GNSS occultation means that electric wave signals transmitted by the GNSS satellite are covered by the earth atmosphere and reach an observation satellite after being refracted by the earth atmosphere and an ionized layer. The refracted GNSS signal received by the observation satellite can be used for inverting the vertical profile of atmospheric temperature, density, air pressure and ionosphere electron density, and the obtained data can be subjected to anabolic treatment in the later period so as to perform accurate weather forecasting service. Meanwhile, meteorological data must be processed in time to ensure good timeliness. Otherwise, the masker data would lose value.
At present, the ground station sends the occultation data sent by the satellite to the inversion system through TCP/IP at certain time intervals, generally 30 minutes. Because massive occultation data caused by multi-channel, multi-system and multi-type occultation events need to be processed in real time, in order to meet the transmission requirement of the occultation data on large data volume, long connection is generally adopted for data transmission, and therefore the problems of packet breaking and packet sticking are inevitably caused. The existing inversion method cannot well process broken packets and sticky packets, data abnormity easily occurs, and the success rate of inversion is influenced.
Disclosure of Invention
The invention provides a ground inversion method of GNSS occultation data aiming at the technical problems in the prior art, which realizes that data abnormity does not occur in the process of processing broken packets and sticky packets, improves the inversion success rate in multiple modes and meets the real-time processing requirement of mass occultation data caused by multi-channel, multi-system and multi-type occultation events.
The technical scheme adopted by the invention is as follows: a ground inversion method of GNSS occultation data comprises the following steps:
s1, receiving the data transmitted by the ground station, and decrypting and restoring the occultation data;
s2, performing sticky packaging processing on the occultation data according to the frame header, and storing each frame of data in a classified manner;
s3, packaging the data, and storing the data belonging to the same occultation event into the same RINEX format file;
s4, repairing missing values and abnormal values in the masker data file in the RINEX format;
s5 inversion calculation is carried out through the ephemeris file and the repaired occultation data file in the RINEX format;
s6 stores the process data and the result data of the inversion calculation.
Preferably, in step S2, the masker data is divided into reference star data, ionosphere data, and neutral atmospheric layer data, each type of data corresponds to a frame header with the same length and different contents, during the sticky packet processing, the frame header is identified, the frame header type and the frame data start position are recorded, and then the frame data are classified and stored according to the frame header type.
Preferably, the frame format of the occultation data is respectively composed of a frame header, a frame count, a data body and a check, wherein the frame header and the frame count both occupy N bytes, and the category of the occultation data is distinguished according to the frame header.
Preferably, during the sticky packet processing, it is first determined whether there are two packets of data spanned by the frame header: searching a occultation data frame header from the last N-1 bytes of the upper packet data and the initial N-1 bytes of the present packet data, and recording the type of the frame header and the initial position of the frame data if the frame header exists; judging whether the data of the packet has a frame header, if so, recording the type of the frame header and the initial position of frame data; recording the last N-1 bytes of the data of the packet; and storing the frame into the corresponding frame buffer area according to the record, and recording the frame buffer area into which the last frame data is stored.
Preferably, the process data and the result data in the step S6 are saved as NC format files.
Compared with the prior art, the invention has the beneficial effects that:
1. the invention disassembles the continuous occultation data into three kinds of data, is convenient for the subsequent data processing, defines the data format for the three kinds of occultation data, adopts the frame headers with the same length and different contents, and distinguishes the occultation data category through the frame headers, thereby effectively solving the problems of packet sticking and packet breaking;
2. the invention saves the last byte data of each packet of data so as to compare with the next packet of data, and searches whether a frame header spans two packets of data, thereby preventing missing data;
3. when the invention is packaged, the data belonging to the same occultation event is stored into the same RINEX format file, which is convenient for the subsequent data processing;
4. the method repairs the missing value and the abnormal value in the occultation data file, and improves the success rate of the later inversion;
5. the method saves the process data and the result data of the inversion calculation, improves the utilization rate of the data, can well apply the process data in many fields, and is convenient for customers to use as the data are saved as standard NC format files.
Drawings
FIG. 1 is a schematic structural diagram of an embodiment of the present invention;
fig. 2 is a schematic view of a package sticking processing flow according to an embodiment of the present 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.
An embodiment of the present invention provides a ground inversion method for GNSS occultation data, as shown in fig. 1-2, which includes the following steps:
s1, receiving the data transmitted by the ground station, and decrypting and restoring the occultation data;
the S2 occultation data comprises reference star data, ionized layer data and neutral atmosphere data, and the content, frequency and format of each data are different. And defining a frame format of the occultation data, wherein the frame format comprises a frame header, a frame count, a data body and a check, and the frame header and the frame count respectively occupy 4 bytes.
When the package is stuck, firstly judging whether frame header spans two package data: searching a occultation data frame header from the last 3 bytes of the upper packet data and the initial 3 bytes of the present packet data, and recording the frame header type and the frame data initial position if the frame header exists; judging whether the data of the packet has a frame header, if so, recording the type of the frame header and the initial position of frame data; finally, the last 3 bytes of the packet data are recorded.
Storing the frames into corresponding frame buffers according to the records: firstly, storing data before a first frame header recorded in the packet data into a last frame data buffer area of a last packet data memory, and indicating that the data at the beginning of the packet belongs to last frame data of the last packet; then, copying data from the data start position corresponding to the recorded frame header to the corresponding frame buffer, if no frame header is found in the packet of data, it is indicated that the packet of data is the last frame of data stored in the previous packet, and placing the packet of data in the corresponding buffer. And finally, recording the frame buffer area into which the last frame data is stored, and judging the frame buffer area into which the foremost data of the next packet of data is stored.
S3, respectively packaging the data in the frame buffer area, storing the binary data into a RINEX format file, and storing the data belonging to the same occultation event into the same RINEX format file; generating a reference star data file, an ionosphere data file and a neutral atmospheric layer data file;
s4, repairing missing values and abnormal values in the masker data file in the RINEX format, and improving the success rate of inversion;
s5 inversion calculation is carried out through the ephemeris file and the repaired occultation data file in the RINEX format;
s6 stores the process data and the result data of the inversion calculation. The process data and the result data are saved as NC format files. The process data includes a coordinate/clock difference file, a neutral atmosphere additive phase file, and an ionosphere additive phase file, and the result data includes an ionosphere profile file and a neutral atmosphere profile file.
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 (5)
1. A ground inversion method of GNSS occultation data is characterized in that: the method comprises the following steps:
s1, receiving the data transmitted by the ground station, and decrypting and restoring the occultation data;
s2, performing sticky packaging processing on the occultation data according to the frame header, and storing each frame of data in a classified manner;
s3, packaging the data, and storing the data belonging to the same occultation event into the same RINEX format file;
s4, repairing missing values and abnormal values in the masker data file in the RINEX format;
s5 inversion calculation is carried out through the ephemeris file and the repaired occultation data file in the RINEX format;
s6 stores the process data and the result data of the inversion calculation.
2. The GNSS occultation data ground inversion method of claim 1, wherein: in step S2, the occultation data is divided into reference star data, ionization layer data, and neutral atmospheric layer data, each type of data corresponds to a frame header with the same length and different contents, during the sticky-wrapping process, the frame header is identified, the frame header type and the frame data start position are recorded, and then the frame data are stored in a classified manner according to the frame header type.
3. The GNSS occultation data ground inversion method of claim 2, wherein: the frame format of the occultation data is respectively composed of a frame head, a frame count, a data body and a check, wherein the frame head and the frame count both occupy N bytes, and the category of the occultation data is distinguished according to the frame head.
4. The GNSS occultation data ground inversion method of claim 3, wherein: when the package is stuck, firstly judging whether frame header spans two package data: searching a occultation data frame header from the last N-1 bytes of the upper packet data and the initial N-1 bytes of the present packet data, and recording the type of the frame header and the initial position of the frame data if the frame header exists; judging whether the data of the packet has a frame header, if so, recording the type of the frame header and the initial position of frame data; recording the last N-1 bytes of the data of the packet; and storing the frame into the corresponding frame buffer area according to the record, and recording the frame buffer area into which the last frame data is stored.
5. The GNSS occultation data ground inversion method of claim 1, wherein: the process data and the result data in the step S6 are saved as NC format files.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010933835.8A CN112087255B (en) | 2020-09-08 | 2020-09-08 | Ground inversion method for GNSS occultation data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010933835.8A CN112087255B (en) | 2020-09-08 | 2020-09-08 | Ground inversion method for GNSS occultation data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112087255A true CN112087255A (en) | 2020-12-15 |
CN112087255B CN112087255B (en) | 2022-04-08 |
Family
ID=73731602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010933835.8A Active CN112087255B (en) | 2020-09-08 | 2020-09-08 | Ground inversion method for GNSS occultation data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112087255B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130016738A1 (en) * | 2010-03-16 | 2013-01-17 | ZTE Coration | Method and apparatus for implementing frame header alignment and multi-frame zeroing |
KR101271006B1 (en) * | 2012-12-03 | 2013-06-04 | 한국항공우주연구원 | Transmission method of gnss data for real time processing |
CN103605140A (en) * | 2013-11-06 | 2014-02-26 | 航天恒星科技有限公司 | Occultation signal rapid recovery method |
CN104486335A (en) * | 2014-12-12 | 2015-04-01 | 北京控制工程研究所 | Real-time intelligent processing system for GNC data of high-complexity spacecraft |
US20160134501A1 (en) * | 2014-11-07 | 2016-05-12 | Nomad Spectrum Limited | Data traffic processing and analysis |
CN108390890A (en) * | 2018-03-27 | 2018-08-10 | 吉旗(成都)科技有限公司 | A kind of packet splicing, the method for packaging of agreement transmission |
CN111596328A (en) * | 2020-06-30 | 2020-08-28 | 天津云遥宇航科技有限公司 | Data time sequence control method of GNSS masker receiver |
-
2020
- 2020-09-08 CN CN202010933835.8A patent/CN112087255B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130016738A1 (en) * | 2010-03-16 | 2013-01-17 | ZTE Coration | Method and apparatus for implementing frame header alignment and multi-frame zeroing |
KR101271006B1 (en) * | 2012-12-03 | 2013-06-04 | 한국항공우주연구원 | Transmission method of gnss data for real time processing |
CN103605140A (en) * | 2013-11-06 | 2014-02-26 | 航天恒星科技有限公司 | Occultation signal rapid recovery method |
US20160134501A1 (en) * | 2014-11-07 | 2016-05-12 | Nomad Spectrum Limited | Data traffic processing and analysis |
CN104486335A (en) * | 2014-12-12 | 2015-04-01 | 北京控制工程研究所 | Real-time intelligent processing system for GNC data of high-complexity spacecraft |
CN108390890A (en) * | 2018-03-27 | 2018-08-10 | 吉旗(成都)科技有限公司 | A kind of packet splicing, the method for packaging of agreement transmission |
CN111596328A (en) * | 2020-06-30 | 2020-08-28 | 天津云遥宇航科技有限公司 | Data time sequence control method of GNSS masker receiver |
Non-Patent Citations (2)
Title |
---|
洪振杰 等: "非球对称电离层掩星数据反演", 《物理学报》 * |
赵莹: "GNSS电离层掩星反演技术及应用研究", 《中国优秀博硕士学位论文全文数据库(博士)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112087255B (en) | 2022-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111949616A (en) | Ground real-time inversion demonstration system for GNSS occultation data | |
CN101470012B (en) | Method for real-time vehicle driving mileage statistics based on wireless network and GPS position information | |
AU3739800A (en) | Fault data synchronization via peer-to-peer communications network | |
CN110191109B (en) | Message sampling method and device | |
CN112087255B (en) | Ground inversion method for GNSS occultation data | |
CN103065466A (en) | Detection method of and detection device for abnormal traffic conditions | |
US20170214454A1 (en) | Efficient aeronautical information download coupled with ads-b function | |
CN111970043A (en) | Text generation method and device | |
CN110058270B (en) | Navigation deception signal generation method based on clock error fitting | |
CN101794510A (en) | Data processing method and device of floating cars | |
US7382780B1 (en) | Method for time coding of asynchronous data transmissions | |
CN109959903B (en) | Device and method for detecting radar data transmission time delay on line | |
US20100049381A1 (en) | System and Method for Simultaneously Processing Telemetry Data | |
CN114545463A (en) | Real-time data quality analysis method and system based on Beidou foundation enhancement system | |
WO2022110280A1 (en) | Method and apparatus for updating local time, device, and storage medium | |
JP4791527B2 (en) | Data processing in distributed communication networks. | |
CN112445156B (en) | Method and device for determining vehicle fault emergency degree and remote monitoring platform | |
JP2987435B2 (en) | Geographic weather mobile monitoring system | |
CN110907953A (en) | Satellite fault identification method and device and software receiver | |
CN114966695B (en) | Digital twin image processing method, device, equipment and medium for radar | |
JP4767314B2 (en) | Data processing in distributed communication networks. | |
KR20070045809A (en) | Method for processing dump data packet of low earth orbit satellite | |
CN112202761B (en) | Fork truck positioning monitoring system | |
CN111103597B (en) | Data processing method and satellite communication system | |
CN106405580A (en) | GNSS continuity assessment method |
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 |