CN103208105A - 一种红外图像细节增强与噪声自适应抑制方法 - Google Patents
一种红外图像细节增强与噪声自适应抑制方法 Download PDFInfo
- Publication number
- CN103208105A CN103208105A CN201310158402XA CN201310158402A CN103208105A CN 103208105 A CN103208105 A CN 103208105A CN 201310158402X A CN201310158402X A CN 201310158402XA CN 201310158402 A CN201310158402 A CN 201310158402A CN 103208105 A CN103208105 A CN 103208105A
- Authority
- CN
- China
- Prior art keywords
- infrared image
- image
- coe
- strengthens
- coordinate
- 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.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种红外图像细节增强与噪声自适应抑制方法,包括以下步骤:(1)获取原始红外图像Forg;(2)使用微分算子计算原始红外图像的滤波分量Fd;(3)计算原始红外图像背景复杂度,确定噪声抑制阈值T1;(4)计算原始红外图像的邻域方差V,根据初始阈值T2将邻域方差图像二值化,用矩阵BW表示二值化图像;(5)计算滤波分量Fd的加权矩阵Coe;(6)输出细节增强图像Fout。本发明与现有技术相比,具有以下显著优点:(1)红外图像细节增强与噪声自适应抑制新方法特别针对信噪比的图像进行处理,在图像细节增强的同时有效抑制了图像的噪声放大。
Description
技术领域
本发明涉及红外图像处理领域,特别涉及一种适合硬件实时实现的图像细节增强与噪声自适应抑制方法。
背景技术
在红外图像摄取、传输和处理过程中,受光的衍射、聚焦不良、景物和取像装置的相对运动的影响,不可避免导致图像质量下降,使图像变模糊,此外,电子***高频性能不好也会损失图像的高频分量,而使图像不清晰。大量研究表明,图像变模糊物理过程的数学模型一般含有求和、平均或积分运算。根据图像变模糊的过程可知,在空间域中运用微分运算或在频率域中用加强信号高频分量的方法可有效增强图像细节。
目前,图像细节增强技术受到了研究者的广泛关注。针对红外成像的特点,研究提出了多种有效的红外图像细节增强处理算法,通过增强场景中目标与背景之间的灰度对比度以及图像细节信息,如边缘、轮廓、纹理等,可提高对弱小目标探测、跟踪及识别的能力,增强对图像内容和关键细节信息理解的准确性。现有的图像细节增强算法存在以下缺点:(1)多数现有的图像细节增强算法仅适用于高信噪比的原始红外图像,随着原始红外图像信噪比下降,增强后的图像视觉效果变差,不利用人或机器对图像的理解和分析;(2)多数现有的噪声抑制算法对整幅图像采用同一滤波参数进行噪声抑制,没有利用图像局部特征信息,导致噪声抑制的同时损失了大量图像细节信息;(3)多数现有的图像细节增强与噪声自适应抑制算法运算量大,不易硬件实时实现。
发明内容
发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供一种红外图像细节增强与噪声自适应抑制方法。
为了解决上述技术问题,本发明公开了一种红外图像细节增强与噪声自适应抑制方法,包括以下步骤:
(1)获取原始红外图像Forg;
(2)运用微分算子计算原始红外图像的滤波分量Fd;
(3)计算原始红外图像背景复杂度,确定噪声抑制阈值T1;
(4)计算原始红外图像的邻域方差V,根据初始阈值T2将邻域方差图像二值化,用矩阵BW表示;
(5)计算滤波分量Fd的加权矩阵Coe;
(6)输出细节增强图像Fout;
本发明中使用无方向性的拉氏算子计算滤波分量Fd。
本发明红外图像细节增强与噪声自适应抑制新方法中,背景复杂度用邻域信息熵特征H表示。红外图像中以坐标(i,j)为中心的P×Q大小的邻域信息熵H(i,j)计算公式为:
其中,Forg(i,j)为原始红外图像中坐标(i,j)位置的灰度值,H(i,j)和p(i,j)分别为坐标(i,j)位置邻域信息熵H和灰度分布概率的输出值,0≤i≤M,0≤j≤N,M和N分别表示原始红外图像的行数和列数。P和Q分别表示邻域的高度和宽度,5≤P≤15,5≤Q≤15。
本发明红外图像细节增强与噪声自适应抑制新方法中,根据原始红外图像平均背景复杂度确定噪声抑制阈值T1,背景复杂度高,阈值T1小,相反背景复杂度低,阈值T1大。
坐标(i,j)位置的噪声抑制阈值T1(i,j)计算公式如下:
其中,T1(i,j)是抑制阈值T1中坐标(i,j)位置的噪声抑制阈值,A和B是常系数。A和B的取值范围为自然数,A取值越大,则噪声抑制能力越强;B取值越大,则噪声抑制能力越强。
本发明红外图像细节增强与噪声自适应抑制新方法中,图像的邻域方差V通过邻域周边像素均值减去中心像素均值方法进行计算。
本发明红外图像细节增强与噪声自适应抑制新方法中,加权系数Coe(i,j)是通过统计二值化图像BW中以坐标(i,j)为中心大小P1×Q1窗口非零值数量,若非零值数量大于噪声抑制阈值T1(i,j),Coe(i,j)=1,否则Coe(i,j)=0,其中,Coe(i,j)是加权矩阵Coe中坐标(i,j)位置的输出值,T1(i,j)是噪声抑制阈值T1中(i,j)位置的输出值,P1和Q1分别是窗口的高度和宽度。P1和Q1取值范围为自然数,3≤P1≤7,3≤Q1≤7。
本发明红外图像细节增强与噪声自适应抑制新方法中,坐标(i,j)位置Fout(i,j)表达式如下所示:
Fout(i,j)=Forg(i,j)-Coe(i,j)×Fd(i,j)
其中,Fout(i,j)、Forg(i,j)和Fd(i,j)分别是坐标(i,j)位置的细节增强图像Fout、原始红外图像Forg(i,j)和滤波分量Fd的灰度值。
本发明与现有技术相比,具有以下显著优点:(1)红外图像细节增强与噪声自适应抑制新方法特别针对信噪比的图像进行处理,在图像细节增强的同时有效抑制了图像的噪声放大;(2)利用背景复杂度局部特征计算滤波系数,可使算法具备对同一图像多尺度抑噪能力,实现平滑区域抑噪作用强,细节区域抑噪作用弱,图像细节有效保留;(3)利用概率统计原理,算法认为区域内图像细节表现为多个高方差点的集合,而噪声则表现为孤立高方差点,它是区分图像细节和噪声的有效手段;(4)本发明提出的红外图像细节增强与噪声自适应抑制新方法不存在高阶运算和复杂结构,算法运算量小,易于硬件实时实现。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
图1为本发明流程图。
图2a为原始图像。
图2b为传统图像细节增强图。
图2c为实施例处理结果图。
具体实施方式
本发明公开了一种红外图像细节增强与噪声自适应抑制方法,包括以下步骤:
(1)获取原始红外图像Forg;
(2)使用微分算子计算原始红外图像的滤波分量Fd;
(3)计算原始红外图像背景复杂度,确定噪声抑制阈值T1;
(4)计算原始红外图像的邻域方差V,根据初始阈值T2将邻域方差图像二值化,用矩阵BW表示二值化图像;
(5)计算滤波分量Fd的加权矩阵Coe;
(6)输出细节增强图像Fout。
使用无方向性的拉氏算子计算滤波分量Fd。
背景复杂度用邻域信息熵特征H表示,红外图像中以坐标(i,j)为中心的P×Q大小的邻域信息熵H(i,j)计算公式为:
其中,Forg(i,j)为原始红外图像中坐标(i,j)位置的灰度值,H(i,j)和p(i,j)分别为坐标(i,j)位置邻域信息熵和灰度分布概率的值,0≤i≤M,0≤j≤N,M和N分别表示原始红外图像的行数和列数,P和Q分别表示邻域的高度和宽度,5≤P≤15,5≤Q≤15。
坐标(i,j)位置的噪声抑制阈值T1(i,j)计算公式如下:
其中,T1(i,j)是抑制阈值T1中坐标(i,j)位置的噪声抑制阈值,A和B是常系数。
步骤(4)中,图像的邻域方差V通过邻域周边像素均值减去中心像素均值方法进行计算。
步骤(5)中,加权矩阵Coe的加权系数Coe(i,j)通过统计矩阵BW中以坐标(i,j)为中心大小P1×Q1窗口非零值数量,若非零值数量大于噪声抑制阈值T1(i,j),令Coe(i,j)=1,否则令Coe(i,j)=0,其中,Coe(i,j)是加权矩阵Coe中(i,j)位置的加权系数,P1和Q1分别是窗口的高度和宽度。
步骤(6)中,坐标(i,j)位置细节增强图像Fout(i,j)采用下式计算:
Fout(i,j)=Forg(i,j)-Coe(i,j)×Fd(i,j),
其中,Fd(i,j)是滤波分量Fd中坐标(i,j)位置的灰度值。
实施例1
结合图1,下面以实例来说明本发明红外图像细节增强与噪声自适应抑制新方法。红外图像的像素个数320×256,帧频25HZ。红外热像仪的数字信号通过光纤传给DSP+FPGA架构的专用图像处理板,红外图像细节增强与噪声自适应抑制新方法在DSP处理器中实现,满足实时处理的需求,具体实施步骤如下:
(1)获取原始红外图像Forg;
原始红外图像Forg是14-bit数字图像。红外热像仪的数字信号通过光纤打包方式传给专用图像处理板的FPGA,FPGA完成光纤解包将数字信号经DSP的EMIF接口传给DSP,得到14-bit原始红外图像Forg。
(2)运用微分算子计算原始红外图像的滤波分量Fd;
由于图像细节边界可能是任意走向的,因此需要微分算子具备无方向性。这里,选用无方向性的8邻域拉氏算子作为微分算子。对离散图像Forg,滤波分量Fd(i,j)表达式如下所示:
(3)计算原始红外图像背景复杂度,确定噪声抑制阈值T1;
背景复杂度用邻域信息熵特征表示。这里,选用10×10计算邻域信息熵,表达式如下所示:
其中,Forg(i,j)为图像中坐标(i,j)位置的灰度值,p(i,j)为坐标(i,j)位置的灰度分布概率,0≤i≤256,0≤j≤320。
坐标(i,j)位置的噪声抑制阈值T1(i,j)表达式如下所示:
其中,T1(i,j)是抑制阈值T1中(i,j)位置的输出值,A和B是常系数,A=3,B=6,抑制阈值T1(i,j)增大则噪声抑制能力变强。
(4)计算原始红外图像的邻域方差V,根据初始阈值T2将邻域方差图像二值化,用矩阵BW表示;
这里,选用5×5邻域计算邻域方差V,坐标(i,j)位置的邻域方差V(i,j)表达式如下所示:
其中,V(i,j)是邻域方差V中(i,j)位置的输出值,Sum3×3和Sum5×5分别是原始红外图像Forg中以(i,j)为中心3×3和5×5邻域像素值总和。
因此,(i,j)位置的二值化矩阵BW(i,j)表达式如下所示:
其中,BW(i,j)是BW中(i,j)位置的输出值,T2是常系数,初始化值T2=8。
(5)计算滤波分量Fd的加权矩阵Coe;
加权系数Coe(i,j)是通过统计BW中以坐标(i,j)为中心P1×Q1窗口中数值1的数量,P1=5,Q1=5,若数值1的数量大于噪声抑制阈值T1(i,j),Coe(i,j)=1,否则Coe(i,j)=0,其中,Coe(i,j)是加权矩阵Coe中(i,j)位置的输出值,T1(i,j)是噪声抑制阈值T1中(i,j)位置的输出值。
(6)输出细节增强图像Fout;
坐标(i,j)位置Fout(i,j)表达式如下所示:
Fout(i,j)=Forg(i,j)-Coe(i,j)×Fd(i,j)。
其中,Fout(i,j)、Forg(i,j)和Fd(i,j)分别是增强图像Fout、原始红外图像Forg(i,j)和滤波分量Fd中(i,j)位置的灰度值,0≤i≤256,0≤j≤320。
实施例2
图2中,图2a表示原始图像,图2b表示传统图像细节增强,图2c是图像细节增强与噪声自适应抑制方法的处理结果。原始图像的细节模糊;传统图像细节增强算法增强了图像细节,但图像噪声被放大;图像细节增强与噪声自适应抑制方法在增强图像细节的同时,有效抑制了点状噪声。
本发明提供了一种红外图像细节增强与噪声自适应抑制方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (7)
1.一种红外图像细节增强与噪声自适应抑制方法,其特征在于,包括以下步骤:
(1)获取原始红外图像Forg;
(2)使用微分算子计算原始红外图像的滤波分量Fd;
(3)计算原始红外图像背景复杂度,确定噪声抑制阈值T1;
(4)计算原始红外图像的邻域方差V,根据初始阈值T2将邻域方差图像二值化,用矩阵BW表示二值化图像;
(5)计算滤波分量Fd的加权矩阵Coe;
(6)输出细节增强图像Fout。
2.根据权利要求1所述的一种红外图像细节增强与噪声自适应抑制方法,其特征在于,使用无方向性的拉氏算子计算滤波分量Fd。
3.根据权利要求1所述的一种红外图像细节增强与噪声自适应抑制方法,其特征在于,背景复杂度用邻域信息熵特征H表示,红外图像中以坐标(i,j)为中心的P×Q大小的邻域信息熵H(i,j)计算公式为:
其中,Forg(i,j)为原始红外图像中坐标(i,j)位置的灰度值,H(i,j)和p(i,j)分别为坐标(i,j)位置邻域信息熵和灰度分布概率的值,0≤i≤M,0≤j≤N,M和N分别表示原始红外图像的行数和列数,P和Q分别表示邻域的高度和宽度,5≤P≤15,5≤Q≤15。
5.根据权利要求4所述的一种红外图像细节增强与噪声自适应抑制方法,其特征在于,步骤(4)中,图像的邻域方差V通过邻域周边像素均值减去中心像素均值方法进行计算。
6.根据权利要求5所述的一种红外图像细节增强与噪声自适应抑制方法,其特征在于,步骤(5)中,加权矩阵Coe的加权系数Coe(i,j)通过统计矩阵BW中以坐标(i,j)为中心大小P1×Q1窗口非零值数量,若非零值数量大于噪声抑制阈值T1(i,j),令Coe(i,j)=1,否则令Coe(i,j)=0,其中,Coe(i,j)是加权矩阵Coe中(i,j)位置的加权系数,P1和Q1分别是窗口的高度和宽度。
7.根据权利要求6所述的一种红外图像细节增强与噪声自适应抑制方法,其特征在于,步骤(6)中,坐标(i,j)位置细节增强图像Fout(i,j)采用下式计算:
Fout(i,j)=Forg(i,j)-Coe(i,j)×Fd(i,j),
其中,Fd(i,j)是滤波分量Fd中坐标(i,j)位置的灰度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310158402.XA CN103208105B (zh) | 2013-05-02 | 2013-05-02 | 一种红外图像细节增强与噪声自适应抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310158402.XA CN103208105B (zh) | 2013-05-02 | 2013-05-02 | 一种红外图像细节增强与噪声自适应抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103208105A true CN103208105A (zh) | 2013-07-17 |
CN103208105B CN103208105B (zh) | 2015-08-19 |
Family
ID=48755322
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310158402.XA Active CN103208105B (zh) | 2013-05-02 | 2013-05-02 | 一种红外图像细节增强与噪声自适应抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103208105B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558754A (zh) * | 2013-10-29 | 2014-02-05 | 国家电网公司 | 一种水电站机组进水控制方法 |
CN104580826A (zh) * | 2015-02-03 | 2015-04-29 | 成都金本华科技股份有限公司 | 一种增强图像锐度的视频信号处理方法 |
CN106228525A (zh) * | 2016-08-25 | 2016-12-14 | 成都市晶林科技有限公司 | 一种抑制红外图像过度增强的方法 |
CN107292834A (zh) * | 2017-05-24 | 2017-10-24 | 杭州天铂红外光电技术有限公司 | 红外图像细节增强方法 |
CN107330863A (zh) * | 2017-05-27 | 2017-11-07 | 浙江大学 | 一种基于噪声估计的图像去噪方法 |
CN107330864A (zh) * | 2017-06-05 | 2017-11-07 | 中国电子科技集团公司第二十八研究所 | 一种基于改进局部纹理特征的红外图像处理方法 |
CN110647854A (zh) * | 2019-09-27 | 2020-01-03 | 华清永安(北京)科技发展有限责任公司 | 一种垃圾分类排放智能管理*** |
CN111652808A (zh) * | 2019-07-29 | 2020-09-11 | 天津大学 | 一种红外图像细节增强及噪声抑制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5818964A (en) * | 1994-12-27 | 1998-10-06 | Texas Instruments Incorporated | Method and apparatus for selecting an adaptive filter for image data |
CN101620727A (zh) * | 2009-08-10 | 2010-01-06 | 电子科技大学 | 红外图像加权直方图自适应增强算法 |
CN101359402B (zh) * | 2008-09-11 | 2010-12-15 | 北京航空航天大学 | 复杂背景下红外背景预测算子的自动选择方法 |
CN102289819A (zh) * | 2011-07-27 | 2011-12-21 | 中国电子科技集团公司第二十八研究所 | 一种背景自适应估计的红外运动目标实时检测方法 |
CN102324021A (zh) * | 2011-09-05 | 2012-01-18 | 电子科技大学 | 一种基于剪切波变换的红外弱小目标检测方法 |
-
2013
- 2013-05-02 CN CN201310158402.XA patent/CN103208105B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5818964A (en) * | 1994-12-27 | 1998-10-06 | Texas Instruments Incorporated | Method and apparatus for selecting an adaptive filter for image data |
CN101359402B (zh) * | 2008-09-11 | 2010-12-15 | 北京航空航天大学 | 复杂背景下红外背景预测算子的自动选择方法 |
CN101620727A (zh) * | 2009-08-10 | 2010-01-06 | 电子科技大学 | 红外图像加权直方图自适应增强算法 |
CN102289819A (zh) * | 2011-07-27 | 2011-12-21 | 中国电子科技集团公司第二十八研究所 | 一种背景自适应估计的红外运动目标实时检测方法 |
CN102324021A (zh) * | 2011-09-05 | 2012-01-18 | 电子科技大学 | 一种基于剪切波变换的红外弱小目标检测方法 |
Non-Patent Citations (2)
Title |
---|
JIAN YANG,ETAL: "Image Texture Feature Extraction Method Based on Regional Average Binary Gray Level Difference Co-occurrence Matrix", 《VIRTUAL REALITY AND VISUALIZATION (ICVRV)》 * |
白俊奇 等: "红外图像噪声滤波对比度增强算法", 《红外与激光工程》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558754B (zh) * | 2013-10-29 | 2016-12-07 | 国家电网公司 | 一种水电站机组进水控制方法 |
CN103558754A (zh) * | 2013-10-29 | 2014-02-05 | 国家电网公司 | 一种水电站机组进水控制方法 |
CN104580826A (zh) * | 2015-02-03 | 2015-04-29 | 成都金本华科技股份有限公司 | 一种增强图像锐度的视频信号处理方法 |
CN106228525A (zh) * | 2016-08-25 | 2016-12-14 | 成都市晶林科技有限公司 | 一种抑制红外图像过度增强的方法 |
CN106228525B (zh) * | 2016-08-25 | 2020-08-21 | 成都市晶林科技有限公司 | 一种抑制红外图像过度增强的方法 |
CN107292834A (zh) * | 2017-05-24 | 2017-10-24 | 杭州天铂红外光电技术有限公司 | 红外图像细节增强方法 |
CN107330863B (zh) * | 2017-05-27 | 2019-10-18 | 浙江大学 | 一种基于噪声估计的图像去噪方法 |
CN107330863A (zh) * | 2017-05-27 | 2017-11-07 | 浙江大学 | 一种基于噪声估计的图像去噪方法 |
CN107330864A (zh) * | 2017-06-05 | 2017-11-07 | 中国电子科技集团公司第二十八研究所 | 一种基于改进局部纹理特征的红外图像处理方法 |
CN107330864B (zh) * | 2017-06-05 | 2019-08-02 | 中国电子科技集团公司第二十八研究所 | 一种基于改进局部纹理特征的红外图像处理方法 |
CN111652808A (zh) * | 2019-07-29 | 2020-09-11 | 天津大学 | 一种红外图像细节增强及噪声抑制方法 |
CN110647854A (zh) * | 2019-09-27 | 2020-01-03 | 华清永安(北京)科技发展有限责任公司 | 一种垃圾分类排放智能管理*** |
CN110647854B (zh) * | 2019-09-27 | 2020-07-28 | 华清永安(北京)科技发展有限责任公司 | 一种垃圾分类排放智能管理*** |
Also Published As
Publication number | Publication date |
---|---|
CN103208105B (zh) | 2015-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103208105B (zh) | 一种红外图像细节增强与噪声自适应抑制方法 | |
CN104657945B (zh) | 复杂背景下多尺度时空联合滤波的红外小目标检测方法 | |
KR101298642B1 (ko) | 영상 잡음 제거 방법 및 장치 | |
CN103020628B (zh) | 一种基于rgb对比图像与目标形状的烟雾检测方法 | |
CN100474337C (zh) | 一种基于径向基神经网络的有噪运动模糊图像复原方法 | |
CN102831591B (zh) | 一种基于高斯滤波的单幅图像的实时去雾方法 | |
Singh et al. | Image dehazing using Moore neighborhood-based gradient profile prior | |
CN103020965B (zh) | 一种基于显著性检测的前景分割方法 | |
CN101489034B (zh) | 一种视频图像噪声估计与去除方法 | |
TWI575503B (zh) | 影像處理裝置及影像處理方法 | |
CN101860667A (zh) | 一种快速去除图像中混合噪声的方法 | |
CN102063704B (zh) | 一种机载视景增强方法及其装置 | |
CN103426148A (zh) | 生成低分辨率输入数据结构的超分辨率版本的方法和设备 | |
CN103295225B (zh) | 弱光照条件下的列车转向架边缘检测方法 | |
CN104657947B (zh) | 针对碱基图像的一种降噪方法 | |
CN104574293A (zh) | 基于有界运算的多尺度Retinex图像清晰化算法 | |
CN103985140B (zh) | 一种基于频域处理的抗噪运动目标检测方法 | |
CN107784639B (zh) | 一种无人机遥感图像改进多边滤波去噪方法 | |
CN110135312B (zh) | 一种基于分级lcm的快速小目标检测方法 | |
CN104036461B (zh) | 一种基于联合滤波的红外复杂背景抑制方法 | |
CN105719251B (zh) | 一种用于大像移线性模糊的压缩降质图像复原方法 | |
CN100367770C (zh) | 一种去除视频孤立噪声点的方法 | |
CN105427255A (zh) | 一种基于grhp的无人机红外图像细节增强方法 | |
CN110351453A (zh) | 一种计算机视频数据处理方法 | |
CN105809633A (zh) | 去除颜色噪声的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |