createiristemplate.m

精确定位 虹膜的内外边缘

function [t, m] = createiristemplate(eyeimage_filename)
% path for writing diagnostic images
global DIAGPATH
DIAGPATH = 'diagnostics';

%normalisation parameters
radial_res = 64;
angular_res =256;
% 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 = eyeimage_filename; 
%eyeimage = rgb2gray(eyeimage1);
[circleiris circlepupil imagewithnoise] = segmentiris(eyeimage);
   
% WRITE NOISE IMAGE
%
imagewithnoise2 = uint8(imagewithnoise);
imagewithcircles = uint8(eyeimage);
%get pixel coords for circle around iris
[x,y] = circlecoords([circleiris(2),circleiris(1)],circleiris(3),size(eyeimage));
ind2 = sub2ind(size(eyeimage),double(y),double(x)); 

%get pixel coords for circle around pupil
[xp,yp] = 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;


[polar_array noise_array] = normaliseiris(imagewithnoise, circleiris(2),...
    circleiris(1), circleiris(3), circlepupil(2), circlepupil(1), circlepupil(3),eyeimage_filename, radial_res, angular_res);
t=polar_array;
m=noise_array;

% WRITE NORMALISED PATTERN, AND NOISE PATTERN
%w = cd;
%cd(DIAGPATH);
%imwrite(polar_array,['D:\虹膜识别\53607937iriscode\diagnostics\','-polar.jpg'],'jpg');
%imwrite(noise_array,['D:\虹膜识别\53607937iriscode\diagnostics\','-polarnoise.
%jpg'],'jpg');
%cd(w);
% perform feature encoding
%[template mask] = encode(polar_array, noise_array, nscales, minWaveLength, mult, sigmaOnf);