Prev: surf function problem
Next: Best approximation
From: Marcus on 18 Dec 2009 10:59
ImageAnalyst <imageanalyst(a)mailinator.com> wrote in message <f7812da0-5d69-4b76-8a8b-0cad671f68c4(a)c34g2000yqn.googlegroups.com>...
> It's just a list of citations. You'd have to order those articles
> that you're interested in. Your university library could probably
> help you with that if you bring them the list of articles you want. I
> don't have any code or articles to send you and I don't work in the
> iris analysis field. Perhaps you could contact this guy who has
> previously shown interest in meeting other people who have an interest
> in irises:
> http://www.mathworks.com/matlabcentral/newsreader/author/96852
>
> If you feel like this is too stretching of a goal, given your current
> skill set, to get done within the alloted time then you should consult
> with your professor to get some project with a higher likelihood of
> success.

Thanks Sir for your advice by the way how to contact this person? He also a brilliant guys i think its hard to find him as well!
From: ImageAnalyst on 18 Dec 2009 11:07
On Dec 18, 10:59 am, "Marcus " <sbeng_...(a)yahoo.com> wrote:
> Thanks Sir for your advice by the way how to contact this person? He also a brilliant guys i think its hard to find him as well!
--------------------------------------------------------------------------------
Did you actually click the link I gave you?
If you did, you'd have your answer.
http://www.limbistraine.com/ro/software/popescu-bodorin-nicolaie-en.html
It has his phone, email, university, everything.
I can't help you beyond that.
From: Marcus on 18 Dec 2009 11:11
ImageAnalyst <imageanalyst(a)mailinator.com> wrote in message <97fd6438-4cad-4f39-ade7-3ef334b90318(a)e27g2000yqd.googlegroups.com>...
> On Dec 18, 10:59 am, "Marcus " <sbeng_...(a)yahoo.com> wrote:
> > Thanks Sir for your advice by the way how to contact this person? He also a brilliant guys i think its hard to find him as well!
> --------------------------------------------------------------------------------
> Did you actually click the link I gave you?
> If you did, you'd have your answer.
> http://www.limbistraine.com/ro/software/popescu-bodorin-nicolaie-en.html
> It has his phone, email, university, everything.
> I can't help you beyond that.

Ops Sorry dint see that page! Ok thanks Sir for ur help!
From: Marcus on 19 Dec 2009 08:08
Hi sir,

I face some difficulties! I can't manage to run this code it appear to come out this error: ??? Maximum recursion limit of 500 reached. Use set(0,'RecursionLimit',N)
to change the limit. Be aware that exceeding your available stack space can
crash MATLAB and/or your computer.

Error in ==> imformats>find_in_registry at 114
varargout{1} = find_in_registry(fmts, varargin{1});

The eye image: http://www.sizlopedia.com/wp-content/uploads/face-recognition.jpg
direct crop the eye part out.

TheCode:

function [y0detect,x0detect,Accumulator] = houghcircle(Imbinary,r,thresh,region)
Imbinary=imread('eyejpg');
houghcircle(i,10,4);
imSize = size(Imbinary);

if nargin > 2 & thresh < 4
error('treshold value must be bigger or equal to 4');
return
end
if nargin < 3 thresh = 4; end
if nargin < 4 region = [1,1,imSize(2),imSize(1)]; end

%Voting
Accumulator = zeros(imSize);
xmin = region(1); xmax = xmin + region(3) - 1;
ymin = region(2); ymax = ymin + region(4) - 1;
[yIndex xIndex] = find(Imbinary);
% loop through edge points
for cnt = 1:length(xIndex)
xCenter = xIndex(cnt); yCenter = yIndex(cnt);
xmx=xCenter-r; xpx=xCenter+r;
ypx=yCenter+r; ymx=yCenter-r;
% skip this point if circle is completely outside region
if xpx<xmin | xmx>xmax | ypx<ymin | ymx>ymax continue; end
if ypx<=imSize(1) & ymx>=1 & xpx<=imSize(2) & xmx>=1 %ypxin & ymxin & xpxin & xmxin % circle is completely inside region
x = 1;
y = r;
d = 5/4 - r;
Accumulator(ypx,xCenter) = Accumulator(ypx,xCenter) + 1;
Accumulator(ymx,xCenter) = Accumulator(ymx,xCenter) + 1;
Accumulator(yCenter,xpx) = Accumulator(yCenter,xpx) + 1;
Accumulator(yCenter,xmx) = Accumulator(yCenter,xmx) + 1;
while y > x
xpx = xCenter+x; xmx = xCenter-x;
ypy = yCenter+y; ymy = yCenter-y;
ypx = yCenter+x; ymx = yCenter-x;
xpy = xCenter+y; xmy = xCenter-y;
Accumulator(ypy,xpx) = Accumulator(ypy,xpx) + 1;
Accumulator(ymy,xpx) = Accumulator(ymy,xpx) + 1;
Accumulator(ypy,xmx) = Accumulator(ypy,xmx) + 1;
Accumulator(ymy,xmx) = Accumulator(ymy,xmx) + 1;
Accumulator(ypx,xpy) = Accumulator(ypx,xpy) + 1;
Accumulator(ymx,xpy) = Accumulator(ymx,xpy) + 1;
Accumulator(ypx,xmy) = Accumulator(ypx,xmy) + 1;
Accumulator(ymx,xmy) = Accumulator(ymx,xmy) + 1;
if d < 0
d = d + x * 2 + 3;
x = x + 1;
else
d = d + (x - y) * 2 + 5;
x = x + 1;
y = y - 1;
end
end
if x == y
xpx = xCenter+x; xmx = xCenter-x;
ypy = yCenter+y; ymy = yCenter-y;
Accumulator(ypy,xpx) = Accumulator(ypy,xpx) + 1;
Accumulator(ymy,xpx) = Accumulator(ymy,xpx) + 1;
Accumulator(ypy,xmx) = Accumulator(ypy,xmx) + 1;
Accumulator(ymy,xmx) = Accumulator(ymy,xmx) + 1;
end
else % circle is partly outside region - need boundary checking
ypxin = ypx>=ymin & ypx<=ymax;
ymxin = ymx>=ymin & ymx<=ymax;
xpxin = xpx>=xmin & xpx<=xmax;
xmxin = xmx>=xmin & xmx<=xmax;
if ypxin Accumulator(ypx,xCenter) = Accumulator(ypx,xCenter) + 1; end
if ymxin Accumulator(ymx,xCenter) = Accumulator(ymx,xCenter) + 1; end
if xpxin Accumulator(yCenter,xpx) = Accumulator(yCenter,xpx) + 1; end
if xmxin Accumulator(yCenter,xmx) = Accumulator(yCenter,xmx) + 1; end
x = 1;
y = r;
d = 5/4 - r;
while y > x
xpx = xCenter+x; xpxin = xpx>=xmin & xpx<=xmax;
xmx = xCenter-x; xmxin = xmx>=xmin & xmx<=xmax;
ypy = yCenter+y; ypyin = ypy>=ymin & ypy<=ymax;
ymy = yCenter-y; ymyin = ymy>=ymin & ymy<=ymax;
ypx = yCenter+x; ypxin = ypx>=ymin & ypx<=ymax;
ymx = yCenter-x; ymxin = ymx>=ymin & ymx<=ymax;
xpy = xCenter+y; xpyin = xpy>=xmin & xpy<=xmax;
xmy = xCenter-y; xmyin = xmy>=xmin & xmy<=xmax;
if ypyin & xpxin Accumulator(ypy,xpx) = Accumulator(ypy,xpx) + 1; end
if ymyin & xpxin Accumulator(ymy,xpx) = Accumulator(ymy,xpx) + 1; end
if ypyin & xmxin Accumulator(ypy,xmx) = Accumulator(ypy,xmx) + 1; end
if ymyin & xmxin Accumulator(ymy,xmx) = Accumulator(ymy,xmx) + 1; end
if ypxin & xpyin Accumulator(ypx,xpy) = Accumulator(ypx,xpy) + 1; end
if ymxin & xpyin Accumulator(ymx,xpy) = Accumulator(ymx,xpy) + 1; end
if ypxin & xmyin Accumulator(ypx,xmy) = Accumulator(ypx,xmy) + 1; end
if ymxin & xmyin Accumulator(ymx,xmy) = Accumulator(ymx,xmy) + 1; end
if d < 0
d = d + x * 2 + 3;
x = x + 1;
else
d = d + (x - y) * 2 + 5;
x = x + 1;
y = y - 1;
end
end
if x == y
xpx = xCenter+x; xpxin = xpx>=xmin & xpx<=xmax;
xmx = xCenter-x; xmxin = xmx>=xmin & xmx<=xmax;
ypy = yCenter+y; ypyin = ypy>=ymin & ypy<=ymax;
ymy = yCenter-y; ymyin = ymy>=ymin & ymy<=ymax;
if ypyin & xpxin Accumulator(ypy,xpx) = Accumulator(ypy,xpx) + 1; end
if ymyin & xpxin Accumulator(ymy,xpx) = Accumulator(ymy,xpx) + 1; end
if ypyin & xmxin Accumulator(ypy,xmx) = Accumulator(ypy,xmx) + 1; end
if ymyin & xmxin Accumulator(ymy,xmx) = Accumulator(ymy,xmx) + 1; end
end
end
end

% Finding local maxima in Accumulator
y0detect = []; x0detect = [];
AccumulatorbinaryMax = imregionalmax(Accumulator);
[Potential_y0 Potential_x0] = find(AccumulatorbinaryMax == 1);
Accumulatortemp = Accumulator - thresh;
for cnt = 1:length(Potential_y0)
if Accumulatortemp(Potential_y0(cnt),Potential_x0(cnt)) >= 0
y0detect = [y0detect;Potential_y0(cnt)];
x0detect = [x0detect;Potential_x0(cnt)];
end
end
First  |  Prev  | 
Pages: 1 2
Prev: surf function problem
Next: Best approximation