function [ending,bf] = ExtractMin(min_im)%is thin image
[r c]=size(min_im);
ending=zeros(200,2);
bf=zeros(100,2);
numending=1;
numbf=1;
%Location=1;% Actual location in structure
for i=20:(r-20)
for j=20:(c-20)
if (min_im(i,j) == 0)
P1=min_im(i,j+1);
P2=min_im(i-1,j+1);
P3=min_im(i-1,j);
P4=min_im(i-1,j-1);
P5=min_im(i,j-1);
P6=min_im(i+1,j-1);
P7=min_im(i+1,j);
P8=min_im(i+1,j+1);
minusP12 = xor(P1,P2);
minusP23 = xor(P2,P3);
minusP34 = xor(P3,P4);
minusP45 = xor(P4,P5);
minusP56 = xor(P5,P6);
minusP67 = xor(P6,P7);
minusP78 = xor(P7,P8);
minusP81 = xor(P8,P1);
NODEID = (double(minusP12) + double(minusP23)...
+ double(minusP34) + double(minusP45)...
+ double(minusP56) + double(minusP67)...
+ double(minusP78) + double(minusP81))/2 ;
if (NODEID == 1) % Ending
ending(numending,1) = i;
ending(numending,2) = j;
%Save in the structure
%MinVector1(Location).CenterRow=i;
%MinVector1(Location).CenterCol=j;
%MinVector1(Location).CenterType=1 ;%1 represent Ending
numending=numending+1;
%Location=Location+1;
end
if (NODEID == 3 ) % Bifurcation
bf(numbf,1) = i;
bf(numbf,2) = j;
%Save in the structure
%MinVector1(Location).CenterRow=i;
%MinVector1(Location).CenterCol=j;
%MinVector1(Location).CenterType=2 ;%2 represent Ending
numbf=numbf+1;
%Location=Location+1;
end
end
end
end
end %function