CN102096909A - 基于对数图像处理模型的改进的反锐化掩模图像增强方法 - Google Patents
基于对数图像处理模型的改进的反锐化掩模图像增强方法 Download PDFInfo
- Publication number
- CN102096909A CN102096909A CN 201010613607 CN201010613607A CN102096909A CN 102096909 A CN102096909 A CN 102096909A CN 201010613607 CN201010613607 CN 201010613607 CN 201010613607 A CN201010613607 A CN 201010613607A CN 102096909 A CN102096909 A CN 102096909A
- Authority
- CN
- China
- Prior art keywords
- image
- formula
- unsharp masking
- gray
- carried out
- 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
一种基于对数图像处理模型的改进的反锐化掩模图像增强方法,属于图像处理领域,方法为:将在传统模型下基于灰度的输入图像转换成对数图像处理模型下基于灰色调的图像;对图像进行改进的反锐化掩模增强:将反锐化掩模增强图像转换为传统模型下基于灰度的图像;本发明在对数图像处理模型下实现反锐化掩模图像增强方法,克服了传统模型下反锐化方法出现灰度“溢出”的不足;本发明利用对数图像处理模型的操作与传统处理模型的操作间的同构关系,通过同构变换和反同构变换,实现更简单、高效;本发明利用图像的梯度信息对图像的边缘区域和平滑区域实施不同程度增强,弥补了传统方法利用常数因子对整个图像进行相同强度增强的不足,所得结果更好。
Description
技术领域
本发明属于图像处理领域,特别是一种基于对数图像处理模型的改进的反锐化掩模图像增强方法。
背景技术
反锐化掩模法(Unsharp Masking,UM)是为了适应处理不同尺寸的图像特征而提出的图像边缘增强方法,该方法最早是应用于摄影技术中,以增强图像的边缘和细节;光学上的操作方法是将聚焦的正片和散焦的负片在底片上进行叠加,结果是增强了正片高频成份,从而增强了轮廓,散焦的负片相当于“模糊”模板(掩模),它与锐化的作用正好相反,因此,该方法被称为反锐化掩模法。
经典的线性反锐化掩模方法的原理如图1所示。首先将原图像低通滤波后产生一个钝化模糊图像,将原图像与这模糊图像相减得到保留高频成份的图像,再将高频图像用参数K放大后与原图像叠加,这就产生一个增强了边缘的图像,传统模型公式如下:
Ium=I+K·(I-Ius)
其中,I为原图像,Ius为低通滤波后图像,这里取均值滤波后的图像,K为正的常数;
尽管这种方法简单,增强效果也相对较好,但它却有2个突出的缺点:
(1)***对噪声非常敏感,由于采用了线性高通滤波器,图像的细节和噪声同时被增强,尤其是在图像的平坦区域,即使很小的噪声也非常明显;
(2)过冲(overshoot)现象,由于图像的高细节区域相对于其它区域增强更大,处理后的图像会呈现明显的人工处理痕迹;
在国内外公开的文献中,为了克服线性反锐化掩模方法的缺点,尤其是对噪声的敏感性,图像领域的研究者提出了各种各样的方法,这些方法大部分是用非线性滤波器代替线性高通滤波器,对减小噪声和增强细节进行折中考虑:S.K.Mitra提出了基于Teager算法的非线性算子,这种非线性算子可近似为局部均值权重高通滤波器;G.Ramponi提出了立方反锐化掩模方法,这种方法的实质是用一个对边缘敏感的平方滤波器算子乘以拉普拉斯算子,仅增强局部亮度变化区域的图像细节,相对较小噪声;Y.H.Lee提出了基于序列统计拉普拉斯方法的算子,这个算子的输出和局部均值与局部中值之差成比例,它能有效地去除高斯白噪声;A.Polesel提出了自适应反锐化掩模方法,这种方法采用自适应滤波器对图像的细节增强程度大一些,而对图像的平坦区域几乎不增强,因而能减小平坦区域的噪声;但是,上面提到的方法虽然相对于线性反锐化掩模方法减小了噪声,但在平坦区域,噪声仍然比较明显,而且,为了使图像中的细节区域达到较好的增强效果,图像的高细节区域往往增强过大,导致过冲现象的出现。
发明内容
为克服上述方法的缺陷,本发明提出一种基于对数图像处理模型的改进的反锐化掩模图像增强方法。
首先,由于经典的图像处理模型(包括传统的加操作和标量乘操作)并不能很好的符合图像叠加和放大的物理法则,而且对于灰度值在一定范围之内的数字图像,传统的加操作常常会产生溢出现象,因此为了改进经典模型的缺点与不足,Jourlin和Pinoli在上世纪八十年代中叶提出了对数图像处理模型:
在对数图像处理模型中,图像以映射的形式表示,称为灰色调函数(gray tone functions)并用字母f表示,对数图像处理模型公式如下:
本发明的技术方案是这样实现的:基于对数图像处理模型的改进的反锐化掩模图像增强方法,包括以下步骤:
步骤1:将在传统模型下基于灰度的输入图像转换成对数图像处理模型下基于灰色调的图像,公式为:
步骤2:对图像进行改进的反锐化掩模增强:利用梯度因子代替传统的常数增强因子,以区别对待平滑区域和边缘区域的增强系数,经步骤1处理后的图像,同时执行步骤2-1和步骤2-2,方法如下:
步骤2-1:计算图像的梯度因子g,首先利用高斯平滑算子对传统模型下基于灰度的图像进行滤波,公式为:
Is=I*Gσ
式中,Is为对传统模型下基于灰度的输入图像进行高斯滤波以后的图像,Gσ为方差为σ的高斯滤波器;
然后,利用差分方法计算图像的梯度信息Gray;
最后,计算图像的梯度因子g,公式为:
对图像进行反锐化掩模增强,公式为:
Ium=I+g×(I-Ius)
式中,Ius表示经低通滤波后的图像,公式为:
式中,N,M为图像的高和宽,I(i,j)表示原始图像,x,y分别表示当前所处理图像像素点的横纵坐标;
最后,利用反同构变换对上述图像进行处理,公式如下:
步骤2-2:对图像进行改进的反锐化掩模增强:首先,根据同构变换对灰色调图像进行处理,变换公式为:
式中,φ(I)表示对灰色调图像进行同构变换后的图像;
然后,利用均值滤波器对传统模型下基于灰度的输入图像进行低通滤波处理,公式为:
式中,N,M为图像的高和宽,I(i,j)表示原始图像;
之后,利用边缘信息作为权重因子对原图像进行反锐化掩模增强处理,此时的操作运算为传统的线性运算,其公式为:
Ium=I+g×(I-Ius)
最后,利用反同构变换对上述图像进行处理,公式如下:
步骤3:将步骤2得到的反锐化掩模增强图像转换为传统模型下基于灰度的图像,公式为:
本发明优点:本发明通过三个关键技术提高反锐化掩模图像增强方法的效果和效率,一是在对数图像处理模型下实现反锐化掩模图像增强方法,克服了传统模型下反锐化方法出现灰度“溢出”的不足,而且对数图像处理模型符合人体视觉感光原理,增强的结果视觉效果更好;二是利用对数图像处理模型的操作与传统处理模型的操作间的同构关系,通过同构变换和反同构变换,实现更简单、高效;三是利用图像的梯度信息定义边缘检测算子来改进反锐化增强的因子,对图像的边缘区域和平滑区域实施不同程度的增强,弥补了传统方法中利用常数因子对整个图像进行相同强度增强的不足,所得结果更好。
附图说明
图1为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法经典的线性反锐化掩模方法的原理图;
图2为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法总流程图;
图3为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法对图像进行改进的反锐化掩模增强流程图;
图4(a)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法原图像示意图;
图4(b)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法斑点噪声图像;
图4(c)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法高斯噪声图像;
图4(d)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法泊松噪声图像;
图4(e)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法传统算法对斑点噪声图像的增强结果示意图;
图4(f)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法传统算法对高斯噪声图像的增强结果示意图;
图4(g)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法传统算法对泊松噪声图像的增强结果示意图;
图4(h)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法改进算法对斑点噪声图像的增强结果示意图;
图4(i)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法改进算法对高斯噪声图像的增强结果示意图;
图4(j)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法改进算法对泊松噪声图像的增强结果示意图;
图5(a)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法医学原图像示意图;
图5(b)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法传统反锐化掩模增强结果示意图;
图5(c)为本发明基于对数图像处理模型的改进的反锐化掩模图像增强方法改进的反锐化掩模增强结果。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。
本实施例采用基于对数图像处理模型的改进的反锐化掩模图像增强方法,流程如图2、图3所示,包括以下步骤:
步骤1:将在传统模型下基于灰度的输入图像转换成对数图像处理模型下基于灰色调的图像,首先,输入一幅传统模型下基于灰度的图像,如图4(a)所示,该图像的一部分用矩阵形式表示为:
步骤2:对图像进行改进的反锐化掩模增强:利用梯度因子代替传统的常数增强因子,经步骤1处理后的图像,同时执行步骤2-1和步骤2-2,方法如下:
步骤2-1:计算图像的梯度因子g,首先利用高斯平滑算子对传统模型下基于灰度的图像进行滤波,公式为:
Is=I*Gσ
式中,Is为对传统模型下基于灰度的输入图像进行高斯滤波以后的图像,Gσ为方差为σ的高斯滤波器;
然后,利用差分方法计算图像的梯度信息Gray;
最后,计算图像的梯度因子g,公式为:
式中,取m=50,计算步骤1所述的传统模型下基于灰度的图像的梯度因子如下:
对图像进行反锐化掩模增强,公式为:
Ium=I+g×(I-Ius)
最后,利用反同构变换对上述图像进行处理,公式如下:
步骤2-2:对图像进行改进的反锐化掩模增强:首先,根据同构变换对灰色调图像进行处理,变换公式为:
式中,φ(I)表示对灰色调图像进行同构变换后的图像;
然后,利用均值滤波器对传统模型下基于灰度的输入图像进行低通滤波处理,公式为:
式中,N,M为图像的高和宽,I(i,j)表示是原始图像;
式中,传统模型下基于灰度的输入图像经过低通滤波后的矩阵为:
之后,利用边缘信息作为权重因子对原图像进行反锐化掩模增强处理,此时的操作运算为传统的线性运算,其公式为:
Ium=I+g×(I-Ius)
最后,利用反同构变换对上述图像进行处理,公式如下:
步骤3:将步骤2得到的反锐化掩模增强图像转换为传统模型下基于灰度的图像,其图像矩阵形式表示为:
图4(a)~4(j),其中图4(a)为原图像,图4(b)、图4(c)和图4(d)为加了噪声的图像(依次添加了斑点噪声、高斯噪声、泊松噪声),图4(e)、图4(f)和图4(g)为传统反锐化掩模算法的增强结果,由图中我们可以看出图像整体都比较亮,这是因为传统方法直接将图像的高频分量与原图像进行线性相加,这样做会造成灰度溢出现象,整幅图像比较亮;图4(h)、图4(i)和图4(j)为改进反锐化掩模算法的增强结果,由图中我们可以看到整幅图像对比很好,并没有出现普遍很亮的结果,尤其在边缘处明显得到了增强效果。
由于人类视觉对图像质量的评价很难用数学方法来表征,因此各种自适应图像增强方法对图像处理的性能也很难进行比较,通常,选择一些准则作为各种增强方法比较的基准,在对增强方法性能进行评估的过程中,本实施例将使用峰值信噪比(PSNR),计算公式为:
式中,N,M为图像的高和宽,I(i,j),Ium(i,j)分别是原始图像和经处理后图像的灰度;
表1为改进的反锐化掩模增强方法与传统反锐化模增强方法对不同噪声峰值信噪比对比表,由表1可知:经改进方法处理后的图像的峰值信噪比都比传统方法的大20%,说明本文方法优于传统方法。
表1.传统与改进反锐化增强算法对不同噪声的峰值信噪比对比表
图5为本发明的另一个应用:将改进的反锐化掩模方法应用于医学图像的处理结果:由图5(a)可知,原图像对比度低,噪声污染严重,骨骼轮廓不明显;由图5(b)可知,传统反锐化掩模增强图像后,对比度更低,整幅图像比较亮,骨骼轮廓没有清晰显示出来;由图5(c)可知,改进的反锐化掩模增强图像后,图像的对比度增大,骨骼的轮廓能清晰地显示出来。
Claims (2)
2.根据权利要求1所述的基于对数图像处理模型的改进的反锐化掩模图像增强方法,其特征在于:步骤2所述的对图像进行改进的反锐化掩模增强,经步骤1处理后的图像,同时执行步骤2-1和步骤2-2,方法如下:
步骤2-1:计算图像的梯度因子g,首先利用高斯平滑算子对传统模型下基于灰度的图像进行滤波,公式为:
Is=I*Gσ
式中,Is为对传统模型下基于灰度的输入图像进行高斯滤波以后的图像,Gσ为方差为σ的高斯滤波器;
然后,利用差分方法计算图像的梯度信息Gray;
最后,计算图像的梯度因子g,公式为:
对图像进行反锐化掩模增强,公式为:
Ium=I+g×(I-Ius)
式中,Ius表示经低通滤波后的图像,公式为:
式中,N,M为图像的高和宽,I(i,j)表示原始图像,x,y分别表示当前所处理图像像素点的横纵坐标;
最后,利用反同构变换对上述图像进行处理,公式如下:
式中,表示对灰色调图像进行反同构变换后的图像;
步骤2-2:对图像进行改进的反锐化掩模增强:首先,根据同构变换对灰色调图像进行处理,变换公式为:
式中,φ(I)表示对灰色调图像进行同构变换后的图像;
然后,利用均值滤波器对传统模型下基于灰度的输入图像进行低通滤波处理,公式为:
式中,N,M为图像的高和宽,I(i,j)表示原始图像;
之后,利用边缘信息作为权重因子对原图像进行反锐化掩模增强处理,此时的操作运算为传统的线性运算,其公式为:
Ium=I+g×(I-Ius)
最后,利用反同构变换对上述图像进行处理,公式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010106136079A CN102096909B (zh) | 2010-12-30 | 2010-12-30 | 基于对数图像处理模型的改进的反锐化掩模图像增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010106136079A CN102096909B (zh) | 2010-12-30 | 2010-12-30 | 基于对数图像处理模型的改进的反锐化掩模图像增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102096909A true CN102096909A (zh) | 2011-06-15 |
CN102096909B CN102096909B (zh) | 2012-07-04 |
Family
ID=44129987
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010106136079A Expired - Fee Related CN102096909B (zh) | 2010-12-30 | 2010-12-30 | 基于对数图像处理模型的改进的反锐化掩模图像增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102096909B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102970461A (zh) * | 2012-11-28 | 2013-03-13 | 信利光电(汕尾)有限公司 | 一种数据处理方法及*** |
CN106238347A (zh) * | 2016-07-27 | 2016-12-21 | 合肥高晶光电科技有限公司 | 一种ccd色选机的图像处理方法 |
CN106534693A (zh) * | 2016-11-25 | 2017-03-22 | 努比亚技术有限公司 | 一种照片处理方法、装置及终端 |
CN106911904A (zh) * | 2015-12-17 | 2017-06-30 | 通用电气公司 | 图像处理方法、图像处理***及成像*** |
CN108428215A (zh) * | 2017-02-15 | 2018-08-21 | 阿里巴巴集团控股有限公司 | 一种图像处理方法、装置和设备 |
CN108632502A (zh) * | 2017-03-17 | 2018-10-09 | 深圳开阳电子股份有限公司 | 一种图像锐化的方法及装置 |
CN108931883A (zh) * | 2017-05-17 | 2018-12-04 | 东方晶源微电子科技(北京)有限公司 | 一种优化掩模版图的方法 |
CN110070508A (zh) * | 2019-04-23 | 2019-07-30 | 西安交通大学 | 一种基于门限值和线性映射的反锐化增强方法 |
CN110580690A (zh) * | 2019-09-02 | 2019-12-17 | 杭州雄迈集成电路技术有限公司 | 一种识别峰值变换非线性曲线的图像增强方法 |
CN112116544A (zh) * | 2020-10-08 | 2020-12-22 | 东南数字经济发展研究院 | 一种抵抗图像平滑滤波的预处理方法 |
CN112150437A (zh) * | 2020-09-23 | 2020-12-29 | 南昌航空大学 | 激光增材制造扩压器裂纹缺陷dr检测图像处理方法 |
CN113042863A (zh) * | 2021-03-31 | 2021-06-29 | 山东齐星铁塔有限公司 | 一种基于激光视觉传感器的焊缝实时跟踪方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1120745A2 (en) * | 1999-12-30 | 2001-08-01 | General Electric Company | Extended dynamic range system for digital X-ray imaging detectors |
CN100563312C (zh) * | 2007-12-25 | 2009-11-25 | 青岛海信信芯科技有限公司 | 一种对比度增强方法 |
CN101878503A (zh) * | 2007-11-30 | 2010-11-03 | 夏普株式会社 | 基于经加权的误差向量的光源光选择的方法、***、程序和介质 |
-
2010
- 2010-12-30 CN CN2010106136079A patent/CN102096909B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1120745A2 (en) * | 1999-12-30 | 2001-08-01 | General Electric Company | Extended dynamic range system for digital X-ray imaging detectors |
CN101878503A (zh) * | 2007-11-30 | 2010-11-03 | 夏普株式会社 | 基于经加权的误差向量的光源光选择的方法、***、程序和介质 |
CN100563312C (zh) * | 2007-12-25 | 2009-11-25 | 青岛海信信芯科技有限公司 | 一种对比度增强方法 |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102970461A (zh) * | 2012-11-28 | 2013-03-13 | 信利光电(汕尾)有限公司 | 一种数据处理方法及*** |
CN106911904A (zh) * | 2015-12-17 | 2017-06-30 | 通用电气公司 | 图像处理方法、图像处理***及成像*** |
CN106911904B (zh) * | 2015-12-17 | 2020-04-21 | 通用电气公司 | 图像处理方法、图像处理***及成像*** |
CN106238347A (zh) * | 2016-07-27 | 2016-12-21 | 合肥高晶光电科技有限公司 | 一种ccd色选机的图像处理方法 |
CN106534693B (zh) * | 2016-11-25 | 2019-10-25 | 努比亚技术有限公司 | 一种照片处理方法、装置及终端 |
CN106534693A (zh) * | 2016-11-25 | 2017-03-22 | 努比亚技术有限公司 | 一种照片处理方法、装置及终端 |
CN108428215A (zh) * | 2017-02-15 | 2018-08-21 | 阿里巴巴集团控股有限公司 | 一种图像处理方法、装置和设备 |
CN108632502A (zh) * | 2017-03-17 | 2018-10-09 | 深圳开阳电子股份有限公司 | 一种图像锐化的方法及装置 |
CN108931883A (zh) * | 2017-05-17 | 2018-12-04 | 东方晶源微电子科技(北京)有限公司 | 一种优化掩模版图的方法 |
CN108931883B (zh) * | 2017-05-17 | 2022-06-21 | 东方晶源微电子科技(北京)有限公司 | 一种优化掩模版图的方法 |
CN110070508A (zh) * | 2019-04-23 | 2019-07-30 | 西安交通大学 | 一种基于门限值和线性映射的反锐化增强方法 |
CN110580690A (zh) * | 2019-09-02 | 2019-12-17 | 杭州雄迈集成电路技术有限公司 | 一种识别峰值变换非线性曲线的图像增强方法 |
CN112150437A (zh) * | 2020-09-23 | 2020-12-29 | 南昌航空大学 | 激光增材制造扩压器裂纹缺陷dr检测图像处理方法 |
CN112150437B (zh) * | 2020-09-23 | 2023-05-09 | 南昌航空大学 | 激光增材制造扩压器裂纹缺陷dr检测图像处理方法 |
CN112116544A (zh) * | 2020-10-08 | 2020-12-22 | 东南数字经济发展研究院 | 一种抵抗图像平滑滤波的预处理方法 |
CN113042863A (zh) * | 2021-03-31 | 2021-06-29 | 山东齐星铁塔有限公司 | 一种基于激光视觉传感器的焊缝实时跟踪方法 |
CN113042863B (zh) * | 2021-03-31 | 2023-01-24 | 山东齐星铁塔有限公司 | 一种基于激光视觉传感器的焊缝实时跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102096909B (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102096909B (zh) | 基于对数图像处理模型的改进的反锐化掩模图像增强方法 | |
CN100568279C (zh) | 一种基于Retinex理论的快速彩色图像增强方法 | |
CN100478993C (zh) | 根据边缘和亮度特征进行自适应调整的图像增强方法 | |
CN104574293B (zh) | 基于有界运算的多尺度Retinex图像清晰化算法 | |
CN103295191B (zh) | 多尺度视觉自适应图像增强方法及评价方法 | |
CN103268598B (zh) | 基于Retinex理论的低照度低空遥感影像增强方法 | |
CN107527333A (zh) | 一种基于伽马变换的快速图像增强方法 | |
CN102289792A (zh) | 一种低照度视频图像增强方法及*** | |
CN104182947A (zh) | 一种低照度图像增强方法和*** | |
CN103530848A (zh) | 一种非均匀光照图像二次曝光的实现方法 | |
CN104240194A (zh) | 一种基于抛物线函数的低照度图像增强算法 | |
CN103295206B (zh) | 一种基于Retinex的微光图像增强方法和装置 | |
CN102222323A (zh) | 基于直方图统计拉伸和梯度滤波的红外图像细节增强方法 | |
CN106127688A (zh) | 一种超分辨率图像重建方法及其*** | |
CN105139366A (zh) | 一种基于空间域的图像增强方法 | |
CN102306378B (zh) | 一种图像增强方法 | |
CN102930517A (zh) | 直方图均衡化图像增强方法 | |
CN107292834A (zh) | 红外图像细节增强方法 | |
CN103106644B (zh) | 克服彩色图像非均匀光照的自适应画质增强方法 | |
Vijayakumar | Selective image enhancement and restoration for skin cancer identification | |
CN102999890B (zh) | 基于环境因素的图像光线强度动态校正方法 | |
CN103700077B (zh) | 一种基于人眼视觉特性的自适应图像增强方法 | |
CN105427255A (zh) | 一种基于grhp的无人机红外图像细节增强方法 | |
CN106981052A (zh) | 基于变分框架的自适应亮度不均变分校正方法 | |
CN101478689A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120704 Termination date: 20211230 |
|
CF01 | Termination of patent right due to non-payment of annual fee |