afm=imread('watershed.tif');
figure,imshow(afm);
se = strel ('disk', 15) ;
Itop = imtophat (afm , se) ;
%figure , imshow( Itop , [ ]) , title ( 'top - hat im2
%age');
Ibot = imbothat (afm , se) ; 
figure , imshow( Ibot , [ ]) , title ('bottom - hat image');
Ienhance = imsubtract (imadd ( Itop , afm) , Ibot) ; 
figure , imshow( Ienhance) , title ('original + top - hat - bottom - hat');
Iec = imcomplement ( Ienhance) ;
figure , imshow( Iec) , title ('complementof enhanced image');
Iemin = imextendedmin( Iec , 22) ; 
figure , imshow( Iemin) , title ('extend2ed minima image’');
Iimpose = imimposemin ( Iec , Iemin) ; 
figure , imshow( Iimpose) , title ( 'imposed minima image');
%18 = bwperim( Iimpose) ; imshow(18);
BW = watershed(Iimpose);
rgb1 = label2rgb(BW);
figure,imshow(rgb1);