CN105549055A - Virtual observation data generation method and device - Google Patents

Virtual observation data generation method and device Download PDF

Info

Publication number
CN105549055A
CN105549055A CN201510822773.2A CN201510822773A CN105549055A CN 105549055 A CN105549055 A CN 105549055A CN 201510822773 A CN201510822773 A CN 201510822773A CN 105549055 A CN105549055 A CN 105549055A
Authority
CN
China
Prior art keywords
network element
base station
baseline
satellite
difference
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
CN201510822773.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.)
Space Star Technology Co Ltd
Original Assignee
Space Star 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 Space Star Technology Co Ltd filed Critical Space Star Technology Co Ltd
Priority to CN201510822773.2A priority Critical patent/CN105549055A/en
Publication of CN105549055A publication Critical patent/CN105549055A/en
Pending legal-status Critical Current

Links

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/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

Abstract

The invention discloses a virtual observation data generation method and device. The method comprises the steps that a base station network element in which a terminal is positioned is determined; a main reference base station, network element global common-view satellites and a reference satellite are determined in the base station network element according to the determined base station network element; double differential atmospheric parameters of a network element baseline in which the main reference base station is positioned are calculated so that the double differential atmospheric parameters of a virtual baseline are obtained through interpolation; and finally the double differential atmospheric parameters are converted into single differential atmospheric parameters through a single differential model, and carrier phase and pseudo-range virtual observation values corresponding to all systems and all frequency points are generated. With application of the method, the problem of redundant matrix transform caused by inconsistence generated by independent selection of all the baseline common-view satellites and the reference satellite in the network element can be solved. Besides, virtual reference base station observation data are constructed by adopting the single differential model equivalent to a double differential model so that result precision is guaranteed and model complexity is reduced.

Description

A kind of generation method of virtual observation data and device
Technical field
The application relates to communication technical field, particularly relates to a kind of generation method and device of virtual observation data.
Background technology
CORS system is (English: ContinuouslyOperatingReferenceStations, be called for short: CORS) conduct is to utilizing GNSS to position one of effective ground enhancements, its application progressively expands, and broadcasting enhancing data in CORS location based on virtual reference station technology to terminal is implementations that this application is comparatively popular at present.
Virtual reference station technology is that CORS data processing centre (DPC) determines its place network element by self rough coordinates that receiving terminal sends, the calculating such as Baselines, error interpolation, dummy observation generation are carried out at network element internal, generate one group of carrier phase and pseudorange observation data corresponding to this position and outwards broadcast, terminal carries out Differential positioning with it after receiving can obtain self high precision position data.
Need in the process of virtual reference base station sight data genaration to resolve the baseline in the network element of terminal place, conventional treatment scheme determines reference star depending on satellite altogether for choosing base station, baseline two ends, set up two difference observation equation based on main website and reference star to resolve respectively every bar baseline, estimate two difference atmospheric parameter last generating virtual reference station observed reading.But, in network element each bar baseline altogether depending on satellite and reference satellite to choose the possibility of result inconsistent, when virtual observed value can be caused to generate needs altogether depending on satellite repeat select, the conversion of reference satellite and correlation matrix operation, increase calculated amount; Virtual reference station observation data generation model conventional in addition obtains based on two poor observation equation, and be proven this model effectively feasible, but its mathematical model desired parameters is many, modeling is more complicated, is unfavorable for programming realization.
Summary of the invention
Embodiments provide a kind of generation method and device of virtual observation data, comprise in order to the problem solved:
1. in network element each bar baseline altogether depending on satellite and reference satellite choose needs when the inconsistent virtual observed value caused generates altogether depending on satellite repeat select, the conversion of reference satellite and correlation matrix operation, increase calculated amount.
2. conventional virtual reference station observation data generation model obtains based on two poor observation equation, and its mathematical model desired parameters is many, and modeling is more complicated, is unfavorable for programming realization.
Its concrete technical scheme is as follows:
A generation method for virtual observation data, described method comprises:
Determine main reference station in base station network element residing for terminal and network element;
According to the base station network element determined, in the network element of described base station, determine that the network element overall situation looks satellite and reference satellite altogether;
Obtain the two difference atmospheric parameters on the network element baseline of described main reference base station place;
According to the two difference atmospheric parameters on network element baseline, obtain the two difference atmospheric parameters on virtual baseline;
Described two difference atmospheric parameter is converted to single poor atmospheric parameter, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation.
Optionally, determine main reference station in base station network element residing for terminal and network element, comprising:
According to the rough coordinates position of terminal, the distance of each base station during acquisition terminal and CORS system CORS net, and filter out the base station nearest with described terminal;
Using the network element residing for described base station as alternative network element;
Obtain the distance of rough coordinates position baseline corresponding to nearest base station to alternative network element of terminal, select the baseline that distance terminal is nearest;
The network element that the base station at the nearest base station of distance terminal and described nearest baseline two ends is formed as the described base station network element at terminal place, using described nearest base station main reference station in network element.
Optionally, according to the base station network element determined, in the network element of described base station, determine that the described network element overall situation looks satellite and reference star altogether, comprising:
Determine in described network element each bar baseline corresponding look satellite altogether;
Each base station corresponding altogether depending on satellite in filter out the overall situation and look satellite altogether;
The satellite with maximum height angle is filtered out as described network element overall situation reference satellite in the described overall situation is altogether depending on satellite.
Optionally, according to the two difference atmospheric parameters on network element baseline, obtain the two difference atmospheric parameters on virtual baseline, comprising:
Calculate the two difference blur level on two baselines at main reference base station place by Kalman filtering and be fixed, then estimating baseline two difference ionosphere delay and two poor tropospheric delay;
Mathematical interpolation is carried out to the two difference ionosphere delays on described two baselines and two poor tropospheric delay, obtains the two difference ionosphere delays on described virtual baseline and two poor tropospheric delay.
Optionally, described two difference atmospheric parameter is converted to single poor atmospheric parameter, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation, be specially:
Carry out mathematical modeling by single differential mode type, the two difference ionosphere delays on the described virtual baseline obtained and two poor tropospheric delay are converted to single poor delay, and generate carrier phase corresponding to each frequency of each system and pseudorange dummy observation.
A generating apparatus for virtual observation data, comprising:
First determination module, for determining main reference station in base station network element residing for terminal and network element; ;
Second determination module, for according to the base station network element determined, determines network element overall situation reference satellite in the network element of described base station;
Acquisition module, for obtaining the two difference atmospheric parameters on the network element baseline of described main reference base station place;
Processing module, for obtaining the two difference atmospheric parameters on virtual baseline according to the two difference atmospheric parameters on described network element baseline;
Generation module, for the two difference atmospheric parameters on described virtual baseline are converted to single poor atmospheric parameter, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation.Optionally, described first determination module, specifically for the rough coordinates position according to terminal, the distance of each base station during acquisition terminal and CORS system CORS net, and filter out the base station nearest with described terminal; Using the network element residing for described base station as alternative network element; Obtain the distance of rough coordinates position baseline corresponding to nearest base station to alternative network element of terminal, select the baseline that distance terminal is nearest; The network element that the base station at the nearest base station of distance terminal and described nearest baseline two ends is formed as the described base station network element at terminal place, and using described nearest base station main reference station in network element.
Optionally, described second determination module, specifically for determine in the network element that described main reference base station is corresponding each base station corresponding look satellite altogether; Each base station corresponding altogether depending on satellite in filter out the overall situation and look satellite altogether; The satellite with maximum height angle is filtered out as described network element overall situation reference satellite in the described overall situation is altogether depending on satellite.
Optionally, described acquisition module, specifically for being calculated the two difference blur leveles on two baselines at main reference base station place by Kalman filtering, calculates two difference ionosphere delays and two poor tropospheric delay after being fixed; Mathematical interpolation is carried out to the two difference ionosphere delays on described two baselines and two poor tropospheric delay, obtains the two difference ionosphere delays on described virtual baseline and two poor tropospheric delay.
Optionally, described generation module, specifically for carrying out mathematical modeling by single differential mode type, the two difference ionosphere delays on the described virtual baseline obtained and two poor tropospheric delay being converted to single poor delay, and generating carrier phase corresponding to each frequency of each system and pseudorange dummy observation.
The invention provides a kind of generation method and device of virtual observation data, the method comprises: determine base station network element residing for terminal, according to the base station network element determined, main reference base station and network element overall situation reference satellite is determined in the network element of base station, calculate the two difference atmospheric parameters on the network element baseline of main reference base station place, and then interpolation obtains the two difference atmospheric parameters on virtual baseline, two poor atmospheric parameter is converted to single poor atmospheric parameter finally by single differential mode type, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation.Solve each baseline in network element by above-mentioned side and independently choose the problem of the inconsistent unnecessary matrixing brought of generation altogether depending on satellite and reference satellite.And adopt the single differential mode type constructing virtual reference base station observation data with two difference model equivalency, while guarantee result precision, reduce model complexity.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of the generation method of a kind of virtual observation data in the embodiment of the present invention;
Fig. 2 is the structural representation that the networking of embodiment of the present invention intermediate cam generates net form;
Fig. 3 is the structural representation of the generating apparatus of a kind of virtual observation data in the embodiment of the present invention.
Embodiment
In order to solve in prior art each bar baseline in network element altogether depending on satellite and reference satellite to choose the possibility of result inconsistent, when virtual observed value can be caused to generate needs altogether depending on satellite repeat select, the conversion of reference satellite and correlation matrix operation, increase the problem of calculated amount, simultaneously many for solving the mathematical model desired parameters that the traditional double differential mode type that uses in terminal dummy observation generative process brings, modeling is more complicated, be unfavorable for the problem of programming realization, embodiments provide a kind of generation method and device of virtual observation data, the method comprises: determine main reference station in base station network element residing for terminal and network element, determine that the network element overall situation looks satellite and reference satellite altogether according to the base station network element determined, obtain the two difference atmospheric parameters on the network element baseline of main reference base station place, according to the two difference atmospheric parameters on network element baseline, obtain the two difference atmospheric parameters on virtual baseline, by single differential mode type, two poor atmospheric parameter is converted to single poor atmospheric parameter, and generate each system and carrier phase corresponding to each frequency and pseudorange dummy observation.Solve each base station in network element by above-mentioned side and independently choose the problem of the inconsistent unnecessary matrixing brought of generation altogether depending on satellite and main reference satellite.And adopt single differential mode type constructing virtual reference base station observation data of two difference model equivalency, while guarantee result precision, reduce model complexity.
Below by accompanying drawing and specific embodiment, technical solution of the present invention is described in detail, be to be understood that, the explanation of concrete technical characteristic in the embodiment of the present invention and embodiment just to technical solution of the present invention, instead of limit, when not conflicting, the concrete technical characteristic in the embodiment of the present invention and embodiment can combine mutually.
Be illustrated in figure 1 the process flow diagram of the generation method of a kind of virtual observation data in the embodiment of the present invention, the method comprises:
S101, determines base station network element residing for terminal and the main reference station of network element;
Specifically, in embodiments of the present invention, first get the rough coordinates position of terminal, according to the rough coordinates position of terminal, the distance of each base station during acquisition terminal and CORS system CORS net, and filter out the base station nearest with terminal.
Using the network element residing for base station as alternative network element, obtain the distance of rough coordinates position baseline corresponding to nearest base station to alternative network element of terminal, select the baseline that distance terminal is nearest, the network element that the base station at the nearest base station of distance terminal and nearest baseline two ends is formed as the base station network element at terminal place, using described nearest base station as the main reference station of network element.
S102, according to the base station network element determined, determines that in the network element of base station the network element overall situation looks satellite and reference satellite altogether;
In S101, determine main reference base station and the base station network element residing for main reference base station, then determine in the network element that main reference base station is corresponding each bar baseline corresponding look satellite altogether, each baseline corresponding altogether depending on satellite in filter out the overall situation and look satellite altogether, in the network element determined, such as contain A, B, E tri-base stations, base station A, B's is a, b, c depending on satellite altogether, i, j, k, base station A, E's is l, m, n depending on satellite altogether, i, j, k, then the overall situation looks satellite is altogether just i, j, k.
The satellite with maximum height angle is filtered out as network element overall situation reference satellite in the overall situation is altogether depending on satellite,
S103, obtains the two difference atmospheric parameters on the network element baseline of main reference base station place;
After determining overall reference satellite, Kalman filtering is used to calculate main reference station place baseline two difference ionosphere delay same epoch and two poor tropospheric delay in network element.
Two difference ionosphere delay herein and two difference tropospheric delay are as the two difference atmospheric parameters on network element baseline.
S104, according to the two difference atmospheric parameters on network element baseline, obtains the two difference atmospheric parameters on virtual baseline;
Calculate the two difference blur level on two baselines at main reference base station place by Kalman filtering and be fixed, then estimating baseline two difference ionosphere delay and two poor tropospheric delay; Mathematical interpolation is carried out to the two difference ionosphere delays on described two baselines and two poor tropospheric delay, obtains the two difference ionosphere delays on described virtual baseline and two poor tropospheric delay.
Carry out interpolation to the two difference atmospheric parameters got, ionosphere interpolation model uses LIM model, and troposphere interpolation model uses LSM model, thus obtains the two difference atmospheric parameters on virtual baseline.
S105, is converted to single poor atmospheric parameter by two poor atmospheric parameter, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation.
Carry out mathematical modeling by single differential mode type, the two difference ionosphere delays on the described virtual baseline obtained and two poor tropospheric delay are converted to single poor delay, and generate carrier phase corresponding to each frequency of each system and pseudorange dummy observation.
Below by concrete implementation, technical solution of the present invention is described in detail.
As shown in Figure 2, there are base station A, B, C, D, E totally 5 places in CORS net, according to Delaunay triangle networking generate rule net form as shown in Figure 2, end flow erect-position is in F place, its rough coordinates position is f, by the method for the invention with the dummy observation of f point for virtual reference station generation correspondence:
(1) determine that terminal place network element is ABE, in network element, main reference base station is A base station;
(2) that selects base station A, B looks satellite altogether: a, b, c, i, j, k;
(3) that selects base station A, E looks satellite altogether: l, m, n, i, j, k;
(4) select AB and AE and look public satellite in satellite altogether: i, j, k;
(5) by select in (4) altogether depending on satellite as stations all in network element ABE (A, B, E, f) look satellite altogether; Then this network element overall situation is i, j, k depending on satellite altogether;
(6) according to selecting star result in (4), selecting the reference star of the maximum satellite of elevation angle as station A, B, E, f, being set to k;
(7) use Kalman filtering to resolve baseline AB, AE and BE, to obtain on baseline same epoch two difference ionosphere delay with two difference tropospheric delay;
Ignore other errors, L1, L2 dual-frequency carrier two difference observation equation is:
L1, L2 double frequency pseudorange two difference observation equation is:
Wherein v is double difference residual error, for two poor operator, for carrier phase observed quantity, P is pseudo range observed quantity, and ρ is for defending the true distance in ground, ε ion-L1for ionosphere delay, ε tropfor tropospheric delay.Because CORS base station accurate coordinates is known, then two difference defends ground very apart from can directly obtain, and in equation, unknown parameter comprises two difference ionospheres of every bar baseline, two poor troposphere and pair difference blur level.
Baseline blur level uses " three-step approach " to resolve, by Householder orthogonal transformation for Kalman filtering provides initial value; Use LAMBDA methods to fix to the two difference float ambiguities drawn after estimating, in network element, each baseline two difference blur level is fixed correctness and is judged according to ratio value, checks further by following formula simultaneously:
Then fixing rear blur level is used to obtain tropospheric delay, the upper two difference ionosphere delay of L1, then according to the relation that L1 and L2 upper ionized layer postpones:
Obtain L2 two difference ionosphere delay.
(8) carry out interpolation to two mistake differences that (7) obtain, ionosphere interpolation model uses LIM model, and troposphere interpolation model uses LSM model, and the two mistakes obtained on virtual baseline Af are poor.
(9) the two mistake differences calculating acquisition in (8) are converted to single mistake by single differential mode type poor, and construct f place virtual observed value, concrete steps are as follows:
First make condition hypothesis, be provided with:
1) error be located in various observed quantities virtual reference base station f observing main reference satellite k is equal to the various errors of main reference satellite k with main reference base station A, non-reference looks satellite altogether except ionosphere and tropospheric delay, and other errors in observed quantity are also equal.
2) f station observes that each to look satellite altogether equal to each depending on satellite blur level and A altogether;
Secondly, single mistake that two difference ionospheres that can obtain in (8) according to assumed condition, tropospheric error are converted to equivalence is poor:
3rd, generating virtual reference station is relative to each Pseudo-range Observations looking satellite altogether.Conventional two differential mode types are:
Need expansion every in formula during calculating:
Wherein P is pseudo range observed quantity, and ρ is for defending the true distance in ground, and Δ ε is the poor summation of list of various error, because f and A two station is equal for the observational error of reference star k, then has:
Δε=0
Replace corresponding two mistake difference by single poor ionosphere delay and tropospheric delay and arrange:
The pseudo range observed quantity computing formula of non-reference satellite relative to virtual reference station can be obtained.
For reference satellite, take condition hypothesis 1 into account, the pseudo range observed quantity of virtual reference base station can generate as follows:
Non-reference can derive by with pseudo range observed quantity similar procedure depending on the generation of satellite carrier carrier phase observable altogether.Non-reference satellite virtual reference base station carrier Phase Double differential mode type:
Above formula need be launched equally during calculating:
From assumed condition Δ ε=0, above formula arranges and is:
λΦ f i = Δρ f A i - Δϵ i o n o ( f - A ) i + Δϵ t r o p ( f - A ) i + λΦ A i
The observed quantity of reference star carrier phase generates and presses following formula:
λΦ f k = Δρ f A k + λΦ A k
Visible according to above derivation, virtual reference station carrier phase and pseudo range observed quantity can by single poor model generations, and its form and the more two differential mode type of Computed-torque control are simplified.
A kind of generation method of virtual observation data in the corresponding embodiment of the present invention, a kind of generating apparatus of virtual observation data is additionally provided in the embodiment of the present invention, be illustrated in figure 3 the structural representation of the generating apparatus of a kind of virtual observation data in the embodiment of the present invention, this device comprises:
First determination module 301, for determining main reference station in base station network element residing for terminal and network element;
Second determination module 302, for according to the base station network element determined, determines that the network element overall situation looks satellite and reference satellite altogether in the network element of described base station;
Acquisition module 303, for obtaining the two difference atmospheric parameters on the network element baseline of described main reference base station place;
Processing module 304, for obtaining the two difference atmospheric parameters on virtual baseline according to the two difference atmospheric parameters on described network element baseline;
Generation module 305, for the two difference atmospheric parameters on described virtual baseline are converted to single poor atmospheric parameter, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation.
Further, in embodiments of the present invention, described first determination module 301, specifically for obtain according to the rough coordinates of terminal terminal and CORS system CORS net in the distance of each base station, and filter out the base station nearest with described terminal; Using the network element residing for described base station as alternative network element; Obtain the distance of rough coordinates position baseline corresponding to nearest base station to alternative network element of terminal, select the baseline that distance terminal is nearest; The network element that the base station at the nearest base station of distance terminal and described nearest baseline two ends is formed as the described base station network element at terminal place, and using described nearest base station as the main reference station of network element.
Further, in embodiments of the present invention, described second determination module 302, specifically for determine in the network element that described main reference base station is corresponding each base station corresponding look satellite altogether; Each base station corresponding altogether depending on satellite in filter out the overall situation and look satellite altogether; Filter out in the described overall situation is altogether depending on satellite there is maximum height angle satellite as reference satellite overall in described network element.
Further, in embodiments of the present invention, described acquisition module 303, specifically for calculating the two difference blur leveles on two baselines at main reference base station place, the two difference ionosphere delays after being fixed on clearing place baseline and two poor tropospheric delay by Kalman filtering; Mathematical interpolation is carried out to the two difference ionosphere delays on described two baselines and two poor tropospheric delay, obtains the two difference ionosphere delays on described virtual baseline and two poor tropospheric delay.
Further, in embodiments of the present invention, described generation module 304, specifically for carrying out mathematical modeling by single differential mode type, two difference ionosphere delays on the described virtual baseline obtained and two poor tropospheric delay are converted to single poor delay, and generate carrier phase corresponding to each frequency of each system and pseudorange dummy observation.
Although described the preferred embodiment of the application, one of ordinary skilled in the art once obtain the basic creative concept of cicada, then can make other change and amendment to these embodiments.So claims are intended to be interpreted as comprising preferred embodiment and falling into all changes and the amendment of the application's scope.
Obviously, those skilled in the art can carry out various change and modification to the application and not depart from the spirit and scope of the application.Like this, if these amendments of the application and modification belong within the scope of the application's claim and equivalent technologies thereof, then the application is also intended to comprise these change and modification.

Claims (10)

1. a generation method for virtual observation data, is characterized in that, described method comprises:
Determine main reference station in base station network element residing for terminal and network element;
According to the base station network element determined, in the network element of described base station, determine that the network element overall situation looks satellite and reference satellite altogether;
Obtain the two difference atmospheric parameters on the network element baseline of described main reference base station place;
According to the two difference atmospheric parameters on network element baseline, obtain the two difference atmospheric parameters on virtual baseline;
Described two difference atmospheric parameter is converted to single poor atmospheric parameter, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation.
2. the method for claim 1, is characterized in that, determines main reference station in base station network element residing for terminal and network element, comprising:
According to the rough coordinates of terminal, the distance of each base station during acquisition terminal and CORS system CORS net, and filter out the base station nearest with described terminal;
Using the network element residing for described base station as alternative network element;
Obtain the distance of rough coordinates position baseline corresponding to nearest base station to alternative network element of terminal, select the baseline that distance terminal is nearest;
The network element that the base station at the nearest base station of distance terminal and described nearest baseline two ends is formed as the described base station network element at terminal place, using described nearest base station main reference station in network element.
3. method as claimed in claim 2, is characterized in that, according to the base station network element determined, determine that the described network element overall situation looks satellite and reference satellite altogether, comprising in the network element of described base station:
Determine in described network element each bar baseline corresponding look satellite altogether;
Each bar baseline corresponding altogether depending on satellite in filter out the overall situation and look satellite altogether;
The satellite with maximum height angle is filtered out as described network element overall situation reference satellite in the described overall situation is altogether depending on satellite.
4. the method for claim 1, is characterized in that, according to the two difference atmospheric parameters on network element baseline, obtains the two difference atmospheric parameters on virtual baseline, comprising:
Calculate the two difference blur level on two baselines at main reference base station place by Kalman filtering and be fixed, then estimating baseline two difference ionosphere delay and two poor tropospheric delay;
Mathematical interpolation is carried out to the two difference ionosphere delays on described two baselines and two poor tropospheric delay, obtains the two difference ionosphere delays on described virtual baseline and two poor tropospheric delay.
5. the method for claim 1, is characterized in that, described two difference atmospheric parameter is converted to single poor atmospheric parameter, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation, is specially:
Carry out mathematical modeling by single differential mode type, the two difference ionosphere delays on the described virtual baseline obtained and two poor tropospheric delay are converted to single poor delay, and generate carrier phase corresponding to each frequency of each system and pseudorange dummy observation.
6. a generating apparatus for virtual observation data, is characterized in that, comprising:
First determination module, for determining main reference station in base station network element residing for terminal and network element;
Second determination module, for according to the base station network element determined, determines that the network element overall situation looks satellite and reference satellite altogether in the network element of described base station;
Acquisition module, for obtaining the two difference atmospheric parameters on the network element baseline of described main reference base station place;
Processing module, for obtaining the two difference atmospheric parameters on virtual baseline according to the two difference atmospheric parameters on described network element baseline;
Generation module, for the two difference atmospheric parameters on described virtual baseline are converted to single poor atmospheric parameter, and generates each system and carrier phase corresponding to each frequency and pseudorange dummy observation.
7. device as claimed in claim 6, it is characterized in that, described first determination module, specifically for obtain according to the rough coordinates of terminal terminal and CORS system CORS net in the distance of each base station, and filter out the base station nearest with described terminal; Using the network element residing for described base station as alternative network element; Obtain the distance of rough coordinates position baseline corresponding to nearest base station to alternative network element of terminal, select the baseline that distance terminal is nearest; The network element that the base station at the nearest base station of distance terminal and described nearest baseline two ends is formed as the described base station network element at terminal place, and using described nearest base station main reference station in network element.
8. device as claimed in claim 6, is characterized in that, described second determination module, specifically for determine in the network element that described main reference base station is corresponding each base station corresponding look satellite altogether; Each base station corresponding altogether depending on satellite in filter out the overall situation and look satellite altogether; Filter out in the described overall situation is altogether depending on satellite there is maximum height angle satellite as reference satellite overall in described network element.
9. device as claimed in claim 6, it is characterized in that, described acquisition module, specifically for being calculated the two difference blur leveles on two baselines at main reference base station place by Kalman filtering, calculates the two difference ionosphere delays on baseline and two poor tropospheric delay after being fixed; Mathematical interpolation is carried out to the two difference ionosphere delays on described two baselines and two poor tropospheric delay, obtains the two difference ionosphere delays on described virtual baseline and two poor tropospheric delay.
10. device as claimed in claim 6, it is characterized in that, described generation module, specifically for carrying out mathematical modeling by single differential mode type, two difference ionosphere delays on the described virtual baseline obtained and two poor tropospheric delay are converted to single poor delay, and generate carrier phase corresponding to each frequency of each system and pseudorange dummy observation.
CN201510822773.2A 2015-11-24 2015-11-24 Virtual observation data generation method and device Pending CN105549055A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510822773.2A CN105549055A (en) 2015-11-24 2015-11-24 Virtual observation data generation method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510822773.2A CN105549055A (en) 2015-11-24 2015-11-24 Virtual observation data generation method and device

Publications (1)

Publication Number Publication Date
CN105549055A true CN105549055A (en) 2016-05-04

Family

ID=55828355

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510822773.2A Pending CN105549055A (en) 2015-11-24 2015-11-24 Virtual observation data generation method and device

Country Status (1)

Country Link
CN (1) CN105549055A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108345017A (en) * 2018-01-04 2018-07-31 千寻位置网络有限公司 New network RTK air interpolating methods
WO2019062030A1 (en) * 2017-09-26 2019-04-04 东南大学 Star network-based bds/gps broadcast network rtk algorithm
CN110618438A (en) * 2019-09-09 2019-12-27 广州市中海达测绘仪器有限公司 Atmospheric error calculation method and device, computer equipment and storage medium
CN110618435A (en) * 2019-11-07 2019-12-27 广东星舆科技有限公司 Observation data generation method and text data generation method for improving VRS stability
CN110674603A (en) * 2019-09-04 2020-01-10 北斗(青岛)导航位置服务有限公司 GNSS observation data simulation method and system
CN111064595A (en) * 2018-10-17 2020-04-24 千寻位置网络有限公司 Method and device for selecting baseline network type and evaluating service quality
CN112630811A (en) * 2020-12-12 2021-04-09 湘潭大学 Real-time PPP-RTK combined positioning method
CN113031037A (en) * 2021-05-25 2021-06-25 腾讯科技(深圳)有限公司 Device positioning method and device, electronic device and computer readable medium
CN113079516A (en) * 2020-08-14 2021-07-06 中移(上海)信息通信科技有限公司 Method, device and equipment for determining base station and computer storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002061449A1 (en) * 2000-12-15 2002-08-08 Trimble Navigation Limited Method and apparatus for correcting the delay error of gps signals
JP2002318273A (en) * 2001-04-20 2002-10-31 Jekku:Kk Mobile station positioning system utilizing unidirectional communication
CN101770033A (en) * 2010-02-08 2010-07-07 东南大学 Fixing method of ambiguity network between CORS and system station
CN101943749A (en) * 2010-09-10 2011-01-12 东南大学 Method for positioning network RTK based on star-shaped virtual reference station
CN102590840A (en) * 2012-02-10 2012-07-18 中国测绘科学研究院 Satellite positioning carrier phase difference method
CN104502935A (en) * 2014-12-29 2015-04-08 东南大学 Network RTK (real-time kinematic) ambiguity resolution method based on un-differential uncombined model

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002061449A1 (en) * 2000-12-15 2002-08-08 Trimble Navigation Limited Method and apparatus for correcting the delay error of gps signals
JP2002318273A (en) * 2001-04-20 2002-10-31 Jekku:Kk Mobile station positioning system utilizing unidirectional communication
CN101770033A (en) * 2010-02-08 2010-07-07 东南大学 Fixing method of ambiguity network between CORS and system station
CN101943749A (en) * 2010-09-10 2011-01-12 东南大学 Method for positioning network RTK based on star-shaped virtual reference station
CN102590840A (en) * 2012-02-10 2012-07-18 中国测绘科学研究院 Satellite positioning carrier phase difference method
CN104502935A (en) * 2014-12-29 2015-04-08 东南大学 Network RTK (real-time kinematic) ambiguity resolution method based on un-differential uncombined model

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019062030A1 (en) * 2017-09-26 2019-04-04 东南大学 Star network-based bds/gps broadcast network rtk algorithm
CN108345017A (en) * 2018-01-04 2018-07-31 千寻位置网络有限公司 New network RTK air interpolating methods
CN111064595A (en) * 2018-10-17 2020-04-24 千寻位置网络有限公司 Method and device for selecting baseline network type and evaluating service quality
CN110674603A (en) * 2019-09-04 2020-01-10 北斗(青岛)导航位置服务有限公司 GNSS observation data simulation method and system
CN110674603B (en) * 2019-09-04 2024-03-01 北斗(青岛)导航位置服务有限公司 GNSS observation data simulation method and system
CN110618438A (en) * 2019-09-09 2019-12-27 广州市中海达测绘仪器有限公司 Atmospheric error calculation method and device, computer equipment and storage medium
CN110618438B (en) * 2019-09-09 2022-05-27 广州市中海达测绘仪器有限公司 Atmospheric error calculation method and device, computer equipment and storage medium
CN110618435B (en) * 2019-11-07 2021-11-23 广东星舆科技有限公司 Observation data generation method and text data generation method for improving VRS stability
CN110618435A (en) * 2019-11-07 2019-12-27 广东星舆科技有限公司 Observation data generation method and text data generation method for improving VRS stability
CN113079516A (en) * 2020-08-14 2021-07-06 中移(上海)信息通信科技有限公司 Method, device and equipment for determining base station and computer storage medium
CN113079516B (en) * 2020-08-14 2021-10-26 中移(上海)信息通信科技有限公司 Method, device and equipment for determining base station and computer storage medium
CN112630811A (en) * 2020-12-12 2021-04-09 湘潭大学 Real-time PPP-RTK combined positioning method
CN113031037B (en) * 2021-05-25 2021-08-06 腾讯科技(深圳)有限公司 Device positioning method and device, electronic device and computer readable medium
CN113031037A (en) * 2021-05-25 2021-06-25 腾讯科技(深圳)有限公司 Device positioning method and device, electronic device and computer readable medium

Similar Documents

Publication Publication Date Title
CN105549055A (en) Virtual observation data generation method and device
Psychas et al. Assessment of ionospheric corrections for PPP-RTK using regional ionosphere modelling
CN105607089B (en) A kind of method of the Big Dipper common-view time transmission based on pseudorange assisted carrier phase
CN110007320B (en) Network RTK resolving method
CN108363084A (en) Utilize the method and apparatus of satellite positioning, satellite navigation receiver, storage medium
CN109683186B (en) Method for eliminating carrier phase time transfer antenna jump of multi-satellite navigation system
CN106125110A (en) Satellite-based augmentation system Improvement Method of Localization Precision corrected based on subregion
CN105699999B (en) A kind of method of the fixed narrow lane ambiguity of Beidou ground strengthening system base station
CN105629279B (en) A kind of wide lane ambiguity fixing means between Fiducial station of the network
CN110208836B (en) GNSS high-adaptability cycle slip detection and restoration method based on Kalman filtering
CN109782313A (en) Ground Static Precise Point Positioning data processing method
JP2010522874A5 (en)
JP2015092185A (en) Method of processing radio-navigation signals
CN103698790B (en) Difference combined method between the Big Dipper and GPS dual system wide lane carrier phase mixing star
CN104656108A (en) Sparse reference station network zenith troposphere delay modeling method considering elevation difference
CN103529459A (en) Precise positioning method through combination of single-frequency GPS and GLONASS and system thereof
CN110058282A (en) A kind of PPP high-precision locating method based on double frequency GNSS smart phone
CN111694030A (en) BDS local difference method and system based on grid virtual observation value
CN105445755A (en) CORS region ionosphere modeling method
CN103197335A (en) Method using improved regularization method to restrain difference global positioning system (DGPS) integer ambiguity ill-condition
CN104483689A (en) Determination method for BDS reference station three-frequency carrier phase whole cycle ambiguities
CN108535749B (en) Positioning enhancement method and system based on CORS and positioning system
CN110824505B (en) Deviation estimation method and system, positioning method and terminal of GNSS satellite receiver
CN104199061A (en) Method for establishing carrier phase frequency standard of GPS (global position system) and BDS (BeiDou Navigation Satellite system)
CN103675858B (en) Dipper system B1 and gps system L1 carrier phase mixing difference method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160504