CN112394376A - Non-differential network parallel processing method for large-scale GNSS network observation data - Google Patents

Non-differential network parallel processing method for large-scale GNSS network observation data Download PDF

Info

Publication number
CN112394376A
CN112394376A CN202011310618.XA CN202011310618A CN112394376A CN 112394376 A CN112394376 A CN 112394376A CN 202011310618 A CN202011310618 A CN 202011310618A CN 112394376 A CN112394376 A CN 112394376A
Authority
CN
China
Prior art keywords
differential
parallel
network
satellite
stations
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011310618.XA
Other languages
Chinese (zh)
Other versions
CN112394376B (en
Inventor
李林阳
吕志平
赵冬青
张勇
邝英才
林家乐
赖路广
杨凯淳
许炜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Information Engineering University of PLA Strategic Support Force
Original Assignee
Information Engineering University of PLA Strategic Support Force
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 Information Engineering University of PLA Strategic Support Force filed Critical Information Engineering University of PLA Strategic Support Force
Priority to CN202011310618.XA priority Critical patent/CN112394376B/en
Publication of CN112394376A publication Critical patent/CN112394376A/en
Application granted granted Critical
Publication of CN112394376B publication Critical patent/CN112394376B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/35Constructional details or hardware or software details of the signal processing chain
    • G01S19/37Hardware or software details of the signal processing chain
    • 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
    • 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
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention relates to a non-differential whole network parallel processing method for large-scale GNSS network observation data, and particularly relates to a non-differential whole network parallel processing method for large-scale GNSS network observation data, which is realized on the basis of a bottom multi-core parallel computing technology. Firstly, screening a certain number of survey stations from a large-scale GNSS network as core stations, estimating satellite orbits and satellite clock errors in parallel, performing parallel non-differential floating solution calculation between stations, and estimating the fractional phase deviation of a wide lane and a narrow lane at a satellite end in parallel; secondly, parallel computing of parallel non-differential floating solutions between stations is carried out on other stations in the GNSS network, and non-differential ambiguities of all stations are fixed in parallel by utilizing fractional phase deviations of wide lanes and narrow lanes estimated by a core station, so that carrier distances of all stations in the GNSS network are generated in parallel; and finally, performing unified solution on the whole network in a non-differential mode by using the pseudo-range observed value and the carrier distance, wherein the parameters to be estimated comprise a survey station coordinate, a satellite orbit, a receiver and satellite clock error, double-difference ambiguity, troposphere delay, decimal phase deviation, earth rotation parameters and the like.

Description

Non-differential network parallel processing method for large-scale GNSS network observation data
Technical Field
The invention relates to a non-differential network parallel processing method for large-scale GNSS network observation data, and belongs to the technical field of large-scale GNSS network observation data processing.
Background
As an important ground infrastructure of a country, a gnss (global Navigation Satellite system) reference station network consisting of hundreds or even thousands of continuously operating reference stations is established in a plurality of countries and regions. Large-scale global GNSS network data provides rich computing resources, but also poses significant challenges in data processing. The time consumed by data processing and resolving is increased in a geometric progression along with the enlargement of the scale of the GNSS network, the global GNSS large network data processing faces the difficult problems of low calculation efficiency and even no resolving, and a corresponding high-efficiency data processing model and method are urgently needed to be used as a support so as to give full play to the advantages of the large network and improve the precision, reliability and real-time performance of large network service. The search for efficient and rapid processing of large-scale GNSS network data becomes one of the current hot spots, and is receiving more and more attention and attention.
Research institutions and researchers at home and abroad have paid extensive attention and research on the aspects of improving models and algorithms, applying parallel computing technology and the like, and published documents mainly comprise: global GPS reference frame definitions of unknown sizes and Parallel resolution of large-scale GNSS networks, Global Navigation Satellite System (GNSS) Global data processing protocols for huge GNSS networks and enhanced Global Navigation Satellite System (GNSS) data processing strategies for GNSS networks, GNSS network Global Navigation Satellite System (GNSS) Global data processing strategies, GNSS network Global Navigation Satellite System (GNSS) Global navigation satellite system (mass navigation satellite system) Global data processing studies, GNSS data Parallel processing studies, large-network double-difference model Parallel fast solving methods and PPP network solution UPD ambiguity fixed large-scale CORS station Global positioning system (GNSS) Global navigation satellite system (GPS) Global surveying and earth dynamics Research abroad, and GNSS data fast processing strategies, GNSS data fast processing and fast processing methods under multi-core environment, of national survey and mapping science report.
In the non-differential whole network computing method of the large-scale GNSS network, the existing solution adopts a serial computing method to synchronously process data of all reference stations in the large-scale GNSS network, serially estimate satellite orbits and clock errors, serially estimate fractional phase deviation of wide lanes and narrow lanes, serially fix non-differential ambiguity, serially generate carrier distance and serially resolve the whole network in sequence. In parallel solution of observation data of a large-scale GNSS network, research has been focused on the aspects of baseline vector network adjustment, double-difference baseline solution, single-station non-difference floating solution, fixed solution calculation and the like.
Disclosure of Invention
The invention aims to provide a non-differential whole network parallel processing method for observation data of a large-scale GNSS network, which is used for solving the problem that the time consumption of data processing and resolving is exponentially increased due to the fact that the current GNSS network is continuously enlarged in size.
In order to achieve the above object, the scheme of the invention comprises:
the invention relates to a non-differential network parallel processing method for observation data of a large-scale GNSS network, which comprises the following steps:
1) selecting a set number of observation stations as core stations, and resolving satellite orbits, satellite clock errors and earth rotation parameters in parallel according to observation data of the core stations;
2) estimating wide lane decimal phase deviation and narrow lane decimal phase deviation of a satellite end by using the satellite orbit, the satellite clock error and the earth rotation parameters;
3) constructing a non-differential ambiguity for each core station, and fixing the combined ambiguity without the ionized layer by using the rounded wide lane decimal phase deviation and narrow lane decimal phase deviation;
for other stations except the core station, performing inter-station parallel non-differential fixed solution calculation by using satellite orbits, satellite clock errors, earth rotation parameters, wide lane decimal phase deviation and narrow lane decimal phase deviation of a satellite end, constructing non-differential ambiguity for each other station, and fixing the non-ionosphere combination ambiguity;
4) according to the carrier phase observed value of each observation station and the fixed ionosphere-free combined mold ambiguity, generating carrier distance observed values of the corresponding observation stations in parallel;
5) and combining the pseudo-range observed values and the carrier distance observed values of all the stations, parallelly constructing a method equation, parallelly inverting the method equation, and carrying out whole-network parallel solution based on a non-differential mode.
The invention provides a novel large-scale GNSS network observation data non-differential whole network multi-core parallel processing method, which realizes the parallel execution of a plurality of links of large-scale GNSS network observation data non-differential whole network calculation, improves the utilization rate of the calculation on a multi-core platform, and improves the whole network resolving efficiency.
Further, the core station is not less than 100 measuring stations which are uniformly distributed in the whole world.
Further, in the step 1), satellite orbit and earth rotation parameters are estimated through Kalman filtering in a non-differential mode, a single core station is used as a parallel granularity to construct a non-differential observation equation, and the non-differential observation equation is constrained through constructing double-differential ambiguity and fixing the ambiguity.
Further, in the step 2), the wide lane decimal phase deviation and the narrow lane decimal phase deviation at the satellite end are subjected to inter-station parallel non-differential floating solution calculation on the core station to realize parallel estimation.
Further, in the step 4), the carrier phase observed value is corrected by using the fixed ionosphere-free combined mold ambiguity and the antenna phase winding effect, so as to obtain a carrier distance observed value.
Further, in step 5), the pseudorange observed value is deducted by satellite-side pseudorange hardware delay.
Further, in step 5), the method for performing whole network parallel computation based on the non-differential mode is as follows: taking the pseudo-range observed value and the carrier distance observed value as observed quantities; establishing a normal equation of each epoch in parallel by taking the epoch as granularity; solving a normal equation by adopting a parameter elimination algorithm based on recursive least squares; and (3) inverting the normal equation by adopting a matrix block inversion method and a Cholesky parallel decomposition method.
At present, the requirement of non-differential network unified solution of observation data of a large-scale GNSS network on real-time performance and computing efficiency is higher and higher, in order to shorten the time of non-differential network computing of the observation data of the large-scale GNSS network and improve the timeliness of the non-differential network computing of the observation data of the large-scale GNSS network, a multi-core parallel computing technology, namely a Task Parallel Library (TPL), is adopted to design the non-differential network computing method of the observation data of the large-scale GNSS network under a multi-core platform, and the parallel solution of satellite orbit, satellite clock difference and earth rotation parameters of a core station is realized; the core station estimates the fractional phase deviation of the wide lane and the fractional phase deviation of the narrow lane in parallel; and performing multi-core parallel execution of whole network resolving in a non-differential mode.
The method is simple and easy to operate, improves the utilization rate of the multi-core platform, shortens the time of calculating the non-differential whole network of the observation data of the large-scale GNSS network, and improves the efficiency of resolving the non-differential whole network of the observation data of the large-scale GNSS network.
Further, the parameters to be estimated obtained by the parallel solution of the whole network include: station coordinates, satellite orbits, double-differenced ambiguities, tropospheric delays, receiver and satellite clock offsets, fractional phase deviations, and earth rotation parameters.
Further, for the non-ionospheric combination model ambiguity which is not fixed by each survey station in the step 3), constructing the non-ionospheric combination model ambiguity of the station-satellite double differences, and fixing the non-ionospheric combination model ambiguity.
For the ambiguity that each station cannot be fixed through a non-differential mode, a double-differential mode is further tried to be adopted for fixing, the success rate of ambiguity fixing is improved, so that more observation data can be utilized in the whole network resolving process, and the resolving precision and reliability of the parameters to be estimated are improved.
Drawings
FIG. 1 is a flow chart of non-differential network parallel computation of observation data of a large-scale GNSS network under a multi-core platform.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings.
As shown in fig. 1, for observation data of a large-scale GNSS network reference station (observation station), the method of the present invention discloses a process for non-differential network parallel processing of observation data of a large-scale GNSS network under a multi-core platform, which comprises the following steps:
step 1: screening a certain number of survey stations which are uniformly distributed globally in a large-scale GNSS network as core stations, and resolving satellite orbits, satellite clock errors and earth rotation parameters in parallel based on task parallelism and data parallelism;
step 2: performing interstation parallel non-differential floating-point solution calculation on a core station by using a satellite orbit, a satellite clock error and an earth rotation parameter, and estimating a wide lane decimal phase deviation and a narrow lane decimal phase deviation at a satellite end in parallel;
and step 3: parallel calculation of non-differential floating solutions of inter-station parallel is carried out on other survey stations in the large-scale GNSS network, and non-differential ambiguity is fixed for each survey station by utilizing fractional phase deviation of a wide lane and a narrow lane estimated by a core station to obtain the fixed non-differential non-ionosphere combined model ambiguity;
and 4, step 4: deducting corresponding fixed non-differential ionosphere-free combined model ambiguity of carrier phase observed values of all stations in a large-scale GNSS network based on a two-stage parallel mechanism between the stations and between satellites, converting the non-differential ionosphere-free combined model ambiguity into ambiguity-free carrier distance observed values, and deducting satellite-side pseudo range hardware delay from pseudo range observed values of all stations;
and 5: combining the pseudo range and the carrier distance obtained in the step 4, parallel constructing a method equation, parallel inverting the method equation, and carrying out whole network parallel resolving based on a non-differential mode.
In the step 1, the number of globally uniform core stations is not less than 100, the strategy of satellite orbit, satellite clock error and earth rotation parameter estimation is Kalman filtering in a non-differential mode, a single survey station is used as parallel granularity to construct a non-differential observation equation, and the non-differential equation is constrained by constructing double-differential ambiguity and fixing.
In step 3), the non-difference fixed solution parallel computing method comprises the following steps:
aiming at the core station, non-differential floating solution calculation is not needed, the wide lane decimal phase deviation and the narrow lane decimal phase deviation obtained in the step 2) are directly utilized, the whole is obtained to obtain the wide lane ambiguity and the narrow lane ambiguity, and the fixed non-differential ionosphere-free combined ambiguity is obtained through the fixed wide lane ambiguity and the fixed narrow lane ambiguity;
and for other stations in the GNSS network, performing parallel non-differential fixed solution calculation between the stations by using the satellite orbit, the satellite clock error, the earth rotation parameter, the wide lane decimal phase deviation and the narrow lane decimal phase deviation of the satellite end.
In the step 4), the carrier distance is generated by utilizing the non-differential ionosphere-free combined model ambiguity fixed in the step 3) and the antenna phase winding effect obtained by model calculation, and the ionosphere-free combined carrier phase observed quantity is corrected to obtain the carrier distance
Figure BDA0002789701390000051
Figure BDA0002789701390000052
In the formula (I), the compound is shown in the specification,
Figure BDA0002789701390000053
the distance is the carrier distance, L is the observed value of the original ionosphere-free combined carrier phase, lambda is the ionosphere-free combined wavelength, N is the fixed non-differential ionosphere-free combined pattern ambiguity, and xi is the antenna phase winding effect.
And (3) further constructing the ionosphere-free combined mold ambiguity of the satellite-satellite double-difference of the stations without realizing the non-difference fixed ionosphere-free combined mold ambiguity for each station in the step 3), trying to fix the double-difference ambiguity, effectively improving the success rate of ambiguity fixation, increasing the effective carrier distance observation value which can be used for the whole network parallel solution in the step 5), and further improving the accuracy and reliability of parameter solution to be estimated.
In step 5), the whole network parallel resolving method in the non-differential mode comprises the following steps:
firstly, the obtained pseudo range and carrier distance are used as observed quantities, when a method equation is constructed, epochs are used as granularity, the method equation of each epoch is constructed in parallel based on data parallelism, and the method equation is solved by adopting a parameter elimination algorithm based on recursive least squares; when the normal equation is used for solving the inverse, matrix block inversion, a Cholesky parallel decomposition method and the like can be adopted.
Parameters to be estimated, which are solved in the whole network in a non-differential mode, comprise station coordinates, satellite orbits, double-difference ambiguity, troposphere delay, receiver and satellite clock errors, decimal phase deviation, earth rotation parameters and the like.
The observation data of the present invention includes pseudorange and carrier phase observations for the GPS, Galileo, BDS and GLONASS systems at multiple frequencies.
The above description is only a preferred embodiment of the present invention, and the scope of the present invention is not limited thereto, and any simple modifications or equivalent substitutions of the technical solutions that can be obviously obtained by those skilled in the art within the technical scope of the present invention are within the scope of the present invention.
The invention discloses a non-differential whole network parallel processing method for large-scale GNSS network observation data, which is characterized by realizing multi-core parallel execution of a plurality of links of non-differential whole network calculation of the large-scale GNSS network observation data. The experiment utilizes the observation data of 500 GNSS reference stations which are distributed globally, and adopts three parallel schemes of dual-core, four-core and six-core respectively. Through tests, compared with the traditional single-core serial method, the large-scale GNSS network observation data non-differential whole network multi-core parallel processing method provided by the invention greatly shortens the calculation time, improves the calculation efficiency, and has the acceleration ratios of the double-core, four-core and six-core parallel reaching 1.60, 3.05 and 4.50 times respectively. The effect of practical application is closely related to the performance of a hardware system, the size of the scale of the GNSS network, the quality of observed data and the like.
Therefore, compared with the prior art, the invention has the following outstanding beneficial technical effects:
(1) the timeliness of non-differential network calculation of observation data of the large-scale GNSS network is improved.
The invention provides a non-differential whole network parallel computing method for observation data of a large-scale GNSS network, which realizes parallel resolving of satellite orbits, satellite clock errors and earth rotation parameters by a certain number of uniformly distributed core stations in the GNSS network under a multi-core platform; the core station estimates decimal phase deviation of a wide lane and a narrow lane in parallel; and the parallel computation of the whole network non-differential ambiguity parallel fixation, the whole network carrier distance parallel generation and the whole network resolving in the non-differential mode shortens the time of the large-scale GNSS network observation data non-differential whole network computation and improves the computation efficiency.
(2) Easy to expand.
The method provided by the invention has wider applicability and stronger expansibility, is suitable for processing the multi-frequency observation values of various satellite navigation systems, and is used for carrying out whole network calculation based on a non-differential mode. Observation data of systems such as GPS, Galileo, BDS, GLONASS and the like can be included in the non-differential whole network parallel computing method, the method provided by the invention can be adopted to carry out whole network parallel computing on multi-frequency observation data of a selected satellite navigation system in a non-differential mode, and the method is not limited to the 2 frequencies, but is still suitable for 3, 4 and 5 frequencies.

Claims (9)

1. A non-differential network parallel processing method for large-scale GNSS network observation data is characterized by comprising the following steps:
1) selecting a set number of observation stations as core stations, and resolving satellite orbits, satellite clock errors and earth rotation parameters in parallel according to observation data of the core stations;
2) estimating wide lane decimal phase deviation and narrow lane decimal phase deviation of a satellite end by using the satellite orbit, the satellite clock error and the earth rotation parameters;
3) constructing a non-differential ambiguity for each core station, and fixing the combined ambiguity without the ionized layer by using the rounded wide lane decimal phase deviation and narrow lane decimal phase deviation;
for other stations except the core station, performing inter-station parallel non-differential fixed solution calculation by using satellite orbits, satellite clock errors, earth rotation parameters, wide lane decimal phase deviation and narrow lane decimal phase deviation of a satellite end, constructing non-differential ambiguity for each other station, and fixing the non-ionosphere combination ambiguity;
4) according to the carrier phase observed value of each observation station and the fixed ionosphere-free combined mold ambiguity, generating carrier distance observed values of the corresponding observation stations in parallel;
5) and combining the pseudo-range observed values and the carrier distance observed values of all the stations, parallelly constructing a method equation, parallelly inverting the method equation, and carrying out whole-network parallel solution based on a non-differential mode.
2. The non-differential network parallel processing method for large-scale GNSS network observation data according to claim 1, wherein the core stations are not less than 100 stations uniformly distributed around the world.
3. The non-differential whole network parallel processing method for large-scale GNSS network observation data according to claim 1, characterized in that in step 1), satellite orbit and earth rotation parameters are estimated through Kalman filtering in a non-differential mode, a single core station is used as parallel granularity to construct a non-differential observation equation, and the non-differential observation equation is constrained through constructing double-differential ambiguity and fixing the ambiguity.
4. The large-scale GNSS network observation data non-differential whole network parallel processing method according to claim 1, wherein in step 2), the wide lane fractional phase bias and the narrow lane fractional phase bias of the satellite end are subjected to inter-station parallel non-differential floating solution calculation to realize parallel estimation.
5. The non-differential global processing method for observation data of the large-scale GNSS network according to claim 1, wherein in the step 4), the carrier phase observation value is corrected to obtain the carrier distance observation value by using the fixed ionosphere-free combined model ambiguity and the antenna phase winding effect.
6. The non-differencing network parallel processing method for large-scale GNSS network observation data according to claim 1, wherein in step 5), satellite-side pseudorange hardware delay is deducted from the pseudorange observation.
7. The non-differential whole network parallel processing method for observation data of the large-scale GNSS network according to claim 1, wherein in the step 5), the method for performing whole network parallel solution based on the non-differential mode comprises: taking the pseudo-range observed value and the carrier distance observed value as observed quantities; establishing a normal equation of each epoch in parallel by taking the epoch as granularity; solving a normal equation by adopting a parameter elimination algorithm based on recursive least squares; and (3) inverting the normal equation by adopting a matrix block inversion method and a Cholesky parallel decomposition method.
8. The non-differential global network parallel processing method for observation data of the large-scale GNSS network according to claim 1, wherein the parameters to be estimated obtained by the parallel computation of the global network comprise: station coordinates, satellite orbits, double-differenced ambiguities, tropospheric delays, receiver and satellite clock offsets, fractional phase deviations, and earth rotation parameters.
9. The non-differential whole-network parallel processing method for observation data of the large-scale GNSS network according to claim 1, wherein the ionospheric-free combined ambiguity free of the satellite double-differential is constructed for the ionospheric-free combined ambiguity free of the unfixed observation station in step 3), and the ionospheric-free combined ambiguity is fixed.
CN202011310618.XA 2020-11-20 2020-11-20 Large-scale GNSS network observation data non-differential whole network parallel processing method Active CN112394376B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011310618.XA CN112394376B (en) 2020-11-20 2020-11-20 Large-scale GNSS network observation data non-differential whole network parallel processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011310618.XA CN112394376B (en) 2020-11-20 2020-11-20 Large-scale GNSS network observation data non-differential whole network parallel processing method

Publications (2)

Publication Number Publication Date
CN112394376A true CN112394376A (en) 2021-02-23
CN112394376B CN112394376B (en) 2023-07-04

Family

ID=74605958

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011310618.XA Active CN112394376B (en) 2020-11-20 2020-11-20 Large-scale GNSS network observation data non-differential whole network parallel processing method

Country Status (1)

Country Link
CN (1) CN112394376B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113282406A (en) * 2021-04-27 2021-08-20 中国人民解放军陆军勤务学院 Collaborative parallel computing method for GNSS large-network observation data
CN113281794A (en) * 2021-04-27 2021-08-20 中国人民解放军陆军勤务学院 Method for parallel estimation of single-difference wide-lane FCB between satellites
CN114355420A (en) * 2021-12-15 2022-04-15 中国科学院国家授时中心 Distributed Beidou position service center PPP product positioning method and device
WO2023019613A1 (en) * 2021-08-19 2023-02-23 中国能源建设集团江苏省电力设计院有限公司 Large-scale gnss network parallel resolution method and system based on dynamic partioning

Citations (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090182502A1 (en) * 2008-01-10 2009-07-16 Riter Bruce D Refining a position estimate of a low earth orbiting satellite
US20090184869A1 (en) * 2008-01-09 2009-07-23 Trimble Navigation Limited, A Corporation Of California Processing Multi-GNSS data from mixed-type receivers
US20100188285A1 (en) * 2009-01-23 2010-07-29 Her Majesty The Queen In Right Of Canada As Represented By The Minister Of Natural Resources Decoupled clock model with ambiguity datum fixing
US20100202300A1 (en) * 2009-02-12 2010-08-12 Zulutime, Llc Systems and methods for space-time determinations with reduced network traffic
US20110090114A1 (en) * 2009-10-15 2011-04-21 Fenton Patrick C Ultra short baseline gnss receiver
US20110267228A1 (en) * 2009-02-22 2011-11-03 Trimble Navigation Limited GNSS Signal Processing Methods and Apparatus with Geometric Filter
JP2012194099A (en) * 2011-03-17 2012-10-11 Seiko Epson Corp Pseudo-distance error estimation method, position calculation method and pseudo-distance error estimation apparatus
US20130141276A1 (en) * 2011-09-16 2013-06-06 Trimble Navigation Limited GNSS Signal Processing Methods and Apparatus
CN103175516A (en) * 2013-02-26 2013-06-26 中国人民解放军信息工程大学 Distributed computing method for adjustment of large-scale geodesic control net
US20130332072A1 (en) * 2008-01-10 2013-12-12 Trimble Navigation Limited Refining a position estimate of a low earth orbiting satellite
WO2014065664A1 (en) * 2012-10-25 2014-05-01 Fugro N.V. Ppp-rtk method and system for gnss signal based position determination
CN104459745A (en) * 2014-12-25 2015-03-25 东南大学 Quick resolving method for multi-constellation long-base-line network RTK partial ambiguity
CN105301619A (en) * 2015-12-02 2016-02-03 武汉大学 Rapid processing method and system for whole large scale GNSS network data
CN105929430A (en) * 2016-07-14 2016-09-07 天津市勘察院 GNSS (global navigation satellite system) zero-baseline inter-reference station ambiguity quick fixation method
CN107656295A (en) * 2017-07-31 2018-02-02 武汉大学 A kind of GNSS high accuracy Baseline Survey methods based on original observed data
CN107942346A (en) * 2017-11-21 2018-04-20 武汉大学 A kind of high-precision GNSS ionized layer TEC observation extracting method
WO2018138904A1 (en) * 2017-01-30 2018-08-02 三菱電機株式会社 Position measurement device and position measurement method
CN108549095A (en) * 2018-04-12 2018-09-18 中国人民解放军战略支援部队信息工程大学 A kind of region CORS nets non-poor Enhancement Method and system parallel
CN109212562A (en) * 2018-08-29 2019-01-15 中国人民解放军61540部队 A method of generating carrier wave pseudo range observed quantity
CN109799521A (en) * 2019-03-14 2019-05-24 苏州工业园区测绘地理信息有限公司 A kind of tri- subtractive combination cycle-slip detection and repair method of BDS/GPS
CN110749911A (en) * 2019-10-28 2020-02-04 中国人民解放军战略支援部队信息工程大学 Method and device for parallel preprocessing of satellite slope path distance of clean station of large-scale GNSS network
CN111025346A (en) * 2019-11-18 2020-04-17 广州南方卫星导航仪器有限公司 Method for rapidly estimating clock error of GNSS precision satellite and storage medium
CN111045042A (en) * 2019-12-20 2020-04-21 西安空间无线电技术研究所 PPP-RTK enhancement method and system based on 'cloud-end' framework
CN111208520A (en) * 2020-01-17 2020-05-29 中国人民解放军战略支援部队信息工程大学 Positioning method and device of submarine acoustic transponder
CN111273327A (en) * 2020-03-20 2020-06-12 中国人民解放军61081部队 Precise single-point positioning method based on combined and non-combined hybrid observation model
CN111290004A (en) * 2020-03-04 2020-06-16 高维时空(北京)网络有限公司 Pseudo-range differential positioning method, pseudo-range differential positioning device, electronic equipment and storage medium

Patent Citations (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090184869A1 (en) * 2008-01-09 2009-07-23 Trimble Navigation Limited, A Corporation Of California Processing Multi-GNSS data from mixed-type receivers
US20090182502A1 (en) * 2008-01-10 2009-07-16 Riter Bruce D Refining a position estimate of a low earth orbiting satellite
US20130332072A1 (en) * 2008-01-10 2013-12-12 Trimble Navigation Limited Refining a position estimate of a low earth orbiting satellite
US20100188285A1 (en) * 2009-01-23 2010-07-29 Her Majesty The Queen In Right Of Canada As Represented By The Minister Of Natural Resources Decoupled clock model with ambiguity datum fixing
US20100202300A1 (en) * 2009-02-12 2010-08-12 Zulutime, Llc Systems and methods for space-time determinations with reduced network traffic
US20110267228A1 (en) * 2009-02-22 2011-11-03 Trimble Navigation Limited GNSS Signal Processing Methods and Apparatus with Geometric Filter
US20110090114A1 (en) * 2009-10-15 2011-04-21 Fenton Patrick C Ultra short baseline gnss receiver
JP2012194099A (en) * 2011-03-17 2012-10-11 Seiko Epson Corp Pseudo-distance error estimation method, position calculation method and pseudo-distance error estimation apparatus
US20130141276A1 (en) * 2011-09-16 2013-06-06 Trimble Navigation Limited GNSS Signal Processing Methods and Apparatus
WO2014065664A1 (en) * 2012-10-25 2014-05-01 Fugro N.V. Ppp-rtk method and system for gnss signal based position determination
CN103175516A (en) * 2013-02-26 2013-06-26 中国人民解放军信息工程大学 Distributed computing method for adjustment of large-scale geodesic control net
CN104459745A (en) * 2014-12-25 2015-03-25 东南大学 Quick resolving method for multi-constellation long-base-line network RTK partial ambiguity
CN105301619A (en) * 2015-12-02 2016-02-03 武汉大学 Rapid processing method and system for whole large scale GNSS network data
CN105929430A (en) * 2016-07-14 2016-09-07 天津市勘察院 GNSS (global navigation satellite system) zero-baseline inter-reference station ambiguity quick fixation method
WO2018138904A1 (en) * 2017-01-30 2018-08-02 三菱電機株式会社 Position measurement device and position measurement method
CN107656295A (en) * 2017-07-31 2018-02-02 武汉大学 A kind of GNSS high accuracy Baseline Survey methods based on original observed data
CN107942346A (en) * 2017-11-21 2018-04-20 武汉大学 A kind of high-precision GNSS ionized layer TEC observation extracting method
CN108549095A (en) * 2018-04-12 2018-09-18 中国人民解放军战略支援部队信息工程大学 A kind of region CORS nets non-poor Enhancement Method and system parallel
CN109212562A (en) * 2018-08-29 2019-01-15 中国人民解放军61540部队 A method of generating carrier wave pseudo range observed quantity
CN109799521A (en) * 2019-03-14 2019-05-24 苏州工业园区测绘地理信息有限公司 A kind of tri- subtractive combination cycle-slip detection and repair method of BDS/GPS
CN110749911A (en) * 2019-10-28 2020-02-04 中国人民解放军战略支援部队信息工程大学 Method and device for parallel preprocessing of satellite slope path distance of clean station of large-scale GNSS network
CN111025346A (en) * 2019-11-18 2020-04-17 广州南方卫星导航仪器有限公司 Method for rapidly estimating clock error of GNSS precision satellite and storage medium
CN111045042A (en) * 2019-12-20 2020-04-21 西安空间无线电技术研究所 PPP-RTK enhancement method and system based on 'cloud-end' framework
CN111208520A (en) * 2020-01-17 2020-05-29 中国人民解放军战略支援部队信息工程大学 Positioning method and device of submarine acoustic transponder
CN111290004A (en) * 2020-03-04 2020-06-16 高维时空(北京)网络有限公司 Pseudo-range differential positioning method, pseudo-range differential positioning device, electronic equipment and storage medium
CN111273327A (en) * 2020-03-20 2020-06-12 中国人民解放军61081部队 Precise single-point positioning method based on combined and non-combined hybrid observation model

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
DE JONG, K: "New Developments in Precise Offshore GNSS Positioning", QUANTITATIVE MONITORING OF THE UNDERWATER ENVIRONMENT, pages 3 - 12 *
HUA CHEN: "An enhanced strategy for GNSS data processing of massive networks", JOURNAL OF GEODESY, pages 857 - 867 *
YAO, YIBI: "An improved approach to model regional ionosphere and accelerate convergence for precise point positioning", ADVANCES IN SPACE RESEARCH, pages 1406 - 1415 *
吕志平: "非差非组合与组合模型定位性能对比及并行设计", 海洋测绘, pages 27 - 30 *
吕浩: "基于分布式计算技术的GNSS数据处理方法研究", 中国优秀硕士学位论文全文数据库 *
崔阳: "GNSS大网双差模型并行快速解算方法", 测绘学报, pages 848 - 856 *
张小红: "GNSS多频多***集成精密单点定位关键技术研究", 科技成果 *
苗亚哲: "非差模型与双差模型的定位精度比较", 测绘工程, pages 6 - 12 *
邹璇: "基于非差观测的网络实时动态定位方法及其在连续运行基准站跨网服务中的应用", 测绘学报, pages 1 - 5 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113282406A (en) * 2021-04-27 2021-08-20 中国人民解放军陆军勤务学院 Collaborative parallel computing method for GNSS large-network observation data
CN113281794A (en) * 2021-04-27 2021-08-20 中国人民解放军陆军勤务学院 Method for parallel estimation of single-difference wide-lane FCB between satellites
WO2023019613A1 (en) * 2021-08-19 2023-02-23 中国能源建设集团江苏省电力设计院有限公司 Large-scale gnss network parallel resolution method and system based on dynamic partioning
CN114355420A (en) * 2021-12-15 2022-04-15 中国科学院国家授时中心 Distributed Beidou position service center PPP product positioning method and device

Also Published As

Publication number Publication date
CN112394376B (en) 2023-07-04

Similar Documents

Publication Publication Date Title
CN112394376B (en) Large-scale GNSS network observation data non-differential whole network parallel processing method
CN110007320B (en) Network RTK resolving method
Hu et al. Multi-GNSS fractional cycle bias products generation for GNSS ambiguity-fixed PPP at Wuhan University
Li et al. Improving the estimation of uncalibrated fractional phase offsets for PPP ambiguity resolution
CN105699999B (en) A kind of method of the fixed narrow lane ambiguity of Beidou ground strengthening system base station
CN108549095B (en) Non-differential parallel enhancement method and system for regional CORS network
CN112462396B (en) Real-time parallel determination method for clock error of navigation satellite with high sampling rate
CN115933356B (en) High-precision time synchronization system and method for virtual atomic clock
CN109212562A (en) A method of generating carrier wave pseudo range observed quantity
Zhao et al. Investigation of real-time kinematic multi-GNSS precise point positioning with the CNES products
CN113253314A (en) Time synchronization method and system between low-earth-orbit satellites
CN110749911B (en) Method and device for parallel preprocessing of satellite slope path distance of clean station of large-scale GNSS network
Zhang et al. Modeling, refinement and evaluation of multipath mitigation based on the hemispherical map in BDS2/BDS3 relative precise positioning
Wang et al. Investigating the contribution of BDS-3 observations to multi-GNSS single-frequency precise point positioning with different ionospheric models
CN117492035A (en) Method for resolving PPP-RTK service integrity information
Li et al. Parallel computation of regional CORS network corrections based on ionospheric-free PPP
CN115308781B (en) BDGIM-assisted phase smoothing pseudo-range high-precision time transfer method
CN110907974B (en) Method and device for quickly fixing PPP ambiguity based on VRS difference
CN115561793A (en) Real-time Beidou phase decimal deviation rapid estimation method based on parallel computation
Zeng et al. Computationally efficient dual-frequency uncombined precise orbit determination based on IGS clock datum
CN114355420B (en) PPP product positioning method and device for distributed Beidou position service center
Li et al. Parallel resolution of large-scale GNSS network un-difference ambiguity
CN114460615A (en) Beidou three-new frequency point positioning method and system with additional virtual observation value
CN114563806A (en) PPP-RTK real-time positioning method and system for Android mobile equipment
Xu et al. BDS precise point positioning ambiguity resolution with high rate data and its application to seismic displacement and marine surveying

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant