指纹图像方向场快速实现代码

指纹图像方向场快速实现代码 matlab实现 卷积实现 效果最佳

function [f_anglef_grad] = FingerprintOrient_convolution_f（fr1）
%%%%%%%%%%%%% function description %%%%%%%%%%%%%%%%%%%%%%%%
%%%%input    f:Array of img  ;  r1:window of average that control
%%%%relationship between the pix and its neighbors
%%%%
%%%%output   f_angle:array that is to show Orientation field   ;
%%%%         f_grad :array that is to show Grad field 
%%%%                            @author  Faith 15.12.22
%%%%                                    v2.0    16.8.3
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%% 读入图片取出接口 %%%%%%%%%%%
f1 = f;
r = r1;

%%%%获取相关图片信息%%%%%%
f2 = im2double（f1）;
% f2 = f1;
[heightwidth] = size（f1）;

EPI	= 57.29578;     %%%%%360/2pi，角度和弧度的换算

%%%%%%%%初始化矩阵%%%%%%%%%%%%%%%%
f_temp = f2;
f_angle = zeros（heightwidth）;
Vx = zeros（heightwidth）;
Vy = zeros（heightwidth）;
f_grad =  zeros（heightwidth）;
%%%%%%%%%%%%%%%%%%%%%%%%%soble算子%%%%%%%%%%%%%%%%%%%%%%%%%%
% soble_x = [-1 0 1; -2 0 2; -1 0 1];
% soble_y = [1 2 1; 0 0 0; -1 -2 -1];

%%%%%%%%%%%%%% 整体像素遍历 %%%%%%%%%%%%%%%

%%%%%%%向量法%%%%%%%%%%%%%%%
f_linex = f_temp（:）;
f_tempy = f_temp‘;
f_liney = f_tempy（:）;

soble_x = [-1 -2 -1 zeros（12*height -3） 1 2 1];
soble_y = [1 2 1 zeros（12*width-3） -1 -2 -1];

f_num = size（f_linex1）;
soblex_num = size（soble_x2）;
sobley_num = size（soble_y2）;

%%%%%%%%%%%%%%%%%%基础乘加%%%%%%%%%%%%%%%%%
% for i = 1:（f_num -soblex_num + 1）
%     f_lineVx（i） =  f_linex（i:i+length（soble_x）-1）‘ * soble_x‘;
% end
% for i = 1:（f_num -sobley_num + 1）
%     f_lineVy（i） =  f_liney（i:i+length（soble_y）-1）‘ * soble_y‘;
% end
% 
% 
% f_x = [zeros（1（soblex_num-1）/2） f_lineVx zeros（1（soblex_num-1）/2）];
% f_y = [zeros（1（sobley_num-1）/2） f_lineVy zeros（1（sobley_num-1）/2）];
%%%%%%%%%%

%%%%%%%%%%%%%%%% conv 卷积%%%%%%%%%%%%%%%%%%
tic

f_lineVx = -conv（f_linex‘soble_x）;
f_lineVy = -conv（f_liney‘soble_y）;

f_x = f_lineVx（1+（soblex_num-1）/2: length（f_lineVx）-（soblex_num-1）/2）;
f_y = f_lineVy（1+（sobley_num-1）/2: length（f_lineVy）-（sobley_num-1）/2）;

%%%%%%%%%%%%

%%%%%%%%%回复图像相应大小%%%%%%%%%%
f_xx = reshape（f_xheightwidth）;
f_yy = reshape（f_ywidthheight）;
f_yy = f_yy‘;

toc
%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%原遍历法%%%%%%%%%%%%%%%%%%%%%%%%%%%
tic

for x = 1:height
    for y = 1: width
         lvx = 0;
         lvy = 0;
         gradsum = 0;
         num = 0;
            %%%%%%%%%%%%%%%%%%%r * r 窗下的处理%%%%%%%%%%%%%%%%%%%
            for i = -r:r
                %%%%%%%%%%%%%%%%%% 边界处理 %%%%%%%%%%%%%%%%%%%%%
                if（x + i -1 < 1 || x + i + 1 > height ）
                    continue;
                end
                    for j = -r:r
                        if （y + j -1 < 1 || y + j + 1 > width）
                            continue;
                        end
                        
                        %%%%%%%%%%%%%%% soble 算子 x方向 和 y方向 的偏导
                        %%%%%%%%%%%%%%% x -1 0 1   y   1  2  1
               

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        4198  2016-08-03 12:43  FingerprintOrient_convolution_f.m