iriscode虹膜识别MATLAB代码Libor Masek写的 2003 matlab代码 修改版本

这是matlab 2003年 比较早了 而且资源效果一般 入门不错 欢迎下载

% createiristemplate - generates a biometric template from an iris in
% an eye image.
%
% Usage: 
% [template mask] = createiristemplate（eyeimage_filename）
%
% Arguments:
%	eyeimage_filename   - the file name of the eye image
%
% Output:
%	template		    - the binary iris biometric template
%	mask			    - the binary iris noise mask
%
% Author: 
% Libor Masek
% masekl01@csse.uwa.edu.au
% School of Computer Science & Software Engineering
% The University of Western Australia
% November 2003


function [template mask] = createiristemplate（eyeimage_filename）

% path for writing diagnostic images
global DIAGPATH
DIAGPATH = ‘diagnostics‘;  %文件夹的名字  需要新建这个文件夹

%normalisation parameters
radial_res = 20;
angular_res = 240;
% with these settings a 9600 bit iris template is
% created

%feature encoding parameters
nscales=1;
minWaveLength=18;
mult=1; % not applicable if using nscales = 1
sigmaOnf=0.5;

eyeimage = imread（eyeimage_filename）; 

savefile = [eyeimage_filename‘-houghpara.mat‘];
[statmess]=fileattrib（savefile）;

if stat == 1
    % if this file has been processed before
    % then load the circle parameters and
    % noise information for that file.
    load（savefile）;
    
else
    
    % if this file has not been processed before
    % then perform automatic segmentation and
    % save the results to a file
    
    [circleiris circlepupil imagewithnoise] = segmentiris（eyeimage）;
    save（savefile‘circleiris‘‘circlepupil‘‘imagewithnoise‘）;
    
end

% WRITE NOISE IMAGE
%

imagewithnoise2 = uint8（imagewithnoise）;
imagewithcircles = uint8（eyeimage）;

%get pixel coords for circle around iris
[xy] = circlecoords（[circleiris（2）circleiris（1）]circleiris（3）size（eyeimage））;
ind2 = sub2ind（size（eyeimage）double（y）double（x））; 

%get pixel coords for circle around pupil
[xpyp] = circlecoords（[circlepupil（2）circlepupil（1）]circlepupil（3）size（eyeimage））;
ind1 = sub2ind（size（eyeimage）double（yp）double（xp））;


% Write noise regions
imagewithnoise2（ind2） = 255;
imagewithnoise2（ind1） = 255;
% Write circles overlayed
imagewithcircles（ind2） = 255;
imagewithcircles（ind1） = 255;
w = cd;
cd（DIAGPATH）;
imwrite（imagewithnoise2[eyeimage_filename‘-noise.jpg‘]‘jpg‘）;
imwrite（imagewithcircles[eyeimage_filename‘-segmented.jpg‘]‘jpg‘）;
cd（w）;

% perform normalisation

[polar_array noise_array] = normaliseiris（imagewithnoise circleiris（2）...
    circleiris（1） circleiris（3） circlepupil（2） circlepupil（1） circlepupil（3）eyeimage_filename radial_res angular_res）;


% WRITE NORMALISED PATTERN AND NOISE PATTERN
w = cd;
cd（DIAGPATH）;
imwrite（polar_array[eyeimage_filename‘-polar.jpg‘]‘jpg‘）;
imwrite（noise_array[eyeimage_filename‘-polarnoise.jpg‘]‘jpg‘）;
cd（w）;

% perform feature encoding
[template mask] = encode（polar_array noise_array nscales minWaveLength mult sigmaOnf）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-04-16 10:18  iriscode虹膜识别MATLAB代码Libor Masek\
     文件       12161  2019-04-15 22:45  iriscode虹膜识别MATLAB代码Libor Masek\001_1_1.jpg
     文件       71070  2019-04-15 23:00  iriscode虹膜识别MATLAB代码Libor Masek\001_1_1.jpg-houghpara.mat
     文件       90678  2003-05-24 22:02  iriscode虹膜识别MATLAB代码Libor Masek\002_1_1.jpg
     文件       92969  2019-04-16 09:21  iriscode虹膜识别MATLAB代码Libor Masek\002_1_1.jpg-houghpara.mat
     文件       90678  2003-05-24 22:02  iriscode虹膜识别MATLAB代码Libor Masek\003_1_1.jpg
     文件       79712  2019-04-16 09:25  iriscode虹膜识别MATLAB代码Libor Masek\003_1_1.jpg-houghpara.mat
     目录           0  2019-04-15 22:25  iriscode虹膜识别MATLAB代码Libor Masek\Matching\
     文件        1733  2003-11-28 15:51  iriscode虹膜识别MATLAB代码Libor Masek\Matching\gethammingdistance.m
     文件        1552  2003-11-28 15:59  iriscode虹膜识别MATLAB代码Libor Masek\Matching\shiftbits.m
     目录           0  2019-04-15 22:25  iriscode虹膜识别MATLAB代码Libor Masek\Normal_encoding\
     文件        1997  2003-11-30 13:01  iriscode虹膜识别MATLAB代码Libor Masek\Normal_encoding\encode.m
     文件        2557  2003-11-30 13:09  iriscode虹膜识别MATLAB代码Libor Masek\Normal_encoding\gaborconvolve.m
     文件        5035  2003-11-30 13:24  iriscode虹膜识别MATLAB代码Libor Masek\Normal_encoding\normaliseiris.m
     文件        1432  2003-11-30 13:44  iriscode虹膜识别MATLAB代码Libor Masek\README.txt
     文件       12282  2006-08-05 08:07  iriscode虹膜识别MATLAB代码Libor Masek\S1007R05.jpg
     文件       70674  2019-04-16 09:44  iriscode虹膜识别MATLAB代码Libor Masek\S1007R05.jpg-houghpara.mat
     目录           0  2019-04-15 23:05  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\
     文件        1472  2003-02-07 09:25  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\addcircle.m
     文件         844  2003-02-11 18:57  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\adjgamma.m
     文件        2219  2003-11-25 11:11  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\canny.m
     文件        1622  2003-11-30 12:20  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\circlecoords.m
     文件        2321  2003-11-25 09:42  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\findcircle.m
     文件        1090  2003-11-25 10:56  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\findline.m
     文件         917  2003-11-25 09:41  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\houghcircle.m
     文件        3304  2003-02-11 18:34  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\hysthresh.m
     文件         812  2003-11-25 09:23  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\linecoords.m
     文件        3683  2003-02-11 18:34  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\nonmaxsup.m
     文件        3261  2003-11-25 10:09  iriscode虹膜识别MATLAB代码Libor Masek\Segmentation\segmentiris.m
     文件        2965  2019-04-16 10:06  iriscode虹膜识别MATLAB代码Libor Masek\createiristemplate.m
     目录           0  2019-04-16 09:44  iriscode虹膜识别MATLAB代码Libor Masek\diagnostics\
............此处省略21个文件信息