指纹增强程序Hong_enhancement

本算法是基于Lin Hong et al 的论文“Fingerprint ImageEnhancement: Algorithm and Performance Evaluation”编写而成。其中一个重要的假设就是：沿脊线垂直方向的灰度变化呈正弦波形。

程序的流程如下所示：

0 Notation

O 是方向场F 是频率场 R 是区域标记

1 Normalization

将原始图像的灰度值的均值和方差调整到所期望的均值和方差，减少沿脊和谷方向上的灰度级的变化。

img=imread('25_2.tif');
img=double(img);
figure,imshow(img./255)
[ht,wt] =   size(img);
m       =   mean(img(:));
v       =   var(img(:));
gmidx   =   find(img > m);
lmidx   =   find(img <=m); m0=100;v0=100;
nimg(gmidx) = m0 + sqrt((v0*(img(gmidx)-m).^2)/v);
nimg(lmidx) = m0 - sqrt((v0*(img(lmidx)-m).^2)/v);
nimg        = reshape(nimg,[ht,wt]);
figure,imshow(nimg./255)

2 Orientation Image

方向场反映了指纹图像上纹线的方向，其准确性直接影响图像增强的效果。

（1）：将图像划分为不重叠的子块（16×16）

（2）：利用Sobel算子计算每个子块的像素点的梯度值。

（3）：利用以下公式计算中心点在（i，j）的子块的脊线的方向值

其中是以像素θ(i,j)为中心的子块的局部脊线方向。由于指纹脊线方向变化缓慢，并在一个小范围内具有相对稳定的变化趋势。因此可采用高斯低通滤波器进行平滑处理。

dbg_show_orientation=1;
img=imread('6_2.tif');
img=double(img);
% figure,imshow(img./255)
hy  = fspecial('sobel');
hx  = hy';
gx = imfilter(img,hx,'same','symmetric');
gy = imfilter(img,hy,'same','symmetric');
g   = gx+i*gy;
gblk=blkproc(g.^2,[10 10],inline('sum(sum(x))'));
oimg    = 0.5*angle(gblk)+pi/2;
[blkht,blkwt]   =   size(oimg)
%smoothen_orientation_field
g =   cos(2*oimg)+i*sin(2*oimg);
g =   imfilter(g,fspecial('gaussian',5));
oimg=   0.5*angle(g);
if(dbg_show_orientation)
quiver(cos(oimg),sin(oimg));
axis ij;axis image;
axis([0 blkwt 0 blkht]);
title('Orientation Image'); pause;
end;

3 Ridge Frequency Image

在指纹图像的局部非奇异区域里，沿着垂直于脊线的方向看，指纹脊线和谷线像素点灰度值大致形成一个二维的正弦波，定义纹线频率为相邻的两个波峰或波谷之间的像素点数的倒数。

（1）：将图像划分为不重叠的子块（16×16）

（2）：以图像子块中心点(i,j)为中心，子块脊线方向为短轴，作一个尺寸为的长方形窗口，如图所示：在窗口中按如下公式计算幅值X[K]