%?ADDCIRCLE
%
%?A?circle?generator?for?adding?（drawing）?weights?into?a?Hough?accumumator
%?array.
%
%?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;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-03-24?10:57??虹膜識別程序\
?????目錄???????????0??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\
?????文件???????90678??2003-05-24?22:15??虹膜識別程序\虹膜識別程序\001_1_1.bmp
?????文件???????77516??2018-03-24?12:58??虹膜識別程序\虹膜識別程序\001_1_1.bmp-houghpara.mat
?????文件???????90678??2003-05-24?22:15??虹膜識別程序\虹膜識別程序\001_1_2.bmp
?????文件???????78641??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\001_1_2.bmp-houghpara.mat
?????文件???????90678??2003-05-24?22:02??虹膜識別程序\虹膜識別程序\002_1_1.bmp
?????文件???????92969??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\002_1_1.bmp-houghpara.mat
?????文件????????1472??2003-02-07?09:25??虹膜識別程序\虹膜識別程序\addcircle.m
?????文件????????1241??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\adjgamma-line.jpg
?????文件????????2157??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\adjgamma.jpg
?????文件?????????844??2003-02-11?18:57??虹膜識別程序\虹膜識別程序\adjgamma.m
?????文件?????????906??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\bottomeyelid.jpg
?????文件????????1686??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\canny-line.jpg
?????文件????????1787??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\canny.jpg
?????文件????????2219??2003-11-25?11:11??虹膜識別程序\虹膜識別程序\canny.m
?????文件????????1622??2003-11-30?12:20??虹膜識別程序\虹膜識別程序\circlecoords.m
?????文件????????3032??2014-06-06?11:59??虹膜識別程序\虹膜識別程序\createiristemplate.m
?????目錄???????????0??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\
?????文件???????12391??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_1.bmp-noise.jpg
?????文件???????21410??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_1.bmp-normal.jpg
?????文件?????????851??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_1.bmp-polar.jpg
?????文件????????1011??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_1.bmp-polarnoise.jpg
?????文件???????13961??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_1.bmp-segmented.jpg
?????文件???????12403??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_2.bmp-noise.jpg
?????文件???????20952??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_2.bmp-normal.jpg
?????文件?????????872??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_2.bmp-polar.jpg
?????文件?????????957??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_2.bmp-polarnoise.jpg
?????文件???????14158??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\001_1_2.bmp-segmented.jpg
?????文件???????11882??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\002_1_1.bmp-noise.jpg
?????文件???????15664??2018-03-24?12:59??虹膜識別程序\虹膜識別程序\diagnostics\002_1_1.bmp-normal.jpg
............此處省略24個文件信息