CN109544494B - 一种人体安检中被动毫米波图像与可见光图像的融合方法 - Google Patents
一种人体安检中被动毫米波图像与可见光图像的融合方法 Download PDFInfo
- Publication number
- CN109544494B CN109544494B CN201811340180.2A CN201811340180A CN109544494B CN 109544494 B CN109544494 B CN 109544494B CN 201811340180 A CN201811340180 A CN 201811340180A CN 109544494 B CN109544494 B CN 109544494B
- Authority
- CN
- China
- Prior art keywords
- image
- millimeter wave
- suspicious
- visible light
- frequency coefficient
- 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
- 238000007500 overflow downdraw method Methods 0.000 title claims description 6
- 239000011159 matrix material Substances 0.000 claims abstract description 48
- 230000004927 fusion Effects 0.000 claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 22
- 230000011218 segmentation Effects 0.000 claims abstract description 14
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 238000007689 inspection Methods 0.000 claims abstract description 4
- 239000002184 metal Substances 0.000 claims abstract description 4
- 239000007788 liquid Substances 0.000 claims abstract description 3
- 210000000746 body region Anatomy 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000010187 selection method Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 230000007547 defect Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 239000003502 gasoline Substances 0.000 description 2
- 238000000889 atomisation Methods 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000001815 facial effect Effects 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000003331 infrared imaging Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种人体安检中被动毫米波图像与可见光图像的融合方法,步骤:确定原始毫米波图像中的金属等低亮温可疑目标,显示为暗色区域;液体等高亮温可疑目标,显示为亮色区域,在最终融合图像中可疑目标均以高亮区域表示;对于暗色可疑目标,将毫米波图像和可见光图像进行小波分解后,融合图像的低频系数矩阵中目标区域来自毫米波图像的低频系数矩阵,非目标区域来自可见光图像的低频系数矩阵,融合图像的高频系数矩阵为原始可见光图像的高频系数矩阵,最终使用小波变换获得融合图像;对于亮色可疑目标,采用固定阈值分割获取毫米波图像可疑目标,通过区域选择融合得到最终的融合图像:目标区域来自毫米波图像,非目标区域来自可见光图像。
Description
技术领域
本发明提供一种用于毫米波图像与可见光图像的融合方法,特别是涉及一种被动毫米波人体安检图像与可见光图像的融合方法,属于毫米波人体安检领域。
背景技术
随着全球范围内的恐怖主义和暴力活动的不断升级,世界各国对人体隐匿违禁物品检测提出了越来越高的要求。相比于传统的X-ray及红外成像,毫米波成像不仅可以检测出隐匿在人体衣物下的金属物品,还可以检测出炸药、汽油、陶瓷刀具等危险品。相比于主动毫米波成像技术,被动毫米波成像具有成像速度快、对被测人员无辐射、不侵犯被测人员的隐私等优点。因此,被动毫米波成像技术在人体安防检测领域得到了广泛的应用和发展。然而,由于受到天线口径、波长尺寸、接收机不理想性以及系统噪声等影响,使得毫米波图像具有分辨率低、噪声水平高、可读性差等不足,图像中的目标及背景的真实边缘将无法得到有效体现,主要依靠亮度和模糊的几何形状来辨别可疑目标。相比而言光学成像是普遍使用的成熟技术,光学图像具有光谱信息丰富、动态范围大、精确直观的特点,却不像毫米波成像能够发现隐匿的目标。单独使用可见光或被动毫米波图像均存在不足之处,对于这两种具有互补性的图像,图像融合技术能够有效地增强场景信息量,突出目标,有利于更精确识别隐藏的危险物品,实现定位,因此提出基于被动毫米波和可见光的图像融合技术。
针对可见光与毫米波图像融合方法,研究者们已开展许多尝试,如基于多小波包的红外/被动毫米波图像融合(何伏春,聂建英.基于多小波包的红外/被动毫米波图像融合[J].电光与控制,2015,22(06):11-16)、基于多尺度变换的无源毫米波图像融合(TangShi,Xiong Jintao,Yang Jianyu,Peng Xinliang.Research on Image Fusion forVisual and PMMW Images Using NSCT and Region Detection.3rd InternationalConference on Multimedia Technology(ICMT-13).Atlantis Press,2013)等,但都没考虑人体安检应用中被动毫米波图像可疑目标区域的特征与其辐射亮温的关系。本发明针对不同辐射亮温的可疑目标在毫米波图像中的特征具有较大差异的问题,将可疑目标分为亮目标与暗目标两类分别处理,提高实际应用中的稳定性与可靠性。
发明内容
本发明的技术解决问题:克服现有技术的不足,提供一种人体安检中被动毫米波图像与可见光图像的融合方法,通过毫米波图像与可见光图像的融合,增强场景信息量,突出目标,以利于更精确识别隐藏的危险物品,实现定位,弥补毫米波图像具有分辨率低、噪声水平高、可读性差等不足;并同时保证融合算法的稳定有效。
本发明技术解决方案:一种人体安检中被动毫米波图像与可见光图像的融合方法,包括以下步骤:
第一步,统一图像显示方式,确定原始毫米波图像中的金属作为低辐射可疑目标,显示为暗色区域;液体作为高辐射目标,显示为亮色区域,在最终融合图像中可疑目标均以高亮区域表示;
第二步,对毫米波图像中的暗色和亮色可疑目标区域分别进行目标区域提取与分割,获取毫米波图像中可疑目标区域中的信息;对于暗色可疑目标,将毫米波图像进行小波分解后,提取可疑暗色目标区域的毫米波图像低频系数矩阵;对于亮色可疑目标,获取亮色可疑目标区域的毫米波图像;
第三步,采用区域选择则的方法对可见光图像和可疑目标区域中的毫米波图像进行融合,得到最终的融合图像;对于暗色可疑目标,通过区域选择获得融合图像的低频系数矩阵:目标区域来自毫米波图像的低频系数矩阵,非目标区域来自可见光图像的低频系数矩阵;融合图像的高频系数矩阵为原始可见光图像的高频系数矩阵;最后通过二维离散小波逆变换得到含有可疑目标的融合图像。对于亮色可疑目标区域,直接通过区域选择融合得到最终的融合图像:目标区域来自毫米波图像,非目标区域来自可见光图像。
所述第二步实现暗色可疑目标提取与图像分割的过程如下:
(1)对毫米波图像和可见光图像进行小波分解,分别得到毫米波图像和可见光图像低频系数矩阵和高频系数矩阵;
(2)采用最大类间方差法,即Otsu算法获取毫米波图像的低频系数矩阵中人体与背景分割阈值;
(3)对毫米波图像人体区域的低频系数矩阵进行基于引导滤波的增强;
(4)获取增强后的毫米波图像人体区域的低频系数矩阵,同时将目标区域以亮色表示;
(5)采用Otsu算法获取增强后的毫米波图像人体区域的低频系数矩阵中目标区域与非目标区域的分割阈值,通过阈值分割得到暗色可疑目标区域的毫米波图像低频系数矩阵。
所述第二步实现亮色可疑目标提取与图像分割的过程如下:
(1)对毫米波源图像和可见光源图像分别进行归一化处理;
(2)对归一化的毫米波图像进行固定阈值分割,获取亮色可疑目标区域的毫米波图像。
本发明与现有技术相比的优点在于:
目前常用的图像融合方法有像素绝对值取最大和像素加权平均。
采用像素绝对值取最大方法进行毫米波图像与可见光图像融合,融合图像很好的继承了可见光图像中的环境信息和毫米波图像中的可疑目标信息,而且可疑目标的显著程度与增强后的毫米波图像中可疑目标的显著程度相当,可以显著地标示可疑目标的位置;但是可见光图像中人的面部信息几乎全部丢失。结果见附图2。
采用像素加权平均方法进行毫米波图像与可见光图像融合,融合图像中人脸可以辨识,但是没有可见光图像中清晰,有类似雾化的现象;可疑目标可以辨识,但是目标的显著程度远低于增强后毫米波图像中可疑目标的显著程度;而且与可见光图像相比,融合图像的亮度明显降低,视觉效果变差。结果见图2。
采用本发明的方法得到的融合图像很好地继承了可见光图像中的环境及人体信息,人脸清晰,图像整体亮度与可见光图像一致;可疑目标区域特征显著,冗余信息少,可以极大提高安检人员对可疑目标的辨别效率。而且经大量实验验证,本发明融合结果和视觉效果非常稳定,可以用于实际应用。结果见图2,图2中从左到右,依次为毫米波原始图像、可见光图像、基于像素绝对值取大获取的融合图像、基于本发明的方法获取的融合图像。
附图说明
图1为毫米波图像与可见光图像的融合算法流程图;(a)为毫米波图像暗色可疑目标区域与可见光图像的融合图;(b)为毫米波图像暗色可疑目标区域与可见光图像的融合图;
图2为含有暗色可疑目标的毫米波图像与可见光图像的融合结果。
具体实施方式
本发明处理的毫米波原始图像中,刀,枪等低辐射可疑目标为暗色区域,汽油等高辐射目标为亮色区域,但是在最终融合图像中可疑目标均以高亮区域表示。
1.毫米波图像暗色可疑目标区域与可见光图像融合算法
如图1中的(a)所示,本发明毫米波图像暗色可疑目标区域与可见光图像的融合算法包括:
第一步,对毫米波图像和可见光图像进行小波分解,分别得到其低频频系数矩阵和高频系数矩阵。假设毫米波图像为A,可见光图像为B,融合图像为F。对毫米波图像A进行二维离散小波变换,得到毫米波图像的低频系数矩阵AL和高频系数矩阵AH,AV,AD;对可见光图像B进行二维离散小波变换,得到低频系数矩阵BL和高频系数矩阵BH,BV,BD。设毫米波图像的低频系数矩阵AL的大小为m*n(m为AL的行数,n为AL的列数),AL的元素表示为AL(i,j)。
第二步,采用Otsu算法(最大类间方差法)获取毫米波图像的低频系数矩阵AL中人体区域低频系数与背景区域低频系数的分割阈值Threshbody,则可由(1)式获得毫米波图像低频系数矩阵中人体区域与背景区域的二值标记模板:
第三步,对毫米波图像人体区域的低频系数矩阵AL进行基于引导滤波的增强。
增强方法如下:
首先,按照式(2)对AL进行归一化处理得到AL01:
(2)式中,ALmin为AL的元素最大值,ALmax为AL的元素最小值。
然后,将AL01作为输入图和引导图进行引导滤波,输出结果记为AL01’,再由式(3)和式(4)得到毫米波图像中人体归一化低频系数矩阵AL01的增强结果AL2:
(3)式中Times是一个表示低频系数矩阵增强程度的常数,参考范围为[5,8]。
第四步,获取增强后的毫米波图像人体区域的低频系数矩阵ALbody,同时将目标区域以亮色表示,如(5)式所示:
ALbody(i,j)=[1-AL2(i,j)]·Maskbody(i,j) (5)
第五步,采用Otsu算法(最大类间方差法)获取ALbody中目标区域与非目标区域的分割阈值Threshtarget,通过阈值分割由(6)式得到目标区域的归一化低频系数矩阵ALtar01:
然后反归一化得到毫米波图像中目标区域的低频系数矩阵ALtarget,如(7)式所示:
ALtarget(i,j)=ALtar01(i,j)*(ALmax-ALmin) (7)
第六步,通过区域选择获得融合图像F的低频系数矩阵为FL:目标区域选择毫米波图像的低频系数加入FL,非目标区域选择可见光图像的低频系数加入FL,如(8)式所示:
第七步,取FL作为融合图像F的低频系数矩阵,BH,BV,BD作为融合图像F的的高频系数矩阵进行二维离散小波逆变换,得到融合图像F。
2.毫米波图像亮色可疑目标区域与可见光图像融合方法
如图1中(b)所示,本发明毫米波图像暗色可疑目标区域与可见光图像的融合算法包括:
第一步,对毫米波源图像和可见光源图像分别进行归一化处理。
假设毫米波图像为A,可见光图像为B,融合图像为图像F。通过(9)式和(10)式对图像进行归一化处理:
其中ALmax和ALmin为A中元素的最大值和最小值,BLmax和BLmin为B中元素的最大值和最小值。
第二步,对归一化的毫米波图像进行阈值分割。经大量实验结果验证,采用固定阈值分割便能很好的分割去亮目标区域,归一化参考阈值0.65.
第三步,通过区域选择融合得到最终的融合图像。区域选择融合规则如(11)式所示:
提供以上实施例仅仅是为了描述本发明的目的,而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改,均应涵盖在本发明的范围之内。
Claims (3)
1.一种人体安检中被动毫米波图像与可见光图像的融合方法,其特征在于,包括以下步骤:
第一步,统一图像显示方式,确定原始毫米波图像中的金属作为低辐射可疑目标,显示为暗色区域;液体作为高辐射目标,显示为亮色区域,在最终融合图像中可疑目标均以高亮区域表示;
第二步,对毫米波图像中的暗色和亮色可疑目标区域分别进行目标区域提取与分割,获取毫米波图像中可疑目标区域中的信息;对于暗色可疑目标,将毫米波图像进行小波分解后,提取可疑暗色目标区域的毫米波图像低频系数矩阵;对于亮色可疑目标,获取亮色可疑目标区域的毫米波图像;
第三步,采用区域选择的方法对可见光图像和可疑目标区域中的毫米波图像进行融合,得到最终的融合图像;对于暗色可疑目标,通过区域选择获得融合图像的低频系数矩阵:目标区域来自毫米波图像的低频系数矩阵,非目标区域来自可见光图像的低频系数矩阵;融合图像的高频系数矩阵为原始可见光图像的高频系数矩阵;最后通过二维离散小波逆变换得到含有可疑目标的融合图像;对于亮色可疑目标区域,直接通过区域选择融合得到最终的融合图像:目标区域来自毫米波图像,非目标区域来自可见光图像。
2.根据权利要求1所述的一种人体安检中被动毫米波图像与可见光图像的融合方法,其特征在于:所述第二步中,提取可疑暗色目标区域的毫米波图像低频系数矩阵的过程如下:
(1)对毫米波图像和可见光图像进行小波分解,分别得到毫米波图像和可见光图像低频系数矩阵和高频系数矩阵;
(2)采用最大类间方差法,即Otsu算法获取毫米波图像的低频系数矩阵中人体与背景分割阈值;
(3)对毫米波图像人体区域的低频系数矩阵进行基于引导滤波的增强;
(4)获取增强后的毫米波图像人体区域的低频系数矩阵,同时将目标区域以亮色表示;
(5)采用Otsu算法获取增强后的毫米波图像人体区域的低频系数矩阵中目标区域与非目标区域的分割阈值,通过阈值分割得到暗色可疑目标区域的毫米波图像低频系数矩阵。
3.根据权利要求1所述的一种人体安检中被动毫米波图像与可见光图像的融合方法,其特征在于:所述第二步中,获取亮色可疑目标区域的毫米波图像的过程如下:
(1)对毫米波源图像和可见光源图像分别进行归一化处理;
(2)对归一化的毫米波图像进行固定阈值分割,获取亮色可疑目标区域的毫米波图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811340180.2A CN109544494B (zh) | 2018-11-12 | 2018-11-12 | 一种人体安检中被动毫米波图像与可见光图像的融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811340180.2A CN109544494B (zh) | 2018-11-12 | 2018-11-12 | 一种人体安检中被动毫米波图像与可见光图像的融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109544494A CN109544494A (zh) | 2019-03-29 |
CN109544494B true CN109544494B (zh) | 2021-08-17 |
Family
ID=65846885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811340180.2A Active CN109544494B (zh) | 2018-11-12 | 2018-11-12 | 一种人体安检中被动毫米波图像与可见光图像的融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109544494B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110211081B (zh) * | 2019-05-24 | 2023-05-16 | 南昌航空大学 | 一种基于图像属性和引导滤波的多模医学图像融合方法 |
CN114581944B (zh) * | 2022-02-18 | 2022-11-04 | 杭州睿影科技有限公司 | 一种毫米波图像处理方法、装置及电子设备 |
CN117495695B (zh) * | 2023-11-10 | 2024-05-03 | 苏州清研浩远汽车科技有限公司 | 基于毫米波和红外图像融合的低光照环境检测系统 |
CN117528274B (zh) * | 2023-11-13 | 2024-05-14 | 北京信息科技大学 | 基于深度学习的红外成像非均匀噪声提取及优化方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1282132C (zh) * | 2004-05-20 | 2006-10-25 | 上海交通大学 | 基于运动目标检测的红外和可见光动态图像融合方法 |
US8102306B2 (en) * | 2010-05-13 | 2012-01-24 | The United States Of America As Represented By The Secretary Of The Navy | Active-radar-assisted passive composite imagery for aiding navigation or detecting threats |
KR20130023930A (ko) * | 2011-08-30 | 2013-03-08 | 대구대학교 산학협력단 | 은닉물체 탐지용 수동형 밀리미터파 영상과 가시 영상의 정합 및 융합방법 |
CN103777252B (zh) * | 2014-02-24 | 2016-08-24 | 黄晓鹏 | 危险品检测方法 |
CN106384341B (zh) * | 2016-09-12 | 2019-11-15 | 南京理工大学 | 一种基于目标极化辐射特性的毫米波被动图像增强方法 |
CN106780392B (zh) * | 2016-12-27 | 2020-10-02 | 浙江大华技术股份有限公司 | 一种图像融合方法及装置 |
CN107202983B (zh) * | 2017-05-19 | 2020-11-13 | 深圳佑驾创新科技有限公司 | 基于图像识别和毫米波雷达融合的自动刹车方法和系统 |
CN108198157A (zh) * | 2017-12-22 | 2018-06-22 | 湖南源信光电科技股份有限公司 | 基于显著目标区域提取和nsst的异源图像融合方法 |
-
2018
- 2018-11-12 CN CN201811340180.2A patent/CN109544494B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109544494A (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109544494B (zh) | 一种人体安检中被动毫米波图像与可见光图像的融合方法 | |
Zhang et al. | A fusion algorithm for infrared and visible images based on saliency analysis and non-subsampled Shearlet transform | |
CN104700381B (zh) | 一种基于显著目标的红外与可见光图像融合方法 | |
Wang et al. | Low-illumination image enhancement for night-time UAV pedestrian detection | |
CN108921803B (zh) | 一种基于毫米波与可见光图像融合的去雾方法 | |
Luo et al. | Infrared and visible image fusion based on visibility enhancement and hybrid multiscale decomposition | |
Chu et al. | Novel fusion method for SAR and optical images based on non-subsampled shearlet transform | |
Wang et al. | Object-oriented change detection approach for high-resolution remote sensing images based on multiscale fusion | |
Kong et al. | Infrared and visible image fusion using structure-transferring fusion method | |
Liu et al. | Research on fusion technology based on low-light visible image and infrared image | |
Yao et al. | Optical remote sensing cloud detection based on random forest only using the visible light and near-infrared image bands | |
Jiang et al. | Automatic space debris extraction channel based on large field of view photoelectric detection system | |
Luo et al. | Infrared and visible image fusion based on visibility enhancement and norm optimization low-rank representation | |
Wu et al. | Fusing optical and synthetic aperture radar images based on shearlet transform to improve urban impervious surface extraction | |
Yue et al. | Low-illumination traffic object detection using the saliency region of infrared image masking on infrared-visible fusion image | |
Dey et al. | A robust FLIR target detection employing an auto-convergent pulse coupled neural network | |
Ducay et al. | Radiometric assessment of four pan-sharpening algorithms as applied to hyperspectral imagery | |
Greenberg et al. | Region-of-interest-based algorithm for automatic target detection in infrared images | |
Zhu et al. | YOLO-SDLUWD: YOLOv7-based small target detection network for infrared images in complex backgrounds | |
Romano et al. | Morphological operators for polarimetric anomaly detection | |
Liu et al. | EWST: an extreme weather scene text detector with dehazing and localization refinement | |
Li et al. | Algorithm for automatic image dodging of unmanned aerial vehicle images using two-dimensional radiometric spatial attributes | |
CN110567886B (zh) | 基于半监督空谱特征的多光谱云检测方法 | |
Lee et al. | Cloud detection of optical satellite images using support vector machine | |
Bayram et al. | Forward looking infrared imagery for landmine detection |
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 |