CN111694021B - GNSS environment model-based single-station landslide deformation monitoring and early warning method - Google Patents
GNSS environment model-based single-station landslide deformation monitoring and early warning method Download PDFInfo
- Publication number
- CN111694021B CN111694021B CN202010460631.7A CN202010460631A CN111694021B CN 111694021 B CN111694021 B CN 111694021B CN 202010460631 A CN202010460631 A CN 202010460631A CN 111694021 B CN111694021 B CN 111694021B
- Authority
- CN
- China
- Prior art keywords
- gnss
- early warning
- monitoring
- environment
- observation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/33—Multimode operation in different systems which transmit time stamped messages, e.g. GPS/GLONASS
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
- G08B21/02—Alarms for ensuring the safety of persons
- G08B21/10—Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L67/00—Network arrangements or protocols for supporting network services or applications
- H04L67/01—Protocols
- H04L67/02—Protocols based on web technology, e.g. hypertext transfer protocol [HTTP]
- H04L67/025—Protocols based on web technology, e.g. hypertext transfer protocol [HTTP] for remote control or remote monitoring of applications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L67/00—Network arrangements or protocols for supporting network services or applications
- H04L67/01—Protocols
- H04L67/12—Protocols specially adapted for proprietary or special-purpose networking environments, e.g. medical networks, sensor networks, networks in vehicles or remote metering networks
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Environmental & Geological Engineering (AREA)
- Emergency Management (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The invention discloses a GNSS environment model-based single-station landslide deformation monitoring and early warning method, which comprises the following steps: deploying terminal equipment, acquiring data, determining the initial position of a monitoring point, calculating environment modeling source data, modeling a sunday environment, solving a correlation coefficient and applying early warning. The monitoring and early warning method mainly utilizes the low-cost consumer GNSS positioning terminal to acquire data and perform environment modeling, early identifies the landslide acceleration state through the change trend of the correlation coefficient of the environment model and effectively early warns, so that the monitoring cost is greatly reduced, and a new thought and a new method are provided for the GNSS landslide monitoring and early warning.
Description
Technical Field
The invention belongs to the technical field of monitoring and early warning, and particularly relates to a method for monitoring and early warning of landslide deformation of a single station based on a GNSS environment model.
Background
China is one of the most serious countries in the world with geological landslide disasters. According to the incomplete statistics of the national geological disaster network, the direct economic loss is 27.7 billion yuan only from the occurrence of geological disasters 6181 in 2019 nationwide; the geological disaster 948 can be successfully forecasted, 2 million persons can be possibly injured, the direct economic loss is avoided by 8.3 hundred million yuan, and the disaster prevention and reduction effects are remarkable. And the landslide disaster loss accounts for 61.2 percent of the total disaster, and the implementation of efficient landslide monitoring and early warning has important practical significance for the disaster prevention and reduction work in China.
The GNSS technology has high-precision, all-time and all-weather monitoring capability in landslide monitoring and early warning, and is an important means for directly acquiring surface three-dimensional vector deformation at present. With the modernization and the gradual improvement of GNSS, currently, there are up to more than 110 global in-orbit GNSS satellites, which become an important means for global high-precision positioning and navigation services.
The GNSS positioning technology has been widely applied to various deformation monitoring fields, in particular to landslide disaster deformation monitoring, by virtue of its unique technical advantages, and mainly uses the RTK positioning technology. The cost of a single GNSS navigation positioning terminal which can be used for high-precision landslide monitoring is generally higher (more than 3 thousands), the cost of a common measurement type double-frequency receiving terminal is at least more than ten thousand yuan, the cost of a GNSS application terminal used for consumer-level positioning is less than one hundred yuan at present, and the positioning precision is low (more than 10m), so that the landslide monitoring and early warning requirements cannot be met.
Therefore, how to utilize a single GNSS positioning terminal with low cost to realize effective early warning of landslide deformation monitoring has important practical significance for geological disaster prevention and reduction work in China.
Disclosure of Invention
Aiming at the problems of high early warning cost, low success rate and the like of the existing GNSS landslide deformation monitoring technology, the invention provides a single-station landslide deformation monitoring and early warning method based on a GNSS environment model.
The technical scheme of the invention is as follows: a GNSS environment model-based single-station landslide deformation monitoring and early warning method comprises the following steps:
s1: deploying terminal equipment and performing data acquisition
The method comprises the steps that a GNSS application terminal is deployed on a monitoring point (landslide deformation characteristic point), and after the GNSS application terminal is started, original observation data of the monitoring point are collected through a 4G or 5G wireless communication network and the like;
s2: determining the initial position of a monitoring point
After the original observation data of several hours (>1 hour) are collected, downloading a real-time broadcast ephemeris file through an Internet IGS service center, and roughly estimating the approximate position of a monitoring point by using a pseudo-range single-point positioning technology;
s3: computing environment modeling source data
Based on the broadcast ephemeris, performing integrity check on the original observation data acquired by S1, calculating satellite azimuth angle and altitude angle information of each epoch through a receiver approximate coordinate, summarizing satellite sight line information, and outputting continuous observation arc-segment observation values, wherein the initial and final observation values of the arc-segment observation values are marked as head and tail points;
s4: sunday environment modeling
Performing gross error elimination on the head and tail point observed values output by the S3, constructing an observation equation by using a plurality of head and tail point observed values, selecting a least square criterion to solve a fitting coefficient in the equation by a fitting system, and establishing an environment model according to the fitting coefficient;
s5: correlation coefficient solving
Setting intervals by taking days as units, and jointly solving correlation coefficients in the environment models according to the corresponding environment models every day;
s6: early warning application
And setting a proper threshold value in practical application, comparing and judging the threshold value and the correlation coefficient calculated in the step S5, if the correlation coefficient is lower than the threshold value, carrying out early warning, and if not, continuing monitoring.
Further, the raw observation data in S1 includes raw pseudorange and ephemeris data.
Further, the GNSS application terminal in S1 is a consumer-grade single frequency GNSS.
Further, the integrity detection method of the original observation data in S3 includes: and comparing the raw observation data collected in the step S1 with the real-time broadcast ephemeris downloaded in the step S2 every 24 hours to determine whether the raw observation data correspond to the real-time broadcast ephemeris or not, so as to detect the integrity of the raw observation data.
Further, the method for constructing the environment model in S4 includes the following steps:
s41: and (2) constructing an observation equation for a plurality of head and tail point observation values by using a formula (1):
in the formula (1), n represents the fitting order, a i 、b i 、c i Denotes the fitting coefficient, e x Representing the measurement height angle, x representing e x A corresponding azimuth angle; n is 360/0.5, i is {0,0.5, …,360 };
s42: selecting a least square criterion to carry out fitting system solution on the observation equation to obtain a fitting coefficient, and establishing an environment model according to the fitting coefficient by taking 0.5 degrees as a step length from 0 to 360 degrees, as shown in a formula (2):
in equation (2), a 'model coefficient vector, e' i Representing the virtual satellite altitude using the model coefficients back-deducted, i ═ 0,0.5, …, 360.
Further, the specific method for solving the correlation coefficient in S5 is as follows:
to set an environmental model a 'of day 1' 0 Reference, environmental model a 'for day k every 24 hours' k And a' 0 Is solved as shown in equation (3):
wherein, Cov (a' k ,a′ 0 ) Is a' k And a' 0 Covariance, Var [ a' k ]Is a' k Variance of (1), Var [ a' 0 ]Is a' 0 The variance of (c).
The invention has the beneficial effects that: the monitoring and early warning method mainly utilizes the low-cost consumer GNSS positioning terminal to acquire data and perform environment modeling, early identifies the landslide acceleration state through the change trend of the correlation coefficient of the environment model and effectively early warns, so that the monitoring cost is greatly reduced, and a new thought and a new method are provided for the GNSS landslide monitoring and early warning.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a view showing modeling results of GNSS environment at monitoring points for 19 consecutive weeks according to embodiment 2 of the present invention;
fig. 3 is a timing chart of three-dimensional position change and correlation coefficient of the monitoring points in embodiment 2 of the present invention.
Detailed Description
Example 1: in this embodiment, a single-station landslide deformation monitoring and early warning method based on a GNSS environment model is designed, and is mainly used for providing a new monitoring and early warning method based on a consumer-grade single-frequency GNSS, which is specifically described based on the method flowchart in fig. 1.
S1: deploying terminal equipment and performing data acquisition
The method comprises the steps that a GNSS application terminal is deployed on a monitoring point (landslide deformation characteristic point), and after the GNSS application terminal is started, original observation data of the monitoring point are collected through a 4G or 5G wireless communication network and the like, wherein the original observation data comprise original pseudo-range and ephemeris data; the GNSS application terminal can be a consumer-grade single-frequency GNSS, or can be a professional type with centimeter-level precision or a common measurement type with sub-meter-level precision, and the GNSS comprises one or more of BDS/GPS/GLONASS/GALILEO.
S2: determining the initial position of a monitoring point
After the original observation data of several hours (>1 hour) are collected, downloading a real-time broadcast ephemeris file through an Internet IGS service center, and roughly estimating the approximate position of a monitoring point by using a pseudo-range single-point positioning technology;
s3: computing environment modeling source data
On the basis of a broadcast ephemeris, comparing original observation data acquired in S1 with the real-time broadcast ephemeris downloaded in S2 every 24 hours to determine whether the original observation data are corresponding or not, performing integrity check on the original observation data acquired in S1, calculating satellite azimuth angle and altitude angle information of each epoch through a receiver approximate coordinate, summarizing satellite sight line information, and outputting continuous observation arc section observation values, wherein the initial observation value and the final observation value of the arc section observation values are marked as head and tail points;
s4: sunday environment modeling
Performing gross error elimination on the head and tail point observed values output by the S3, constructing an observation equation by using a plurality of head and tail point observed values, selecting a least square criterion to solve a fitting coefficient in the equation by a fitting system, and establishing an environment model according to the fitting coefficient;
the construction method of the environment model comprises the following steps:
s41: the formula (1) is utilized to construct an observation equation for a plurality of head and tail point observation values:
in the formula (1), n represents the fitting order, a i 、b i 、c i Denotes the fitting coefficient, e x Representing the measurement height angle, x representing e x A corresponding azimuth angle; n is 360/0.5, i is {0,0.5, …,360 };
s42: selecting a least square criterion to carry out fitting system solution on the observation equation to obtain a fitting coefficient, and establishing an environment model according to the fitting coefficient by taking 0.5 degrees as a step length from 0 to 360 degrees, as shown in a formula (2):
in equation (2), a 'model coefficient vector, e' i Representing the virtual satellite altitude using the model coefficients back-deducted, i ═ 0,0.5, …, 360.
S5: correlation coefficient solving
Setting intervals by taking days as units, and jointly solving correlation coefficients in the environment models according to the corresponding environment models every day;
the concrete method for solving the correlation coefficient comprises the following steps:
to set an environmental model a 'of day 1' 0 Reference, environmental model a 'for day k every 24 hours' k And a' 0 Is solved as shown in equation (3):
wherein, Cov (a' k ,a′ 0 ) Is a' k And a' 0 Covariance, Var [ a' k ]Is a' k Variance of (1), Var [ a' 0 ]Is a' 0 The variance of (c).
S6: early warning application
And setting a proper threshold value in practical application, comparing and judging the threshold value and the correlation coefficient calculated in the step S5, if the correlation coefficient is lower than the threshold value, carrying out early warning, and if not, continuing monitoring.
Example 2: in this embodiment, the method in embodiment 1 is used for actual measurement, data of a certain landslide monitoring point in a creep stage is selected for an experiment, the monitoring point gradually sinks along with the time, original observation data are collected and environment modeling is performed, the result is shown in fig. 2, the direction indicated by an arrow in fig. 2 represents the modeling result of 19 consecutive weeks from the 1 st week to the 19 th week, and it can be seen that the satellite height angle of the environment model in the same azimuth gradually increases along with the creep movement of the monitoring point.
In the time chart of the three-dimensional position change and the correlation coefficient of the monitoring point given in conjunction with fig. 3, the horizontal axis in fig. 3 represents the time axis, the left vertical axis represents the 3D vector displacement, and the right vertical axis represents the correlation coefficient. As can be seen from the displacement sequence, the monitoring point passes through a creep stage- > an acceleration stage- > and recovers the creep; the correlation coefficient is correspondingly restored to be stable from a stable stage- > an acceleration stage- >; therefore, it can be seen that the trend of the correlation coefficient can be consistent with the position change trend of the monitoring point, and therefore effective landslide early warning can be carried out through setting a reasonable threshold value.
Claims (6)
1. A GNSS environment model-based single-station landslide deformation monitoring and early warning method is characterized by comprising the following steps:
s1: deploying terminal equipment and performing data acquisition
The GNSS application terminal is deployed on a monitoring point, and the original observation data of the monitoring point is collected through a wireless communication network after the GNSS application terminal is started;
s2: determining the initial position of a monitoring point
After the original observation data are collected, downloading a real-time broadcast ephemeris file, and roughly estimating the rough position of the monitoring point by utilizing a pseudo-range single-point positioning technology;
s3: computing environment modeling source data
Based on the broadcast ephemeris, performing integrity check on the original observation data acquired by S1, calculating satellite azimuth angle and altitude angle information of each epoch through a receiver approximate coordinate, summarizing satellite sight line information, and outputting continuous observation arc-segment observation values, wherein the initial and final observation values of the arc-segment observation values are marked as head and tail points;
s4: sunday environment modeling
Performing gross error elimination on the head and tail point observed values output by the S3, constructing an observation equation by using a plurality of head and tail point observed values, selecting a least square criterion to solve a fitting coefficient in the equation by a fitting system, and establishing an environment model according to the fitting coefficient;
s5: correlation coefficient solving
Setting intervals by taking days as units, and jointly solving correlation coefficients in the environment models according to the corresponding environment models every day;
s6: early warning application
And setting a proper threshold value in practical application, comparing and judging the threshold value and the correlation coefficient calculated in the step S5, if the correlation coefficient is lower than the threshold value, carrying out early warning, and if not, continuing monitoring.
2. The GNSS environment model-based single-station landslide deformation monitoring and early warning method according to claim 1, wherein the raw observation data in S1 comprises raw pseudoranges and ephemeris data.
3. The GNSS environment model-based single-station landslide deformation monitoring and early warning method according to claim 1, wherein the GNSS application terminal in S1 is a consumer-grade single-frequency GNSS.
4. The GNSS environment model-based single-station landslide deformation monitoring and early warning method of claim 1, wherein the integrity detection method of the original observation data in S3 comprises: and comparing the raw observation data collected in the step S1 with the real-time broadcast ephemeris downloaded in the step S2 every 24 hours to determine whether the raw observation data correspond to the real-time broadcast ephemeris or not, so as to detect the integrity of the raw observation data.
5. The GNSS environment model-based single-station landslide deformation monitoring and early warning method according to claim 1, wherein the construction method of the environment model in S4 comprises the following steps:
s41: and (2) constructing an observation equation for a plurality of head and tail point observation values by using a formula (1):
in the formula (1), n represents the fitting order, a i 、b i 、c i Denotes the fitting coefficient, e x Representing the measurement height angle, x representing e x A corresponding azimuth angle; n is 360/0.5, i is {0,0.5, …,360 };
s42: selecting a least square criterion to carry out fitting system solution on the observation equation to obtain a fitting coefficient, and establishing an environment model according to the fitting coefficient by taking 0.5 degrees as a step length from 0 to 360 degrees, as shown in a formula (2):
6. The GNSS environment model-based single-station landslide deformation monitoring and early warning method as claimed in claim 1, wherein the specific method for solving the correlation coefficient in S5 is as follows:
to set an environmental model a 'of day 1' 0 Reference, environmental model a 'for day k every 24 hours' k And a' 0 Is solved as shown in equation (3):
wherein, Cov (a' k ,a′ 0 ) Is a' k And a' 0 Covariance, Var [ a' k ]Is a' k Variance of (1), Var [ a' 0 ]Is a' 0 The variance of (c).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010460631.7A CN111694021B (en) | 2020-05-26 | 2020-05-26 | GNSS environment model-based single-station landslide deformation monitoring and early warning method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010460631.7A CN111694021B (en) | 2020-05-26 | 2020-05-26 | GNSS environment model-based single-station landslide deformation monitoring and early warning method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111694021A CN111694021A (en) | 2020-09-22 |
CN111694021B true CN111694021B (en) | 2022-08-05 |
Family
ID=72478602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010460631.7A Active CN111694021B (en) | 2020-05-26 | 2020-05-26 | GNSS environment model-based single-station landslide deformation monitoring and early warning method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111694021B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114089055B (en) * | 2021-09-30 | 2024-08-09 | 安徽继远软件有限公司 | Method and system for monitoring safety state of power grid limited space operation personnel |
CN115659497B (en) * | 2022-10-18 | 2023-07-21 | 武汉理工大学 | Ship body modeling method based on direct control free deformation |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8457882B2 (en) * | 2011-04-19 | 2013-06-04 | Honeywell International Inc. | Systems and methods for navigation in a GPS-denied environment |
CN104597471B (en) * | 2015-01-13 | 2017-04-19 | 华东师范大学 | Orientation attitude determination method oriented to clock synchronization multi-antenna GNSS receiver |
CN106405589A (en) * | 2016-06-24 | 2017-02-15 | 西安科技大学 | Method and device for determining global ionized layer grid model |
CN106500928B (en) * | 2016-11-02 | 2018-09-04 | 济南大学 | A kind of reservoir dam Seepage Safety online auto monitoring system |
CN108332649B (en) * | 2018-02-07 | 2020-04-24 | 桂林电子科技大学 | Landslide deformation comprehensive early warning method and system |
CN109977181A (en) * | 2019-02-28 | 2019-07-05 | 广东赛诺科技股份有限公司 | A method of based on bridge stereoscopic monitoring |
CN111025357B (en) * | 2019-12-31 | 2021-11-19 | 中国科学院国家授时中心 | High-precision GNSS virtual SNR weighted positioning method |
-
2020
- 2020-05-26 CN CN202010460631.7A patent/CN111694021B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN111694021A (en) | 2020-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP4206743A1 (en) | Method for vehicle positioning, related apparatus, device, and storage medium | |
CN110515097B (en) | GNSS satellite observation gross error rejection method and device applied to reference station | |
CN107765275B (en) | Wide-area differential positioning method, device, terminal and computer readable storage medium | |
CN112327340B (en) | Terminal positioning accuracy evaluation method, device, equipment and medium | |
CN108490459B (en) | Method and system for applying precision and risk balance to GNSS position service | |
CN1670542A (en) | Method and system for using altitude information in a satellite positioning system | |
CN104483691B (en) | A kind of GNSS combines accurate one-point positioning method | |
CN111694021B (en) | GNSS environment model-based single-station landslide deformation monitoring and early warning method | |
CN103900580A (en) | Compass/GPS (global positioning system) and INS (inertial navigation system) combination vehicle navigation positioning system based on GIS (geographic information system) technology | |
CN103364803A (en) | Satellite selection method and satellite navigation positioning method applying the satellite selection method | |
CN111290004A (en) | Pseudo-range differential positioning method, pseudo-range differential positioning device, electronic equipment and storage medium | |
CN117739797B (en) | Beidou/GNSS-based multi-time-scale deformation monitoring method | |
CN107807373A (en) | GNSS high-precision locating methods based on mobile intelligent terminal | |
CN117826200A (en) | PPP-B2B-based marine real-time precise positioning method, system and medium | |
CN115902968A (en) | PPP terminal positioning method based on Beidou third GEO broadcast enhancement information | |
KR101617033B1 (en) | Global standard point positioning apparatus using multi global satellite positioning systems and the method thereof | |
CN108491616B (en) | Ionized layer vertical total electron content modeling method based on ellipsoid harmonic function theory | |
CN112230254B (en) | Correction method and device for GPS carrier phase multipath error | |
Tarig | Positioning with wide-area GNSS networks: Concept and application | |
Raghuvanshi et al. | Precise positioning of smartphones using a robust adaptive Kalman filter | |
CN114355420B (en) | PPP product positioning method and device for distributed Beidou position service center | |
Blomenhofer et al. | Sensitivity analysis of the Galileo integrity performance dependent on the ground sensor station network | |
CN113866801A (en) | Beidou satellite positioning accuracy evaluation improvement method and system based on vertical projection | |
CN114167469A (en) | Vehicle navigation information monitoring method and device based on 5G/GNSS combination | |
KR100264977B1 (en) | Method for the real-time estimation of user position in dgps |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |