虹膜识别+源代码+matlab+霍夫变换+hough变换+2

虹膜识别+源代码+matlab+霍夫变换+hough变换2

% ADDCIRCLE
%
% A circle generator for adding （drawing） weights into a Hough accumulator
% array.（添加一个圆到hough累加阵列中）
%
% Usage:  h = addcircle（h c radius weight）
% 
% Arguments:
%            h      - 2D accumulator array.
%            c      - [xy] coords of centre of circle.
%            radius - radius of the circle
%            weight - optional weight of values to be added to the
%                     accumulator array （defaults to 1）
%
% Returns:   h - Updated accumulator array.

% Peter Kovesi
% Department of Computer Science & Software Engineering
% The University of Western Australia
% April 2002

function h = addcircle（h c radius weight）

    [hr hc] = size（h）;
    
    if nargin == 3
	weight = 1;
    end
    
    % c and radius must be integers
    if any（c-fix（c））
	error（‘Circle centre must be in integer coordinates‘）;
    end
    
    if radius-fix（radius）
	error（‘Radius must be an integer‘）;
    end
    
    x = 0:fix（radius/sqrt（2））;
    costheta = sqrt（1 - （x.^2 / radius^2））;
    y = round（radius*costheta）;
    
    % Now fill in the 8-way symmetric points on a circle given coords 
    % [px py] of a point on the circle.
    
    px = c（2） + [x  y  y  x -x -y -y -x];
    py = c（1） + [y  x -x -y -y -x  x  y];

    % Cull points that are outside limits
    validx = px>=1 & px<=hr;
    validy = py>=1 & py<=hc;    
    valid = find（validx & validy）;

    px = px（valid）;
    py = py（valid）;
    
    ind = px+（py-1）*hr;
    h（ind） = h（ind） + weight;

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

     文件     308278  2004-11-08 22:29  0033_006.bmp

     文件     261919  2008-12-19 19:54  0023_006.bmp-houghpara.mat

     文件     251373  2008-12-19 10:03  0025_001.bmp-houghpara.mat

     文件     250228  2008-12-19 19:53  0025_006.bmp-houghpara.mat

     文件     240812  2008-12-19 10:04  0025_002.bmp-houghpara.mat

     文件     240603  2008-12-19 20:55  0033_006.bmp-houghpara.mat

     文件      41058  2008-12-16 19:45  testP.jpg

     文件       5407  2008-12-20 11:31  normaliseiris.m

     文件       3885  2008-12-19 19:42  createiristemplate.m

     文件       3809  2008-12-12 21:06  hough.txt

     文件       3699  2008-12-09 16:58  nonmaxsup.m

     文件       3684  2008-12-19 10:03  segmentiris.m

     文件       3340  2008-12-09 19:39  hysthresh.m

     文件       2817  2008-12-19 16:00  writeoriginal.m

     文件       2618  2008-12-17 15:56  gaborconvolve.m

     文件       2499  2008-12-19 09:59  findcircle.m

     文件       2490  2008-12-17 10:42  encode.m

     文件       2350  2008-12-19 16:57  Untitled.m

     文件       2180  2008-12-09 21:45  canny.m

     文件       1780  2008-12-09 19:45  gethammingdistance.m

     文件       1709  2008-12-20 09:57  final2.m

     文件       1691  2008-12-09 19:59  circlecoords.m

     文件       1622  2008-12-09 16:56  shiftbits.m

     文件       1503  2008-12-09 20:13  addcircle.m

     文件       1432  2003-11-30 13:44  README.txt

     文件       1420  2008-12-19 15:40  findline.m

     文件       1202  2009-02-25 10:09  说明.txt

     文件        979  2008-12-09 19:40  houghcircle.m

     文件        864  2008-12-09 20:00  adjgamma.m

     文件        847  2008-12-09 19:27  linecoords.m

............此处省略7个文件信息