[go: up one dir, main page]

CN119689532B - Global post-processing positioning method and system based on various types of Beidou enhanced services - Google Patents

Global post-processing positioning method and system based on various types of Beidou enhanced services Download PDF

Info

Publication number
CN119689532B
CN119689532B CN202510203989.4A CN202510203989A CN119689532B CN 119689532 B CN119689532 B CN 119689532B CN 202510203989 A CN202510203989 A CN 202510203989A CN 119689532 B CN119689532 B CN 119689532B
Authority
CN
China
Prior art keywords
positioning
post
modes
mobile station
mode
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.)
Active
Application number
CN202510203989.4A
Other languages
Chinese (zh)
Other versions
CN119689532A (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.)
Leador Spatial Information Technology Co ltd
Original Assignee
Leador Spatial Information 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 Leador Spatial Information Technology Co ltd filed Critical Leador Spatial Information Technology Co ltd
Priority to CN202510203989.4A priority Critical patent/CN119689532B/en
Publication of CN119689532A publication Critical patent/CN119689532A/en
Application granted granted Critical
Publication of CN119689532B publication Critical patent/CN119689532B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention discloses a global post-processing positioning method based on Beidou various types of enhancement services, which comprises the steps of acquiring three-dimensional position information of a mobile station at different positions by using various positioning modes, calculating quality factors of the mobile station at different positions in various positioning modes, acquiring high-precision three-dimensional positions of the mobile station at different positions by using reference equipment, classifying and training a random forest model by using the three-dimensional position information and the quality factors of the mobile station at different positions in various positioning modes and the high-precision three-dimensional positions acquired by the reference equipment, inputting the quality factors of the mobile station at a certain position in the various positioning modes into the trained random forest classification model to obtain position variances of the mobile station at various positioning modes, and fusing the three-dimensional position information of the mobile station at various positioning modes by using self-adaptive weighting to obtain optimal position information. According to the invention, different types of positioning results are fused, and the availability, reliability and positioning accuracy of Beidou positioning are improved.

Description

Global post-processing positioning method and system based on Beidou various types of enhanced services
Technical Field
The invention belongs to the technical field of satellite navigation positioning, and particularly relates to a global post-processing positioning method and system based on Beidou various types of enhancement services.
Background
The Beidou satellite navigation system (Beidou system for short) is built and developed according to a three-step strategy. The Beidou I system is put into use in 2000, and adopts an active positioning system to provide pseudo-range single-point positioning, time service, wide area difference and short message communication service for Chinese users. The Beidou No. two system is put into use in 2012, and a passive positioning system is added on the basis of being compatible with the Beidou No. one system technical system, so that positioning, speed measurement, time service and short message communication service are provided for users in the asia-Tai area. The Beidou No. three system 2020 is put into use, and on the basis of the Beidou No. two system, the performance and the expansion function are further improved, and various services are provided, including global-oriented positioning navigation time service, global short message communication and international search and rescue service, satellite-based enhancement, satellite-based precise single-point positioning, foundation enhancement and regional short message communication service are provided in China and peripheral areas, and post differential positioning and post precise single-point positioning service are supported.
The pseudo-range single-point positioning (Single Point Positioning, SPP) is a mode for realizing positioning by using pseudo-range observation values and broadcast ephemeris of the Beidou system, and is a basic service provided by the Beidou system. The coverage range of the system is global, the positioning accuracy is within 5m, the system can be obtained in real time and calculated afterwards, and the positioning convergence time is 10 seconds.
The satellite-based precise single point positioning (Precise Point Positioning, PPP) is realized by using the Beidou system GEO satellite to broadcast PPP-B2B satellite-based differential correction, integrity information and other information, and dual-frequency observation value enhancement service is provided for users in China and peripheral areas. The coverage range of the system is China and surrounding areas, the positioning accuracy is within 0.4m, the system can be obtained in real time, and the positioning convergence time is 20 minutes.
The foundation enhancement positioning is based on a network RTK (Real Time Kinematic) positioning technology, and is characterized in that a data processing center processes synchronous observation data of a plurality of Beidou system reference stations covered in a certain range, differential data are generated, the differential data are broadcast through mobile communication or the Internet, and users in the area receive satellite signals and differential signals, so that high-precision real-time dynamic positioning is realized. The coverage area is the area covered by mobile communication in China and surrounding areas, the positioning accuracy is within 0.05m, the positioning accuracy can be obtained in real time, and the positioning convergence time is 30 seconds.
The post differential positioning (Post Process Kinematic, PPK) is that a user uses multi-frequency pseudo-range and carrier phase data synchronously observed by a single Beidou system reference station and a mobile station to carry out double-difference ionosphere-free combination to weaken the influence of atmospheric errors, thereby realizing high-precision post processing positioning. The coverage range of the system is 20km range with a single Beidou reference station as the center, the positioning accuracy is 0.02m+1ppm multiplied by distance kilometers, the system can be calculated afterwards, and no positioning convergence time exists.
The post-precision single point positioning (PPP, precise Point Positioning) is a high-precision post-processing positioning mode realized by correcting a series of precision error correction models by using a non-differential non-combination observation model which is carried out by a user by utilizing the multi-frequency pseudo-range and carrier phase observation value of a Beidou system mobile station and the simultaneous IGS precision ephemeris and precision satellite clock difference. The coverage range of the method is global, the positioning accuracy is within 0.2m, the method can be solved afterwards, and the positioning convergence time is not needed.
Based on pseudo-range single-point positioning, satellite-based precise single-point positioning, foundation enhanced positioning, post differential positioning and post precise single-point positioning, the Beidou system can provide enhanced services and positioning capabilities of meter-level, decimeter-level, centimeter-level and other types. However, as coverage ranges, positioning accuracy, convergence time and external dependence conditions of different types of enhanced services are different, how to make a flexible fusion strategy, the best comprehensive position information is obtained by taking the advantages and the advantages of the best comprehensive position information, and the availability, reliability and accuracy of Beidou positioning service are improved.
Disclosure of Invention
Aiming at the problems of insufficient availability, reliability and precision of Beidou positioning service caused by the difference of the existing positioning technology, the invention provides a global post-processing positioning method based on Beidou various types of enhancement services. The method has strong practicability and high flexibility, improves the availability, reliability and positioning accuracy of Beidou positioning, and can be widely applied to the field of Beidou navigation and positioning.
In order to achieve the above purpose, the technical scheme provided by the invention is a global post-processing positioning method based on Beidou various types of enhancement services, which comprises the following steps:
Step 1, acquiring three-dimensional position information of a mobile station at different positions by using multiple positioning modes based on a Beidou system, calculating and recording quality factors of the mobile station at the different positions in various positioning modes, and acquiring high-precision three-dimensional positions of the mobile station at the different positions by using reference equipment based on the Beidou system;
Step2, classifying and training the random forest model by using the three-dimensional position information and quality factors of the mobile station obtained in the step1 under a plurality of positioning modes at different positions and the high-precision three-dimensional positions of the mobile station at different positions obtained by the reference equipment to obtain a trained random forest classification model;
And 3, inputting quality factors of the mobile station in a certain position in a plurality of positioning modes into the random forest classification model trained in the step 2 to obtain position variances in the plurality of positioning modes, and fusing three-dimensional position information in the plurality of positioning modes by using an adaptive weighting method to obtain an optimal positioning result in the certain position.
Further, the positioning modes in the step 1 include two types of real-time positioning and post-processing positioning, wherein the real-time positioning includes three modes of pseudo-range single-point positioning, satellite-based precise single-point positioning and foundation enhancement positioning, and the post-processing positioning includes two modes of post-differential positioning and post-precise single-point positioning. The quality factors of the pseudo-range single-point positioning mode, the satellite-based precise single-point positioning mode and the foundation enhanced positioning mode comprise satellite quantity, position error standard deviation, a positioning resolving mode, differential time delay and precision factors. The quality factors of the post differential positioning mode include the use of satellite particles, standard deviation of position error, positioning solution mode, precision factor, signal to noise ratio, baseline length, observed data integrity, multipath error estimate, carrier phase noise and fixed check ratio of ambiguity. The quality factors of the post-precise single-point positioning mode include the use of satellite numbers, standard deviation of position errors, positioning solution mode, precision factors, signal to noise ratio, observed data integrity, multipath error estimates, carrier phase noise and fixed check ratio of ambiguity.
Further, in the step 2, the three-dimensional positions of the mobile station at different positions in the plurality of positioning modes obtained in the step 1 and the high-precision three-dimensional positions obtained by using the reference device at the positions are subjected to difference to obtain a plurality of position deviations, the position deviations are classified according to the set precision classification grades, the corresponding position variances are obtained, and the grading rule is as follows:
(7)
In the formula, Representing the variance of the position of the rover station in the ith position mode at position p,Representing the difference between the three-dimensional position of the rover station in the ith position mode at position p and the high-accuracy three-dimensional position obtained by the reference apparatus,All are set position deviation thresholds, and n is the set precision classification grade number.
Inputting the quality factors of the mobile station in the various positioning modes at different positions obtained in the step 1 and the calculated precision classification grades corresponding to the position deviations of the mobile station in the various positioning modes at different positions into a random forest model, training the random forest classifier, and evaluating the classification accuracy of the random forest classifier by using a five-fold cross validation method.
Further, in the step3, the quality factors of the mobile station in different types of positioning modes at a certain position are input into the random forest classification model trained in the step 2 to obtain the precision classification level of the mobile station in different types of positioning modes at the position, and further obtain the position variance corresponding to the three-dimensional position information in various positioning modes at the position, wherein the self-adaptive weight of the three-dimensional position information of the mobile station in the ith positioning mode at the position p is assumed to beThree-dimensional position information in various positioning modes are independent of each other and are unbiased estimates of optimal position information, then:
(8)
(9)
In the formula, Representing the mathematical expectation that the data will be,Three-dimensional position information of the mobile station in the ith and jth positioning modes at the position p is respectively represented, N represents the number of types of positioning modes,An optimal positioning result of the mobile station to be solved;
In order to fuse variances The minimum value is obtained by carrying out multi-element function extremum calculation through Lagrangian multiplier method, the minimum self-adaptive weight value of different kinds of positioning modes and the minimum value of fusion variance are obtained, the three-dimensional position information under various positioning modes and the minimum self-adaptive weight value corresponding to the three-dimensional position information are multiplied, and the optimal positioning result of the mobile station is obtained through accumulationThe specific calculation mode is as follows:
(10)
(11)
(12)
In the formula, The minimum adaptive weight representing the ith positioning pattern at position p,The variance of the optimal positioning result is represented,The position variance corresponding to the three-dimensional position information in the ith positioning mode at the position p.
The invention also provides a global post-processing positioning system based on the Beidou various types of enhancement services, which is used for realizing the global post-processing positioning method based on the Beidou various types of enhancement services.
And the system comprises a processor and a memory, wherein the memory is used for storing program instructions, and the processor is used for calling the stored instructions in the memory to execute the global post-processing positioning method based on the Beidou various types of enhancement services.
Or comprises a readable storage medium, wherein the readable storage medium is stored with a computer program, and the computer program realizes the global post-processing positioning method based on the Beidou various types of enhancement services when being executed.
Compared with the prior art, the invention has the following advantages:
1) Aiming at pseudo-range single-point positioning, star-based precise single-point positioning, foundation-enhanced three-dimensional position information and precise single-point positioning and differential positioning three-dimensional position information obtained by post-processing, the application range, positioning precision, convergence time and the like of various positioning modes are comprehensively analyzed, precision classification is carried out by machine learning according to the three-dimensional position information of various types and quality factors thereof, and optimal comprehensive three-dimensional position information is obtained according to self-adaptive weighted fusion.
2) The method can realize flexible fusion and comprehensive application of different types of positioning modes, improves the availability, reliability and positioning precision of Beidou positioning, and can be widely applied to the fields of Beidou navigation positioning such as land measurement, aviation mapping, ocean navigation surveying and the like.
Drawings
In order to more clearly illustrate the invention or the technical solutions of the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the invention, and other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a flowchart of a global post-processing positioning method based on Beidou various types of enhanced services according to an embodiment of the present invention.
FIG. 2 is a diagram of a random forest classification modeling process according to an embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings and examples of the present invention, and it is apparent that the described examples are some, but not all, examples of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1
As shown in fig. 1, the embodiment of the invention provides a global post-processing positioning method based on Beidou various types of enhancement services, which comprises the following steps:
Step 1, based on different real-time positioning modes (pseudo-range single-point positioning, satellite-based precise single-point positioning and foundation enhancement positioning) set by a Beidou system according to a mobile station receiver, three-dimensional position information and quality factors of the mobile station at different positions in different modes are obtained in real time, and meanwhile, based on the Beidou system, high-precision three-dimensional positions of the mobile station at different positions are obtained by using reference equipment.
The quality factors include the number of satellites used, standard deviation of position errors, positioning solution mode, differential time delay, precision factors, etc. The satellite number is used as the total number of satellites participating in the positioning calculation. The standard deviation of the position error is the measurement error of latitude, longitude and altitude positioning solution, if the least square positioning solution is used, the measurement error is represented by the least square residual error, and if the Kalman filtering positioning solution is used, the measurement error is represented by the position state variance. The positioning resolving mode of the pseudo-range single-point positioning mode is a pseudo-range single-point positioning mode. The positioning resolving mode of the star-based precise single-point positioning comprises a pseudo-range single-point positioning mode, a fixed solution mode of the star-based precise single-point positioning and a floating solution mode of the star-based precise single-point positioning. The positioning resolving mode of the foundation enhanced positioning comprises a pseudo-range single-point positioning mode, a foundation enhanced positioning pseudo-range differential mode, a fixed solution mode of the foundation enhanced positioning and a floating solution mode of the foundation enhanced positioning. The differential time delay is the difference between the differential correction time used and the current observation time. The precision factor is calculated by a horizontal precision factor, an elevation precision factor and a time precision factor, and the calculation formula is as follows:
(1)
In the formula, GDOP represents a precision factor, HDOP represents a horizontal precision factor, VDOP represents an elevation precision factor, and TDOP represents a time precision factor.
And 2, when the coverage range of the erected reference station comprises the position of the mobile station, performing post PPK calculation according to the acquired original observation values of the reference station and the Beidou system of the mobile station to obtain three-dimensional position information and quality factors of the mobile station at different positions, and acquiring high-precision three-dimensional positions of the mobile station at different positions by using reference equipment.
The quality factors include the number of satellites used, standard deviation of position error, positioning solution mode, accuracy factor, signal to noise ratio, baseline length, observed data integrity, multipath error estimate, carrier phase noise, fixed ambiguity rate, etc. The satellite number is used as the total number of satellites participating in the positioning calculation. The standard deviation of the position error is the measurement error of latitude, longitude and altitude positioning solution, if the least square positioning solution is used, the measurement error is represented by the least square residual error, and if the Kalman filtering positioning solution is used, the measurement error is represented by the position state variance. The positioning resolving mode of the post differential positioning mode comprises a post pseudo-range single-point positioning mode, a post differential positioning pseudo-range differential mode, a post differential positioning fixed solution mode and a post differential positioning floating solution mode. The precision factor is calculated by a horizontal precision factor, an elevation precision factor and a time precision factor, and the calculation mode is the same as the formula (1). The signal-to-noise ratio is the ratio of the carrier signal power to the noise power spectral density. The baseline length is the relative distance calculated from the known reference station coordinates and the resolved rover coordinates.
The integrity of the observed data comprises the integrity rate of the observed data of the single frequency pointAnd single system observation data integrity rateThe calculation is performed by the following formula:
(2)
(3)
In the formula, Represents the complete rate of single-frequency point observation data, n represents the number of satellites observed in an observation period,Representing the actual total number of observation epochs of the jth satellite at a certain frequency point in the observation period,Represents the theoretical total number of epochs of the jth satellite at a certain frequency point in the observation period,Representing the overall rate of observed data for a single system,The epoch number representing valid observation data for all frequency points of the jth satellite during the observation period,Representing the theoretical total number of epochs for the j-th satellite during the observation period.
Multipath error estimationThe calculation is performed by the following formula:
(4)
In the formula, Indicating that satellite is observedThe multipath error estimate over the frequency is determined,The number of epochs representing the sliding window,Expressed in calendar elementObserve satellite presenceThe frequency includes the amount of calculation of the multipath error and the integer ambiguity information.
Carrier phase noiseThe calculation is performed by the following formula:
(5)
In the formula, The number of three differences of the observed quantity of the carrier phases of the adjacent epochs of the satellite at a certain frequency point is shown,Expressed in calendar elementThe observed quantity of the phase carrier phase of the satellite at a certain frequency point is observed,And the group difference value (third difference value) representing the second difference value of the carrier phase observed quantity of the adjacent epoch of a certain frequency point.
The fixed ambiguity resolution ratio is calculated using the following formula:
(6)
Where ratio represents the fixed check ratio of ambiguity, Representing the sum of squares of the next smallest residuals in the fixed solution,Representing the sum of squares of the smallest residuals in the fixed solution.
And 3, performing post PPP calculation according to the acquired original observation value of the Beidou system of the mobile station, the precise satellite ephemeris, the clock error file and the correction to acquire three-dimensional position information and quality factors of the mobile station at different positions, and acquiring high-precision three-dimensional positions of the mobile station at different positions by using reference equipment by the Beidou system.
The quality factors include the number of satellites used, standard deviation of position error, positioning solution mode, accuracy factor, signal to noise ratio, observed data integrity, multipath error estimate, carrier phase noise, fixed check ratio of ambiguity, etc. The satellite number is used as the total number of satellites participating in the positioning calculation. The standard deviation of the position error is the measurement error of latitude, longitude and altitude positioning solution, if the least square positioning solution is used, the measurement error is represented by the least square residual error, and if the Kalman filtering positioning solution is used, the measurement error is represented by the position state variance. The positioning resolving mode of the post-precision single-point positioning mode comprises a post-pseudo-range single-point positioning mode, a post-precision single-point positioning fixed solution mode and a post-precision single-point positioning floating solution mode. The precision factor is calculated by a horizontal precision factor, an elevation precision factor and a time precision factor, and the calculation mode is the same as the formula (1). The signal-to-noise ratio is the ratio of the carrier signal power to the noise power spectral density. The observation data integrity evaluation method is the same as formulas (2) and (3). The multipath error evaluation is the same as in equation (4). The carrier phase noise calculation method is the same as the formula (5). The fuzzy degree fixed check ratio evaluation mode is the same as the formula (6).
And 4, classifying and training the random forest model by using the three-dimensional position information and quality factors of the mobile station in the different positions in a plurality of positioning modes obtained in the step 1-step 3 and the high-precision three-dimensional positions of the mobile station in the different positions obtained by the reference equipment to obtain a trained random forest classification model.
And (3) carrying out difference between three-dimensional position information of the five positioning modes (pseudo-range single-point positioning, satellite-based precise single-point positioning, foundation enhancement positioning, post differential positioning and post precise single-point positioning) obtained in the steps (1) to (3) and the high-precision three-dimensional position obtained by the reference equipment to obtain a plurality of position deviations. The position deviation is classified into 8 precision classification classes according to the precision class classification table set in table 1, wherein the precision classification class of class 1 is highest, the precision classification class of class 8 is lowest, and each precision classification class has a corresponding position variance.
TABLE 1 precision class Classification List
The random forest is a supervised machine learning integration algorithm taking the decision tree as a base learner, and attribute random selection is added in the decision tree construction process, so that the problem of over-fitting can be effectively solved, and the anti-noise performance is improved. The random forest classification method is characterized in that when a plurality of decision trees are constructed, characteristic variables and labeled samples of a training data set are randomly sampled, a decision tree is obtained by sampling results each time, and decision rules and classification results which accord with the attribute of the random forest classification algorithm can be generated by each tree, and the random forest classification algorithm is realized by integrating the decision rules and the classification results of all the trees.
And (3) inputting the quality factors of the mobile station in the various positioning modes at different positions obtained in the steps (1) to (3) and the calculated precision classification grades corresponding to the position deviations of the mobile station in the various positioning modes at different positions into a random forest model, and training the random forest classifier. Meanwhile, the accuracy of classification of the random forest classifier is evaluated by using a five-fold cross validation method, namely, five-fold repeated wheel flows are used for model training by using four equal parts of data, and the classification accuracy of the random forest classifier is tested by using one part of data, as shown in figure 2, so that the accuracy of the overall random forest classification model is improved.
And 5, inputting quality factors of the mobile station in a plurality of positioning modes at a certain position into the random forest classification model trained in the step 4 to obtain the position variance in the various positioning modes, and fusing the three-dimensional position information in the various positioning modes by using an adaptive weighting method to obtain an optimal positioning result at the certain position.
The redundant measurement information of the multi-source data can reduce the measurement error of a single data source, and the self-adaptive weighted fusion method is based on the multi-source measurement data and the corresponding characteristic standard deviation thereof, and takes the fusion variance as the principle of minimum mean square error, the weighting factors of the multi-source measurement data are automatically estimated and obtained, so that the fusion measurement data result is the optimal result.
And (3) inputting quality factors of the mobile station in different types of positioning modes at a certain position into the random forest classification model trained in the step (4) to obtain precision classification grades of the position in the different types of positioning modes, and further obtaining position variances corresponding to three-dimensional position information in the different types of positioning modes at the position. Assume that the adaptive weight of three-dimensional position information of the mobile station in the ith positioning mode at position p isThree-dimensional position information in various positioning modes are independent of each other and are unbiased estimates of optimal position information, then:
(8)
(9)
In the formula, Representing the mathematical expectation that the data will be,Three-dimensional position information of the mobile station in the ith and jth positioning modes at the position p is respectively represented, N represents the number of types of positioning modes,And (5) obtaining an optimal positioning result of the mobile station to be solved.
In order to fuse variancesThe minimum value is obtained by carrying out multi-element function extremum calculation through Lagrangian multiplier method, the minimum self-adaptive weight value of various positioning modes and the minimum value of fusion variance are obtained, the three-dimensional position information under various positioning modes and the corresponding minimum self-adaptive weight value are multiplied, and the optimal positioning result of the mobile station is obtained by accumulationThe specific calculation mode is as follows:
(9)
(10)
(11)
In the formula, The minimum adaptive weight representing the ith positioning pattern at position p,The variance of the optimal positioning result is represented,The position variance corresponding to the three-dimensional position information in the ith positioning mode at the position p.
The method can realize flexible fusion of various Beidou positioning modes by utilizing the self-adaptive weighting method, improves the availability, reliability and positioning precision of Beidou positioning, and expands the comprehensive application of Beidou.
Example 2
Based on the same inventive concept, the invention also provides a global post-processing positioning system based on the Beidou various types of enhancement services, which comprises a processor and a memory, wherein the memory is used for storing program instructions, and the processor is used for calling the program instructions in the memory to execute the global post-processing positioning method based on the Beidou various types of enhancement services.
Example 3
Based on the same inventive concept, the invention also provides a global post-processing positioning system based on the Beidou various types of enhancement services, which comprises a readable storage medium, wherein a computer program is stored on the readable storage medium, and the global post-processing positioning method based on the Beidou various types of enhancement services is realized when the computer program is executed.
In particular, the method according to the technical solution of the present invention may be implemented by those skilled in the art using computer software technology to implement an automatic operation flow, and a system apparatus for implementing the method, such as a computer readable storage medium storing a corresponding computer program according to the technical solution of the present invention, and a computer device including the operation of the corresponding computer program, should also fall within the protection scope of the present invention.
The specific embodiments described herein are offered by way of example only to illustrate the spirit of the invention. Those skilled in the art may make various modifications or additions to the described embodiments or substitutions thereof without departing from the spirit of the invention or exceeding the scope of the invention as defined in the accompanying claims.

Claims (9)

1.一种基于北斗各类型增强服务的全域后处理定位方法,其特征在于,包括以下步骤:1. A global post-processing positioning method based on various types of Beidou enhanced services, characterized in that it includes the following steps: 步骤1,基于北斗系统使用多种定位模式获取流动站在不同位置处的三维位置信息,计算并记录不同位置处各种定位模式下的质量因子,同时基于北斗系统使用基准设备获取流动站在不同位置处的高精度三维位置;Step 1, using multiple positioning modes based on the BeiDou system to obtain the three-dimensional position information of the mobile station at different positions, calculating and recording the quality factors under various positioning modes at different positions, and using reference equipment based on the BeiDou system to obtain the high-precision three-dimensional position of the mobile station at different positions; 定位模式包括实时定位和后处理定位两类,实时定位包括伪距单点定位、星基精密单点定位和地基增强定位三种模式,后处理定位包括事后差分定位和事后精密单点定位两种模式;伪距单点定位、星基精密单点定位和地基增强定位三种模式的质量因子包括使用卫星颗数、位置误差标准差、定位解算模式、差分时延和精度因子,事后差分定位模式的质量因子包括使用卫星颗数、位置误差标准差、定位解算模式、精度因子、信噪比、基线长度、观测数据完整性、多路径误差估计值、载波相位噪声和模糊度固定检核比率,事后精密单点定位模式的质量因子包括使用卫星颗数、位置误差标准差、定位解算模式、精度因子、信噪比、观测数据完整性、多路径误差估计值、载波相位噪声和模糊度固定检核比率;Positioning modes include real-time positioning and post-processing positioning. Real-time positioning includes three modes: pseudo-range single-point positioning, satellite-based precise single-point positioning and ground-based enhanced positioning. Post-processing positioning includes two modes: post-differential positioning and post-precision single-point positioning. The quality factors of the three modes, pseudo-range single-point positioning, satellite-based precise single-point positioning and ground-based enhanced positioning, include the number of satellites used, standard deviation of position error, positioning solution mode, differential delay and precision factor. The quality factor of the post-differential positioning mode includes the number of satellites used, standard deviation of position error, positioning solution mode, precision factor, signal-to-noise ratio, baseline length, observation data integrity, multipath error estimate, carrier phase noise and ambiguity fixed check ratio. The quality factor of the post-precision single-point positioning mode includes the number of satellites used, standard deviation of position error, positioning solution mode, precision factor, signal-to-noise ratio, observation data integrity, multipath error estimate, carrier phase noise and ambiguity fixed check ratio. 步骤2,使用步骤1得到的流动站在不同位置处多种定位模式下的三维位置信息、质量因子,以及基准设备获取的流动站在不同位置处的高精度三维位置对随机森林模型进行分类训练,得到训练好的随机森林分类模型;Step 2, using the three-dimensional position information and quality factor of the mobile station at different positions in multiple positioning modes obtained in step 1, and the high-precision three-dimensional position of the mobile station at different positions obtained by the reference equipment to perform classification training on the random forest model to obtain a trained random forest classification model; 步骤3,将流动站在某一位置处多种定位模式下的质量因子输入到步骤2训练好的随机森林分类模型中,得到多种定位模式下的位置方差,利用自适应加权方法将各种定位模式下的三维位置信息进行融合,得到某一位置处的最优定位结果。Step 3: Input the quality factors of the mobile station in multiple positioning modes at a certain location into the random forest classification model trained in step 2 to obtain the position variance in multiple positioning modes. Use the adaptive weighting method to fuse the three-dimensional position information in various positioning modes to obtain the optimal positioning result at a certain location. 2.如权利要求1所述的一种基于北斗各类型增强服务的全域后处理定位方法,其特征在于:步骤1中使用卫星颗数为参与定位解算的卫星总颗数;位置误差标准差为纬度、经度、高度定位解算的量测误差,若使用最小二乘定位解算,则采用最小平方残差表示量测误差,若使用卡尔曼滤波定位解算,则采用位置状态方差表示量测误差;伪距单点定位模式的定位解算模式为伪距单点定位模式;星基精密单点定位的定位解算模式包括伪距单点定位模式、星基精密单点定位的固定解模式和星基精密单点定位的浮点解模式;地基增强定位的定位解算模式包括伪距单点定位模式、地基增强定位伪距差分模式、地基增强定位的固定解模式和地基增强定位的浮点解模式;事后差分定位模式的定位解算模式包括事后伪距单点定位模式、事后差分定位伪距差分模式、事后差分定位固定解模式和事后差分定位浮点解模式;事后精密单点定位模式的定位解算模式包括事后伪距单点定位模式、事后精密单点定位固定解模式和事后精密单点定位浮点解模式;差分时延为所使用的差分改正数时刻与当前观测值时刻之差。2. A global post-processing positioning method based on various types of Beidou enhanced services as described in claim 1, characterized in that: the number of satellites used in step 1 is the total number of satellites involved in the positioning solution; the standard deviation of the position error is the measurement error of the latitude, longitude, and altitude positioning solution. If the least squares positioning solution is used, the least square residual is used to represent the measurement error. If the Kalman filter positioning solution is used, the position state variance is used to represent the measurement error; the positioning solution mode of the pseudo-range single-point positioning mode is the pseudo-range single-point positioning mode; the positioning solution mode of the satellite-based precise single-point positioning includes the pseudo-range single-point positioning mode, the fixed solution mode of the satellite-based precise single-point positioning, and the satellite-based precise single-point positioning mode. The positioning solution modes of post differential positioning mode include post pseudorange single point positioning mode, post differential positioning pseudorange difference mode, post differential positioning fixed solution mode and post differential positioning floating point solution mode; the positioning solution modes of post precise point positioning mode include post pseudorange single point positioning mode, post precise single point positioning fixed solution mode and post precise single point positioning floating point solution mode; the differential delay is the difference between the time of the differential correction used and the time of the current observation value. 3.如权利要求1所述的一种基于北斗各类型增强服务的全域后处理定位方法,其特征在于:步骤1中精度因子由水平精度因子、高程精度因子、时间精度因子计算得到,计算公式为:3. A global post-processing positioning method based on various types of Beidou enhanced services as claimed in claim 1, characterized in that: the precision factor in step 1 is calculated by the horizontal precision factor, the elevation precision factor, and the time precision factor, and the calculation formula is: (1) (1) 式中,GDOP表示精度因子,HDOP表示水平精度因子,VDOP表示高程精度因子,TDOP表示时间精度因子;Where GDOP represents the Dilution of Precision, HDOP represents the Horizontal Dilution of Precision, VDOP represents the Vertical Dilution of Precision, and TDOP represents the Temporal Dilution of Precision. 信噪比为载波信号功率与噪声功率谱密度之比;基线长度为根据已知基准站坐标与解算流动站坐标计算得到的相对距离。The signal-to-noise ratio is the ratio of the carrier signal power to the noise power spectral density; the baseline length is the relative distance calculated based on the known coordinates of the base station and the resolved mobile station coordinates. 4.如权利要求1所述的一种基于北斗各类型增强服务的全域后处理定位方法,其特征在于:步骤1中观测数据完整性包括单频点观测数据完整率和单系统观测数据完整率,通过下列公式进行计算:4. A global post-processing positioning method based on various types of Beidou enhanced services as claimed in claim 1, characterized in that: the observation data integrity in step 1 includes the single frequency observation data integrity rate and single system observation data completeness rate , calculated by the following formula: (2) (2) (3) (3) 式中,表示单频点观测数据完整率,n表示在观测时段内观测的卫星颗数,表示在观测时段内第j颗卫星在某频点的实际观测历元总数,表示在观测时段内第j颗卫星在某频点的理论历元总数,表示单系统观测数据完整率,表示在观测时段内第j颗卫星所有频点均有有效观测数据的历元数,表示在观测时段内第j颗卫星的理论历元总数;In the formula, represents the completeness rate of single-frequency observation data, n represents the number of satellites observed during the observation period, represents the total number of actual observation epochs of the jth satellite at a certain frequency point during the observation period, represents the total number of theoretical epochs of the jth satellite at a certain frequency point during the observation period, represents the completeness rate of single system observation data, It indicates the number of epochs during the observation period when all the frequencies of the jth satellite have valid observation data. represents the total number of theoretical epochs of the jth satellite during the observation period; 多路径误差估计值通过下列公式进行计算:Multipath Error Estimate Calculated using the following formula: (4) (4) 式中,表示观测到卫星在频率上多路径误差估计值,表示滑动窗口的历元个数,表示在历元观测到卫星在频率上包含多路径误差和整周模糊度信息的计算量。In the formula, Indicates that the satellite is observed Multipath error estimate in frequency, represents the number of epochs of the sliding window, Indicates the epoch Satellites were observed The amount of calculation that includes multipath error and integer ambiguity information on frequency. 5.如权利要求1所述的一种基于北斗各类型增强服务的全域后处理定位方法,其特征在于:步骤1中载波相位噪声通过下列公式进行计算:5. A global post-processing positioning method based on various types of Beidou enhanced services as claimed in claim 1, characterized in that: in step 1, the carrier phase noise Calculated using the following formula: (5) (5) 式中,表示观测到卫星在某频点相邻历元载波相位观测量的三次差值的个数,表示在历元观测到卫星在某频点的相载波相位观测量,表示某频点相邻历元载波相位观测量二次差值的组差值;In the formula, It indicates the number of cubic differences of the satellite carrier phase observations at adjacent epochs at a certain frequency point. Indicates the epoch The satellite phase observation at a certain frequency is observed. It represents the group difference of the secondary difference of the carrier phase observations of adjacent epochs at a certain frequency point; 模糊度固定检核比率采用下列公式计算:The fuzzy fixed check ratio is calculated using the following formula: (6) (6) 式中,ratio表示模糊度固定检核比率,表示固定解中次最小残差平方和,表示固定解中最小残差平方和。In the formula, ratio represents the fuzzy fixed check ratio, represents the sub-minimum residual sum of squares in the fixed solution, represents the minimum residual sum of squares among the stationary solutions. 6.如权利要求1所述的一种基于北斗各类型增强服务的全域后处理定位方法,其特征在于:步骤2中将步骤1得到的多种定位模式下流动站在不同位置处的三维位置与该位置处使用基准设备获取的高精度三维位置作差得到多个位置偏差,将位置偏差根据设定的精度分类等级进行划分,得到对应的位置方差,等级划分规则如下:6. A global post-processing positioning method based on various types of Beidou enhanced services as described in claim 1, characterized in that: in step 2, the three-dimensional position of the mobile station at different positions in the multiple positioning modes obtained in step 1 is subtracted from the high-precision three-dimensional position obtained at the position using the reference equipment to obtain multiple position deviations, and the position deviations are divided according to the set accuracy classification level to obtain the corresponding position variance, and the level division rules are as follows: (7) (7) 式中,表示位置p处第i种定位模式下流动站的位置方差,表示位置p处第i种定位模式下流动站三维位置与基准设备获取的高精度三维位置的差值,均为设定的位置偏差阈值,n为设定的精度分类等级数;In the formula, represents the position variance of the mobile station in the ith positioning mode at position p , It represents the difference between the three-dimensional position of the mobile station in the i -th positioning mode at position p and the high-precision three-dimensional position obtained by the reference equipment, , , , , are the set position deviation thresholds, and n is the set accuracy classification level number; 将步骤1得到的流动站在不同位置处各种定位模式下的质量因子和计算得到的流动站在不同位置处各种定位模式下的位置偏差对应的精度分类等级输入随机森林模型里,对随机森林分类器进行训练,同时利用五折交叉验证法评估随机森林分类器分类的准确度。The quality factors of the mobile station in various positioning modes at different locations obtained in step 1 and the accuracy classification levels corresponding to the position deviations of the mobile station in various positioning modes at different locations obtained in step 1 are input into the random forest model to train the random forest classifier. At the same time, the classification accuracy of the random forest classifier is evaluated using the five-fold cross-validation method. 7.如权利要求1所述的一种基于北斗各类型增强服务的全域后处理定位方法,其特征在于:步骤3中将流动站在某一位置处不同种类定位模式下的质量因子输入到步骤2训练好的随机森林分类模型中,得到该位置不同种类定位模式下的精度分类等级,进而得到该位置处不同种类定位模式下三维位置信息对应的位置方差;假设在位置p处第i种定位模式下流动站三维位置信息的自适应权值为,各种定位模式下的三维位置信息彼此互相独立并为最优位置信息的无偏估计,则:7. A global post-processing positioning method based on various types of Beidou enhanced services as described in claim 1, characterized in that: in step 3, the quality factor of the mobile station in different types of positioning modes at a certain position is input into the random forest classification model trained in step 2 to obtain the accuracy classification level under different types of positioning modes at the position, and then obtain the position variance corresponding to the three-dimensional position information under different types of positioning modes at the position; assuming that the adaptive weight of the three-dimensional position information of the mobile station in the i- th positioning mode at position p is , the three-dimensional position information under various positioning modes is independent of each other and is an unbiased estimate of the optimal position information, then: (8) (8) (9) (9) 式中,表示数学期望,分别表示位置p处流动站在第ij种定位模式下的三维位置信息,N表示定位模式的种类数,为待求解流动站最优定位结果;In the formula, represents the mathematical expectation, They represent the three-dimensional position information of the mobile station at position p in the i -th and j -th positioning modes, respectively. N represents the number of positioning modes. is the optimal positioning result of the mobile station to be solved; 为了使得融合方差 值最小,通过拉格朗日乘数法进行多元函数极值解算,得到各种定位模式的最小自适应权值及融合方差的最小值,将各种定位模式下的三维位置信息和其对应的最小自适应权值相乘,累加得到流动站最优定位结果,具体计算方式如下:In order to make the fusion variance The value is minimized, and the Lagrange multiplier method is used to solve the extreme value of the multivariate function to obtain the minimum adaptive weights and the minimum value of the fusion variance of various positioning modes. The three-dimensional position information under various positioning modes is multiplied by the corresponding minimum adaptive weights, and the optimal positioning result of the mobile station is obtained by accumulation. , the specific calculation method is as follows: (10) (10) (11) (11) (12) (12) 式中,表示位置p处第i种定位模式的最小自适应权值,表示最优定位结果的方差,为位置p处第i种定位模式下三维位置信息对应的位置方差。In the formula, represents the minimum adaptive weight of the i- th positioning mode at position p , represents the variance of the optimal positioning result, is the position variance corresponding to the three-dimensional position information in the i -th positioning mode at position p . 8.一种基于北斗各类型增强服务的全域后处理定位系统,其特征在于,包括处理器和存储器,存储器用于存储程序指令,处理器用于调用存储器中的程序指令执行如权利要求1-7任一项所述的一种基于北斗各类型增强服务的全域后处理定位方法。8. A global post-processing positioning system based on various types of Beidou enhanced services, characterized in that it includes a processor and a memory, the memory is used to store program instructions, and the processor is used to call the program instructions in the memory to execute a global post-processing positioning method based on various types of Beidou enhanced services as described in any one of claims 1-7. 9.一种基于北斗各类型增强服务的全域后处理定位系统,其特征在于,包括可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序执行时,实现如权利要求1-7任一项所述的一种基于北斗各类型增强服务的全域后处理定位方法。9. A global post-processing positioning system based on various types of Beidou enhanced services, characterized in that it includes a readable storage medium, on which a computer program is stored. When the computer program is executed, it implements a global post-processing positioning method based on various types of Beidou enhanced services as described in any one of claims 1-7.
CN202510203989.4A 2025-02-24 2025-02-24 Global post-processing positioning method and system based on various types of Beidou enhanced services Active CN119689532B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202510203989.4A CN119689532B (en) 2025-02-24 2025-02-24 Global post-processing positioning method and system based on various types of Beidou enhanced services

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202510203989.4A CN119689532B (en) 2025-02-24 2025-02-24 Global post-processing positioning method and system based on various types of Beidou enhanced services

Publications (2)

Publication Number Publication Date
CN119689532A CN119689532A (en) 2025-03-25
CN119689532B true CN119689532B (en) 2025-07-11

Family

ID=95033498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202510203989.4A Active CN119689532B (en) 2025-02-24 2025-02-24 Global post-processing positioning method and system based on various types of Beidou enhanced services

Country Status (1)

Country Link
CN (1) CN119689532B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731489A (en) * 2020-12-11 2021-04-30 国网辽宁省电力有限公司朝阳供电公司 High-precision positioning method based on seamless fusion of BDS (brain-based distributed system) satellite-based foundation enhancement system
CN115685287A (en) * 2022-10-27 2023-02-03 国网山西省电力公司阳泉供电公司 Method and device for substation regional elastic fusion positioning based on deep learning

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8059028B2 (en) * 2008-08-14 2011-11-15 Trueposition, Inc. Hybrid GNSS and TDOA wireless location system
CN115326066B (en) * 2022-08-16 2025-04-29 华东交通大学 A high-precision positioning method and system for high-speed trains
CN119355764A (en) * 2024-12-23 2025-01-24 山东理工大学 GNSS/5G fusion location information service performance monitoring method and system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731489A (en) * 2020-12-11 2021-04-30 国网辽宁省电力有限公司朝阳供电公司 High-precision positioning method based on seamless fusion of BDS (brain-based distributed system) satellite-based foundation enhancement system
CN115685287A (en) * 2022-10-27 2023-02-03 国网山西省电力公司阳泉供电公司 Method and device for substation regional elastic fusion positioning based on deep learning

Also Published As

Publication number Publication date
CN119689532A (en) 2025-03-25

Similar Documents

Publication Publication Date Title
Zhang et al. 3D mapping database aided GNSS based collaborative positioning using factor graph optimization
Wang et al. Smartphone shadow matching for better cross-street GNSS positioning in urban environments
CN115373005B (en) High-precision product conversion method between satellite navigation signals
US20140306842A1 (en) Accuracy and performance of the hybrid positioning system
US9291716B2 (en) Method and apparatus for determining physical characteristics of a receiving device in a navigation system
CN106054216A (en) Multi-mode GNSS satellite selection method based on GDOP and UERE
Zhang et al. A novel GNSS based V2V cooperative localization to exclude multipath effect using consistency checks
JP2024543946A (en) Positioning method, device, equipment and program
Dai et al. Innovative algorithms to improve long range RTK reliability and availability
Zhang et al. Ins-aided gnss pseudo-range error prediction using machine learning for urban vehicle navigation
Zhong et al. Outlier detection for 3D-mapping-aided GNSS positioning
Luo Precise relative positioning of multiple moving platforms using GPS carrier phase observables
CN112596088B (en) High-precision positioning method, device and storage medium applied to land measurement
CN119689532B (en) Global post-processing positioning method and system based on various types of Beidou enhanced services
CN111323748B (en) Differential positioning method and system
CN117630996A (en) Positioning method, device, equipment and storage medium
CN113917509B (en) Double-difference ambiguity fixing method, device and readable storage medium
Uaratanawong et al. Optimization technique for pseudorange multipath mitigation using different signal selection methods
Luo et al. Multi-Epoch Grid-Based 3DMA Positioning in Dense Urban Canyons by Optimizing Reflection Modelling
EP4490548A1 (en) Methods and systems for enhanced ransac selection of gnss signals
Wielgosz et al. Research on GNSS positioning and applications in Poland in 2015–2018
Alves et al. The effects of network geometry on network RTK using simulated GPS data
Dawod et al. Accuracy assessment of the PGM17 global geopotential model: a case study of Egypt and Northeast Africa
Próchniewicz et al. Quality indicator for ionospheric biases interpolation in the Network RTK
Roongpiboonsopit et al. Integrated global navigation satellite system (iGNSS) QoS prediction

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