division.rar_matlab 灰度图像 分割_图像阈值分割

afm = imread('E:\3.bmp'); figure, imshow(afm); se = strel('disk', 15); Itop = imtophat(afm, se); % 高帽变换 Ibot = imbothat(afm, se); % 低帽变换 figure, imshow(Itop, []); title('高帽变换，体现原始图像的灰度峰值');% 高帽变换，体现原始图像的灰度峰值 figure, imshow(Ibot, []); title('低帽变换，体现原始图像的灰度谷值');% 低帽变换，体现原始图像的灰度谷值 Ienhance = imsubtract(imadd(Itop, afm), Ibot);% 高帽图像与低帽图像相减，增强图像 figure, imshow(Ienhance); title('高帽图像与低帽图像相减，增强图像'); Iec = imcomplement(Ienhance); % 进一步增强图像 Iemin = imextendedmin(Iec, 20); figure,imshow(Iemin);% 搜索Iec中的谷值 title('进一步增强图像'); Iimpose = imimposemin(Iec, Iemin); wat = watershed(Iimpose); % 分水岭分割 figure, imshow(wat); title('分水岭分割');