Content of the invention
In view of the shortcomings of the prior art, the invention provides a kind of distributed GNSS real-time data processing method and be
System, the resolving state before the method can make service end can continue after restarting to restart is run it is ensured that the real-time of difference product and one
Cause property, thus ensure to provide the user with the positioning service of continuous precise.
For solving above-mentioned technical problem, the present invention adopts the following technical scheme that:
First, distributed GNSS real-time data processing method, including:
Step 1, carries out Detection of Cycle-slip in receiver end to the original observed data receiving, and obtains cycle slip temporal information, week
Jumping temporal information is the time that each satellite occurs cycle slip for the last time;
Step 2, sends to server end after receiver end is by cycle slip temporal information and original observed data coding;
Step 3, is decoded to the coded data receiving in server end, obtains observation data and cycle slip temporal information;
Step 4, server end by compare current epoch decoding obtain cycle slip temporal information and on one epoch decoding
Whether the cycle slip temporal information obtaining is identical, judges under current epoch, whether satellite occurs cycle slip;
Step 5, is resolved according to the observation data that cycle slip judged result and decoding obtain in server end, and in service
Device end stores current epoch and resolves the filter status parameter information obtaining;
Step 6, if resolve interrupting, after server end is restarted, the filter status Parameter reconstruction using storage interrupts epoch
Under wave filter, using recover wave filter continue resolve interrupt before resolving state.
In receiver end, Detection of Cycle-slip is carried out to the original observed data receiving using TurboEdit method in step 1.
In step 2, cycle slip temporal information and original observed data are compiled using RTCM coded system in receiver end
Code;Meanwhile, in step 3, it is decoded obtaining observation data and cycle slip temporal information using RTCM decoding process in server end.
Step 4 is specially:
The second cycle slip time of relatively the first cycle slip temporal information server end storage that current epoch decoding obtains believes
Breath, the initial value of the second cycle slip temporal information is the cycle slip temporal information that initial epoch decoding obtains;If the first cycle slip time letter
Breath is different with the second cycle slip temporal information, then respective satellite occurs cycle slip, using first cycle slip time information renovation the second cycle slip
Temporal information;Conversely, there is not cycle slip in respective satellite.
Preferably, after step 5 is resolved, using the observation data of residual error editing and processing decoding acquisition, rejecting and not visiting
Measure the Satellite Observations of cycle slip.
Being resolved according to the observation data that cycle slip judged result and decoding obtain in server end described in step 5,
It is specially:
If satellite does not occur cycle slip, directly resolved using the observation data that decoding obtains;
If satellite occurs cycle slip, resolve the fuzziness parameter of initialization cycle slip satellite by pseudorange, then using decoding
The observation data obtaining is resolved.
Resolved using Kalman filtering method in step 5.
2nd, distributed GNSS real-time data processing system, including:
Detection of Cycle-slip module, for carrying out Detection of Cycle-slip in receiver end to the original observed data receiving, obtains cycle slip
Temporal information, cycle slip temporal information is the time that each satellite occurs cycle slip for the last time;
Coding and transport module, for receiver end by cycle slip temporal information and original observed data coding after send to
Server end;
Decoder module, for being decoded to the coded data receiving in server end, when obtaining observation data and cycle slip
Between information;
Cycle slip judge module, for server end by compare current epoch decoding obtain cycle slip temporal information and on
Whether the cycle slip temporal information that the decoding of one epoch obtains is identical, judges under current epoch, whether satellite occurs cycle slip;
Resolve module, for being resolved according to the observation data that cycle slip judged result and decoding obtain in server end,
And store, in server end, the filter status parameter information that current epoch resolves acquisition;
Wave filter recovery module, if interrupting for resolving, after server end is restarted, using the filter status parameter of storage
Recover to interrupt the wave filter under epoch, continue to resolve the resolving state before interrupting using the wave filter recovering.
Compared to the prior art, the present invention has following features and beneficial effect:
1st, Detection of Cycle-slip module is moved on to receiver end, overcome server end Detection of Cycle-slip when data sampling rate is higher
Effect is poor and cannot be carried out the defect of Detection of Cycle-slip after interrupting restarting;And different samplings can be broadcast according to server end demand
The observation data of rate.
2nd, server end passes through real-time to preserve filter status parameter information so that having no progeny in server end and can utilizing
Before disconnected, the state parameter information of storage recovers wave filter, need not reinitialize the resolving state fortune before can continuing at once to interrupt
OK, it is to avoid longer convergence time, real-time and the reliability of service system are substantially increased, thus can the company of providing the user
Continuous, consistent precision differential product is it is ensured that hi-Fix service.
Specific embodiment
Further illustrate technical solution of the present invention below in conjunction with the drawings and specific embodiments.
" receiver " i.e. appearing below " reference receiver ", " server " i.e. " receiver resolving server ".
The present invention includes receiver end data prediction server end post processing two parts.When being embodied as, receiver
First pass through Detection of Cycle-slip module and the original observed data receiving is carried out with Detection of Cycle-slip acquisition cycle slip temporal information, then by week
Jump temporal information and original observed data to send to server end by RTCM coding.Server end first passes through RTCM decoding and obtains
Obtain observation data and the cycle slip temporal information of each satellite;Then, judged by cycle slip and filter estimation and carry out precision differential product
Resolve, resolve and preserve the state parameter information of current epoch respective filter after terminating to storage disk.Service end solver
Once interrupting restarting, recover wave filter so that wave filter can continue immediately by reading the state parameter information interrupting front preservation
State before continuous interruption is run.
First, receiver end data prediction.
The Detection of Cycle-slip module of real-time resolving system is located at server end at present.For low sampled data, difference between epoch
Residual error is larger afterwards, and Detection of Cycle-slip reliability is relatively low;Once running into network blockage or server end interrupting leading to adjacent epoch
Time interval is larger, and data preprocessing module need to reinitialize.
Detection of Cycle-slip module is built into receiver end by the present invention:
1) the commonly used dual-frequency receiver of base station and sample rate higher (not less than 0.2Hz), through epoch between each after difference
Item error concealment is all more clean, highly beneficial to Detection of Cycle-slip.The Cycle Slips Detection commonly used at present has TurboEdit method,
The method combines MW combination and the pluses and minuses of GF combination, can obtain preferable Detection of Cycle-slip effect.
(1) MW combined detection cycle slip:
MW combination observation lMWFor:
In formula (1),For wide lane wavelength, f1And f2It is carrier wave L respectively1And L2Frequency, Φ and P be respectively with
Rice is carrier wave and the Pseudo-range Observations of unit.
MW combination eliminates the shadow of the error relevant with geometric distance, tropospheric delay error and ionosphere delay error
Ring, only remaining width lane ambiguity, in the case of not having cycle slip, wide lane ambiguity is constant, can be used to accordingly detect cycle slip;But MW
Combination introduce Pseudo-range Observations, its noise larger it is difficult to detect little cycle slip, simultaneously MW combination also cannot detect L1 and L2 frequency
The observation of identical cycle slip is occurred on rate.
(2) Geometry-Free combined detection cycle slip:
Free-geometry combination can be expressed as:
In formula (2), ρ is the geometric distance between receiver for the satellite, and A is the ionosphere factor, N1And N2For what L2 frequency of L1
Fuzziness parameter, orderThen (2) can be reduced to:
lGF=Δ dIon+ (λ1N1-λ2N2) (3)
GF combination eliminates the relevant error of geometric distance and the impact of the error term relevant with frequency, only remaining remnants
Ionosphere delay impact and fuzziness parameter, because reference receiver sample rate is higher, between epoch, Ionospheric variability is more flat
Slow, therefore GF combination observation is slowly varying and more smooth curve.
2) through Detection of Cycle-slip, the time that satellite occurs cycle slip for the last time, i.e. cycle slip temporal information are obtained, receiver end will
The cycle slip temporal information of each satellite and original observed data encoded along after, send to server end, its coded format is shown in Table 1.
Table 1 coded format
Data content |
Bit/s |
Explanation |
Original observed data |
— |
Follow corresponding RTCM agreement |
The cycle slip time |
20 |
Zhou Miao |
3) according to actual needs coded data is broadcast to service end using different sample rates.
2nd, server end post processing
Server end obtains observing data and cycle slip temporal information by decoding, no longer carries out Detection of Cycle-slip, only carries out essence
The real-time resolving of close difference product, can be using Kalman filter or other routine filtering methods of estimation.Resolve knot in each epoch
Preserve the state parameter information of this epoch respective filter after bundle, after service end meets abnormal restarting, deposit before being interrupted by reading
The filter status parameter information of storage builds wave filter so that server end can carry out real-time resolving without convergence process.
Server end post processing of the present invention mainly includes:
S1, cycle slip judge:
A () server end passes through the cycle slip temporal information that decoding obtains each every satellite of base station;
B () compares the second cycle slip time of the first cycle slip temporal information server end record that current epoch decoding obtains
Information judges whether satellite occurs cycle slip:If hop-information was different with second week hop-information in first week, illustrate that corresponding satellite is sent out
Raw cycle slip;Otherwise then there is not cycle slip.Second cycle slip temporal information initial value of server end record is that initial epoch decoding obtains
The cycle slip temporal information obtaining.If it is determined that generation cycle slip, using the second cycle slip of the first cycle slip information updating server end record
Information, the cycle slip for next epoch judges.
When in this step, the comparison of cycle slip temporal information is by the first cycle slip temporal information of same satellite and the second cycle slip
Between information compare, to judge whether this satellite occurs cycle slip.
C () receiver end data prediction there may be the cycle slip that very small amount not yet detects, can pass through residual error Editorial Services
Reason observation data, after the cycle slip not detected with rejecting, is resolved again and eliminates remaining cycle slip impact, renewal simultaneously is recorded
All hop-informations.
The ultimate principle of residual error editor is:If certain satellite observation has the cycle slip not detected or rough error, it is fixed
Position result residual error should be big compared with other satellites.The purpose of residual error editor is that algorithm positions the larger satellite of residual error and rejects this satellite.
S2, resolving obtain filter status parameter:
Obtain observation data and the cycle slip temporal information of each satellite of current epoch by decoding;Satellite is judged according to step S1
Whether there is cycle slip, if there is cycle slip, resolving the fuzziness parameter of initialization cycle slip satellite by pseudorange;Then, using solution
The observation data that code obtains carries out real-time resolving, thus realizing positioning.If there is not cycle slip, directly using the sight of decoding acquisition
Survey data and carry out real-time resolving.
S3, the storage of filter status parameter information:
After resolving terminates, under storage current epoch, the state parameter information of wave filter, is shown in Table 2:
Table 2 filter status parameter information
Furthermore it is also possible to store following data according to the needs of the real-time resolving algorithm being adopted in server end:
(1) the cycle slip temporal information of each satellite that decoding obtains;
(2) the observation data of nearest n epoch, n carries out value with the computation being adopted according to actual needs;
(3) current epoch temporal information;
(4) other information, such as:In data delay information in real-time resolving, adaptive-filtering resolving, each base station is many
Routing information, data quality information of each satellite etc..
To taking Kalman filter method as a example illustrate the solution process of server end below, step is as follows:
(1) the filter status parameter being resolved according to a upper epoch and its variance-covariance matrix information are carried out
Kalman filter is predicted, if initial epoch, then needs given filter status initial parameter values.
(2) server end decoding obtains observation data and the cycle slip temporal information of each satellite of each base station, during according to cycle slip
Between information judge respective satellite the need of reinitializing unknown parameter:
If 2.1 current epoch are initial epoch Tk(k=0), record the cycle slip temporal information of each satellite of base station, continue step
Suddenly (3).
2.2, by comparing first week hop-information and the second week hop-information recording that current epoch decoding obtains, judge to defend
Whether star there is cycle slip:If hop-information was different with second week hop-information in first week, illustrate that satellite occurs cycle slip;Conversely, not sending out
Raw cycle slip.In the event of cycle slip, using the second week hop-information of the first cycle slip information updating record, for the cycle slip of next epoch
Judge.
2.3 pass through the very small amount cycle slip that residual error editor rejecting does not detect, and the cycle slip temporal information of more new record.
If not having satellite that cycle slip occurs, it is directly entered step (3);If there is satellite cycle slip occurs, reinitialize
Step (3) is entered after the relevant parameter of the satellite of this generation cycle slip.
(3) the observation data real-time resolving according to decoding the filter status parameter information under obtaining current epoch, adopt
Update wave filter with filter status parameter information and preserve.
(4) repeat step (1)~(3) can carry out real-time resolving.
Server end post-processing function mainly uses observation data and the cycle slip temporal information that decoding obtains, and carries out in real time
Resolve, and each epoch resolve terminate after preserve the state parameter information of wave filter under this epoch;When server end runs into exception
After interrupting restarting, by reading the filter status parameter information reconfigurable filter interrupting front storage so that server end need not
Convergence process can carry out real-time resolving.
S4, server end interrupt restarting postfilter structure:
If server end collapse or aborted, the filter status parameter information restarting the storage of rear reading service device end is extensive
Wave filter before multiple interruption, and repeat step S2 can continue the resolving state operation before interruption.
Assume server end in epoch TkCollapse or aborted, in epoch Tk+mRestart the storage of rear reading service device end
Epoch TkLower filter status parameter information, for recovering epoch TkUnder wave filter, then repeat step S2 can make server
Resolving state before end continues directly to without convergence time interrupt is run.