CN111694021B - GNSS environment model-based single-station landslide deformation monitoring and early warning method - Google Patents

GNSS environment model-based single-station landslide deformation monitoring and early warning method Download PDF

Info

Publication number
CN111694021B
CN111694021B CN202010460631.7A CN202010460631A CN111694021B CN 111694021 B CN111694021 B CN 111694021B CN 202010460631 A CN202010460631 A CN 202010460631A CN 111694021 B CN111694021 B CN 111694021B
Authority
CN
China
Prior art keywords
gnss
early warning
monitoring
environment
observation
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
CN202010460631.7A
Other languages
Chinese (zh)
Other versions
CN111694021A (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.)
National Time Service Center of CAS
Original Assignee
National Time Service Center of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National Time Service Center of CAS filed Critical National Time Service Center of CAS
Priority to CN202010460631.7A priority Critical patent/CN111694021B/en
Publication of CN111694021A publication Critical patent/CN111694021A/en
Application granted granted Critical
Publication of CN111694021B publication Critical patent/CN111694021B/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/14Receivers specially adapted for specific applications
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/33Multimode operation in different systems which transmit time stamped messages, e.g. GPS/GLONASS
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02Alarms for ensuring the safety of persons
    • G08B21/10Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L67/00Network arrangements or protocols for supporting network services or applications
    • H04L67/01Protocols
    • H04L67/02Protocols based on web technology, e.g. hypertext transfer protocol [HTTP]
    • H04L67/025Protocols based on web technology, e.g. hypertext transfer protocol [HTTP] for remote control or remote monitoring of applications
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L67/00Network arrangements or protocols for supporting network services or applications
    • H04L67/01Protocols
    • H04L67/12Protocols specially adapted for proprietary or special-purpose networking environments, e.g. medical networks, sensor networks, networks in vehicles or remote metering networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Environmental & Geological Engineering (AREA)
  • Emergency Management (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Business, Economics & Management (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

The invention discloses a GNSS environment model-based single-station landslide deformation monitoring and early warning method, which comprises the following steps: deploying terminal equipment, acquiring data, determining the initial position of a monitoring point, calculating environment modeling source data, modeling a sunday environment, solving a correlation coefficient and applying early warning. The monitoring and early warning method mainly utilizes the low-cost consumer GNSS positioning terminal to acquire data and perform environment modeling, early identifies the landslide acceleration state through the change trend of the correlation coefficient of the environment model and effectively early warns, so that the monitoring cost is greatly reduced, and a new thought and a new method are provided for the GNSS landslide monitoring and early warning.

Description

GNSS environment model-based single-station landslide deformation monitoring and early warning method
Technical Field
The invention belongs to the technical field of monitoring and early warning, and particularly relates to a method for monitoring and early warning of landslide deformation of a single station based on a GNSS environment model.
Background
China is one of the most serious countries in the world with geological landslide disasters. According to the incomplete statistics of the national geological disaster network, the direct economic loss is 27.7 billion yuan only from the occurrence of geological disasters 6181 in 2019 nationwide; the geological disaster 948 can be successfully forecasted, 2 million persons can be possibly injured, the direct economic loss is avoided by 8.3 hundred million yuan, and the disaster prevention and reduction effects are remarkable. And the landslide disaster loss accounts for 61.2 percent of the total disaster, and the implementation of efficient landslide monitoring and early warning has important practical significance for the disaster prevention and reduction work in China.
The GNSS technology has high-precision, all-time and all-weather monitoring capability in landslide monitoring and early warning, and is an important means for directly acquiring surface three-dimensional vector deformation at present. With the modernization and the gradual improvement of GNSS, currently, there are up to more than 110 global in-orbit GNSS satellites, which become an important means for global high-precision positioning and navigation services.
The GNSS positioning technology has been widely applied to various deformation monitoring fields, in particular to landslide disaster deformation monitoring, by virtue of its unique technical advantages, and mainly uses the RTK positioning technology. The cost of a single GNSS navigation positioning terminal which can be used for high-precision landslide monitoring is generally higher (more than 3 thousands), the cost of a common measurement type double-frequency receiving terminal is at least more than ten thousand yuan, the cost of a GNSS application terminal used for consumer-level positioning is less than one hundred yuan at present, and the positioning precision is low (more than 10m), so that the landslide monitoring and early warning requirements cannot be met.
Therefore, how to utilize a single GNSS positioning terminal with low cost to realize effective early warning of landslide deformation monitoring has important practical significance for geological disaster prevention and reduction work in China.
Disclosure of Invention
Aiming at the problems of high early warning cost, low success rate and the like of the existing GNSS landslide deformation monitoring technology, the invention provides a single-station landslide deformation monitoring and early warning method based on a GNSS environment model.
The technical scheme of the invention is as follows: a GNSS environment model-based single-station landslide deformation monitoring and early warning method comprises the following steps:
s1: deploying terminal equipment and performing data acquisition
The method comprises the steps that a GNSS application terminal is deployed on a monitoring point (landslide deformation characteristic point), and after the GNSS application terminal is started, original observation data of the monitoring point are collected through a 4G or 5G wireless communication network and the like;
s2: determining the initial position of a monitoring point
After the original observation data of several hours (>1 hour) are collected, downloading a real-time broadcast ephemeris file through an Internet IGS service center, and roughly estimating the approximate position of a monitoring point by using a pseudo-range single-point positioning technology;
s3: computing environment modeling source data
Based on the broadcast ephemeris, performing integrity check on the original observation data acquired by S1, calculating satellite azimuth angle and altitude angle information of each epoch through a receiver approximate coordinate, summarizing satellite sight line information, and outputting continuous observation arc-segment observation values, wherein the initial and final observation values of the arc-segment observation values are marked as head and tail points;
s4: sunday environment modeling
Performing gross error elimination on the head and tail point observed values output by the S3, constructing an observation equation by using a plurality of head and tail point observed values, selecting a least square criterion to solve a fitting coefficient in the equation by a fitting system, and establishing an environment model according to the fitting coefficient;
s5: correlation coefficient solving
Setting intervals by taking days as units, and jointly solving correlation coefficients in the environment models according to the corresponding environment models every day;
s6: early warning application
And setting a proper threshold value in practical application, comparing and judging the threshold value and the correlation coefficient calculated in the step S5, if the correlation coefficient is lower than the threshold value, carrying out early warning, and if not, continuing monitoring.
Further, the raw observation data in S1 includes raw pseudorange and ephemeris data.
Further, the GNSS application terminal in S1 is a consumer-grade single frequency GNSS.
Further, the integrity detection method of the original observation data in S3 includes: and comparing the raw observation data collected in the step S1 with the real-time broadcast ephemeris downloaded in the step S2 every 24 hours to determine whether the raw observation data correspond to the real-time broadcast ephemeris or not, so as to detect the integrity of the raw observation data.
Further, the method for constructing the environment model in S4 includes the following steps:
s41: and (2) constructing an observation equation for a plurality of head and tail point observation values by using a formula (1):
Figure GDA0002614093530000031
in the formula (1), n represents the fitting order, a i 、b i 、c i Denotes the fitting coefficient, e x Representing the measurement height angle, x representing e x A corresponding azimuth angle; n is 360/0.5, i is {0,0.5, …,360 };
s42: selecting a least square criterion to carry out fitting system solution on the observation equation to obtain a fitting coefficient, and establishing an environment model according to the fitting coefficient by taking 0.5 degrees as a step length from 0 to 360 degrees, as shown in a formula (2):
Figure GDA0002614093530000032
in equation (2), a 'model coefficient vector, e' i Representing the virtual satellite altitude using the model coefficients back-deducted, i ═ 0,0.5, …, 360.
Further, the specific method for solving the correlation coefficient in S5 is as follows:
to set an environmental model a 'of day 1' 0 Reference, environmental model a 'for day k every 24 hours' k And a' 0 Is solved as shown in equation (3):
Figure GDA0002614093530000033
wherein, Cov (a' k ,a′ 0 ) Is a' k And a' 0 Covariance, Var [ a' k ]Is a' k Variance of (1), Var [ a' 0 ]Is a' 0 The variance of (c).
The invention has the beneficial effects that: the monitoring and early warning method mainly utilizes the low-cost consumer GNSS positioning terminal to acquire data and perform environment modeling, early identifies the landslide acceleration state through the change trend of the correlation coefficient of the environment model and effectively early warns, so that the monitoring cost is greatly reduced, and a new thought and a new method are provided for the GNSS landslide monitoring and early warning.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a view showing modeling results of GNSS environment at monitoring points for 19 consecutive weeks according to embodiment 2 of the present invention;
fig. 3 is a timing chart of three-dimensional position change and correlation coefficient of the monitoring points in embodiment 2 of the present invention.
Detailed Description
Example 1: in this embodiment, a single-station landslide deformation monitoring and early warning method based on a GNSS environment model is designed, and is mainly used for providing a new monitoring and early warning method based on a consumer-grade single-frequency GNSS, which is specifically described based on the method flowchart in fig. 1.
S1: deploying terminal equipment and performing data acquisition
The method comprises the steps that a GNSS application terminal is deployed on a monitoring point (landslide deformation characteristic point), and after the GNSS application terminal is started, original observation data of the monitoring point are collected through a 4G or 5G wireless communication network and the like, wherein the original observation data comprise original pseudo-range and ephemeris data; the GNSS application terminal can be a consumer-grade single-frequency GNSS, or can be a professional type with centimeter-level precision or a common measurement type with sub-meter-level precision, and the GNSS comprises one or more of BDS/GPS/GLONASS/GALILEO.
S2: determining the initial position of a monitoring point
After the original observation data of several hours (>1 hour) are collected, downloading a real-time broadcast ephemeris file through an Internet IGS service center, and roughly estimating the approximate position of a monitoring point by using a pseudo-range single-point positioning technology;
s3: computing environment modeling source data
On the basis of a broadcast ephemeris, comparing original observation data acquired in S1 with the real-time broadcast ephemeris downloaded in S2 every 24 hours to determine whether the original observation data are corresponding or not, performing integrity check on the original observation data acquired in S1, calculating satellite azimuth angle and altitude angle information of each epoch through a receiver approximate coordinate, summarizing satellite sight line information, and outputting continuous observation arc section observation values, wherein the initial observation value and the final observation value of the arc section observation values are marked as head and tail points;
s4: sunday environment modeling
Performing gross error elimination on the head and tail point observed values output by the S3, constructing an observation equation by using a plurality of head and tail point observed values, selecting a least square criterion to solve a fitting coefficient in the equation by a fitting system, and establishing an environment model according to the fitting coefficient;
the construction method of the environment model comprises the following steps:
s41: the formula (1) is utilized to construct an observation equation for a plurality of head and tail point observation values:
Figure GDA0002614093530000051
in the formula (1), n represents the fitting order, a i 、b i 、c i Denotes the fitting coefficient, e x Representing the measurement height angle, x representing e x A corresponding azimuth angle; n is 360/0.5, i is {0,0.5, …,360 };
s42: selecting a least square criterion to carry out fitting system solution on the observation equation to obtain a fitting coefficient, and establishing an environment model according to the fitting coefficient by taking 0.5 degrees as a step length from 0 to 360 degrees, as shown in a formula (2):
Figure GDA0002614093530000052
in equation (2), a 'model coefficient vector, e' i Representing the virtual satellite altitude using the model coefficients back-deducted, i ═ 0,0.5, …, 360.
S5: correlation coefficient solving
Setting intervals by taking days as units, and jointly solving correlation coefficients in the environment models according to the corresponding environment models every day;
the concrete method for solving the correlation coefficient comprises the following steps:
to set an environmental model a 'of day 1' 0 Reference, environmental model a 'for day k every 24 hours' k And a' 0 Is solved as shown in equation (3):
Figure GDA0002614093530000053
wherein, Cov (a' k ,a′ 0 ) Is a' k And a' 0 Covariance, Var [ a' k ]Is a' k Variance of (1), Var [ a' 0 ]Is a' 0 The variance of (c).
S6: early warning application
And setting a proper threshold value in practical application, comparing and judging the threshold value and the correlation coefficient calculated in the step S5, if the correlation coefficient is lower than the threshold value, carrying out early warning, and if not, continuing monitoring.
Example 2: in this embodiment, the method in embodiment 1 is used for actual measurement, data of a certain landslide monitoring point in a creep stage is selected for an experiment, the monitoring point gradually sinks along with the time, original observation data are collected and environment modeling is performed, the result is shown in fig. 2, the direction indicated by an arrow in fig. 2 represents the modeling result of 19 consecutive weeks from the 1 st week to the 19 th week, and it can be seen that the satellite height angle of the environment model in the same azimuth gradually increases along with the creep movement of the monitoring point.
In the time chart of the three-dimensional position change and the correlation coefficient of the monitoring point given in conjunction with fig. 3, the horizontal axis in fig. 3 represents the time axis, the left vertical axis represents the 3D vector displacement, and the right vertical axis represents the correlation coefficient. As can be seen from the displacement sequence, the monitoring point passes through a creep stage- > an acceleration stage- > and recovers the creep; the correlation coefficient is correspondingly restored to be stable from a stable stage- > an acceleration stage- >; therefore, it can be seen that the trend of the correlation coefficient can be consistent with the position change trend of the monitoring point, and therefore effective landslide early warning can be carried out through setting a reasonable threshold value.

Claims (6)

1. A GNSS environment model-based single-station landslide deformation monitoring and early warning method is characterized by comprising the following steps:
s1: deploying terminal equipment and performing data acquisition
The GNSS application terminal is deployed on a monitoring point, and the original observation data of the monitoring point is collected through a wireless communication network after the GNSS application terminal is started;
s2: determining the initial position of a monitoring point
After the original observation data are collected, downloading a real-time broadcast ephemeris file, and roughly estimating the rough position of the monitoring point by utilizing a pseudo-range single-point positioning technology;
s3: computing environment modeling source data
Based on the broadcast ephemeris, performing integrity check on the original observation data acquired by S1, calculating satellite azimuth angle and altitude angle information of each epoch through a receiver approximate coordinate, summarizing satellite sight line information, and outputting continuous observation arc-segment observation values, wherein the initial and final observation values of the arc-segment observation values are marked as head and tail points;
s4: sunday environment modeling
Performing gross error elimination on the head and tail point observed values output by the S3, constructing an observation equation by using a plurality of head and tail point observed values, selecting a least square criterion to solve a fitting coefficient in the equation by a fitting system, and establishing an environment model according to the fitting coefficient;
s5: correlation coefficient solving
Setting intervals by taking days as units, and jointly solving correlation coefficients in the environment models according to the corresponding environment models every day;
s6: early warning application
And setting a proper threshold value in practical application, comparing and judging the threshold value and the correlation coefficient calculated in the step S5, if the correlation coefficient is lower than the threshold value, carrying out early warning, and if not, continuing monitoring.
2. The GNSS environment model-based single-station landslide deformation monitoring and early warning method according to claim 1, wherein the raw observation data in S1 comprises raw pseudoranges and ephemeris data.
3. The GNSS environment model-based single-station landslide deformation monitoring and early warning method according to claim 1, wherein the GNSS application terminal in S1 is a consumer-grade single-frequency GNSS.
4. The GNSS environment model-based single-station landslide deformation monitoring and early warning method of claim 1, wherein the integrity detection method of the original observation data in S3 comprises: and comparing the raw observation data collected in the step S1 with the real-time broadcast ephemeris downloaded in the step S2 every 24 hours to determine whether the raw observation data correspond to the real-time broadcast ephemeris or not, so as to detect the integrity of the raw observation data.
5. The GNSS environment model-based single-station landslide deformation monitoring and early warning method according to claim 1, wherein the construction method of the environment model in S4 comprises the following steps:
s41: and (2) constructing an observation equation for a plurality of head and tail point observation values by using a formula (1):
Figure FDA0003667104720000021
in the formula (1), n represents the fitting order, a i 、b i 、c i Denotes the fitting coefficient, e x Representing the measurement height angle, x representing e x A corresponding azimuth angle; n is 360/0.5, i is {0,0.5, …,360 };
s42: selecting a least square criterion to carry out fitting system solution on the observation equation to obtain a fitting coefficient, and establishing an environment model according to the fitting coefficient by taking 0.5 degrees as a step length from 0 to 360 degrees, as shown in a formula (2):
Figure FDA0003667104720000022
in the formula (2), the first and second groups,
Figure FDA0003667104720000023
vector of model coefficients, e' i Representing the virtual satellite altitude using the model coefficients back-deducted, i ═ 0,0.5, …, 360.
6. The GNSS environment model-based single-station landslide deformation monitoring and early warning method as claimed in claim 1, wherein the specific method for solving the correlation coefficient in S5 is as follows:
to set an environmental model a 'of day 1' 0 Reference, environmental model a 'for day k every 24 hours' k And a' 0 Is solved as shown in equation (3):
Figure FDA0003667104720000024
wherein, Cov (a' k ,a′ 0 ) Is a' k And a' 0 Covariance, Var [ a' k ]Is a' k Variance of (1), Var [ a' 0 ]Is a' 0 The variance of (c).
CN202010460631.7A 2020-05-26 2020-05-26 GNSS environment model-based single-station landslide deformation monitoring and early warning method Active CN111694021B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010460631.7A CN111694021B (en) 2020-05-26 2020-05-26 GNSS environment model-based single-station landslide deformation monitoring and early warning method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010460631.7A CN111694021B (en) 2020-05-26 2020-05-26 GNSS environment model-based single-station landslide deformation monitoring and early warning method

Publications (2)

Publication Number Publication Date
CN111694021A CN111694021A (en) 2020-09-22
CN111694021B true CN111694021B (en) 2022-08-05

Family

ID=72478602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010460631.7A Active CN111694021B (en) 2020-05-26 2020-05-26 GNSS environment model-based single-station landslide deformation monitoring and early warning method

Country Status (1)

Country Link
CN (1) CN111694021B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114089055B (en) * 2021-09-30 2024-08-09 安徽继远软件有限公司 Method and system for monitoring safety state of power grid limited space operation personnel
CN115659497B (en) * 2022-10-18 2023-07-21 武汉理工大学 Ship body modeling method based on direct control free deformation

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8457882B2 (en) * 2011-04-19 2013-06-04 Honeywell International Inc. Systems and methods for navigation in a GPS-denied environment
CN104597471B (en) * 2015-01-13 2017-04-19 华东师范大学 Orientation attitude determination method oriented to clock synchronization multi-antenna GNSS receiver
CN106405589A (en) * 2016-06-24 2017-02-15 西安科技大学 Method and device for determining global ionized layer grid model
CN106500928B (en) * 2016-11-02 2018-09-04 济南大学 A kind of reservoir dam Seepage Safety online auto monitoring system
CN108332649B (en) * 2018-02-07 2020-04-24 桂林电子科技大学 Landslide deformation comprehensive early warning method and system
CN109977181A (en) * 2019-02-28 2019-07-05 广东赛诺科技股份有限公司 A method of based on bridge stereoscopic monitoring
CN111025357B (en) * 2019-12-31 2021-11-19 中国科学院国家授时中心 High-precision GNSS virtual SNR weighted positioning method

Also Published As

Publication number Publication date
CN111694021A (en) 2020-09-22

Similar Documents

Publication Publication Date Title
EP4206743A1 (en) Method for vehicle positioning, related apparatus, device, and storage medium
CN110515097B (en) GNSS satellite observation gross error rejection method and device applied to reference station
CN107765275B (en) Wide-area differential positioning method, device, terminal and computer readable storage medium
CN112327340B (en) Terminal positioning accuracy evaluation method, device, equipment and medium
CN108490459B (en) Method and system for applying precision and risk balance to GNSS position service
CN1670542A (en) Method and system for using altitude information in a satellite positioning system
CN104483691B (en) A kind of GNSS combines accurate one-point positioning method
CN111694021B (en) GNSS environment model-based single-station landslide deformation monitoring and early warning method
CN103900580A (en) Compass/GPS (global positioning system) and INS (inertial navigation system) combination vehicle navigation positioning system based on GIS (geographic information system) technology
CN103364803A (en) Satellite selection method and satellite navigation positioning method applying the satellite selection method
CN111290004A (en) Pseudo-range differential positioning method, pseudo-range differential positioning device, electronic equipment and storage medium
CN117739797B (en) Beidou/GNSS-based multi-time-scale deformation monitoring method
CN107807373A (en) GNSS high-precision locating methods based on mobile intelligent terminal
CN117826200A (en) PPP-B2B-based marine real-time precise positioning method, system and medium
CN115902968A (en) PPP terminal positioning method based on Beidou third GEO broadcast enhancement information
KR101617033B1 (en) Global standard point positioning apparatus using multi global satellite positioning systems and the method thereof
CN108491616B (en) Ionized layer vertical total electron content modeling method based on ellipsoid harmonic function theory
CN112230254B (en) Correction method and device for GPS carrier phase multipath error
Tarig Positioning with wide-area GNSS networks: Concept and application
Raghuvanshi et al. Precise positioning of smartphones using a robust adaptive Kalman filter
CN114355420B (en) PPP product positioning method and device for distributed Beidou position service center
Blomenhofer et al. Sensitivity analysis of the Galileo integrity performance dependent on the ground sensor station network
CN113866801A (en) Beidou satellite positioning accuracy evaluation improvement method and system based on vertical projection
CN114167469A (en) Vehicle navigation information monitoring method and device based on 5G/GNSS combination
KR100264977B1 (en) Method for the real-time estimation of user position in dgps

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