CN108629744B - 一种图像增强方法 - Google Patents
一种图像增强方法 Download PDFInfo
- Publication number
- CN108629744B CN108629744B CN201810311065.6A CN201810311065A CN108629744B CN 108629744 B CN108629744 B CN 108629744B CN 201810311065 A CN201810311065 A CN 201810311065A CN 108629744 B CN108629744 B CN 108629744B
- Authority
- CN
- China
- Prior art keywords
- image
- iris
- sub
- subblocks
- iris image
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005070 sampling Methods 0.000 claims abstract description 18
- 238000001228 spectrum Methods 0.000 claims description 16
- 230000000903 blocking effect Effects 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000012163 sequencing technique Methods 0.000 claims description 4
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 5
- 238000000638 solvent extraction Methods 0.000 abstract 2
- 210000000554 iris Anatomy 0.000 description 107
- 230000004927 fusion Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 239000006002 Pepper Substances 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
-
- 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
-
- 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/30004—Biomedical image processing
- G06T2207/30041—Eye; Retina; Ophthalmic
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)
Abstract
本发明公开了一种虹膜图像增强方法,解决现有方法无法对低质虹膜图像进行增强的问题。所述方法,包含:获取归一化虹膜图像序列,选取基准图像、伴生图像;将伴生图像序列对照基准图像全局对齐,得到对齐伴生图像序列;将预处理图像序列进行上采样和分块,并将上采样后分块位置相同的子块按照图像质量高低存入堆栈,得到虹膜图像子块;将各堆栈中虹膜图像子块进行局部配准,得到虹膜图像配准子块;对各堆栈中的虹膜图像配准子块进行超分辨率多阶局部信息加权重建;将各堆栈中的虹膜图像重建子块进行归一化拼接、下采样,得到重建图像;将重建图像进行去噪和锐化处理,得到图像增强后的归一化虹膜图像。本发明有助于提高系统的准确率和鲁棒性。
Description
技术领域
本发明涉及虹膜图像识别领域,尤其涉及一种虹膜图像增强方法。
背景技术
目前的虹膜图像增强方法主要以对比度扩展、形态学修正、选择高通滤波和矩阵填充等手段为主,从提高图像的清晰度,去噪和复原虹膜细节信息等方面对单帧虹膜图像的全局信息进行增强,从而提升图像的质量,存在的缺点是,单帧低质虹膜图像的有效信息较少,对其进行去噪、去模糊等操作容易形成斑块效应,而图像复原操作填充的信息大部分为冗余信息,因此增强后的虹膜图像的信息熵并没有显著增加,无法实现有效的特征提取与识别,同时低质虹膜图像的全局变形和局部变形是异质的,对图像整体进行校准无法使各局部的变形得到最优的抑制,同时会使局部信息均质化,因此对感兴趣的虹膜区域的调谐具有不确定性,难以增强虹膜局部细节信息。
发明内容
本发明提供一种虹膜图像增强方法,解决现有方法对远距离动态复杂环境成像条件下的低质虹膜图像进行增强的问题。
一种虹膜图像增强方法,包含以下步骤:获取归一化虹膜图像序列,选取图像质量最高的图像作为基准图像,其余按照图像质量高低排序作为伴生图像;将所述伴生图像序列根据多方向汉明距离判别准则对照所述基准图像进行全局对齐,得到对齐伴生图像序列;将所述基准图像、对齐伴生图像组成的预处理图像序列中的每帧图像进行上采样和分块,并依次将上采样后图像序列中每帧图像的分块位置相同的子块按照图像质量高低的顺序存入堆栈中,得到虹膜图像子块;采用归一化互信息熵判别准则对各堆栈中的所述虹膜图像子块的基准子块和伴生子块进行局部配准,得到虹膜图像配准子块;对各堆栈中的所述虹膜图像配准子块进行超分辨率多阶局部信息加权重建,得到虹膜图像重建子块;将各堆栈中的所述虹膜图像重建子块按照原始相对位置进行归一化拼接,并下采样,得到重建图像;将所述重建图像进行去噪和锐化处理,得到图像增强后的归一化虹膜图像。
进一步地,所述图像质量高低由图像质量分数的高低决定,所述图像质量分数为:
QI=HFE/MFE
其中,QI为所述图像质量分数,HFE为在所述归一化虹膜图像序列的频谱中去除噪声频谱后的高频能量,MFE为在所述归一化虹膜图像序列的频谱中去除直流频谱后的低频能量。
优选地,所述将伴生图像序列根据多方向汉明距离判别准则对照所述基准图像进行全局对齐,得到对齐伴生图像序列的步骤,进一步包含:将所述伴生图像序列中的每帧图像沿四周方向扩展,得到扩展伴生图像序列;在所述扩展伴生图像序列的中开一个窗口,所述窗口大小与所述伴生图像序列中的每帧图像大小相同;对所述扩展伴生图像序列的每帧图像,逐行、逐列滑动所述窗口,计算所述窗口内图像与所述基准图像的汉明距离,将所述汉明距离最小的所述窗口内图像记为对齐图像,根据每帧图像的所述对齐图像得到所述对齐伴生图像序列。
优选地,所述上采样采用双立方插值法。
优选地,所述采用归一化互信息熵判别准则对各堆栈中的所述虹膜图像子块的基准子块和伴生子块进行局部配准,得到虹膜图像配准子块的步骤,进一步包含:计算各堆栈中所述虹膜图像子块中基准子块的信息熵;计算各堆栈中所述虹膜图像子块中每个伴生子块的信息熵;计算各堆栈中所述虹膜图像子块中基准子块与伴生子块联合熵;计算各堆栈中所述基准子块与所述伴生子块的归一化互信息熵为:
NMI(BB,CBm)=(H(BB)+H(CBm))/H(BB,CBm)
其中,BB为所述基准子块,CBm为第m个伴生子块,NMI(BB,CBm)为所述归一化互信息熵,H(BB)为所述基准子块的信息熵,H(CBm)为所述第m个伴生子块的信息熵,H(BB,CBm)为所述联合熵;将各堆栈中所述虹膜图像子块中的所述伴生子块与所述基准子块重合,并在沿水平方向、垂直方向滑动的范围内,找到所述归一化互信息熵最大的重叠区域为所述伴生子块的配准子块,并求得所述基准子块各像素点的累计配准次数矩阵。
进一步地,所述虹膜图像重建子块为:
RB=BB+diff(BB)+diff2(BB)+∑m=1→M-1wm*(ABm+diff(ABm)+diff2(ABm))
其中,RB为所述虹膜图像重建子块,BB为所述基准子块,ABm为质量排序为m的配准子块,wm为所述质量排序为m的配准子块的权值,diff为图像一阶导数函数,diff2为图像二阶导数函数。
优选地,所述将各堆栈中的所述虹膜图像重建子块按照原始相对位置进行归一化拼接,进一步包含:对各堆栈中的所述虹膜图像重建子块,利用归一化因子平滑拼接边缘,得到拼接归一化重建子块为:
RBe’=RBe/NRe
其中,RBe’为所述拼接归一化重建子块,RBe为所述虹膜图像重建子块,NRe为所述归一化因子。
优选地,所述去噪处理采用中值联合高斯低通滤波,所述锐化处理采用拉普拉斯算子。
优选地,所述下采样采用双线性插值法。
优选地,所述归一化虹膜图像序列的图像帧数为10,单帧图像分别率为288*64。
本发明有益效果包括:针对远距离动态虹膜识别系统成像易受用户形态影响,获取的虹膜图像往往存在变形、模糊和噪声等污染的特点,通过对序列低质虹膜图像进行分块超分辨率局部信息加权融合来增强虹膜的细节信息,解决了常规图像增强方法难以对远距离动态低质虹膜图像进行有效增强的问题;本发明方法能在抵抗虹膜图像变形、模糊和噪声污染的同时有效复原虹膜局部细节信息,从而有助于提高系统的准确率和鲁棒性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为一种虹膜图像增强方法流程实施例;
图2为一种包含全局对其的虹膜图像增强方法流程实施例。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
虹膜识别以其准确性、稳定性、安全性和非接触性等显著优势已成为生物识别领域的重点研究方向和发展趋势。当前的虹膜识别产品以近距离静止识别为主,需要用户进行较好的配合以实现精准的识别,这在一定程度上限制了虹膜识别技术的普及和推广。为了提高用户使用时的自由度和认可度,虹膜识别技术和产品必须向远距离动态虹膜识别应用场景扩展。而远距离动态虹膜识别的难点在于,由于用户在识别过程中形态各异,行进速度时快时慢,因此远距离动态虹膜识别系统的成像模块很难捕获到清晰的高质量的虹膜图像用于识别,获取的原始虹膜图像往往存在变形、模糊和噪声等污染。所以,由远距离动态虹膜识别系统采集的虹膜图像必须经过图像增强满足一定质量要求才能实现有效的识别,图像增强是远距离动态虹膜识别的关键。
以下结合附图,详细说明本发明各实施例提供的技术方案。
图1为一种虹膜图像增强方法流程实施例,本发明实施例提供一种虹膜图像增强方法,包含以下步骤:
步骤101,获取归一化虹膜图像序列,选取图像质量最高的图像作为基准图像,其余按照图像质量高低排序作为伴生图像。
进一步地,所述图像质量高低由图像质量分数的高低决定,所述图像质量分数为:
QI=HFE/MFE (1)
其中,QI为所述图像质量分数,HFE为在所述归一化虹膜图像序列的频谱中去除噪声频谱后的高频能量,MFE为在所述归一化虹膜图像序列的频谱中去除直流频谱后的低频能量。
在步骤101中,所述归一化虹膜图像序列为从远距离动态虹膜识别系统拍摄的视频流中选取连续多帧含有虹膜的图像,并对每帧图像进行虹膜定位和极坐标归一化处理,得到归一化虹膜图像序列。
在步骤101中,假设所述归一化虹膜图像序列的傅立叶变换频谱为F(u,v),其中心点为原点,假设RL={(u,v)|-x≤u≤x,-y≤v≤y,(u,v)≠(0,0)}为去除直流频谱的低频分量区,其系数的幅度平方和为低频能量LFE,RH={(u,v)|-x-x’≤u<-x∪x<u≤x+x’,-y-y’≤v≤y+y’}∪{(u,v)|-x≤u≤x,-y-y’≤v<-y∪y<v≤y+y’}为去除噪声频谱的高频分量区,其系数的幅度平方和为高频能量HFE,其中x、y为低频分量区的门限,x’、y’分别为高频分量区的门限,典型地,x=48,y=16,x’=88,y’=12,这里不做特别限定。由于傅立叶频谱的高频分量代表图像的细节信息,因此所述图像质量分数越高表示图像越清晰,质量越好。
需要说明的是,所述归一化虹膜图像序列的图像帧数可以为10,也可以为其他数值,所述图像的单帧分别率可以为288*64,也可以为其他数值,这里不做特别限定。
步骤102,将所述伴生图像序列根据多方向汉明距离判别准则对照所述基准图像进行全局对齐,得到对齐伴生图像序列。
在步骤102中,为了将伴生图像序列与基准图像进行融合以增强虹膜区域的细节,首先以基准图像作为参照,将伴生图像逐帧进行全局对齐,以减少由相邻帧图像的全局变形带来的干扰,利用多方向滑动汉明距离窗检测的方法将伴生图像与基准图像进行对齐,将汉明距离最小的窗内图像作为对齐后的图像。
步骤103,将所述基准图像、对齐伴生图像组成的预处理图像序列中的每帧图像进行上采样和分块,并依次将上采样后图像序列中每帧图像的分块位置相同的子块按照图像质量高低的顺序存入堆栈中,得到虹膜图像子块。
优选地,所述上采样采用双立方插值法。
在步骤103中,由于低质序列虹膜图像存在局部变形和失真,因此分块信息融合更有利于细节信息增强,为了将多帧图像进行分块超分辨率信息融合,首先对每帧图像进行上采样将图像分辨率扩展,典型地所述上采用为1:2上采样,假设所述图像的分别率为X×Y,则上采样后的图像分辨率为2X×2Y,需要说明的是,所述上采样也可以采用1:3上采样,或其他采样率上采样,这里不做特别限定。
在步骤103中,假设分块的大小为x×y的子块,因此单帧图像包含(2X/x)×(2Y/y)个子块。
假设所述预处理图像的分辨率为288*64,为了将多帧图像进行分块超分辨率信息融合,对每帧图像进行1:2上采样将图像分辨率扩展至576×128,然后将其划分成大小为48×32的子块,因此单帧图像包含12×4=48个子块,依次将上采样后的序列中每帧图像的相同位置的子块按照质量高低的顺序存入堆栈中,得到48个堆栈,每个堆栈中含有10个子块。
其中,图像上采样采用双立方插值法:
I’=imresize(I,[128 576],‘bicubic’) (2)
其中,imresize为图像插值函数,I表示分辨率为288×64的原图像,I’表示分辨率为576×128的1:2上采样图像。
需要说明的是,所述上采样的采样率、分块的子块分辨率可以是本发明实施例中的数值,也可以是其他数值,这里不做特别限定。
步骤104,采用归一化互信息熵判别准则对各堆栈中的所述虹膜图像子块的基准子块和伴生子块进行局部配准,得到虹膜图像配准子块。
步骤105,对各堆栈中的所述虹膜图像配准子块进行超分辨率多阶局部信息加权重建,得到虹膜图像重建子块。
进一步地,所述虹膜图像重建子块为:
RB=BB+diff(BB)+diff2(BB)+∑m=1→M-1wm*(ABm+diff(ABm)+diff2(ABm)) (3)
其中,RB为所述虹膜图像重建子块,BB为所述基准子块,ABm为所述质量排序为m的配准子块,wm为所述质量排序为m的配准子块的权值,diff为图像一阶导数函数,diff2为图像二阶导数函数。
在步骤105中,由于图像的一阶导数可以提取较粗的边缘信息,二阶导数能提取更为精细的纹理信息,因此通过叠加所述基准子块、配准伴生子块及其一阶、二阶导数实现超分辨率局部信息融合,以增强图像的局部细节信息。考虑到各配准伴生子块图像质量的高低,通过加权叠加重建局部细节信息,计算各堆栈的所述虹膜图像配准子块的权值为:
wm=e(8M-7m-9)/(M-2)/∑m=1→M-1e(8M-7m-9)/(M-2) (4)
其中,wm为所述质量排序为m的配准子块的权值,M为所述归一化虹膜图像序列的图像帧数,m为所述各堆栈中配准子块的质量序数。
步骤106,将各堆栈中的所述虹膜图像重建块按照原始相对位置进行归一化拼接,并下采样,得到重建图像。
进一步地,对各堆栈中的所述虹膜图像重建子块,利用归一化因子平滑拼接边缘,得到拼接归一化重建子块为:
RBe’=RBe/NRe (5)
其中,RBe’为所述拼接归一化重建子块,RBe为所述虹膜图像重建子块,NRe为所述归一化因子。
在步骤106中,假设所述虹膜图像重建子块的分辨率为(2X/x)×(2Y/y),为了消除子块拼接时产生的边界灰度跳跃现象,利用归一化因子平滑拼接边缘,并进行下采样,假设所述下采样为2:1下采样,则所述重建图像的分辨率为X×Y。优选地,所述下采样的方法为双线性插值法。
假设在各堆栈中的所述虹膜图像中,有48个堆栈,每个堆栈有10个虹膜图像,每帧图像的分辨率为576×128,将48个重建子块按照原始相对位置进行拼接,并进行2:1下采样,得到一幅分辨率大小为288×64的超分辨率信息融合重建归一化虹膜图像。其中,图像下采样采用双线性插值法:
P’=imresize(P,[64 288],‘bilinear’) (6)
其中,imresize为图像插值函数,P表示分辨率为576×128的原图像,P’表示分辨率为288×64的2:1下采样图像。
步骤107,将所述重建图像进行去噪和锐化处理,得到图像增强后的归一化虹膜图像。
在步骤107中,将重建图像进行去噪和锐化等后处理的目的是进一步提升图像的对比度和清晰度,最终得到一幅增强后的归一化虹膜图像。
在步骤107中,去噪采用中值联合高斯低通滤波,以有效去除白噪声和椒盐噪声:
RI’=imfilter(RI,‘median’*‘Gaussian_low’) (7)
其中,imfilter表示二维中值联合高斯低通滤波函数,RI’表示去噪后的图像。
在步骤107中,锐化采用拉普拉斯算子,以有效增强局部细节像素的灰度对比度:
EI=del2(RI’) (8)
其中,RI’表示去噪后的图像,del2表示离散拉普拉斯变换函数,EI为增强后的归一化虹膜图像。
本发明实施例提供了一种虹膜图像增强方法,所述虹膜图像序列为连续10帧低质的归一化虹膜图像序列,通过分块超分辨率局部信息加权融合的方法,得到一幅增强后的归一化虹膜图像。可以看出,本发明所述技术方案,能在抵抗虹膜图像变形、模糊和噪声污染的同时有效复原虹膜局部细节信息,解决了常规图像增强方法难以对远距离动态低质虹膜图像进行有效增强的问题,从而有助于提高系统的准确率和鲁棒性。
图2为一种包含全局对其的虹膜图像增强方法流程实施例,具体包含以下步骤:
步骤101,获取归一化虹膜图像序列,选取图像质量最高的图像作为基准图像,其余按照图像质量高低排序作为伴生图像。
步骤201,将所述伴生图像序列中的每帧图像沿四周方向扩展,得到扩展伴生图像序列。
在步骤201中,假设所述伴生图像分别沿四周的方向扩展p列和q行,其分辨率大小变为(X+2p)×(Y+2q)。
步骤202,在所述扩展伴生图像序列的中开一个窗口,所述窗口大小与所述伴生图像序列中的每帧图像大小相同。
步骤203,对所述扩展伴生图像序列的每帧图像,逐行、逐列滑动所述窗口,计算所述窗口内图像与所述基准图像的汉明距离,将所述汉明距离最小的所述窗口内图像记为对齐图像,根据每帧图像的所述对齐图像得到所述对齐伴生图像序列。
在步骤203中,在(X+2p)×(Y+2q)大小的伴生图像中开一个X×Y大小的窗口,逐行、逐列滑动窗口并计算窗内图像与基准图像的汉明距离,将汉明距离最小的窗内图像记为对齐图像,逐帧对伴生图像序列进行全局对齐,得到对齐伴生图像序列。
假设所述伴生图像的分辨率大小为288×64,将伴生图像分别沿四周的方向扩展8列和8行,则扩展后的伴生图像分辨率大小变为304×80。在304×80大小的伴生图像中开一个288×64大小的窗口,逐行、逐列滑动窗口并计算窗内图像与基准图像的汉明距离,将汉明距离最小的窗内图像记为对齐图像。逐帧对伴生图像序列进行全局对齐,得到对齐伴生图像序列。
步骤103,将所述基准图像、对齐伴生图像组成的预处理图像序列中的每帧图像进行上采样和分块,并依次将上采样后图像序列中每帧图像的分块位置相同的子块按照图像质量高低的顺序存入堆栈中,得到虹膜图像子块。
步骤204,计算各堆栈中所述虹膜图像子块中基准子块的信息熵。
所述基准子块的信息熵为:
H(BB)=-∑i(silogsi) (9)
其中,BB为所述基准子块,H(BB)为所述基准子块的信息熵,si表示灰度级i在所述基准子块中出现的概率,i=0,1,…,255。
步骤205,计算各堆栈中所述虹膜图像子块中每个伴生子块的信息熵。
所述每个伴生子块的信息熵为:
H(CBm)=-∑i(tilogti) (10)
其中,CBm为第m个伴生子块,H(CBm)为所述第m个伴生子块的信息熵,ti表示灰度级i在所述伴生子块CBm中出现的概率,i=0,1,…,255。
步骤206,计算各堆栈中所述虹膜图像子块中基准子块与伴生子块联合熵。
所述联合熵为:
H(BB,CBm)=-∑j,k(z(j,k)logz(j,k)) (11)
其中,BB为所述基准子块,CBm为所述第m个伴生子块,H(BB,CBm)为所述联合熵,z(j,k)为灰度级对(j,k)在所述基准子块BB和所述第m个伴生子块CBm相应位置上出现的概率,j,k=0,1,…,255。
步骤207,计算各堆栈中所述基准子块与所述伴生子块的归一化互信息熵。
所述基准子块与所述伴生子块的归一化互信息熵为:
NMI(BB,CBm)=(H(BB)+H(CBm))/H(BB,CBm) (12)
其中,BB为所述基准子块,CBm为所述第m个伴生子块,NMI(BB,CBm)为所述归一化互信息熵,H(BB)为所述基准子块的信息熵,H(CBm)为所述第m个伴生子块的信息熵,H(BB,CBm)为所述联合熵。
在步骤207中,所述基准子块与所述伴生子块的归一化互信息熵越大,表示所述基准子块与所述伴生子块两幅图像越相似。
步骤208,将各堆栈中所述虹膜图像子块中的所述伴生子块与所述基准子块重合,并在沿水平方向、垂直方向滑动的范围内,找到所述归一化互信息熵最大的重叠区域为所述伴生子块的配准子块,并求得所述基准子块各像素点的累计配准次数矩阵。
在步骤208中,所述累计配准次数矩阵的个数为0,1,…,(2X/x)×(2Y/y),其中X、Y为所述归一化虹膜图像序列的图像分辨率,x、y为所述分块的子块分辨率。
步骤105,对各堆栈中的所述虹膜图像配准子块进行超分辨率多阶局部信息加权重建,得到虹膜图像重建子块。
步骤106,将各堆栈中的所述虹膜图像重建子块按照原始相对位置进行归一化拼接,并下采样,得到重建图像。
步骤107,将所述重建图像进行去噪和锐化处理,得到图像增强后的归一化虹膜图像。
需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上所述仅为本发明的实施例而已,并不用于限制本发明。对于本领域技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (8)
1.一种虹膜图像增强方法,其特征在于,包含以下步骤:
获取归一化虹膜图像序列,选取图像质量最高的图像作为基准图像,其余按照图像质量高低排序作为伴生图像;
将所述伴生图像序列根据多方向汉明距离判别准则对照所述基准图像进行全局对齐,得到对齐伴生图像序列;
将所述基准图像、对齐伴生图像组成的预处理图像序列中的每帧图像进行上采样和分块,并依次将上采样后图像序列中每帧图像的分块位置相同的子块按照图像质量高低的顺序存入堆栈中,得到虹膜图像子块;
计算各堆栈中所述虹膜图像子块中基准子块的信息熵;
计算各堆栈中所述虹膜图像子块中每个伴生子块的信息熵;
计算各堆栈中所述虹膜图像子块中基准子块与伴生子块联合熵;
计算各堆栈中所述基准子块与所述伴生子块的归一化互信息熵为:
NMI(BB,CBm)=(H(BB)+H(CBm))/H(BB,CBm)
其中,BB为所述基准子块,CBm为第m个伴生子块,NMI(BB,CBm)为所述归一化互信息熵,H(BB)为所述基准子块的信息熵,H(CBm)为所述第m个伴生子块的信息熵,H(BB,CBm)为所述联合熵;
将各堆栈中所述虹膜图像子块中的所述伴生子块与所述基准子块重合,并在沿水平方向、垂直方向滑动的范围内,找到所述归一化互信息熵最大的重叠区域为所述伴生子块的配准子块,并求得所述基准子块各像素点的累计配准次数矩阵;
对各堆栈中的所述虹膜图像配准子块进行重建,得到虹膜图像重建子块为:
RB=BB+diff(BB)+diff2(BB)+∑m=1→M-1wm*(ABm+diff(ABm)+diff2(ABm))
其中,RB为所述虹膜图像重建子块,BB为所述基准子块,ABm为质量排序为m的配准子块,em为所述质量排序为m的配准子块的权值,diff为图像一阶导数函数,diff2为图像二阶导数函数;
将各堆栈中的所述虹膜图像重建子块按照原始相对位置进行归一化拼接,并下采样,得到重建图像;
将所述重建图像进行去噪和锐化处理,得到图像增强后的归一化虹膜图像。
2.如权利要求1所述的虹膜图像增强方法,其特征在于,所述图像质量高低由图像质量分数的高低决定,所述图像质量分数为:
QI=HFE/MFE
其中,QI为所述图像质量分数,HFE为在所述归一化虹膜图像序列的频谱中去除噪声频谱后的高频能量,MFE为在所述归一化虹膜图像序列的频谱中去除直流频谱后的低频能量。
3.如权利要求1所述的虹膜图像增强方法,其特征在于,所述将所述 伴生图像序列根据多方向汉明距离判别准则对照所述基准图像进行全局对齐,得到对齐伴生图像序列的步骤,进一步包含:
将所述伴生图像序列中的每帧图像沿四周方向扩展,得到扩展伴生图像序列;
在所述扩展伴生图像序列的中开一个窗口,所述窗口大小与所述伴生图像序列中的每帧图像大小相同;
对所述扩展伴生图像序列的每帧图像,逐行、逐列滑动所述窗口,计算所述窗口内图像与所述基准图像的汉明距离,将所述汉明距离最小的所述窗口内图像记为对齐图像,根据每帧图像的所述对齐图像得到所述对齐伴生图像序列。
4.如权利要求1所述的虹膜图像增强方法,其特征在于,所述上采样采用双立方插值法。
5.如权利要求1所述的虹膜图像增强方法,其特征在于,所述将各堆栈中的所述虹膜图像重建子块按照原始相对位置进行归一化拼接,进一步包含:
对各堆栈中的所述虹膜图像重建子块,利用归一化因子平滑拼接边缘,得到拼接归一化重建子块为:
RBe’=RBe/NRe
其中,RBe’为所述拼接归一化重建子块,RBe为所述虹膜图像重建子块,NRe为所述归一化因子。
6.如权利要求1~5任一项所述的虹膜图像增强方法,其特征在于,所述去噪处理采用中值联合高斯低通滤波,所述锐化处理采用拉普拉斯算子。
7.如权利要求1~5任一项所述的虹膜图像增强方法,其特征在于,所述下采样采用双线性插值法。
8.如权利要求1~5任一项所述的虹膜图像增强方法,其特征在于,所述归一化虹膜图像序列的图像帧数为10,单帧图像分别率为288*64。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810311065.6A CN108629744B (zh) | 2018-04-09 | 2018-04-09 | 一种图像增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810311065.6A CN108629744B (zh) | 2018-04-09 | 2018-04-09 | 一种图像增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108629744A CN108629744A (zh) | 2018-10-09 |
CN108629744B true CN108629744B (zh) | 2021-07-30 |
Family
ID=63704887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810311065.6A Active CN108629744B (zh) | 2018-04-09 | 2018-04-09 | 一种图像增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108629744B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020107145A1 (zh) * | 2018-11-26 | 2020-06-04 | 温州生物材料与工程研究所 | 一种虹膜识别装置和方法 |
CN110148121B (zh) * | 2019-05-09 | 2022-11-22 | 腾讯科技(深圳)有限公司 | 一种皮肤图像处理方法、装置、电子设备及介质 |
CN110674737B (zh) * | 2019-09-23 | 2022-04-19 | 北京无线电计量测试研究所 | 一种虹膜识别的增强方法 |
CN112887758B (zh) * | 2019-11-29 | 2023-04-14 | 北京百度网讯科技有限公司 | 视频处理方法和装置 |
CN111091566A (zh) * | 2020-02-04 | 2020-05-01 | 浙江省人民医院 | 一种糖尿病并发症视网膜病变数据检测方法 |
CN112967207B (zh) * | 2021-04-23 | 2024-04-12 | 北京恒安嘉新安全技术有限公司 | 一种图像处理方法、装置、电子设备及存储介质 |
CN114418927B (zh) * | 2021-11-09 | 2023-06-16 | 四川大学 | 一种基于空间关系特征匹配的闭环检测方法及系统 |
CN115049636A (zh) * | 2022-07-12 | 2022-09-13 | 郑州工商学院 | 一种计算机图像处理方法 |
CN117635506B (zh) * | 2024-01-24 | 2024-04-05 | 成都航天凯特机电科技有限公司 | 一种基于AI赋能Mean Shift算法的图像增强方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6084989A (en) * | 1996-11-15 | 2000-07-04 | Lockheed Martin Corporation | System and method for automatically determining the position of landmarks in digitized images derived from a satellite-based imaging system |
CN1426760A (zh) * | 2001-12-18 | 2003-07-02 | 中国科学院自动化研究所 | 基于活体虹膜的身份识别方法 |
EP1516493B1 (en) * | 2002-06-14 | 2006-10-11 | Koninklijke Philips Electronics N.V. | A method and system for optimizing image sharpness during coding |
CN102324031A (zh) * | 2011-09-07 | 2012-01-18 | 江西财经大学 | 老龄用户多生物特征身份认证中的隐式语义特征提取方法 |
CN102521227A (zh) * | 2011-10-18 | 2012-06-27 | 中国科学院自动化研究所 | 基于用户信息建模的图像标注增强方法 |
CN103025229A (zh) * | 2010-04-29 | 2013-04-03 | 麻省理工学院 | 适用于光学相干断层扫描技术的移动修正和图像增强的方法和装置 |
CN103020898A (zh) * | 2012-12-03 | 2013-04-03 | 北京理工大学 | 序列虹膜图像超分辨率重建方法 |
CN104599260A (zh) * | 2015-02-02 | 2015-05-06 | 天津三英精密仪器有限公司 | 一种基于双能谱和小波融合的x射线图像增强方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2016138608A (ru) * | 2016-09-29 | 2018-03-30 | Мэджик Лип, Инк. | Нейронная сеть для сегментации изображения глаза и оценки качества изображения |
-
2018
- 2018-04-09 CN CN201810311065.6A patent/CN108629744B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6084989A (en) * | 1996-11-15 | 2000-07-04 | Lockheed Martin Corporation | System and method for automatically determining the position of landmarks in digitized images derived from a satellite-based imaging system |
CN1426760A (zh) * | 2001-12-18 | 2003-07-02 | 中国科学院自动化研究所 | 基于活体虹膜的身份识别方法 |
EP1516493B1 (en) * | 2002-06-14 | 2006-10-11 | Koninklijke Philips Electronics N.V. | A method and system for optimizing image sharpness during coding |
CN103025229A (zh) * | 2010-04-29 | 2013-04-03 | 麻省理工学院 | 适用于光学相干断层扫描技术的移动修正和图像增强的方法和装置 |
CN102324031A (zh) * | 2011-09-07 | 2012-01-18 | 江西财经大学 | 老龄用户多生物特征身份认证中的隐式语义特征提取方法 |
CN102521227A (zh) * | 2011-10-18 | 2012-06-27 | 中国科学院自动化研究所 | 基于用户信息建模的图像标注增强方法 |
CN103020898A (zh) * | 2012-12-03 | 2013-04-03 | 北京理工大学 | 序列虹膜图像超分辨率重建方法 |
CN104599260A (zh) * | 2015-02-02 | 2015-05-06 | 天津三英精密仪器有限公司 | 一种基于双能谱和小波融合的x射线图像增强方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108629744A (zh) | 2018-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108629744B (zh) | 一种图像增强方法 | |
JP7542156B2 (ja) | 人物画像修復方法、装置、電子機器、記憶媒体及びプログラム製品 | |
CN108765343B (zh) | 图像处理的方法、装置、终端及计算机可读存储介质 | |
Zeyde et al. | On single image scale-up using sparse-representations | |
US20170365046A1 (en) | Algorithm and device for image processing | |
KR101795271B1 (ko) | 영상의 선명화를 위한 전처리를 수행하는 영상 처리 장치 및 방법 | |
CN111275626A (zh) | 一种基于模糊度的视频去模糊方法、装置及设备 | |
CN103020898B (zh) | 序列虹膜图像超分辨率重建方法 | |
Teng et al. | Adaptive morphological filtering method for structural fusion restoration of hyperspectral images | |
Li et al. | Image/video restoration via multiplanar autoregressive model and low-rank optimization | |
KR101428728B1 (ko) | 차량 번호 인식 장치 및 방법 | |
Peng et al. | Efficient image resolution enhancement using edge-directed unsharp masking sharpening for real-time ASIC applications | |
WO2009047208A1 (en) | Method of generating a multiscale contrast enhanced image | |
Rakhshanfar et al. | Low-frequency image noise removal using white noise filter | |
Xiong et al. | Image hallucination with feature enhancement | |
Zhao et al. | Single depth image super-resolution with multiple residual dictionary learning and refinement | |
KR101515686B1 (ko) | 주파수 성분과 영역 분할을 이용한 고해상도 얼굴 영상 복원 장치 및 방법 | |
CN110400270B (zh) | 一种利用图像分解和多重校正融合的车牌去雾方法 | |
He et al. | Joint motion deblurring and superresolution from single blurry image | |
Subathro et al. | Detecting digital image forgeries using resampling by automatic Region of Interest (ROI) | |
Arif et al. | Super-resolution using edge modification through stationary wavelet transform | |
Zhang et al. | Digital image forensics of non-uniform deblurring | |
Toan | Super-resolution method for reconstructing street images from surveillance system based on real-ESRGAN | |
KR101650897B1 (ko) | 저해상도 컨텐츠를 위한 윈도우 사이즈 주밍 방법 및 그 장치 | |
CN111652808A (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 |