function [circles,cImg] = findCircles(img)

edges = edgeDetect(img);

[rows, cols] = size(edges);

Rm = round(min(rows, cols) / 2);
Count = zeros(rows, cols, Rm);

for colIndex = 1 : cols
    for rowIndex = 1 : rows
        if(edges(rowIndex, colIndex) == 1)
            for Cx = 1 : cols
                for Cy = 1 : rows
                    radius = round(sqrt((colIndex - Cx)^2 + (rowIndex - Cy)^2));
                    if(radius > 0 && radius < Rm)
                        Count(Cx, Cy, radius) = Count(Cx, Cy, radius) + 1;
                    end
                end
            end
        end
    end
end

threshold = max(max(max(Count))) * 0.85; % there is a good result when it's 60

result = [];

for Cx = 1 : cols
    for Cy = 1 : rows
        for radius = 1 : min(rows, cols)

            %threshold = pi * radius;

            if(radius > 0 && radius < Rm && Count(Cx, Cy, radius) > threshold)
                if(isMaximal(Count, Cx, Cy, radius))
                    result = [result ; Cx, Cy, radius];
                end
            end
        end
    end
end

circles = filterCircles(result);

%Print the list into the screen.
printList(circles);

%Mark the circles in the original image.
cImg = addCircles(img, circles);

end


function filteredList = filterCircles(circles)

[rows, cols] = size(circles);
filteredList = [];
counter = 1;

for r1 = 1 : rows
   for r2 = 1 : rows
       if(r1 ~= r2)
           x1 = circles(r1, 1);
           y1 = circles(r1, 2);
           x2 = circles(r2, 1);
           y2 = circles(r2, 2);
           distance = sqrt((x1 - x2)^2 + (y1 - y2)^2);
           if(distance < 5)
             %  if(y1 > y2)
                %filteredList = [filteredList ; x1, y1, circles(r1, 3)];
                %counter = counter + 1;
                    circles(r2, 1) = 0;
                    circles(r2, 2) = 0;
                    circles(r2, 3) = 0;
             %  else
             %       circles(r1, 1) = 0;
            %        circles(r1, 2) = 0;
            %        circles(r1, 3) = 0;
            %   end
           end
       end
   end
end

for r = 1 : size(circles, 1)
    if(circles(r, 3) ~= 0)
        filteredList = [filteredList ; circles(r, 1), circles(r, 2), circles(r, 3)];
    end
end

end


function [] = printList(circles)

rows = size(circles, 1);

for r = 1 : rows
    fprintf('Circle %d: %d, %d, %d\n',r ,circles(r, 1), circles(r, 2), circles(r, 3));
end
end

function answer = isMaximal(Count, Cx, Cy, radius)

currentValue = Count(Cx, Cy, radius);
[rows, cols, depth] = size(Count);

if(Cx > 1 &&Count(Cx - 1, Cy, radius) > currentValue)
    answer = 0;
elseif (Cx < cols && Count(Cx + 1, Cy, radius) > currentValue)
    answer = 0;
elseif(Cy > 1 && Count(Cx, Cy - 1, radius) > currentValue)
    answer = 0;
elseif(Cy < rows && Count(Cx, Cy + 1, radius) > currentValue)
    answer = 0;
elseif(radius > 1 && Count(Cx, Cy, radius - 1) > currentValue)
    answer = 0;
elseif(radius < depth && Count(Cx, Cy, radius + 1) > currentValue)
    answer = 0;
else
    answer = 1;
end
end


function imageWithCircles = addCircles(img, circles)

rows = (size(circles, 1));

if(rows > 0)

    circleParameters = circles(1, :);
    imageWithCircles = insertShape(img,'circle',circleParameters,'LineWidth',2);

    for r = 2 : rows
        circleParameters = circles(r, :);
        imageWithCircles = insertShape(imageWithCircles,'circle',circleParameters,'LineWidth',2);
    end
end
end