CN105488494B - 一种虹膜内边缘的快速精确的定位方法 - Google Patents
一种虹膜内边缘的快速精确的定位方法 Download PDFInfo
- Publication number
- CN105488494B CN105488494B CN201511006255.XA CN201511006255A CN105488494B CN 105488494 B CN105488494 B CN 105488494B CN 201511006255 A CN201511006255 A CN 201511006255A CN 105488494 B CN105488494 B CN 105488494B
- Authority
- CN
- China
- Prior art keywords
- image
- edge
- iris
- inner boundary
- coefficient
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Ophthalmology & Optometry (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种虹膜内边缘的快速精确定位方法。现有的虹膜内边缘定位方法要么受到定位精度的限制，要么受到定位速度的限制，在或者需要预先满足许多前提条件的问题等情况下，很难满足图像生物识别技术的需求。本发明对采集的虹膜图像先进行边缘点提取；然后对前面提取的两个边缘图像进行边缘强度自适应融合，获得最终边缘图像融合图；之后对融合的新边缘图深度进行去除多余的杂边缘点和干扰点；最后再进行重新确定边缘点，将其进行数据存储以及显示于虹膜原图中。本发明对虹膜内边缘的定位在速度和精确度上都有很大程度的提高，能够满足图像生物识别技术的需求。
Description
技术领域
本发明属于图像生物识别技术领域，涉及一种虹膜内边缘的快速精确的定位方法。
背景技术
随着信息技术、网络技术的快速发展，以及网络信息化时代的一大特征就是身份的数字化和隐性化，且传统的身份鉴别方法是基于鉴别一些标识个人身份的事物，主要包括两个方面：(1).身份标识物品；(2).身份标识知识。这些传统的身份识别存在着严重的缺陷，使人类安全面临重大的挑战。而虹膜识别是当前应用最为方便和精确的一种。因此，如何准确鉴定一个人的身份，保护信息安全是当今信息化时代必须解决的一个关键社会问题，具有十分重大的价值与市场应用前景非常广阔。
虹膜内边缘的准确快速定位是虹膜定位的重要关键一部分，而虹膜定位是虹膜识别技术中的关键基础环节，对后期的特征提取和模式匹配步骤有很大的作用。虹膜边缘定位发展至今，产生了许多的定位算法，但绝大多数算法要么受到定位精度的限制，要么受到定位速度的限制，在或者需要预先满足许多前提条件的问题等，都存在着一些缺陷，而被产业界广泛接受的始终还是那些简单、高效且无需过多前提条件的虹膜定位方法。所以，研究一种内边缘的简单快速精确的定位算法仍然是一个富有挑战性的课题。
发明内容
本发明针对现有技术的不足，提出了一种虹膜内边缘的快速精确的定位方法。
本发明解决技术问题所采取的技术方案步骤如下:
步骤(1).虹膜图像采集。
步骤(2).对提取的虹膜灰度图像进行边缘点提取，这里用的是基于改进Canny算子进行检测边缘。
步骤(3).对提取的虹膜灰度图像利用小波变换提取边缘局部有效信息。
步骤(4).对步骤(2)和步骤(3)提取的边缘图像进行边缘强度自适应融合，获得最终边缘图像融合图。
步骤(5).对融合的新边缘图深度进行去除多余的杂边缘点和干扰点，这里用的是基于边缘信息的四向扫描法提取目标物。
步骤(6).之后再进行重新确定边缘点，将其进行数据存储以及显示于虹膜原图中。
本发明的有益效果：
(1).本发明不仅降低算法的运算量，而且避免了不必要的纹理细节，从而使后期的图像分割结果更快且更准确。
(2).四向区域扫描法是一种基于边缘几何信息的图像分割，所有参与的边缘点都可完全独立完成算法，无需与其他的边缘点进行信息交互，具备很好的并行特性，对算法的速度有很大程度的提升。
(3).可以通过对整幅矩形图像的四条边界线选择性去除和添加，以不同颜色区分，而且扫描长度参数是可控的，可筛选出不同面积的闭合区域；这样以利于瞳孔边缘的准确检测与定位。
附图说明
图1为本发明的主流程图；
图2为本发明步骤(2)的流程图；
图3为本发明步骤(5)四向区域填充算法的流程图。
具体的实施方式：
以下结合附图对本发明作进一步的描述。
如图1所示，本发明具体的技术方案为：
步骤(1).虹膜图像采集。由于虹膜的区分主要在于纹理细节的不同，因此首要的任务是获取高质量的虹膜图像，但用普通的CCD摄像头在正常光照下是很难获得清晰的虹膜图像且现有的虹膜取像装置仍价格昂贵，这里使用的是中科院提供的虹膜库，里面包含四万七千多张各式各样条件下的虹膜图，很具有研究价值。
步骤(2).对中科院提供的虹膜灰度图像进行边缘点提取，这里用的是基于改进Canny算子进行检测边缘。由于虹膜图像的获取过程是着力于高质量的图像获取，因此，从信息理论的角度来看，最好的预处理是没有预处理。对采集到的虹膜图像直接进行分析处理。
首先，根据几万张的虹膜图片实验检测结果分析，设置它的空间尺度系数和高阈值，以控制图像平滑程度和筛除次要边缘；然后根据图2做进一步的说明，具体如下：
第一步：用高斯滤波器平滑图像；
其高斯平滑函数为：原图像与高斯平滑函数卷积：H(x,y)＝f(x,y)*G(x,y)；
第二步：先对高斯平滑函数进行求一阶偏导，分别为：
用一阶偏导与图像卷积，然后计算每个像素点的梯度幅值和梯度方向，其分别为：P(x,y)＝f(x,y)*Gx，Q(x,y)＝f(x,y)*Gy，
θ(x,y)＝arctan(Q(x,y)/P(x,y))；
第三步：对梯度幅值进行非极大值抑制；
第四步：进行高阈值边缘检测生成强边缘图。
步骤(3).由于小波变换具有多分辨率分析的特点，在时频两域都具有表征信号局部特征的能力，且它的窗口面积固定不变，而窗口长度及宽度可以自动变化，能根据需要对信号高、低频部分具体分析，所以这里使用小波变换进行提取虹膜图像边缘的局部特征，然后与步骤(2)中的强边缘图进行信息融合，目的是为了尽量弥补瞳孔的遗漏边缘点，使瞳孔的边缘是一个完整的闭合区域。
对于尺度参数j,平移参数k,将yj,k(t)中的连续变量j和k取做整数离散形式，将其表示为：
yj,k(t)＝2j/2y(2jt-k)
得到离散小波，将其表示为：
Wf(j,k)＝(f(t),yj,k(t))
小波变换通过参数j、k的调节，从而实现在频域与时域方面的局部化作用。
步骤(4).由于边缘图像只保留了图像的边缘信息，而且为了使融合后图像的边缘效果达到最佳，选取基于小波变换的窗口内边缘强度自适应融合算法。
基本思想：选定窗口，计算窗内各像素点的边缘强度，并以归一化的边缘强度作为权值对两幅图像的高频系数进行加权求和。
设F1(x,y)、F2(x,y)分别为待融合的两幅图像，F(x,y)为融合后的图像。D1(x,y)、D2(x,y)分别为两幅图像矩阵中点(x,y)处的像素值，D(x,y)为融合后图像点(x,y)处的像素值。为第m幅图像在尺度系数为j、方向系数为k时，点(x,y)处经小波分解得到的系数；为第m幅图像在尺度系数为j、方向系数为k时，以(x,y)点为中心的窗口区域中小波系数的均值。W为窗口区域，设所选区域是边长为a的方形区域，则W＝a*a。
点(x,y)处的边缘强度定义为：
两幅图像中，的权系数取为：
将小波边缘图与强边缘图融合后，得到的最终边缘图像用F(x,y)来表示，则各像素点的像素值D(x,y)表示为：
步骤(5).对获得的融合边缘图进行深度彻底的去除多余的杂边缘点和干扰点，这里用的是基于边缘信息的四向扫描法。根据图3所示。
用四个方向(从上到下，从左到右，从下到上，从右到左)以不同颜色进行扫描填充，首先设置扫描长度参数和颜色参数，再是进行第一次的从上到下的区域扫描填充，根据图3流程依次进行，其中上一次的填充色将成为下一次的扫描检测色，最后一个方向的填充色将是唯一的保留色，即分割出了所需要的提取目标物。
步骤(6).之后再对提取的目标物进行重新确定边缘点，将其进行数据存储以及显示于虹膜原图中。
Claims (3)
1.一种虹膜内边缘的快速精确的定位方法，其特征在于该方法包括以下步骤：
步骤(1).虹膜图像采集；
步骤(2).对提取的虹膜灰度图像进行边缘点提取；
步骤(3).对提取的虹膜灰度图像利用小波变换提取边缘局部有效信息；
步骤(4).对步骤(2)和步骤(3)提取的边缘图像进行边缘强度自适应融合，获得最终边缘图像融合图；
步骤(5).对融合的新边缘图深度进行去除多余的杂边缘点和干扰点；
步骤(6).之后再进行重新确定边缘点，将其进行数据存储以及显示于虹膜原图中；
步骤(2)中使用了改进Canny算子进行检测边缘，具体是：
首先，根据几万张的虹膜图片实验检测结果分析，设置它的空间尺度系数和高阈值，以控制图像平滑程度和筛除次要边缘；然后如下：
第一步：用高斯滤波器平滑图像；
其高斯平滑函数为：原图像与高斯平滑函数卷积：H(x,y)＝f(x,y)*G(x,y)；
第二步：先对高斯平滑函数进行求一阶偏导，分别为：
用一阶偏导与图像卷积，然后计算每个像素点的梯度幅值和梯度方向，其分别为：P(x,y)＝f(x,y)*Gx，Q(x,y)＝f(x,y)*Gy，
第三步：对梯度幅值进行非极大值抑制；
第四步：进行高阈值边缘检测生成强边缘图；
步骤(4)具体是：
设F1(x,y)、F2(x,y)分别为待融合的两幅图像，F(x,y)为融合后的图像；D1(x,y)、D2(x,y)分别为两幅图像矩阵中点(x,y)处的像素值，D(x,y)为融合后图像点(x,y)处的像素值；为第m幅图像在尺度系数为j、方向系数为k时，点(x,y)处经小波分解得到的系数；为第m幅图像在尺度系数为j、方向系数为k时，以(x,y)点为中心的窗口区域中小波系数的均值；W为窗口区域，设所选区域是边长为a的方形区域，则W＝a*a；
点(x,y)处的边缘强度定义为：
两幅图像中，的权系数取为：
将小波边缘图与强边缘图融合后，得到的最终边缘图像用F(x,y)来表示，则各像素点的像素值D(x,y)表示为：
2.根据权利要求1所述的虹膜内边缘的快速精确定位方法，其特征在于：步骤(3)中利用小波变换提取边缘局部有效信息，具体是：
对于尺度参数j,平移参数k,将yj,k(t)中的连续变量j和k取做整数离散形式，将其表示为：
yj,k(t)＝2j/2y(2jt-k)
得到离散小波，将其表示为：
Wf(j,k)＝(f(t),yj,k(t))
小波变换通过参数j、k的调节，从而实现在频域与时域方面的局部化作用。
3.根据权利要求1所述的虹膜内边缘的快速精确定位方法，其特征在于：步骤(5)使用的是基于边缘信息的四向扫描法提取目标物，具体是：
采用从上到下、从左到右、从下到上、从右到左这四个方向，以不同颜色进行扫描填充，首先设置扫描长度参数和颜色参数，再是进行第一次的从上到下的区域扫描填充，其中上一次的填充色将成为下一次的扫描检测色，最后一个方向的填充色将是唯一的保留色，即分割出了所需要的提取目标物。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511006255.XA CN105488494B (zh) | 2015-12-29 | 2015-12-29 | 一种虹膜内边缘的快速精确的定位方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511006255.XA CN105488494B (zh) | 2015-12-29 | 2015-12-29 | 一种虹膜内边缘的快速精确的定位方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105488494A CN105488494A (zh) | 2016-04-13 |
| CN105488494B true CN105488494B (zh) | 2019-01-08 |
Family
ID=55675466
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511006255.XA Active CN105488494B (zh) | 2015-12-29 | 2015-12-29 | 一种虹膜内边缘的快速精确的定位方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105488494B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008091278A3 (en) * | 2006-09-25 | 2008-09-25 | Retica Systems Inc | Iris data extraction |
| CN101894256A (zh) * | 2010-07-02 | 2010-11-24 | 西安理工大学 | 基于奇对称2D Log-Gabor滤波器的虹膜识别方法 |
| CN103824061A (zh) * | 2014-03-03 | 2014-05-28 | 山东大学 | 基于光源反光区域检测和改进Hough变换的虹膜定位方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004288156A (ja) * | 2002-11-20 | 2004-10-14 | Stmicroelectronics Sa | 眼の虹彩イメージの定義の評価 |
-
2015
- 2015-12-29 CN CN201511006255.XA patent/CN105488494B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008091278A3 (en) * | 2006-09-25 | 2008-09-25 | Retica Systems Inc | Iris data extraction |
| CN101894256A (zh) * | 2010-07-02 | 2010-11-24 | 西安理工大学 | 基于奇对称2D Log-Gabor滤波器的虹膜识别方法 |
| CN103824061A (zh) * | 2014-03-03 | 2014-05-28 | 山东大学 | 基于光源反光区域检测和改进Hough变换的虹膜定位方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105488494A (zh) | 2016-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | Multifocus image fusion based on NSCT and focused area detection | |
| CN106296666B (zh) | 一种彩色图像去阴影方法和应用 | |
| CN104809734B (zh) | 一种基于引导滤波的红外图像与可见光图像融合的方法 | |
| CN105182350B (zh) | 一种应用特征跟踪的多波束声呐目标探测方法 | |
| CN106898015B (zh) | 一种基于自适应子块筛选的多线索视觉跟踪方法 | |
| CN102542275B (zh) | 一种证件照背景自动识别方法及其系统 | |
| CN109584204A (zh) | 一种图像噪声强度估计方法、存储介质、处理及识别装置 | |
| CN109242888A (zh) | 一种结合图像显著性和非下采样轮廓波变换的红外与可见光图像融合方法 | |
| CN105046206B (zh) | 基于视频中运动先验信息的行人检测方法及装置 | |
| CN105761214A (zh) | 基于轮廓波变换与引导滤波的遥感图像融合方法 | |
| CN111462027B (zh) | 基于多尺度梯度和抠图的多聚焦图像融合方法 | |
| CN109523506A (zh) | 基于视觉显著图像特征增强的全参考立体图像质量客观评价方法 | |
| CN106447640B (zh) | 基于字典学习、旋转引导滤波的多聚焦图像融合方法及装置 | |
| CN113592782B (zh) | 一种复合材料碳纤维芯棒x射线图像缺陷提取方法及系统 | |
| CN108830856B (zh) | 一种基于时间序列sd-oct视网膜图像的ga自动分割方法 | |
| CN110728302A (zh) | 一种基于HSV和Lab颜色空间的色纺织物组织识别的方法 | |
| CN109492647A (zh) | 一种电网机器人障碍物识别方法 | |
| CN109726681A (zh) | 一种基于机器学习识别与图像分割的盲道识别定位算法 | |
| CN106599891B (zh) | 一种基于尺度相位谱显著性的遥感图像兴趣区快速提取方法 | |
| CN109559273A (zh) | 一种面向车底图像的快速拼接方法 | |
| CN105513060A (zh) | 一种视觉感知启发的高分辨率遥感图像分割方法 | |
| CN104637060B (zh) | 一种基于邻域主成分分析-拉普拉斯的图像分割方法 | |
| CN106415606B (zh) | 一种基于边缘的识别、系统和方法 | |
| CN103679726A (zh) | 一种提高岩屑图像成像质量的方法 | |
| Sofou et al. | Generalized flooding and multicue PDE-based image segmentation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |