CN104867139B - A kind of remote sensing image clouds and shadow detection method based on radiation field - Google Patents
A kind of remote sensing image clouds and shadow detection method based on radiation field Download PDFInfo
- Publication number
- CN104867139B CN104867139B CN201510239763.6A CN201510239763A CN104867139B CN 104867139 B CN104867139 B CN 104867139B CN 201510239763 A CN201510239763 A CN 201510239763A CN 104867139 B CN104867139 B CN 104867139B
- Authority
- CN
- China
- Prior art keywords
- numerical value
- remote sensing
- sensing image
- radiation field
- detected
- 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
Abstract
The invention discloses a kind of remote sensing image clouds and shadow detection method based on radiation field.This method includes：Obtain atmospheric radiation transmission corresponding to remote sensing image to be detected, utilize the model, using the atmospheric temperature moisture profile under the conditions of earth surface albedo product and default clear sky and aerosol optical depth as driving data, simulation radiation field corresponding to remote sensing image to be detected is generated；According to remote sensing image Ratio index image corresponding with simulation computation of radiation field remote sensing image；Numerical value corresponding to clear sky unshadowed area on Ratio index image is extracted, obtains numerical value set；The first numerical value and second value in numerical value set are obtained, wherein, the first numerical value is more than second value；Density slice is carried out using the first numerical value and second value reduced value index image, to obtain the classification of cloud and shade.By means of the invention it is possible to provide a kind of detection method of cloud versatile, precision is high and shade.
Description
Technical field
The present invention relates to technical field of remote sensing image processing, specifically, Yun Jiyin more particularly on remote sensing image
Shadow universal test method, namely a kind of remote sensing image clouds and shadow detection method based on radiation field.
Background technology
Remote sensing image processing and a key link of analysis are to produce the cloud mask (binary map that mark cloud whether there is
Picture), to determine the position of cloud and shade, these information go further analysis to have cloud sector domain or clear sky area image special for user
Sign, extract required information.
The detection method of cloud and shade mainly includes following two methods on current remote sensing images：The first：Believe by spectrum
The detection method of breath；Second, the geometrical relationship based on cloud-sensor-sun determines the detection method of cloud and shade.Wherein,
In first method, mainly enter the detection racked with shade using specific wave band, but and not all sensor have
Available for the wave band of cloud detection, for example, current most typical cloud detection algorithm is MODIS cloud mask algorithms, it makes use of a lot
The cooperation of wave band.And in general image, do not possess these typical wave bands if the satellite datas such as TM, SPOT, therefore, believe by spectrum
The method of cloud and shade is not general on breath detection remote sensing images, in addition, this method is also easy to cause mixed point of snow-cloud, water-shade
Mixed point, so as to cause the accuracy of detection low.In the second approach, it is necessary to which the parameter such as cloud level is as necessary input, but
These parameters are but difficult to directly obtain, and therefore, the geometrical relationship based on cloud-sensor-sun determines the detection side of cloud and shade
Method is in the case of can not obtain the parameters such as the cloud level and does not apply to, and the versatility of this method is also bad, further, since the cloud obtained
High parameter accuracy is poor, so cloud and shadow Detection have a very big uncertainty, the determination at cloud and Shadow edge position also compared with
Hardly possible is accurately portrayed.
On the detection method of cloud and shade does not possess the problem of versatility on remote sensing images in the prior art, at present not yet
It is proposed effective solution method.
The content of the invention
It is a primary object of the present invention to provide a kind of remote sensing image clouds and shadow detection method based on radiation field, with solution
Certainly the detection method of cloud and shade does not possess the problem of versatility on remote sensing images in the prior art.
According to one aspect of the present invention, there is provided a kind of remote sensing image clouds and shadow detection method based on radiation field,
This method includes：Obtain atmospheric radiation transmission corresponding to remote sensing image to be detected；With earth surface albedo product and
The aerosol optical depth under the conditions of atmospheric temperature-moisture profile, default clear sky under the conditions of default clear sky is driving number
According to being generated using atmospheric radiation transmission and simulate radiation field corresponding to remote sensing image to be detected；According to be detected
Remote sensing image Ratio index image corresponding with the remote sensing image that simulation computation of radiation field is to be detected；Extraction ratio refers to
Numerical value corresponding to clear sky unshadowed area on number image, obtains numerical value set；Obtain the first numerical value and second in numerical value set
Numerical value, wherein, the first numerical value is more than second value；Line density point is entered using the first numerical value and second value reduced value index image
Cut, to obtain the classification of cloud and shade.
Further, according to the remote sensing image to be detected remote sensing image to be detected with simulation computation of radiation field
The step of corresponding Ratio index image, specifically includes：By remote sensing image to be detected and each wave band of simulation radiation field
Radiance seeks ratio；The each ratio tried to achieve is subjected to product calculation, to obtain corresponding to remote sensing image to be detected
The Ratio index image of cloud and shade.
Further, numerical value corresponding to clear sky unshadowed area on Ratio index image is extracted, obtains the step of numerical value set
Suddenly specifically include：Remote sensing image to be detected is calculated with simulating the similarity of radiation field；On Ratio index image, it is determined that
The region that similarity is more than default similarity is clear sky unshadowed area；Numerical value corresponding to extracting clear sky unshadowed area, is obtained
Numerical value set.
Further, calculate remote sensing image to be detected and simulate radiation field similarity the step of be specially：Adopt
Remote sensing image to be detected and the similarity of simulation radiation field are calculated with structural similarity texture recognition algorithm.
Further, obtain the first numerical value in numerical value set and second value specifically includes：Number in logarithm value set
Value is ranked up according to size；Remove the numerical value for front and back end 2% of being sorted in numerical value set；In numerical value set, acquisition removes 2%
Numerical value after maximum as the first numerical value, minimum value is as second value.
Further, density slice is carried out using the first numerical value and second value reduced value index image, with obtain cloud and
The classification of shade specifically includes：First numerical value and second value are carried out close as initial segmentation threshold value, reduced value index image
Degree segmentation, to obtain the preliminary classification of cloud and shade；Initial segmentation threshold value is adjusted according to the classification results of preliminary classification, using tune
Segmentation threshold reduced value index image after whole carries out density slice, to obtain cloud mask.
Further, remote sensing image to be detected can be MODIS images, TM/ETM+ images, VIIRS images,
SPOT images, ASTER images, AVHRR images or HJ-1A/B images etc..
By the present invention, using the atmospheric radiation transmission corresponding with remote sensing image to be detected, and with ground
The aerosol light under the conditions of atmospheric temperature-moisture profile, default clear sky under the conditions of table albedo product and default clear sky
Thickness is driving data, and simulation generates simulation radiation field corresponding to remote sensing image to be detected, then according to be detected
The remote sensing image remote sensing image corresponding Ratio index image to be detected with simulation computation of radiation field, extract ratio
Numerical value corresponding to clear sky unshadowed area on index image, obtains numerical value set, obtains the first numerical value in numerical value set and the
Two numerical value, wherein, the first numerical value is more than second value, is finally carried out using the first numerical value and second value reduced value index image
Density slice, to obtain the classification of cloud and shade, in this method, the radiation field of simulation is directly used come for reference, Ke Yitong
When cloud and shade are detected, and not by earth's surface snow lid, water body influenceed, there are versatility, nearly all remote sensing images
It can directly apply, not limited by wave band, the detection method of cloud and shade does not possess on remote sensing images in the prior art for solution
The problem of versatility.
Described above is only the general introduction of technical solution of the present invention, in order to better understand the technological means of the present invention,
And can be practiced according to the content of specification, and in order to allow above and other objects of the present invention, feature and advantage can
Become apparent, below especially exemplified by the embodiment of the present invention.
Brief description of the drawings
By reading the detailed description of hereafter preferred embodiment, it is various other the advantages of and benefit it is common for this area
Technical staff will be clear understanding.Accompanying drawing is only used for showing the purpose of preferred embodiment, and is not considered as to the present invention
Limitation.And in whole accompanying drawing, identical process is denoted by the same reference numerals.In the accompanying drawings：
Fig. 1 is the flow of remote sensing image clouds based on radiation field and shadow detection method according to a first embodiment of the present invention
Figure；
Fig. 2 is the flow of remote sensing image clouds based on radiation field and shadow detection method according to a second embodiment of the present invention
Figure.
Embodiment
The present invention will be further described with reference to the accompanying drawings and detailed description.It is pointed out that do not conflicting
In the case of, the feature in embodiment and embodiment in the application can be mutually combined.
Embodiment one
First, the embodiment of the present invention one provides a kind of remote sensing image clouds and shadow detection method based on radiation field, ginseng
See Fig. 1, this method may comprise steps of：
Step S102：Obtain atmospheric radiation transmission corresponding to remote sensing image to be detected.
In this step, according to the species of remote sensing image to be detected, obtain corresponding atmospheric radiation and pass
Defeated model, for example, remote sensing image to be detected is MODIS images, then the atmospheric radiation transmission obtained is MODTRAN
Radiative transfer model.
Step S104：With atmospheric temperature-moisture profile under the conditions of earth surface albedo product and default clear sky, default
Aerosol optical depth under the conditions of clear sky is driving data, and optical remote sensing to be detected is generated using atmospheric radiation transmission
Simulation radiation field corresponding to image.
Wherein, earth surface albedo product preferably uses MCD43B3 earth surface albedo product.Under the conditions of default clear sky
Aerosol optical depth under the conditions of atmospheric temperature-moisture profile, default clear sky is the experience ginseng under the conditions of typical clear sky
Number.
Step S106：According to the remote sensing image that remote sensing image to be detected is to be detected with simulation computation of radiation field
Corresponding Ratio index image.
Specifically, remote sensing image to be detected and each wave band radiance of simulation radiation field are sought into ratio, so
Product calculation is carried out to each ratio afterwards, the result of product calculation is Ratio index corresponding to remote sensing image to be detected
Image.
Step S108：Numerical value corresponding to clear sky unshadowed area on Ratio index image is extracted, obtains numerical value set.
Specifically, it is first determined clear sky unshadowed area：Calculated using structural similarity texture recognition algorithm to be detected
Remote sensing image and the similarity of simulation radiation field, on Ratio index image, determine that similarity is more than default similarity
Region is clear sky unshadowed area, and then numerical value corresponding to extraction clear sky unshadowed area, obtains numerical value set.
Step S110：The first numerical value and second value in numerical value set are obtained, wherein, the first numerical value is more than the second number
Value.
Specifically, handled using following steps logarithm value set, and obtain the first numerical value and second value：It is right first
Numerical value in numerical value set is ranked up according to size, is then removed the numerical value for front and back end 2% of being sorted in numerical value set, is finally existed
The maximum after the numerical value for removing 2% is obtained in numerical value set as the first numerical value, minimum value is as second value.
Step S112：Density slice is carried out using the first numerical value and second value reduced value index image, with obtain cloud and
The classification of shade.
Specifically, the first numerical value and second value can be entered line density as initial segmentation threshold value, reduced value index image
Segmentation, to obtain the preliminary classification of cloud and shade, initial segmentation threshold value then is adjusted according to the classification results of preliminary classification, is used
Segmentation threshold reduced value index image after adjustment carries out density slice again, to obtain cloud mask.
Because cloud in reflected waveband shows high reflection characteristic, and cloud shade shows low reflection characteristic in reflected waveband,
The two is the extreme case of reflectance spectrum, if existing without cloud or cloud shade, wave band radiation is in therebetween.Therefore, exist
In the embodiment, it is used as reference by simulating the radiation field under the conditions of clear sky, cloud and cloud shade can be distinguished, and can be simultaneously right
Cloud and shade are detected, and are not influenceed by earth's surface snow lid, water body, have a versatility, nearly all remote sensing images, including
MODIS images, TM/ETM+ images, VIIRS images, SPOT images, ASTER images, AVHRR images or HJ-1A/B images, all
Can directly apply, not limited by wave band, solve in the prior art on remote sensing images the detection method of cloud and shade do not possess it is logical
With property the problem of.
Embodiment two
The embodiment is the further preferred remote sensing image clouds and shade based on radiation field on the basis of embodiment one
Detection method, in the method, using remote sensing image to be detected as MODIS image examples, referring to Fig. 2, what this method included
Step is described as follows：
(1) using MODTRAN radiative transfer models, MCD43B3 earth surface albedo product, typical Cloudless atmosphere temperature-
Moisture profile and aerosol optical depth etc. are driving data, 7 ripples of the MODIS shortwave radiations under the conditions of simulation generation clear sky
Section radiance.
(2) using 7 wave band radiances of MODIS as ambient field, the cloud detection index of short-wave band is built, MODIS is schemed
As (to be detected) and the radiation field subrane of simulation seek ratio, then carry out seizing the opportunity computing obtaining the ratio of cloud and shade to each ratio
Value index number image；
(3) structural similarity texture recognition algorithm (SSIM) is used to the radiation field of MODIS images (to be detected) and simulation
The similarity of the two is calculated, similarity is considered clear sky unshadowed area more than 0.9.
(4) be extracted on locus and the above-mentioned clear sky unshadowed area obtained using texture recognition algorithm corresponding to compare
Numerical value on value index number image, is ranked up according to size, removes the numerical value of front and back end 2%, then obtains minimum and maximum threshold
Value, the threshold value is then as subsequently entering the preliminary threshold that racks with shadow Detection.
(5) density slice is carried out using preliminary threshold reduced value index image, obtains the preliminary classification of cloud and shade, then
According to classification situation adjustment threshold value, the reliable cloud mask of precision, shadow image are obtained, while realizing remote sensing image clouds and shade
Detection.
Although here by taking MODIS images as an example, this method is applicable to any remote sensing image, including, TM/ETM
+, VIIRS, SPOT, a series of classical satellite images such as ASTER, AVHRR, HJ-1A/B.
The method provided using the embodiment, it is possible to achieve the high precision test of cloud and shade, so as to be remote sensing images point
Class research, remote sensing images inverting, which are studied and carries out atmospheric research using remote sensing images, provides key technology and authentic data branch
Hold.
From the description of various embodiments above, it can be seen that the embodiment of the present invention realizes following technique effect：
(1) the cloud detection technology based on radiomimesis field, spectral information and observation geological information are eliminated the reliance on, is a kind of complete
New cloud detection thinking；
(2) cloud detection technology different from the past, this method detect while can realizing cloud and shade, not by underlying surface height
Reflect the limitation of earth's surface (e.g., ice and snow, desert etc.) and low reflection earth's surface (e.g., water body) etc.；
(3) compared with being only applicable to the detection method of specific remote sensing images in the prior art, this method goes for appointing
What remote sensing image, not by wave band any restrictions, has extremely strong versatility, and precision is reliable.
Each embodiment in this specification is described by the way of progressive, what each embodiment stressed be with
The difference of other embodiment, between each embodiment identical similar part mutually referring to.
The foregoing is only a preferred embodiment of the present invention, but protection scope of the present invention be not limited thereto,
Any people for being familiar with the technology disclosed herein technical scope in, the change or replacement that can readily occur in should all be covered
Within protection scope of the present invention.Therefore, protection scope of the present invention should be defined by scope of the claims.
Claims (5)
1. a kind of remote sensing image clouds and shadow detection method based on radiation field, it is characterised in that including：
Obtain atmospheric radiation transmission corresponding to remote sensing image to be detected；
Under the conditions of atmospheric temperature-moisture profile under the conditions of earth surface albedo product and default clear sky, default clear sky
Aerosol optical depth is driving data, and the remote sensing image to be detected is generated using the atmospheric radiation transmission
Corresponding simulation radiation field；
The remote sensing image to be detected according to the remote sensing image to be detected with the simulation computation of radiation field
Corresponding Ratio index image；Wherein, by the remote sensing image to be detected and each wave band of the simulation radiation field
Radiance seeks ratio；The each ratio tried to achieve is subjected to product calculation, to obtain the remote sensing image pair to be detected
The Ratio index image answered；
Numerical value corresponding to clear sky unshadowed area on the Ratio index image is extracted, obtains numerical value set；
The first numerical value and second value in the numerical value set are obtained, wherein, first numerical value is more than the second value；
The first numerical value and second value in the numerical value set are obtained, is specifically included：Numerical value in the numerical value set is pressed
It is ranked up according to size；Remove the numerical value for front and back end 2% of being sorted in the numerical value set；In the numerical value set, acquisition is removed
Maximum after 2% numerical value is as first numerical value, and minimum value is as the second value；
Density slice is carried out to the Ratio index image using first numerical value and the second value, to obtain Yun Heyin
The classification of shadow.
2. remote sensing image clouds and shadow detection method according to claim 1 based on radiation field, it is characterised in that extraction
Numerical value corresponding to clear sky unshadowed area on the Ratio index image, the step of obtaining numerical value set, specifically include：
Calculate the remote sensing image to be detected and the similarity of the simulation radiation field；
On the Ratio index image, determine that the region that the similarity is more than default similarity is clear sky unshadowed area；
Numerical value corresponding to extracting the clear sky unshadowed area, obtain the numerical value set.
3. remote sensing image clouds and shadow detection method according to claim 2 based on radiation field, it is characterised in that calculate
The remote sensing image to be detected and it is described simulation radiation field similarity the step of be specially：
The remote sensing image to be detected and the simulation radiation field are calculated using structural similarity texture recognition algorithm
Similarity.
4. remote sensing image clouds and shadow detection method according to claim 1 based on radiation field, it is characterised in that use
First numerical value and the second value carry out density slice to the Ratio index image, to obtain the classification of cloud and shade
Specifically include：
Using first numerical value and the second value as initial segmentation threshold value, line density point is entered to the Ratio index image
Cut, to obtain the preliminary classification of cloud and shade；
The initial segmentation threshold value is adjusted according to the classification results of the preliminary classification, using the segmentation threshold after adjustment to described
Ratio index image carries out density slice, to obtain cloud mask.
5. remote sensing image clouds and shadow detection method according to claim 1 based on radiation field, it is characterised in that described
Remote sensing image to be detected can be MODIS images, TM/ETM+ images, VIIRS images, SPOT images, ASTER images,
AVHRR images or HJ-1A/B images.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510239763.6A CN104867139B (en) | 2015-05-12 | 2015-05-12 | A kind of remote sensing image clouds and shadow detection method based on radiation field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510239763.6A CN104867139B (en) | 2015-05-12 | 2015-05-12 | A kind of remote sensing image clouds and shadow detection method based on radiation field |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104867139A CN104867139A (en) | 2015-08-26 |
| CN104867139B true CN104867139B (en) | 2018-02-09 |
Family
ID=53912954
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510239763.6A Active CN104867139B (en) | 2015-05-12 | 2015-05-12 | A kind of remote sensing image clouds and shadow detection method based on radiation field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104867139B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105426903A (en) * | 2015-10-27 | 2016-03-23 | 航天恒星科技有限公司 | Cloud determination method and system for remote sensing satellite images |
| CN105678777B (en) * | 2016-01-12 | 2018-03-13 | 武汉大学 | A kind of optical satellite image cloud of multiple features combining and cloud shadow detection method |
| CN107944413A (en) * | 2017-12-04 | 2018-04-20 | 中国科学院南京地理与湖泊研究所 | Aquatic vegetation Classification in Remote Sensing Image threshold value calculation method based on spectral index ranking method |
| AU2019362019A1 (en) * | 2018-10-19 | 2021-05-20 | Climate Llc | Machine learning techniques for identifying clouds and cloud shadows in satellite imagery |
| CN110599488B (en) * | 2019-09-27 | 2022-04-29 | 广西师范大学 | Cloud detection method based on Sentinel-2 aerosol wave band |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101034477A (en) * | 2007-03-29 | 2007-09-12 | 上海大学 | Method for eliminating shadow on remote sensing digital image and recovering picture element remote sensing value in shadow |
| CN103901420A (en) * | 2014-04-18 | 2014-07-02 | 山东科技大学 | Method for dynamic threshold method remote sensing data cloud identification supported by prior surface reflectance |
| CN104156923A (en) * | 2014-08-12 | 2014-11-19 | 西北工业大学 | Multispectral remote sensing image cloud removing method based on sparse representation |
-
2015
- 2015-05-12 CN CN201510239763.6A patent/CN104867139B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101034477A (en) * | 2007-03-29 | 2007-09-12 | 上海大学 | Method for eliminating shadow on remote sensing digital image and recovering picture element remote sensing value in shadow |
| CN103901420A (en) * | 2014-04-18 | 2014-07-02 | 山东科技大学 | Method for dynamic threshold method remote sensing data cloud identification supported by prior surface reflectance |
| CN104156923A (en) * | 2014-08-12 | 2014-11-19 | 西北工业大学 | Multispectral remote sensing image cloud removing method based on sparse representation |
Non-Patent Citations (2)
| Title |
|---|
| 基于 MODIS 数据的高分辨率气溶胶光学厚度反演;姬大彬;《中国优秀硕士学位论文全文数据库 信息科技辑》;20111215;I140-836 * |
| 基于MPHD算法的高分辨率遥感影像对海洋上空的云及其阴影的识别与匹配;陈晓东 等;《海洋学研究》;20090615;第27卷(第2期);51-57 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104867139A (en) | 2015-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104867139B (en) | A kind of remote sensing image clouds and shadow detection method based on radiation field | |
| Makarau et al. | Haze detection and removal in remotely sensed multispectral imagery | |
| CN108051371B (en) | A kind of shadow extraction method of ecology-oriented environment parameter remote-sensing inversion | |
| KR101922645B1 (en) | cloud area detection device and cloud area detection method | |
| CN109101955A (en) | Industrial heat anomaly area recognizing method based on Multi-sensor satellite remote sensing | |
| CN110648347A (en) | Coastline extraction method and device based on remote sensing image | |
| JP6885509B2 (en) | Image processing device, image processing method and image processing program | |
| CN107203742B (en) | Gesture recognition method and device based on significant feature point extraction | |
| CN107421892A (en) | A kind of hyperspectral data processing method for water body information | |
| CN104408705A (en) | Anomaly detection method of hyperspectral image | |
| Zhang et al. | Extraction of tree crowns damaged by Dendrolimus tabulaeformis Tsai et Liu via spectral-spatial classification using UAV-based hyperspectral images | |
| CN103389255A (en) | Method for predicting water content of tea on basis of near-infrared hyperspectral textural feature modeling | |
| CN103839267A (en) | Building extracting method based on morphological building indexes | |
| CN111339948A (en) | Automatic identification method for newly-added buildings of high-resolution remote sensing images | |
| WO2020027167A1 (en) | System, method, and non-transitory, computer-readable medium containing instructions for image processing | |
| WO2019184269A1 (en) | Landsat 8 snow-containing image-based cloud detection method | |
| CN114049566B (en) | Method and device for detecting cloud and cloud shadow of land satellite image in step-by-step refinement manner | |
| CN109580497B (en) | Hyperspectral mineral abnormal information extraction method based on singularity theory | |
| CN111310771A (en) | Road image extraction method, device and equipment of remote sensing image and storage medium | |
| Ali et al. | A Modified Built-up Index (MBI) for automatic urban area extraction from Landsat 8 Imagery | |
| CN110738134B (en) | Soil information extraction method and device for visible light image of unmanned aerial vehicle | |
| CN104966273B (en) | Haze method is gone suitable for the DCM-HTM of optical remote sensing image | |
| Ahmad et al. | Haze reduction from remotely sensed data | |
| CN110717413B (en) | Unmanned aerial vehicle visible light image-oriented water body information extraction method and device | |
| CN110378294B (en) | Hyperspectral target detection method and system based on local energy constraint and feature fusion |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |