CN107088069A - 基于人体ppg信号分段的身份识别方法 - Google Patents
基于人体ppg信号分段的身份识别方法 Download PDFInfo
- Publication number
- CN107088069A CN107088069A CN201710195051.8A CN201710195051A CN107088069A CN 107088069 A CN107088069 A CN 107088069A CN 201710195051 A CN201710195051 A CN 201710195051A CN 107088069 A CN107088069 A CN 107088069A
- Authority
- CN
- China
- Prior art keywords
- mrow
- waveform
- collection
- sub
- msubsup
- 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 42
- 238000012360 testing method Methods 0.000 claims abstract description 64
- 238000012549 training Methods 0.000 claims abstract description 52
- 230000004927 fusion Effects 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims description 54
- 239000013598 vector Substances 0.000 claims description 46
- 239000000203 mixture Substances 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 238000005192 partition Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 238000012706 support-vector machine Methods 0.000 claims description 4
- 230000007812 deficiency Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 239000012141 concentrate Substances 0.000 claims 2
- 210000000056 organ Anatomy 0.000 abstract 1
- 230000003278 mimic effect Effects 0.000 description 6
- 244000287680 Garcinia dulcis Species 0.000 description 5
- 239000000284 extract Substances 0.000 description 3
- 108010024605 neurotensin mimic 2 Proteins 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000003542 behavioural effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000035479 physiological effects, processes and functions Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 125000001967 indiganyl group Chemical group [H][In]([H])[*] 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000004816 latex Substances 0.000 description 1
- 229920000126 latex Polymers 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/117—Identification of persons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
Abstract
本发明公开了一种基于人体PPG信号分段的身份识别方法，主要解决现有基于PPG信号身份识别率较低的问题。实现步骤：1)获取训练数据库和测试数据；2)对训练数据预处理，建立单周期波形库；3)对单周期波形分段，获取各子波形集及权重因子；4)利用DNMF方法对各子波形集分解，获取基空间和训练子特征集；5)利用权重因子对训练子特征融合，获取训练模板库；6)对测试数据预处理并分段，获取子波形集及权重因子；7)将测试数据子波形在基空间上投影，获取测试子特征集；8)对测试子特征加权融合，获取测试特征集，利用SVM分类器完成身份识别。本发明识别率在99.62％以上，可应用于金融、政府机构等领域的个人身份识别。
Description
技术领域
本发明属于信息处理技术领域，具体涉及一种身份识别方法，可作为金融、政府机构等领域中维护个人信息安全的一种手段。
背景技术
当今社会，安全问题越来越突出，人们不得不记忆复杂的密码或者携带额外的电子密码器，这使得传统的身份识别方法越来越失去它的实用性和可靠性，这一现状使人们对生物识别的需求越来越大。今天，大多数的系统比如金融交易、计算机网络和对安全领域的访问系统仍是通过身份证或口令进行识别授权的。这样的系统并不足够安全，因为身份证或口令信息很容易被窃取或者遗忘。生物识别系统可以提供更多的可靠性和隐秘性，因为它是根据个人的生理信号和行为特征进行身份认证的，这种生理信号或行为特征是个人独有的，并可以用来区分开不同个体。人体或行为属性独特的特性，如指纹、人脸、声音、脑电图和心电图等被用来进行身份识别。基于这些特征的应用提供了一种有发展前景和不可取代的识别方法。然而，指纹可以用乳胶提取特征，人脸识别可以用伪造的照片进行欺骗，声音可以被模仿，基于脑电信号或者心电信号的方法因需要各种各样的电极采集信号而不能广泛使用。
光电容积脉搏波PPG信号是一种非入侵式的光电方法，通过靠近皮肤测试身体的某一部位，获取关于血管中血液流动体积变化的信息。PPG信号作为人体固有的一种生理特征，具有难以被复制和模仿的特点，具有较高的安全性，且采集简单。目前基于PPG信号的时域身份识别方法，识别率不足够高，难以满足实际应用需求。
目前已提出的基于PPG信号的身份识别方法有：
A.Kemal Polatb等人2014年在“Computers in Biology&medicine”期刊上发表的“Anovel feature ranking algorithm for biometricrecognition with PPG signals”一文，提出了一种利用PPG信号波形特征点进行身份识别的方法，该方法首先从PPG信号单周期波形、一阶差分波形和二阶差分波形中提取了40维时域特征，然后根据这40维特征对正确识别率的贡献度，选取贡献度大的特征组成特征向量作为分类器的输入，并利用K最近邻分类器进行分类识别，正确识别率最大时为94.44％，不能满足识别率要求高的场合。
NI MohammedNadzr，M Sulaimi，LF Umadi，KA Sidek等人2016年在“IndianJournal of Science and Technology”期刊上发表的文章“PhotoplethysmogramBasedBiometric Recognition forTwins”中，研究了一种利用PPG信号对双胞胎的身份进行识别的方法，该方法首先利用低通滤波器对原始PPG信号进行去噪，然后对PPG信号波形进行分割，提取单周期波形，再利用径向基函数网络和朴素贝叶斯分类器分别对单周期波形进行识别分类，最终身份正确识别率达到97％以上，该方法验证了PPG信号的单周期波形特征对个体身份识别的有效性，但身份识别率仍有待进一步提升。
发明内容
本发明的目的在于针对上述已有技术的不足，提出一种基于人体PPG信号分段的身份识别方法，以提高身份识别的正确率。
本发明的技术方案是通过对人体PPG信号单周期波形进行分段处理，再利用鉴别式非负矩阵分解DNMF方法获取各子波段的特征向量，最后将各子波段的特征向量加权融合，生成融合的特征向量，进行身份识别，其实现步骤如下：
(1)获取训练数据库和测试数据。采集M个人在规定时间段内的光电容积脉搏波PPG信号，组成训练数据库S；再采集其中一人在另一时间段内的PPG信号，作为被鉴定者的测试数据xG；
(2)对训练数据库S依次进行去噪，归一化处理，波峰检测，波形分割，插值，去除差异性大的波形和进行波形平均，得到单周期平均波形数据库W；
(3)将单周期平均波形数据库W中每一个单周期波形进行分段处理，得到训练数据的上子波形集V1、中子波形集V2和下子波形集V3，并计算各子波形集中同一个人两两子波形之间的相似度，得到训练数据上子波形集的权重因子d1、中子波形集的权重因子d2和下子波形集的权重因子d3；
(4)利用鉴别式非负矩阵分解DNMF方法分别对训练数据的上子波形集V1、中子波形集V2和下子波形集V3进行分解，得到训练数据上子波形集的基空间Z1、中子波形集的基空间Z2、下子波形集的基空间Z3、训练数据上子波形特征集H1、中子波形特征集H2和下子波形特征集H3；
(5)利用训练数据上子波形集的权重因子d1，中子波形集的权重因子d2和下子波形集的权重因子d3分别对训练数据上子波形特征集H1，中子波形特征集H2下子波形特征集H3中相应的子特征进行加权融合，得到训练模板库H；
(6)对被xG鉴定者的测试数据依次进行步骤(2)-(3)操作，得到测试数据xG的上子波形集Α1，中子波形集Α2，下子波形集Α3和测试数据上子波形集的权重因子a1，测试数据中子波形集的权重因子a2和测试数据下子波形集的权重因子a3；
(7)将测试数据xG的上子波形集Α1，中子波形集Α2，下子波形集Α3分别在训练数据上子波形集的基空间Z1，中子波形集的基空间Z2，下子波形集的基空间Z3上进行投影，获得测试数据上子波形特征集F1，中子波形特征集F2和下子波形特征集F3；
(8)利用测试数据xG的各子波形集的权重因子a1，a2和a3，将测试数据xG的上子波形特征集F1，中子波形特征集F2和下子波形特征集F3进行加权融合，得到测试特征集F；
(9)利用训练模板库H和被鉴定者的测试特征集F，对被鉴定者的身份进行识别。
本发明与现有技术相比具有以下优点：
第一，本发明充分挖掘并利用PPG信号单周期波形的特征，通过将被鉴定者PPG信号的单周期波形进行分段处理，利用余弦相似公式计算各段子波形的权重因子，并利用权重因子对各段子波形进行加权融合，利用获得的融合特征向量进行身份识别，提高了被鉴定者身份正确识别率。
第二，本发明利用鉴别式非负矩阵分解方法提取PPG信号单周期波形的主要特征，使得相同个体的单周期波形特征之间的差异性变小，不同个体单周期波形特征之间的差异性增大，从而提高了被鉴定者身份的正确识别率。
附图说明
图1为本发明的实现总流程图；
图2为MIMIC数据库的身份识别率结果图；
图3为MIMIC2数据库的身份识别率结果图；
图4为CapnoBase数据库的身份识别率结果图。
具体实施方式
下面结合附图对本发明的实施及效果作进一步详细描述。
参照图1，本发明的实现如下：
步骤1.采集PPG信号，得到训练数据库和测试数据。
采集M个人在规定时间段内的PPG信号，设采集的每个人的PPG信号采样点数为N，将采集的每个人的PPG信号作为一个行向量，构造一个M×N大小的矩阵，作为训练数据库S；再采集其中一人在另一时间段内的PPG信号，作为被鉴定者的测试数据，用符号xG表示，则测试数据xG是一个包含多个采样点的向量。
本发明以MIMIC数据库中的PPG信号作为实验数据，模拟从人体采集到的PPG信号，从MIMIC数据库中随机选取其中50个个体的PPG数据文件，读取每个人的PPG数据文件的前200秒的PPG信号，组成训练数据库S；再从该50个个体中随机读取其中一个人的PPG数据文件后200秒的PPG信号，作为被鉴定者的测试数据xG；MIMIC数据库中的PPG信号的采样频率f为125Hz，所以训练数据库S是一个50×25000大小的矩阵，测试数据xG是一个1×25000维的向量。
步骤2.对训练数据库S进行预处理，获取单周期波形数据库。
(2a)对训练数据库S每行PPG信号进行去噪处理，由去噪后的所有行PPG信号组成去噪后的训练数据库S1，其中，训练数据库S1中每行的PPG信号S1i由N个采样点组成，表示为：S1i＝{S1(i,j)|j∈[1,N]}，i∈[1,M]，S1(i,j)表示去噪后的训练数据库S1第i行第j列的采样点，N表示每人PPG信号的采样点数；
常用的去噪方法有低通滤波器、小波去噪、傅里叶分析等，本实例采用2014年胡广书编著的清华大学出版社出版的“现代信号处理教程第2版”的第12.4节“小波去噪”；
(2b)对去噪后训练数据库S1的每行PPG信号的每个采样点进行归一化处理，使归一化后的所有采样点的取值都在区间[0,1]之内，得到归一化后的训练数据库S2，其中，归一化后的训练数据库S2中第i行第j列的采样点S2(i,j)计算公式为：
(2c)对归一化后的训练数据库S2的每行PPG信号进行收缩期波峰检测，获取所有收缩期波峰的位置，去除第一个位置和最后一个位置，由剩余的所有位置组成每行PPG信号收缩期波峰位置的集合其中，loci表示第i行PPG信号的所有收缩峰位置集合，表示位置集合loci的第nl个元素，表示第i行PPG信号所有收缩峰的位置集合loci的元素个数；
PPG信号是周期性信号，PPG信号的一个周期包括两个波峰，即收缩期波峰和舒张期波峰，收缩期波峰的幅度值高于舒张期波峰的幅度值。常用的波峰检测方法有极大值检测法、差分阈值法、自适应阈值法等，本实例使用王黎，韩清鹏编著的2011年科学出版社出版的“人体生理信号的非线性分析方法”一书的第4.4.1节“P波波峰点的提取方法”；
(2d)波形分割，获取单周期波形；
以第i行PPG信号的位置集合loci中的所有元素为分割点，将相邻两个分割点之间的波形作为一个单周期波形，对归一化后的训练数据库S2的第i行PPG信号S2i进行波形分割，使PPG信号S2i分割后变成个单周期波形，组成第i类单周期集合：其中，表示第i类单周期集合的第nw个单周期，表示第i类单周期集合中的波形个数；
(2e)对所有类单周期集合中每一个单周期波形进行插值，使插值后的每个单周期波形的采样点数均为n，得到插值后的单周期波形为由第i类所有插值后的单周期波形组成插值后的第i类单周期波形集合其中，n为正整数，表示第i类的第nw个插值后的单周期波形；常见的插值方法主要有拉格朗日插值、分段线性插值和三次样条插值等，本实例采用但不限于三次样条插值方法；
(2f)去除差异性较大的单周期波形；
计算第i类所有插值后的单周期波形的平均周期波形，作为参考波形；将每个单周期波形作为//随机变量，计算第i类的每个单周期波形与参考波形的皮氏积矩相关系数；然后将相关系数与设定的阈值Th进行比较，若相关系数小于设定的阈值Th，则删除相应的单周期波形，否则保留相应的单周期波形；由保留下来的所有单周期组合成第i类去除杂波的单周期波形集合其中，Th为一正实数，表示第i类去除杂波的单周期波形集合中的第nc个单周期波形，表示第i类去除杂波的单周期波形集合中的单周期波形个数；
(2g)对第i类去除杂波的单周期波形集合中的所有单周期波形进行分组，每连续m个单周期波形划分为一组，剔除剩余的未能凑成组的波形，然后对每组中的所有单周期波形进行平均，得到平均周期波形；将每一个平均周期波形作为一个列向量，由所有类的所有平均周期波形组合成单周期平均波形数据库单周期平均波形数据库W是一个n×K大小的矩阵，其中，表示单周期平均波形数据库W中第nm个单周期波形，K表示单周期平均波形数据库W的所有单周期波形个数，ki表示第i类单周期平均波形的个数，m为大于零的整数，表示向下取整操作。
步骤3.对单周期平均波形数据库W的单周期波形进行分段处理。取单周期波形数据库W的第1行到第行的所有数据组成训练数据的上子波形集V1，取数据库W的第行到第行的所有数据组成训练数据的中子波形集V2，取数据库W的第行到第行的所有数据组成训练数据的下子波形集V3，其中， 表示实u×K矩阵空间Ru×K中全体非负矩阵构成的子集合，u表示单个子波形的采样点数。这里的分段处理也可以采用其他不均衡的分段方式。
步骤4.计算训练数据的上子波形集V1、中子波形集V2和下子波形集V3各自的相似因子，并根据相似因子，得到各子波形集的权重因子。
(4a)将两个向量夹角的余弦值作为两个子波形之间的相似度，设vi,θ和vi,χ分别为上子波形集V1的第i类的第θ个子波形和第χ个子波形，则上子波形集V1的第i类所有子波形之间的相似度simi计算公式如下：
其中，θ,χ∈[1,ki]，(·)T表示向量或矩阵的转置，||·||L2表示向量的L2范数；
(4b)按照步骤(4a)的计算公式，计算训练数据的上子波形集V1的每一类子波形的相似度，再利用下列公式，得到上子波形集V1的相似因子s1：
(4c)按照步骤(4a)-(4b)，求出训练数据的中子波形集V2的相似因子s2和下子波形集V3的相似因子s3；
(4d)利用训练数据的上子波形集的相似因子s1，中子波形集的相似因子s2和下子波形集的相似因子s3，按照下列公式，求得各子波形集的权重因子d1，d2和d3：
其中，d1，d2和d3分别为训练数据的上子波形集的权重因子，中子波形集权重因子和下子波形集权重因子，d1，d2和d3都是0到1之间的实数。
步骤5.获取基空间和训练模板库。
利用鉴别式非负矩阵分解DNMF方法分别对训练数据的上子波形集V1，中子波形集V2和下子波形集V3进行分解，得到上子波形集的基空间Z1，中子波形集的基空间Z2，下子波形集的基空间Z3和各子波形的子特征集H1，H2，H3；并对各子波形的子特征集H1，H2，H3进行融合，得到训练模板库H。
此处的鉴别式非负矩阵分解方法是采用的2006年Stefanos Zafeiriou,Anastasios Tefas等人在“IEEE TRANSACTIONS ON NEURAL NETWORKS”期刊发表的“Exploiting Discriminant Information in Nonnegative Matrix Factorization WithApplication to Frontal Face Verification”一文中描述的鉴别式非负矩阵分解方法DNMF。
(5a)利用鉴别式非负矩阵分解方法对训练数据的上子波形集V1进行分解：
(5a1)随机初始化基矩阵Z(0)和系数矩阵H(0)，使基矩阵Z(0)中的任意元素满足ρ∈[1,u]，φ∈[1,r]，系数矩阵H(0)的任意元素满足τ∈[1,K]，其中，是基矩阵Z(0)中的第ρ行第φ列元素，r表示分解维数，表示系数矩阵H(0)中的第φ行第τ列元素；
(5a2)根据如下公式，对基矩阵Z(t)中的元素进行更新：
首先，按照如下公式更新，得到中间变量值
然后，对中间变量值归一化处理得到
将作为本次对基矩阵Z(t)中元素的迭代更新结果，其中，为迭代t-1次的基矩阵Z(t-1)的第ρ行第φ列元素，t∈[1,iter]，iter为预先定义的最大迭代次数，为迭代t-1次的系数矩阵H(t-1)的第φ行第τ列元素，V1(ρ,τ)表示训练数据的上子波形集V1的第ρ行第τ列元素；
(5a3)根据步骤(5a2)得到的迭代t次后的基矩阵Z(t)，按如下迭代规则更新系数矩阵H(t)中的元素
其中，γ，δ分别为类内散度约束项和类间散度约束项的约束因子，μφ表示系数矩阵H(t-1)中所有列向量的均值向量μ中的第φ个元素；
(5a4)采用预定义的最大迭代次数iter作为停止迭代条件，当迭代次数t达到iter次后，停止迭代，输出基矩阵Z(iter)和系数矩阵H(iter)；否则，返回步骤(5a2)；
(5b)将基矩阵Z(iter)作为上子波形集V1的基空间Z1，将系数矩阵H(iter)的每列作为一个子特征向量，组成上子波形的子特征集其中，表示实u×r矩阵空间Ru×r中全体非负矩阵构成的子集合，表示上子波形的子特征集H1中第fe个子特征向量，fe＝1,2,…,K，表示实r维向量空间Rr中全体非负向量构成的子集合；
(5c)按照步骤(5a)，分别对训练数据的中子波形集V2和下子波形集V3进行分解，得到中子波形集的基空间Z2和子特征集下子波形集的基空间Z3和子特征集其中，表示中子波形的子特征集H2中第fe个子特征向量；表示下子波形的子特征集H3中第fe个子特征向量；
(5d)利用训练数据的上子波形集的权重因子d1，中子波形集的权重因子d2和下子波形集的权重因子d3对上子波形的子特征集H1中的子特征向量中子波形的子特征集H2中的子特征向量和下子波形的子特征集H3中的子特征向量进行加权串联，得到融合特征模板然后由所有的融合特征模板组成训练模板库H＝{h1,h2,…,hfe,…,hK}，其中，hfe表示训练模板库H中的第fe个融合特征模板，表示实3r维向量空间R3r中全体非负向量构成的子集合。
步骤6.对被鉴定者的测试数据xG进行处理，得到测试数据xG的各子波形集及各子波形集的权重因子。
(6a)对测试数据xG进行步骤2-3操作，得到测试数据xG的上子波形集Α1，中子波形集Α2和下子波形集A3；其中，表示实u×G矩阵空间Ru×G中全体非负矩阵构成的子集合，G表示测试数据子波形集的子波形个数；
(6b)对测试数据xG的上子波形集Α1，中子波形集Α2和下子波形集A3进行步骤4操作，得到测试数据上子波形集的权重因子a1，中子波形集的权重因子a2，下子波形集的权重因子a3，其中，a1,a2,a3∈(0,1)。
步骤7.按照下式，将测试数据xG的上子波形集Α1，中子波形集Α2和下子波形集A3分别在训练数据的上子波形集基空间Z1，中子波形集基空间Z2和下子波形集基空间Z3上进行投影，获得测试数据上子波形的子特征集F1，中子波形的子特征集F2和下子波形的子特征集F3：
F1＝inv((Z1)T×Z1)×(Z1)T×A1，
F2＝inv((Z2)T×Z2)×(Z2)T×A2，
F3＝inv((Z3)T×Z3)×(Z3)T×A3，
其中，表示实r×G矩阵空间Rr×G中全体非负矩阵构成的子集合，inv(·)表示矩阵求逆运算。
步骤8.获取测试特征集。按照步骤(5d)，利用权重因子a1，a2和a3，对F1，F2和F3的对应列加权融合，得到测试特征集F＝{ξ1,ξ2,…,ξg,…,ξG}，其中，ξg表示测试特征集F的第g个测试特征向量，g＝1,2,…,G。
步骤9.利用支持向量机SVM对被鉴定者的身份进行识别。
(9a)将训练模板库H中所有模板输入到支持向量机中进行训练，得出支持向量机模型；
(9b)将被鉴定者的测试特征集F中的所有测试特征向量依次输入到训练好的支持向量机模型中进行类别预测，得到预测类别信息；
(9c)根据所有测试特征向量的预测类别信息，分别统计各类中测试特征向量的个数，将测试特征向量个数最多的类预测为被鉴定者的身份。
本发明的效果可通过以下仿真做进一步说明。
1.仿真条件
本发明的仿真实验使用三个公开的PPG信号数据库MIMIC，MIMIC2和Capnobase数据库，模拟从人体采集到的PPG信号，仿真实验在Intel Pentium E58003.2GHz CPU、内存2GB的计算机上进行。
2.仿真内容
首先，分别从MIMIC数据库，MIMIC2数据库和Capnobase数据库中随机选取50个人，50个人和42个人的PPG信号，使用本发明分别对数据库中的每个人进行身份预测，计算每个人的识别率：
身份识别率＝类别预测正确的测试特征数目/被鉴定者的测试特征总数；
然后，取数据库中所有人身份识别率的平均值作为数据库的身份识别率，得到每个库身份识别率的结果图，如图2、图3和图4。
从图2、图3和图4可以看出，每个库的正确身份识别率均达到99.62％以上，充分说明了本发明的有效性和高识别率。
Claims (8)
1.基于人体PPG信号分段的身份识别方法，包括如下步骤：
(1)获取训练数据库和测试数据。采集M个人在规定时间段内的光电容积脉搏波PPG信号，组成训练数据库S；再采集其中一人在另一时间段内的PPG信号，作为被鉴定者的测试数据xG；
(2)对训练数据库S依次进行去噪，归一化处理，波峰检测，波形分割，插值，去除差异性大的波形和进行波形平均，得到单周期平均波形数据库W；
(3)将单周期平均波形数据库W中每一个单周期波形进行分段处理，得到训练数据的上子波形集V1、中子波形集V2和下子波形集V3，并计算各子波形集中同一个人两两子波形之间的相似度，得到训练数据上子波形集的权重因子d1、中子波形集的权重因子d2和下子波形集的权重因子d3；
(4)利用鉴别式非负矩阵分解DNMF方法分别对训练数据的上子波形集V1、中子波形集V2和下子波形集V3进行分解，得到训练数据上子波形集的基空间Z1、中子波形集的基空间Z2、下子波形集的基空间Z3、训练数据上子波形特征集H1、中子波形特征集H2和下子波形特征集H3；
(5)利用训练数据上子波形集的权重因子d1，中子波形集的权重因子d2和下子波形集的权重因子d3分别对训练数据上子波形特征集H1，中子波形特征集H2下子波形特征集H3中相应的子特征进行加权融合，得到训练模板库H；
(6)对被鉴定者的测试数据xG依次进行步骤(2)-(3)操作，得到测试数据xG的上子波形集Α1，中子波形集Α2，下子波形集Α3和测试数据上子波形集的权重因子a1，测试数据中子波形集的权重因子a2和测试数据下子波形集的权重因子a3；
(7)将测试数据xG的上子波形集Α1，中子波形集Α2，下子波形集Α3分别在训练数据上子波形集的基空间Z1，中子波形集的基空间Z2，下子波形集的基空间Z3上进行投影，获得测试数据上子波形特征集F1，中子波形特征集F2和下子波形特征集F3；
(8)利用测试数据xG的各子波形集的权重因子a1，a2和a3，将测试数据xG的上子波形特征集F1，中子波形特征集F2和下子波形特征集F3进行加权融合，得到测试特征集F；
(9)利用训练模板库H和被鉴定者的测试特征集F，对被鉴定者的身份进行识别。
2.根据权利要求1所述的方法，其中步骤(2)中的归一化处理，采用的计算公式为：其中S2(i,j)表示为归一化后的训练数据库S2中第i行第j列的采样点，i∈[1,M]，j∈[1,N]，N表示每人PPG信号的采样点数，S1(i,j)表示去噪后的训练数据库S1中第i行第j列的采样点，S1i表示去噪后的训练数据库S1中第i行所有采样点。
3.根据权利要求1所述的方法，其中步骤(2)中的波形分割，是以第i行PPG信号的位置集合loci中的所有元素为分割点，将相邻两个分割点之间的波形作为一个单周期波形，对训练数据库S2的第i行PPG信号S2i进行波形分割，使PPG信号S2i分割后变成个单周期波形，组成第i类单周期集合：其中，表示第i类单周期集合的第nw个单周期，表示第i类单周期集合中的波形个数，i＝1,2,…,M。
4.根据权利要求1所述的方法，其中步骤(3)中计算各子波形集中同一人两两子波形之间的相似度，是将子波形看作为一个向量，计算两个子波形向量夹角的余弦值作为两个子波形之间的相似度。
5.根据权利要求1所述的方法，其中步骤(4)中利用鉴别式非负矩阵分解方法DNMF分别对训练数据的上子波形集V1，中子波形集V2和下子波形集V3进行分解，具体步骤如下：
(4a)利用鉴别式非负矩阵分解方法对训练数据的上子波形集V1进行分解：
(4a1)随机初始化基矩阵Z(0)和系数矩阵H(0)，使基矩阵Z(0)中的任意元素满足ρ∈[1,u]，φ∈[1,r]，系数矩阵H(0)的任意元素满足τ∈[1,K]，其中，是基矩阵Z(0)中的第ρ行第φ列元素，表示系数矩阵H(0)中的第φ行第τ列元素，u表示单个子波形的采样点数，r表示分解维数，K表示单周期平均波形数据库W中的波形总数；
(4a2)根据如下公式，对基矩阵Z(t)中的元素进行更新：
首先，按照如下公式更新，得到中间变量值
<mrow>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&phi;</mi>
</mrow>
<mrow>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&phi;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>H</mi>
<mrow>
<mi>&phi;</mi>
<mo>,</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mfrac>
<msub>
<mi>V</mi>
<mrow>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>r</mi>
</munderover>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&eta;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>H</mi>
<mrow>
<mi>&eta;</mi>
<mo>,</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mfrac>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>H</mi>
<mrow>
<mi>&phi;</mi>
<mo>,</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
然后，对中间变量值归一化处理得到
<mrow>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&phi;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&phi;</mi>
</mrow>
<mrow>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>u</mi>
</munderover>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>&eta;</mi>
<mo>,</mo>
<mi>&phi;</mi>
</mrow>
<mrow>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
将作为本次对基矩阵Z(t)中元素的迭代更新结果，其中，为迭代t-1次的基矩阵Z(t-1)的第ρ行第φ列元素，t∈[1,iter]，iter为预先定义的最大迭代次数，为迭代t-1次的系数矩阵H(t-1)的第φ行第τ列元素，V1(ρ,τ)表示训练数据的上子波形集V1的第ρ行第τ列元素；
(4a3)根据步骤(4a2)得到的迭代t次后的基矩阵Z(t)，按如下迭代规则更新系数矩阵H(t)中的元素
<mrow>
<msubsup>
<mi>H</mi>
<mrow>
<mi>&phi;</mi>
<mo>,</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msqrt>
<mrow>
<msubsup>
<mi>T</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mn>4</mn>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&gamma;</mi>
<mo>-</mo>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&gamma;</mi>
<mo>+</mo>
<mn>2</mn>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msubsup>
<mi>H</mi>
<mrow>
<mi>&phi;</mi>
<mo>,</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>u</mi>
</munderover>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&phi;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mfrac>
<msub>
<mi>V</mi>
<mrow>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>r</mi>
</munderover>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>&rho;</mi>
<mo>,</mo>
<mi>&eta;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>H</mi>
<mrow>
<mi>&eta;</mi>
<mo>,</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&gamma;</mi>
<mo>-</mo>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&gamma;</mi>
<mo>+</mo>
<mn>2</mn>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&gamma;</mi>
<mo>+</mo>
<mn>2</mn>
<mi>&delta;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>h</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>&lambda;</mi>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>&NotEqual;</mo>
<mi>&eta;</mi>
</mrow>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
</munderover>
<msubsup>
<mi>H</mi>
<mrow>
<mi>&phi;</mi>
<mo>,</mo>
<mi>&lambda;</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&delta;&mu;</mi>
<mi>&phi;</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
</mrow>
其中，γ，δ分别为类内散度约束项和类间散度约束项的约束因子，ki表示第i类单周期平均波形的个数，μφ表示系数矩阵H(t-1)中所有列向量的均值向量μ中的第φ个元素；
(4a4)采用预定义的最大迭代次数iter作为停止迭代条件，当迭代次数t达到iter次后，停止迭代，输出基矩阵Z(iter)和系数矩阵H(iter)；否则，返回步骤(4a2)；
(4b)将基矩阵Z(iter)作为上子波形集V1的基空间Z1，将系数矩阵H(iter)的每列作为一个子特征向量，组成上子波形的子特征集其中， 表示实u×r矩阵空间Ru×r中全体非负矩阵构成的子集合，表示上子波形的子特征集H1中第fe个子特征向量，fe＝1,2,…,K， 表示实r维向量空间Rr中全体非负向量构成的子集合；
(4c)按照步骤(4a)，分别对训练数据的中子波形集V2和下子波形集V3进行分解，得到中子波形集的基空间Z2和子特征集下子波形集的基空间Z3和子特征集其中， 表示中子波形的子特征集H2中第fe个子特征向量； 表示下子波形的子特征集H3中第fe个子特征向量。
6.根据权利要求1所述的方法，其中步骤(5)中得到训练模板库H，是利用训练数据的上子波形集的权重因子d1，中子波形集的权重因子d2和下子波形集的权重因子d3对上子波形的子特征集H1中的子特征向量中子波形的子特征集H2中的子特征向量和下子波形的子特征集H3中的子特征向量进行加权串联，得到融合特征模板然后由所有的融合特征模板组成训练模板库H＝{h1,h2,…,hfe,…,hK}，其中，hfe表示训练模板库H中的第fe个融合特征模板，表示实r维向量空间Rr中全体非负向量构成的子集合，表示实3r维向量空间R3r中全体非负向量构成的子集合，K表示单周期平均波形数据库W中的波形总数。
7.根据权利要求1所述的方法，其中步骤(7)中获得测试数据上子波形特征集F1，中子波形特征集F2和下子波形特征集F3，是按照下列公式得出：
F1＝inv((Z1)T×Z1)×(Z1)T×A1，
F2＝inv((Z2)T×Z2)×(Z2)T×A2，
F3＝inv((Z3)T×Z3)×(Z3)T×A3，
其中， 表示实r×G矩阵空间Rr×G中全体非负矩阵构成的子集合，r表示分解维数，G表示测试数据子波形集的子波形个数，inv(·)表示矩阵求逆运算。
8.根据权利要求1所述的方法，其中步骤(9)中利用训练模板库H和被鉴定者的测试特征集F，对被鉴定者的身份进行识别，具体步骤如下：
(9a)将训练模板库H中所有模板输入到支持向量机中进行训练，得出支持向量机模型；
(9b)将被鉴定者的测试特征集F中的所有测试特征向量依次输入到训练好的支持向量机模型中进行类别预测，得到预测类别信息；
(9c)根据所有测试特征向量的预测类别信息，分别统计各类中测试特征向量的个数，将个数最多的类预测为被鉴定者的身份。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710195051.8A CN107088069B (zh) | 2017-03-29 | 2017-03-29 | 基于人体ppg信号分段的身份识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710195051.8A CN107088069B (zh) | 2017-03-29 | 2017-03-29 | 基于人体ppg信号分段的身份识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107088069A true CN107088069A (zh) | 2017-08-25 |
| CN107088069B CN107088069B (zh) | 2019-10-11 |
Family
ID=59649000
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710195051.8A Active CN107088069B (zh) | 2017-03-29 | 2017-03-29 | 基于人体ppg信号分段的身份识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107088069B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109557520A (zh) * | 2018-12-26 | 2019-04-02 | 苏州天必佑科技有限公司 | 一种基于多方法融合的人体微弱呼吸信号增强方法 |
| CN109692005A (zh) * | 2018-12-03 | 2019-04-30 | 南京邮电大学 | 基于ppg信号稀疏分解的身份识别方法 |
| CN109919050A (zh) * | 2019-02-21 | 2019-06-21 | 天津惊帆科技有限公司 | 身份识别方法和装置 |
| CN110123289A (zh) * | 2019-04-08 | 2019-08-16 | 清华大学深圳研究生院 | 一种基于脉搏波的生物识别方法及相关装置 |
| CN110811582A (zh) * | 2019-12-09 | 2020-02-21 | 北京明略软件系统有限公司 | 心率检测方法、装置和设备 |
| CN111222498A (zh) * | 2020-03-19 | 2020-06-02 | 桂林电子科技大学 | 一种基于光电容积脉搏波的身份识别方法 |
| CN111914655A (zh) * | 2020-07-06 | 2020-11-10 | 电子科技大学 | 一种基于显著信号子段提取的脉象分类方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5954664A (en) * | 1995-04-06 | 1999-09-21 | Seegobin; Ronald D. | Noninvasive system and method for identifying coronary disfunction utilizing electrocardiography derived data |
| US20060161065A1 (en) * | 2005-01-18 | 2006-07-20 | Heartlab, Inc. | Similarity scores for electrocardiography |
| CN103714281A (zh) * | 2013-12-12 | 2014-04-09 | 深圳先进技术研究院 | 一种基于心电信号的身份识别方法 |
| CN105286876A (zh) * | 2014-07-17 | 2016-02-03 | 原相科技股份有限公司 | 应用生理特征的个人化控制系统和其操作方法 |
| CN105740680A (zh) * | 2014-12-30 | 2016-07-06 | 三星电子株式会社 | 用于基于生物信号对用户进行认证的方法和设备 |
| CN105787420A (zh) * | 2014-12-24 | 2016-07-20 | 北京三星通信技术研究有限公司 | 用于生物认证的方法、装置以及生物认证系统 |
| CN106473750A (zh) * | 2016-10-08 | 2017-03-08 | 西安电子科技大学 | 基于光电容积脉搏波最佳周期波形的身份识别方法 |
-
2017
- 2017-03-29 CN CN201710195051.8A patent/CN107088069B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5954664A (en) * | 1995-04-06 | 1999-09-21 | Seegobin; Ronald D. | Noninvasive system and method for identifying coronary disfunction utilizing electrocardiography derived data |
| US20060161065A1 (en) * | 2005-01-18 | 2006-07-20 | Heartlab, Inc. | Similarity scores for electrocardiography |
| CN103714281A (zh) * | 2013-12-12 | 2014-04-09 | 深圳先进技术研究院 | 一种基于心电信号的身份识别方法 |
| CN105286876A (zh) * | 2014-07-17 | 2016-02-03 | 原相科技股份有限公司 | 应用生理特征的个人化控制系统和其操作方法 |
| CN105787420A (zh) * | 2014-12-24 | 2016-07-20 | 北京三星通信技术研究有限公司 | 用于生物认证的方法、装置以及生物认证系统 |
| CN105740680A (zh) * | 2014-12-30 | 2016-07-06 | 三星电子株式会社 | 用于基于生物信号对用户进行认证的方法和设备 |
| CN106473750A (zh) * | 2016-10-08 | 2017-03-08 | 西安电子科技大学 | 基于光电容积脉搏波最佳周期波形的身份识别方法 |
Non-Patent Citations (2)
| Title |
|---|
| GANG ZHENG,ET AL.: ""Study of Human Identificaiton by Electrocardiogram Waveform Morph"", 《BIOMETRIC RECOGNITION》 * |
| SHIH-CHIN FANG: "Human identification by quantifying similarity and dissimilarity in electrocardiogram phase space", 《PATTERN RECOGNITION》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109692005A (zh) * | 2018-12-03 | 2019-04-30 | 南京邮电大学 | 基于ppg信号稀疏分解的身份识别方法 |
| CN109557520A (zh) * | 2018-12-26 | 2019-04-02 | 苏州天必佑科技有限公司 | 一种基于多方法融合的人体微弱呼吸信号增强方法 |
| CN109919050A (zh) * | 2019-02-21 | 2019-06-21 | 天津惊帆科技有限公司 | 身份识别方法和装置 |
| CN109919050B (zh) * | 2019-02-21 | 2021-07-20 | 天津惊帆科技有限公司 | 身份识别方法和装置 |
| CN110123289A (zh) * | 2019-04-08 | 2019-08-16 | 清华大学深圳研究生院 | 一种基于脉搏波的生物识别方法及相关装置 |
| CN110811582A (zh) * | 2019-12-09 | 2020-02-21 | 北京明略软件系统有限公司 | 心率检测方法、装置和设备 |
| CN111222498A (zh) * | 2020-03-19 | 2020-06-02 | 桂林电子科技大学 | 一种基于光电容积脉搏波的身份识别方法 |
| CN111914655A (zh) * | 2020-07-06 | 2020-11-10 | 电子科技大学 | 一种基于显著信号子段提取的脉象分类方法 |
| CN111914655B (zh) * | 2020-07-06 | 2022-07-12 | 电子科技大学 | 一种基于显著信号子段提取的脉象分类方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107088069B (zh) | 2019-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107088069B (zh) | 基于人体ppg信号分段的身份识别方法 | |
| CN101866421B (zh) | 基于离散度约束非负稀疏编码的自然图像特征提取方法 | |
| Yuan et al. | Fingerprint liveness detection using an improved CNN with image scale equalization | |
| Hejazi et al. | ECG biometric authentication based on non-fiducial approach using kernel methods | |
| Rashidi et al. | Feature extraction based DCT on dynamic signature verification | |
| CN103268497B (zh) | 一种人脸姿态检测方法与在人脸识别中的应用 | |
| MageshKumar et al. | Gabor features and LDA based face recognition with ANN classifier | |
| Fong et al. | A biometric authentication model using hand gesture images | |
| CN107122643B (zh) | 基于ppg信号和呼吸信号特征融合的身份识别方法 | |
| Zhang | Off‐line signature verification and identification by pyramid histogram of oriented gradients | |
| CN103366182B (zh) | 基于全监督非负矩阵分解的人脸识别方法 | |
| Gutierrez et al. | Modular neural network integrator for human recognition from ear images | |
| Dong et al. | Feature extraction through contourlet subband clustering for texture classification | |
| Khalifa et al. | Wavelet, gabor filters and co-occurrence matrix for palmprint verification | |
| Chelali et al. | Face recognition using MLP and RBF neural network with Gabor and discrete wavelet transform characterization: a comparative study | |
| Pratama et al. | Face recognition for presence system by using residual networks-50 architecture | |
| Bodini et al. | Single sample face recognition by sparse recovery of deep-learned lda features | |
| Rezaei et al. | Persian signature verification using fully convolutional networks | |
| CN110929243B (zh) | 一种基于手机惯性传感器的行人身份识别方法 | |
| CN104573728A (zh) | 一种基于极端学习机的纹理分类方法 | |
| Boukhris et al. | Artificial human face recognition via Daubechies wavelet transform and SVM | |
| Zalasiński et al. | Fuzzy-genetic approach to identity verification using a handwritten signature | |
| Laskaris et al. | Use of random time-intervals (RTIs) generation for biometric verification | |
| Adedeji et al. | Comparative Analysis of Feature Selection Techniques For Fingerprint Recognition Based on Artificial Bee Colony and Teaching Learning Based Optimization | |
| Alsaidi et al. | Hierarchical sparse method with applications In vision and speech recognition |
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 |