afm = imread('afmsurf.tif');figure, imshow(afm);
se = strel('disk', 15);
Itop = imtophat(afm, se);  % 高帽變換
Ibot = imbothat(afm, se);  % 低帽變換
figure, imshow(Itop, []);   % 高帽變換，體現原始圖像的灰度峰值
figure, imshow(Ibot, []);   % 低帽變換，體現原始圖像的灰度谷值
Ienhance = imsubtract(imadd(Itop, afm), Ibot);% 高帽圖像與低帽圖像相減，增強圖像
figure, imshow(Ienhance);
Iec = imcomplement(Ienhance); % 進一步增強圖像
Iemin = imextendedmin(Iec, 20); figure,imshow(Iemin) % 搜索Iec中的谷值
Iimpose = imimposemin(Iec, Iemin);
wat = watershed(Iimpose);  % 分水嶺分割
rgb = label2rgb(wat); figure, imshow(rgb); %  用不同的顏色表示分割出的不同區域