CN112087255A - Ground inversion method for GNSS occultation data - Google Patents

Ground inversion method for GNSS occultation data Download PDF

Info

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
Application number
CN202010933835.8A
Other languages
Chinese (zh)
Other versions
CN112087255B (en
Inventor
李一路
李峰辉
李兴国
王鹏程
温凯
付乃锋
刘永成
黄满义
高阳
张小飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin Yunyao Aerospace Technology Co ltd
Original Assignee
Tianjin Yunyao Aerospace Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin Yunyao Aerospace Technology Co ltd filed Critical Tianjin Yunyao Aerospace Technology Co ltd
Priority to CN202010933835.8A priority Critical patent/CN112087255B/en
Publication of CN112087255A publication Critical patent/CN112087255A/en
Application granted granted Critical
Publication of CN112087255B publication Critical patent/CN112087255B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/185Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851Systems using a satellite or space-based relay
    • H04B7/18513Transmission in a satellite or space-based system
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L69/00Network arrangements, protocols or services independent of the application payload and not provided for in the other groups of this subclass
    • H04L69/22Parsing 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

Ground inversion method for GNSS occultation data
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.
CN202010933835.8A 2020-09-08 2020-09-08 Ground inversion method for GNSS occultation data Active CN112087255B (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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