KR101407008B1 - System and method for processing GNSS data automatically - Google Patents

Abstract

The present invention relates to a system and method for automatically processing GNSS data, including a satellite orbit acquisition definition unit, a GNSS observation station definition unit, a processing condition definition unit, a satellite orbit information data unit, a GNSS observation data unit, A database; A GNSS data acquisition program, a GNSS data acquisition program, a GNSS data acquisition program, a GNSS data acquisition program, a GNSS data acquisition program, and a GNSS data processing program.

Description

Technical Field [0001] The present invention relates to a system and method for automatically processing data,

The present invention relates to a system and method for automatically processing GNSS data, and more particularly, to a system and method for automatically processing GNSS data, And directions and the like, and displays it on an online map, thereby intuitively recognizing the crustal variation information.

 In general, the Global Navigation Satellite System (GNSS) is a satellite navigation positioning system that is used extensively from vehicle navigation systems used by the general public to geodetic surveying, surveying, and crustal monitoring using high performance GNSS receivers.

GNSS stations (or GNSS stations) are distributed throughout the country, and GNSS data is collected and stored in real time for various purposes such as surveying, monitoring of crustal movement, weather utilization and so on.

The GNSS observation data can be used to monitor the crustal variation with the coordinates obtained through the separate processing. We use a high-precision GNSS data processing program such as Bernese GNSS Software (Bernese), GIPSY or GAMIT.

 However, if the precise coordinates of each GNSS station are calculated and accumulated over a long period of time using the above program, the magnitude and direction of the crustal fluctuation at each GNSS station can be obtained.

However, it is very troublesome to process GNSS data by manpower over a long period of time.

In addition, the high-precision GNSS data processing program is running on a LINUX operating system that is not familiar to the general public, and it is very difficult for the general public to use because of the unfamiliar and difficult concept.

 Currently, programs for automatic acquisition of GNSS data and processing of high-precision data have been developed and operated on the basis of LINUX OS, but systems and programs for collecting, managing and processing data based on Windows have not been developed yet, There is a demand for a high-precision GNSS data processing system that can easily and conveniently process GNSS data.

SUMMARY OF THE INVENTION The present invention has been made to solve the conventional problems as described above, and its object is to automate the process of acquiring highly precise GNSS data based on Windows and automating high-precision data processing, not only systematically managing GNSS data, The present invention provides a system and method for automatic processing of GNSS data, which can easily calculate crustal fluctuation amount, speed and direction, and display it on an online map to intuitively recognize crustal variation information.

The system includes a satellite GNSS data processing system database, a data processing unit, and a control unit. The database includes a satellite orbit acquisition definition unit for defining a precise orbital force of the GNSS satellite, a GNSS station definition unit for defining the GNSS station information, A processing condition defining unit for defining setting conditions for processing GNSS data in the data processing unit, a satellite orbit information data unit for storing the orbit information defined and acquired by the satellite orbit acquisition defining unit, A GNSS observation data section in which GNSS observation data measured from an observing station defined in a definition section is stored, and a data processing result section that stores coordinates and a displacement of a GNSS observing point where processing after data processing of the data processing section is completed, Orbit Acquisition of Orbit Information at the Time Defined by the Orbit Acquisition Definition Part A GNSS observation data acquisition unit for acquiring a GNSS data file at an observing station defined by the observing station defining unit; a determination unit for determining a setting condition for GNSS data processing defined by the processing condition defining unit, A coordinate change mapping unit for displaying coordinates or displacements of the observing stations defined as a result of the data processing, and a user interface unit for user operation of the data processing program operation of the system.

Further, in the GNSS automatic data processing system according to the present invention, the coordinate variation mapping unit displays the result of the coordinate / displacement of the defined GNSS station on the online map as a speed or a vector.

Meanwhile, a method for automatically processing GNSS data according to the present invention includes the steps of: acquiring predefined satellite orbit information data and GNSS observation data in a GNSS data automatic processing system and storing them in a database; Determining a setting condition for processing GNSS data in an automatic GNSS data processing system; Processing the satellite information data and the GNSS observation data defined in the automatic GNSS data processing system and copying the satellite information data and the GNSS observation data to the data processing unit; Processing the GNSS data processing program of the data processing unit to process the satellite orbital information data and the GNSS observation data to calculate information on the coordinate displacement of the observatory; Storing a processing result of a data processing unit in a database in an automatic GNSS data processing system; And the GNSS data automatic processing system are characterized by calculating the result of the coordinate / displacement of the defined GNSS station at the speed of the GNSS station.

Further, in the method for automatically processing GNSS data according to the present invention, the step of calculating the result of the coordinate / displacement of the defined GNSS station at the speed of the GNSS observing station may be such that the GNSS data automatic processing system selects the defined observing station and the period date Receiving an input; The automatic processing system of GNSS data processes the processing result and the period date of the inputted GNSS observing station, and removes an abnormal value among the inquired results; And the GNSS data automatic processing system includes a step of calculating a velocity for each coordinate component.

Further, in the method for automatically processing GNSS data according to the present invention, the GNSS data automatic processing system includes a step of calculating the coordinates of the processed GNSS station, the result of the displacement, at the speed of the GNSS station, .

Further, in the method of automatically processing GNSS data according to the present invention, the GNSS data automatic processing system includes a step of displaying the coordinates of the processed GNSS station and the result of the displacement as error ellipses.

According to the system and method for automatically processing GNSS data according to the present invention, GNSS data is automatically or manually acquired and managed through a DB, thereby enabling effective management of GNSS data.

In addition, by automating the process of the high-precision GNSS data processing program, the user can automatically process the GNSS observation data without any additional manipulation and obtain the coordinates and displacements. By cumulatively managing the coordinates and displacements of the processed GNSS stations, It is possible to calculate the coordinate variation (speed and direction).

On the other hand, expressing the speed and direction of GNSS station on an online electronic map such as Google Map, it is possible to intuitively understand perceptual fluctuation.

1 is a schematic diagram of an automatic GNSS data processing system according to the present invention.
2 is a detailed block diagram of a database of an automatic GNSS data processing system according to the present invention.
3 is a detailed block diagram of a control unit of the GNSS data automatic processing system according to the present invention.
4 is a data processing execution process of the GNSS data automatic processing system according to the present invention.
5 is a coordinate mapping process after data processing according to the automatic GNSS data processing method according to the present invention.

Hereinafter, embodiments of the present invention will be described in more detail with reference to the accompanying drawings.

1 is a schematic diagram of an automatic GNSS data processing system according to the present invention.

Referring to FIG. 1, the GNSS data automatic processing system according to the present invention includes three parts including a database 100, a data processing unit 200, and a control unit 500.

The database 100 includes data on satellite orbits, observation stations, and GNSS observations.

The data processing unit 200 includes a program for processing the precise GNSS data based on various data provided in the database 100.

The control unit 500 is for outputting speed and coordinates / displacement on an on-line map together with the orbit information and the GNSS observation data acquired from the database 100 and the data processed by the data processing unit 200 together.

2 is a detailed block diagram of a database of an automatic GNSS data processing system according to the present invention.

2, the database 100 includes a satellite orbit acquisition defining unit 110, a station defining unit 120, a processing condition defining unit 130, a satellite orbit information data unit 140, a GNSS observation data unit 150, and a data processing result unit 160.

The satellite trajectory acquisition definition unit 110 defines a precision force to use the precise orbital force (hereinafter referred to as precision power) of the GNSS satellite for data processing in order to calculate high-precision coordinates of the GNSS station.

IGS (International GNSS Service) and Bern University's CODE (Center for Orbit Determination in Europe) are representative institutions that provide precision power and provide them through Internet.

  The precision force files provided by IGS and CODE are as follows.

  - Orbit file: Includes precise orbital information of GNSS satellite

  - Earth Rotation Pole File: Includes information about the Earth's polar motion

  - Clock file: Includes error information of high precision atomic clock mounted on GNSS satellite.

  - Ion file: Including ionospheric activity information

  - DCB file: Includes hardware bias information for each GNSS satellite

  In addition, the precision force is provided in three types at regular intervals, and the current status and characteristics are as follows.

  - Ultra Rapid Orbit: Provides GNSS satellite precision data for the day, 12 hours ahead of that day, and 12 hours of satellite orbit prediction data from the date and time. This precision force is provided based on the day (D + 0).

  - Rapid Orbit: Rapid Orbit: usually provided one day after satellite orbit observation.

  - Final Orbit: The most precise final precision, usually provided 21 days after satellite orbital observation.

The table below summarizes the file statuses provided by CODE and IGS.

Type Orbit Earth Rotation Pole Clock ION DCB Ultra Rapid O O X O O Rapid O O O O O Final O O O O O

<CODE>

Type Orbit Earth Rotation Pole Clock ION DCB Ultra Rapid O O X X X Rapid O O O X X Final O O O X X

<IGS>

In the processing system of the present invention, the precision power provided by the IGS and the CODE is automatically downloaded four times a day at a time designated by the user, and the current status is inputted into the DB by precision date, type and type, do.

Meanwhile, the GNSS station definition unit 120 is for defining accurate information of the observation station to process the GNSS observation data through the data processing unit 200.

In order to process the GNSS observation data through the data processing unit 200, it is necessary to input accurate information of the observing station.

Further, in order to automatically acquire GNSS observation data via the Internet or a network, communication information (IP), connection ID, and password of an observing station providing GNSS observation data are required.

  The following information is essential for this course.

  - Observation site name: It should normally have a unique name that does not overlap with other observer sites.

  - GNSS antenna information: Type of antenna and type of dome covered by antenna

  - GNSS receiver information: the type of receiver used for data acquisition

  - Station coordinates / speed information: 3-D centering angle and velocity information

  - Reference Frame: Information about the reference frame

  - Reference Epoch: Base date for calculating coordinates

  - Data information: Access information for downloading data from Internet

Meanwhile, the processing condition defining unit 130 defines a setting condition for processing GNSS data in the data processing unit 200.

The information that must be input in order to process the GNSS data in the data processing unit 200 is as follows.

 - Data processing method: Select high-precision GNSS data processing program processing method (selection: double difference method, precision single positioning method)

  - Precision force type: Select the type of precision force to be used. (Optional: final precision force, rapid trajectory force, super rapid trajectory force)

  - Process Observation Site: Select the observation site that you want to process data on online map based on predefined observatory information and select it.

  - Select the reference coordinate system and reference time: Select the reference coordinate system and reference time at which the calculation result coordinates will be expressed.

  - Select the tide information file: Select the file that contains the tide information of the observatory to be processed.

The satellite orbit information data unit 140 stores the orbit information defined and acquired by the satellite orbit acquisition definition unit 110. [

In addition, the GNSS observation data part 150 stores GNSS observation data measured from an observing station defined in the observing station defining part 120.

The data processing result section 160 stores coordinates and displacements of the GNSS station where post-processing of the data of the satellite orbit information data and the GNSS observed data is completed by driving the data processing section 200.

3 is a detailed block diagram of a control unit of the GNSS data automatic processing system according to the present invention.

The orbit information acquiring unit 510 downloads the orbit information to the system through the Internet at a precision power providing agency at a time defined by the user based on the information input from the satellite orbit acquisition defining unit 110, The status of the file is stored in the satellite orbit information data unit 140 and is inquired by the high-precision GNSS data processing program of the data processing unit 200.

The GNSS observation data acquisition unit 520 sequentially downloads and acquires GNSS data files through an Internet connection at an observation station defined at a time defined by the user based on the information input from the station definition unit 120, The current state of the file is stored in the GNSS observation data unit 150 and is inquired by the high-precision GNSS data processing program of the data processing unit 200. [

The condition judgment execution unit 530 judges a setting condition for processing the GNSS data defined in the processing condition defining unit 130 and causes the data processing unit 200 to execute data processing when the condition is satisfied.

4 is a data processing execution process of the GNSS data automatic processing system according to the present invention.

If both the satellite orbit information data of the database and the GNSS observation data are acquired, the condition judgment execution unit 530 judges the setting condition for processing the GNSS data (S10).

The condition judgment execution unit 530 sequentially copies the satellite information data defined in the processing condition definition unit 130 and the GNSS observation data to the data processing unit 200 in turn (S20, S30).

At this time, the high-precision GNSS data processing program of the data processing unit 200 is operated to process the satellite information data and the GNSS observation data to calculate information on the coordinate displacement of the observatory (S40).

The processing result of the data processing unit 200 is stored in the data processing result unit 160 (S50).

The coordinate variation mapping unit 540 is used to display the coordinate / displacement of the processed observation station on an online map or the like.

In addition, the user interface unit 550 is a user operation portion for operating a program for data processing of the JNSES automatic data processing system according to the present invention.

Meanwhile, the result of the coordinate / displacement of the GNSS station defined by the coordinate variation mapping unit 540 is expressed through the calculation of the speed and the expression of the speed.

The coordinates of the GNSS station defined as a result of the data processing according to the present invention can be expressed by the concept of speed together with time.

5 is a coordinate mapping process after data processing according to the automatic GNSS data processing method according to the present invention.

First, the user selects and inputs an observation station defined in the system and a period date (S60).

The automatic GNSS data processing system inquires the processing result and the period date of the GNSS observing station input from the data processing result unit 160, and removes an abnormal value among the inquired results (S70, S80).

The automatic GNSS data processing system calculates the velocity for each coordinate component (S90).

The coordinates and the velocity of the observation point of the data calculated in the automatic GNSS data processing system are expressed in a chart (S100).

In addition, the GNSS data automatic processing system displays the coordinates and the speed of the observation point of the data calculated in the GNSS data automatic processing system on an online electronic map (S110).

In addition, the GNSS data automatic processing system expresses an observation point displacement of data calculated in the GNSS data automatic processing system by an error ellipse (S120).

The moving speed of the GNSS station can be calculated using the coordinates or the displacement amount of each station stored in the data processing result unit 160. The calculation process at this time is as follows.

 First, the processing result of the GNSS observation station and the inquiry period are inquired by the data processing result section 160.

Next, an outlier among the retrieved results is removed. The elimination of an ideal value uses the average and standard deviation calculated from the coordinates retrieved, and removes coordinate values outside the 1-sigma, 2-sigma, or 3-sigma region centered on the average.

  Next, a straight line equation representing the most suitable speed for the period searched using the least squares method is obtained as a result of the process in which the ideal value is removed.

  Finally, the slope is calculated from the derived straight line equation to calculate the moving speed of the corresponding station.

On the other hand, the present invention relates to a method for indicating the speed of the corresponding observation station.

The speed calculated in the above procedure can be expressed by the following two methods.

The first is a method of displaying the calculated speed together with the coordinates or the displacement of each observation station retrieved by the processing result section 160. The second is a method of displaying the size and direction of movement speed of each observation station in the form of a vector on an online electronic map Method.

  The first method expresses the coordinates or displacements of each station by date in graph form and expresses the calculated straight line (speed) on it. This method is suitable for confirming the magnitude of the absolute velocity of crustal variation.

  In the second method, the velocity vector is divided into the horizontal component and the vertical component on the on-line electronic map.

This method is suitable for comparing the direction of the speed and the relative magnitude of the speed between the various stations. Since it is a relative magnitude comparison rather than absolute magnitude magnitude comparison, magnification can be applied to enlarge or reduce the velocity vector.

Then, the coordinate change of the observatory can express the distribution of the displacement with the error ellipse.

 Error ellipses are the best way to represent the spatial distribution of the coordinates of the processed observations. An error ellipse is drawn by calculating the elements (length / short axis length and rotation angle) that can draw an error ellipse by using the coordinates of each observation station stored in the processing result section 160. The calculation process is as follows.

A. The user selects the display period of the GNSS station and the error ellipse (start date, end date).

B. The two-dimensional longitudinal latitude coordinates or displacements corresponding to the selected station and period are inquired in the data processing result section (160).

 C. Calculate the space center, rotation angle, major axis, and minor axis components of the error ellipse using coordinates (latitude (dN), hardness (dE)) or displacement among the queried elements.

D. The error ellipse is displayed on the online electronic map using the calculated element or it is created as the image using the plotting program.

E. Express the image generated by day or time by animation.

 While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is clearly understood that the same is by way of illustration and example only and is not to be construed as limiting the scope of the invention as defined by the appended claims. .

100: Database 110: Satellite Orbit Acquisition Definition section
120: Observation site defining section 130: Processing condition defining section
140: satellite orbital information data part 150: GNSS observation data part
160: Data processing result section 200: Data processing section
500: control unit 510: trajectory information acquisition unit
520: GNSS observation data acquisition unit 530: condition judgment execution unit
540: Coordinate variation mapping unit 550: User interface unit

Claims (6)

delete delete In an automatic GNSS data processing method,
Acquiring predefined satellite orbit information data and GNSS observation data in a GNSS data automatic processing system and storing the satellite orbit information data and the GNSS observation data in a database;
Determining a setting condition for processing GNSS data in an automatic GNSS data processing system;
Searching the satellite orbit information data and GNSS observation data defined in the automatic GNSS data processing system and copying them to a data processing unit;
Processing the GNSS data processing program of the data processing unit to process the satellite orbital information data and the GNSS observation data to calculate information on the coordinate / displacement of the observatory;
Storing a processing result of a data processing unit in a database in an automatic GNSS data processing system; And
The automatic processing system of GNSS data includes calculating the result of the coordinate / displacement of the defined GNSS station at the speed of the GNSS station,
Calculating the result of the coordinate / displacement of the defined GNSS station at the speed of the GNSS station,
The automatic GNSS data processing system includes a step of selecting a defined observation station and a period date,
The automatic processing system of GNSS data processes the processing result and the period date of the inputted GNSS observing station, and removes an abnormal value among the inquired results; And
Wherein the GNSS data automatic processing system includes a step of calculating a velocity for each coordinate component. delete The method of claim 3,
Wherein the GNSS data automatic processing system includes a step of calculating the coordinates of the processed GNSS station, the result of the displacement, at the speed of the GNSS station, and displaying it on an on-line electronic map. The method of claim 3,
Wherein the GNSS data automatic processing system includes a step of displaying the coordinates of the processed GNSS station and the result of the displacement as an error ellipse of the coordinate results.

Images