POC 10_1

close all;
clearvars;
clc;

knee=imread("knee.png");

figure(1);
imshow(knee);

kx=363.0;
ky=163.0;

visited=zeros(779,693);
segmented=zeros(779,693);
stackk=zeros(2,10000);
stackk(1,1)=kx;
stackk(2,1)=ky;
visited(kx,ky)=1;
segmented(kx,ky)=1;

istack=1;

Q=4;

while(istack>0)

    x=stackk(1,istack);
    y=stackk(2,istack);
    istack=istack-1;
    if(x>1 && x<779 && y>1 && y<693)

    for i=(-1):2
        for j=(-1):2
            px=x+i;
            py=y+j;
            if(px>1 && px<779 && py>1 && py<693 && x~=px && y~=py && visited(px,py)==0)
                visited(px,py)=1;
                if(abs(knee(x,y)-knee(px,py))<Q)
                    segmented(px,py)=1;
                    istack=istack+1;
                    px=stackk(1,istack);
                    py=stackk(2,istack);
                end
            end



        end
    end

    end
end

figure(2);
imshow(segmented);