CN116383191A - Method for cleaning, treating and quality checking and evaluating mass data of CORS station network - Google Patents

Method for cleaning, treating and quality checking and evaluating mass data of CORS station network Download PDF

Info

Publication number
CN116383191A
CN116383191A CN202310648390.2A CN202310648390A CN116383191A CN 116383191 A CN116383191 A CN 116383191A CN 202310648390 A CN202310648390 A CN 202310648390A CN 116383191 A CN116383191 A CN 116383191A
Authority
CN
China
Prior art keywords
data
frequency
observation
beidou
station network
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.)
Pending
Application number
CN202310648390.2A
Other languages
Chinese (zh)
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.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute of CAS
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 Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202310648390.2A priority Critical patent/CN116383191A/en
Publication of CN116383191A publication Critical patent/CN116383191A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/33Multimode operation in different systems which transmit time stamped messages, e.g. GPS/GLONASS
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/10File systems; File servers
    • G06F16/17Details of further file system functions
    • G06F16/174Redundancy elimination performed by the file system
    • G06F16/1748De-duplication implemented within the file system, e.g. based on file segments
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/21Design, administration or maintenance of databases
    • G06F16/215Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/25Integrating or interfacing systems involving database management systems
    • G06F16/258Data format conversion from or to a database
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Quality & Reliability (AREA)
  • Data Exchanges In Wide-Area Networks (AREA)
  • General Factory Administration (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method for cleaning, managing and quality checking and evaluating mass data of a CORS station network, which is oriented to mass GNSS monitoring and sensing data generated by a continuously running reference station system (Continuous Operational Reference System, CORS). The method specifically comprises the following steps: monitoring station network data file format standardization and normalization; monitoring station network data primary cleaning; monitoring station network data management and quality check; and (5) monitoring the quality evaluation and grading evaluation of the station network data assets. The invention aims at mass data characteristics and data flow attributes of the CORS station network, carries out data cleaning, treatment processing, hierarchical management and control and the like in a layered manner, forms a unified and executable data cleaning treatment standard, realizes a repeatable data management flow and has transparent flow management, thereby ensuring the data quality of the CORS station network and promoting the creation of application value of data assets close to business scenes.

Description

Method for cleaning, treating and quality checking and evaluating mass data of CORS station network
Technical Field
The invention relates to the technical field of satellite data processing, in particular to a method for cleaning, managing and quality checking and evaluating mass data of a (Continuous Operational Reference System, CORS) station network.
Background
The city continuous operation reference station system (CORS) can acquire the position and time information of various spaces and relevant dynamic changes, is one of hot spots of modern GNSS, and is established or is being established in more developed countries in the world at present.
With the establishment of more and more CORS stations, the generated data volume is larger and larger, and as the key of the whole CORS system to play a role, the evaluation of the data quality of the CORS is a key problem. In addition, the Beidou/GNSS receivers of different manufacturers have different design standards at present, so that the adopted received file formats are different, and the data processing of different receivers does not have format uniformity and universal applicability, so that the method has important significance for unifying the received data formats.
Disclosure of Invention
In view of the foregoing, the present invention has been developed to provide a method for cleaning and quality inspection and assessment of mass data of a CORS site network that overcomes or at least partially solves the foregoing problems.
According to one aspect of the invention, a method for cleaning, managing and evaluating mass data and quality check of a CORS station network is provided, which comprises the following steps:
s1, monitoring station network data file format standardization and normalization processing;
s2, monitoring primary cleaning of station network data;
s3, monitoring station network data management and quality check;
and S4, monitoring the quality evaluation and the grading evaluation of the station network data assets.
The step S1 specifically includes:
according to different Beidou/GNSS high-precision measurement type receiver manufacturers and different parameter settings of CORS station network receivers, one version of Rinex2.X series and Rinex3.X series can be selected as a data file storage format; the data format of multiple Rinex versions necessarily introduces the trouble of data import processing, and file format normalization and normalization processing are necessary.
The file format normalization and normalization process mainly comprises the following steps: a format checking and repairing and format converting function;
format checking and repairing: the method mainly comprises the steps of performing normalization inspection on RINEX format files, attempting automatic repair on incorrect format contents, and updating header file information according to inspection results to ensure consistency of the format and the contents;
format conversion: converting the format of any data which accords with the RINEX format standard, and normalizing the Rinex2.X format file into the RINEX2.11 format standard; the RINEX3.x format file is normalized to the RINEX3.04 format standard.
Further, the step S2 specifically includes:
the CORS station network data file is inevitably subjected to various conditions of data repetition, deletion, error and data redundancy, and special service software operation can be supported by cleaning, otherwise, the service software operation thread is extremely easy to terminate and even crash and the like.
The primary cleaning function of the monitoring station network data mainly comprises the following steps: redundant data deletion, error data correction, duplicate/invalid data deletion, logical missing data filling, and the like, specifically:
redundant data deletion: and deleting unnecessary satellite navigation system observation data in the observation file according to the subsequent service data requirement, such as only keeping Beidou/GPS observation data, deleting unnecessary navigation text files and the like.
Error data correction: checking the type of a receiver, the type of an antenna, the height of the antenna and the like in a file by using CORS station network element data, and correcting error information possibly existing; and checking the rationality of the observation epoch, satellite number and the like in the file body and attempting to correct.
Duplicate/invalid data deletion: and performing traversal checking on the data in the file according to the file header information, and identifying and deleting repeated and invalid data.
Logic missing data padding: according to the Rinex file standard, padding and filling are attempted according to logic or rules for partially missing data.
The step S3 specifically includes: data integrity rate analysis, data efficiency analysis, cycle slip rate analysis, multipath error analysis, pseudo-range noise analysis, carrier phase noise analysis and carrier-to-noise ratio analysis.
Further, the data integrity rate analysis specifically includes:
calculating any system according to a formula (1)
Figure SMS_1
Any frequency signal->
Figure SMS_2
Calculating the observation data integrity rate of any system according to a formula (2);
Figure SMS_3
(1)
Figure SMS_4
(2)
in the method, in the process of the invention,
Figure SMS_10
for Beidou/GNSS system->
Figure SMS_8
Frequency signal +.>
Figure SMS_17
For observing the data integrity rate, the unit is; />
Figure SMS_9
The total number of satellites observed in the observation period; />
Figure SMS_15
For observation guardStar sequence number->
Figure SMS_16
;/>
Figure SMS_25
For the observation period, beidou/GNSS system +.>
Figure SMS_7
Satellite->
Figure SMS_13
In frequency signal->
Figure SMS_12
Is a total number of actual observation epochs; />
Figure SMS_22
For the observation period, beidou/GNSS system +.>
Figure SMS_5
Satellite->
Figure SMS_26
In frequency signal->
Figure SMS_11
Is a theoretical epoch count of (2); />
Figure SMS_14
For system->
Figure SMS_19
Observing the data integrity rate in units of;
Figure SMS_23
for the observation period, beidou/GNSS system +.>
Figure SMS_21
Satellite->
Figure SMS_24
Actual observation epoch Total with all observation frequencies having observation dataA number; />
Figure SMS_6
For the observation period, beidou/GNSS system +.>
Figure SMS_18
Satellite->
Figure SMS_20
Is a theoretical epoch count of (a).
Further, the data efficiency analysis specifically includes:
the effective rate can be defined by equation (3):
Figure SMS_27
(3)
wherein,,
Figure SMS_28
representing a high cut-off angle->
Figure SMS_29
Observing the data quantity above the degree;
Figure SMS_30
representing the observed data volume of unhealthy satellites, wherein +.>
Figure SMS_31
The signal to noise ratio of the data is less than the specified threshold epoch number.
Further, cycle slip analysis specifically includes:
any Beidou/GNSS system in the observation period is calculated according to a formula (4)
Figure SMS_32
Cycle slip ratio of (c):
Figure SMS_33
(4)
in the method, in the process of the invention,
Figure SMS_34
for the observation period, beidou/GNSS system +.>
Figure SMS_35
Cycle slip ratio of (2); />
Figure SMS_36
For Beidou/GNSS system->
Figure SMS_37
The total number of actual epochs observed during the observation period; />
Figure SMS_38
For Beidou/GNSS system->
Figure SMS_39
The total number of cycle slip epochs occurring during the observation period is also known as Zhou Tiaoshu. Zhou Tiaoshu->
Figure SMS_40
A method of MW (Melbourne-Mubbena) combined detection and GF (Geometry-Free) combined detection is adopted.
Further, the multipath error analysis specifically includes:
calculating the multipath error value RMS of any Beidou/GNSS system, any frequency and any satellite in the observation period according to a formula (5):
Figure SMS_41
(5)
in the method, in the process of the invention,
Figure SMS_44
for observation frequency in observation period->
Figure SMS_46
Multipath error values RMS in meters (m);
Figure SMS_49
the total number of the calendar elements is observed in the observation period; />
Figure SMS_43
For epoch number,/->
Figure SMS_47
;/>
Figure SMS_50
For observing frequency +.>
Figure SMS_51
In epoch->
Figure SMS_42
Multipath calculations (containing integer ambiguity effects) at time in meters (m); />
Figure SMS_45
For observing frequency +.>
Figure SMS_48
The average value in meters (m) is calculated for the multipath over the observation period.
Calendar element
Figure SMS_52
Multipath calculation value +.>
Figure SMS_53
I.e. < ->
Figure SMS_54
Calculated according to the formula (6):
Figure SMS_55
(6)
in the method, in the process of the invention,
Figure SMS_56
is->
Figure SMS_60
Multipath calculation value of any observation epoch of frequency, the unit is meter (m); />
Figure SMS_62
The unit is meter (m) for the pseudo-range observation of the epoch corresponding to the first frequency; />
Figure SMS_57
A carrier frequency that is the first frequency in megahertz (MHz);
Figure SMS_61
a carrier frequency in megahertz (MHz) that is the second frequency; />
Figure SMS_63
The unit is meter (m) for the observed carrier phase of the epoch corresponding to the first frequency; />
Figure SMS_64
The unit is meter (m) for the observed carrier phase of the epoch corresponding to the second frequency; />
Figure SMS_58
Is that
Figure SMS_59
Multipath calculation value of any observation epoch of frequency, the unit is meter (m); />
Figure SMS_65
The unit is meter (m) for the pseudo-range observation of the epoch corresponding to the second frequency;
Figure SMS_66
calculated according to the formula (7):
Figure SMS_67
(7)
in the method, in the process of the invention,
Figure SMS_68
the total number of the calendar elements is observed in the observation period; />
Figure SMS_69
For epoch number,/->
Figure SMS_70
Specifically, the pseudo-range noise analysis specifically includes:
calculating pseudo-range noise average values of all observation satellites of any frequency of any Beidou/GNSS system in the observation period according to a formula (8):
Figure SMS_71
(8)
in the method, in the process of the invention,
Figure SMS_72
in the observation period, the unit of any frequency pseudo-range noise average value of any Beidou/GNSS system is meter (m); />
Figure SMS_73
In the observation period, the total number of satellites observed by any frequency of any Beidou/GNSS system; />
Figure SMS_74
In order to observe the satellite serial number,
Figure SMS_75
;/>
Figure SMS_76
for the observation period, any frequency (signal) of any Beidou/GNSS system is +.>
Figure SMS_77
Pseudo-range noise RMS for a satellite is in meters (m).
Figure SMS_78
Calculated according to the formula (9):
Figure SMS_79
(9)
wherein:
Figure SMS_80
in order to observe the period, any frequency of any Beidou/GNSS system is +.>
Figure SMS_81
The total number of observation epochs of the satellites; />
Figure SMS_82
For epoch number,/->
Figure SMS_83
;/>
Figure SMS_84
For any frequency of any Beidou/GNSS system +.>
Figure SMS_85
Satellite, in epoch->
Figure SMS_86
Pseudo-range noise estimation at time in meters (m);
Figure SMS_87
calculated according to the formula (10):
Figure SMS_88
(10)
in the method, in the process of the invention,
Figure SMS_89
for any frequency of any Beidou/GNSS system +.>
Figure SMS_90
Satellite, in epoch->
Figure SMS_91
The unit of the pseudo-range observation value of the moment is meter (m); />
Figure SMS_92
For any frequency of any Beidou/GNSS system +.>
Figure SMS_93
Satellite, in epoch->
Figure SMS_94
The pseudorange quadratic polynomial fit value of time is given in meters (m).
Further: performing pseudo-range quadratic polynomial fitting according to a formula (11); after fitting calculation to obtain a quadratic polynomial coefficient, a pseudo-range quadratic polynomial fitting value of each epoch can be obtained
Figure SMS_95
The method comprises the steps of carrying out a first treatment on the surface of the In the observation period, starting from the initial epoch, each 120 epochs are a fitting window, and the fitting windows are not overlapped; fitting according to the number of the remaining calendar elements when the number of the remaining calendar elements is not less than 3 calendar element observation data near the end of the observation period; when the number of the remaining epoch is smaller than 3 epoch observation data, the method can be abandoned; when the observation data has an interruption phenomenon, the method can be used for respectively processing the following steps:
Figure SMS_96
(11)
in the method, in the process of the invention,
Figure SMS_97
fitting a function for a pseudo-range quadratic polynomial; />
Figure SMS_98
To fit epoch number within window, quadratic polynomial argument, ++>
Figure SMS_99
;/>
Figure SMS_100
Is a quadratic polynomial coefficient; />
Figure SMS_101
A first order term coefficient which is a quadratic polynomial;
Figure SMS_102
is a quadratic polynomial constant term.
Specifically, the carrier phase noise analysis specifically includes:
calculating the carrier-phase noise average value of all observation satellites at any frequency of any Beidou/GNSS system in the observation period according to a formula (12):
Figure SMS_103
(12)
in the method, in the process of the invention,
Figure SMS_104
in the observation period, the average value of phase noise of any frequency carrier of any Beidou/GNSS system is given in units of weeks; />
Figure SMS_105
In the observation period, the total number of satellites observed by any frequency of any Beidou/GNSS system; />
Figure SMS_106
In order to observe the satellite serial number,
Figure SMS_107
;/>
Figure SMS_108
in order to observe the period, any frequency of any Beidou/GNSS system is +.>
Figure SMS_109
Carrier phase noise of a satellite is in units of weeks.
Further, according to the formula (13)
Figure SMS_110
Figure SMS_111
(13)
Wherein:
Figure SMS_112
in order to observe the period, any frequency of any Beidou/GNSS system is +.>
Figure SMS_113
The total number of observation epochs of the satellites;
Figure SMS_114
for epoch number,/->
Figure SMS_115
;/>
Figure SMS_116
For any frequency of any Beidou/GNSS system +.>
Figure SMS_117
And the noise estimation of three differences of the carrier phase observation values of adjacent epochs of the satellites is given in units of weeks.
Calculated according to the formula (14)
Figure SMS_118
Figure SMS_119
(14)
In the method, in the process of the invention,
Figure SMS_120
for epoch->
Figure SMS_121
Moment, any satellite->
Figure SMS_122
The carrier phase observations at any frequency are in units of weeks.
Specifically, the carrier-to-noise ratio analysis specifically includes:
calculating the average value of the carrier-to-noise ratio statistics of all satellites at any frequency of any Beidou/GNSS system according to a formula (15):
Figure SMS_123
(15)
in the method, in the process of the invention,
Figure SMS_124
the average value of the carrier-to-noise ratio statistics values of all satellites at any frequency of any Beidou/GNSS system is expressed in decibel hertz (dBHz); />
Figure SMS_125
Observing the total number of satellites in an observation period; />
Figure SMS_126
For observing satellite serial number>
Figure SMS_127
;/>
Figure SMS_128
In order to be within the observation period, any frequency of any Beidou/GNSS system is +.>
Figure SMS_129
The average carrier-to-noise ratio of a satellite is expressed in decibel hertz (dBHz).
Further, calculate according to equation (16)
Figure SMS_130
Figure SMS_131
(16)
Wherein:
Figure SMS_132
the total number of the calendar elements is observed in the observation period; />
Figure SMS_133
For epoch number,/->
Figure SMS_134
;/>
Figure SMS_135
For epoch->
Figure SMS_136
Moment, satellite->
Figure SMS_137
The observed amount of carrier-to-noise ratio at any frequency is in decibel hertz (dBHz).
Further, the step S4 specifically includes:
monitoring station network data asset quality assessment and grading assessment, wherein indexes such as data integrity rate, data effective rate, pseudo-range multipath, cycle slip ratio and signal to noise ratio are mainly used as indexes of a comprehensive evaluation model of a sequencing method TOPSIS (Technique for Order Preference by Similarity to Ideal Solution) approaching ideal points; the comprehensive evaluation process is to 'synthesize' a plurality of evaluation indexes into an integral comprehensive evaluation index through a mathematical model, and combine the special navigation time-frequency service type to determine the weight coefficient by forward and dimensionless methods and entropy value methods on the indexes, and determine the final comprehensive evaluation model so as to obtain the final evaluation result; and (3) carrying out grade assessment on the observation data file by utilizing the comprehensive evaluation result of the observation data and referring to a set threshold value, wherein the grade assessment is respectively four grades of excellent, good, qualified and unavailable.
According to the resources and technical characteristics utilized by the invention, the invention has the following technical advantages:
A. constructing a set of transparent, standard and standardized CORS station network data cleaning treatment flow
Aiming at mass data characteristics and data flow attributes of the CORS station network, data cleaning, treatment processing, hierarchical management and control and the like are carried out in a layered manner, a unified and executable data cleaning treatment standard is formed, a repeatable data management flow is realized, and flow management is transparent, so that the data quality of the CORS station network is ensured, and the creation of application value of data assets close to business scenes is promoted.
B. Forming a comprehensive data quality comprehensive evaluation system
The method for realizing the multidimensional evaluation of data quality by monitoring the service demand of the perceived data scene facing the CORS station network comprises the following steps: observation data integrity rate analysis, availability analysis, pseudo-range observation noise analysis, carrier phase observation noise analysis, signal-to-noise ratio analysis, cycle slip ratio analysis, multipath influence analysis and the like, and comprehensive multidimensional evaluation conclusion realizes hierarchical screening and management and control of data/site resources, and well corresponds to multi-type service requirements.
C. Ensuring balance of data governance benefits and costs
The CORS station network monitors and perceives data to generate massive data at real time, and the management of the massive data does not need to pursue extremely high data quality.
Drawings
In order that the above-recited objects, features and advantages of the present application will become more apparent and fully apparent from the following detailed description of embodiments of the invention, it should be read in connection with the accompanying drawings. Based on the embodiments of the present invention, those skilled in the art may implement other embodiments without making any inventive effort, which fall within the scope of the present invention.
FIG. 1 shows a flow chart of a method for cleaning and treating mass data and evaluating quality check of a CORS station network.
Detailed Description
The invention provides a method for cleaning, treating and quality checking and evaluating mass data of a CORS station network, which aims to overcome the defects of various formats and different data quality of the received data of the CORS station and solve the increasing demands for standardization and analysis and treatment of the CORS data.
The method realizes data standardization and cleaning, data quality analysis and evaluation and the like of the CORS station network observation data with the characteristics of high updating frequency, strong real-time performance, large data quantity, uncertain quality and the like, and forms the data quality checking capability and deep excavation data potential which cover comprehensively, have transparent and standard flows and are unified and feasible.
According to the embodiment of the invention, the method for cleaning, treating and quality checking and evaluating mass data of the CORS station network comprises the following steps:
s1, monitoring station network data file format standardization and normalization processing;
s2, monitoring primary cleaning of station network data;
s3, monitoring station network data management and quality check;
and S4, monitoring the quality evaluation and the grading evaluation of the station network data assets.
As shown in fig. 1, the detailed technical scheme of the steps is as follows:
1. data standardization and primary cleaning module
From the purpose of unifying (or reducing) the format version of the observation data file of the CORS station network, carrying out standardized normalization processing on the observation data files of different types, wherein the standardized normalization processing mainly comprises format checking and repairing, normalization format conversion and the like; the primary cleaning of the content in the observation data file mainly comprises repeated data deletion, logic missing data filling, error data identification and correction, invalid data elimination and the like in the file, and the readability and usability of the data in the file are ensured.
1.1 normalization and normalization treatments
The format checking and repairing functions are aimed at Beidou/GNSS observation data files (Obs), broadcast ephemeris files (Nav) and the like, and due to the fact that Beidou/GNSS receiver manufacturers are different and parameter settings are different, file format versions are inconsistent, and a plurality of inconveniences are introduced for later data processing and use.
The file format normalization and normalization process mainly comprises the following steps: and (5) a format checking and repairing and format conversion function. Format checking and repairing: the method mainly comprises the steps of performing normalization inspection on RINEX format files, attempting automatic repair on incorrect format contents, and updating header file information according to inspection results to ensure consistency of the format and the contents; format conversion: converting the format of any data which accords with the RINEX format standard, and normalizing the Rinex2.X format file into the RINEX2.11 format standard; the RINEX3.x format file is normalized to the RINEX3.04 format standard.
The step S2 specifically includes:
1.2 data Primary cleaning
The primary cleaning function of the monitoring station network data mainly comprises the following steps: redundant data deletion, error data correction, duplicate/invalid data deletion, logical missing data filling, and the like, specifically:
redundant data deletion: and deleting unnecessary satellite navigation system observation data in the observation file according to the subsequent service data requirement, such as only keeping Beidou/GPS observation data, deleting unnecessary navigation text files and the like.
Error data correction: checking the type of a receiver, the type of an antenna, the height of the antenna and the like in a file by using CORS station network element data, and correcting error information possibly existing; and checking the rationality of the observation epoch, satellite number and the like in the file body and attempting to correct.
Duplicate/invalid data deletion: and performing traversal checking on the data in the file according to the file header information, and identifying and deleting repeated and invalid data.
Logic missing data padding: according to the Rinex file standard, padding and filling are attempted according to logic or rules for partially missing data.
2 monitoring station network data management and quality inspection
2.1 data integrity Rate analysis
Calculating any system according to a formula (1)
Figure SMS_138
Any frequency (signal)>
Figure SMS_139
The observation data integrity rate of any system is calculated according to the formula (2).
Figure SMS_140
(1)
Figure SMS_141
(2)
In the method, in the process of the invention,
Figure SMS_158
for Beidou/GNSS system->
Figure SMS_143
Frequency (signal)/(frequency)>
Figure SMS_161
For observing the data integrity rate, the unit is; />
Figure SMS_145
The total number of satellites observed in the observation period; />
Figure SMS_155
For observing satellite serial number>
Figure SMS_146
;/>
Figure SMS_151
For the observation period, beidou/GNSS system +.>
Figure SMS_160
Satellite->
Figure SMS_163
At frequency (signal)/(frequency)>
Figure SMS_142
Is a total number of actual observation epochs; />
Figure SMS_156
For the observation period, beidou/GNSS system +.>
Figure SMS_147
Satellite->
Figure SMS_153
At frequency (signal)/(frequency)>
Figure SMS_159
Is a theoretical epoch count of (2); />
Figure SMS_162
For system->
Figure SMS_148
Observing the data integrity rate in units of; />
Figure SMS_154
For the observation period, beidou/GNSS system +.>
Figure SMS_149
Satellite->
Figure SMS_152
The total number of actual observation epochs of the observation data exists in all the observation frequencies (signals); />
Figure SMS_144
For the observation period, beidou/GNSS system +.>
Figure SMS_150
Satellite->
Figure SMS_157
Is a theoretical epoch count of (a).
2.2 data efficient analysis
The effective rate can be defined by equation (3):
Figure SMS_164
(3)
wherein,,
Figure SMS_165
indicating the height cut-off angle
Figure SMS_166
Observing the data quantity above the degree;
Figure SMS_167
representing the observed data volume for unhealthy satellites, wherein,
Figure SMS_168
the signal to noise ratio of the data is less than the specified threshold epoch number.
2.3 cycle skip analysis
Any Beidou/GNSS system in the observation period is calculated according to a formula (4)
Figure SMS_169
Cycle slip ratio of (c):
Figure SMS_170
(4)
in the method, in the process of the invention,
Figure SMS_171
for the observation period, beidou/GNSS system +.>
Figure SMS_172
Cycle slip ratio of (2); />
Figure SMS_173
For Beidou/GNSS system->
Figure SMS_174
The total number of actual epochs observed during the observation period; />
Figure SMS_175
For Beidou/GNSS system->
Figure SMS_176
The total number of cycle slip epochs occurring during the observation period (Zhou Tiaoshu). Zhou Tiaoshu->
Figure SMS_177
A method of MW (Melbourne-Mubbena) combined detection and GF (Geometry-Free) combined detection is adopted.
2.4 multipath error analysis
And calculating the multipath error value (RMS) of any Beidou/GNSS system, any frequency (signal) and any satellite in the observation period according to the formula (5).
Figure SMS_178
(5)
In the method, in the process of the invention,
Figure SMS_181
for the observation period, the observation frequency (signal)>
Figure SMS_184
Multipath error values (RMS) in meters (m); />
Figure SMS_186
The total number of the calendar elements is observed in the observation period; />
Figure SMS_179
For epoch number,/->
Figure SMS_182
;/>
Figure SMS_185
For observing frequency (signal)>
Figure SMS_188
In epoch->
Figure SMS_180
Multipath calculations (containing integer ambiguity effects) at time in meters (m); />
Figure SMS_183
For observing frequency (signal)>
Figure SMS_187
The average value (containing integer ambiguity effects) is calculated in multipath over the observation period in meters (m).
Calendar element
Figure SMS_189
Multipath calculation value +.>
Figure SMS_190
I.e. < ->
Figure SMS_191
Calculated according to the formula (6):
Figure SMS_192
(6)
in the method, in the process of the invention,
Figure SMS_194
is->
Figure SMS_198
Multipath calculation value (including integer ambiguity influence) of any observation epoch of frequency (signal) in meters (m); />
Figure SMS_201
The unit is meter (m) for the pseudo-range observation of the epoch corresponding to the first frequency (signal); />
Figure SMS_195
A carrier frequency in megahertz (MHz) that is the first frequency (signal); />
Figure SMS_197
A carrier frequency in megahertz (MHz) which is the second frequency (signal); />
Figure SMS_199
The unit is meter (m) for the carrier phase observed quantity of the epoch corresponding to the first frequency (signal); />
Figure SMS_202
The unit is meter (m) for the carrier phase observed quantity of the epoch corresponding to the second frequency (signal); />
Figure SMS_193
Is->
Figure SMS_196
Multipath calculation of any observation epoch for frequency (signal) (including integer ambiguity effect) in units ofRice (m); />
Figure SMS_200
The unit is meter (m) for the pseudorange observations of the epoch corresponding to the second frequency (signal).
Figure SMS_203
Calculated according to the formula (7):
Figure SMS_204
(7)
in the method, in the process of the invention,
Figure SMS_205
the total number of the calendar elements is observed in the observation period; />
Figure SMS_206
For epoch number,/->
Figure SMS_207
The method comprises the steps of carrying out a first treatment on the surface of the In the process of multi-path error calculation, firstly, the influence of coarse error, clock jump and cycle slip of a repair receiver and the like should be removed from observed data.
2.5 pseudo-range noise analysis
And (3) calculating pseudo-range noise (average value) of all observation satellites of any frequency (signal) of any Beidou/GNSS system in the observation period according to a formula (8).
Figure SMS_208
(8)
In the method, in the process of the invention,
Figure SMS_209
in order to observe the time interval, the pseudo-range noise (average value) of any frequency (signal) of any Beidou/GNSS system is expressed in meters (m); />
Figure SMS_210
The total number of satellites observed by any frequency (signal) of any Beidou/GNSS system in an observation period; />
Figure SMS_211
For observing satellite serial number>
Figure SMS_212
;/>
Figure SMS_213
In order to observe the period, any frequency (signal) of any Beidou/GNSS system
Figure SMS_214
Pseudo-range noise (RMS) of a satellite in meters (m).
Figure SMS_215
Calculated according to the formula (9):
Figure SMS_216
(9)
wherein:
Figure SMS_217
for the observation period, any frequency (signal) of any Beidou/GNSS system is +.>
Figure SMS_218
The total number of observation epochs of the satellites; />
Figure SMS_219
For epoch number,/->
Figure SMS_220
;/>
Figure SMS_221
For any frequency (signal) of any Beidou/GNSS system +.>
Figure SMS_222
Satellite, in epoch->
Figure SMS_223
Pseudo-range noise estimates of the time of day,the unit is meter (m).
Figure SMS_224
Calculated according to the formula (10):
Figure SMS_225
(10)
in the method, in the process of the invention,
Figure SMS_226
for any frequency (signal) of any Beidou/GNSS system +.>
Figure SMS_227
Satellite, in epoch->
Figure SMS_228
The pseudorange observations at time are in meters (m). />
Figure SMS_229
For any frequency (signal) of any Beidou/GNSS system +.>
Figure SMS_230
Satellite, in epoch->
Figure SMS_231
The pseudorange quadratic polynomial fit value of time is given in meters (m).
Performing pseudo-range quadratic polynomial fitting according to a formula (11); after fitting calculation to obtain a quadratic polynomial coefficient, a pseudo-range quadratic polynomial fitting value of each epoch can be obtained
Figure SMS_232
. During the observation period, starting from the starting epoch, every 120 epochs is a fit window, which does not overlap. Fitting according to the number of the remaining calendar elements when the number of the remaining calendar elements is not less than 3 calendar element observation data near the end of the observation period; when the number of remaining epochs is less than 3 epoch observation data, it may be discarded. When the observation data has an interruption phenomenon, the observation data can be processed in a sectionalized way.
Figure SMS_233
(11)
In the method, in the process of the invention,
Figure SMS_234
fitting a function for a pseudo-range quadratic polynomial; />
Figure SMS_235
To fit epoch number within window, quadratic polynomial argument, ++>
Figure SMS_236
;/>
Figure SMS_237
Is a quadratic polynomial coefficient; />
Figure SMS_238
A first order term coefficient which is a quadratic polynomial;
Figure SMS_239
is a quadratic polynomial constant term.
2.6 Carrier phase noise analysis
And (3) calculating carrier phase noise (average value) of all observation satellites of any frequency (signal) of any Beidou/GNSS system in the observation period according to a formula (12).
Figure SMS_240
(12)
In the method, in the process of the invention,
Figure SMS_241
in order to observe the period, the carrier phase noise (average value) of any frequency (signal) of any Beidou/GNSS system is given in units of weeks; />
Figure SMS_242
The total number of satellites observed by any frequency (signal) of any Beidou/GNSS system in an observation period; />
Figure SMS_243
For observing satellite serial number>
Figure SMS_244
;/>
Figure SMS_245
For the observation period, any frequency (signal) of any Beidou/GNSS system is +.>
Figure SMS_246
Carrier phase noise of a satellite is in units of weeks.
Calculated according to the formula (13)
Figure SMS_247
Figure SMS_248
(13)
Wherein:
Figure SMS_249
for the observation period, any frequency (signal) of any Beidou/GNSS system is +.>
Figure SMS_250
The total number of observation epochs of the satellites; />
Figure SMS_251
For epoch number,/->
Figure SMS_252
;/>
Figure SMS_253
For any frequency (signal) of any Beidou/GNSS system +.>
Figure SMS_254
And the noise estimation of three differences of the carrier phase observation values of adjacent epochs of the satellites is given in units of weeks.
Calculated according to the formula (14)
Figure SMS_255
Figure SMS_256
(14)
In the method, in the process of the invention,
Figure SMS_257
for epoch->
Figure SMS_258
Moment, any satellite->
Figure SMS_259
Carrier phase observations of any frequency (signal) are in cycles. In carrier phase noise calculation, the influence of coarse difference, clock jump and cycle slip of the repair receiver should be removed from the observed data. When the observation data has an interruption phenomenon, the observation data can be processed in a sectionalized way.
2.7 analysis of the Carrier to noise ratio
And calculating the carrier-to-noise ratio statistic (average value) of all satellites at any frequency of any Beidou/GNSS system according to a formula (15).
Figure SMS_260
(15)
In the method, in the process of the invention,
Figure SMS_261
the carrier-to-noise ratio statistics (average) of all satellites at any frequency of any Beidou/GNSS system is in decibel hertz (dBHz); />
Figure SMS_262
Observing the total number of satellites in an observation period; />
Figure SMS_263
For observing satellite serial number>
Figure SMS_264
;/>
Figure SMS_265
In order to be within the observation period, any frequency of any Beidou/GNSS system is +.>
Figure SMS_266
The average carrier-to-noise ratio of a satellite is expressed in decibel hertz (dBHz).
Calculated according to the formula (16)
Figure SMS_267
Figure SMS_268
(16)
Wherein:
Figure SMS_269
the total number of the calendar elements is observed in the observation period; />
Figure SMS_270
For epoch number,/->
Figure SMS_271
;/>
Figure SMS_272
For epoch->
Figure SMS_273
Moment, satellite->
Figure SMS_274
The observed amount of carrier-to-noise ratio of any frequency (signal) is in decibel hertz (dBHz).
3. Monitoring station network data asset quality assessment and grading assessment
The comprehensive evaluation analysis of the quality of the observed data adopts TOPSIS to integrate information of a plurality of indexes to comprehensively evaluate the observed data, and the index values of all the evaluated objects are examined to be close to ideal points and sequenced in sequence, wherein the ideal points are the optimal index values in the measuring station. Because the pseudo-range noise, the carrier phase noise and the carrier phase multipath have smaller values, the pseudo-range multipath contains the information of the pseudo-range noise, the carrier phase multipath and the like, and the pseudo-range multipath is used for replacing the complex degree of the comprehensive evaluation model. And taking indexes such as data integrity rate, data effective rate, pseudo-range multipath, cycle slip ratio, signal-to-noise ratio and the like as indexes of a comprehensive evaluation model, determining weight coefficients by forward and dimensionless methods and entropy value methods, and grading observation data, wherein the four grades are excellent, good, qualified and unavailable.
While preferred embodiments of the present embodiments have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. It is therefore intended that the following claims be interpreted as including the preferred embodiments and all such alterations and modifications as fall within the scope of the embodiments of the present application.
The method for cleaning, managing and evaluating mass data and quality check of the CORS station network provided by the application is described in detail, and specific examples are applied to the explanation of the principle and the implementation mode of the application, and the explanation of the examples is only used for helping to understand the method and the core idea of the application; meanwhile, as those skilled in the art will have modifications in the specific embodiments and application scope in accordance with the ideas of the present application, the present description should not be construed as limiting the present application in view of the above.

Claims (17)

1. A mass data cleaning treatment and quality check evaluation method for a CORS station network is characterized by comprising the following steps:
s1, monitoring station network data file format standardization and normalization processing;
s2, monitoring primary cleaning of station network data;
s3, monitoring station network data management and quality check;
and S4, monitoring the quality evaluation and the grading evaluation of the station network data assets.
2. The mass data cleaning treatment and quality check evaluation method for the CORS station network as claimed in claim 1, wherein the mass data cleaning treatment and quality check evaluation method is characterized by comprising the following steps of: s1 specifically comprises:
according to different Beidou/GNSS high-precision measurement type receiver manufacturers, the parameter settings of the CORS station network receiver are different, the data file storage format has a plurality of Rinex version data formats, and file format standardization and normalization processing are needed;
the file format normalization and normalization process comprises the following steps: a format checking and repairing and format converting function;
format checking and repairing: performing normalization inspection on the RINEX format file, attempting automatic repair on incorrect format content, and updating header file information according to an inspection result to ensure consistency of the format and the content;
format conversion: performing format conversion on any data meeting the RINEX format standard, and normalizing the RINEX2.X format file into the RINEX2.11 format standard; the RINEX3.x format file is normalized to the RINEX3.04 format standard.
3. The mass data cleaning treatment and quality check evaluation method for the CORS station network as claimed in claim 1, wherein the mass data cleaning treatment and quality check evaluation method is characterized by comprising the following steps of: s2 specifically comprises:
monitoring station network data primary cleaning functions includes: redundant data deletion, error data correction, repeated data/invalid data deletion, and logical missing data filling, specifically:
redundant data deletion: deleting unnecessary satellite navigation system observation data in the observation file according to the subsequent service data requirement;
error data correction: checking the type of a receiver, the type of an antenna and the height of the antenna in a file by using CORS station network element data, and correcting error information possibly existing; rationality checking and correcting the observation epoch and satellite number in the file body;
duplicate/invalid data deletion: performing traversal inspection on the data in the file according to the file header information, and identifying and deleting repeated and invalid data;
logic missing data padding: and filling the partially missing data according to logic or rules according to the Rinex file standard.
4. The mass data cleaning treatment and quality check evaluation method for the CORS station network as claimed in claim 1, wherein the mass data cleaning treatment and quality check evaluation method is characterized by comprising the following steps of: s3 specifically comprises: data integrity rate analysis, data efficiency analysis, cycle slip rate analysis, multipath error analysis, pseudo-range noise analysis, carrier phase noise analysis and carrier-to-noise ratio analysis.
5. The mass data cleaning treatment and quality check evaluation method for the CORS station network according to claim 4, which is characterized in that: the data integrity rate analysis specifically comprises the following steps:
calculating any system according to a formula (1)
Figure QLYQS_1
Any frequency signal->
Figure QLYQS_2
Calculating the observation data integrity rate of any system according to a formula (2);
Figure QLYQS_3
(1)
Figure QLYQS_4
(2)
in the method, in the process of the invention,
Figure QLYQS_17
for Beidou/GNSS system->
Figure QLYQS_7
Frequency signal +.>
Figure QLYQS_14
For observing the data integrity rate, the unit is; />
Figure QLYQS_8
The total number of satellites observed in the observation period; />
Figure QLYQS_21
For observing satellite serial number>
Figure QLYQS_11
;/>
Figure QLYQS_24
For the observation period, beidou/GNSS system +.>
Figure QLYQS_12
Satellite->
Figure QLYQS_16
In frequency signal->
Figure QLYQS_9
Is a total number of actual observation epochs; />
Figure QLYQS_18
For the observation period, beidou/GNSS system +.>
Figure QLYQS_5
Satellite->
Figure QLYQS_13
In frequency signal->
Figure QLYQS_10
Is a theoretical epoch count of (2); />
Figure QLYQS_19
For system->
Figure QLYQS_22
Observing the data integrity rate in units of; />
Figure QLYQS_25
For the observation period, beidou/GNSS system +.>
Figure QLYQS_23
Satellite->
Figure QLYQS_26
All the observation frequencies have the total number of actual observation epochs of the observation data; />
Figure QLYQS_6
For the observation period, beidou/GNSS system +.>
Figure QLYQS_15
Satellite->
Figure QLYQS_20
Is a theoretical epoch count of (a).
6. The mass data cleaning treatment and quality check evaluation method for the CORS station network according to claim 4, which is characterized in that: the data efficiency analysis specifically comprises:
the data efficiency is defined by equation (3):
Figure QLYQS_27
(3)
wherein,,
Figure QLYQS_28
representing a high cut-off angle->
Figure QLYQS_29
Observing the data quantity above the degree; />
Figure QLYQS_30
Representing the observed data volume of unhealthy satellites, wherein +.>
Figure QLYQS_31
The signal to noise ratio of the data is less than the specified threshold epoch number.
7. The mass data cleaning treatment and quality check evaluation method for the CORS station network according to claim 4, which is characterized in that: the cycle slip analysis specifically includes:
any Beidou/GNSS system in the observation period is calculated according to a formula (4)
Figure QLYQS_32
Cycle slip ratio of (c):
Figure QLYQS_33
(4)
in the formula (4), the amino acid sequence of the compound,
Figure QLYQS_34
for the observation period, beidou/GNSS system +.>
Figure QLYQS_35
Cycle slip ratio of (2); />
Figure QLYQS_36
For Beidou/GNSS system->
Figure QLYQS_37
The total number of actual epochs observed during the observation period; />
Figure QLYQS_38
For Beidou/GNSS system->
Figure QLYQS_39
The total number of cycle slip epochs occurring during the observation period, also known as Zhou Tiaoshu; zhou Tiaoshu->
Figure QLYQS_40
Using the first MW groupAnd (3) a method for combined detection and then GF combined detection.
8. The mass data cleaning treatment and quality check evaluation method for the CORS station network according to claim 4, which is characterized in that: the multipath error analysis specifically includes:
calculating multipath RMS (Root Mean Square) error values of any Beidou/GNSS system, any frequency and any satellite in the observation period according to a formula (5):
Figure QLYQS_41
(5)
in the formula (5), the amino acid sequence of the compound,
Figure QLYQS_43
for observation frequency in observation period->
Figure QLYQS_47
Multipath RMS error values in meters (m);
Figure QLYQS_49
the total number of the calendar elements is observed in the observation period; />
Figure QLYQS_44
For epoch number,/->
Figure QLYQS_45
;/>
Figure QLYQS_50
For observing frequency +.>
Figure QLYQS_51
In the epoch
Figure QLYQS_42
The multipath calculation value of the moment contains the influence of integer ambiguity, and the unit is meter (m); />
Figure QLYQS_46
For observing frequency +.>
Figure QLYQS_48
The average value in meters is calculated for multiple paths during the observation period.
9. The mass data cleaning treatment and quality check assessment method for a CORS station network according to claim 8, wherein the mass data cleaning treatment and quality check assessment method is characterized by comprising the following steps of:
calendar element
Figure QLYQS_52
Multipath calculation value +.>
Figure QLYQS_53
I.e. +.>
Figure QLYQS_54
Calculated according to the formula (6):
Figure QLYQS_55
(6)
in the formula (6), the amino acid sequence of the compound,
Figure QLYQS_57
、/>
Figure QLYQS_60
corresponding +.>
Figure QLYQS_63
、/>
Figure QLYQS_58
Multipath calculation value of any observation epoch of frequency, the unit is meter; />
Figure QLYQS_59
And
Figure QLYQS_62
carrier frequency points corresponding to the first frequency and the second frequency respectively are in megahertz; />
Figure QLYQS_65
And->
Figure QLYQS_56
Pseudo-range observed quantity corresponding to the epoch of the first frequency and pseudo-range observed quantity corresponding to the epoch of the second frequency are respectively measured in meters; />
Figure QLYQS_61
And->
Figure QLYQS_64
The unit is meter for the observed quantity of the carrier phase of the epoch corresponding to the first frequency and the observed quantity of the carrier phase of the epoch corresponding to the second frequency;
the said
Figure QLYQS_66
Calculated according to the formula (7):
Figure QLYQS_67
(7)
in the formula (7), the amino acid sequence of the compound,
Figure QLYQS_68
the total number of the calendar elements is observed in the observation period; />
Figure QLYQS_69
For epoch number,/->
Figure QLYQS_70
10. The mass data cleaning treatment and quality check evaluation method for the CORS station network according to claim 4, which is characterized in that: the pseudo-range noise analysis specifically comprises:
calculating pseudo-range noise average values of all observation satellites of any frequency of any Beidou/GNSS system in the observation period according to a formula (8):
Figure QLYQS_71
(8)
in the formula (8), the amino acid sequence of the compound,
Figure QLYQS_72
in the observation period, the unit of any frequency pseudo-range noise average value of any Beidou/GNSS system is meter (m); />
Figure QLYQS_73
In the observation period, the total number of satellites observed by any frequency of any Beidou/GNSS system; />
Figure QLYQS_74
In order to observe the satellite serial number,
Figure QLYQS_75
;/>
Figure QLYQS_76
in order to observe the period, any frequency of any Beidou/GNSS system is +.>
Figure QLYQS_77
Pseudo-range noise RMS for a satellite is in meters.
11. The mass data cleaning and quality check and assessment method for a CORS station network of claim 10, further comprising:
Figure QLYQS_78
calculated according to the formula (9):
Figure QLYQS_79
(9)
in the formula (9):
Figure QLYQS_80
in order to observe the period, any frequency of any Beidou/GNSS system is +.>
Figure QLYQS_81
The total number of observation epochs of the satellites; />
Figure QLYQS_82
For epoch number,/->
Figure QLYQS_83
;/>
Figure QLYQS_84
For any frequency of any Beidou/GNSS system +.>
Figure QLYQS_85
Satellite, in epoch->
Figure QLYQS_86
Pseudo-range noise estimation of time, the unit is meter;
Figure QLYQS_87
calculated according to the formula (10):
Figure QLYQS_88
(10)
in the formula (10), the amino acid sequence of the compound,
Figure QLYQS_89
for any frequency of any Beidou/GNSS system +.>
Figure QLYQS_90
Satellite, in epoch->
Figure QLYQS_91
The unit of the pseudo-range observation value of the moment is meter; />
Figure QLYQS_92
For any frequency of any Beidou/GNSS system +.>
Figure QLYQS_93
Satellite, in epoch->
Figure QLYQS_94
And the unit of the pseudo-range quadratic polynomial fitting value of the moment is meter.
12. The mass data cleaning and inspection and assessment method for a CORS station network of claim 11, further comprising:
performing pseudo-range quadratic polynomial fitting according to a formula (11); after fitting calculation to obtain a quadratic polynomial coefficient, obtaining a pseudo-range quadratic polynomial fitting value of each epoch
Figure QLYQS_95
The method comprises the steps of carrying out a first treatment on the surface of the In the observation period, starting from the initial epoch, each 120 epochs are a fitting window, and the fitting windows are not overlapped; fitting according to the number of the remaining calendar elements when the number of the remaining calendar elements is not less than 3 calendar element observation data near the end of the observation period; discarding the data when the number of the remaining epoch is less than 3 epoch observation data; when the observation data has an interruption phenomenon, the segmentation is respectively processed:
Figure QLYQS_96
(11)
in the formula (11), the amino acid sequence of the compound,
Figure QLYQS_97
fitting a function for a pseudo-range quadratic polynomial; />
Figure QLYQS_98
To fit epoch number within window, quadratic polynomial argument, ++>
Figure QLYQS_99
;/>
Figure QLYQS_100
Is a quadratic polynomial coefficient; />
Figure QLYQS_101
A first order term coefficient which is a quadratic polynomial;
Figure QLYQS_102
is a quadratic polynomial constant term.
13. The mass data cleaning treatment and quality check evaluation method for a CORS station network as set forth in claim 4, wherein the carrier phase noise analysis specifically includes:
calculating the carrier-phase noise average value of all observation satellites at any frequency of any Beidou/GNSS system in the observation period according to a formula (12):
Figure QLYQS_103
(12)
in the formula (12), the amino acid sequence of the compound,
Figure QLYQS_104
in the observation period, the average value of phase noise of any frequency carrier of any Beidou/GNSS system is given in units of weeks; />
Figure QLYQS_105
In the observation period, the total number of satellites observed by any frequency of any Beidou/GNSS system; />
Figure QLYQS_106
In order to observe the satellite serial number,
Figure QLYQS_107
;/>
Figure QLYQS_108
in order to observe the period, any frequency of any Beidou/GNSS system is +.>
Figure QLYQS_109
Carrier phase noise of a satellite is in units of weeks.
14. The mass data cleaning and inspection and assessment method for a CORS station network according to claim 13, further comprising:
calculated according to the formula (13)
Figure QLYQS_110
Figure QLYQS_111
(13)
In the formula (13):
Figure QLYQS_112
in order to observe the period, any frequency of any Beidou/GNSS system is +.>
Figure QLYQS_113
The total number of observation epochs of the satellites; />
Figure QLYQS_114
For epoch number,/->
Figure QLYQS_115
;/>
Figure QLYQS_116
For any frequency of any Beidou/GNSS system +.>
Figure QLYQS_117
Noise estimation of three differences of carrier phase observation values of adjacent epochs of satellites is given in units of weeks;
calculated according to the formula (14)
Figure QLYQS_118
Figure QLYQS_119
(14)
In the formula (14), the amino acid sequence of the compound,
Figure QLYQS_120
for epoch->
Figure QLYQS_121
Moment, any satellite->
Figure QLYQS_122
The carrier phase observations at any frequency are in units of weeks.
15. The mass data cleaning treatment and quality check evaluation method for a CORS station network as set forth in claim 4, wherein the carrier-to-noise ratio analysis specifically includes:
calculating the average value of the carrier-to-noise ratio statistics of all satellites at any frequency of any Beidou/GNSS system according to a formula (15):
Figure QLYQS_123
(15)
in the formula (15), the amino acid sequence of the compound,
Figure QLYQS_124
the average value of the carrier-to-noise ratio statistics values of all satellites at any frequency of any Beidou/GNSS system is expressed in decibel hertz; />
Figure QLYQS_125
Observing the total number of satellites in an observation period; />
Figure QLYQS_126
For observing satellite serial number>
Figure QLYQS_127
;/>
Figure QLYQS_128
In order to be within the observation period, any frequency of any Beidou/GNSS system is +.>
Figure QLYQS_129
The average value of the carrier-to-noise ratio of the satellite is expressed in decibel hertz.
16. The mass data cleaning and inspection and assessment method for a CORS station network according to claim 15, further comprising:
calculating the said according to formula (16)
Figure QLYQS_130
Figure QLYQS_131
(16)
In formula (16):
Figure QLYQS_132
the total number of the calendar elements is observed in the observation period; />
Figure QLYQS_133
For epoch number,/->
Figure QLYQS_134
;/>
Figure QLYQS_135
For epoch->
Figure QLYQS_136
Moment, satellite->
Figure QLYQS_137
The observed amount of carrier-to-noise ratio at any frequency is in dB hertz.
17. The mass data cleaning treatment and quality check evaluation method for the CORS station network as claimed in claim 1, wherein the mass data cleaning treatment and quality check evaluation method is characterized by comprising the following steps of: s4 specifically comprises the following steps:
the quality evaluation and grading evaluation of the data assets of the monitoring station network are used for carrying out comprehensive quality evaluation/grading on the data assets formed by the data of the monitoring station after treatment and quality verification; the specific method comprises the following steps: taking the indexes of the data integrity rate, the data effective rate, the pseudo-range multipath, the cycle slip ratio and the signal to noise ratio as the indexes of a comprehensive evaluation model of the TOPSIS of the sequencing method approaching ideal points; the comprehensive evaluation process is to synthesize a plurality of evaluation indexes into an integrated comprehensive evaluation index through a mathematical model, combine the special navigation time-frequency service types, determine weight coefficients by forward and dimensionless methods and entropy methods on the indexes, and determine a final comprehensive evaluation model so as to obtain a final evaluation result; and (3) carrying out grade assessment on the observation data file by utilizing the comprehensive evaluation result of the observation data and referring to a set threshold value, wherein the grade assessment is respectively four grades of excellent, good, qualified and unavailable.
CN202310648390.2A 2023-06-02 2023-06-02 Method for cleaning, treating and quality checking and evaluating mass data of CORS station network Pending CN116383191A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310648390.2A CN116383191A (en) 2023-06-02 2023-06-02 Method for cleaning, treating and quality checking and evaluating mass data of CORS station network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310648390.2A CN116383191A (en) 2023-06-02 2023-06-02 Method for cleaning, treating and quality checking and evaluating mass data of CORS station network

Publications (1)

Publication Number Publication Date
CN116383191A true CN116383191A (en) 2023-07-04

Family

ID=86971434

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310648390.2A Pending CN116383191A (en) 2023-06-02 2023-06-02 Method for cleaning, treating and quality checking and evaluating mass data of CORS station network

Country Status (1)

Country Link
CN (1) CN116383191A (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103901440A (en) * 2014-03-14 2014-07-02 中国测绘科学研究院 GNSS data signal quality monitor method
CN110967719A (en) * 2019-12-25 2020-04-07 北斗天地股份有限公司 Processing method and device of multi-navigation system based on Beidou navigation
WO2021237804A1 (en) * 2020-05-29 2021-12-02 湖南联智科技股份有限公司 Infrastructure structure deformation monitoring method based on beidou high-precision positioning

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103901440A (en) * 2014-03-14 2014-07-02 中国测绘科学研究院 GNSS data signal quality monitor method
CN110967719A (en) * 2019-12-25 2020-04-07 北斗天地股份有限公司 Processing method and device of multi-navigation system based on Beidou navigation
WO2021237804A1 (en) * 2020-05-29 2021-12-02 湖南联智科技股份有限公司 Infrastructure structure deformation monitoring method based on beidou high-precision positioning

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"BD 420022-2019,北斗/全球卫星导航系统(GNSS)测量型接收机观测数据质量评估方法", 中国第二代卫星导航系统重大专项标准 *
郭亮亮 等: "GNSS测量型接收机观测数据质量评估", 第八届中国卫星导航学术年会论文集——S08测试评估技术 *
郭亮亮: "GNSS数据质量评估软件研制与应用", 中国优秀硕士学位论文全文数据库 基础科学辑, no. 06 *

Similar Documents

Publication Publication Date Title
CN109359270B (en) Threshold model establishing method for integrity risk monitoring of Beidou foundation enhancement system
TWI425238B (en) Method of position determination in a global navigation satellite system receiver
CN104965207A (en) Method for acquiring area troposphere zenith delay
CN101950024B (en) Code carrier consistency detection method applied to local area augmentation system
CN115993623B (en) Adaptive star selection method, device, equipment and readable storage medium
CN111435167B (en) Receiver data quality analysis method and device based on Beidou III
KR101152399B1 (en) DGNSS Reference Station and method of estimating a User Differential Range Error thereof
CN116931026B (en) Abnormality determination method for satellite navigation signals
CN112213742A (en) Signal quality monitoring method for satellite navigation system
CN114485655B (en) GNSS/INS combined navigation data quality control method
CN112130177A (en) Foundation reinforcement system integrity monitoring method based on stable distribution
CN113325446A (en) Multi-mode common-frequency GNSS carrier phase time transfer method and system
Xiao et al. Data quality check and visual analysis of CORS station based on ANUBIS software
CN114384565A (en) Dynamic positioning coordinate sequence abnormal value identification method based on VMD iterative decomposition
CN111007541B (en) Simulation performance evaluation method for satellite navigation foundation enhancement system
CN116383191A (en) Method for cleaning, treating and quality checking and evaluating mass data of CORS station network
CN116719073B (en) GNSS (Global navigation satellite System) solution domain-oriented coarse difference detection and rejection method
KR100884611B1 (en) Apparatus of computing the pseudorange measurement noise of reference station receiver for GNSS augmentation systems and method thereof
CN112880633A (en) Sea surface height measuring method based on Berger algorithm
CN113077121A (en) Signal quality abnormity reason positioning method based on multi-dimensional attribute group decision
CN114397679A (en) Observation data processing method
CN113532588B (en) Water level acquisition method, device, equipment and storage medium
CN115267845A (en) Method for acquiring difference code deviation, computer equipment and readable storage medium
CN106405580B (en) A kind of GNSS continuity appraisal procedure
CN111239779B (en) Blind-spot-free GNSS tri-frequency combined cycle slip detection and repair 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