KR101290085B1 - Method for monitoring tropospheric delay irregularity in multi-reference stations environment and system using monitoring method - Google Patents

Abstract

PURPOSE: A method for monitoring delay anomaly phenomenon in the convection zone is provided to detect a sudden convection zone delay anomaly phenomenon by comparing the differences and the threshold values between the reference stations as the method uses weather parameters from each reference station as the key parameter necessary for monitoring delay anomaly phenomenon. CONSTITUTION: A method for monitoring delay anomaly phenomenon in multiple reference stations comprises the steps of: collecting weather parameters in multiple reference stations (S101), setting up a low performance point (S102), producing an delay estimation value (eTztd) in the convection zone ceiling (S103) , producing an delay calculation value (cTztd) in the convection zone ceiling by using the collected weather parameters (S104), selecting the maximum value as the delta delay value in the convection zone ceiling (ΔTztd) among the difference values from the delay calculation values (cTztd) and the delay estimation values (eTztd) (S105), producing a delta delay value in the delta convection zone in line-of-sight (S106), and comparing a delta delay value in the delta convection zone in line-of-sight and the setting threshold value. [Reference numerals] (AA) Start; (BB) No; (CC) Yes; (DD) End; (S101) Collecting weather parameters in multiple reference stations; (S102) Setting low performance point; (S103) Calculating estimated convection zone ceiling delay calculation value (eTztd); (S104) Calculating convection zone ceiling delay calculation value (cTztd); (S105) Calculating delta convection zone ceiling delay calculation value (ΔTztd); (S106) Calculating delta convection zone in line-of-sight delay calculation value (ΔTsd); (S108) Output flag'1'; (S109) Output flag'0'

Description

METHOD FOR MONITORING TROPOSPHERIC DELAY IRREGULARITY IN MULTI-REFERENCE STATIONS ENVIRONMENT AND SYSTEM USING MONITORING METHOD}

The present invention relates to a method for monitoring a convective delay anomaly and a system using the same. More particularly, the present invention relates to the correction of information due to an unbalanced convective environment between reference stations when generating correction information using a network RTK (Real Time Kinematic). The present invention relates to a method for monitoring a convective delay anomaly in a multi-reference station environment capable of improving reliability of correction information by detecting a satellite deteriorating and providing a detected satellite to a user.

Global Positioning System (GPS) is a global radio navigation satellite system developed and promoted by the US Department of Defense (DOD), which is essential for the safe operation of land, marine and aviation.

The GPS is also referred to as NAVSTAR / GPS in the sense of a system using a navigation system with Time And Ranging (NAVSTAR).

The GPS consists of a total of 24 navigating satellites, each of which was launched on six circular orbits with an altitude of approximately 20,000 km, a period of about 12 hours, and an orbital inclination of 55 degrees, a ground control station that manages satellites, and a user's mobile station.

Service interruption due to the failure of the Global Navigation Satellite System (GNSS) can disrupt transportation services in the ocean, air and land, resulting in economic losses. Therefore, IM (Integrity Monitoring) needs to be executed for the user during the GNSS service.

In other words, the factors that can cause an abnormality in the GNSS service include the standby factors such as satellite clock error, satellite orbit error, and navigation data error, and the signal waiting such as multipath error, convective delay and ionospheric delay error. It can be classified as a factor due to irregularities in the process of passing through.

In relation to the prior art of receiving a signal from a GPS satellite and correcting a position error, Korean Patent Laid-Open No. 10-2008-0065040 relates to an apparatus and a method for satellite navigation signal correction, and includes an α-β parameter included in a satellite navigation message. A method of minimizing the cost value by calculating the ionospheric delay of each satellite and measuring the ionospheric delay from the actual satellite navigation signal by using the difference between the ionospheric delay calculation value and the measured value of each satellite.

However, local rainfall and typhoons can cause weather differences between reference countries. As shown in FIG. 1, when the gas phase of the convective layer, in particular, the water vapor pressure is irregular, a difference in the amount of delay in the convective ceiling of the common satellite between reference stations may occur.

Therefore, there is a problem in that the convective delay anomaly occurs due to the weather difference of the convective layers between the reference stations, thereby reducing the quality of the correction information.

SUMMARY OF THE INVENTION The present invention has been made to solve the above-mentioned problems, and it is possible to detect a sudden convective delay anomaly between reference stations through comparison between the reference stations and thresholds of the main weather parameters necessary for monitoring the convective delay anomaly. To provide a method for monitoring convective delay anomalies in multiple base station environments and a system using the same.

In addition, the present invention provides a user with information on the detected satellites by detecting satellites causing sudden convective delay anomaly between reference stations, thereby improving the reliability of correction information due to the convective delay. It aims to provide a method of monitoring delay anomaly and a system using the same.

However, the object of the present invention is not limited to the above-mentioned objects, and other objects not mentioned can be clearly understood by those skilled in the art from the following description.

In order to achieve the above object, in the multiple reference station environment according to the present invention, the method of monitoring the convective delay anomaly, each weather signal receiving a satellite signal from a common satellite consisting of a plurality of satellites, and measuring weather parameters, respectively A first step of collecting weather parameters in a multiple reference station comprising a plurality of reference stations, and a low performance point (LPP) which is a point located farthest from four reference stations among the multiple reference stations A second step of calculating a convective ceiling delay amount eTztd in proportion to the distance between the first reference station selected from the multiple reference stations and the second reference station in the LPP, the first reference station; A fourth step of calculating the convective ceiling delay amount cTztd using the weather parameters collected at the second reference station, respectively; The difference between the estimated convective ceiling delay amount eTztd and the calculated convective ceiling delay amount cTztd calculated in the fourth step is calculated, and the maximum value of the calculated difference values is the delta convective ceiling delay amount ( A fifth step of selecting ΔTztd), and a sixth step and a convective layer of calculating the delta convective bed direction delay amount projected in the eye direction using a mapping function of the maximum delta convective ceiling delay amount selected in the fifth step. And a seventh step of outputting the presence or absence of convective delay anomalies for each reference station by comparing the delayed anomaly determination unit with the set threshold value calculated in the sixth step. do.

In addition, the method of monitoring the convective delay anomaly in a multiple reference station environment according to the present invention is characterized in that, in the second step, LPP is established through the following relational expression.

Here, n represents the number of reference stations, and x, y, and z represent coordinates.

In addition, the method for monitoring the convective delay anomaly in a multiple reference station environment according to the present invention, the weather parameters collected in the first step, the pressure (P), temperature (T), humidity (H) information It is characterized by.

In addition, the method of monitoring the convective delay anomaly in a multiple reference station environment according to the present invention is characterized in that the first reference station and the second reference station selected in the third step are located on a line near the LPP. It is done.

In addition, the method of monitoring the convective delay anomaly in a multi-reference station environment according to the present invention, in the fourth step, the calculated convective ceiling delay amount (cTztd) is the wet delay amount (ZHD, Zenith Hydrostatic Delay) and the drying delay amount. (ZWD, Zenith Wet Delay) is calculated.

In addition, the method for monitoring the convective delay anomaly in a multiple reference station environment according to the present invention is characterized in that the wet delay amount (ZHD) is calculated through the following relational expression.

Here, P ₀ represents the barometric pressure [hPa], φ the latitude [deg], and H the altitude [km] of the sensor provided in the reference station.

In addition, the method of monitoring the convective delay anomaly in a multiple reference station environment according to the present invention is characterized in that the drying delay amount (ZWD) is calculated through the following relationship.

Here, e ₀ represents the water vapor pressure [hPa] and T ₀ represents the absolute temperature [K].

In addition, the method of monitoring the convective delay anomaly in a multiple reference station environment according to the present invention is characterized in that, in the third step, the convective ceiling delay amount eTztd is calculated through the following relational expression.

Where Tztd ₁ is the convective ceiling delay amount of the first reference station, Tztd ₂ is the convective ceiling delay amount of the second reference station, d1 is the distance between the LPP and the first reference station, d2 is the LPP and the second reference station. Indicates the distance of.

In addition, in the fifth step, the method for monitoring the convective delay anomaly in the multiple reference station environment according to the present invention is characterized in that the delta convective ceiling delay amount ΔTztd is obtained through the following relational expression.

Here, cTztd ₁ represents the calculated convective ceiling delay amount using the weather parameter of the first reference station at the LPP point, and cTztd ₂ represents the calculated convective ceiling delay amount using the weather parameter of the second reference station at the LPP point.

In addition, in the sixth step, the method for monitoring the convective delay anomaly in the multiple reference station environment according to the present invention is characterized in that the delta convective line direction delay amount ΔTsd is calculated through the following relational expression.

ΔTztd is the delta convective ceiling delay, m _i (E) is the mapping function, and E is the satellite angle of the satellite, which is common to both reference stations near the LPP and the straight line.

Further, in the multi-reference station environment monitoring method according to the present invention, in the seventh step, the convective delay anomaly determination unit outputs a flag indicating the presence or absence of convective delay anomaly. The output flag is defined by the following relational expression.

Here, flag 1 represents an abnormal state and flag 0 represents a normal state.

In the multiple reference station environment according to the present invention, a convective delay anomaly monitoring system comprises: a multiple reference station comprising a plurality of reference stations each receiving a GPS signal from a common satellite consisting of a plurality of satellites, and having a weather sensor; A convective delay anomaly monitoring apparatus for monitoring a convective delay anomaly due to meteorological differences between reference stations, wherein the convective delay anomaly monitoring apparatus includes a weather parameter for collecting weather parameters of each reference station from the weather sensor. An LPP setting unit for setting a LPP (Low Performance Point) which is a point located farthest from four reference stations among the multiple reference stations, and a first reference station selected from the multiple reference stations in the LPP; Convective ceiling delay amount that calculates the convective ceiling delay amount eTztd in proportion to the distance of the second reference station A convective ceiling delay calculation unit for calculating a convective ceiling delay amount (cTztd) at the LPP point using the weather parameters collected from the government, the first reference station and the second reference station, and the convective ceiling Calculate the difference between the delay estimate (eTztd) and the convection ceiling delay calculation (cTztd), and calculate the delta convection ceiling delay that selects the maximum value as the delta convection ceiling delay (ΔTztd) Delta convective line direction delay calculation unit and delta convective line direction delay amount and setting for calculating the delta convective visual direction delay projected in the eye direction using the mapping function with the maximum delta convective ceiling delay amount And a convective delay anomaly determination unit for outputting the presence or absence of a convective delay anomaly for each reference station by comparing the threshold values.

According to the present invention, there is provided a recording medium having recorded thereon a program for executing a method for monitoring a convective delay anomaly in a multi-reference station environment as described above.

According to the method of monitoring the convective delay anomaly in a multi-reference station environment of the present invention and a system using the same, the difference between the reference stations and the threshold value is used by using the weather parameters of each reference station as the main parameters required for monitoring the convective delay anomaly. Compared with, there is an advantage that can detect sudden convective delay anomaly between reference stations.

In addition, according to the method for monitoring the convective delay anomaly in a multi-reference station environment of the present invention and a system using the same, the satellites causing the rapid convective delay anomaly between the reference stations are identified and the detected satellite information is provided to the user. As a result, the reliability of correction information due to convective delay anomaly may be improved.

1 is an exemplary diagram illustrating a difference in convective ceiling delay of a common satellite due to a weather difference between reference stations.
2 is a block diagram schematically illustrating a monitoring system for monitoring a convective delay anomaly in a multiple reference station environment according to the present invention.
FIG. 3 is an exemplary diagram showing a Low Performance Point (LPP), which is a point located farthest from four reference stations in a multiple reference station environment.
4 is a block diagram showing the configuration of the convective layer delay anomaly monitoring device constituting the monitoring system according to the present invention.
5 is a flowchart illustrating a method for monitoring a convective delay anomaly in a multiple reference station environment in accordance with the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, a detailed description of preferred embodiments of the present invention will be given with reference to the accompanying drawings. In the following description of the present invention, detailed description of known functions and configurations incorporated herein will be omitted when it may make the subject matter of the present invention rather unclear.

Embodiments in accordance with the concepts of the present invention can make various changes and have various forms, so that specific embodiments are illustrated in the drawings and described in detail in this specification or application. However, this is not intended to limit the embodiments in accordance with the concept of the present invention to a particular disclosed form, it should be understood to include all changes, equivalents, and substitutes included in the spirit and scope of the present invention.

It is to be understood that when an element is referred to as being "connected" or "connected" to another element, it may be directly connected or connected to the other element, . On the other hand, when an element is referred to as being "directly connected" or "directly connected" to another element, it should be understood that there are no other elements in between. Other expressions that describe the relationship between components, such as "between" and "between" or "neighboring to" and "directly adjacent to" should be interpreted as well.

The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. Singular expressions include plural expressions unless the context clearly indicates otherwise. In this specification, the terms "comprises ",or" having ", or the like, specify that there is a stated feature, number, step, operation, , Steps, operations, components, parts, or combinations thereof, as a matter of principle.

2 is a block diagram schematically illustrating a monitoring system for monitoring a convective delay anomaly in a multi-reference station environment according to the present invention, and FIG. 3 is located at the longest distance from four reference stations in a multi-reference station environment. 4 is a block diagram illustrating a convective delay anomaly monitoring device constituting a monitoring system according to the present invention.

As shown, the system 100 for monitoring the convective delay anomaly in a multi-reference station environment of the present invention comprises a common satellite (10) consisting of a plurality of satellites (10-1 to 10-n) for transmitting GPS signals. And a weather sensor 21 for receiving GPS signals transmitted from the common satellite 10 and measuring weather parameters such as pressure (P), temperature (T), and humidity (H), respectively. Convective delay anomaly monitoring the presence or absence of convective delay anomalies for GPS signals received from multiple reference stations 20 and multiple reference stations 20 consisting of a plurality of reference stations 20-1 to 20-n. It may be configured as a phenomenon monitoring device 30.

In such a multiple reference station environment for monitoring the convective delay anomaly, as shown in FIG. 3, the first reference station 20-1, the second reference station 20-2, and the first reference station are shown as four reference stations. A low performance point (LPP) 40, which is a point located farthest from the third reference station 20-3 and the fourth reference station 20-4, can be set.

The LPP 40 may be located on a straight line or near a straight line between the first reference station 20-1 and the second reference station 20-2.

As shown in FIG. 4, the convective delay anomaly monitoring apparatus 30 includes a weather parameter collection unit 310, a low performance point (LPP) setting unit 320, a convective ceiling delay amount estimation unit 330, The convective ceiling delay amount calculating unit 340, the delta convective ceiling delay amount calculating unit 350, the delta convective layer visual delay amount calculating unit 360, and the convective layer delay anomaly determining unit 370 may be configured. have.

The meteorological parameter collection unit 310 is a meteorological parameter required for calculating the convective layer delay at each reference station 20-1 to 20-n measured from the weather sensor 21 provided in each of the multiple reference stations 20. For example, pressure (P), temperature (T) and humidity (H), etc.) can be collected.

The LPP setting unit 320 sets the LPP 40, which is a point located farthest from the four reference stations among the multiple reference stations 20, and the convective ceiling delay amount estimator 330 is the LPP 40. In Equation 20, an estimated convective ceiling delay amount eTztd may be calculated in proportion to the distance between the first reference station 20-1 and the second reference station 20-2 selected from the multiple reference stations 20.

The convective ceiling delay calculation unit 340 uses the weather parameters collected through the first reference station 20-1 and the second reference station 20-2 to delay the convective ceiling at the LPP 40 point. The amount calculation value cTztd can be calculated respectively.

The delta convective ceiling delay calculation unit 350 calculates the convective ceiling delay amount estimate eTztd calculated by the convective ceiling delay amount estimator 330 and the convective ceiling delay amount calculation unit 340. After calculating the difference value of the ceiling delay amount cTztd, the maximum value may be selected from the calculated difference values.

The delta convection visual delay amount calculating unit 360 projects the maximum delta convective ceiling delay amount selected by the delta convective ceiling delay calculation unit 350 in a visual direction using a mapping function. The layer line direction delay amount ΔTsd can be calculated.

In addition, the convective layer delay anomaly determination unit 370 determines whether there is a convective layer delay anomaly for the multiple reference stations 20 by comparing the delta convective line direction delay amount ΔTsd and the set threshold value Th. In case of a normal state, a flag of "0" and an abnormal state of "1" can be generated and output.

Therefore, when generating correction information using network RTK, by identifying satellites with irregular convective delay due to meteorological differences between reference stations, it is possible to provide users with information on the quality of correction information to improve satisfaction with correction information. have.

5 is a flowchart illustrating a method for monitoring a convective delay anomaly in a multiple reference station environment according to the present invention.

2 to 5, first, a plurality of reference stations 20 forming a multiple reference station 20 for transmitting and receiving satellite signals to and from a common satellite 10 including a plurality of satellites 10-1 to 10-n. In the weather sensor 21 provided in each of -1 to 20-n, information on weather parameters, for example, pressure (P), temperature (T), and humidity (H), is measured in real time, and the weather parameter collection unit ( In 320, the real-time measured weather parameters are collected (S101).

The low performance point (LPP) setting unit 320 sets the LPP 40, which is a point located farthest from four reference stations among the multiple reference stations 20 (S102). In this case, the LPP 40 may be set through the following relational expression.

[Formula 1]

Here, n represents the number of reference stations, and x, y, z represent coordinates. The LPP 40 may be located in a straight line between the first reference station 20-1 and the second reference station 20-2 of the four reference stations.

Subsequently, through the convective ceiling delay amount estimator 330, the LPP 40 estimates a convective ceiling delay amount in proportion to the distance between the first reference station 20-1 and the second reference station 20-2. eTztd) is calculated through the following relational expression (S103).

[Formula 2]

Where Tztd ₁ is the convective ceiling delay amount of the first reference station, Tztd ₂ is the convective ceiling delay amount of the second reference station, d1 is the distance between the LPP and the first reference station, d2 is the LPP and the second reference station. Indicates the distance of.

The convective ceiling delay amount calculation unit 340 calculates the convective ceiling delay amount cTztd based on the sum of the wet delay amount (ZHD, Zenith Hydrostatic Delay) and the dry delay amount (ZWD, Zenith Wet Delay) (S104). .

At this time, the wet delay amount (ZHD) is calculated through the following relationship,

[Equation 3]

Here, P ₀ represents the barometric pressure [hPa], φ the latitude [deg], and H the altitude [km] of the sensor provided in the reference station.

In addition, the drying delay amount ZWD may be calculated through the following relationship.

[Formula 4]

Here, e ₀ represents the water vapor pressure [hPa] and T ₀ represents the absolute temperature [K].

Next, the delta convective ceiling delay calculation unit 350 calculates a difference value between the convective ceiling delay amount estimated value eTztd and the convective ceiling delay amount calculated value cTztd, as shown in the following relational expression. The maximum value of the calculated difference values is selected as the delta convective ceiling delay amount ΔTztd (S105).

[Formula 5]

Here, cTztd ₁ represents the calculated convective ceiling delay amount using the weather parameter of the first reference station at the LPP point, and cTztd ₂ represents the calculated convective ceiling delay amount using the weather parameter of the second reference station at the LPP point.

The delta convex line direction delay amount calculating unit 360 calculates the delta convective line direction delay amount ΔTsd through the following relational expression (S106).

[Formula 6]

E denotes the satellite's bottom angle, which is common to two reference stations forming the closest straight line of the LPP. That is, the amount of delta convection line direction direction projected to the elevation angle E may be calculated by using the mapping function of the delta convective ceiling delay ΔTztd. In this case, the mapping function m _i (E) may be represented in the form of soft fractions.

Subsequently, the convective layer delay anomaly determination unit 370 compares the delta convective layer line direction delay amount ΔTsd calculated by the delta convective bed line direction delay amount calculation unit 360 and compares the preset threshold Th. It is determined whether there is an error (S107), and a flag is generated as shown in the following relational expression.

[Equation 7]

That is, when the visual direction delay amount ΔTsd exceeds the set threshold Th, the convective delay anomaly determination unit 370 determines that the abnormal state is an abnormal state and generates the abnormal flag “1” to determine the PRN number of the abnormal satellite. If the visual direction delay amount ΔTsd is equal to or less than the set threshold Th, the output signal is determined together (S108), and the normal flag "0" is generated and output (S109).

As described above, by identifying the satellite causing the sudden convective delay anomaly between the reference stations due to the weather difference and providing the user with the detected satellite information, the correction information reliability due to the convective delay anomaly is generated when the correction information is generated. There are features that can be improved.

Convective delay anomalies can be affected by meteorological differences in specific reference stations between base stations in a multi-reference station environment, and there is a need to improve real-time monitoring of convective delay anomalies. In contrast, the present invention provides a monitoring method and system for the detection of convective delay anomalies based on multiple reference stations.

In addition, the method of monitoring the convective delay anomaly in a multiple reference station environment according to the present invention may be implemented in the form of program instructions that can be executed by various computer means and recorded in a computer readable medium. Computer-readable media may include program instructions, data files, data structures, etc. alone or in combination. Program instructions recorded on the media may be those specially designed and constructed for the present invention, or they may be of the kind well-known and available to those having skill in the computer software arts.

Examples of computer-readable recording media include magnetic media such as hard disks, floppy disks, and magnetic tape, optical recording media such as CD-ROMs, DVDs, magnetic-optical media such as floppy disks, and ROM, RAM, flash memory, and the like. Hardware devices specifically configured to store and execute program instructions are included.

Examples of program instructions include not only machine code generated by a compiler, but also high-level language code that can be executed by a computer using an interpreter or the like. Hardware devices may be configured to operate as one or more software modules to perform the operations of the present invention, and vice versa.

While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is to be understood that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims. will be. Therefore, the true technical protection scope of the present invention will be defined by the technical spirit of the appended claims.

10: common satellite
20: multiple reference stations
21: weather sensor
30: convective delay anomaly monitoring device
40: Low Performance Point (LPP)
100: convective delay anomaly monitoring system
310: weather parameter collection unit
320: LPP setting section
330: convection ceiling delay amount estimation unit
340: convective ceiling delay amount calculation unit
350: delta convection ceiling delay calculation unit
360: Delta convection layer visual delay delay calculation unit
370: convective delay anomaly determination unit

Claims (20)

Receiving a satellite signal from a common satellite consisting of a plurality of satellites and collecting weather parameters from a plurality of reference stations each comprising a plurality of reference stations each having a weather sensor for measuring weather parameters;
A second step of establishing a low performance point (LPP), which is a point located farthest from four reference stations among the multiple reference stations;
A third step of calculating a convective ceiling delay amount (eTztd) in proportion to the distance between the first reference station and the second reference station selected from the multiple reference stations in the LPP;
A fourth step of calculating a convective ceiling delay amount (cTztd) using the weather parameters collected at the first reference station and the second reference station, respectively;
The difference between the estimated convective ceiling delay amount eTztd calculated in the third step and the calculated convective ceiling delay amount cTztd calculated in the fourth step is calculated, and the maximum value of the calculated difference values is delta. A fifth step of selecting the convective ceiling delay amount [Delta] Tztd;
A sixth step of calculating a delta convective bed line direction delay amount projected in the line of sight direction using a mapping function of the maximum delta convective ceiling delay amount selected in the fifth step; And
A seventh step of outputting the presence or absence of convective delay anomaly for each reference station by comparing the delta convective line direction delay amount calculated in the sixth step with a set threshold value; A method of monitoring convective delay anomalies in a multi-reference station environment.
The method of claim 1,
The second step is a convective delay anomaly monitoring method in a multiple reference station environment, characterized in that for setting the LPP through the following relationship.

(Where n is the number of reference stations and x, y, z are the coordinates)
The method of claim 1,
The meteorological parameter collected in the first step includes pressure (P), temperature (T), and humidity (H) information.
The method of claim 1,
The first reference station and the second reference station selected in the third step is located on the line close to the straight line to the LPP convective delay anomaly monitoring method in a multiple reference station environment.
The method of claim 1,
In the third step,
Convective ceiling delay estimate (eTztd) is calculated by the following equation.

Where Tztd ₁ is the convective ceiling delay of the first reference station, Tztd ₂ is the convective ceiling delay of the second reference station, d1 is the distance between the LPP and the first reference station, and d2 is the LPP and second reference. Indicates the distance to the station)
The method of claim 1,
In the fourth step, the convective ceiling delay amount (cTztd) is calculated as the sum of the wet delay amount (ZHD, Zenith Hydrostatic Delay) and the drying delay amount (ZWD, Zenith Wet Delay). A method of monitoring convective delay anomalies.
The method according to claim 6,
Wet delay (ZHD) is a method of monitoring a convective delay anomaly in a multiple reference station environment, characterized by the following relationship.

Where P ₀ is the barometric pressure [hPa], φ is the latitude [deg], and H is the altitude of the sensor in the reference station [km].
The method according to claim 6,
Drying delay (ZWD) is a method of monitoring a convective delay anomaly in a multiple reference station environment, characterized by the following relationship.

(Where e0 represents water vapor pressure [hPa] and T0 represents absolute temperature [K])
The method of claim 1,
In the fifth step,
Delta convective ceiling delay amount (ΔTztd) is obtained by the following relationship.

(Where cTztd ₁ is the convective ceiling delay calculation using the weather parameters of the first reference station at the LPP point, and cTztd ₂ is the convection ceiling delay calculation using the weather parameters of the second reference station at the LPP point)
The method of claim 1,
In the sixth step,
Delta convective line direction delay amount ΔTsd is calculated by the following relationship.

(Where ΔTztd is the delta convective ceiling delay, m _i (E) is the mapping function, and E is the lower angle of the satellite commonly seen in LPP and two reference stations near the straight line)
The method of claim 1,
In the seventh step,
The convective delay anomaly determination unit outputs a flag indicating the presence or absence of the convective delay anomaly, and the output flag is defined by the following relational expression. Monitoring method.

(Where flag 1 is anomalous and flag 0 is normal)
Multiple reference stations each receiving a GPS signal from a common satellite consisting of a plurality of satellites, the plurality of reference stations comprising a weather sensor;
A convective delay anomaly monitoring device for monitoring the convective delay anomaly due to meteorological differences between the reference stations;
The convective layer delay anomaly monitoring device,
A weather parameter collection unit for collecting weather parameters of each reference station from the weather sensor;
An LPP setting unit for setting an LPP (Low Performance Point), which is a point located farthest from four reference stations among the multiple reference stations;
A convective ceiling delay amount estimator for estimating a convective ceiling delay amount eTztd in proportion to the distance between a first reference station selected from the multiple reference stations and a second reference station;
A convective ceiling delay calculation unit for calculating a convective ceiling delay amount cTztd at the LPP point using the weather parameters collected at the first reference station and the second reference station, respectively;
A delta convection ceiling that calculates the difference between the convection ceiling delay estimate (eTztd) and the convection ceiling delay calculation (cTztd) and selects the maximum value of the calculated difference as the delta convection ceiling delay (ΔTztd). A delay calculation unit;
A delta convective bed line direction delay calculation unit configured to calculate a delta convective bed line direction delay amount projected in the line of sight using a mapping function of the maximum delta convective ceiling delay; And
Convection layer in a multi-reference station environment comprising a; convective delay anomaly determination unit for outputting the presence or absence of convective delay anomaly for each reference station by comparing the delta convection line direction delay amount and the set threshold value. A system for monitoring delay anomalies.
13. The method of claim 12,
The LPP setting unit, the system for monitoring the convective delay anomaly in a multiple reference station environment, characterized in that for setting the LPP through the following relationship.

(Where n is the number of reference stations and x, y, z are the coordinates)
13. The method of claim 12,
The meteorological parameter is a system for monitoring the convective delay anomaly in a multiple reference station environment, characterized in that the pressure (P), temperature (T), humidity (H) information.
13. The method of claim 12,
The convective ceiling delay amount estimator may calculate a convective ceiling delay amount (cTztd) based on a sum of a wet delay amount (ZHD, Zenith Hydrostatic Delay) and a dry delay amount (ZWD, Zenith Wet Delay). A system for monitoring convective delay anomalies in a reference station environment.
13. The method of claim 12,
The convective ceiling delay amount estimating unit calculates the convective ceiling delay amount (eTztd) through the following relational expression.

Where Tztd ₁ is the convective ceiling delay of the first reference station, Tztd ₂ is the convective ceiling delay of the second reference station, d1 is the distance between the LPP and the first reference station, and d2 is the LPP and second reference. Indicates the distance to the station)
13. The method of claim 12,
And the delta convective ceiling delay calculation unit calculates a delta convective ceiling delay amount ΔTztd through the following relational expression.

(Where cTztd ₁ is the convective ceiling delay calculation using the weather parameters of the first reference station at the LPP point, and cTztd ₂ is the convection ceiling delay calculation using the weather parameters of the second reference station at the LPP point)
13. The method of claim 12,
And the delta convective visual delay amount calculating unit calculates a delta convective visual direction delay amount ΔTsd through the following relational expression.

(Where ΔTztd is the delta convective ceiling delay, m _i (E) is the mapping function, and E is the lower angle of the satellite commonly seen in LPP and two reference stations near the straight line)
13. The method of claim 12,
The convective delay anomaly determination unit outputs a flag indicating the presence or absence of the convective delay anomaly, and the output flag is defined by the following relational expression. System to monitor them.

(Where flag 1 is anomalous and flag 0 is normal)
12. A computer-readable recording medium having recorded thereon a computer program for executing the method of monitoring the convective delay anomaly in the multi-reference station environment according to any one of claims 1 to 11.

Images