虹膜滤波归一化特征编码

可以实现对原始图像的去噪声、增强和归一化处理，生成多种gabor滤波器，用词滤波器作用于归一化虹膜，可以提取其特征编码。

% 批量计算Gabor编码，并保存在recd中
% 相关文件：mygabor.m，gabor_encode.mget_all_gaborcoef.mgabor_eer.m
function [recd result2 frr far]= gabor_all（）;
gabor = get_all_gaborcoef（）;
tic;
n=1;
for i=0:59
    for j=0:19
%         if （n<=600）%只试验前600幅图
            if （i<=9&j<=9）
                name = strcat（‘000‘int2str（i） ‘_00‘ int2str（j）‘.bmp‘）;
                code = gabor_encode（namegabor）;
                %recd（1n） = {name};
                recd（i+1j+1） = {code};
                n=n+1;
            elseif （i<=9&j>=10）
                name = strcat（‘000‘int2str（i） ‘_0‘ int2str（j） ‘.bmp‘）;
                code = gabor_encode（namegabor）;
                %recd（1n） = {name};
                recd（i+1j+1） = {code};
                n=n+1;
            elseif （i>=10&j<=9） 
                name = strcat（‘00‘int2str（i） ‘_00‘ int2str（j） ‘.bmp‘）;
                code = gabor_encode（namegabor）;
                %recd（1n） = {name};
                recd（i+1j+1） = {code};
                n=n+1;
            elseif （i>=10&j>=10）
                name = strcat（‘00‘int2str（i） ‘_0‘ int2str（j） ‘.bmp‘）;
                code = gabor_encode（namegabor）;
                %recd（1n） = {name};
                recd（i+1j+1） = {code};
                n=n+1;
            end
%         end
    end
end
average_encode_time = toc/1200
%----------------------------------------------
%hamming
for j=1:60      %第j人
    for i=j:60
        for k=1:20  %第k图
            for l=k:20
                template = recd{jk};   %第j人，第k图
                putin = recd{il};      %第i人，第l图
                result（kl）=mean（xor（templateputin））;                       %Hamming distance
%                 result（kl） = mean2（abs（template - putin））;                     %Manhatten distance
%                 [mm nn] = size（putin）;
%                 result（kl） = sqrt（sum（sum（（template - putin）.^2）））/（mm*nn）;  %Euclidian distance
            end
        end
        result2{ji}=result;
    end
end
%--------------------------------------------
%eer
thr=0:0.01:1;
numthr = numel（thr）;
for t=1:numthr;
thrd=thr（t）;
[frr（t）far（t）]=gabor_eer（thrdresult2）;
end
plot（thrfrrthrfar）
figure;plot（farfrr）
ee=0:0.01:1;hold on;plot（eeee）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      90678  2004-08-13 00:00  虹膜的特征提取算法\001_1_1.bmp

     文件       2312  2005-12-09 11:53  虹膜的特征提取算法\gabor_all.m

     文件      43195  2005-12-08 12:54  虹膜的特征提取算法\gabor_coef_data.m

     文件       1323  2005-12-09 12:25  虹膜的特征提取算法\gabor_eer.m

     文件       1045  2005-12-09 07:55  虹膜的特征提取算法\gabor_encode.m

     文件        371  2005-12-08 14:10  虹膜的特征提取算法\gabor说明.txt

     文件        278  2005-12-09 07:52  虹膜的特征提取算法\get_all_gaborcoef.m

     文件       1650  2006-04-10 02:52  虹膜的特征提取算法\mygabor.m

     目录          0  2009-02-03 16:11  虹膜的特征提取算法

----------- ---------  ---------- -----  ----

               140852                    9