KR101738384B1 - System for integrity checking of measured location of DGNSS and method for integrity checking using for thereof - Google Patents

System for integrity checking of measured location of DGNSS and method for integrity checking using for thereof Download PDF

Info

Publication number
KR101738384B1
KR101738384B1 KR1020160162790A KR20160162790A KR101738384B1 KR 101738384 B1 KR101738384 B1 KR 101738384B1 KR 1020160162790 A KR1020160162790 A KR 1020160162790A KR 20160162790 A KR20160162790 A KR 20160162790A KR 101738384 B1 KR101738384 B1 KR 101738384B1
Authority
KR
South Korea
Prior art keywords
error
dgnss
module
correction
user
Prior art date
Application number
KR1020160162790A
Other languages
Korean (ko)
Inventor
김영기
서기열
황태현
박상현
Original Assignee
한국해양과학기술원
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 한국해양과학기술원 filed Critical 한국해양과학기술원
Priority to KR1020160162790A priority Critical patent/KR101738384B1/en
Application granted granted Critical
Publication of KR101738384B1 publication Critical patent/KR101738384B1/en

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/20Integrity monitoring, fault detection or fault isolation of space segment
    • 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/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/08Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing integrity information, e.g. health of satellites or quality of ephemeris data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position

Abstract

The present invention relates to an integrity checking system for a DGNSS measuring position and an integrity checking method using the same. The integrity checking system of a DGNSS measurement position according to an embodiment of the present invention includes a global navigation satellite system (GNSS) receiving module; an acquisition module for generating GNSS information and supplying the GNSS information to an error estimation unit and a correction positioning module; a correction information receiving module; a correction information acquisition module; a correction positioning module for supplying the position data of a satellite; an error estimation module composed of a pseudo range error estimation module, a non-common error estimation module and a correction value error estimation module; a generation module for calculating a protection level which a user wants; and an integrity checking module. Therefore, the present invention checks the integrity of the measurement position by only the DGNSS receiver.

Description

DGNSS 측정 위치의 무결성 검사 시스템 및 이를 이용한 무결성 검사방법{System for integrity checking of measured location of DGNSS and method for integrity checking using for thereof}TECHNICAL FIELD [0001] The present invention relates to an integrity checking system for a DGNSS measuring position and an integrity checking method using the same.

본 발명은 DGNSS 측정 위치의 무결성 검사 시스템 및 이를 이용한 무결성 검사방법에 관한 것으로서, DGNSS 이용자들이 수신기가 측정한 위치 값을 받아 항법에 활용할때 매번 측정되는 위치 값이 포함한 오차의 크기가 이용자가 원하는 수준보다 작은지 여부를 검사하고 그 결과를 제공함으로써 이용자의 안전한 항해를 도모할 수 있는 DGNSS 측정 위치의 무결성 검사 시스템 및 이를 이용한 무결성 검사방법에 관한 것이다.The present invention relates to an integrity checking system for a DGNSS measuring position and an integrity checking method using the same. More specifically, when the DGNSS users receive a position value measured by a receiver, The present invention relates to an integrity checking system for a DGNSS measurement position and an integrity checking method using the same.

일반적으로, 위성항법시스템(GNSS, Global Navigation Satellite System)은 위성을 이용한 전파항법시스템으로 미국의 GPS, 러시아의 GLONASS, 유럽연합의 Galileo, 중국의 BeiDou 등이 운영 중이다. In general, the Global Navigation Satellite System (GNSS) is a radio navigation system using satellites such as US GPS, GLONASS in Russia, Galileo in EU and BeiDou in China.

이러한 위성항법시스템 이용자는 위성항법시스템만을 이용한 경우보다 더욱 정확한 측정 결과를 획득하기 위하여 위성항법시스템과 함께 위성항법보정시스템을 이용하는데 이러한 위성항법보정시스템(DGNSS, Differential Global Navigation Satellite System)은 GNSS의 신호를 수신하고 이에 대한 보정정보를 생성하여 방송함으로써 사용자들이 보다 정확한 위치를 측정하도록 하는 시스템이다. The satellite navigation system (DGNSS, Differential Global Navigation Satellite System) is used by GNSS to obtain more accurate measurement results than satellite navigation system. A system for receiving a signal and generating correction information therefrom and broadcasting it so that users can measure a more accurate position.

이러한 육상/해양/항공 항법에 사용되는 DGNSS 수신기는 정확한 위치 값을 제공하지만 때때로 오차가 매우 큰 위치 값을 제공하는 경우도 있다.The DGNSS receivers used in these land / marine / air navigation systems provide accurate position values, but occasionally they provide very large position values.

만약 DGNSS 수신기가 오차가 큰 위치 값을 제공하면 사용자는 항로를 이탈하거나 장애물과의 관계해석에 실패하여 위험에 빠지게 될 수 있다.If the DGNSS receiver provides a positional value with a large error, the user may be in danger because he or she fails to interpret the relationship with the obstacle.

이와 관련하여, 국제해사기구(IMO, International Maritime Organization)는 미래 위성항법시스템에 대한 해양분야의 요구성능을 권고하는데, 그 중 하나로 위성항법시스템이 항법에 이용되면 안 될 경우에 경보를 제공하는 무결성 요건이 있다.In this regard, the International Maritime Organization (IMO) recommends the performance requirements of the maritime sector for future satellite navigation systems, one of which is the integrity of the satellite navigation system providing alarms There is a requirement.

또한, 국제민간항공기구(ICAO, International Civil Aviation Organization)도 항공기의 항법이나 이착륙 안전을 도모하기 위한 무결성 성능을 요구하고 있다.In addition, the International Civil Aviation Organization (ICAO) also calls for integrity performance to ensure aircraft navigation and landing and landing safety.

그러나 현재 DGNSS 사용자는 매번 측정된 위치 값이 얼마나 부정확한지 확인할 방법이 없는 문제가 있다.However, the present DGNSS user has a problem in that there is no way to check the accuracy of the measured position value each time.

선행특허는 GPS 위치오차 보정정보 제공시스템에서 IGS(International GPS Service)를 이용하여 GPS 위치오차 보정정보를 생성하여 제공하기 위한 GPS 위치오차 보정정보 제공방법 및 상기 GPS 위치오차 보정정보를 이용하여 GPS 수신기에서 GPS 위치오차를 보정하기 위한 GPS 위치오차 보정방법과 상기 방법들을 실현시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체를 제공하는 구성이 기재되어 있다.The prior art discloses a method for providing GPS position error correction information for generating and providing GPS position error correction information using IGS (International GPS Service) in a GPS position error correction information providing system and a method for providing GPS position error correction information using a GPS receiver A GPS position error correction method for correcting a GPS position error in a vehicle, and a computer-readable recording medium on which a program for realizing the methods are recorded.

그러나, 선행특허는 기준국의 보정정보를 이용한 보정 후에도 존재하는 잔여오차가 사용자 위치에 미치는 영향을 나타내는 것이 아닌 위치오차를 보정하기만 하는 것으로 지속적으로 DGNSS 사용자가 매번 측정된 위치 값이 얼마나 부정확한지 확인할 방법이 없는 기존의 문제점을 그대로 가지고 있는 문제가 있다.However, the prior patent does not show the influence of the residual error existing on the user position even after the correction using the correction information of the reference station, but merely corrects the position error, and it is difficult for the DGNSS user There is a problem that the existing problem which does not have a check method is inherited.

선행특허 : 한국 등록특허 공보 제10-0498185호Prior Patent: Korean Patent Registration No. 10-0498185

본 발명은 DGNSS 수신기가 때때로 제공하는 매우 큰 오차를 포함한 위치 측정치를 이용자가 사용하기 전에 제거하여 이용자가 원하는 수준의 품질을 가진 측정치만을 활용함으로써 이용자의 안전을 도모할 수 있도록 하는 DGNSS 측정 위치의 무결성 검사 시스템 및 이를 이용한 무결성 검사방법을 제공하는데 그 목적이 있다.The present invention eliminates position measurements, including very large errors that are sometimes provided by DGNSS receivers, before they are used by the user, thereby ensuring the safety of the user by utilizing only measurements having a desired level of quality. And an integrity checking method using the same.

상기 과제를 달성하기 위한 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템은, GNSS 위성군에서 송신되는 GNSS 신호를 수신하고 해석하여 GNSS 원시데이터를 제공하는 GNSS(위성항법 시스템: Global Navigation Satellite System) 수신모듈; 상기 GNSS 수신모듈에서 제공한 GNSS 원시데이터를 받아 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보를 획득하여 오차 추정부 및 보정측위모듈에 제공하는 획득모듈; DGNSS 기준국에서 송신되는 DGNSS 신호를 수신하고 해석하여 DGNSS 원시데이터를 제공하는 보정정보 수신모듈; 상기 보정정보 수신모듈에서 제공한 DGNSS 원시데이터를 받아 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 획득하여 오차 추정부 및 보정측위모듈에 제공하는 보정정보 획득모듈; 상기 획득모듈에서 제공한 GNSS 정보와 보정정보 획득모듈에서 제공한 DGNSS 정보를 이용하여 차분보정(Differential correction) 방법으로 사용자의 DGNSS 위치를 측정하고, 측정된 사용자의 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치를 오차 추정부 및 생성모듈에 제공하는 보정 측위모듈; 상기 획득모듈에서 제공한 GNSS 정보를 이용하여 추정한 위성별 의사거리 오차 추정치를 제공하는 의사거리 오차 추정모듈과, 보정정보 획득모듈에서 제공한 DGNSS 정보와 보정 측위모듈에서 제공한 사용자의 DGNSS 위치를 이용하여 추정한 위성별 비공통 오차 추정치를 제공하는 비공통 오차 추정모듈 및 보정정보 획득모듈에서 제공한 DGNSS 정보를 이용하여 추정한 위성별 보정치 오차 추정치를 제공하는 보정치 오차 추정모듈로 이루어지는 오차 추정부; 상기 보정 측위모듈 및 오차 추정부를 통해 제공되는 사용자의 DGNSS 위치, DGNSS 위치 측정에 사용된 각 위성의 위치, 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치, 위성별 보정치 오차 추정치를 이용하여 사용자가 원하는 수준의 보호수준(Protection Level)을 계산하여 제공하는 생성모듈 및 상기 생성모듈을 통해 제공되는 보호수준을 받아 사용자가 원하는 경보한계(Alert Limit)와 비교하여 사용자의 해당 DGNSS 위치 측정치를 항법에 활용할 것인지 여부를 결정하는 무결성 검사모듈을 포함하는 것을 특징으로 한다.According to another aspect of the present invention, there is provided a system for checking the integrity of a DGNSS measurement location, comprising: a GNSS (Global Navigation Satellite System) receiving and analyzing a GNSS signal transmitted from a GNSS satellite group and providing GNSS raw data; module; And receives GNSS raw data provided by the GNSS receiving module to obtain measurement time, receiver clock error, pseudo distance by each level, standard deviation of each pseudo distance, position of each satellite, and GNSS information, which is a clock error of each satellite An acquisition module for providing the error estimation unit and the corrected positioning module; A correction information receiving module for receiving and interpreting a DGNSS signal transmitted from a DGNSS reference station and providing DGNSS raw data; (DGNSS) data provided by the correction information receiving module to obtain correction value creation time (MZC, Modified-Z Count), position correction value, UDRE (User Differential Range Error) A correction information acquisition module provided to the government and the corrected positioning module; The DGNSS position of the user is measured by a differential correction method using the GNSS information provided by the acquisition module and the DGNSS information provided by the correction information acquisition module, and the DGNSS position of the user and the DGNSS position A correction positioning module for providing a position of each satellite to an error estimation unit and a generation module; A pseudo distance error estimation module for estimating a pseudorange distance error estimated value using the GNSS information provided by the acquisition module; a DGNSS information providing module for providing the correction information and a DGNSS position of the user provided by the correction positioning module; And a correction value error estimation module for providing a gross-gender-by-gender-by-gender-error estimation value using the DGNSS information provided by the correction information acquisition module, ; A DGNSS position, a position of each satellite used for DGNSS position measurement, a pseudorange error estimation value for each gender, a non-common error estimation value for each gender, and a gender-specific correction error estimation value provided through the corrected positioning module and error estimation unit, And a protection level provided through the generation module and compares the protection level with a desired alert limit to thereby measure the corresponding DGNSS position measurement value of the user. And an integrity checking module that determines whether or not to utilize the information.

구체적으로, 상기 보정정보 수신모듈(20)은 DGNSS 기준국에서 송신되는 DGNSS 신호를 중파, NTRIP, Eurofix, AIS(Auto Identification System) 또는 DMB(Digital Multimedia Broadcasting) 중 어느 하나 또는 다수를 통하여 수신받는 것을 특징으로 한다.Specifically, the correction information receiving module 20 receives a DGNSS signal transmitted from a DGNSS reference station through one or more of a medium wave, an NTRIP, a Eurofix, an AIS (Auto Identification System) or a DMB (Digital Multimedia Broadcasting) .

상기 과제를 달성하기 위한 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법은, DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법은, GNSS 수신모듈에서 GNSS 위성군의 GNSS 신호를 수신하여 해석한 GNSS 원시데이터 및 보정정보 수신모듈에서 DGNSS 기준국의 DGNSS 신호를 수신하여 해석한 DGNSS 원시데이터를 획득모듈과 보정정보 획득모듈에서 제공받는 단계; 상기 제공받는 단계에서 GNSS 수신모듈을 통해 GNSS 원시데이터를 제공받아 획득모듈에서 생성한 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보 및 상기 제공받는 단계에서 보정정보 수신모듈을 통해 DGNSS 원시데이터를 제공받아 보정정보 획득모듈에서 생성한 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 오차 추정부 및 보정측위모듈에 제공하는 단계; 상기 제공하는 단계를 통해 생성된 GNSS 정보와 DGNSS 정보를 이용하여 보정 측위모듈에서 차분보정(Differential correction) 방법으로 사용자의 DGNSS 위치를 측정하는 단계; 상기 사용자의 DGNSS 위치를 측정하는 단계를 통해 측정된 사용자의 DGNSS 위치, 상기 제공하는 단계를 통해 획득모듈에서 생성한 GNSS 정보 및 상기 보정정보 획득모듈에서 생성한 DGNSS 정보를 통해 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치 및 위성별 보정치 오차 추정치를 오차 추정부를 통해 각각 추정하는 단계; 상기 오차 추정부를 통해 각각 추정하는 단계에서 추정된 위성별 오차 추정치 및 보정 측위모듈에서 측정된 사용자의 DGNSS 위치를 통해 생성모듈에서 사용자가 원하는 수준의 보호수준을 계산하는 단계 및 상기 보호수준을 계산하는 단계를 통해 계산하여 생성한 보호수준을 사용자가 원하는 경보한계(Alert Limit)와 비교하여 해당 DGNSS 위치 측정치를 항법에 활용할 것인지 여부를 무결성 검사모듈을 통해 결정하는 단계를 포함하는 것을 특징으로 한다.
상기 오차 추정부를 통해 각각 추정하는 단계는, 상기 획득모듈에서 제공한 GNSS 정보를 이용하여 의사거리 오차 추정모듈을 통해 위성별 의사거리 오차를 추정하는 단계; 상기 보정정보 획득모듈에서 제공한 DGNSS 정보 및 보정 측위모듈에서 제공한 사용자의 DGNSS 위치를 이용하여 비공통 오차 추정모듈을 통해 기준국과 사용자 간의 거리를 구하고 이에 따른 위성별 비공통 오차를 추정하는 단계 및 상기 보정정보 획득모듈에서 제공한 DGNSS 정보를 이용하여 보정치 오차 추정모듈을 통해 위성별 보정치 오차를 추정하는 단계로 이루어지는 것을 특징으로 한다.
상기 의사거리 오차 추정모듈을 통해 위성별 의사거리 오차를 추정하는 단계는, 상기 획득모듈에서 오차 추정부 및 보정측위모듈에 전달하는 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정하며, 상기 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정할 때 사용할 것으로 결정되면, 관계식 1을 통해 각 위성별 의사거리의 표준편차(

Figure 112017036602278-pat00193
)를 취하여 위성별 의사거리 오차 추정치(
Figure 112017036602278-pat00194
)를 구하고, 상기 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정할 때 사용하지 않을 것으로 결정되면, 관계식 2를 통해 각 위성별 의사거리(
Figure 112017036602278-pat00195
)에서 각위성의 시계 오차(
Figure 112017036602278-pat00196
)와 수신기 시계 오차(
Figure 112017036602278-pat00197
)를 보상한 각 위성별 조정된 의사거리 오차(
Figure 112017036602278-pat00198
)를 구하고 관계식 3을 통해 각 위성별 조정된 의사거리 오차(
Figure 112017036602278-pat00199
)의 표준편차(
Figure 112017036602278-pat00200
)를 취하여 위성별 의사거리 오차 추정치(
Figure 112017036602278-pat00201
)를 구하는 것을 특징으로 한다.
상기 비공통 오차 추정모듈을 통해 위성별 비공통 오차를 추정하는 단계는, 상기 DGNSS 정보의 기준국 위치와 상기 보정 측위모듈(40)에서 측정된 사용자의 DGNSS 위치를 이용하여 기준국과 사용자 사이의 거리(
Figure 112017036602278-pat00202
)를 구하고, 이를 관계식 4을 통해 비공통 오차(
Figure 112017036602278-pat00203
)를 추정하는 것을 특징으로 한다.According to another aspect of the present invention, there is provided an integrity checking method using an integrity checking system for a DGNSS measurement location, the method comprising: receiving a GNSS signal of a GNSS satellite group from a GNSS receiving module; Receiving DGNSS primitive data obtained by analyzing the GNSS primitive data and correction information receiving module and receiving the DGNSS signal of the DGNSS reference station and receiving the DGNSS primitive data from the obtaining module and the correction information obtaining module; In the providing and receiving step, GNSS source data is received through the GNSS receiving module, and the measurement time, receiver clock error, standard deviation of each pseudo range, satellite position, satellite (MZC, Modified-Z Count) generated by the correction information acquiring module, DGNSS raw data through the correction information receiving module in the receiving step, (User Differential Range Error) and DGNSS information, which is a reference station position, to the error estimation unit and the correction positioning module; Measuring a DGNSS position of a user using a differential correction method in a correction positioning module using GNSS information and DGNSS information generated through the providing step; A DGNSS position of the user measured through the step of measuring the DGNSS position of the user, GNSS information generated by the acquisition module through the providing step, and DGNSS information generated by the correction information acquisition module, , A common non-common error estimate for each gender and an error estimate for each gender by an error estimator; Calculating a protection level of a desired level by a user in a generation module through a positional error estimation value estimated by the error estimating unit and a DGNSS position of a user measured by the correction positioning module, Comparing the generated protection level with an alert limit desired by the user and determining whether the DGNSS position measurement value is to be used for navigation through the integrity checking module.
Estimating the pseudorange error by the pseudorange error estimating module using the GNSS information provided by the acquiring module; Calculating a distance between the reference station and the user through the non-common error estimation module using the DGNSS information provided by the correction information acquisition module and the DGNSS position of the user provided by the correction positioning module, and estimating a non- And estimating a correction value for each sex by using a correction value error estimation module using the DGNSS information provided by the correction information acquisition module.
Wherein the step of estimating the pseudo distance error by the pseudo distance using the pseudo distance error estimation module determines whether or not to use the standard deviation of the pseudo distance by each degree transmitted to the error estimation unit and the correction positioning module in the acquisition module, wherein when each of the above when deciding whether to use the standard deviation of the sex pseudoranges determined to be used, the standard deviation of each of the above gender pseudorange through a relational expression 1 (
Figure 112017036602278-pat00193
) And the estimated pseudorange error
Figure 112017036602278-pat00194
) For obtaining, when each of the above decision that the gender is not used when determining whether or not to use the standard deviation of a distance, each of the above gender pseudorange through a second relational expression (
Figure 112017036602278-pat00195
), The clock error of each satellite
Figure 112017036602278-pat00196
) And receiver clock error (
Figure 112017036602278-pat00197
) Adjusted pseudorange error for each position compensated (
Figure 112017036602278-pat00198
) To obtain the above, each sex adjusted by a relational expression 3 pseudorange error (
Figure 112017036602278-pat00199
) Of the standard deviation
Figure 112017036602278-pat00200
) And the estimated pseudorange error
Figure 112017036602278-pat00201
) Is obtained.
The step of estimating non-common error for each gender through the non-common error estimation module may include estimating a non-common error between the reference station and the user using the reference station position of the DGNSS information and the DGNSS position of the user measured by the correction positioning module Street(
Figure 112017036602278-pat00202
), And this is calculated as a non-common error (
Figure 112017036602278-pat00203
).

삭제delete

삭제delete

삭제delete

상기 보정치 오차 추정모듈을 통해 위성별 보정치 오차를 추정하는 단계는, 상기 보정정보 획득모듈(35)에서 오차 추정부 및 보정측위모듈에 전달하는 각 위성별 UDRE를 사용할 것인지의 여부를 결정하며, 상기 각 위성별 UDRE를 사용할 것인지의 여부를 결정할 때 사용할 것으로 결정되면, 관계식 5을 통해 각 위성별 UDRE(

Figure 112017036602278-pat00204
)를 취하여 위성별 보정치 오차 추정치(
Figure 112017036602278-pat00205
)를 구하고, 상기 각 위성별 UDRE를 사용할 것인지의 여부를 결정할 때 사용하지 않을 것으로 결정되면, 관계식 6을 통해 각 위성별 보정치(
Figure 112017036602278-pat00206
)의 표준편차를 취하여 위성별 보정치 오차 추정치(
Figure 112017036602278-pat00207
)로 구하는 것을 특징으로 한다.
상기 사용자가 원하는 수준의 보호수준을 계산하는 단계는, 상기 오차 추정부를 통해 각각 추정하는 단계에서 추정되어 제공되는 DGNSS 정보를 통해 위성별 의사거리 오차 추정치(
Figure 112017036602278-pat00208
), 위성별 비공통 오차 추정치(
Figure 112017036602278-pat00209
) 및 위성별 보정치 오차 추정치(
Figure 112017036602278-pat00210
)를 식 1에 대입하여 오차 모델(
Figure 112017036602278-pat00211
)을 구성하는 단계; 상기 제공하는 단계의 획득모듈에서 획득한 각 위성의 위치 및 사용자의 DGNSS 위치를 측정하는 단계의 보정 측위모듈에서 측정된 사용자의 DGNSS 위치로 ENU 좌표계에서 사용자 위치로부터 위성의 방향에 대한 각 앙각(
Figure 112017036602278-pat00212
)과 방위각(
Figure 112017036602278-pat00213
)을 구하고, 식 2를 통해 위성 배치를 나타내는 단위벡터 행렬(
Figure 112017036602278-pat00214
)을 구성하는 단계; 상기 오차 모델(
Figure 112017036602278-pat00215
)를 통해 구성된 위성별 오차모델(
Figure 112017036602278-pat00216
) 값을 통해 오차모델 값의 역수값(
Figure 112017036602278-pat00217
)을 식 3으로 구하고, 상기 역수값(
Figure 112017036602278-pat00218
)을 식 4에 대입하여 식 4를 통해 가중치 행렬(
Figure 112017036602278-pat00219
)를 구하는 단계; 상기 단위벡터 행렬(
Figure 112017036602278-pat00220
)을 구성하는 단계를 통해 구해진 단위벡터 행렬(
Figure 112017036602278-pat00221
)과 가중치 행렬(
Figure 112017036602278-pat00222
)를 구하는 단계를 통해 구해진 가중치 행렬(
Figure 112017036602278-pat00223
)을 가중최소제곱법인 식 5에 대입하여 거리영역에서 위치영역으로 투영한 투영행렬(
Figure 112017036602278-pat00224
)을 구하는 단계; 상기 투영행렬(
Figure 112017036602278-pat00225
)을 구하는 단계를 통해 구해진 투영행렬(
Figure 112017036602278-pat00226
)의 일부 요소를 취하여 상기 오차 모델(
Figure 112017036602278-pat00227
)을 구성하는 단계에서 구한 위성별 오차 모델을 곱한 후 더하는 방법으로 식 6, 식 7, 식 8식 9를 통해 동쪽 방향으로의 분산 오차(
Figure 112017036602278-pat00228
), 북쪽 방향으로의 분산 오차(
Figure 112017036602278-pat00229
), 동쪽과 북쪽 방향으로의 공분산 오차(
Figure 112017036602278-pat00230
) 및 수직 방향으로의 분산 오차(
Figure 112017036602278-pat00231
)를 구하는 단계 및 상기 분산 오차를 구하는 단계를 통해 구해진 분산 오차 요소와 수평보호수준(
Figure 112017036602278-pat00232
, Horizontal Protection Level)에 대하여 사용자가 요구하는 무결성 리스크를 만족시키는 표준점수(
Figure 112017036602278-pat00233
)와 수직보호수준(
Figure 112017036602278-pat00234
, Vertical Protection Level)에 대하여 사용자가 요구하는 무결성 리스크를 만족시키는 표준점수(
Figure 112017036602278-pat00235
)를 이용하여 수평보호수준과 수직보호수준을 식 10식 11을 통해 구하는 단계를 포함하는 것을 특징으로 한다.The step of estimating the correction value error by the correction value error estimation module determines whether or not to use the UDRE of each position transmitted to the error estimation unit and the correction positioning module in the correction information acquisition module 35, If it is decided to use to determine whether to use each UDRE of each rank, the relation UDRE (
Figure 112017036602278-pat00204
) To obtain the correction value error estimate (
Figure 112017036602278-pat00205
) For obtaining, when each of the above is determined that is not used when determining whether or not to use the gender UDRE, each of the above gender correction value through equations 6 (
Figure 112017036602278-pat00206
), And calculate the error correction value
Figure 112017036602278-pat00207
). ≪ / RTI >
The step of calculating the level of protection desired by the user may further include the step of estimating the degree of protection by using the DGNSS information estimated by the estimating step through the error estimating unit,
Figure 112017036602278-pat00208
), Non-common error estimates by sex
Figure 112017036602278-pat00209
) And the gender-specific correction value error estimate
Figure 112017036602278-pat00210
) To the equation (1 ) to obtain the error model
Figure 112017036602278-pat00211
); The angle of the satellite from the user position in the ENU coordinate system to the DGNSS position of the user measured by the calibration positioning module in the step of measuring the position of each satellite acquired by the acquisition module of the providing step and the DGNSS position of the user,
Figure 112017036602278-pat00212
) And azimuth angle
Figure 112017036602278-pat00213
) To obtain a unit vector matrix representing the satellite deployment through a formula (2) (
Figure 112017036602278-pat00214
); The error model (
Figure 112017036602278-pat00215
), Which is composed of
Figure 112017036602278-pat00216
) Value, the reciprocal value of the error model value (
Figure 112017036602278-pat00217
) Is obtained by the equation (3 ), and the reciprocal value
Figure 112017036602278-pat00218
), The weight matrix through formula (4) are substituted in equation (4) (
Figure 112017036602278-pat00219
); The unit vector matrix (
Figure 112017036602278-pat00220
) Obtained by the step of constructing the unit vector matrix (
Figure 112017036602278-pat00221
) And a weighting matrix (
Figure 112017036602278-pat00222
) Obtained from the weight matrix
Figure 112017036602278-pat00223
) Is substituted into the weighted least squares equation (5 ), and a projection matrix (
Figure 112017036602278-pat00224
); The projection matrix (
Figure 112017036602278-pat00225
) Obtained by obtaining the projection matrix (
Figure 112017036602278-pat00226
) To obtain the error model (
Figure 112017036602278-pat00227
), And then multiplying the result by the error model obtained by the above step. Then, the variance error in the east direction is calculated by the equation 6, the equation 7, the equation 8 and the equation 9
Figure 112017036602278-pat00228
), The dispersion error in the north direction (
Figure 112017036602278-pat00229
), The covariance error in the east and north directions
Figure 112017036602278-pat00230
) And the dispersion error in the vertical direction (
Figure 112017036602278-pat00231
) And a variance error element obtained through the step of obtaining the dispersion error and a horizontal protection level
Figure 112017036602278-pat00232
, Horizontal Protection Level) that satisfies the user's required integrity risk.
Figure 112017036602278-pat00233
) And vertical protection level
Figure 112017036602278-pat00234
, Vertical Protection Level), which satisfies the user's required integrity risk.
Figure 112017036602278-pat00235
) To obtain the horizontal protection level and the vertical protection level through equations (10) and (11 ).

삭제delete

본 발명은 GNSS 위성신호와 DGNSS 보정신호를 수신하여 항해하는 사용자가 외부의 지원없이 DGNSS 수신기 자체만으로 측정 위치의 무결성을 검사할 수 있는 효과가 있다.The present invention has an effect that a user who receives and navigates a GNSS satellite signal and a DGNSS correction signal can check the integrity of a measurement location with only the DGNSS receiver itself without external assistance.

또한, 무결성 검사를 수행함에 있어서 사용자가 원하는 수준의 무결성 리스크를 충족시키는 보호수준을 계산할 수 있고, 사용자가 원하는 크기의 경보 한계와 비교할 수 있는 효과가 있다.Also, in performing the integrity check, the user can calculate the protection level satisfying the desired level of integrity risk, and the user can compare with the alarm limit of the desired size.

도 1은 기준국의 의사거리 보정치(PRC, Pseudo-Range Correction)를 적용하여 생성된 보정된 의사거리(CPR, Corrected Pseudo-Range)에 포함된 잔여오차가 최종적인 DGNSS 사용자 위치의 오차 원인으로 작용함을 나타내는 도면이다.
도 2는 GNSS 수신기를 사용하는 사용자의 항로와 GNSS 수신기가 제공한 위치값 사이에 거리 오차가 발생되는 것을 나타낸 도면이다.
도 3은 본 발명의 일 실시예에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 나타낸 도면이다.
도 4는 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법을 나타낸 순서도이다.
도 5는 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 오차를 추정하는 단계의 순서를 나탄낸 도면이다.
도 6은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 기준국과 사용자 수신기 사이의 거리에 따라 증가하는 오차를 나타낸 도면이다.
도 7은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 UDRE의 스케일 팩터(Scale Factor)를 나타낸 도면이다.
도 8은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 1시그마 오차에 대한 UDRE(User Differential Range Error) 수치를 나타낸 도면이다.
도 9는 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 보호수준을 결정하는 과정에 대한 순서도이다.
도 10은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법에 따른 실험 결과를 나타낸 그래프이다.
도 11은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법에 따른 실험결과에서 나타난 활용 횟수에 따른 가용성을 나타낸 도면이다.
FIG. 1 is a view illustrating a case where a residual error included in a corrected pseudo-range (CPR) generated by applying a pseudo-range correction (PRC) of a reference station is an error cause of a final DGNSS user position Fig.
2 is a diagram showing that a distance error is generated between the route of the user using the GNSS receiver and the position value provided by the GNSS receiver.
FIG. 3 is a diagram illustrating an integrity checking system for a DGNSS measurement location according to an exemplary embodiment of the present invention. Referring to FIG.
4 is a flowchart illustrating an integrity checking method using an integrity checking system for a DGNSS measurement location according to the present invention.
5 is a flowchart illustrating a procedure of estimating an error of an integrity checking method using an integrity checking system for a DGNSS measuring position according to the present invention.
FIG. 6 is a diagram illustrating an error that increases according to a distance between a reference station and a user receiver in an integrity checking method using an integrity checking system for a DGNSS measurement location according to the present invention.
FIG. 7 is a diagram illustrating a scale factor of a UDRE of an integrity checking method using an integrity checking system for a DGNSS measurement position according to the present invention.
FIG. 8 is a diagram showing a UDRE (User Differential Range Error) value for a 1-sigma error of an integrity checking method using an integrity checking system for a DGNSS measurement position according to the present invention.
FIG. 9 is a flowchart illustrating a procedure for determining a protection level of an integrity checking method using an integrity checking system for a DGNSS measurement location according to the present invention.
FIG. 10 is a graph showing an experimental result according to the integrity checking method using the integrity checking system of the DGNSS measuring position according to the present invention.
FIG. 11 is a diagram illustrating availability according to the number of times of use, which is an experiment result according to an integrity checking method using an integrity checking system for a DGNSS measurement position according to the present invention.

이하 본 발명의 실시를 위한 구체적인 실시예를 도면을 참고하여 설명한다. 본 발명의 실시예는 하나의 발명을 설명하기 위한 것으로서 권리범위는 예시된 실시예에 한정되지 아니하고, 예시된 도면은 발명의 명확성을 위하여 핵심적인 내용만 확대 도시하고 부수적인 것을 생략하였으므로 도면에 한정하여 해석하여서는 아니 된다.Hereinafter, embodiments of the present invention will be described with reference to the drawings. It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory and are not intended to limit the invention to the particular forms disclosed. And shall not interpret it.

도 3은 본 발명의 일 실시예에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 나타낸 도면으로서, GNSS 위성군에서 송신되는 GNSS 신호를 수신하고 해석하여 GNSS 원시데이터를 제공하는 GNSS(위성항법 시스템: Global Navigation Satellite System) 수신모듈(10); 상기 GNSS 수신모듈(10)에서 제공한 GNSS 원시데이터를 받아 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보를 획득하여 오차 추정부(50) 및 보정측위모듈(40)에 제공하는 획득모듈(30); DGNSS 기준국에서 송신되는 DGNSS 신호를 수신하고 해석하여 DGNSS 원시데이터를 제공하는 보정정보 수신모듈(20); 상기 보정정보 수신모듈(20)에서 제공한 DGNSS 원시데이터를 받아 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 획득하여 오차 추정부(50) 및 보정측위모듈(40)에 제공하는 보정정보 획득모듈(35); 상기 획득모듈(30)에서 제공한 GNSS 정보와 보정정보 획득모듈(35)에서 제공한 DGNSS 정보를 이용하여 차분보정(Differential correction) 방법으로 사용자의 DGNSS 위치를 측정하고, 측정된 사용자의 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치를 오차 추정부(50) 및 생성모듈(60)에 제공하는 보정 측위모듈(40); 상기 획득모듈(30)에서 제공한 GNSS 정보를 이용하여 추정한 위성별 의사거리 오차 추정치를 제공하는 의사거리 오차 추정모듈(51)과, 보정정보 획득모듈(35)에서 제공한 DGNSS 정보와 보정 측위모듈(40)에서 제공한 사용자의 DGNSS 위치를 이용하여 추정한 위성별 비공통 오차 추정치를 제공하는 비공통 오차 추정모듈(52) 및 보정정보 획득모듈(35)에서 제공한 DGNSS 정보를 이용하여 추정한 위성별 보정치 오차 추정치를 제공하는 보정치 오차 추정모듈(53)로 이루어지는 오차 추정부(50); 상기 보정 측위모듈(40) 및 오차 추정부(50)를 통해 제공되는 사용자의 DGNSS 위치, DGNSS 위치 측정에 사용된 각 위성의 위치, 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치, 위성별 보정치 오차 추정치를 이용하여 사용자가 원하는 수준의 보호수준(Protection Level)을 계산하여 제공하는 생성모듈(60) 및 상기 생성모듈(60)을 통해 제공되는 보호수준을 받아 사용자가 원하는 경보한계(Alert Limit)와 비교하여 사용자의 해당 DGNSS 위치 측정치를 항법에 활용할 것인지 여부를 결정하는 무결성 검사모듈(70)을 포함한다.FIG. 3 is a diagram illustrating a system for checking the integrity of a DGNSS measurement location according to an exemplary embodiment of the present invention. Referring to FIG. 3, a Global Navigation (GNSS) system, which receives and analyzes a GNSS signal transmitted from a GNSS satellite group and provides GNSS raw data, Satellite System) receiving module 10; And receives GNSS raw data provided by the GNSS receiving module 10, and receives the GNSS raw data provided by the GNSS receiving module 10, such as a measurement time, a receiver clock error, a pseudorange by each level, a standard deviation of pseudoranges by each level, An acquisition module (30) for acquiring and providing to the error estimation unit (50) and the correction positioning module (40); A correction information receiving module (20) for receiving and interpreting a DGNSS signal transmitted from a DGNSS reference station and providing DGNSS raw data; And receives DGNSS raw data provided by the correction information receiving module 20 to obtain correction value generation time (MZC, Modified-Z Count), statistical correction value, UDRE (User Differential Range Error) and DGNSS information A correction information acquisition module (35) for providing the correction information to the error estimation unit (50) and the correction positioning module (40); The DGNSS position of the user is measured by a differential correction method using the GNSS information provided by the acquisition module 30 and the DGNSS information provided by the correction information acquisition module 35, A correction positioning module 40 for providing the position of each satellite used for the DGNSS position measurement to the error estimation unit 50 and the generation module 60; A pseudo distance error estimation module 51 for providing a pseudorange distance error estimation value by using the GNSS information provided by the acquisition module 30, a DGNSS information providing module Common error estimation module 52 for providing a gender-specific non-common error estimation value estimated using the DGNSS position of the user provided by the module 40 and a DGNSS information provided by the correction information acquisition module 35, An error estimator (50) comprising a correction value error estimation module (53) for providing a correction value error estimate for the upper gender; A DGNSS position of the user provided through the correction positioning module 40 and the error estimation unit 50, a position of each satellite used for the DGNSS position measurement, a pseudorange error estimate by sex, a non-common error estimate by sex, A generation module 60 for calculating and providing a protection level of a desired level by using a correction value error estimation value and a protection level provided through the generation module 60 and receiving an alert limit And an integrity checking module 70 for determining whether to use the user's corresponding DGNSS position measurements for navigation.

도 3을 참조하면, 상기 GNSS(위성항법 시스템: Global Navigation Satellite System) 수신모듈(10)은 GNSS 위성군에서 송신되는 GNSS 신호를 수신하고 해석하여 GNSS 원시데이터를 제공한다.Referring to FIG. 3, the Global Navigation Satellite System (GNSS) receiving module 10 receives and interprets a GNSS signal transmitted from a GNSS satellite group and provides GNSS raw data.

상기 보정정보 수신모듈(20)은 DGNSS(Differential Global Navigation Satellite System) 기준국에서 송신되는 DGNSS 신호를 수신하고 해석하여 DGNSS 원시데이터를 제공한다.The correction information receiving module 20 receives and analyzes a DGNSS signal transmitted from a DGNSS (Differential Global Navigation Satellite System) reference station and provides DGNSS source data.

상기 보정정보 수신모듈(20)은 기준국에서 전송되는 DGNSS 보정정보를 중파, NTRIP, Eurofix, AIS(Auto Identification System) 또는 DMB(Digital Multimedia Broadcasting) 중 어느 하나 또는 다수를 통하여 DGNSS 보정정보를 수신하게 된다.The correction information receiving module 20 receives the DGNSS correction information through one or more of the medium wave, NTRIP, Eurofix, AIS (Auto Identification System), or DMB (Digital Multimedia Broadcasting) do.

상기 획득모듈(30)은 상기 GNSS 수신모듈(10)에서 제공한 GNSS 원시데이터를 받아 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보를 생성하여 보정 측위모듈(40) 및 오차 추정부(50)로 제공한다.The acquisition module 30 receives the GNSS raw data provided by the GNSS reception module 10, and receives the measurement time, the receiver clock error, the pseudorange by each level, the standard deviation of the pseudoranges by each level, Generates GNSS information, which is a clock error of the satellite, and provides it to the correction positioning module 40 and the error estimation unit 50. [

상기 획득모듈(30)에서 GNSS 정보를 생성할 때 동시에 보정정보 획득모듈(35)에서 DGNSS 정보를 생성한다.At the same time when generating the GNSS information in the acquisition module 30, the correction information acquisition module 35 generates DGNSS information.

상기 보정정보 획득모듈(35)은 보정정보 수신모듈(20)에서 제공한 DGNSS 원시데이터를 받아 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 생성하여 보정 측위모듈(40) 및 오차 추정부(50)로 제공한다.The correction information acquisition module 35 receives the DGNSS raw data provided by the correction information reception module 20 and receives the correction value generation time (MZC, Modified-Z Count), the statistical correction value, the UDRE (User Differential Range Error) And provides the DGNSS information to the correction positioning module 40 and the error estimating unit 50. [

상기 획득모듈(30)에서 제공받은 GNSS 정보와 보정정보 획득모듈(35)에서 제공받은 DGNSS 정보를 이용하여 보정 측위모듈(40)에서 차분보정(Differential correction) 방법으로 사용자의 DGNSS 위치를 측정하고, 측정된 사용자의 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치를 오차 추정부(50) 및 생성모듈(50)에 제공한다.The DGNSS position of the user is measured by the correction positioning module 40 using the differential correction method using the GNSS information provided by the acquisition module 30 and the DGNSS information provided by the correction information acquisition module 35, And provides the error estimation unit 50 and the generation module 50 with the DGNSS position of the measured user and the position of each satellite used for the DGNSS position measurement.

상기 차분보정(Differential correction) 방법은 위성항법시스템을 이용한 상대측위 방법의 하나인 단일차분(Single difference) 방법 중 수신기간 단일차분에 해당하는 방법으로 DGNSS 기준국에서 생성한 보정정보를 이용하여 사용자 수신기의 거리측정치를 보정함으로써 사용자 측정 위치의 정확도를 향상시키는 방법이다.The differential correction method is a method corresponding to a single difference of reception period among single differential methods, which is one of relative positioning methods using a satellite navigation system, using a correction information generated by a DGNSS reference station, Thereby improving the accuracy of the user measurement position.

상기 보정 측위모듈(40)에서 측정된 사용자의 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치가 제공되면, 오차 추정부(50)에서 획득모듈(30)에서 제공한 GNSS 정보, 보정정보 획득모듈(35)에서 제공한 DGNSS 정보 및 보정 측위모듈(40)을 통해 측정된 사용자의 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치를 이용하여 기준국과 사용자 간의 거리를 통해 의사거리 오차 추정치, 비공통 오차 추정치 및 보정치 오차 추정치를 추정한다.When the DGNSS position of the user measured by the corrected positioning module 40 and the position of each satellite used for the DGNSS position measurement are provided, the error estimator 50 obtains the GNSS information provided by the acquisition module 30, A pseudorange error estimate is obtained through the distance between the reference station and the user using the DGNSS information provided by the module 35 and the DGNSS position of the user measured through the correction positioning module 40 and the position of each satellite used for the DGNSS position measurement , A non-common error estimate and a correction value error estimate.

상기 오차 추정부(50)는 획득모듈(30)에서 제공한 GNSS 정보를 이용하여 추정한 위성별 의사거리 오차 추정치를 제공하는 의사거리 오차 추정모듈(51), 보정정보 획득모듈(35)에서 제공한 DGNSS 정보와 보정 측위모듈(40)에서 제공한 사용자의 DGNSS 위치를 이용하여 추정한 위성별 비공통 오차 추정치를 제공하는 비공통 오차 추정모듈(52) 및 보정정보 획득모듈(35)에서 제공한 DGNSS 정보를 이용하여 추정한 위성별 보정치 오차 추정치를 제공하는 보정치 오차 추정모듈(53)로 이루어진다.The error estimator 50 is provided by a pseudorange error estimation module 51 and a correction information acquisition module 35 that provide estimates of pseudorange error distances by using the GNSS information provided by the acquisition module 30 Common error estimation module 52 and correction information acquisition module 35 that provide a non-common error estimation value for each gender estimated using the DGNSS information and the DGNSS position of the user provided by the correction positioning module 40 And a correction value error estimation module 53 for providing a gender-specific correction value error estimation value estimated using the DGNSS information.

상기 의사거리 측정치 오차 추정모듈(51)은 획득모듈(30)이 제공한 측정 시각, 수신기 시계 오차, 각 위성에 대한 의사거리, 각 위성에 대한 의사거리의 표준편차, 각 위성의 시계 오차를 획득하고 이를 이용하여 각 위성에 대한 의사거리에 포함된 1 시그마(

Figure 112016118091315-pat00042
) 오차를 추정하여 추정된 값을 위성별 의사거리 오차 추정치로 취한다.The pseudo distance measurement error estimation module 51 acquires the measurement time provided by the acquisition module 30, the receiver clock error, the pseudorange for each satellite, the standard deviation of the pseudorange for each satellite, (1 sigma) included in pseudoranges for each satellite
Figure 112016118091315-pat00042
) Is taken as the estimated value of the pseudorange error for each sex.

상기 1시그마(

Figure 112016118091315-pat00043
) 오차는 분산 형태의 오차를 말하는 것으로, 일반적으로 거리 또는 위치의 오차를 표현할때 사용되는 용어이다.The 1-sigma (
Figure 112016118091315-pat00043
) Error refers to the variance type error, and is generally used to describe the error of distance or position.

상기 획득모듈(30)로부터 각 위성에 대한 의사거리의 표준편차를 획득하였다면, 이 값을 위성별 의사거리 오차 추정치로 사용한다.If the standard deviation of the pseudorange for each satellite is obtained from the acquisition module 30, this value is used as the pseudorange error estimate for each satellite.

상기 획득모듈(30)로부터 각 위성에 대한 의사거리의 표준편차를 획득하지 못하였다면, 각 위성에 대한 의사거리에서 수신기 시계 오차와 각 위성의 시계 오차에 해당하는 거리 오차를 보상하여 각 위성별 조정된 의사거리를 구하고 측정 시각을 이용하여 이동 표준편차 방법으로 각 위성별 조정된 의사거리의 표준편차를 구한다. 구해진 각 위성별 조정된 의사거리의 표준편차는 각 위성별 조정된 의사거리 오차의 표준편차가 된다. 각 위성별 조정된 의사거리 오차의 표준편차를 위성별 의사거리 오차 추정치로 정하여 제공한다.If the standard deviation of the pseudorange for each satellite is not obtained from the acquisition module 30, the distance error corresponding to the receiver clock error and the clock error of each satellite in the pseudo range for each satellite is compensated for, And the standard deviation of the pseudo distance adjusted for each position is obtained by the moving standard deviation method using the measurement time. The standard deviation of the adjusted pseudoranges for each position obtained is the standard deviation of the adjusted pseudorange errors for each position. The standard deviation of the adjusted pseudorange error for each position is provided as the estimated value of the pseudorange error.

상기 비공통 오차 추정모듈(52)은 보정정보 획득모듈(35)에서 제공한 DGNSS 정보와 보정 측위모듈(40)에서 제공한 사용자의 DGNSS 위치를 이용하여 기준국과 사용자 간의 거리 값을 구하고, 이 거리 값에

Figure 112016118091315-pat00044
을 곱하여 구한 값을 비공통 오차 추정치로 정하고 모든 위성에 대해 동일한 값을 적용하여 위성별 비공통 오차의 추정치로 취하여 제공한다.The non-common error estimation module 52 obtains the distance value between the reference station and the user by using the DGNSS information provided by the correction information acquisition module 35 and the DGNSS position of the user provided by the correction positioning module 40, To the distance value
Figure 112016118091315-pat00044
And the same value is applied to all satellites, and it is taken as an estimate of the common error of each sex.

이때, 비공통 오차는 기준국과 사용자 수신기 사이의 이격이 증가함에 따라 동일 위성이 송신한 신호가 이온층/대류층을 통과하면서 지연되는 량의 상관도가 하락하여 발생하는 오차이며 이로 인하여 기준국과 사용자 사이의 거리가 증가함에 따라 증가하는 거리 측정치의 1 시그마 오차는 도 6을 참고하여 기준국과 사용자 사이의 거리 1

Figure 112016118091315-pat00045
당 각 의사거리 측정치에 최대 약 0.57
Figure 112016118091315-pat00046
만큼의 1 시그마 오차가 발생하여 DGNSS 측위 결과에 영향을 미치게 된다.In this case, the non-common error is an error caused by a decrease in the degree of correlation of the amount of delay that the signal transmitted from the same satellite passes through the ionosphere / convection layer as the distance between the reference station and the user receiver increases, 1-sigma error of the range measurement to increase as the distance is increased between the user reference to FIG. 6, the distance between the reference station and the user 1
Figure 112016118091315-pat00045
Up to about 0.57 per pseudorange
Figure 112016118091315-pat00046
1 sigma error occurs, which affects the DGNSS positioning result.

상기 보정치 오차 추정모듈(53)은 보정정보 획득모듈(35)에서 제공한 DGNSS 정보를 이용하여 위성별 UDRE(User Differential Range Error)를 해석하여 보정치에 들어있는 1시그마(

Figure 112016118091315-pat00047
) 오차를 추정한 위성별 보정치 오차 추정치를 구한다.The correction value error estimation module 53 interprets the UDRE (User Differential Range Error) by using the DGNSS information provided by the correction information acquisition module 35, and outputs the 1-sigma
Figure 112016118091315-pat00047
) Error is estimated.

만약 상기 위성별 UDRE가 제공되지 않았다면 보정치 생성 시각과 위성별 보정치를 이용하여 위성별 보정치의 표준편차를 구하고 이 값을 위성별 보정치 오차 추정치로 정하여 제공한다.If the above UDRE per URI is not provided, the standard deviation of the correction value for each sex is obtained using the correction value generation time and the correction value for each sex, and this value is provided as the correction value error estimation value for each sex.

상기 생성모듈(60)은 보정 측위모듈(40) 및 오차 추정부(50)를 통해 제공되는 사용자의 DGNSS 위치, DGNSS 위치 측정에 사용된 각 위성의 위치, 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치, 위성별 보정치 오차 추정치를 이용하여 사용자가 원하는 수준의 보호수준(Protection Level)을 계산하여 제공한다.The generation module 60 calculates the DGNSS position of the user provided through the correction positioning module 40 and the error estimation unit 50, the position of each satellite used for the DGNSS position measurement, the estimated value of the pseudo- A common error estimate, and a correction value error estimate for each gender, and provides the calculated protection level.

상기 생성모듈(60)은 보정 측위모듈(40)이 제공한 사용자 DGNSS 위치, DGNSS 위치 측정에 사용된 각 위성의 위치, 의사거리 오차 추정모듈(51)이 제공한 위성별 의사거리 오차 추정치, 비공통 오차 추정모듈(52)이 제공한 위성별 비공통 오차 추정치 및 보정치 오차 추정모듈(53)이 제공한 보정치 오차 추정치를 이용하여 사용자가 원하는 수준의 보호수준을 계산한다.The generation module 60 receives the user DGNSS position provided by the correction positioning module 40, the position of each satellite used for the DGNSS position measurement, the estimated value of the pseudo distance error provided by the pseudo range error estimation module 51, Common non-common error estimates provided by the common error estimation module 52 and the correction value error estimates provided by the correction value error estimation module 53. [

상기 생성모듈(60)을 통해 계산된 사용자가 원하는 수준의 보호수준을 무결성 검사모듈(70)에서 생성모듈(60)을 통해 제공되는 보호수준을 받아 사용자가 원하는 경보한계(Alert Limit)와 비교하여 사용자의 해당 DGNSS 위치 측정치를 항법에 활용할 것인지 여부를 결정하게 된다.The protection level provided by the user through the generation module 60 is provided to the integrity checking module 70 through the generation module 60 and is compared with an alert limit desired by the user And determines whether to use the user's corresponding DGNSS position measurements for navigation.

이와 같이 구성된 본 발명의 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법은 다음과 같다.The integrity checking method using the integrity checking system of the DGNSS measuring position according to the present invention configured as described above is as follows.

도 4는 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법을 나타낸 순서도이다.4 is a flowchart illustrating an integrity checking method using an integrity checking system for a DGNSS measurement location according to the present invention.

도 4를 참조하면, DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법은, GNSS 수신모듈(10)에서 GNSS 위성군의 GNSS 신호를 수신하여 해석한 GNSS 원시데이터 및 보정정보 수신모듈(20)에서 DGNSS 기준국의 DGNSS 신호를 수신하여 해석한 DGNSS 원시데이터를 획득모듈(30)과 보정정보 획득모듈(35)에서 제공받는 단계(S10); 상기 제공받는 단계(S10)에서 GNSS 수신모듈(10)을 통해 GNSS 원시데이터를 제공받아 획득모듈(30)에서 생성한 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보 및 상기 제공받는 단계(S10)에서 보정정보 수신모듈(20)을 통해 DGNSS 원시데이터를 제공받아 보정정보 획득모듈(35)에서 생성한 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 오차 추정부(50) 및 보정측위모듈(40)에 제공하는 단계(S20); 상기 제공하는 단계(S20)를 통해 생성된 GNSS 정보와 DGNSS 정보를 이용하여 보정 측위모듈(40)에서 차분보정(Differential correction) 방법으로 사용자의 DGNSS 위치를 측정하는 단계(S30); 상기 사용자의 DGNSS 위치를 측정하는 단계(S30)를 통해 측정된 사용자의 DGNSS 위치, 상기 제공하는 단계(S20)를 통해 획득모듈(30)에서 생성한 GNSS 정보 및 상기 보정정보 획득모듈(35)에서 생성한 DGNSS 정보를 통해 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치 및 위성별 보정치 오차 추정치를 오차 추정부(50)를 통해 각각 추정하는 단계(S40); 상기 오차 추정부(50)를 통해 각각 추정하는 단계(S40)에서 추정된 위성별 오차 추정치 및 보정 측위모듈(40)에서 측정된 사용자의 DGNSS 위치를 통해 생성모듈(60)에서 사용자가 원하는 수준의 보호수준을 계산하는 단계(S50) 및 상기 보호수준을 계산하는 단계(S50)를 통해 계산하여 생성한 보호수준을 사용자가 원하는 경보한계(Alert Limit)와 비교하여 해당 DGNSS 위치 측정치를 항법에 활용할 것인지 여부를 무결성 검사모듈(70)을 통해 결정하는 단계(S60)를 포함한다.Referring to FIG. 4, the integrity checking method using the integrity checking system of the DGNSS measuring position is performed by receiving the GNSS signal of the GNSS satellite group from the GNSS receiving module 10 and analyzing the GNSS source data and the correction information receiving module 20 (S10) receiving DGNSS raw data, which is received and analyzed by the DGNSS signal of the DGNSS reference station, from the acquisition module (30) and the correction information acquisition module (35); In the receiving step S10, GNSS source data is received through the GNSS receiving module 10, and the measurement time generated by the acquisition module 30, the receiver clock error, the pseudo distance by each level, The GNSS information which is the clock error of each satellite and the DGNSS raw data through the correction information receiving module 20 in the receiving step S10 and received by the correction information obtaining module 35, Providing the error estimation unit 50 and the correction positioning module 40 with the generation time (MZC, Modified-Z Count), the statistical correction value, the UDRE (User Differential Range Error) S20); A step (S30) of measuring a DGNSS position of a user by a differential correction method in the correction positioning module 40 using the GNSS information and the DGNSS information generated through the providing step S20; The DGNSS position of the user measured through the DGNSS position measurement of the user (S30), the GNSS information generated by the acquisition module 30 through the providing step S20, and the GNSS information generated by the correction information acquisition module 35 (S40) estimating the pseudo distance error estimate, the gender-by-common non-common error estimation value, and the gender-specific correction value error estimation value through the error estimation unit 50 through the generated DGNSS information, respectively; The generation module 60 estimates the user's desired level of error by using the estimated error by the estimated error in step S40 and the DGNSS position of the user measured by the corrected positioning module 40 through the error estimator 50, (S50) calculating the protection level and calculating the protection level (S50), compares the generated protection level with the desired alert limit (Alert Limit), and determines whether to use the corresponding DGNSS position measurement for navigation (Step S60) through the integrity checking module 70. [0050]

상기 제공받는 단계(S10)는 GNSS 위성군에서 송신한 다수 위성의 신호를 GNSS 수신모듈(10)에서 수신받아 해석하여 제공한 GNSS 원시데이터를 획득모듈(20)에서 제공받고 DGNSS 기준국에서 송신한 DGNSS 신호를 보정정보 수신모듈(20)에서 수신받아 해석하여 제공한 DGNSS 원시데이터를 보정정보 획득모듈(35)에서 제공받는 단계이다.In step S10, the GNSS reception module 10 receives the GNSS satellite group signal from the GNSS satellite group. The GNSS source data is received from the acquisition module 20 and transmitted from the DGNSS reference station. The DGNSS signal is received by the correction information receiving module 20 and the DGNSS raw data provided by the analysis is supplied from the correction information obtaining module 35. [

상기 제공하는 단계(S20)는 제공받는 단계(S10)에서 GNSS 수신모듈(10)을 통해 GNSS 원시데이터를 제공받은 획득모듈(30)에서 생성한 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보 및 제공받는 단계(S10)에서 보정정보 수신모듈(20)을 통해 DGNSS 원시데이터를 제공받은 보정정보 획득모듈(35)에서 생성한 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 보정 측위모듈(40) 및 오차 추정부(50)로 제공하는 단계(S20)이다.The providing step S20 may include measuring time, receiver clock error, pseudo distance according to each level, generated by the acquisition module 30 provided with the GNSS raw data through the GNSS receiving module 10 in step S10, The GNSS information, which is a clock error of each satellite, and the correction information acquisition module (DGNSS), which are provided with the DGNSS raw data through the correction information receiving module 20 in the receiving step (S10) (DGNSS), which is a reference station position, to the correction positioning module 40 and the error estimator 50 (Step S20).

상기 보정정보 수신모듈(20)은 기준국에서 전송되는 DGNSS 보정정보를 중파, NTRIP, Eurofix, AIS(Auto Identification System) 또는 DMB(Digital Multimedia Broadcasting)를 통하여 DGNSS 보정정보를 수신하게 된다.The correction information receiving module 20 receives the DGNSS correction information through the medium wave, the NTRIP, the Eurofix, the AIS (Auto Identification System), or the DMB (Digital Multimedia Broadcasting) through the DGNSS correction information transmitted from the reference station.

상기 보정정보 수신모듈(20)이 중파, NTRIP, Eurofix, AIS(Auto Identification System) 또는 DMB(Digital Multimedia Broadcasting) 중 어느 하나 또는 다수를 통해 DGNSS 보정정보를 수신하는 이유는 다양한 매체를 통하여 DGNSS 보정정보를 수신받는 DGNSS 사용자를 고려하기 위함이다.The reason why the correction information receiving module 20 receives the DGNSS correction information through one or more of medium frequency, NTRIP, Eurofix, AIS (Auto Identification System) or DMB (Digital Multimedia Broadcasting) To the DGNSS user.

상기 제공받는 단계에서 GNSS 수신모듈(10)을 통해 제공된 GNSS 원시데이터 및 상기 제공받는 단계에서 보정정보 수신모듈(20)을 통해 제공된 DGNSS 원시데이터를 획득모듈(30)과 보정정보 획득모듈(35)에서 제공받으면, 제공하는 단계(S20)를 통해 획득모듈(30)에서 GNSS 원시데이터를 이용하여 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보를 생성하고 보정정보 획득모듈(35)에서 DGNSS 원시데이터를 이용하여 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 생성하고 상기의 GNSS 정보와 DGNSS 정보를 하위에 연결되는 오차 추정부(50) 및 보정 측위모듈(40)에 제공한다.The GNSS raw data provided through the GNSS receiving module 10 and the DGNSS raw data provided through the correction information receiving module 20 in the providing and receiving step are supplied to the obtaining module 30 and the correction information obtaining module 35, The receiver clock error, the pseudorange by each level, the standard deviation of the pseudoranges by each level, and the position of each satellite using the GNSS raw data in the providing module 30 through the providing step S20 , GNSS information which is a clock error of each satellite, and the correction information acquisition module 35 uses the DGNSS raw data to generate a correction value (MZC, Modified-Z Count), a gender-specific correction value, a user differential range error And the DGNSS information, and supplies the GNSS information and the DGNSS information to the error estimation unit 50 and the correction positioning module 40 which are connected to each other.

상기 제공하는 단계(S20)를 통해 제공되는 GNSS 정보와 DGNSS 정보를 이용하여 보정측위 모듈(40)에서 차분보정(Differential correction) 방법으로 사용자의 DGNSS 위치를 측정하고 측정된 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치를 제공한다.The DGNSS position of the user is measured by the differential correction method in the correction positioning module 40 using the GNSS information and the DGNSS information provided through the providing step S20 and the measured DGNSS position and the DGNSS position measurement Provide the location of each satellite used.

도 5는 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 오차를 추정하는 단계(S40)의 순서를 나타낸 도면이다.FIG. 5 is a flowchart illustrating a procedure of estimating an error of an integrity checking method using the integrity checking system of the DGNSS measurement position (S40) according to the present invention.

도 5를 참조하면, 상기 오차 추정부(50)를 통해 각각 추정하는 단계(S40)는 상기 사용자의 DGNSS 위치를 측정하고 측정된 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치를 제공하는 단계(S30)를 통해 측정된 사용자의 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치, 상기 제공하는 단계(S20)를 통해 획득 모듈(30)에서 생성한 GNSS 정보 및 상기 보정정보 획득 모듈(35)에서 생성한 DGNSS 정보를 통해 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치 및 위성별 보정치 오차 추정치를 오차 추정부(50)를 통해 각각 추정하게 된다.Referring to FIG. 5, the estimating step S40 through the error estimating unit 50 may include measuring the DGNSS position of the user and providing the measured DGNSS position and the position of each satellite used in the DGNSS position measurement, The position of each satellite used for DGNSS position measurement, the GNSS information generated by the acquisition module 30 through the providing step S20, and the correction information acquisition module 35 The error estimator 50 estimates the pseudorange error estimate, the gender-specific non-common error estimate, and the gender-specific correction error estimation value, respectively, using the DGNSS information generated by the DGNSS.

상기 오차 추정부(50)를 통해 각각 추정하는 단계(S40)를 통해 각각 추정되는 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치 및 위성별 보정치 오차 추정치는 오차 추정부(50)의 의사거리 오차 추정모듈(51), 비공통 오차 추정모듈(52) 및 보정치 오차 추정모듈(53)을 통해 각각 추정하게 된다.The gender-specific non-common error estimates and the gender-specific correction error estimates estimated through the estimating step S40 through the error estimating unit 50 are calculated by the error estimating unit 50, Common error estimation module 52, and correction value error estimation module 53, respectively.

상기 오차 추정부(50)를 통해 각각 추정하는 단계(S40)는, 획득 모듈(30)에서 제공한 GNSS 정보를 이용하여 의사거리 오차 추정 모듈(51)을 통해 위성별 의사거리 오차를 추정하는 단계(S41), 보정정보 획득 모듈(35)에서 제공한 DGNSS 정보 및 보정측위 모듈(40)에서 제공한 사용자의 DGNSS 위치를 이용하여 비공통 오차 추정모듈(52)을 통해 기준국과 사용자 간의 거리를 구하고 이에 따른 위성별 비공통 오차를 추정하는 단계(S42) 및 보정정보 획득 모듈(35)에서 제공한 DGNSS 정보를 이용하여 보정치 오차 추정 모듈(53)을 통해 위성별 보정치 오차를 추정하는 단계(S43)로 이루어진다.The estimating step S40 through the error estimating unit 50 may include estimating the pseudorange error for each position through the pseudo range error estimating module 51 using the GNSS information provided by the acquiring module 30 (S41). The DGNSS information provided by the correction information acquisition module 35 and the DGNSS position of the user provided by the correction positioning module 40 are used to calculate the distance between the reference station and the user through the non-common error estimation module 52 (S42) of estimating the gender-specific uncorrelated error according to the error, and estimating the gender-specific correction value error using the correction value error estimation module 53 using the DGNSS information provided by the correction information acquisition module 35 (S43 ).

상기 의사거리 오차 추정모듈(51)을 통해 위성별 의사거리 오차를 추정하는 단계(S40), 비공통 오차 추정모듈(52)을 통해 위성별 비공통 오차를 추정하는 단계(S42) 및 위성별 보정치 오차 추정치를 보정치 오차 추정모듈(53)을 통해 구하는 단계는 순차적으로 이루어질 수 있고, 동시에 이루어 질 수도 있다.A step S40 of estimating a pseudorange error by sex using the pseudo range error estimation module 51, a step S42 of estimating a non-common error by sex using the non-common error estimation module 52, The step of obtaining the error estimate value through the correction value error estimation module 53 may be performed sequentially or simultaneously.

상기 의사거리 오차 추정 모듈(51)을 통해 위성별 의사거리 오차를 추정하는 단계(S41)는, 획득 모듈(30)에서 보정 측위모듈(40) 및 오차 추정부(50)에 전달하는 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정한다.The step S41 of estimating the pseudo distance error by sex through the pseudo distance error estimation module 51 is performed by the acquisition module 30 in accordance with the respective distances to be transmitted to the corrected positioning module 40 and the error estimator 50 And determines whether to use the standard deviation of the pseudo range.

이때, 상기 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정할 때 사용할 것으로 결정되면, 관계식 1을 통해 각 위성별 의사거리의 표준편차(

Figure 112016118091315-pat00048
)를 취하여 위성별 의사거리 오차 추정치(
Figure 112016118091315-pat00049
)를 구한다.At this time, if it is decided to use the standard deviation of the pseudoranges for each position, it is determined through use of the relational expression 1 that the standard deviation of the pseudoranges
Figure 112016118091315-pat00048
) And the estimated pseudorange error
Figure 112016118091315-pat00049
).

관계식 1Relationship 1

Figure 112016118091315-pat00050
Figure 112016118091315-pat00050

상기 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정할 때 사용하지 않을 것으로 결정되면, 관계식 2를 통해 각 위성별 의사거리(

Figure 112016118091315-pat00051
)에서 각위성의 시계 오차(
Figure 112016118091315-pat00052
)와 수신기 시계 오차(
Figure 112016118091315-pat00053
)를 보상한 각 위성별 조정된 의사거리 오차(
Figure 112016118091315-pat00054
)를 구한다.If it is determined that the each of the above not used to determine whether to use the standard deviation of the sex pseudoranges, each of the above gender pseudorange through a second relational expression (
Figure 112016118091315-pat00051
), The clock error of each satellite
Figure 112016118091315-pat00052
) And receiver clock error (
Figure 112016118091315-pat00053
) Adjusted pseudorange error for each position compensated (
Figure 112016118091315-pat00054
).

관계식 2Relation 2

Figure 112016118091315-pat00055
Figure 112016118091315-pat00055

이렇게 상기 의사거리 오차(

Figure 112016118091315-pat00056
)가 구해지면, 관계식 3을 통해 각 위성별 조정된 의사거리 오차(
Figure 112016118091315-pat00057
)의 표준편차(
Figure 112016118091315-pat00058
)를 취하여 위성별 의사거리 오차 추정치(
Figure 112016118091315-pat00059
)를 구하게 된다.Thus, the pseudorange error (
Figure 112016118091315-pat00056
) Is obtained when, with each of the above sex adjusted by a relational expression 3 pseudorange error (
Figure 112016118091315-pat00057
) Of the standard deviation
Figure 112016118091315-pat00058
) And the estimated pseudorange error
Figure 112016118091315-pat00059
).

관계식 3Relation 3

Figure 112016118091315-pat00060
Figure 112016118091315-pat00060

상기 비공통 오차 추정 모듈(52)을 통해 위성별 비공통 오차를 추정하는 단계(S42)는, DGNSS 정보의 기준국 위치와 보정측위 모듈(40)에서 측정된 사용자의 DGNSS 위치를 이용하여 기준국과 사용자 사이의 거리(

Figure 112016118091315-pat00061
)를 구하고, 이를 관계식 4을 통해 비공통 오차(
Figure 112016118091315-pat00062
)를 추정하게 된다.The step of estimating the non-common error for each gender using the non-common error estimation module 52 is performed by using the reference station position of the DGNSS information and the DGNSS position of the user measured by the correction positioning module 40, And the distance between the user (
Figure 112016118091315-pat00061
) To obtain, this non-common error by the relation 4 (
Figure 112016118091315-pat00062
).

관계식 4Relation 4

Figure 112016118091315-pat00063
Figure 112016118091315-pat00063

이때, 기준국과 사용자 수신기 사이의 이격이 증가함에 따라 동일 위성이 송신한 신호가 이온층/대류층을 통과하면서 지연되는 량의 상관도가 하락하며, 이로 인하여 기준국과 사용자 사이의 거리 1

Figure 112016118091315-pat00064
당 각 의사거리 측정치에 약 0.57
Figure 112016118091315-pat00065
만큼의 1 시그마 오차가 발생하여 DGNSS 측위 결과에 영향을 미치게 된다.At this time, as the distance between the reference station and the user receiver increases, the degree of correlation of the amount of delay that the signal transmitted by the same satellite passes through the ionosphere / convection layer decreases. As a result,
Figure 112016118091315-pat00064
Approximately 0.57 for each pseudorange measurement
Figure 112016118091315-pat00065
1 sigma error occurs, which affects the DGNSS positioning result.

상기 보정치 오차 추정모듈(53)을 통해 위성별 보정치 오차를 추정하는 단계(S43)는, 보정정보 획득 모듈(35)에서 보정 측위모듈(40) 및 오차 추정부(50)에 전달하는 각 위성별 UDRE를 사용할 것인지의 여부를 결정하며, 각 위성별 UDRE를 사용할 것인지의 여부를 결정할 때 사용할 것으로 결정되면, 관계식 5을 통해 각 위성별 UDRE(

Figure 112017036602278-pat00066
)를 취하여 위성별 보정치 오차 추정치(
Figure 112017036602278-pat00067
)를 구한다.The step S43 of estimating the correction value for each sex by the correction value error estimation module 53 is performed by the correction information acquisition module 35 in accordance with each of the above- when determining whether to use the UDRE and determined to use when deciding whether or not to use each of the above gender UDRE, each of the above gender relations through 5 UDRE (
Figure 112017036602278-pat00066
) To obtain the correction value error estimate (
Figure 112017036602278-pat00067
).

관계식 5Relation 5

Figure 112016118091315-pat00068
Figure 112016118091315-pat00068

상기 각 위성별 UDRE를 사용할 것인지의 여부를 결정할 때 사용하지 않을 것으로 결정되면, 관계식 6을 통해 각 위성별 보정치(

Figure 112016118091315-pat00069
)의 표준편차를 취하여 위성별 보정치 오차 추정치(
Figure 112016118091315-pat00070
)로 구하게 된다.If it is determined that the each of the above not used to determine whether or not to use the gender UDRE, each of the above gender correction value through equations 6 (
Figure 112016118091315-pat00069
), And calculate the error correction value
Figure 112016118091315-pat00070
).

관계식 6Relation 6

Figure 112016118091315-pat00071
Figure 112016118091315-pat00071

이때, 상기 위성별 UDRE는 (

Figure 112017036602278-pat00236
)는 의사거리 보정치에 포함되었을 것으로 추정되는 1시그마(
Figure 112017036602278-pat00237
) 오차로서, DGNSS를 서비스하기 위한 RTCM의 표준문서(SC104-STD)에 따라 기준국의 보정정보에 포함되어 위성별 보정치와 함께 DGNSS 기준국에서 방송된다.At this time, the UDRE by state is (
Figure 112017036602278-pat00236
) Is the one-sigma (< RTI ID = 0.0 >
Figure 112017036602278-pat00237
) Is included in the correction information of the reference station according to the RTCM standard document (SC104-STD) for service of the DGNSS, and is broadcast on the DGNSS reference station together with the correction value of the station.

도 6은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 기준국과 사용자 수신기 사이의 거리에 따라 증가하는 오차를 나타낸 도면이며, 도 7은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 UDRE의 스케일 팩터(Scale Factor)를 나타낸 도면이고, 도 8은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 1시그마 오차에 대한 UDRE(User Differential Range Error) 수치를 나타낸 도면이다.FIG. 6 is a diagram illustrating an error that increases according to the distance between the reference station and the user receiver in the integrity checking method using the integrity checking system of the DGNSS measurement position according to the present invention. FIG. FIG. 8 is a diagram illustrating a scale factor of the UDRE of the integrity checking method using the inspection system. FIG. 8 is a diagram illustrating a scale factor of the UDRE of the integrity checking method using the inspection system. Range Error).

도 6 내지 도 8을 참조하면, 상기 보정치 오차 추정치를 보정치 오차 추정모듈(53)을 통해 추정하는 단계(S40)에서 구해지는 상기 위성별 UDRE(

Figure 112016118091315-pat00074
)는 의사거리 보정치에 포함되었을 것으로 추정되는 1시그마(
Figure 112016118091315-pat00075
) 오차로서, 기준국의 보정정보에 포함되어 보정치와 함께 UDRE(User Differntial Range Error)라는 요소로서 기준국에서 방송된다.Referring to FIGS. 6 to 8, the correction value error estimating module 53 estimates the correction value error estimation value in step S40,
Figure 112016118091315-pat00074
) Is the one-sigma (< RTI ID = 0.0 >
Figure 112016118091315-pat00075
), Which is included in the correction information of the reference station and is broadcast on the reference station as an element called UDRE (User Differential Range Error) together with the correction value.

이때, 상기 위성별 UDRE는 도 7 및 도 8의 Scale Factor에 One-Sigma Differential Error 한계 값을 곱하여 사용한다.At this time, the above-mentioned UDRE by UDRE is used by multiplying the scale factor of FIG. 7 and FIG. 8 by the One-Sigma Differential Error limit value.

상기 오차 추정부(50)를 통해 각각 추정하는 단계(S40)에서 추정된 위성별 오차 추정치 및 보정 측위모듈(40)에서 측정된 사용자의 DGNSS 위치를 통해 생성모듈(60)에서 보호수준을 계산하는 단계(S50)를 통해 사용자가 원하는 수준의 보호수준을 계산하게 된다.The generation module 60 calculates the protection level through the error estimation values estimated in step S40 and the DGNSS position measured by the correction positioning module 40 through the error estimation unit 50, In step S50, the user calculates a desired level of protection.

도 9는 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법의 보호수준을 결정하는 과정에 대한 순서도이다.FIG. 9 is a flowchart illustrating a procedure for determining a protection level of an integrity checking method using an integrity checking system for a DGNSS measurement location according to the present invention.

이때, 상기 보호수준을 계산하는 단계(S50)는, 도 9와 같이 오차 모델(

Figure 112016118091315-pat00076
)을 구성하는 단계(S51), 단위벡터 행렬(G)을 구성하는 단계(S52), 가중치 행렬(W)를 구하는 단계(S53), 투영행렬(S)을 구하는 단계(S54), 동쪽 분산오차, 북쪽 분산오차, 동북쪽 공분산 오차 및 수직방향 분산 오차를 구하는 단계(S55) 및 수평보호수준(HPL, Horizontal Protection Level)과 수직보호수준(VPL, Vertical Protection Level)을 식 10식 11을 통해 구하는 단계(S56)로 이루어진다.At this time, the step of calculating the protection level (S50)
Figure 112016118091315-pat00076
A step S52 of constructing a unit vector matrix G, a step S53 of obtaining a weighting matrix W, a step S54 of obtaining a projection matrix S, (S55), a horizontal protection level (HPL), and a vertical protection level (VPL) are obtained through equations (10 ) and (11 ) Step S56 is performed.

상기 오차 모델(

Figure 112016118091315-pat00077
)을 구성하는 단계(S51)는, 오차 추정부(50)를 통해 각각 추정하는 단계(S40)에서 추정되어 제공되는 DGNSS 정보를 통해 위성별 의사거리 오차 추정치(
Figure 112016118091315-pat00078
), 위성별 비공통 오차 추정치(
Figure 112016118091315-pat00079
) 및 위성별 보정치 오차 추정치(
Figure 112016118091315-pat00080
)를 식 1에 대입하여 오차 모델(
Figure 112016118091315-pat00081
)을 구하게 된다.The error model (
Figure 112016118091315-pat00077
(Step S51) is performed by using the DGNSS information estimated in step S40, which is estimated through the error estimation unit 50,
Figure 112016118091315-pat00078
), Non-common error estimates by sex
Figure 112016118091315-pat00079
) And the gender-specific correction value error estimate
Figure 112016118091315-pat00080
) To the equation (1 ) to obtain the error model
Figure 112016118091315-pat00081
).

식 1Equation 1

Figure 112016118091315-pat00082
Figure 112016118091315-pat00082

상기 오차 모델(

Figure 112017036602278-pat00083
)이 구해지면, 단위벡터 행렬(G)을 구성하는 단계(S52)를 제공하는 단계(S20)의 획득모듈(30)에서 제공한 각 위성의 위치 및 사용자의 DGNSS 위치를 측정하는 단계(S30)의 보정측위 모듈(40)에서 측정된 사용자의 DGNSS 위치로 ENU 좌표계에서 사용자 위치로부터 위성의 방향에 대한 각 앙각(
Figure 112017036602278-pat00084
)과 방위각(
Figure 112017036602278-pat00085
)을 구하고, 식 2를 통해 위성 배치를 나타내는 단위벡터 행렬(
Figure 112017036602278-pat00086
)을 구성하게 된다.The error model (
Figure 112017036602278-pat00083
A step S30 of measuring the position of each satellite and the DGNSS position of the user provided by the acquisition module 30 of step S20 of providing a step S52 of constructing a unit vector matrix G, To the DGNSS position of the user measured by the calibration positioning module 40 in the ENU coordinate system,
Figure 112017036602278-pat00084
) And azimuth angle
Figure 112017036602278-pat00085
) To obtain a unit vector matrix representing the satellite deployment through a formula (2) (
Figure 112017036602278-pat00086
).

식 2Equation 2

Figure 112016118091315-pat00087
Figure 112016118091315-pat00087

상기 단위벡터 행렬(G)을 구성하면, 가중치 행렬(W)를 구하는 단계(S53)에서 오차 모델을 구성하는 단계(S51)를 통해 구성된 위성별 오차모델(

Figure 112016118091315-pat00088
) 값을 통해 오차모델 값의 역수값(
Figure 112016118091315-pat00089
)을 식 3으로 구하고, 상기 역수값(
Figure 112016118091315-pat00090
)을 식 4에 대입하여 식 4를 통해 가중치 행렬(
Figure 112016118091315-pat00091
)를 구하게된다.The unit vector matrix G is constructed so that the weight matrix W is obtained through the step S51 of constructing the error model in the step S53 of obtaining the weight matrix W
Figure 112016118091315-pat00088
) Value, the reciprocal value of the error model value (
Figure 112016118091315-pat00089
) Is obtained by the equation (3 ), and the reciprocal value
Figure 112016118091315-pat00090
), The weight matrix through formula (4) are substituted in equation (4) (
Figure 112016118091315-pat00091
).

식 3Equation 3

Figure 112016118091315-pat00092
Figure 112016118091315-pat00092

식 4Equation 4

Figure 112016118091315-pat00093
Figure 112016118091315-pat00093

상기 가중치 행렬(W)이 구해지면, 투영행렬(S)을 구하는 단계(S54)에서 단위벡터 행렬()을 구성하는 단계(S52)를 통해 구해진 단위벡터 행렬(

Figure 112016118091315-pat00095
)과 가중치 행렬(
Figure 112016118091315-pat00096
)를 구하는 단계(S53)를 통해 구해진 가중치 행렬(
Figure 112016118091315-pat00097
)을 가중최소제곱법인 식 5에 대입하여 거리영역에서 위치영역으로 투영한 투영행렬(
Figure 112016118091315-pat00098
)을 구하게 된다.When the weight matrix W is obtained, the unit matrix W (k) is calculated in step S54 of obtaining the projection matrix S (S52) of the unit vector matrix < RTI ID = 0.0 >
Figure 112016118091315-pat00095
) And a weighting matrix (
Figure 112016118091315-pat00096
(Step S53) of obtaining a weight matrix
Figure 112016118091315-pat00097
) Is substituted into the weighted least squares equation (5 ), and a projection matrix (
Figure 112016118091315-pat00098
).

식 5Equation 5

Figure 112016118091315-pat00099
Figure 112016118091315-pat00099

상기 투영행렬(S)이 구해지면, 분산 오차를 구하는 단계(S55)에서 투영행렬(

Figure 112016118091315-pat00100
)을 구하는 단계(S54)를 통해 구해진 투영행렬(
Figure 112016118091315-pat00101
)의 일부 요소를 취하여 상기 오차 모델을 구성하는 단계(S51)에서 구한 위성별 오차 모델을 곱한 후 더하는 방법으로 식 6, 식 7, 식 8식 9를 통해 동쪽 방향으로의 분산 오차(
Figure 112016118091315-pat00102
), 북쪽 방향으로의 분산 오차(
Figure 112016118091315-pat00103
), 동쪽과 북쪽 방향으로의 공분산 오차(
Figure 112016118091315-pat00104
) 및 수직 방향으로의 분산 오차(
Figure 112016118091315-pat00105
)를 구하게 된다.When the projection matrix S is obtained, in the step S55 of obtaining the dispersion error, the projection matrix (
Figure 112016118091315-pat00100
(S54) of obtaining the projection matrix (
Figure 112016118091315-pat00101
) Is calculated by multiplying the error model obtained in the step S51 of constructing the error model by a factor of 6, 7, 8 and 9 to calculate the variance error in the east direction
Figure 112016118091315-pat00102
), The dispersion error in the north direction (
Figure 112016118091315-pat00103
), The covariance error in the east and north directions
Figure 112016118091315-pat00104
) And the dispersion error in the vertical direction (
Figure 112016118091315-pat00105
).

식 6Equation 6

Figure 112016118091315-pat00106
Figure 112016118091315-pat00106

식 7Equation 7

Figure 112016118091315-pat00107
Figure 112016118091315-pat00107

식 8Equation 8

Figure 112016118091315-pat00108
Figure 112016118091315-pat00108

식 9Equation 9

Figure 112016118091315-pat00109
Figure 112016118091315-pat00109

상기 분산오차들이 구해지면, 수평보호수준(HPL, Horizontal Protection Level)과 수직보호수준(VPL, Vertical Protection Level)을 식 10식 11을 통해 구하는 단계(S56)에서 분산 오차를 구하는 단계(S55)를 통해 구해진 분산 오차 요소와 수평보호수준(

Figure 112016118091315-pat00110
, Horizontal Protection Level)에 대하여 사용자가 요구하는 무결성 리스크를 만족시키는 표준점수(
Figure 112016118091315-pat00111
)와 수직보호수준(
Figure 112016118091315-pat00112
, Vertical Protection Level)에 대하여 사용자가 요구하는 무결성 리스크를 만족시키는 표준점수(
Figure 112016118091315-pat00113
)를 이용하여 (HPL, Horizontal Protection Level)과 수직보호수준(VPL, Vertical Protection Level)을 식 10식 11을 통해 구하게 된다.When the dispersion errors are obtained, a dispersion error is obtained (S55) in a step S56 of obtaining a horizontal protection level (HPL) and a vertical protection level (VPL) through Equations 10 and 11 , And the horizontal protection level (
Figure 112016118091315-pat00110
, Horizontal Protection Level) that satisfies the user's required integrity risk.
Figure 112016118091315-pat00111
) And vertical protection level
Figure 112016118091315-pat00112
, Vertical Protection Level), which satisfies the user's required integrity risk.
Figure 112016118091315-pat00113
(HPL) and vertical protection level (VPL) using Equation (10 ) and Equation (11 ).

식 10Equation 10

Figure 112016118091315-pat00114
Figure 112016118091315-pat00114

수평보호수준(HPL)은 타원의 장축을 구하는 식으로 구해짐.The horizontal protection level (HPL) is obtained by obtaining the long axis of the ellipse.

식 11Equation 11

Figure 112016118091315-pat00115
Figure 112016118091315-pat00115

이처럼 상기 보호수준을 계산하는 단계(S50)를 통해 계산하여 생성한 보호수준은 무결성 검사모듈(70)을 통해 결정하는 단계(S60)를 통해 사용자가 원하는 경보한계(Alert Limit)와 비교하고, 경보한계와 비교된 보호수준을 해당 DGNSS 위치 측정치를 항법에 활용할 것인지 여부를 무결성 검사모듈(70)을 통해 결정하게 된다.The protection level calculated through the calculation of the protection level S50 is compared with the desired alarm limit through the integrity checking module 70 in step S60, The integrity check module 70 determines whether the protection level compared to the limit is to be utilized for navigation in the DGNSS location measurement.

상기 무결성 검사모듈(70)을 통해 결정하는 단계(S60)의 경보한계와 위치오차 및 보호수준의 상관 관계를 아래 표 1을 통해 나타낸다.The correlation between the alarm limit, the position error and the protection level of the step S60 of determining through the integrity checking module 70 is shown in Table 1 below.

위치오차(PE) < 보호수준(PL) (99.999%)Position error (PE) <Protection level (PL) (99.999%) 보호수준(PL) < 위치오차(PE) (00.001%)Protection level (PL) <Position error (PE) (00.001%) 관계relation 경보Warning 문제Problem 관계relation 경보Warning 설명Explanation 1One PE < PL < ALPE <PL <AL -- -- 44 PL < PE < ALPL <PE <AL -- -- 22 PE < AL < PLPE <AL <PL 경보Warning 가용성 저하Low availability 55 PL < AL < PEPL <AL <PE -- 무결성상실Loss of Integrity 33 AL < PE < PLAL <PE <PL 경보Warning -- 66 AL < PL < PEAL <PL <PE 경보Warning --

PE-Position Error, AL-Alert Limit, PL-Protection LevelPE-Position Error, AL-Alert Limit, PL-Protection Level

상기 표 1은 0.00001의 무결성 리스크를 만족시키는 표준점수 4.42를

Figure 112016118091315-pat00116
Figure 112016118091315-pat00117
에 적용한 경우 무결성 분석 결과가 발생할 수 있는 경우를 나타낸다.Table 1 above shows a standard score of 4.42 satisfying the integrity risk of 0.00001
Figure 112016118091315-pat00116
and
Figure 112016118091315-pat00117
, The integrity analysis result can be generated.

상기 표 1의 전체 6가지의 경우 중 1, 2, 3번의 경우가 합하여 99.999%의 확률로 발생하며, 4, 5, 6번의 경우가 합하여 0.001% 확률로 발생한다.Of the six cases shown in Table 1, the cases 1, 2, and 3 occur together at a probability of 99.999%, and the cases 4, 5, and 6 occur at a probability of 0.001%.

상기 표 1을 설명하면, 1번의 경우는 위치오차와 보호수준이 모두 경보한계 미만으로 경보할 필요가 없고, 경보가 울리지 않는다.As shown in Table 1, in the case of No. 1, it is unnecessary to warn the position error and the protection level to be below the alarm limit, and the alarm does not sound.

상기 2번의 경우는 위치오차가 경보한계 미만인데 보호수준이 경보한계 이상인 경우로, 경보를 발생할 필요가 없었지만 발생한 경우를 의미한다. 이 경우는 무결성 성능과는 상관이 없으나 불필요한 경우에 경보를 발생시켜 시스템 전체 가용성이 저하 된다.In case 2, the position error is below the alarm limit, but the protection level is above the alarm limit, which means that the alarm was not generated but occurred. In this case, there is no relation to the integrity performance, but when it is unnecessary, an alarm is generated and the overall system availability is degraded.

상기 3번의 경우는 위치오차가 경보한계를 넘어서 경보를 발생시켜야 하며, 이 때 보호수준도 경보한계를 넘어 경보를 잘 발생시킨 경우를 나타낸다.In case of above 3, the position error must exceed the alarm limit, and the alarm level should exceed the alarm limit.

상기 4번의 경우는 보호수준이 위치오차보다 낮게 계산되었지만, 결과적으로 경보를 발생시키지 않아야 할 때 발생시키지 않은 정상상황이다.In case 4, the protection level is calculated to be lower than the position error, but it is a normal situation that does not occur when the alarm should not be generated.

상기 5번의 경우는 보호수준이 위치오차보다 낮게 계산되어, 경보를 발생시켜야 하는데 발생시키지 못하여 무결성이 상실된 경우이다.In case 5, the protection level is calculated to be lower than the position error, and the alarm is generated, but the integrity is lost due to failure.

상기 6번의 경우는 보호수준이 위치오차보다 낮게 계산되었지만, 결과적으로 경보를 발생시켜야 할 때 경보를 발생한 정상상황이다. In case of 6 above, the protection level is calculated to be lower than the position error, but as a result, when the alarm should be generated, the alarm is generated in a normal state.

도 10은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법에 따른 실험 결과를 나타낸 그래프이며, 도 11은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법에 따른 실험결과에서 나타난 활용 횟수에 따른 가용성을 나타낸 도면이다.FIG. 10 is a graph showing an experimental result according to the integrity checking method using the integrity checking system of the DGNSS measurement position according to the present invention. FIG. 11 is a graph showing the results of the integrity checking method using the integrity checking system of the DGNSS measuring position according to the present invention. And the availability according to the number of utilization shown in the result.

상기 본 발명의 구성을 통해 나타나는 실험 결과는 도 10에 도시된 바와 같이 나타난다.The result of the experiment through the configuration of the present invention is shown in FIG.

상기 도 10은 거제도에 정밀하게 측정된 위치에서 GPS 사용자 수신기르 설치하고 영도 기준국의 보정정보를 받아 DGPS 측위를 하고 보호수준을 계산하였다.In FIG. 10, a GPS user receiver is installed at a precisely measured position in a geographical diagram, DGPS positioning is performed by receiving the correction information of the reference station, and the level of protection is calculated.

상기 도 10의 그래프는 계산된 보호수준과 경보한계 그리고 실제 위치의 오차를 나타내며, 파란점이 위치의 오차이고 검정색 선이 제안된 방법으로 계산된 보호수준이다. The graph of FIG. 10 shows the calculated protection level, alarm limit and error of the actual position, the blue dot is the position error and the black line is the protection level calculated by the proposed method.

이러한 도 10의 그래프에서 붉은 색 선은 사용자가 요구하는 경보한계를 나타내고, 대부분의 경우 보호수준이 실제 오차를 커버하는 것을 확인할 수 있다.In the graph of FIG. 10, the red line indicates the alarm limit required by the user, and in most cases, the protection level covers the actual error.

그리고, 붉은색 경보 한계 위에 나타난 두 개의 파란 점은 실제 위치 오차가 25m 이상인 경우를 나타내며, 이 때 계산된 보호수준도 25m를 넘어 잘 경보를 해낸 것을 확인할 수 있다.The two blue dots on the red alarm limit indicate that the actual position error is greater than 25m, and the calculated protection level is also over 25m, indicating a good alarm.

상기 도 10의 그래프를 보면 작은 붉은 원은 2번 발생하였는데, 이 경우는 무결성을 상실하지는 않았지만 보호수준이 위치오차보다 크게 계산된 경우로 0.001%의 확률로 발생할 것으로 기대된다. 실험은 342835 epoch 실행되었으며, 이 때 기대값은 약 3.4회가 되며, 실제 2회 발생하였다.Referring to the graph of FIG. 10, a small red circle occurs twice. In this case, the degree of protection is estimated to be greater than the positional error, which is expected to occur at a probability of 0.001% although the integrity is not lost. Experiments were performed on 342835 epochs, with an expectation value of about 3.4 times, which actually occurred twice.

도 11은 본 발명에 따른 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법에 따른 실험결과에서 나타난 활용 횟수에 따른 가용성을 나타낸 도면이다.FIG. 11 is a diagram illustrating availability according to the number of times of use, which is an experiment result according to an integrity checking method using an integrity checking system for a DGNSS measurement position according to the present invention.

도 11에 도시된 바와 같이 경보를 하여 사용자가 활용할 수 없었던 경우는 342835회 중 208회 발생하였으며, 이 결과는 약 0.06%의 가용성 저하를 보여준다.As shown in FIG. 11, when the user could not utilize the alarm, 208 times out of 342835 times occurred, and the result shows a decrease in availability of about 0.06%.

이와 같이 이루어지는 본 발명은 GNSS 위성신호와 DGNSS 보정신호를 수신하여 항해하는 사용자가 외부의 지원없이 DGNSS 수신기 자체만으로 측정 위치의 무결성을 검사할 수 있는 효과가 있다.According to the present invention as described above, a user who receives a GNSS satellite signal and a DGNSS correction signal and sails can check the integrity of a measurement location with only the DGNSS receiver itself without external support.

또한, 무결성 검사를 수행함에 있어서 사용자가 원하는 수준의 무결성 리스크를 충족시키는 보호수준을 계산할 수 있고, 사용자가 원하는 크기의 경보 한계와 비교할 수 있는 효과가 있다.Also, in performing the integrity check, the user can calculate the protection level satisfying the desired level of integrity risk, and the user can compare with the alarm limit of the desired size.

상기와 같은 DGNSS 측정 위치의 무결성 검사 시스템은 위에서 설명된 실시예들의 구성과 작동 방식에 한정되는 것이 아니다. 상기 실시예들은 각 실시예들의 전부 또는 일부가 선택적으로 조합되어 다양한 변형이 이루어질 수 있도록 구성될 수도 있다.The integrity checking system of the DGNSS measurement position is not limited to the configuration and operation of the embodiments described above. The embodiments may be configured so that all or some of the embodiments may be selectively combined so that various modifications may be made.

10 : DGNSS 수신모듈 20 : 보정정보 수신모듈
30 : 획득모듈 35 : 보정정보 획득모듈
40 : 보정 측위모듈
50 : 오차 추정부 51 : 의사거리 오차 추정보듈
52 : 비공통 오차 추정모듈 53 : 보정치 오차 추정모듈
60 : 생성모듈 70 : 무결성 검사모듈
10: DGNSS receiving module 20: correction information receiving module
30: Acquisition module 35: Calibration information acquisition module
40: Calibration positioning module
50: error estimation unit 51: pseudo range error estimation module
52: non-common error estimation module 53: correction value error estimation module
60: Generation module 70: Integrity check module

Claims (8)

GNSS 위성군에서 송신되는 GNSS 신호를 수신하고 해석하여 GNSS 원시데이터를 제공하는 GNSS(위성항법 시스템: Global Navigation Satellite System) 수신모듈(10);
상기 GNSS 수신모듈(10)에서 제공한 GNSS 원시데이터를 받아 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보를 획득하여 오차 추정부(50) 및 보정측위모듈(40)에 제공하는 획득모듈(30);
DGNSS 기준국에서 송신되는 DGNSS 신호를 수신하고 해석하여 DGNSS 원시데이터를 제공하는 보정정보 수신모듈(20);
상기 보정정보 수신모듈(20)에서 제공한 DGNSS 원시데이터를 받아 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 획득하여 오차 추정부(50) 및 보정측위모듈(40)에 제공하는 보정정보 획득모듈(35);
상기 획득모듈(30)에서 제공한 GNSS 정보와 보정정보 획득모듈(35)에서 제공한 DGNSS 정보를 이용하여 차분보정(Differential correction) 방법으로 사용자의 DGNSS 위치를 측정하고, 측정된 사용자의 DGNSS 위치와 DGNSS 위치 측정에 사용된 각 위성의 위치를 오차 추정부(50) 및 생성모듈(60)에 제공하는 보정 측위모듈(40);
상기 획득모듈(30)에서 제공한 GNSS 정보를 이용하여 추정한 위성별 의사거리 오차 추정치를 제공하는 의사거리 오차 추정모듈(51)과, 보정정보 획득모듈(35)에서 제공한 DGNSS 정보와 보정 측위모듈(40)에서 제공한 사용자의 DGNSS 위치를 이용하여 추정한 위성별 비공통 오차 추정치를 제공하는 비공통 오차 추정모듈(52) 및 보정정보 획득모듈(35)에서 제공한 DGNSS 정보를 이용하여 추정한 위성별 보정치 오차 추정치를 제공하는 보정치 오차 추정모듈(53)로 이루어지는 오차 추정부(50);
상기 보정 측위모듈(40) 및 오차 추정부(50)를 통해 제공되는 사용자의 DGNSS 위치, DGNSS 위치 측정에 사용된 각 위성의 위치, 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치, 위성별 보정치 오차 추정치를 이용하여 사용자가 원하는 수준의 보호수준(Protection Level)을 계산하여 제공하는 생성모듈(60) 및
상기 생성모듈(60)을 통해 제공되는 보호수준을 받아 사용자가 원하는 경보한계(Alert Limit)와 비교하여 사용자의 해당 DGNSS 위치 측정치를 항법에 활용할 것인지 여부를 결정하는 무결성 검사모듈(70)을 포함하는 것을 특징으로 하는 DGNSS 측정 위치의 무결성 검사 시스템.
A GNSS (Global Navigation Satellite System) receiving module 10 for receiving and interpreting GNSS signals transmitted from the GNSS satellite group and providing GNSS raw data;
And receives GNSS raw data provided by the GNSS receiving module 10, and receives the GNSS raw data provided by the GNSS receiving module 10, such as a measurement time, a receiver clock error, a pseudorange by each level, a standard deviation of pseudoranges by each level, An acquisition module (30) for acquiring and providing to the error estimation unit (50) and the correction positioning module (40);
A correction information receiving module (20) for receiving and interpreting a DGNSS signal transmitted from a DGNSS reference station and providing DGNSS raw data;
And receives DGNSS raw data provided by the correction information receiving module 20 to obtain correction value generation time (MZC, Modified-Z Count), statistical correction value, UDRE (User Differential Range Error) and DGNSS information A correction information acquisition module (35) for providing the correction information to the error estimation unit (50) and the correction positioning module (40);
The DGNSS position of the user is measured by a differential correction method using the GNSS information provided by the acquisition module 30 and the DGNSS information provided by the correction information acquisition module 35, A correction positioning module 40 for providing the position of each satellite used for the DGNSS position measurement to the error estimation unit 50 and the generation module 60;
A pseudo distance error estimation module 51 for providing a pseudorange distance error estimation value by using the GNSS information provided by the acquisition module 30, a DGNSS information providing module Common error estimation module 52 for providing a gender-specific non-common error estimation value estimated using the DGNSS position of the user provided by the module 40 and a DGNSS information provided by the correction information acquisition module 35, An error estimator (50) comprising a correction value error estimation module (53) for providing a correction value error estimate for the upper gender;
A DGNSS position of the user provided through the correction positioning module 40 and the error estimation unit 50, a position of each satellite used for the DGNSS position measurement, a pseudorange error estimate by sex, a non-common error estimate by sex, A generation module 60 for calculating and providing a protection level of a desired level by using a correction value error estimation value,
And an integrity checking module 70 for determining whether to use the user's corresponding DGNSS position measurement value in navigation by comparing the user's desired alert limit with a desired protection level provided by the creation module 60 The integrity check system of the DGNSS measurement location.
제1항에 있어서,
상기 보정정보 수신모듈(20)은 DGNSS 기준국에서 송신되는 DGNSS 신호를 중파, NTRIP, Eurofix, AIS(Auto Identification System) 또는 DMB(Digital Multimedia Broadcasting) 중 어느 하나 또는 다수를 통하여 수신받는 것을 특징으로 하는 DGNSS 측정 위치의 무결성 검사 시스템.
The method according to claim 1,
The correction information receiving module 20 receives a DGNSS signal transmitted from a DGNSS reference station through one or more of a medium wave, an NTRIP, a Eurofix, an AIS (Auto Identification System), or a DMB (Digital Multimedia Broadcasting) DGNSS Integrity check system of measurement location.
DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법은,
GNSS 수신모듈(10)에서 GNSS 위성군의 GNSS 신호를 수신하여 해석한 GNSS 원시데이터 및 보정정보 수신모듈(20)에서 DGNSS 기준국의 DGNSS 신호를 수신하여 해석한 DGNSS 원시데이터를 획득모듈(30)과 보정정보 획득모듈(35)에서 제공받는 단계(S10);
상기 제공받는 단계(S10)에서 GNSS 수신모듈(10)을 통해 GNSS 원시데이터를 제공받아 획득모듈(30)에서 생성한 측정 시각, 수신기 시계 오차, 각 위성별 의사거리, 각 위성별 의사거리의 표준편차, 각 위성의 위치, 각 위성의 시계 오차인 GNSS 정보 및 상기 제공받는 단계(S10)에서 보정정보 수신모듈(20)을 통해 DGNSS 원시데이터를 제공받아 보정정보 획득모듈(35)에서 생성한 보정치 생성 시각(MZC, Modified-Z Count), 위성별 보정치, 위성별 UDRE(User Differential Range Error) 및 기준국 위치인 DGNSS 정보를 오차 추정부(50) 및 보정측위모듈(40)에 제공하는 단계(S20);
상기 제공하는 단계(S20)를 통해 생성된 GNSS 정보와 DGNSS 정보를 이용하여 보정 측위모듈(40)에서 차분보정(Differential correction) 방법으로 사용자의 DGNSS 위치를 측정하는 단계(S30);
상기 사용자의 DGNSS 위치를 측정하는 단계(S30)를 통해 측정된 사용자의 DGNSS 위치, 상기 제공하는 단계(S20)를 통해 획득모듈(30)에서 생성한 GNSS 정보 및 상기 보정정보 획득모듈(35)에서 생성한 DGNSS 정보를 통해 위성별 의사거리 오차 추정치, 위성별 비공통 오차 추정치 및 위성별 보정치 오차 추정치를 오차 추정부(50)를 통해 각각 추정하는 단계(S40);
상기 오차 추정부(50)를 통해 각각 추정하는 단계(S40)에서 추정된 위성별 오차 추정치 및 보정 측위모듈(40)에서 측정된 사용자의 DGNSS 위치를 통해 생성모듈(60)에서 사용자가 원하는 수준의 보호수준을 계산하는 단계(S50) 및
상기 보호수준을 계산하는 단계(S50)를 통해 계산하여 생성한 보호수준을 사용자가 원하는 경보한계(Alert Limit)와 비교하여 해당 DGNSS 위치 측정치를 항법에 활용할 것인지 여부를 무결성 검사모듈(70)을 통해 결정하는 단계(S60)를 포함하는 것을 특징으로 하는 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법.
The integrity check method using the integrity check system of the DGNSS measurement location,
The GNSS source data obtained by receiving the GNSS signal of the GNSS satellite group from the GNSS receiving module 10 and the DGNSS source signal of the DGNSS reference station received from the correction information receiving module 20 and interpreting the DGNSS source data, And a step (S10) provided by the correction information acquisition module 35;
In the receiving step S10, GNSS source data is received through the GNSS receiving module 10, and the measurement time generated by the acquisition module 30, the receiver clock error, the pseudo distance by each level, The GNSS information which is the clock error of each satellite and the DGNSS raw data through the correction information receiving module 20 in the receiving step S10 and received by the correction information obtaining module 35, Providing the error estimation unit 50 and the correction positioning module 40 with the generation time (MZC, Modified-Z Count), the statistical correction value, the UDRE (User Differential Range Error) S20);
A step (S30) of measuring a DGNSS position of a user by a differential correction method in the correction positioning module 40 using the GNSS information and the DGNSS information generated through the providing step S20;
The DGNSS position of the user measured through the DGNSS position measurement of the user (S30), the GNSS information generated by the acquisition module 30 through the providing step S20, and the GNSS information generated by the correction information acquisition module 35 (S40) estimating the pseudo distance error estimate, the gender-by-common non-common error estimation value, and the gender-specific correction value error estimation value through the error estimation unit 50 through the generated DGNSS information, respectively;
The generation module 60 estimates the user's desired level of error by using the estimated error by the estimated error in step S40 and the DGNSS position of the user measured by the corrected positioning module 40 through the error estimator 50, Calculating a protection level (S50) and
The protection level calculated and generated through the calculation of the protection level (S50) is compared with a desired alert limit (Alert Limit), and whether the DGNSS position measurement is used for navigation is checked through the integrity checking module 70 And a step (S60) of determining the integrity of the DGNSS measurement location.
제3항에 있어서,
상기 오차 추정부(50)를 통해 각각 추정하는 단계(S40)는,
상기 획득모듈(30)에서 제공한 GNSS 정보를 이용하여 의사거리 오차 추정모듈(51)을 통해 위성별 의사거리 오차를 추정하는 단계(S41);
상기 보정정보 획득모듈(35)에서 제공한 DGNSS 정보 및 보정 측위모듈(40)에서 제공한 사용자의 DGNSS 위치를 이용하여 비공통 오차 추정모듈(52)을 통해 기준국과 사용자 간의 거리를 구하고 이에 따른 위성별 비공통 오차를 추정하는 단계(S42) 및
상기 보정정보 획득모듈(35)에서 제공한 DGNSS 정보를 이용하여 보정치 오차 추정모듈(53)을 통해 위성별 보정치 오차를 추정하는 단계(S43)로 이루어지는 것을 특징으로 하는 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법.
The method of claim 3,
The estimating step S40 through the error estimating unit 50, respectively,
A step (S41) of estimating a pseudorange error by sex using the pseudo range error estimation module 51 using the GNSS information provided by the acquisition module 30;
The distance between the reference station and the user is obtained through the non-common error estimation module 52 using the DGNSS information provided by the correction information acquisition module 35 and the DGNSS position of the user provided by the correction positioning module 40, A step (S42) of estimating non-common error by sex and
And a step (S43) of estimating a threshold correction value error through a correction value error estimation module (53) using the DGNSS information provided by the correction information acquisition module (35). Integrity check method used.
제4항에 있어서,
상기 의사거리 오차 추정모듈(51)을 통해 위성별 의사거리 오차를 추정하는 단계(S41)는,
상기 획득모듈(30)에서 오차 추정부(50) 및 보정측위모듈(40)에 전달하는 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정하며,
상기 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정할 때 사용할 것으로 결정되면, 관계식 1을 통해 각 위성별 의사거리의 표준편차(
Figure 112017036602278-pat00118
)를 취하여 위성별 의사거리 오차 추정치(
Figure 112017036602278-pat00119
)를 구하고,
상기 각 위성별 의사거리의 표준편차를 사용할 것인지의 여부를 결정할 때 사용하지 않을 것으로 결정되면, 관계식 2를 통해 각 위성별 의사거리(
Figure 112017036602278-pat00120
)에서 각위성의 시계 오차(
Figure 112017036602278-pat00121
)와 수신기 시계 오차(
Figure 112017036602278-pat00122
)를 보상한 각 위성별 조정된 의사거리 오차(
Figure 112017036602278-pat00123
)를 구하고 관계식 3을 통해 각 위성별 조정된 의사거리 오차(
Figure 112017036602278-pat00124
)의 표준편차(
Figure 112017036602278-pat00125
)를 취하여 위성별 의사거리 오차 추정치(
Figure 112017036602278-pat00126
)를 구하는 것을 특징으로 하는 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법.
관계식 1
Figure 112017036602278-pat00127

관계식 2
Figure 112017036602278-pat00128

관계식 3
Figure 112017036602278-pat00129

5. The method of claim 4,
The step (S41) of estimating the pseudorange distance error by the pseudo range error estimation module (51)
The acquisition module 30 determines whether or not to use the standard deviation of the pseudo distance for each position transmitted to the error estimation unit 50 and the correction positioning module 40,
Wherein when each of the above when deciding whether to use the standard deviation of the sex pseudoranges determined to be used, the standard deviation of each of the above gender pseudorange through a relational expression 1 (
Figure 112017036602278-pat00118
) And the estimated pseudorange error
Figure 112017036602278-pat00119
),
If it is determined that the each of the above not used to determine whether to use the standard deviation of the sex pseudoranges, each of the above gender pseudorange through a second relational expression (
Figure 112017036602278-pat00120
), The clock error of each satellite
Figure 112017036602278-pat00121
) And receiver clock error (
Figure 112017036602278-pat00122
) Adjusted pseudorange error for each position compensated (
Figure 112017036602278-pat00123
) To obtain the above, each sex adjusted by a relational expression 3 pseudorange error (
Figure 112017036602278-pat00124
) Of the standard deviation
Figure 112017036602278-pat00125
) And the estimated pseudorange error
Figure 112017036602278-pat00126
The integrity check method using the integrity check system of the DGNSS measurement position is characterized by determining the integrity check method.
Relationship 1
Figure 112017036602278-pat00127

Relation 2
Figure 112017036602278-pat00128

Relation 3
Figure 112017036602278-pat00129

제4항에 있어서,
상기 비공통 오차 추정모듈(52)을 통해 위성별 비공통 오차를 추정하는 단계(S42)는,
상기 DGNSS 정보의 기준국 위치와 상기 보정 측위모듈(40)에서 측정된 사용자의 DGNSS 위치를 이용하여 기준국과 사용자 사이의 거리(
Figure 112017036602278-pat00130
)를 구하고, 이를 관계식 4을 통해 비공통 오차(
Figure 112017036602278-pat00131
)를 추정하는 것을 특징으로 하는 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법.
관계식 4
Figure 112017036602278-pat00132

이때, 기준국과 사용자 수신기 사이의 이격이 증가함에 따라 동일 위성이 송신한 신호가 이온층/대류층을 통과하면서 지연되는 량의 상관도가 하락하며, 이로 인하여 기준국과 사용자 사이의 거리 1
Figure 112017036602278-pat00133
당 각 의사거리 측정치에 0.57
Figure 112017036602278-pat00134
만큼의 1 시그마 오차가 발생하여 DGNSS 측위 결과에 영향을 미침.
5. The method of claim 4,
The step (S42) of estimating the non-common error for each sex by the non-common error estimation module (52)
The distance between the reference station and the user using the reference station position of the DGNSS information and the DGNSS position of the user measured by the correction positioning module 40
Figure 112017036602278-pat00130
), And this is calculated as a non-common error (
Figure 112017036602278-pat00131
) Of the DGNSS measurement position is estimated.
Relation 4
Figure 112017036602278-pat00132

At this time, as the distance between the reference station and the user receiver increases, the degree of correlation of the amount of delay that the signal transmitted by the same satellite passes through the ionosphere / convection layer decreases. As a result,
Figure 112017036602278-pat00133
0.57 for each pseudo distance measure
Figure 112017036602278-pat00134
1 sigma error occurs and affects DGNSS positioning result.
제4항에 있어서,
상기 보정치 오차 추정모듈(53)을 통해 위성별 보정치 오차를 추정하는 단계(S43)는,
상기 보정정보 획득모듈(35)에서 오차 추정부(50) 및 보정측위모듈(40)에 전달하는 각 위성별 UDRE를 사용할 것인지의 여부를 결정하며,
상기 각 위성별 UDRE를 사용할 것인지의 여부를 결정할 때 사용할 것으로 결정되면, 관계식 5을 통해 각 위성별 UDRE(
Figure 112017036602278-pat00135
)를 취하여 위성별 보정치 오차 추정치(
Figure 112017036602278-pat00136
)를 구하고,
상기 각 위성별 UDRE를 사용할 것인지의 여부를 결정할 때 사용하지 않을 것으로 결정되면, 관계식 6을 통해 각 위성별 보정치(
Figure 112017036602278-pat00137
)의 표준편차를 취하여 위성별 보정치 오차 추정치(
Figure 112017036602278-pat00138
)로 구하는 것을 특징으로 하는 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법.
이때, 상기 위성별 UDRE는 (
Figure 112017036602278-pat00139
)는 의사거리 보정치에 포함되었을 것으로 추정되는 1시그마(
Figure 112017036602278-pat00140
) 오차로서, DGNSS를 서비스하기 위한 RTCM의 표준문서(SC104-STD)에 따라 기준국의 보정정보에 포함되어 위성별 보정치와 함께 DGNSS 기준국에서 방송된다.
관계식 5
Figure 112017036602278-pat00141

관계식 6
Figure 112017036602278-pat00142

5. The method of claim 4,
The step (S43) of estimating the correction value for each sex by the correction value error estimation module (53)
The correction information acquisition module 35 determines whether or not to use each UDRE for each position to be transmitted to the error estimation unit 50 and the correction positioning module 40,
If it is determined that the use in determining whether or not each of the above to use the gender UDRE, each of the above equations 5 through gender UDRE (
Figure 112017036602278-pat00135
) To obtain the correction value error estimate (
Figure 112017036602278-pat00136
),
If it is determined that the each of the above not used to determine whether or not to use the gender UDRE, each of the above gender correction value through equations 6 (
Figure 112017036602278-pat00137
), And calculate the error correction value
Figure 112017036602278-pat00138
The integrity check method using the integrity check system of the DGNSS measurement position is characterized by that the integrity check method is used.
At this time, the UDRE by state is (
Figure 112017036602278-pat00139
) Is the one-sigma (&lt; RTI ID = 0.0 &gt;
Figure 112017036602278-pat00140
) Is included in the correction information of the reference station according to the RTCM standard document (SC104-STD) for service of the DGNSS, and is broadcast on the DGNSS reference station together with the correction value of the station.
Relation 5
Figure 112017036602278-pat00141

Relation 6
Figure 112017036602278-pat00142

제3항에 있어서,
상기 사용자가 원하는 수준의 보호수준을 계산하는 단계(S50)는,
상기 오차 추정부(50)를 통해 각각 추정하는 단계(S40)에서 추정되어 제공되는 DGNSS 정보를 통해 위성별 의사거리 오차 추정치(
Figure 112017036602278-pat00143
), 위성별 비공통 오차 추정치(
Figure 112017036602278-pat00144
) 및 위성별 보정치 오차 추정치(
Figure 112017036602278-pat00145
)를 식 1에 대입하여 오차 모델(
Figure 112017036602278-pat00146
)을 구성하는 단계(S51);
상기 제공하는 단계(S20)의 획득모듈(30)에서 획득한 각 위성의 위치 및 사용자의 DGNSS 위치를 측정하는 단계(S30)의 보정 측위모듈(40)에서 측정된 사용자의 DGNSS 위치로 ENU 좌표계에서 사용자 위치로부터 위성의 방향에 대한 각 앙각(
Figure 112017036602278-pat00147
)과 방위각(
Figure 112017036602278-pat00148
)을 구하고, 식 2를 통해 위성 배치를 나타내는 단위벡터 행렬(
Figure 112017036602278-pat00149
)을 구성하는 단계(S52);
상기 오차 모델(
Figure 112017036602278-pat00150
)를 통해 구성된 위성별 오차모델(
Figure 112017036602278-pat00151
) 값을 통해 오차모델 값의 역수값(
Figure 112017036602278-pat00152
)을 식 3으로 구하고, 상기 역수값(
Figure 112017036602278-pat00153
)을 식 4에 대입하여 식 4를 통해 가중치 행렬(
Figure 112017036602278-pat00154
)를 구하는 단계(S53);
상기 단위벡터 행렬(
Figure 112017036602278-pat00155
)을 구성하는 단계(S52)를 통해 구해진 단위벡터 행렬(
Figure 112017036602278-pat00156
)과 가중치 행렬(
Figure 112017036602278-pat00157
)를 구하는 단계(S53)를 통해 구해진 가중치 행렬(
Figure 112017036602278-pat00158
)을 가중최소제곱법인 식 5에 대입하여 거리영역에서 위치영역으로 투영한 투영행렬(
Figure 112017036602278-pat00159
)을 구하는 단계(S54);
상기 투영행렬(
Figure 112017036602278-pat00160
)을 구하는 단계(S54)를 통해 구해진 투영행렬(
Figure 112017036602278-pat00161
)의 일부 요소를 취하여 상기 오차 모델(
Figure 112017036602278-pat00162
)을 구성하는 단계(S51)에서 구한 위성별 오차 모델을 곱한 후 더하는 방법으로 식 6, 식 7, 식 8식 9를 통해 동쪽 방향으로의 분산 오차(
Figure 112017036602278-pat00163
), 북쪽 방향으로의 분산 오차(
Figure 112017036602278-pat00164
), 동쪽과 북쪽 방향으로의 공분산 오차(
Figure 112017036602278-pat00165
) 및 수직 방향으로의 분산 오차(
Figure 112017036602278-pat00166
)를 구하는 단계(S55) 및
상기 분산 오차를 구하는 단계(S55)를 통해 구해진 분산 오차 요소와 수평보호수준(
Figure 112017036602278-pat00167
, Horizontal Protection Level)에 대하여 사용자가 요구하는 무결성 리스크를 만족시키는 표준점수(
Figure 112017036602278-pat00168
)와 수직보호수준(
Figure 112017036602278-pat00169
, Vertical Protection Level)에 대하여 사용자가 요구하는 무결성 리스크를 만족시키는 표준점수(
Figure 112017036602278-pat00170
)를 이용하여 수평보호수준과 수직보호수준을 식 10식 11을 통해 구하는 단계(S56)를 포함하는 것을 특징으로 하는 DGNSS 측정 위치의 무결성 검사 시스템을 이용한 무결성 검사방법.
식 1
Figure 112017036602278-pat00171

식 2
Figure 112017036602278-pat00172

식 3
Figure 112017036602278-pat00173

식 4
Figure 112017036602278-pat00174

식 5
Figure 112017036602278-pat00175

식 6
Figure 112017036602278-pat00176

식 7
Figure 112017036602278-pat00177

식 8
Figure 112017036602278-pat00178

식 9
Figure 112017036602278-pat00179

식 10
Figure 112017036602278-pat00180

수평보호수준(HPL)은 타원의 장축을 구하는 식으로 구해짐.
식 11
Figure 112017036602278-pat00181

The method of claim 3,
The step S50 of calculating the level of protection desired by the user,
The DGNSS information estimated in step S40, which is estimated through the error estimation unit 50,
Figure 112017036602278-pat00143
), Non-common error estimates by sex
Figure 112017036602278-pat00144
) And the gender-specific correction value error estimate
Figure 112017036602278-pat00145
) To the equation (1 ) to obtain the error model
Figure 112017036602278-pat00146
(S51);
The DGNSS position of the user measured by the calibration positioning module 40 of the step S30 of measuring the position of each satellite obtained by the acquisition module 30 of the providing step S20 and the DGNSS position of the user is calculated in the ENU coordinate system Each elevation angle relative to the direction of the satellite from the user location (
Figure 112017036602278-pat00147
) And azimuth angle
Figure 112017036602278-pat00148
) To obtain a unit vector matrix representing the satellite deployment through a formula (2) (
Figure 112017036602278-pat00149
(S52);
The error model (
Figure 112017036602278-pat00150
), Which is composed of
Figure 112017036602278-pat00151
) Value, the reciprocal value of the error model value (
Figure 112017036602278-pat00152
) Is obtained by the equation (3 ), and the reciprocal value
Figure 112017036602278-pat00153
), The weight matrix through formula (4) are substituted in equation (4) (
Figure 112017036602278-pat00154
(S53);
The unit vector matrix (
Figure 112017036602278-pat00155
(S52) of the unit vector matrix &lt; RTI ID = 0.0 &gt;
Figure 112017036602278-pat00156
) And a weighting matrix (
Figure 112017036602278-pat00157
(Step S53) of obtaining a weight matrix
Figure 112017036602278-pat00158
) Is substituted into the weighted least squares equation (5 ), and a projection matrix (
Figure 112017036602278-pat00159
(S54);
The projection matrix (
Figure 112017036602278-pat00160
(S54) of obtaining the projection matrix (
Figure 112017036602278-pat00161
) To obtain the error model (
Figure 112017036602278-pat00162
(S51), and then adding the result of the multiplication to the error in the east direction ( Equation 6, Equation 7, Equation 8 and Equation 9 )
Figure 112017036602278-pat00163
), The dispersion error in the north direction (
Figure 112017036602278-pat00164
), The covariance error in the east and north directions
Figure 112017036602278-pat00165
) And the dispersion error in the vertical direction (
Figure 112017036602278-pat00166
(Step S55); and
The dispersion error component obtained through the step S55 of obtaining the dispersion error and the horizontal protection level (
Figure 112017036602278-pat00167
, Horizontal Protection Level) that satisfies the user's required integrity risk.
Figure 112017036602278-pat00168
) And vertical protection level
Figure 112017036602278-pat00169
, Vertical Protection Level), which satisfies the user's required integrity risk.
Figure 112017036602278-pat00170
And a step (S56) of obtaining a horizontal protection level and a vertical protection level using Equation 10 and Equation 11 using the integrity check system of the DGNSS measurement position.
Equation 1
Figure 112017036602278-pat00171

Equation 2
Figure 112017036602278-pat00172

Equation 3
Figure 112017036602278-pat00173

Equation 4
Figure 112017036602278-pat00174

Equation 5
Figure 112017036602278-pat00175

Equation 6
Figure 112017036602278-pat00176

Equation 7
Figure 112017036602278-pat00177

Equation 8
Figure 112017036602278-pat00178

Equation 9
Figure 112017036602278-pat00179

Equation 10
Figure 112017036602278-pat00180

The horizontal protection level (HPL) is obtained by obtaining the long axis of the ellipse.
Equation 11
Figure 112017036602278-pat00181

KR1020160162790A 2016-12-01 2016-12-01 System for integrity checking of measured location of DGNSS and method for integrity checking using for thereof KR101738384B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020160162790A KR101738384B1 (en) 2016-12-01 2016-12-01 System for integrity checking of measured location of DGNSS and method for integrity checking using for thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020160162790A KR101738384B1 (en) 2016-12-01 2016-12-01 System for integrity checking of measured location of DGNSS and method for integrity checking using for thereof

Publications (1)

Publication Number Publication Date
KR101738384B1 true KR101738384B1 (en) 2017-05-22

Family

ID=59050053

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020160162790A KR101738384B1 (en) 2016-12-01 2016-12-01 System for integrity checking of measured location of DGNSS and method for integrity checking using for thereof

Country Status (1)

Country Link
KR (1) KR101738384B1 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190050554A (en) 2017-11-03 2019-05-13 세종대학교산학협력단 Multi-constellation gnss positioning system and method by correcting the inter-system time difference
KR101970239B1 (en) * 2017-12-18 2019-08-27 한국과학기술원 Method and System of Optimal Protection Level for Local-Area Differential GNSS to Support UAV Navigation
CN112596080A (en) * 2020-11-25 2021-04-02 中国人民解放军93216部队 Method for testing integrity index of differential Beidou take-off and landing guide system of unmanned aerial vehicle
CN113302452A (en) * 2019-01-16 2021-08-24 罗伯特·博世有限公司 Method for providing an integrity range for parameter estimation
WO2022033440A1 (en) * 2020-08-11 2022-02-17 维沃移动通信有限公司 Integrity protection method and system
KR20230056895A (en) * 2021-10-21 2023-04-28 (주)답스 Pseudo-range measurement and apparatus about Continuous wave signal in the MF R-mode System
KR102645996B1 (en) * 2022-11-02 2024-03-12 한국해양과학기술원 Geospatial information platform service system and method for maritime safety based on smart devices

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006112822A (en) 2004-10-12 2006-04-27 Mitsubishi Electric Corp Mobile station, dgps basis station, dgps center station and dgps positioning system
KR101151670B1 (en) 2011-12-14 2012-06-08 한국항공우주연구원 Standard deviation establishing method for gnss pseudorange correction considering reference receiver clock adjustment effect
KR101477041B1 (en) 2013-08-13 2015-01-02 한국해양과학기술원 Satellite Signal Anomaly Monitoring System for DGNSS Reference Station and Its Monitoring Method
JP5778036B2 (en) 2008-11-17 2015-09-16 クゥアルコム・インコーポレイテッドQualcomm Incorporated DGNSS correction for positioning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006112822A (en) 2004-10-12 2006-04-27 Mitsubishi Electric Corp Mobile station, dgps basis station, dgps center station and dgps positioning system
JP5778036B2 (en) 2008-11-17 2015-09-16 クゥアルコム・インコーポレイテッドQualcomm Incorporated DGNSS correction for positioning
KR101151670B1 (en) 2011-12-14 2012-06-08 한국항공우주연구원 Standard deviation establishing method for gnss pseudorange correction considering reference receiver clock adjustment effect
KR101477041B1 (en) 2013-08-13 2015-01-02 한국해양과학기술원 Satellite Signal Anomaly Monitoring System for DGNSS Reference Station and Its Monitoring Method

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190050554A (en) 2017-11-03 2019-05-13 세종대학교산학협력단 Multi-constellation gnss positioning system and method by correcting the inter-system time difference
KR101970239B1 (en) * 2017-12-18 2019-08-27 한국과학기술원 Method and System of Optimal Protection Level for Local-Area Differential GNSS to Support UAV Navigation
CN113302452A (en) * 2019-01-16 2021-08-24 罗伯特·博世有限公司 Method for providing an integrity range for parameter estimation
WO2022033440A1 (en) * 2020-08-11 2022-02-17 维沃移动通信有限公司 Integrity protection method and system
CN114125846A (en) * 2020-08-11 2022-03-01 维沃移动通信有限公司 Integrity protection method and system
CN114125846B (en) * 2020-08-11 2023-09-12 维沃移动通信有限公司 Integrity protection method and system
CN112596080A (en) * 2020-11-25 2021-04-02 中国人民解放军93216部队 Method for testing integrity index of differential Beidou take-off and landing guide system of unmanned aerial vehicle
CN112596080B (en) * 2020-11-25 2024-04-16 中国人民解放军93216部队 Method for testing integrity index of unmanned aerial vehicle differential Beidou lifting guide system
KR20230056895A (en) * 2021-10-21 2023-04-28 (주)답스 Pseudo-range measurement and apparatus about Continuous wave signal in the MF R-mode System
KR102583747B1 (en) 2021-10-21 2023-09-27 (주)답스 Pseudo-range measurement and apparatus about Continuous wave signal in the MF R-mode System
KR102645996B1 (en) * 2022-11-02 2024-03-12 한국해양과학기술원 Geospatial information platform service system and method for maritime safety based on smart devices

Similar Documents

Publication Publication Date Title
KR101738384B1 (en) System for integrity checking of measured location of DGNSS and method for integrity checking using for thereof
US8014948B2 (en) Navigation system with apparatus for detecting accuracy failures
US8976064B2 (en) Systems and methods for solution separation for ground-augmented multi-constellation terminal area navigation and precision approach guidance
US8742984B2 (en) Integrity method for differential corrections
US10739471B2 (en) GNSS receiver with a capability to resolve ambiguities using an uncombined formulation
CN110007317B (en) Star-selection optimized advanced receiver autonomous integrity monitoring method
KR101761782B1 (en) System and method for determining the protection level
Angrisano et al. GNSS Reliability Testing in Signal-Degraded Scenario.
CN111007552B (en) LDACS-based air-ground cooperative positioning and integrity monitoring method
US20110231038A1 (en) Aircraft landing system using relative gnss
JP4723932B2 (en) Positioning system
CA2891363A1 (en) Method for estimating the level of error in satellite geolocation measurements and for monitoring the reliability of said estimations and associated device
CN110208835A (en) A kind of cross-system tight integration Differential positioning method based on iono-free combination
KR101874974B1 (en) Apparatus and method for generating differential global navigation satellite system pseudo range correction information
CN103983986A (en) Particle filter-based improved RAIM (Receiver Autonomous Integrity Monitoring) anti-deception jamming method
US20190004181A1 (en) Method and apparatus applicable to positioning in nlos environment
CN114235007A (en) Method and system for positioning and integrity monitoring of APNT service
JP5566598B2 (en) Navigation system comprising a device for detecting inaccuracy
RU2577846C1 (en) Method of determining integrity of high-precision navigation determinations of consumer and system therefor
Weinbach et al. Integrity of the Trimble® centerpoint RTX correction service
EP2367023A1 (en) Aircraft landing system using relative GNSS
KR20150110398A (en) Geopositioning method with trust index, and associated terminal
KR101040053B1 (en) An integrity monitoring method to detect and identify the gnss satellite clock anomaly by monitoring the pseudorange measurement noise
KR101040054B1 (en) An integrity monitoring method to detect and identify the gnss satellite clock anomaly by monitoring the receiver clock
JP2009175134A (en) Navigation system with apparatus for detecting accuracy failure

Legal Events

Date Code Title Description
A201 Request for examination
GRNT Written decision to grant