CN107991676B - Troposphere error correction method of satellite-borne single-navigation-pass InSAR system - Google Patents
Troposphere error correction method of satellite-borne single-navigation-pass InSAR system Download PDFInfo
- Publication number
- CN107991676B CN107991676B CN201711250005.XA CN201711250005A CN107991676B CN 107991676 B CN107991676 B CN 107991676B CN 201711250005 A CN201711250005 A CN 201711250005A CN 107991676 B CN107991676 B CN 107991676B
- Authority
- CN
- China
- Prior art keywords
- error
- troposphere
- coordinate system
- dimensional
- scene
- 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
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
Abstract
The invention provides a troposphere error correction method of a satellite-borne single-navigation InSAR system. The technical scheme is as follows: firstly, calculating a two-dimensional measurement error introduced by a troposphere in a distance-height plane by utilizing a troposphere influence effect model of a satellite-borne single-navigation InSAR system, an auxiliary radar incident angle and a priori troposphere zenith delay; then, decomposing the two-dimensional error into a three-dimensional ground-fixed coordinate system for representation through projection transformation according to the spaceborne SAR imaging geometric model; and finally, correcting the troposphere error influence in the detected scene point by point. The method has high correction precision, can be used for directly processing the interference positioning result, and has simple and efficient operation process.
Description
Technical field
The invention belongs to the interleaving techniques field of space remote sensing and atmospheric effect, navigated InSAR for spaceborne list
The tropospheric error correction side of (Interferometric Synthetic Aperture Radar, interference synthetic aperture radar)
Method.
Background technique
Spaceborne list InSAR system of navigating is single using a pair of the same area obtained under main and auxiliary radar difference observation visual angle
The three-dimensional information of earth's surface is obtained depending on the interference phase difference inverting between complex pattern, complete the mapping task to global landform and is generated
High-precision DEM (Digital Elevation Model, digital elevation model).
Atmospheric effect is to restrict spaceborne list to navigate a key factor of InSAR system intervention measurement accuracy.Wherein, convection current
Layer is the bottom of atmosphere, connects with earth's surface downwards, extends upwardly at about earth's surface overhead 16km, constituent contains a variety of
Gas (nitrogen, oxygen, carbon dioxide etc.) and steam, numerous weather phenomena occur in troposphere.The orbit altitude of satellite
It is above convection current layer height, will receive tropospheric influence when the electromagnetic wave of radar emission passes through troposphere, gives satellite-borne SAR figure
As introducing certain time delay and phase error, and then influence the performance of InSAR system.Studies have shown that conventional convection layer is to spaceborne
The influence for InSAR system of singly navigating up to several meters of magnitudes, must usually be corrected in high-precision interferometry.Therefore have
Necessity apply tropospheric radio propagation mechanism, propose for spaceborne list navigated InSAR system high-precision tropospheric error correction
Method.Currently, it is not yet found that closing the related data of bearing calibration.
Summary of the invention
The present invention in order to effectively solve the problems, such as spaceborne list navigated InSAR system tropospheric error influence, propose a kind of base
It navigated InSAR system tropospheric error bearing calibration in the spaceborne list of priori tropospheric zenith delay.This method correction accuracy is high,
Treatment process is relatively easy, navigates InSAR system suitable for the spaceborne list under known priori tropospheric zenith delay distribution situation
Tropospheric error correction.
Basic ideas of the invention are: firstly, navigated InSAR system troposphere influential effect model using spaceborne list, benefit
With auxiliary radar incidence angle and priori tropospheric zenith delay, calculates the two-dimensional measurement that troposphere introduces in distance-height plane and miss
Difference;Then, two-dimentional error is decomposed under three-dimensional body-fixed coordinate system by projective transformation according to Space-borne SAR Imaging geometrical model
It indicates;Finally, pointwise correction, which is tested the tropospheric error in scene, to be influenced.
The technical scheme is that
It is known at any time under the orbit coordinate A (x that navigated under the auxiliary star body-fixed coordinate system of InSAR system of spaceborne listA, yA,
zA), auxiliary radar incidence angle θ is tested the corresponding tropospheric zenith time delay τ of scenezenith;It is known to be influenced by tropospheric
It is tested coordinate B (x of any point under body-fixed coordinate system in scene afterwardsB, yB, zB), which is fixed by directly reading InSAR
Position result obtains.
It is B (x for coordinate in tested sceneB, yB, zB) arbitrary point, using following steps complete scene three-dimensional coordinate
Tropospheric error influences correction:
Step 1: navigating InSAR system troposphere influential effect model according to spaceborne list, distance-height is calculated
Horizontal measurement error and vertical measurement error in plane.
According to auxiliary radar incidence angle θ tropospheric zenith delay τ corresponding with tested scenezenithTroposphere is calculated to spaceborne
The horizontal measurement error that InSAR system of singly navigating introducesWith vertical measurement error
In above formula, c indicates the light velocity.
Step 2: two-dimensional measurement error pro is transformed to be indicated under three-dimensional body-fixed coordinate system.
Firstly, calculating vertical measurement error using following formulaError vector under body-fixed coordinate system
Wherein xBC、yBCAnd zBCIt respectively indicates are as follows:
Then, horizontal measurement error is calculatedError vector under under body-fixed coordinate system
In above formula, O is that coordinate of the earth's core under body-fixed coordinate system is (0,0,0), because A, B coordinate are it is known that vectorWithIt can ask.
Finally, calculating the body-fixed coordinate system three-dimensional coordinate offset vector that troposphere introduces according to the following formula
Step 3: being tested the three-dimensional body-fixed coordinate system coordinate of scene after obtaining correction tropospheric error influence.
The body-fixed coordinate system three-dimensional coordinate offset vector introduced according to troposphereIt is B (x to coordinateB, yB, zB) appoint
Meaning point correction tropospheric error influences, and the three-dimensional body-fixed coordinate system B'(x of measured point in scene is tested after being correctedB', yB', zB'),
Wherein xB'、yB'And zB'It respectively indicates are as follows:
xB'=xB+xBE
yB'=yB+yBE
zB'=zB+zBE
Following technical effect can use to obtain using the present invention:
It navigated InSAR system tropospheric error bearing calibration the present invention provides a kind of spaceborne list, this method is based on priori
Tropospheric zenith delay distribution, navigated InSAR system troposphere influential effect model using spaceborne list, by will theoretical two dimension
Measurement error is transformed under three-dimensional body-fixed coordinate system the troposphere for indicating the tested scene of InSAR system, it can be achieved that spaceborne list navigated
Error accurately corrects point by point.The bearing calibration fully considered spaceborne list navigated InSAR basic principle and tropospheric propagation effect
The space-variant answered, correction accuracy is high, and can be used for directly handling interference positioning result, and operating process is simple and efficient.
Detailed description of the invention
Fig. 1 is the principle of the present invention flow diagram;
Fig. 2 is the radar system set in emulation and troposphere basic parameter;
Fig. 3 is 1 spatial distribution map of priori tropospheric zenith delay in emulation experiment;
Fig. 4 is the three dimensional effects error space distribution map that troposphere introduces InSAR positioning result in emulation experiment;
Fig. 5 be the embodiment of the present invention in spaceborne list navigated InSAR system tropospheric error correct residual error experimental result picture.
Specific embodiment
It elaborates in the following with reference to the drawings and specific embodiments to the present invention.Wherein, the spaceborne list in embodiment navigated
The method that InSAR system and tested scene are all made of emulation generates.Specific principle of simulation and process refer to Min Wang,
Diannong Liang, which is equal to 2007, is published in IEEE International Geoscience and Remote
Article SBRAS-An Advanced Simulator of Spaceborne on Sensing Symposium.
Fig. 1 is the principle of the present invention flow diagram, and whole flow process is divided into three big steps.It is flat to calculate distance-height for the first step
Troposphere horizontal measurement error and vertical measurement error in face;Second step, using satellite-borne SAR geometrical relationship, by two-dimensional measurement
Error pro transforms to be indicated under three-dimensional body-fixed coordinate system, and it is inclined to obtain the three-dimensional body-fixed coordinate system that troposphere introduces tested scene
It moves;Third step, point-by-point accurate correction are tested the measurement error that measured point is introduced by tropospheric zenith delay in scene, are corrected
Tropospheric error is tested the three-dimensional body-fixed coordinate system coordinate of measured point in scene after influencing.
Fig. 2 is the radar system set in emulation and troposphere basic parameter.Spaceborne list navigated the work of InSAR system in X
Wave band (radar signal centre frequency is 9.65GHz), it is residual in order to compare the correction of priori tropospheric zenith delay of different accuracy
Difference, the priori tropospheric zenith delay that two kinds of precision are respectively adopted in the embodiment of the present invention carry out tropospheric error correction experiment.
Wherein, priori tropospheric zenith delay 1 is at a time to be tested the corresponding true tropospheric zenith delay of scene, in scene
For space-variant, zenith delay variation range is 2.27 meters to 2.54 meters;Priori tropospheric zenith delay 2 is experience tropospheric zenith
Delay, takes 2.3 meters of constant value.
Fig. 3 is 1 spatial distribution map of priori tropospheric zenith delay in emulation experiment.The sky of priori tropospheric zenith delay 1
Between coordinate grid it is identical as tested scene coordinate grid, it is 1117 that horizontal axis, which represents distance to, pixel number, in figure, and the longitudinal axis represents orientation
It is 1446 to, pixel number, distance is to being 3km with orientation size.
Gray value indicates the tropospheric zenith delay value in tested scene in figure, and gray scale gets over superficial and shows tropospheric zenith in figure
Length of delay is bigger.From figure 3, it can be seen that the troposphere distribution in tested scene is distributed close to true space troposphere.
Fig. 4 is the three dimensional effects error space distribution map that introduces to InSAR positioning result of troposphere in emulation experiment, in figure
It is 1117 that horizontal axis, which represents distance to, pixel number, and the longitudinal axis represents orientation, and pixel number is 1446, and distance is to equal with orientation size
For 3km.Figure (a), which represents X-direction troposphere under body-fixed coordinate system, influences error, and whole scene mean value is about -2.17 meters;Scheme (b)
Representing Y direction troposphere under body-fixed coordinate system influences error, and whole scene mean value is about -2.8 meters;Figure (c) represents body-fixed coordinate system
It is that lower Z-direction troposphere influences error, whole scene mean value is about 0.51 meter.Gray value indicates that troposphere influences error in figure
Size, Fig. 4 (a) and gray scale in (b) get over superficial and show that error influence is smaller, and gray scale, which is more deeply felt, in Fig. 4 (c) shows that error influence is smaller.
As seen from the figure, troposphere spaceborne list navigated error caused by InSAR system influence it is very significant, usually up to several meters of magnitudes,
And more significant space-variant is presented.
Fig. 5 be the embodiment of the present invention in spaceborne list navigated InSAR system tropospheric error correct residual error experimental result picture, figure
It is 1117 that middle horizontal axis, which represents distance to, pixel number, and the longitudinal axis represents orientation, and pixel number is 1446, distance to orientation size
It is 3km.Gray value indicates the residual values that tropospheric error correction in scene is tested after being corrected using the present invention, figure in figure
5 (a), (c), (e) and gray scale in (f), which are more deeply felt, shows that correction residual error is smaller, and Fig. 5 (b) gets over superficial with gray scale in (d) and shows correction residual error
It is smaller.Wherein, Fig. 5 (a), (c) and (e) are respectively to use the three-dimensional correction that priori tropospheric zenith delay 1 obtains in Fig. 2 residual
Difference, the whole scene mean value of the tropospheric correction residual error of body-fixed coordinate system X-axis, Y-axis and Z-direction respectively may be about 4.4,2.6 and
14.4mm, correction accuracy respectively may be about 0.20%, 0.09% and 2.82%, illustrate to work as the corresponding accurate convection current of known tested scene
The correction residual error that layer zenith delay carries out tropospheric error correction is smaller, correction accuracy is higher.Fig. 5 (b), (d) and (f) are respectively
The three-dimensional correction residual error obtained using priori tropospheric zenith delay 2 in Fig. 2, pair of body-fixed coordinate system X-axis, Y-axis and Z-direction
The whole scene mean value of tropospheric correction residual error respectively may be about -85.2, -113.3 and 35.4mm, correction accuracy respectively may be about 3.93%,
4.04% and 6.94%, correction residual result has the characteristic of space-variant.Wherein, the calculation formula of correction accuracy is that correction residual error is removed
Size is influenced with initial error.The experimental results showed that, convection current is carried out using the priori tropospheric zenith delay of two kinds of precision above
Layer error correction, correction residual error is smaller, in the case where inaccurate known tropospheric zenith delay spatial distribution, using experience
2.3 meters of tropospheric zenith delay carry out the InSAR system tropospheric error correction of navigating of spaceborne lists and good correction also can be obtained tie
Fruit.Using the present invention can effectively correct tropospheric propagation to spaceborne list navigated InSAR system positioning result introducing three-dimensional position
Offset, correction accuracy can meet the application demand of high-precision mapping.
Claims (1)
1. a kind of spaceborne list navigated, InSAR system tropospheric error bearing calibration, InSAR refer to interference synthetic aperture radar,
Spaceborne list navigated orbit coordinate A (x of the auxiliary star under body-fixed coordinate system of InSAR system under knowing at any timeA,yA,zA), it is auxiliary
Radar incidence angle θ;The tested corresponding tropospheric zenith time delay τ of scenezenith;Any point is in body-fixed coordinate system in tested scene
Coordinate B (x under systemB,yB,zB), which is characterized in that
It is B (x for coordinate in tested sceneB,yB,zB) arbitrary point, using following steps complete scene three-dimensional coordinate convection current
Layer error influences correction:
Step 1: calculating the horizontal measurement error and vertical measurement error in distance-height plane;
Using following formula calculate troposphere to spaceborne list navigated InSAR system introducing horizontal measurement errorIt is missed with vertical measurement
Difference
In above formula, c indicates the light velocity;
Step 2: two-dimensional measurement error pro is transformed to be indicated under three-dimensional body-fixed coordinate system:
Firstly, calculating vertical measurement error using following formulaError vector under body-fixed coordinate system
Then, horizontal measurement error is calculatedError vector under body-fixed coordinate system
In above formula, O is that coordinate of the earth's core under body-fixed coordinate system is (0,0,0);
Finally, calculating the body-fixed coordinate system three-dimensional coordinate offset vector that troposphere introduces according to the following formula
Step 3: correction tropospheric error is tested the three-dimensional body-fixed coordinate system coordinate of scene after influencing;
The three-dimensional body-fixed coordinate system B'(x that measured point in scene is tested after correcting is calculated using following formulaB',yB',zB'):
xB'=xB+xBE
yB'=yB+yBE
zB'=zB+zBE
Correct result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711250005.XA CN107991676B (en) | 2017-12-01 | 2017-12-01 | Troposphere error correction method of satellite-borne single-navigation-pass InSAR system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711250005.XA CN107991676B (en) | 2017-12-01 | 2017-12-01 | Troposphere error correction method of satellite-borne single-navigation-pass InSAR system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107991676A CN107991676A (en) | 2018-05-04 |
CN107991676B true CN107991676B (en) | 2019-09-13 |
Family
ID=62035042
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711250005.XA Active CN107991676B (en) | 2017-12-01 | 2017-12-01 | Troposphere error correction method of satellite-borne single-navigation-pass InSAR system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107991676B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109471104B (en) * | 2018-11-15 | 2020-08-04 | 首都师范大学 | Method for acquiring three-dimensional movement amount of earth surface from SAR data of two parallel tracks |
CN109583132B (en) * | 2018-12-18 | 2023-04-07 | 中国人民解放军国防科技大学 | Hypersonic aircraft landing probability calculation method based on radar detection error |
CN109782282B (en) * | 2019-03-13 | 2023-04-18 | 武汉大学 | Time series InSAR analysis method integrating troposphere atmospheric delay correction |
CN111239736B (en) * | 2020-03-19 | 2022-02-11 | 中南大学 | Single-baseline-based surface elevation correction method, device, equipment and storage medium |
CN112344847B (en) * | 2020-11-20 | 2022-03-08 | 中国有色金属长沙勘察设计研究院有限公司 | Ground-based synthetic aperture radar data noise reduction method |
CN116859391B (en) * | 2023-09-05 | 2023-11-10 | 中国科学院空天信息创新研究院 | InSAR troposphere delay correction method based on high-frequency texture cross correlation |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004095055A8 (en) * | 2003-04-17 | 2005-03-03 | Secr Defence | Correction of troposphere induced errors in global positioning systems |
CN101770027A (en) * | 2010-02-05 | 2010-07-07 | 河海大学 | Ground surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion |
CN102243300A (en) * | 2011-04-15 | 2011-11-16 | 中国人民解放军国防科学技术大学 | Method for radio frequency interference suppression and error correction of low-frequency synthetic aperture radar |
CN103033833A (en) * | 2012-12-13 | 2013-04-10 | 中国航天科工信息技术研究院 | Method of correcting troposphere delaying errors |
CN103364766B (en) * | 2012-04-01 | 2015-02-04 | 中国科学院电子学研究所 | External calibration method of satellite-borne interferometric synthetic aperture radar (InSAR) system |
CN106600551A (en) * | 2016-12-03 | 2017-04-26 | 中国人民解放军国防科学技术大学 | High-precision geometry geometric correction method for large-scene spaceborne SAR image |
-
2017
- 2017-12-01 CN CN201711250005.XA patent/CN107991676B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004095055A8 (en) * | 2003-04-17 | 2005-03-03 | Secr Defence | Correction of troposphere induced errors in global positioning systems |
CN101770027A (en) * | 2010-02-05 | 2010-07-07 | 河海大学 | Ground surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion |
CN102243300A (en) * | 2011-04-15 | 2011-11-16 | 中国人民解放军国防科学技术大学 | Method for radio frequency interference suppression and error correction of low-frequency synthetic aperture radar |
CN103364766B (en) * | 2012-04-01 | 2015-02-04 | 中国科学院电子学研究所 | External calibration method of satellite-borne interferometric synthetic aperture radar (InSAR) system |
CN103033833A (en) * | 2012-12-13 | 2013-04-10 | 中国航天科工信息技术研究院 | Method of correcting troposphere delaying errors |
CN106600551A (en) * | 2016-12-03 | 2017-04-26 | 中国人民解放军国防科学技术大学 | High-precision geometry geometric correction method for large-scene spaceborne SAR image |
Also Published As
Publication number | Publication date |
---|---|
CN107991676A (en) | 2018-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107991676B (en) | Troposphere error correction method of satellite-borne single-navigation-pass InSAR system | |
US20200103530A1 (en) | Method for extracting elevation control point with assistance of satellite laser altimetry data | |
CN110058236B (en) | InSAR and GNSS weighting method oriented to three-dimensional surface deformation estimation | |
CN107367716B (en) | A kind of high-precision satellite-borne SAR geometric calibration method | |
CN109917356B (en) | Airborne laser scanning system error calibration method | |
CN107272039B (en) | A kind of positioning survey attitude positioning method based on double antenna GPS | |
CN103323855B (en) | A kind of precision acquisition methods of baseline dynamic measurement system | |
CN106556822B (en) | Spaceborne Sliding spotlight SAR pointing accuracy Orbital detection method | |
CN101876701B (en) | Positioning method of remote sensing image of side-looking radar | |
CN107202582A (en) | The in-orbit real time imagery navigation of fixed statellite and method for registering | |
CN110031877B (en) | GRNN model-based regional NWP troposphere delay correction method | |
CN112305510A (en) | DEM matching-based synthetic aperture radar image geometric calibration method | |
Tang et al. | Verification of ZY-3 satellite imagery geometric accuracy without ground control points | |
CN102243299B (en) | Image orthographic correction device of unmanned airborne SAR (Synthetic Aperture Radar) | |
CN110907932A (en) | Distributed InSAR satellite height measurement precision influence factor analysis method and system | |
CN103995260A (en) | Synthetic aperture radar SAR imaging method and device | |
CN108562900B (en) | SAR image geometric registration method based on elevation correction | |
Liu et al. | Accurate mapping method for UAV photogrammetry without ground control points in the map projection frame | |
CN109239710A (en) | Acquisition methods and device, the computer readable storage medium of radar elevation information | |
Wang et al. | Large-scale orthorectification of GF-3 SAR images without ground control points for China’s land area | |
CN108008367B (en) | Ionosphere error correction method for satellite-borne single-navigation-pass InSAR system | |
Li et al. | Improve the ZY-3 height accuracy using ICESat/GLAS laser altimeter data | |
CN108205128B (en) | Passive distance measurement method based on long baseline interferometer | |
CN111127334B (en) | SAR image real-time geometric correction method and system based on RD plane pixel mapping | |
CN109946682B (en) | GF3 data baseline estimation method based on ICESat/GLAS |
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 |