CN106778631B - 一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法 - Google Patents
一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法 Download PDFInfo
- Publication number
- CN106778631B CN106778631B CN201611196794.9A CN201611196794A CN106778631B CN 106778631 B CN106778631 B CN 106778631B CN 201611196794 A CN201611196794 A CN 201611196794A CN 106778631 B CN106778631 B CN 106778631B
- Authority
- CN
- China
- Prior art keywords
- iris
- image
- edge
- analysis
- sample
- 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.)
- Active
Links
- 210000000554 iris Anatomy 0.000 title claims abstract description 238
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000001914 filtration Methods 0.000 title claims abstract description 23
- 238000007781 pre-processing Methods 0.000 title claims abstract description 15
- 238000013461 design Methods 0.000 title claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 57
- 238000011156 evaluation Methods 0.000 claims abstract description 38
- 230000008569 process Effects 0.000 claims abstract description 20
- 238000003708 edge detection Methods 0.000 claims abstract description 16
- 230000008859 change Effects 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000009499 grossing Methods 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 4
- 101100398248 Arabidopsis thaliana KIN10 gene Proteins 0.000 claims description 3
- 101000835860 Homo sapiens SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily B member 1 Proteins 0.000 claims description 3
- 102100025746 SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily B member 1 Human genes 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims 1
- 238000012854 evaluation process Methods 0.000 abstract description 2
- 230000002146 bilateral effect Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000000605 extraction Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
- G06V40/18—Eye characteristics, e.g. of the iris
- G06V40/193—Preprocessing; Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/285—Selection of pattern recognition techniques, e.g. of classifiers in a multi-classifier system
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
- G06V40/18—Eye characteristics, e.g. of the iris
- G06V40/197—Matching; Classification
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Multimedia (AREA)
- Human Computer Interaction (AREA)
- Ophthalmology & Optometry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法，此分类器的设计包括改进的伪造虹膜数据库的设计、基于MSE分析的虹膜图像评估、基于SNR分析的虹膜图像评估、边缘检测、基于改进的边缘锐度分析的虹膜图像评估流程。在虹膜识别领域创新性的采用了双边滤波、MSE分析以及SNR分析相结合的虹膜图像评估体系针对无关类别虹膜攻击、散焦虹膜图像、运动模糊虹膜图像、美瞳虹膜图像的灰度特征及噪声分布特征实现对其高效且稳定的滤除。进一步针对仿真度较高的合成虹膜的边缘锐度特征创新性的设计基于改进的边缘锐度分析的虹膜图像评估体系。最终实现虹膜识别预处理过程中伪造虹膜快速、稳定的滤除。
Description
技术领域
本发明涉及计算机图像领域，具体涉及基于中科院自动化所CAISA虹膜数据库和C++图像处理函数、图像中指定灰度提取、降噪以及处理结果的可视化显示。
背景技术
随着网络和信息化技术的飞速发展，信息安全问题日益凸显。软硬件技术的提升使计算机及智能电子设备得以普及和使用，加速了网络技术对社会生活各个方面的渗透。通过网络技术非法获取用户的个人信息的事情时有发生，防不胜防，这使得人们更多的关注个人信息安全。传统的身份验证手段依据密码口令、证件信息，但是，由于其易于丢失、窃取、遗忘、伪造的缺陷，弊端凸显，也愈来愈满足不了当今信息化社会的安全需求，跟不上人们安全观念的提升。这就需要更加可靠的身份识别手段提高信息安全级别，增加身份识别的难度。而虹膜特征以其唯一性、稳定性、可采集性、非侵犯性、防规避性等特性，使得虹膜特征识别技术成为当前身份认证技术领域的热点。随着虹膜识别的推广，大量的异质虹膜(伪造虹膜和低质虹膜)也不断充斥着市场，对虹膜识别系统的识别性能和分类性能提出了更高的要求。
目前主要的伪造虹膜攻击及其典型特征：
(1)无关类别虹膜攻击：非虹膜图像检测；
(2)打印虹膜攻击：高清彩色打印虹膜图像检测；
(3)彩色隐形眼镜：美瞳；
(4)合成虹膜：取本人真实虹膜图像的有效虹膜区域与另外其他人虹膜进行合成。
目前主要的低质虹膜的典型特征：
(1)过度噪声遮挡；
(2)光线明暗不稳导致的散焦模糊；
(3)运动引起的运动模糊。
目前市场上的伪造虹膜和低质虹膜鱼目混珠的现象愈演愈烈，使得异质虹膜分类器技术受到学术界与企业界的关注。本文基于目前市场上主流的伪造虹膜攻击及低质虹膜攻击，从图像质量评估角度，摆脱了传统异质虹膜分类器依赖虹膜纹理特征进行分类的特性，设计了一种在虹膜识别系统实现定位之前即可滤除绝大多数异质虹膜的快速分类器。
发明内容
本发明的目的是针对传统的基于虹膜纹理特征的异质虹膜分类器需要在特征提取环节才能进行异质虹膜的滤除操作导致整个虹膜预处理、定位与分割及归一化操作环节几乎作废，因而浪费了大量的计算资源同时也降低虹膜识别系统的整体执行效率的问题。设计了基于MSE分析的虹膜图像评估、基于SNR分析的虹膜图像评估和基于改进的边缘锐度分析的虹膜图像评估三种评估方式为核心的异质虹膜分类器。该分类器分析了目前市场上主流的伪造虹膜攻击方式以及低质虹膜的攻击方式，有针对性的采用了上述的评估方式，可以在虹膜预处理环节过滤掉绝大多数的伪造及低质虹膜攻击，节约了计算机的计算资源，有效提升了效率。实验结果表明，在相同硬件环境、相同的综合虹膜样本集合中传统的虹膜分类器平均耗时875.33毫秒，而本文设计的异质虹膜分类器平均耗时433.79毫秒；传统分类器异质虹膜接收率为7％，而本文的分类器异质虹膜接收率为2％。因此，相比传统的基于虹膜纹理特征分类器，本文设计的分类器过滤效果更加全面。
本发明采用的技术方案如下：一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计，包括基于SNR分析的虹膜图像评估、基于MSE分析的虹膜图像评估和基于改进的边缘锐度分析的虹膜图像评估流程；具体的步骤如下：
步骤1：SNR分析的虹膜图像评估；
引入虹膜样本图像，生成虹膜图像灰度直方图，提取标准虹膜样本和伪造虹膜样本的图像像素灰度值，分别对其进行高斯平滑，再次计算标准虹膜样本和伪造虹膜样本的图像像素灰度值，带入SNR公式计算图像信噪比值。
步骤2：基于MSE分析的虹膜图像评估；
引入通过基于SNR分析的图像评估的虹膜样本图像，获取SNR图像评估获取的标准虹膜样本图像和伪造虹膜样本图像的像素灰度值，带入MSE均方误差计算公式，通过原始图像的灰度信息，计算出伪造虹膜图像相对于原始图像的失真程度。
步骤3：基于改进的边缘锐度分析的虹膜图像评估；
引入通过基于MSE分析的图像评估的虹膜样本图像，分析统计出图像某一边缘法向方向的灰度变化情况对图像进行评价。首先利用Sobel边缘检测算子提取出边缘，提取出边缘两侧像素灰度值，求出边缘法向的灰度变化率；提取边缘两侧总体像素灰度，求出该方向上的总体灰度变化，带入改进的边缘锐度分析公式计算出点锐度。
进一步，所述步骤1中，生成引入的虹膜样本图像的灰度直方图，提取标准虹膜样本和伪造虹膜样本的图像像素灰度值分别记为α(i,j)和β(i,j)；对标准虹膜样本和伪造虹膜样本分别进行高斯平滑降噪，再提取标准虹膜样本和伪造虹膜样本的图像像素灰度值分别记为α′(i,j)和β′(i，j)，带入图像信噪比的计算公式。如式(1)、式(2)所示。
式(1)、式(2)中的M和N分别代表图像长度和宽度上的像素个数。求出的信噪比值为SNR1和SNR2，图像信噪比值越大，图像质量越高。而打印虹膜和无关类别虹膜经过人为处理，噪声相对标准虹膜图像较多，基于SNR分析可以有效的实现滤除。
进一步，所述步骤2的具体过程为：引入通过基于SNR分析的图像评估的虹膜样本图像，获取SNR图像评估获取的标准虹膜样本图像和伪造虹膜样本图像的像素灰度值α(i,j)和β(i,j)，带入MSE均方误差计算公式如式(3)所示。
式(3)中的M和N分别代表图像长度和宽度上的像素个数。根据标准虹膜样本图像与伪造虹膜样本图像的灰度值，计算出MSE均方误差值，该值反映出伪造虹膜图像相对于原始图像的失真程度，图像失真越大(灰度变化较大)均方误差值越高，超过一定的误差范围即可认定为伪造虹膜。该方法有效的针对彩色隐形眼镜、散焦模糊、运动模糊等灰度变化较大的异质虹膜过滤效果很好，而且该方法相对于传统的纹理特征分析计算简单，执行效率高。
进一步，所选步骤(3)的具体过程为：引入通过基于MSE分析的图像评估的虹膜样本图像，分析统计出图像某一边缘法向方向的灰度变化情况对图像进行评价。首先利用Sobel边缘检测算子提取出边缘，提取出边缘两侧像素灰度值，求出边缘法向的灰度变化率记为df/dx；提取边缘两侧总体像素灰度，求出该方向上的总体灰度变化记为f(b)-f(a)。
由于传统的锐度分析算法只对特定边缘区域统计，不能代表整个图像的锐度，此外，进行边缘锐度计算前需要人为确定边缘，不方便实现系统的自动化。本发明在传统的锐度分析算法的基础上做出改进：
1)使用逐个像素邻域梯度计算代替了传统算法中边缘梯度计算，便于图像的整体评估。
2)对方格像素8个邻域灰度变化进行距离加权，0°与90°相邻的像素权值为1，45°和90°相邻的像素权值为
将上述提取的参数带入改进的边缘锐度分析公式计算出点锐度，公式如式(4)所示。
式(4)可以描述为：逐个对图像中的每点取8邻域点与之相减，先求8个差值的加权和(权的大小取决于距离,距离远，则权小，如45°和135°方向的差值需乘以1/2)，再将所有点所得值相加除以像素总个数。
本发明的有益效果是：
(1)对比目前市场上主流的伪造虹膜和低质虹膜样本图像与标准虹膜样本图像，分析出几种主流的伪造、低质虹膜的特征，并在此基础上有针对性的采用基于SNR分析、基于MSE分析、基于边缘锐度分析的虹膜图像评估方式，按照阶梯式的过滤结构，实现异质虹膜的滤除。
(2)针对传统的基于虹膜纹理特征的异质虹膜分类器需要在特征提取环节才能进行异质虹膜的滤除操作导致整个虹膜预处理、定位与分割及归一化操作环节几乎作废，浪费大量的计算资源的问题，设计基于SNR、MSE、边缘锐度分析的图像评估方法为核心的分类器，成功在虹膜预处理环节将异质虹膜的滤除，有效提升了执行的效率。
(3)针对传统的异质虹膜分类器针对虹膜纹理特征进行分类导致过滤异质虹膜种类有限，无法对合成红膜进行有效滤出的问题做出改进。针对合成虹膜边缘锐度变化的特征，将锐度分析引入到异质虹膜图像评估中，提出基于改进边缘锐度分析的虹膜图像评估方法。
(4)针对传统的点锐度算法只对特定边缘区域统计，不能代表整个图像的锐度，且进行边缘锐度计算前需要人为确定边缘，不方便实现系统的自动化的问题进行改进。
1)使用逐个像素邻域梯度计算代替了传统算法中边缘梯度计算，便于图像的整体评估。
2)对方格像素8个邻域灰度变化进行距离加权，0°与90°相邻的像素权值为1,45°和90°相邻的像素权值为
(5)采用如4所述的基于改进锐度分析的图像评估方法结合Sobel边缘检测算子检测边缘，可以相互取长补短，Sobel算子简单快速，效率极高但未能将图像主体与背景严格分开(非常适用于采集方式严格，背景相对简单的虹膜图像)，由于本文所述分类器采用如1所述的阶梯式分类结构，SNR与MSE分析图像评估会将明显低质的虹膜去除，保证Sobel算子处理的虹膜样本图像相对稳定，一定程度上也保证了Sobel边缘检测的效果。因此Sobel算子结合基于改进锐度分析的图像评估方法能够稳定、快速地对合成虹膜进行滤除。
(6)为了能够进一步适应现代计算机硬件多核化发展，实现编码的体现多核并行运算的设计思想，改变CPU性能过剩的特征。从硬件角度配合算法优化，实现了虹膜分类器执行效率的进一步改进。
附图说明
图1是本发明的方法实现原理流程图。
图2为本发明的标准虹膜与异质虹膜高斯平滑处理图。
其中，(a)表示标准虹膜图像，(b)表示标准虹膜图像高斯平滑处理图，(c)表示合成虹膜图像，(d)表示合成虹膜高斯平滑处理图。
图3为本发明的标准虹膜与异质虹膜Sobel边缘检测图。
其中，(a)表示，标准虹膜图像，(b)表示，标准虹膜图像Sobel边缘检测图，(c)表示，异质虹膜图像，(d)表示，异质虹膜图像Sobel边缘检测图。
图4为本发明的标准虹膜与异质虹膜SNR分析图。
图5为本发明的标准虹膜与异质虹膜MSE分析图。
图6为本发明的标准虹膜与异质虹膜点锐度分析图。
图7为本发明的异质虹膜分类器整体设计模型图。
具体实施方式
以下结合附图对发明的实施做进一步说明。
如图1所示，本发明的一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计，包括以下步骤实现：
步骤1：SNR分析的虹膜图像评估；
引入虹膜样本图像，生成虹膜图像灰度直方图，提取标准虹膜样本和伪造虹膜样本的图像像素灰度值，分别对其进行高斯平滑，再次计算标准虹膜样本和伪造虹膜样本的图像像素灰度值，带入SNR公式计算图像信噪比值。
步骤2：基于均方误差分析MSE的虹膜图像评估；
引入通过基于SNR分析的图像评估的虹膜样本图像，获取SNR图像评估获取的标准虹膜样本图像和伪造虹膜样本图像的像素灰度值，带入MSE均方误差计算公式，通过原始图像的灰度信息，计算出伪造虹膜图像相对于原始图像的失真程度。
步骤3：基于改进的边缘锐度分析的虹膜图像评估；
引入通过基于MSE分析的图像评估的虹膜样本图像，分析统计出图像某一边缘法向方向的灰度变化情况对图像进行评价。首先利用Sobel边缘检测算子提取出边缘，提取出边缘两侧像素灰度值，求出边缘法向的灰度变化率；提取边缘两侧总体像素灰度，求出该方向上的总体灰度变化，带入改进的边缘锐度分析公式计算出点锐度。
所述步骤1的具体过程为：生成引入的虹膜样本图像的灰度直方图，提取标准虹膜样本和伪造虹膜样本的图像像素灰度值分别记为α(i,j)和β(i,j)；对标准虹膜样本和伪造虹膜样本分别进行高斯平滑降噪，如图2所示，再提取标准虹膜样本和伪造虹膜样本的图像像素灰度值分别记为α′(i,j)和β′(i,j)，带入图像信噪比的计算公式如式(1)、式(2)所示。
式(1)、式(2)中的M和N分别代表图像长度和宽度上的像素个数。求出的信噪比值为SNR1和SNR2，如图4所示，图像信噪比值越大，图像质量越高。而打印虹膜和无关类别虹膜经过人为处理，噪声相对标准虹膜图像较多，基于SNR分析可以有效的实现滤除。
所述步骤2的具体过程为：引入通过基于SNR分析的图像评估的虹膜样本图像，获取SNR图像评估获取的标准虹膜样本图像和伪造虹膜样本图像的像素灰度值α(i,j)和β(i,j)，带入MSE均方误差计算公式如式(3)所示。
式(3)中的M和N分别代表图像长度和宽度上的像素个数。根据标准虹膜样本图像与伪造虹膜样本图像的灰度值，计算出MSE均方误差值，该值反映出出伪造虹膜图像相对于原始图像的失真程度，图像失真越大(灰度变化较大)均方误差值越高，如图5所示，超过一定的误差范围即可认定为伪造虹膜。该方法有效的针对彩色隐形眼镜、散焦模糊、运动模糊等灰度变化较大的异质虹膜过滤效果很好，而且该方法相对于传统的纹理特征分析计算简单，执行效率高。本发明采用中科院的虹膜CASIA V1.0版本的虹膜图像样本，阈值设定为83，针对不同的虹膜图像数据库，此阈值的划分也是有区别。
所选步骤3的具体过程为：引入通过基于MSE分析的图像评估的虹膜样本图像，分析统计出图像某一边缘法向方向的灰度变化情况对图像进行评价。首先利用Sobel边缘检测算子提取出边缘，提取出边缘两侧像素灰度值，求出边缘法向的灰度变化率记为df/dx；提取边缘两侧总体像素灰度，求出该方向上的总体灰度变化记为f(b)-f(a)。
Sobel边缘测算子提取边缘过程如下：
滑动图像窗口与Sobel边缘检测算子模板卷积，求得当前窗口中心像素的梯度值。不断滑动卷积窗口求得完整虹膜样本图像像素点的梯度值。生成原图的梯度图像，可以有效的凸显出图像的边缘部分。该发明中Sobel算子卷积模板表示为：
其中Sx表示为水平方向的卷积模板，Sy表示为垂直方向的卷积模板，K[i,j]表示为滤波窗口，ai表示像素点，i＝1、2、3…。图像像素点梯度值计算公式如式(7)：
其中Gx[i,j]代表水平方向梯度，Gy[i,j]表示垂直方向梯度，G[i,j]表示总的梯度。*代表卷积，h、v分别代表水平、垂直方向梯度权值，取值0或1。该发明中内边缘检测参数(h＝1，v＝1)效果如图5所示，外边缘检测参数(h＝0,v＝1)效果如图3所示。
由于传统的锐度分析算法只对特定边缘区域统计，不能代表整个图像的锐度，此外，进行边缘锐度计算前需要人为确定边缘，不方便实现系统的自动化。本发明在传统的锐度分析算法的基础上做出改进：
a)使用逐个像素邻域梯度计算代替了传统算法中边缘梯度计算，便于图像的整体评估；
b)对方格像素8个邻域灰度变化进行距离加权，0°与90°相邻的像素权值为1,45°和90°相邻的像素权值为
按照如上a)、b)的方法改进传统的锐度分析算法便于图像的整体评估，算法的自动化实现；同时对邻域像素的灰度加权，更好的表现出图像特征。
将上述a)、b)中提取的参数带入改进的边缘锐度分析公式计算出点锐度，公式如式(8)所示。
式(8)可以描述为：逐个对图像中的每点取8邻域点与之相减，先求8个差值的加权和(权的大小取决于距离，距离远，则权小，如45°和135°方向的差值需乘以
)，再将所有点所得值相加除以像素总个数，类似于对图像点扩展函数能量分布情况的统计。合成虹膜由于需要你对虹膜区域分割，因此往往边缘锐度受损，因此可以通过对比标准虹膜图像与异质虹膜图像的边缘点锐度，区分标准虹膜与异质虹膜，结果如图6所示。
最后发明的异质虹膜分类器整体设计模型图如图7所示。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。
Claims (6)
1.一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法，其特征在于，包括如下步骤：
步骤1：SNR分析的虹膜图像评估；
引入虹膜样本图像，生成虹膜图像灰度直方图，提取标准虹膜样本和伪造虹膜样本的图像像素灰度值，分别对其进行高斯平滑，再次计算标准虹膜样本和伪造虹膜样本的图像像素灰度值，代入SNR公式计算图像信噪比SNR的值；
步骤2：基于MSE分析的虹膜图像评估；
引入通过基于SNR分析的图像评估的虹膜样本图像，获取SNR图像评估获取的标准虹膜样本图像和伪造虹膜样本图像的像素灰度值，代入MSE均方误差计算公式，通过原始图像的灰度信息，计算出伪造虹膜图像相对于原始图像的失真程度；
步骤3：基于改进的边缘锐度分析的虹膜图像评估；
引入通过基于MSE分析的图像评估的虹膜样本图像，分析统计出图像某一边缘法向方向的灰度变化情况对图像进行评价；首先利用Sobel边缘检测算子提取出边缘，提取出边缘两侧像素灰度值，求出边缘法向的灰度变化率；提取边缘两侧总体像素灰度，求出该方向上的总体灰度变化，代入改进的边缘锐度分析算法公式计算出点锐度；
所述步骤1的实现包括：生成引入的虹膜样本图像的灰度直方图，提取标准虹膜样本和伪造虹膜样本的图像像素灰度值分别记为α(i,j)和β(i,j)；对标准虹膜样本和伪造虹膜样本分别进行高斯平滑降噪，再提取标准虹膜样本和伪造虹膜样本的图像像素灰度值分别记为α′(i,j)和β′(i,j)，代入图像信噪比的计算公式：
式中的M和N分别代表图像长度和宽度上的像素个数；
求出的信噪比值为SNR1和SNR2，图像信噪比值越大，图像质量越高；
所选步骤(3)的具体过程为：
引入通过基于MSE分析的图像评估的虹膜样本图像，分析统计出图像某一边缘法向方向的灰度变化情况对图像进行评价；首先利用Sobel边缘检测算子提取出边缘，提取出边缘两侧像素灰度值，求出边缘法向的灰度变化率记为df/dx；提取边缘两侧总体像素灰度，求出该方向上的总体灰度变化记为f(b)-f(a)；
所述Sobel边缘测算子提取边缘的具体过程如下：
滑动图像窗口与Sobel边缘检测算子模板卷积，求得当前窗口中心像素的梯度值；不断滑动卷积窗口求得完整虹膜样本图像像素点的梯度值，生成原图的梯度图像，有效凸显出图像的边缘部分；
所述锐度分析算法采用：a)使用逐个像素邻域梯度计算代替了传统算法中边缘梯度计算，便于图像的整体评估；b)对方格像素8个邻域灰度变化进行距离加权，0°与90°相邻的像素权值为1,45°和90°相邻的像素权值为
c)将上述提取的参数代入改进的边缘锐度分析公式计算出点锐度。
2.根据权利要求1所述的一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法，其特征在于，所述步骤2的具体过程为：
引入通过基于SNR分析的图像评估的虹膜样本图像，获取SNR图像评估获取的标准虹膜样本图像和伪造虹膜样本图像的像素灰度值α(i,j)和β(i,j)，代入MSE均方误差计算公式：
式中的M和N分别代表图像长度和宽度上的像素个数；
根据标准虹膜样本图像与伪造虹膜样本图像的灰度值，计算出MSE均方误差值，该值反映出伪造虹膜图像相对于原始图像的失真程度，图像失真越大，即灰度变化较大，均方误差值越高，超过误差范围即可认定为伪造虹膜。
3.根据权利要求1所述的一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法，其特征在于，所述Sobel算子卷积模板表示为：
其中Sx表示为水平方向的卷积模板，Sy表示为垂直方向的卷积模板，K[i,j]表示为滤波窗口，ai表示像素点。
4.根据权利要求1所述的一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法，其特征在于，所述图像像素点梯度值计算表达式为：
其中Gx[i,j]代表水平方向梯度，Gy[i,j]表示垂直方向梯度，G[i,j]表示总的梯度，*代表卷积，h、v分别代表水平、垂直方向梯度权值，Sx表示为水平方向的卷积模板，Sy表示为垂直方向的卷积模板，K[i,j]表示为滤波窗口。
5.根据权利要求4所述的一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法，其特征在于，所述h、所述v的取值为0或1。
6.根据权利要求2所述的一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法，其特征在于，所述误差范围设为83。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611196794.9A CN106778631B (zh) | 2016-12-22 | 2016-12-22 | 一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611196794.9A CN106778631B (zh) | 2016-12-22 | 2016-12-22 | 一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106778631A CN106778631A (zh) | 2017-05-31 |
| CN106778631B true CN106778631B (zh) | 2020-11-20 |
Family
ID=58899306
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611196794.9A Active CN106778631B (zh) | 2016-12-22 | 2016-12-22 | 一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106778631B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108470170B (zh) * | 2018-07-25 | 2018-11-02 | 上海聚虹光电科技有限公司 | 美瞳检测方法 |
| CN110619272A (zh) * | 2019-08-14 | 2019-12-27 | 中山市奥珀金属制品有限公司 | 一种虹膜图像分割方法 |
| CN112528720B (zh) * | 2020-04-03 | 2024-04-05 | 西安钗瑞信息科技有限公司 | 一种基于深度学习的红外体温测量系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201893A (zh) * | 2006-09-30 | 2008-06-18 | 电子科技大学中山学院 | 一种基于灰度信息的虹膜识别预处理方法 |
| CN101923640A (zh) * | 2010-08-04 | 2010-12-22 | 中国科学院自动化研究所 | 基于鲁棒纹理特征和机器学习对伪造虹膜图像判别的方法 |
| CN102567996A (zh) * | 2012-01-05 | 2012-07-11 | 中国地质大学（武汉） | 一种遥感卫星图像质量评价的方法 |
| CN103824293A (zh) * | 2014-02-28 | 2014-05-28 | 北京中科虹霸科技有限公司 | 虹膜采集设备成像质量评测系统 |
-
2016
- 2016-12-22 CN CN201611196794.9A patent/CN106778631B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201893A (zh) * | 2006-09-30 | 2008-06-18 | 电子科技大学中山学院 | 一种基于灰度信息的虹膜识别预处理方法 |
| CN101923640A (zh) * | 2010-08-04 | 2010-12-22 | 中国科学院自动化研究所 | 基于鲁棒纹理特征和机器学习对伪造虹膜图像判别的方法 |
| CN102567996A (zh) * | 2012-01-05 | 2012-07-11 | 中国地质大学（武汉） | 一种遥感卫星图像质量评价的方法 |
| CN103824293A (zh) * | 2014-02-28 | 2014-05-28 | 北京中科虹霸科技有限公司 | 虹膜采集设备成像质量评测系统 |
Non-Patent Citations (1)
| Title |
|---|
| 基于稀疏表示的伪造虹膜检测算法研究;韩晋;《中国优秀硕士学位论文全文数据库信息科技辑（月刊）》;20160415;第21-33页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106778631A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106407917B (zh) | 基于动态尺度分配的视网膜血管提取方法及系统 | |
| EP3455782B1 (en) | System and method for detecting plant diseases | |
| Raghunandan et al. | Riesz fractional based model for enhancing license plate detection and recognition | |
| CN110334706B (zh) | 一种图像目标识别方法及装置 | |
| US10558841B2 (en) | Method and apparatus for recognizing fingerprint ridge point | |
| Savkare et al. | Automatic system for classification of erythrocytes infected with malaria and identification of parasite's life stage | |
| CN107909081B (zh) | 一种深度学习中图像数据集的快速获取和快速标定方法 | |
| Liang et al. | Multi-spectral fusion based approach for arbitrarily oriented scene text detection in video images | |
| CN100373397C (zh) | 一种虹膜图像预处理方法 | |
| CN101201893A (zh) | 一种基于灰度信息的虹膜识别预处理方法 | |
| CN109636824A (zh) | 一种基于图像识别技术的多目标计数方法 | |
| CN108765465A (zh) | 一种无监督sar图像变化检测方法 | |
| CN110120056A (zh) | 基于自适应直方图阈值和轮廓检测的血液白细胞分割方法 | |
| CN110188614A (zh) | 一种基于皮裂纹分割的nlm滤波指静脉去噪方法 | |
| Roy et al. | Fractional poisson enhancement model for text detection and recognition in video frames | |
| CN110659649A (zh) | 一种基于近红外光成像的图像处理与字符识别算法 | |
| Pamplona Segundo et al. | Pore-based ridge reconstruction for fingerprint recognition | |
| CN111611907A (zh) | 一种图像增强的红外目标检测方法 | |
| WO2021057395A1 (zh) | 一种鞋跟型号识别方法、装置及存储介质 | |
| CN106778631B (zh) | 一种虹膜识别预处理过程中快速滤除伪造虹膜的异质虹膜分类器设计方法 | |
| CN110222647B (zh) | 一种基于卷积神经网络的人脸活体检测方法 | |
| WO2013075295A1 (zh) | 低分辨率视频的服装识别方法及系统 | |
| CN107038445A (zh) | 一种针对中文字符验证码的二值化和分割方法 | |
| CN117649606A (zh) | 一种基于本征表示模型的高光谱图像阴影去除方法与系统 | |
| CN116524269A (zh) | 一种视觉识别检测系统 |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210107 Address after: Room 508, building 3, 99 business Avenue, Huaqiao Town, Kunshan City, Suzhou City, Jiangsu Province Patentee after: JIANGSU START DIMA DATA PROCESSING Co.,Ltd. Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301 Patentee before: JIANGSU University |
|
| TR01 | Transfer of patent right |