I = imread('pout.tif');
figure
nexttile
imshow(I)%
se = strel('disk',3);
J=imbothat(I,se);%
nexttile
imshow(J)
J3=imtophat(I,se);%
nexttile
imshow(J3)
JJ=imsubtract(I,J);%
nexttile
imshow(JJ);
J2=imadd(I,imtophat(I,se));%
nexttile
imshow(J2)
J4 = imsubtract(imadd(I,imtophat(I,se)),imbothat(I,se));%
nexttile
imshow(J4)