CN105631816A - 一种虹膜图像噪声分类检测方法 - Google Patents
一种虹膜图像噪声分类检测方法 Download PDFInfo
- Publication number
- CN105631816A CN105631816A CN201510968095.0A CN201510968095A CN105631816A CN 105631816 A CN105631816 A CN 105631816A CN 201510968095 A CN201510968095 A CN 201510968095A CN 105631816 A CN105631816 A CN 105631816A
- Authority
- CN
- China
- Prior art keywords
- iris
- eyelashes
- region
- noise
- eyelid
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 46
- 210000000720 eyelash Anatomy 0.000 claims abstract description 78
- 210000000744 eyelid Anatomy 0.000 claims abstract description 59
- 210000001747 pupil Anatomy 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 19
- 230000011218 segmentation Effects 0.000 claims description 15
- 238000009826 distribution Methods 0.000 claims description 11
- 238000002372 labelling Methods 0.000 claims description 10
- 238000006424 Flood reaction Methods 0.000 claims description 2
- 229910052705 radium Inorganic materials 0.000 claims description 2
- HCWPIIXVSYCSAN-UHFFFAOYSA-N radium atom Chemical compound [Ra] HCWPIIXVSYCSAN-UHFFFAOYSA-N 0.000 claims description 2
- 238000000605 extraction Methods 0.000 abstract description 11
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000009795 derivation Methods 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Ophthalmology & Optometry (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种虹膜图像噪声分类检测方法,该方法的步骤包括检测虹膜区域内的光源像点噪声S1、检测虹膜区域内的眼睑噪声S2、检测虹膜区域内的睫毛噪声S3、和基于步骤S1至S3中检测得到的光源像点、眼睑和睫毛的噪声信息,标记虹膜区域内的所有噪声点,并剔除该像素点S4。本发明所述技术方案能够高效的、准确的剔除噪声对虹膜定位和特征提取的干扰。
Description
技术领域
本发明涉及一种图像噪声检测方法,特别是涉及一种用于虹膜生物识别系统的虹膜图像噪声检测方法。
背景技术
虹膜识别具有准确性高、稳定性强、防伪性好和非接触检测等显著优势,是目前已知最先进的生物识别技术。在虹膜识别系统中,获取高质量的虹膜图像是提高识别效率的前提。然而,在实际应用中,虹膜图像中的虹膜区域会被光源像点、眼睑和睫毛等噪声干扰;如果不对这些噪声进行检测并将其干扰予以剔除,势必会影响虹膜定位的准确性,被提取的无效虹膜特征也会增加,从而使虹膜识别的准确率大大降低。因此,在虹膜图像预处理过程中需要检测虹膜噪声并将其干扰剔除,以增强虹膜定位的准确性和虹膜特征提取的有效性,从而保证虹膜识别的准确率。
虹膜图像中的噪声主要包括那些出现在虹膜区域内的光源像点、眼睑和睫毛部分,它们易于遮挡虹膜的边缘和纹理,干扰虹膜定位和特征提取,且难以消除。在目前的虹膜识别系统中,常用的虹膜图像噪声检测方法主要是通过阈值分割和模型匹配或粗略或精细地筛选噪声点并加以标记,剔除其对虹膜图像预处理和特征编码的干扰,从而减小虹膜图像噪声对虹膜识别准确率的影响。这类方法虽然能对虹膜图像中的各类噪声进行检测,但是存在以下明显的缺点:
1、阈值分割法主要用于光源像点和睫毛噪声的检测,这两类噪声的灰度分布没有统一的规律,不同区域的噪声点对选取的阈值较为敏感,并且虹膜图像的整体灰度分布受背景光照影响严重。单一的阈值选取仅能粗略地检测噪声;自适应的多阈值选取较为困难,且影响检测速度。因此,利用一般的阈值分割法检测噪声,效率较低。
2、模型匹配法主要用于眼睑和睫毛噪声的检测,这两类噪声的区域分布不具有规则的形状,且因人而异。想要建立合适的模型去匹配噪声点的分布比较困难,粗略的模型在检测噪声的同时会将部分有效的虹膜区域视为噪声,使可用的虹膜特征点减少,从而影响虹膜识别的效率;精细的模型设计复杂,噪声检测速度慢,同样会影响虹膜识别的效率。因此,单纯地利用模型匹配法检测虹膜图像的噪声,很难兼顾检测的速度与精度。
发明内容
本发明要解决的技术问题是提供一种虹膜图像噪声分类检测方法,以解决虹膜识别系统中虹膜图像干扰噪声影响系统识别率的问题。
为解决上述技术问题,本发明采用下述技术方案:
一种虹膜图像噪声分类检测方法,该方法的步骤包括
S1、检测虹膜区域内的光源像点噪声;
S2、检测虹膜区域内的眼睑噪声;
S3、检测虹膜区域内的睫毛噪声;
S4、基于步骤S1至S3中检测得到的光源像点、眼睑和睫毛的噪声信息,标记虹膜区域内的所有噪声点,并剔除该像素点。
优选的,所述步骤S1包括:
S11、确定虹膜区域中瞳孔的位置参数,所述位置参数包括瞳孔半径、中心点的行坐标和列坐标;
S12、对虹膜区域中瞳孔区域的像素灰度进行统计,得到分布曲线中小于高灰度峰值的波谷中对应的最大灰度值gssa:
并确定虹膜区域内光源像点的灰度分割阈值gsst:gsst=gssa+α,其中,α<0为安全系数;
S13、利用阈值分割结合连通域筛选的方法对虹膜区域内的光源像点进行检测,将Ψ内灰度值小于gsst,且聚集在n-连通域内的像素点作为光源像点SD:
SD={(x,y)|[x,y]=find(bwlabel(gs(Ψ)>gsst,n))}。
优选的,所述步骤S2包括:
S21、划定上、下眼睑的待搜索区域Φt和Φb:
其中,x,y分别表示虹膜图像中的行、列坐标变量;
S22、对下眼睑的待搜索区域Φb进行恢复对比扩展:Φb=imadjust(Φb);
S23、将Prewitt边缘提取与误差平方和最小化抛物线拟合方法相结合,分别对上眼睑和下眼睑的边缘进行检测,获得虹膜区域内的眼睑像素点集:
ED=EDt∪EDb,其中,EDt为虹膜区域内的上眼睑,EDb为虹膜区域内的下眼睑。
优选的,所述步骤S3包括:
S31、划定睫毛的待搜索区域Θ:
Θ={(x,y)|xp-rim≤x≤xp+rp,yp-rim≤y≤yp+rim}\Ω;
S32、将上眼睑边缘附件的睫毛灰度较低且成簇状的部分作为浓密睫毛的检测区域Θd:Θd={(x,y)|xp-rim≤x≤max(x′),y=y′,(x′,y′)∈Θ∩EDt};
将上眼睑边缘远离眼睑边缘的睫毛零落分散的部分作为稀疏睫毛的检测区域Θs:Θs=Θ\Θd;
S33、利用灰度直方图统计自适应地确定浓密睫毛的分割阈值gshd:
其中,λ>0为安全系数;并确定虹膜区域内的浓密睫毛点集EHd为:
EHd={(x,y)|gs(x,y)<gshd,(x,y)∈Θd};
S34、选取邻近瞳孔边缘的两个虹膜区域小像素块Λ1和Λ2作为参考,根据Λ1和Λ2的灰度分布确定稀疏睫毛的分割阈值gshs:
其中,ω、η是位置参数,σ是经验参数;并确定虹膜区域内系数睫毛点集EHs为:EHs={(x,y)|gs(x,y)<gshs,(x,y)∈Θs};
S35、根据虹膜区域内的浓密睫毛点集和稀疏睫毛点集,获得虹膜区域内的睫毛像素点集EH:EH=EHd∪EHs。
优选的,所述步骤S4包括:将噪声淹没图像定义为M:
其中,size(M)=size(I),I为原始虹膜图像,SD光源像点,ED眼睑像素点集,EH睫毛像素点集。
本发明的有益效果如下:
本发明所述技术方案从影响虹膜图像质量最主要的几个因素着手,快速而有效地剔除噪声干扰对虹膜定位与特征提取的影响,同时可兼顾虹膜图像预处理算法与特征提取匹配算法的适应性,从而提高虹膜识别的效率。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明;
图1示出本发明所述一种虹膜图像噪声分类检测方法的示意图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
如图1所示,本发明公开了一种虹膜图像噪声分类检测方法的具体步骤是:
第一步,检测虹膜区域内的光源像点噪声
由于虹膜识别系统的成像模块采用红外光源照明,在虹膜图像中的人眼部分会产生光源像点,其中位于虹膜区域的光源像点会遮挡虹膜纹理而影响识别的准确率,因此应当被视为噪声检测出来。一般地,光源像点的灰度值远远大于其他像素点,且分布集中,可以利用阈值分割结合连通域筛选的方法进行检测。由于瞳孔内的像素灰度分布对比明显,分为瞳孔低灰度像素区域和光源像点高灰度像素区域两部分,因此可以利用瞳孔内的光源像点灰度分布自适应地确定虹膜区域内光源像点的灰度分割阈值。
在检测噪声前,已对虹膜图像进行了瞳孔定位,已知瞳孔的位置参数为(xp,yp,rp),其中,rp、xp和yp分别表示瞳孔的半径和中心点的行、列坐标。设瞳孔区域为Ω,其像素灰度矩阵为gs(Ω),统计Ω的像素灰度分布,得到分布曲线中小于高灰度峰值的波谷中对应的最大灰度值gssa:
因此,得到虹膜区域内光源像点的灰度分割阈值gsst:
gsst=gssa+α(2)
其中,α<0为安全系数。
利用阈值分割结合连通域筛选的方法对虹膜区域内的光源像点进行检测。设待搜索虹膜区域为以(xp,yp)为中心点,rim为半径的圆形区域Ψ,其中rim为虹膜区域半径的最大限定值。将Ψ内灰度值小于gsst,且聚集在n-连通域(n一般取值为4或8,表示在n个方向上像素的连通性)内的像素点视作光源像点SD:
SD={(x,y)|[x,y]=find(bwlabel(gs(Ψ)>gsst,n))}(3)
第二步,检测虹膜区域内的眼睑噪声
虹膜区域易受到上、下眼睑的遮挡,尤其受上眼睑遮挡较为严重。眼睑遮挡虹膜纹理,使无效特征点增多,大大影响虹膜定位和识别的准确性,因此必须将虹膜区域内的眼睑噪声检测出来,剔除其造成的干扰。
首先,划定上、下眼睑的待搜索区域Φt和Φb:
其中,x,y分别表示虹膜图像中的行、列坐标变量。由于下眼睑边缘像素点的灰度变化不明显,因此对Φb进行灰度对比度扩展:
Φb=imadjust(Φb)(5)
由于眼睑的边缘近似抛物线形状,因此利用Prewitt边缘提取结合误差平方和最小化抛物线拟合的方法分别检测上、下眼睑的边缘。设眼睑边缘的待搜索区域为Φ,对其进行水平方向的Prewitt边缘提取,得到眼睑边缘的样本图像Γ:
Γ={(x,y)|[x,y]=find(edge(Φ,′prewitt′,′horizontal′))}(6)
设眼睑边缘形状满足:
y=ax2+bx+c(7)
其中,P=[a,b,c]T为抛物线参数向量。对Γ上的点(xn,yn)进行抛物线拟合,误差平方和δ(a,b,c)为:
δ(a,b,c)=∑(yn-(axn 2+bxn+c))2(8)
根据误差平方和最小化准则,令
求得P:
其中,Xn、Yn分别为:
因此,得到眼睑边缘L:
L={(x,y)∈Γ|y=ax2+bx+c}(12)
对应的眼睑区域ED为:
ED=Φ(L)(13)
根据以上推导,令Φ=Φt,检测到虹膜区域内的上眼睑为EDt;令Φ=Φb,检测到虹膜区域内的下眼睑为EDb。因此,得到虹膜区域内的眼睑像素点集ED:
ED=EDt∪EDb(14)
第三步,检测虹膜区域内的睫毛噪声
虹膜区域易受到睫毛的遮挡,尤其是上眼睑的睫毛。睫毛遮挡虹膜纹理,干扰虹膜定位和特征提取,使识别准确率降低,因此必须将虹膜区域内的睫毛噪声检测出来,剔除其造成的干扰。睫毛多而杂乱,形状各异,很难采用固定的阈值和模板分割检测,因此根据睫毛的分布特点,利用自适应的阈值分别对浓密睫毛和稀疏睫毛进行检测。
首先,划定睫毛的待搜索区域Θ:
Θ={(x,y)|xp-rim≤x≤xp+rp,yp-rim≤y≤yp+rim}\Ω(15)
由于下眼睑的睫毛对虹膜区域的干扰很少,因此只对上眼睑的睫毛进行检测。其中,位于上眼睑边缘附近的睫毛灰度较低且成簇状,被视作浓密睫毛;离上眼睑边缘较远的睫毛零落分散,被视作稀疏睫毛。浓密睫毛的检测区域为Θd:
Θd={(x,y)|xp-rim≤x≤max(x′),y=y′,(x′,y′)∈Θ∩EDt}(16)
稀疏睫毛的检测区域Θs为:
Θs=Θ\Θd(17)
浓密睫毛像素点的灰度值较低,且与邻近的上眼睑像素点灰度对比明显,因此利用灰度直方图统计自适应地确定浓密睫毛的分割阈值gshd:
其中,λ>0为安全系数。因此,检测到虹膜区域内的浓密睫毛点集EHd为:
EHd={(x,y)|gs(x,y)<gshd,(x,y)∈Θd}(19)
其中,gs(x,y)表示坐标(x,y)处像素点的灰度值。
稀疏睫毛像素点的灰度值虽然也较低,但与其邻近虹膜区域的像素点灰度对比并不明显,因此选取邻近瞳孔边缘的两个虹膜区域小像素块Λ1和Λ2作为参考,根据Λ1和Λ2的灰度分布确定稀疏睫毛的分割阈值gshs:
其中,ω、η是位置参数,σ是经验参数。因此,检测到虹膜区域内的稀疏睫毛点集EHs为:
EHs={(x,y)|gs(x,y)<gshs,(x,y)∈Θs}(22)
根据以上推导,得到虹膜区域内的睫毛像素点集EH:
EH=EHd∪EHs(23)
第四步,标记虹膜区域内的所有噪声点。
由以上三步分别检测到虹膜图像中的光源像点、眼睑和睫毛,将这些干扰噪声点在虹膜图像中的坐标位置用掩膜图像标记,以剔除其对虹膜定位和特征提取的影响。将噪声掩膜图像定义为M:
其中,size(M)=size(I),I为原始虹膜图像。
经过以上四步,将虹膜图像中的噪声点检测出来并加以标记,能有效抑制其对虹膜定位和特征提取的干扰,从而大大提高虹膜识别的准确率。
下面通过一组实施例对本发明做进一步说明:
一种虹膜图像噪声分类检测方法的具体步骤是:
第一步,检测虹膜区域内的光源像点噪声。
在检测噪声前,已对虹膜图像进行了瞳孔定位,已知瞳孔的位置参数为(xp,yp,rp),其中rp、xp和yp分别表示瞳孔的半径和中心点的行、列坐标。设瞳孔区域为Ω,其像素灰度矩阵为gs(Ω),统计Ω的像素灰度分布,得到分布曲线中小于高灰度峰值的波谷中对应的最大灰度值gssa:
因此,得到虹膜区域内光源像点的灰度分割阈值gsst:
gsst=gssa+α(2)
其中,α<0为安全系数,这里取α=(255-gssa)/4。
利用阈值分割结合连通域筛选的方法对虹膜区域内的光源像点进行检测。设待搜索虹膜区域为以(xp,yp)为中心点,rim为半径的圆形区域Ψ,其中rim为虹膜区域半径的最大限定值。将Ψ内灰度值小于gsst,且聚集在n-连通域内的像素点视作光源像点SD:
SD={(x,y)|[x,y]=find(bwlabel(gs(Ψ)>gsst,n))}(3)
式中,取n=8。
第二步,检测虹膜区域内的眼睑噪声。
首先,划定上、下眼睑的待搜索区域Φt和Φb:
其中,x,y分别表示虹膜图像中的行、列坐标变量。由于下眼睑边缘像素点的灰度变化不明显,因此对Φb进行灰度对比度扩展:
Φb=imadjust(Φb)(5)
由于眼睑的边缘近似抛物线形状,因此利用Prewitt边缘提取结合误差平方和最小化抛物线拟合的方法分别检测上、下眼睑的边缘。设眼睑边缘的待搜索区域为Φ,对其进行水平方向的Prewitt边缘提取,得到眼睑边缘的样本图像Γ:
Γ={(x,y)|[x,y]=find(edge(Φ,′prewitt′,′horizontal′))}(6)
设眼睑边缘形状满足:
y=ax2+bx+c(7)
其中,P=[a,b,c]T为抛物线参数向量。对Γ上的点(xn,yn)进行抛物线拟合,误差平方和δ(a,b,c)为:
δ(a,b,c)=∑(yn-(axn 2+bxn+c))2(8)
根据误差平方和最小化准则,令
求得P:
其中,Xn、Yn分别为:
因此,得到眼睑边缘L:
L={(x,y)∈Γ|y=ax2+bx+c}(12)
对应的眼睑区域ED为:
ED=Φ(L)(13)
根据以上推导,令Φ=Φt,检测到虹膜区域内的上眼睑为EDt;令Φ=Φb,检测到虹膜区域内的下眼睑为EDb。因此,得到虹膜区域内的眼睑像素点集ED:
ED=EDt∪EDb(14)
第三步,检测虹膜区域内的睫毛噪声。
首先,划定睫毛的待搜索区域Θ:
Θ={(x,y)|xp-rim≤x≤xp+rp,yp-rim≤y≤yp+rim}\Ω(15)
由于下眼睑的睫毛对虹膜区域的干扰很少,因此只对上眼睑的睫毛进行检测。其中,位于上眼睑边缘附近的睫毛灰度较低且成簇状,被视作浓密睫毛;离上眼睑边缘较远的睫毛零落分散,被视作稀疏睫毛。浓密睫毛的检测区域为Θd:
Θd={(x,y)|xp-rim≤x≤max(x′),y=y′,(x′,y′)∈Θ∩EDt}(16)
稀疏睫毛的检测区域Θs为:
Θs=Θ\Θd(17)
浓密睫毛像素点的灰度值较低,且与邻近的上眼睑像素点灰度对比明显,因此利用灰度直方图统计自适应地确定浓密睫毛的分割阈值gshd:
其中,λ>0为安全系数,本实例中,λ=gshd’/5,其中,gshd’是gshd的辅助计算量,是计算gshd的中间值。因此,检测到虹膜区域内的浓密睫毛点集EHd为:
EHd={(x,y)|gs(x,y)<gshd,(x,y)∈Θd}(19)
其中,gs(x,y)表示坐标(x,y)处像素点的灰度值。
稀疏睫毛像素点的灰度值虽然也较低,但与其邻近虹膜区域的像素点灰度对比并不明显,因此选取邻近瞳孔边缘的两个虹膜区域小像素块Λ1和Λ2作为参考,根据Λ1和Λ2的灰度分布确定稀疏睫毛的分割阈值gshs:
其中,ω、η是位置参数,本实例中,ω=rp/3,η=rp/5,σ是经验参数,本实例中,σ=0.75。因此,检测到虹膜区域内的稀疏睫毛点集EHs为:
EHs={(x,y)|gs(x,y)<gshs,(x,y)∈Θs}(22)
根据以上推导,得到虹膜区域内的睫毛像素点集EH:
EH=EHd∪EHs(23)
第四步,标记虹膜区域内的所有噪声点。
由以上三步分别检测到虹膜图像中的光源像点、眼睑和睫毛,将这些干扰噪声点在虹膜图像中的坐标位置用掩膜图像标记,以剔除其对虹膜定位和特征提取的影响。将噪声掩膜图像定义为M:
其中,size(M)=size(I),I为原始虹膜图像。
经过以上四步,将虹膜图像中的噪声点检测出来并加以标记,能有效抑制其对虹膜定位和特征提取的干扰,从而大大提高虹膜识别的准确率。
综上所述,本发明所述技术方案着重检测光源像点、眼睑和睫毛等影响虹膜图像质量最主要的几种噪声,结合虹膜定位算法和特征提取算法的适应性,快速而有效地剔除噪声干扰。在实际应用中,本发明所述技术方案大大缩短了虹膜图像预处理的时间,同时保证了虹膜定位的准确性和特征提取的稳定性,从而有效地提高了虹膜识别的效率。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (5)
1.一种虹膜图像噪声分类检测方法,其特征在于,该方法的步骤包括
S1、检测虹膜区域内的光源像点噪声;
S2、检测虹膜区域内的眼睑噪声;
S3、检测虹膜区域内的睫毛噪声;
S4、基于步骤S1至S3中检测得到的光源像点、眼睑和睫毛的噪声信息,标记虹膜区域内的所有噪声点,并剔除该像素点。
2.根据权利要求1所述的检测方法,其特征在于,所述步骤S1包括:
S11、确定虹膜区域中瞳孔的位置参数,所述位置参数包括瞳孔半径、中心点的行坐标和列坐标;
S12、对虹膜区域中瞳孔区域的像素灰度进行统计,得到分布曲线中小于高灰度峰值的波谷中对应的最大灰度值gssa:
并确定虹膜区域内光源像点的灰度分割阈值gsst:gsst=gssa+α,其中,α<0为安全系数;
S13、利用阈值分割结合连通域筛选的方法对虹膜区域内的光源像点进行检测,将Ψ内灰度值小于gsst,且聚集在n-连通域内的像素点作为光源像点SD:SD={(x,y)|[x,y]=find(bwlabel(gs(Ψ)>gsst,n))}。
3.根据权利要求1所述的检测方法,其特征在于,所述步骤S2包括:
S21、划定上、下眼睑的待搜索区域Φt和Φb:
其中,x,y分别表示虹膜图像中的行、列坐标变量;
S22、对下眼睑的待搜索区域Φb进行恢复对比扩展:Φb=imadjust(Φb);
S23、将Prewitt边缘提取与误差平方和最小化抛物线拟合方法相结合,分别对上眼睑和下眼睑的边缘进行检测,获得虹膜区域内的眼睑像素点集:
ED=EDt∪EDb,其中,EDt为虹膜区域内的上眼睑,EDb为虹膜区域内的下眼睑。
4.根据权利要求1所述的检测方法,其特征在于,所述步骤S3包括:
S31、划定睫毛的待搜索区域Θ:
Θ={(x,y)|xp-rim≤x≤xp+rp,yp-rim≤y≤yp+rim}\Ω;
S32、将上眼睑边缘附件的睫毛灰度较低且成簇状的部分作为浓密睫毛的检测区域Θd:Θd={(x,y)|xp-rim≤x≤max(x′),y=y′,(x′,y′)∈Θ∩EDt};
将上眼睑边缘远离眼睑边缘的睫毛零落分散的部分作为稀疏睫毛的检测区域Θs:Θs=Θ\Θd;
S33、利用灰度直方图统计自适应地确定浓密睫毛的分割阈值gshd:
其中,λ>0为安全系数;并确定虹膜区域内的浓密睫毛点集EHd为:
EHd={(x,y)|gs(x,y)<gshd,(x,y)∈Θd};
S34、选取邻近瞳孔边缘的两个虹膜区域小像素块Λ1和Λ2作为参考,根据Λ1和Λ2的灰度分布确定稀疏睫毛的分割阈值gshs:
其中,ω、η是位置参数,σ是经验参数;并确定虹膜区域内系数睫毛点集EHs为:EHs={(x,y)|gs(x,y)<gshs,(x,y)∈Θs};
S35、根据虹膜区域内的浓密睫毛点集和稀疏睫毛点集,获得虹膜区域内的睫毛像素点集EH:EH=EHd∪EHs。
5.根据权利要求1所述的检测方法,其特征在于,所述步骤S4包括:将噪声淹没图像定义为M:
其中,size(M)=size(I),I为原始虹膜图像,SD光源像点,ED眼睑像素点集,EH睫毛像素点集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510968095.0A CN105631816B (zh) | 2015-12-22 | 2015-12-22 | 一种虹膜图像噪声分类检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510968095.0A CN105631816B (zh) | 2015-12-22 | 2015-12-22 | 一种虹膜图像噪声分类检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105631816A true CN105631816A (zh) | 2016-06-01 |
CN105631816B CN105631816B (zh) | 2018-04-03 |
Family
ID=56046704
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510968095.0A Active CN105631816B (zh) | 2015-12-22 | 2015-12-22 | 一种虹膜图像噪声分类检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105631816B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107451546A (zh) * | 2017-07-14 | 2017-12-08 | 广东欧珀移动通信有限公司 | 虹膜识别方法及相关产品 |
CN109446935A (zh) * | 2018-10-12 | 2019-03-08 | 北京无线电计量测试研究所 | 一种用于远距离行进中虹膜识别的虹膜定位方法 |
CN110210357A (zh) * | 2019-05-24 | 2019-09-06 | 浙江大学 | 一种基于静态照片面部识别的上睑下垂图像测量方法 |
CN111950376A (zh) * | 2020-07-15 | 2020-11-17 | 杭州电子科技大学 | 强方向加权边缘检测及rnl拟合优度的眼睑定位方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101447025A (zh) * | 2008-12-30 | 2009-06-03 | 东南大学 | 一种大型动物虹膜识别方法 |
US20090252382A1 (en) * | 2007-12-06 | 2009-10-08 | University Of Notre Dame Du Lac | Segmentation of iris images using active contour processing |
CN104268527A (zh) * | 2014-09-26 | 2015-01-07 | 北京无线电计量测试研究所 | 一种基于径向梯度检测的虹膜定位方法 |
-
2015
- 2015-12-22 CN CN201510968095.0A patent/CN105631816B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090252382A1 (en) * | 2007-12-06 | 2009-10-08 | University Of Notre Dame Du Lac | Segmentation of iris images using active contour processing |
CN101447025A (zh) * | 2008-12-30 | 2009-06-03 | 东南大学 | 一种大型动物虹膜识别方法 |
CN104268527A (zh) * | 2014-09-26 | 2015-01-07 | 北京无线电计量测试研究所 | 一种基于径向梯度检测的虹膜定位方法 |
Non-Patent Citations (2)
Title |
---|
周俊: "一种用于虹膜识别的眼睑和睫毛检测算法", 《武汉理工大学学报 信息与管理工程》 * |
田启川等: "虹膜识别中噪声的检测和处理方法", 《计算机工程》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107451546A (zh) * | 2017-07-14 | 2017-12-08 | 广东欧珀移动通信有限公司 | 虹膜识别方法及相关产品 |
US10607076B2 (en) | 2017-07-14 | 2020-03-31 | Guangdong Oppo Mobile Telecommunications Corp., Ltd. | Method for iris recognition and related products |
CN109446935A (zh) * | 2018-10-12 | 2019-03-08 | 北京无线电计量测试研究所 | 一种用于远距离行进中虹膜识别的虹膜定位方法 |
CN109446935B (zh) * | 2018-10-12 | 2021-06-29 | 北京无线电计量测试研究所 | 一种用于远距离行进中虹膜识别的虹膜定位方法 |
CN110210357A (zh) * | 2019-05-24 | 2019-09-06 | 浙江大学 | 一种基于静态照片面部识别的上睑下垂图像测量方法 |
CN110210357B (zh) * | 2019-05-24 | 2021-03-23 | 浙江大学 | 一种基于静态照片面部识别的上睑下垂图像测量方法 |
CN111950376A (zh) * | 2020-07-15 | 2020-11-17 | 杭州电子科技大学 | 强方向加权边缘检测及rnl拟合优度的眼睑定位方法 |
CN111950376B (zh) * | 2020-07-15 | 2023-09-29 | 杭州电子科技大学 | 强方向加权边缘检测及rnl拟合优度的眼睑定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105631816B (zh) | 2018-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102043950B (zh) | 基于canny算子和边缘点统计的车辆轮廓识别方法 | |
CN104268872B (zh) | 基于一致性的边缘检测方法 | |
CN105046250A (zh) | 人脸识别的眼镜消除方法 | |
CN110443807A (zh) | 一种基于亮度均衡的光照不均匀图像阈值分割方法 | |
CN105631816A (zh) | 一种虹膜图像噪声分类检测方法 | |
CN109410230A (zh) | 一种可抗噪声的改进Canny图像边缘检测方法 | |
CN105117705B (zh) | 一种虹膜图像质量级联式评价方法 | |
CN106570447B (zh) | 基于灰度直方图匹配的人脸照片太阳镜自动去除方法 | |
CN103440476A (zh) | 一种人脸视频中瞳孔的定位方法 | |
CN102521836A (zh) | 一种基于特定类灰度图像的边缘检测方法 | |
CN103077386A (zh) | 一种级联的视频流虹膜图像质量检测方法 | |
CN106096491B (zh) | 一种眼底彩色照相图像中微动脉瘤自动化识别方法 | |
CN105095907B (zh) | 一种基于rbf神经网络的棉花异质纤维鉴别方法 | |
CN109886935A (zh) | 一种基于深度学习的道面异物检测方法 | |
CN106651888A (zh) | 基于多特征融合的彩色眼底图像视杯分割方法 | |
CN110245626B (zh) | 一种在虹膜图像中精确检测眼睫毛图像的方法 | |
CN110728185B (zh) | 一种判别驾驶人存在手持手机通话行为的检测方法 | |
CN117197140B (zh) | 基于机器视觉的不规则金属扣成型检测方法 | |
CN109377450A (zh) | 一种边缘保护的去噪方法 | |
CN108710832B (zh) | 一种无参考虹膜图像清晰度检测方法 | |
CN101615241B (zh) | 一种证件照片的筛选方法 | |
CN114119603A (zh) | 基于图像处理的快餐盒短射缺陷检测方法 | |
CN109840484A (zh) | 一种基于边缘滤波、椭圆评价和瞳孔验证的瞳孔检测方法 | |
CN104598914A (zh) | 一种肤色检测的方法及装置 | |
CN111476804A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |