%Part 1:Compare various edge detectors
x=imread('fp.bmp');

e1=edge(x,'sobel');
e2=edge(x,'canny');
figure(1);

subplot(1,2,1)
imshow(e1);
title 'sobel edge detector';
subplot(1,2,2)
imshow(e2);
title 'canny edge detector';


%Part 2:applying morphological filtering to binarized image
th=graythresh(x);
y=im2bw(x,th);
b1=bwmorph(y,'thin');
b2=bwareaopen(y,10);
figure(2);
subplot(1,2,1)
imshow(x);
subplot(1,2,2)
imshow(b2);

%Part 3:Minutiae Extraction
b3=minutiae_extraction(b2);
figure(3)
 imshow(b3);
% check the validity of your results
figure(4)
imshow(x);hold on;
[i,j]=find(b3==1)
plot(j,i,'o');