CN101046847A - 一种基于二次多项式光照模型的人脸光照对齐方法 - Google Patents
一种基于二次多项式光照模型的人脸光照对齐方法 Download PDFInfo
- Publication number
- CN101046847A CN101046847A CN 200710027817 CN200710027817A CN101046847A CN 101046847 A CN101046847 A CN 101046847A CN 200710027817 CN200710027817 CN 200710027817 CN 200710027817 A CN200710027817 A CN 200710027817A CN 101046847 A CN101046847 A CN 101046847A
- Authority
- CN
- China
- Prior art keywords
- illumination
- human face
- model
- component
- quadratic polynomial
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005286 illumination Methods 0.000 claims abstract description 106
- 238000012549 training Methods 0.000 claims abstract description 19
- 230000001815 facial effect Effects 0.000 claims description 22
- 238000010606 normalization Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 abstract 1
- 238000000513 principal component analysis Methods 0.000 description 8
- 230000000007 visual effect Effects 0.000 description 8
- 230000006872 improvement Effects 0.000 description 5
- 101100042371 Caenorhabditis elegans set-3 gene Proteins 0.000 description 4
- 101150104646 SET4 gene Proteins 0.000 description 4
- 101150117538 Set2 gene Proteins 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 210000001938 protoplast Anatomy 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于二次多项式光照模型的人脸光照对齐方法,包括以下步骤:(1)建立光照模型;(2)对所有训练图像以及目标图像均进行形状规一化;(3)训练系数矩阵Ai和Bi,其中i=2,…,64;(4)对目标人脸图像进行分解;(5)对目标人脸图像进行光照类型估计;(6)对目标人脸图像的光照分量进行对齐矫正;(7)目标人脸图像重构。本发明采用二次多项式描述了光照变化引起的图像变化,从而对基于LTV模型的光照处理算法中丢弃的人脸图像光照分量进行了对齐矫正再次利用,在获得更高的人脸识别率同时达到人脸视觉效果的改善。
Description
技术领域
本发明涉及一种人脸光照对齐方法,尤其涉及一种基于二次多项式光照模型的人脸光照对齐方法。
背景技术
人脸技术在公共安全系统、身份鉴别、虚拟游戏等方面均有广泛应用,但是光照问题却是长期困扰该项技术实用化的主要原因之一。光照不均匀不仅影响视觉效果,而且严重影响人脸识别率。近二十年来,各种各样的光照处理技术被提出来用于人脸技术,但是其中大多数尚达不到实用要求或者不能适应多方面的应用要求。
Terrence Chen等人于2006年提出基于LTV模型的光照处理算法,即在对数域用全变分模型提取人脸图像的反射分量用于人脸识别。由于反射分量对光照变化不敏感,所以基于LTV模型的光照处理算法能取得很高的人脸识别结果。但是实际上,用LTV模型分解得到的光照分量中仍然含有大量对人脸识别有用的信息,这些成分在基于LTV模型的光照处理算法中并没有用于人脸识别。另外该方法仅仅是提取光照不变量用于人脸识别,并没有达到人脸可视效果上的真正改善,这也限制了该方法的应用范围。
发明内容
本发明的目的在于克服现有技术的不足,提供一种实用性强、应用范围广、可明显提高人脸识别率的基于二次多项式光照模型的人脸光照对齐方法。本方法对基于LTV模型的光照处理算法中丢弃的光照成分进行了对齐矫正后利用,在获得更高的人脸识别率同时达到视觉效果的改善。
本发明的目的通过下述技术方案实现:一种基于二次多项式光照模型的人脸光照对齐方法,包括下述步骤:
(1)建立光照模型。根据Yale B人脸库的采集条件,将人脸光照条件划分为64种,其中第1种为正面标准光照,对标准光照条件下的光照分量Lstandrad和任意一种光照条件下的光照分量Lobject,可以建立以下模型:
(2)对所有训练图像以及目标图像均进行形状规一化处理。即对每张图像,先检测若干个特征点的坐标,通过旋转使得每张人脸的两只眼睛处于水平位置,再运用双插值算法拉伸图像,使得这些特征点位于图像的固定位置,最后将图像裁剪为同一大小;
(3)训练系数矩阵Ai和Bi,i=2Λ,64。用Yale B人脸库中的正面人脸采用最小二乘法对模型①进行训练,其中每张人脸图像的光照分量L用LTV模型分解获得(其分解过程见步骤(4))。第i种光照条件对应的系数矩阵Ai和Bi求解模型如下:
其中m为训练样本组数,L为标准光照条件下的人脸光照分量,L′为某特定光照条件下的人脸光照分量。a,b标记标准光照条件下光照分量的范围,即a≤L(x,y)≤b,c为一个大于0的数。另外
L′(x,y)为训练样本平均值,即
(4)对目标人脸图像进行分解。对每张待进行光照对齐的人脸图像Iobject,根据光照—反射模型Iobject(x,y)=Lobject(x,y)Robject(x,y),用LTV模型对其进行分解可以得到光照分量Lobject和反射分量Robject。分解过程如下:
f(x,y)=logIobject(x,y)=log Robject(x,y)+log Lobject(x,y) ④,
求解以下TV模型:
便可得到Lobject≈exp(u),Robject≈exp(f-u) ⑥,
实际上步骤(3)中也是用同样的方法对所有训练样本进行分解。
(5)对目标人脸图像进行光照类型估计。对步骤(4)中分解得到的光照分量Lobject,直接用主成分分析(principal component analysis,PCA)识别方法估计出Iobject的光照类型。即Iobject的光照条件与Yale B人脸库中的64种光照条件中的第i种最接近时,便记其光照类型为i,若光照类型为i=1,则不用进行光照对齐。
(6)对目标人脸图像的光照分量Lobject进行对齐矫正。根据步骤(5)中估计得到的光照类型i相应选取步骤(3)中训练得到的系数矩阵Ai和Bi,用以下式子对Lobject进行对齐矫正:
其中如果i=1则不必再进行调整。
(7)目标人脸图像重构。用步骤(4)中分解得到的反射分量Robject和步骤(6)中光照对齐后得到的光照分量Laligned重构最终图像:
Ialigned(x,y)=Laligned(x,y)Robject(x,y) ⑧。
本发明与现有技术相比具有如下优点和有益效果:
1、本发明提出的省掉一次项的二次多项式模型能很好描述光照变化引起人脸图像的光照分量变化,用以进行光照对齐的结果真正达到了视觉效果的改善。
2、本发明方法对进行光照对齐的人脸图像不必进行严格的形状对齐,实际上,只需选取三个特征点进行简单的人脸形状对齐便可达到很好的光照对齐效果。而且提出的模型形式简单,非常容易实现。与基于LTV模型的光照处理算法相比,本发明方法有较大的改进,对用LTV模型提取的光照不变量Robject,本发明在整个算法过程中并没有对其进行任何不必要的调整。
3、本发明利用到光照分量中对人脸识别有用的信息,在基于LTV模型的光照处理算法的基础上再次提高识别率。最重要的是,利用到光照分量,本方法进行光照对齐的结果真正达到了视觉效果的改善。
附图说明
图1为本发明方法的操作流程图。
图2为本发明的2D训练图像以及目标图像形状规一化状态图;左为原人脸图;中为裁剪为112×92象素的人脸图;右为裁剪为112×100象素的人脸图。
图3为本发明的LTV模型分解结果图;左为原图;中为光照分量图;右为反射分量图。
图4为本发明的光照对齐结果图;左为原光照分量图;右为对齐后的光照分量图。
图5为本发明的光照对齐整体流程图。
图6为视觉效果恢复比较图。第一行为原图;第二行为基于LTV模型的光照处理算法的结果图;最后一行为本发明的光照对齐结果图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例
图1示出了本发明的操作过程,由图1可见,本基于二次多项式光照模型的人脸光照对齐方法,包括下述步骤:
(1)对所有训练图像以及目标图像均进行形状规一化,即对每张图像,先人工点取三个特征点(两只眼睛的中心点和嘴巴的中心点)的坐标,通过旋转使得每张人脸的两只眼睛处于水平位置,再运用双插值算法拉伸图像,使得这三个特征点位于图像的固定位置,最后将图像裁剪为同一大小(见图2)。
(2)训练系数矩阵Ai和Bi,i=2…,64。用Yale B图库中的正面人脸采用最小二乘法模型②进行训练,其中每张图像的光照分量L用LTV模型分解获得(见图3)。
(3)对目标人脸图像进行分解。对每张待进行光照对齐的人脸图像Iobject,用LTV模型对其进行分解可以得到光照分量Lobject和反射分量Robject(见图3)。
(4)对目标人脸图像进行光照类型估计。直接用(3)中分解得到的光照分量Lobject通过PCA识别方法估计出Iobject的光照类型ii。
(5)对目标人脸图像的光照成分Lobject进行对齐矫正。根据(4)中估计得到的光照类型ii相应选取(2)中训练得到的系数矩阵Ai和Bi,用式⑦对Lobject进行对齐矫正(见图4),矫正结果为Laligned。
(6)目标人脸图像重构。用(3)中分解得到的反射分量Robject和(5)中光照对齐后得到的光照分量Laligned按⑧式重构最终图像。
其中步骤(3)~步骤(6)显示了本发明对具体某张人脸图像的光照对齐整体流程(见图5)。
图6为本发明方法与其他几种光照处理方法在视觉效果方面的比较。第一行为原图,第二行为LTV模型的结果,最后一行为我们的方法处理结果。由图6对比可知,用本发明提出的方法,恢复得到的图像与标准图像整体灰度比较接近,阴影部分已经基本消除,同时能保留大量人脸细节信息,真正达到了视觉效果的恢复,能符合各种实际应用要求。
本发明通过以下实验对本发明的效果进行说明:人脸识别实验分别在100×100和112×92象素大小的两种图库(对齐后的Yale B人脸库)上进行,所用的识别方法为主成分分析(PCA)识别方法,用欧氏距离描述图像的相似度。Yale B人脸库的图像根据光照类型可以分为5个子集,用第1个子集里面每个人若干张图片进行训练,对其余4个子集进行识别,用各种光照处理方法处理后的图像进行人脸识别对应识别率如下(其中线性光照模型即用线性关系来描述光照变化引起人脸图像的变化,具体见“Xudong Xie,Kin-Man Lam,’Facerecognition under varying illumination based on 2D face shape model,’PatternRecognition 38(2005)221-230”,而LTV就是基于LTV模型的光照恢复算法,也就是用LTV模型分解得到反射分量进行识别):
表1人脸识别结果(%),其中选取30个主成分,用第1个子集里面每个人7张图片做训练样本(100×100大小图片)
Set2 | Set3 | Set4 | Set5 | |
未进行光照恢复 | 97.50 | 73.33 | 25.71 | 10.53 |
LTV | 95.83 | 93.33 | 80.00 | 83.68 |
线性光照模型 | 99.17 | 100.00 | 73.57 | 21.58 |
本文方法 | 98.33 | 100.00 | 88.57 | 90.53 |
表2人脸识别结果(%),其中选取20个主成分,仅用第1个子集里面每个人3张图片做训练样本(100×100大小图片)
Set2 | Set3 | Set4 | Set5 | |
未进行光照恢复 | 96.67 | 70.83 | 26.43 | 12.10 |
LTV | 99.17 | 91.67 | 79.29 | 81.58 |
线性光照模型 | 99.17 | 100.00 | 69.29 | 26.84 |
本文方法 | 98.33 | 97.50 | 87.14 | 89.47 |
表3人脸识别结果(%),其中选取30个主成分,用第1个子集里面每个人7张图片做训练样本(112×92大小图片)
Set2 | Set3 | Set4 | Set5 | |
未进行光照恢复 | 95.83 | 69.17 | 38.00 | 16.32 |
LTV | 96.67 | 97.50 | 87.33 | 75.26 |
线性光照模型 | 99.17 | 95.83 | 62.14 | 26.84 |
本文方法 | 98.33 | 99.17 | 84.67 | 80.53 |
表4人脸识别结果(%),其中选取20个主成分,仅用第1个子集里面每个人3张图片做训练样本(112×92大小图片)
Set2 | Set3 | Set4 | Set5 | |
未进行光照恢复 | 96.67 | 70.83 | 35.71 | 17.37 |
LTV | 99.17 | 95.83 | 86.00 | 77.36 |
线性光照模型 | 99.17 | 95.83 | 62.14 | 25.26 |
本文方法 | 98.33 | 100 | 82.86 | 80.53 |
从实验结果可以看出,采用本发明进行光照恢复,在YaleB各个光照类型子集中均能较大地提高识别率。从识别结果可以得出,LTV模型分解得到的光照分量L中的确还存在可用的人脸细节信息,本文提出的方法利用这些信息,能再次提高人脸识别率。然而,对光照分量,仅仅用线性模型进行对齐矫正尚存在很大的误差。从上面识别实验可以看出,对光照变化比较大的图像,如第四、五个子集,用线性模型进行对齐,识别率并没有得到非常显著的提高。而本文选取二次多项式模型,从识别率可以看出确实二次多项式模型能较好地描述光照变化引起的人脸图像变换。
所述实施例为本发明较佳的实施方式,但本发明的实施方式并不受所述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (7)
1、一种基于二次多项式光照模型的人脸光照对齐方法,其特征在于包括以下步骤:
(1)建立光照模型;
(2)对所有训练图像以及目标图像均进行形状规一化处理;
(3)训练系数矩阵Ai和Bi,其中i=2,Λ,64;
(4)对目标人脸图像进行分解;
(5)对目标人脸图像进行光照类型估计;
(6)对目标人脸图像的光照分量进行对齐矫正;
(7)目标人脸图像重构。
2、根据权利要求1所述的一种基于二次多项式光照模型的人脸光照对齐方法,其特征在于:所述步骤(1)中光照模型是
3、根据权利要求2所述的一种基于二次多项式光照模型的人脸光照对齐方法,其特征在于:所述步骤(3)是采用最小二乘法对光照模型中系数矩阵Ai和Bi进行训练。
4、根据权利要求3所述的一种基于二次多项式光照模型的人脸光照对齐方法,其特征在于:所述步骤(4)具体是对每张待进行光照对齐的目标人脸图像Iobject,根据光照-反射模型Iobject(x,y)=Lobject(x,y)Robject(x,y),用LTV模型对Iobject进行分解,得到光照分量Lobject和反射分量Robject。
5、根据权利要求4所述的一种基于二次多项式光照模型的人脸光照对齐方法,其特征在于,所述步骤(5)具体是对光照分量Lobject,直接用PCA识别方法估计出Iobject的光照类型。
6、根据权利要求5所述的一种基于二次多项式光照模型的人脸光照对齐方法,其特征在于,所述步骤(6)具体是根据Iobject的光照类型相应选取训练得到的系数矩阵Ai和Bi,用 i=2,…,64对Lobject进行对齐矫正,并得到矫正结果Laligned。
7、根据权利要求6所述的一种基于二次多项式光照模型的人脸光照对齐方法,其特征在于,所述步骤(7)具体是用分解得到的反射分量Robject和光照对齐后得到的光照分量Laligned重构最终图像:Ialigned(x,y)=Laligned(x,y)Robject(x,y)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100278178A CN100444193C (zh) | 2007-04-29 | 2007-04-29 | 一种基于二次多项式光照模型的人脸光照对齐方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100278178A CN100444193C (zh) | 2007-04-29 | 2007-04-29 | 一种基于二次多项式光照模型的人脸光照对齐方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101046847A true CN101046847A (zh) | 2007-10-03 |
CN100444193C CN100444193C (zh) | 2008-12-17 |
Family
ID=38771445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007100278178A Expired - Fee Related CN100444193C (zh) | 2007-04-29 | 2007-04-29 | 一种基于二次多项式光照模型的人脸光照对齐方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100444193C (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101916384A (zh) * | 2010-09-01 | 2010-12-15 | 汉王科技股份有限公司 | 一种人脸图像重建方法、装置及人脸识别系统 |
CN101916371A (zh) * | 2010-09-01 | 2010-12-15 | 北京工业大学 | 对图像进行光照归一化处理的方法及采用其的图像识别方法 |
CN102279925A (zh) * | 2011-08-25 | 2011-12-14 | 三峡大学 | 链式处理人脸识别方法及系统 |
CN102346857A (zh) * | 2011-09-14 | 2012-02-08 | 西安交通大学 | 人脸图像光照参数和去光照图的高精度同时估计方法 |
CN101425179B (zh) * | 2008-11-18 | 2012-03-28 | 清华大学 | 一种人脸图像重光照的方法及装置 |
CN104424483A (zh) * | 2013-08-21 | 2015-03-18 | 中移电子商务有限公司 | 一种人脸图像的光照预处理方法、装置及终端 |
WO2017024963A1 (zh) * | 2015-08-11 | 2017-02-16 | 阿里巴巴集团控股有限公司 | 图像识别方法、度量学习方法、图像来源识别方法及装置 |
CN109040729A (zh) * | 2018-08-16 | 2018-12-18 | Oppo广东移动通信有限公司 | 图像白平衡校正方法、装置、存储介质及终端 |
CN110647859A (zh) * | 2019-09-29 | 2020-01-03 | 浙江商汤科技开发有限公司 | 人脸图像分解方法和装置、电子设备及存储介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100375105C (zh) * | 2005-04-28 | 2008-03-12 | 中国科学院自动化研究所 | 一种基于图像比的光照变化下的人脸识别方法 |
-
2007
- 2007-04-29 CN CNB2007100278178A patent/CN100444193C/zh not_active Expired - Fee Related
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101425179B (zh) * | 2008-11-18 | 2012-03-28 | 清华大学 | 一种人脸图像重光照的方法及装置 |
CN101916371A (zh) * | 2010-09-01 | 2010-12-15 | 北京工业大学 | 对图像进行光照归一化处理的方法及采用其的图像识别方法 |
CN101916384A (zh) * | 2010-09-01 | 2010-12-15 | 汉王科技股份有限公司 | 一种人脸图像重建方法、装置及人脸识别系统 |
CN101916371B (zh) * | 2010-09-01 | 2012-11-21 | 北京工业大学 | 对图像进行光照归一化处理的方法及采用其的图像识别方法 |
CN101916384B (zh) * | 2010-09-01 | 2012-11-28 | 汉王科技股份有限公司 | 一种人脸图像重建方法、装置及人脸识别系统 |
CN102279925B (zh) * | 2011-08-25 | 2013-09-11 | 三峡大学 | 链式处理人脸识别方法及系统 |
CN102279925A (zh) * | 2011-08-25 | 2011-12-14 | 三峡大学 | 链式处理人脸识别方法及系统 |
CN102346857A (zh) * | 2011-09-14 | 2012-02-08 | 西安交通大学 | 人脸图像光照参数和去光照图的高精度同时估计方法 |
CN102346857B (zh) * | 2011-09-14 | 2014-01-15 | 西安交通大学 | 人脸图像光照参数和去光照图的高精度同时估计方法 |
CN104424483A (zh) * | 2013-08-21 | 2015-03-18 | 中移电子商务有限公司 | 一种人脸图像的光照预处理方法、装置及终端 |
WO2017024963A1 (zh) * | 2015-08-11 | 2017-02-16 | 阿里巴巴集团控股有限公司 | 图像识别方法、度量学习方法、图像来源识别方法及装置 |
CN109040729A (zh) * | 2018-08-16 | 2018-12-18 | Oppo广东移动通信有限公司 | 图像白平衡校正方法、装置、存储介质及终端 |
CN109040729B (zh) * | 2018-08-16 | 2020-04-07 | Oppo广东移动通信有限公司 | 图像白平衡校正方法、装置、存储介质及终端 |
CN110647859A (zh) * | 2019-09-29 | 2020-01-03 | 浙江商汤科技开发有限公司 | 人脸图像分解方法和装置、电子设备及存储介质 |
CN110647859B (zh) * | 2019-09-29 | 2023-01-31 | 浙江商汤科技开发有限公司 | 人脸图像分解方法和装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN100444193C (zh) | 2008-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101046847A (zh) | 一种基于二次多项式光照模型的人脸光照对齐方法 | |
CN101673338B (zh) | 基于多角度投影的模糊车牌识别方法 | |
CN1908960A (zh) | 一种基于特征分组的多分类器组合人脸识别方法 | |
CN113392711B (zh) | 一种基于高层语义与噪声抑制的烟雾语义分割方法及系统 | |
CN1932882A (zh) | 基于目标检测的红外与可见光序列图像特征级融合方法 | |
CN1885341A (zh) | 一种用于图像认证的可逆水印方法 | |
CN102147812A (zh) | 基于三维点云模型的地标建筑图像分类方法 | |
CN103646381B (zh) | 一种线阵ccd的行进畸变校正方法 | |
CN1438603A (zh) | 基于通用办公设备的二维条码自动识读方法 | |
CN1878319A (zh) | 基于平面单应矩阵特征线的摄像机标定方法 | |
CN101261678A (zh) | 在不同尺度特征图像上进行人脸光照规一化的方法 | |
CN106127763A (zh) | 一种具有广泛适应性的图像二值化方法 | |
CN1921562A (zh) | 基于变换域数学形态学的图像降噪方法 | |
CN107545259A (zh) | 一种基于大津法的二维码重构方法 | |
CN1917577A (zh) | 一种图像组合降噪方法 | |
CN1920881A (zh) | 一种Contourlet变换域的图像降噪方法 | |
CN110121109A (zh) | 面向监控系统数字视频实时溯源方法、城市视频监控系统 | |
CN1885339A (zh) | 一种基于图像特征的电子印章认证方法 | |
CN101039371A (zh) | 用于印刷品版权保护的数字水印新方法 | |
CN106296599A (zh) | 一种自适应图像增强方法 | |
CN1975762A (zh) | 一种皮肤检测方法 | |
CN117853725A (zh) | 一种木料区域分割方法、系统及存储介质 | |
CN1246779C (zh) | 小波域混沌半脆弱数字水印嵌入与提取方法 | |
CN1525387A (zh) | 检测图像模糊性的设备和方法 | |
CN113724156B (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20081217 Termination date: 20100429 |