CN101833654A - 基于约束采样的稀疏表示人脸识别方法 - Google Patents

基于约束采样的稀疏表示人脸识别方法 Download PDF

Info

Publication number
CN101833654A
CN101833654A CN 201010140797 CN201010140797A CN101833654A CN 101833654 A CN101833654 A CN 101833654A CN 201010140797 CN201010140797 CN 201010140797 CN 201010140797 A CN201010140797 A CN 201010140797A CN 101833654 A CN101833654 A CN 101833654A
Authority
CN
China
Prior art keywords
eyedistance
face
training set
nosetip
sampling
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
Application number
CN 201010140797
Other languages
English (en)
Other versions
CN101833654B (zh
Inventor
苏光大
王晶
熊英
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JIANGSU TSINGDA VISION TECHNOLOGY Co Ltd
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN2010101407977A priority Critical patent/CN101833654B/zh
Publication of CN101833654A publication Critical patent/CN101833654A/zh
Application granted granted Critical
Publication of CN101833654B publication Critical patent/CN101833654B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)

Abstract

本发明涉及基于约束采样的稀疏表示人脸识别方法,属于图像处理领域,该方法包括:对训练集中所有人脸图像分别基于约束区域进行采样进行特征提取,得到的特征矩阵,特征矩阵的一个类别对应训练集的一个人的多幅人脸图像;对待识别人的人脸图像基于约束区域进行采样进行特征提取,得到待识别图像的特征向量;该特征向量用训练集人脸图像的特征向量线性系数;对于训练集的每一类别计算对应于该类别的残差;残差的最小值对应的训练集的类别作为待识别人的识别结果;该方法通过对约束区域进行采样,使其具有更高的人脸识别率。

Description

基于约束采样的稀疏表示人脸识别方法
技术领域
本发明属于图像处理、计算机视觉、模式识别技术领域,特别涉及人脸识别方法。
背景技术
生物特征识别技术是身份识别的有效技术,近来发展最快的是人脸识别技术以及与人脸识别技术相融合的生物特征识别技术。
目前已有的人脸识别方法主要是对整个人脸进行识别的,而在诸多的识别方法中,主要采用主分量分析(PCA-Principal Component Analysis)、弹性匹配、神经网络、几何特征等方法。
同时,人脸识别的难点还在于:
(1)表情引起的人脸塑性变形
(2)姿态引起的人脸多样性
(3)年龄引起的人脸变化
(4)发型、胡须、眼镜、化装等因素引起的人脸模式的多重性
(5)光照的角度、强度以及传感器特性等因素引起的人脸图像的差异性
诸多的因素使得人脸识别成为一项棘手而富挑战性的课题,也因此在近年成为科研的热点。
目前已有的一种基于稀疏表示的人脸识别方法:其主要特点是将稀疏信号表示这一新理论应用于人脸识别。该方法首先假设待识别的人脸可以由训练集中该类别的人脸线性表示,利用这一假设,通过计算L1范数最小值求得稀疏表达。该方法具体包括以下步骤:
1)对于训练集中每一幅人脸图像进行处理,提取特征向量,将训练集的所有特征向量排列形成特征矩阵A,将特征矩阵中的一行或多行作为训练集的一个类别,一个类别对应训练集的一个人的多幅人脸图像;
2)对于待识别的人脸图像按照同样的方法提取特征向量y;
3)将待识别人脸图像的特征向量用训练集图像的特征向量线性表示,求出线性表示的系数;(由于这种线性表示的系数具有多样性)计算如下L1范数最小值问题,选择其中L1范数最小的系数作为线性表示的系数
Figure GSA00000075020700011
x ^ = arg min | | x | | 1 , 满足||y-Ax||2≤ε    (7)
其中ε为预设的误差,||y-Ax||2为向量y-Ax的L2范数;
4)对于训练集的每一类别i=1,…,k(k为训练集的类别数),根据步骤3)得到的系数
Figure GSA00000075020700022
按照公式(8)计算对应于该类别的残差ri(y);
r i ( y ) = | | y - Aδ i ( x ^ ) | | 2 - - - ( 8 )
其中,
Figure GSA00000075020700024
对于第i类,δi(x)=x;对于其他类,δi(x)=0;
5)比较步骤4)求出每一类的对应残差,残差的最小值对应的训练集的类别作为待识别人的识别结果。
该方法对人脸光照和表情的变化具有鲁棒性,但是对图像的配准有较高的要求。当图像之间的配准有几个像素的偏差时,它的性能就会受到很大的影响。
与本发明相关的已有技术说明如下:
人脸检测:利用Real Adaboost方法与新的特征定义形式以及构造方式得到层次型检测器,并选择计算次数少的特征作为候选特征以提高检测速度。
人脸关键特征点的提取:其主要特点是利用SUSAN算子提取边缘,并与方向积分投影方法结合,获得眼球、鼻尖、下巴等一系列关键特征点位置。实验结果表明,该方法平均定位准确率在95%以上。
发明内容
本发明的目的是为了提高人脸识别算法的鲁棒性,利用人脸关键点的精确定位结果,提出一种基于约束采样的稀疏表示人脸识别方法,该方法通过对约束区域进行采样,使其具有更高的人脸识别率。
本发明提出的基于约束采样的稀疏表示人脸识别方法,其特征在于,该方法包括以下步骤:
1)对训练集中所有人脸图像分别进行特征提取,得到训练集人脸图像的特征向量,将训练集的所有人脸图像的特征向量排列形成特征矩阵,将特征矩阵中的一行或多行作为训练集的一个类别,一个类别对应训练集的一个人的多幅人脸图像;
2)对待识别人的人脸图像进行特征提取,得到待识别图像的特征向量;
3)将待识别人脸图像的特征向量用训练集人脸图像的特征向量线性表示,求出线性表示的系数;选择其中L1范数最小的系数作为线性表示的系数;
4)对于训练集的每一类别,根据步骤3)得到的系数计算对应于该类别的残差;
5)按照步骤4)求出训练集的每一类的对应残差,残差的最小值对应的训练集的类别作为待识别人的识别结果;
所述步骤1)或步骤2)中的进行特征提取,具体包括:
(a)对于给定的一幅人脸图像,通过人脸检测定位其人脸区域;定位两个眼球点,将图像尺度归一化,并进行旋转,使两眼呈水平状;
(b)选取人脸区域中的两个以上的关键特征点进行定位;
(c)对约束区域进行采样:根据所述关键特征点的位置选定眼睛、鼻尖和嘴巴区域;对这些约束区域进行均匀采样,得到人脸图像的纹理特征;
(d)将约束采样得到的人脸图像纹理特征排列成向量,并进行L2范数归一化,得到人脸图像的特征向量。
本发明的特点及效果
本发明根据人脸的关键特征点,划定眼睛、鼻尖、嘴巴区域,再通过约束采样提取出纹理特征。接着根据稀疏表示的理论进行识别的矩阵运算,满足了原算法对配准的要求,具有更高的人脸识别率。
具体实施方式
本发明提出的基于约束采样的稀疏表示人脸识别方法结合实施例详细说明如下:本实施例方法包括以下步骤:
1)对训练集中所有人脸图像分别进行特征提取,得到训练集人脸图像的特征向量,将训练集的所有人脸图像的特征向量排列形成特征矩阵A,将特征矩阵中的一行或多行作为训练集的一个类别,一个类别对应训练集的一个人的多幅人脸图像;
2)对待识别人的人脸图像进行特征提取,得到待识别图像的特征向量y;
3)将待识别图像的特征向量用训练集图像的特征向量线性表示,求出线性表示的系数;(由于这种表示具有多样性)选择其中L1范数最小的系数作为线性表示的系数;步骤1)-2)得到训练矩阵A,特征向量y,步骤3)就是求解如下线性方程:
Ax=y    (9)
(由于通常情况下,训练矩阵A的行数小于列数,方程(9)为不定方程,有多个非零解,必须增加约束条件)定义x的L1范数||x||1为:
||x||1=|x1|+|x2|+...+|xn|    (10)
其中xj(j=1,2,...,n)为向量x的一个分量,|xj|为xj的绝对值,n为向量x的维数。
根据公式(7),计算L1范数最小的系数作为线性表示的系数。
4)对于训练集的每一类别,根据3)得到的系数计算对应于该类别的残差;
对于每一类别i=1,…,k,提取出对应于该类的系数,与该类的特征向量线性组合,再与待识别的特征向量比较,分别按照(8)式计算其残差ri。
6)按照4)求出每一类对应的残差ri(i=1,…,k),残差的最小值对应的训练集的类别作为待识别人的识别结果;
上述步骤1)或步骤2)中的进行特征提取,具体包括:
a)对于储存在计算机的人脸图像,根据人脸检测定位其人脸区域;定位两眼球点,并将图像尺度归一化到360×480像素,使人脸位于图像中间,且下颌点到两眼中心连线的距离大约为200个像素(该尺度只作为举例说明,并不对图像大小进行约束),通过旋转,使双眼在同一水平线上;
b)定位双眼球、鼻尖、下巴四个关键特征点(四个点的选取只作为举例说明,实际应用中也可以选取两个以上的其他部位的关键点,随着关键点选取的不同,性能将有微小差别);
本实施例采用SUSAN算子提取边缘,并与方向积分投影方法结合,获得眼球、鼻尖、下巴四个关键特征点的位置;
c)对约束区域进行采样:根据所述关键特征点的位置选定眼睛、鼻尖和嘴巴区域;对这些约束区域进行均匀采样,得到人脸图像的纹理特征;
本实施例具体的做法是:
用LeftEye_x,LeftEye_y分别表示左眼球的横、纵坐标值;RightEye_x,RightEye_y分别表示右眼球的横、纵坐标值;NoseTip_x、NoseTip_y表示鼻尖的横、纵坐标值;Chin_x、Chin_y分别表示下颌点的横、纵坐标值。按照公式(1)-(6),取参数a1=0.3,a2=0,4,a3=0.3,a4=0.2,b1=0.35,b2=0,1,b3=0.3,c1=0.5,c2=0,2,c3=0.8,分别得到左眼区域ImgEyeLeft,右眼区域ImgEyeRight,鼻尖区域ImgNose以及EyeDistance嘴巴区域。
接着根据预设的采样间隔r=6(r≥1),对这四个区域进行固定点均匀采样,得到人脸图像的纹理特征(属于常规技术);
d)将约束采样得到的人脸图像纹理特征排列成向量,并进行L2范数归一化,得到人脸图像的特征向量;
向量x的L2范数||x||2定义为:
| | x | | 2 = ( x 1 2 + x 2 3 + . . . + x n 2 ) - - - ( 11 )
其中xj(j=1,2,...,n)为向量x的一个分量,n为向量x的维数。
基于约束采样的稀疏表示人脸识别方法在AR库以及中科院人脸表情库进行测试。AR库选取90个不同类别的人脸,其中每一类包含14张不同光照和表情的正面图像,实验中7张训练另外7张用来测试,所达到的识别率为:98.1%,(纹理特征维数为432);中科院人脸库选取202类不同人脸,其中每一类包含6张不同表情的正面图像,实验中任意选择4张训练,另外2张用来测试,所达到识别率为:98.51%(纹理特征维数768)。

Claims (2)

1.一种基于约束采样的稀疏表示人脸识别方法,其特征在于,该方法包括以下步骤:
1)对训练集中所有人脸图像分别进行特征提取,得到训练集人脸图像的特征向量,将训练集的所有人脸图像的特征向量排列形成特征矩阵,将特征矩阵中的一行或多行作为训练集的一个类别,一个类别对应训练集的一个人的多幅人脸图像;
2)对待识别人的人脸图像进行特征提取,得到待识别图像的特征向量;
3)将待识别人脸图像的特征向量用训练集人脸图像的特征向量线性表示,求出线性表示的系数;选择其中L1范数最小的系数作为线性表示的系数;
4)对于训练集的每一类别,根据步骤3)得到的系数计算对应于该类别的残差;
5)按照步骤4)求出训练集的每一类的对应残差,残差的最小值对应的训练集的类别作为待识别人的识别结果;
所述步骤1)或步骤2)中的进行特征提取,具体包括:
(a)对于给定的一幅人脸图像,通过人脸检测定位其人脸区域;定位两个眼球点,将图像尺度归一化,并进行旋转,使两眼呈水平状;
(b)选取人脸区域中的两个以上的关键特征点进行定位;
(c)对约束区域进行采样:根据所述关键特征点的位置选定眼睛、鼻尖和嘴巴区域;对这些约束区域进行均匀采样,得到人脸图像的纹理特征;
(d)将约束采样得到的人脸图像纹理特征排列成向量,并进行L2范数归一化,得到人脸图像的特征向量。
2.如权利要求1所述的方法,其特征在于,所述步骤(c)对约束区域进行采样,具体包括:
用LeftEye_x,LeftEye_y分别表示左眼球的横、纵坐标值;RightEye_x,RightEye_y分别表示右眼球的横、纵坐标值;NoseTip_x、NoseTip_y表示鼻尖的横、纵坐标值;Chin_x、Chin_y分别表示下颌点的横、纵坐标值;定义EyeDistance为双眼距离:
EyeDistance=RightEye_x-LeftEye_x                            (1)
左眼区域ImgEyeLeft定义为一个矩形区域,它表示为:
ImgEyeLeft=LeftEye_x-a1×EyeDistance,LeftEye_y-a2×EyeDistance
                                                                    (2)
LeftEye_x+a3×EyeDistance,LeftEye_y+a4×EyeDistance
其中等式右边分别为矩形区域左上角横坐标,纵坐标,右下角的横坐标,纵坐标。系数a1,a3的取值范围为[0.2,0.4],a2,a4取值范围为[0.3,0.5]。同理,右眼区域ImgEyeRight表示为:
ImgEyeRight=RightEye_x-a1×EyeDistance,RightEye_y-a2×EyeDistance
                                                                      (3)
RighttEye_x+a3×EyeDistance,RightEye_y+a4×EyeDistance
其中系数a1,a2,a3,a4取值同式(2);
鼻尖区域ImgNose表示为:
ImgNose=NoseTip_x-b1×EyeDistance,NoseTip_y-b2×EyeDistance
                                                                      (4)
NoseTip_x+b1×EyeDistance,NoseTip_y+b3×EyeDistance
其中系数b1的取值范围为[0.3,0.4],b2,b3取值范围为[0,0.3]。定义MouthDistance为嘴巴区域的垂直距离:
MouthDistance=Chin_y-NoseTip_y                                       (5)
嘴巴区域ImgMouth定义为一矩形区域,它表示为:
ImgMouth=NoseTip_x-c1×EyeDistance,NoseTip_y+c2×MouthDistance
                                                                      (6)
NoseTip_x+c1×EyeDistance,NoseTip_y+c3×MouthDistance
其中系数c1的取值范围为[0.4,0.6],c2取值范围为[0.2,0.4],c3的取值范围为[0.7,0.9]。接着根据预设的采样间隔r(r≥1),对这四个区域进行固定点均匀采样,得到人脸图像的纹理特征。
CN2010101407977A 2010-04-02 2010-04-02 基于约束采样的稀疏表示人脸识别方法 Expired - Fee Related CN101833654B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101407977A CN101833654B (zh) 2010-04-02 2010-04-02 基于约束采样的稀疏表示人脸识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101407977A CN101833654B (zh) 2010-04-02 2010-04-02 基于约束采样的稀疏表示人脸识别方法

Publications (2)

Publication Number Publication Date
CN101833654A true CN101833654A (zh) 2010-09-15
CN101833654B CN101833654B (zh) 2012-01-11

Family

ID=42717718

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101407977A Expired - Fee Related CN101833654B (zh) 2010-04-02 2010-04-02 基于约束采样的稀疏表示人脸识别方法

Country Status (1)

Country Link
CN (1) CN101833654B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102073880A (zh) * 2011-01-13 2011-05-25 西安电子科技大学 利用稀疏表示进行人脸识别的集成方法
CN102404581A (zh) * 2011-11-02 2012-04-04 清华大学 基于插值与近红外的彩色图像处理的方法和装置
CN103167606A (zh) * 2013-03-12 2013-06-19 钱钢 基于稀疏表示的wlan室内定位方法
CN103345621A (zh) * 2013-07-09 2013-10-09 东南大学 一种基于稀疏浓度指数的人脸分类方法
CN103530657A (zh) * 2013-09-26 2014-01-22 华南理工大学 一种基于加权l2抽取深度学习人脸识别方法
CN103984922A (zh) * 2014-04-30 2014-08-13 苏亚 一种基于稀疏表示和形状约束的人脸识别方法
CN104143090A (zh) * 2014-07-30 2014-11-12 哈尔滨工业大学深圳研究生院 一种基于人脸识别的汽车开门方法
CN104166860A (zh) * 2014-07-25 2014-11-26 哈尔滨工业大学深圳研究生院 基于约束的面向单个测试样本的人脸识别方法
CN105320964A (zh) * 2015-10-26 2016-02-10 中国矿业大学(北京) 基于稀疏表示的煤岩识别方法
CN107563328A (zh) * 2017-09-01 2018-01-09 广州智慧城市发展研究院 一种基于复杂环境下的人脸识别方法及系统
CN110334587A (zh) * 2019-05-23 2019-10-15 北京市威富安防科技有限公司 人脸关键点定位模型的训练方法、装置及关键点定位方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101667246A (zh) * 2009-09-25 2010-03-10 西安电子科技大学 基于核稀疏表示的人脸识别方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101667246A (zh) * 2009-09-25 2010-03-10 西安电子科技大学 基于核稀疏表示的人脸识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《北方交通大学学报》 20010430 李华胜, 杨 桦, 袁保宗 人脸识别系统中的特征提取 第25卷, 第2期 2 *
《西安电子科技大学硕士学位论文》 20090701 刘彬 基于随机投影和稀疏表征的红外人脸识别方法 , 2 *
《计算机工程》 20040505 汪 宁,丁晓青 结合整体与局部信息的人脸识别方法 第30卷, 第5期 2 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102073880A (zh) * 2011-01-13 2011-05-25 西安电子科技大学 利用稀疏表示进行人脸识别的集成方法
CN102404581A (zh) * 2011-11-02 2012-04-04 清华大学 基于插值与近红外的彩色图像处理的方法和装置
CN103167606B (zh) * 2013-03-12 2015-12-23 钱钢 基于稀疏表示的wlan室内定位方法
CN103167606A (zh) * 2013-03-12 2013-06-19 钱钢 基于稀疏表示的wlan室内定位方法
CN103345621A (zh) * 2013-07-09 2013-10-09 东南大学 一种基于稀疏浓度指数的人脸分类方法
CN103345621B (zh) * 2013-07-09 2016-05-18 东南大学 一种基于稀疏浓度指数的人脸分类方法
CN103530657A (zh) * 2013-09-26 2014-01-22 华南理工大学 一种基于加权l2抽取深度学习人脸识别方法
CN103530657B (zh) * 2013-09-26 2016-07-06 华南理工大学 一种基于加权l2抽取深度学习人脸识别方法
CN103984922A (zh) * 2014-04-30 2014-08-13 苏亚 一种基于稀疏表示和形状约束的人脸识别方法
CN103984922B (zh) * 2014-04-30 2017-04-26 苏亚 一种基于稀疏表示和形状约束的人脸识别方法
CN104166860A (zh) * 2014-07-25 2014-11-26 哈尔滨工业大学深圳研究生院 基于约束的面向单个测试样本的人脸识别方法
CN104166860B (zh) * 2014-07-25 2017-06-13 哈尔滨工业大学深圳研究生院 基于约束的面向单个测试样本的人脸识别方法
CN104143090A (zh) * 2014-07-30 2014-11-12 哈尔滨工业大学深圳研究生院 一种基于人脸识别的汽车开门方法
CN104143090B (zh) * 2014-07-30 2018-02-02 哈尔滨工业大学深圳研究生院 一种基于人脸识别的汽车开门方法
CN105320964A (zh) * 2015-10-26 2016-02-10 中国矿业大学(北京) 基于稀疏表示的煤岩识别方法
CN105320964B (zh) * 2015-10-26 2019-06-04 中国矿业大学(北京) 基于稀疏表示的煤岩识别方法
CN107563328A (zh) * 2017-09-01 2018-01-09 广州智慧城市发展研究院 一种基于复杂环境下的人脸识别方法及系统
CN110334587A (zh) * 2019-05-23 2019-10-15 北京市威富安防科技有限公司 人脸关键点定位模型的训练方法、装置及关键点定位方法

Also Published As

Publication number Publication date
CN101833654B (zh) 2012-01-11

Similar Documents

Publication Publication Date Title
CN101833654B (zh) 基于约束采样的稀疏表示人脸识别方法
CN101819628B (zh) 结合形状特征的稀疏表示人脸识别方法
CN101833672B (zh) 基于约束采样与形状特征的稀疏表示人脸识别方法
CN108549873B (zh) 三维人脸识别方法和三维人脸识别系统
US10049262B2 (en) Method and system for extracting characteristic of three-dimensional face image
CN107169455B (zh) 基于深度局部特征的人脸属性识别方法
US7512255B2 (en) Multi-modal face recognition
CN107403142B (zh) 一种微表情的检测方法
Al-Osaimi et al. An expression deformation approach to non-rigid 3D face recognition
TWI383325B (zh) 臉部表情辨識
CN103218609B (zh) 一种基于隐最小二乘回归的多姿态人脸识别方法及其装置
CN102013011B (zh) 基于正脸补偿算子的多姿态人脸识别方法
CN103136516B (zh) 可见光与近红外信息融合的人脸识别方法及系统
Cheng et al. Image-to-class dynamic time warping for 3D hand gesture recognition
CN108182397B (zh) 一种多姿态多尺度的人脸验证方法
CN108629336A (zh) 基于人脸特征点识别的颜值计算方法
CN106096517A (zh) 一种基于低秩矩阵与特征脸的人脸识别方法
CN101840509A (zh) 人眼观察视角的测量方法及装置
CN109409298A (zh) 一种基于视频处理的视线追踪方法
CN107480586A (zh) 基于人脸特征点位移的生物识别照片仿冒攻击检测方法
CN110110603A (zh) 一种基于面部生理信息的多模态唇读方法
CN103093237A (zh) 一种基于结构化模型的人脸检测方法
CN112232128A (zh) 基于视线追踪的老年残障人士照护需求识别方法
Kare et al. Using bidimensional regression to assess face similarity
CN107862246A (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
ASS Succession or assignment of patent right

Owner name: JIANGSU TSINGDA VISION TECHNOLOGY CO., LTD.

Free format text: FORMER OWNER: TSINGHUA UNIVERSITY

Effective date: 20130508

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100084 HAIDIAN, BEIJING TO: 210000 NANJING, JIANGSU PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20130508

Address after: 210000 Jiangsu Baixia District of Nanjing City Shi Yang Road No. 56 A District No. 2 building 019 room

Patentee after: Jiangsu Tsingda Vision Technology Co., Ltd.

Address before: 100084 Haidian District Tsinghua Yuan Beijing No. 1

Patentee before: Tsinghua University

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120111

Termination date: 20210402