matlab-分水岭算法处理图像分割(用分水岭算法对米粒进行分割)

个人邮箱:570664186@qq.com,如有疑问请备注CSDN

代码实现：

f=imread('C:\Users\lenovo\Desktop\3(1).bmp');     % 读入图像
figure,imshow(f);        % 显示原始图像
title('原始图像'); 
bw=imbinarize(f);    % 转换为黑白二值图像
figure;
imshow(bw);       
title('二值图像');      
bwc=~bw;    % 图像反色，黑白反转
dst=bwdist(bwc);     % 图像距离
ws=watershed(-dst);   % 分水岭算法
figure;
imshow(ws);                                      
title('分水岭算法后距离变换');                  
w=ws==0;
rf=bw&~w;  %原图灰度图与分水岭算法处理后的图进行重叠
figure,imshow(rf);                        
title('分水岭算法和原始图像重叠后的图像');         

h=fspecial('sobel');    % 获得纵方向的sobel算子
fd=im2double(f);  % 数据类型转换-》双精度浮点型
%使用梯度算子进行计算
sq=sqrt(imfilter(fd,h,'replicate').^2+imfilter(fd,h','replicate').^2);  
sqoc=imclose(imopen(sq,ones(3,3)),ones(3,3));  % 图像闭合运算
wsd=watershed(sqoc);     %分水岭运算
figure;
imshow(wsd);       
title('梯度图像');            

wg=wsd==0;
rfg=f;
rfg(wg)=255;
figure;
imshow(rfg);                            
title('分水岭图像和原始图像叠加后图像');      
im=imextendedmin(f,18);   %用来修改图像，使其只是在特定的要求位置有局部极小。
%这里可以使用imimposemin来修改梯度幅值图像
Lim=watershed(bwdist(im));
figure,imshow(Lim);                          
title('分水岭标志图像');                 
em=Lim==0;
rfmin=imimposemin(sq,im|em);
wsdmin=watershed(rfmin);
figure,imshow(rfmin);                  
title('流域-梯度和标记控制');         
rfgm=f;
rfgm(wsdmin==0)=255;
figure,imshow(rfgm);       
title('分水岭和原始图像叠加后图像'); 

参考博主：

基于Matlab的标记分水岭分割算法
分水岭算法的原理、改进及matlab实现
分水岭算法 matlab的三种实现方法

感觉文章不错的同学麻烦动动小手点点关注订阅呗,您的肯定是对我持续更新最大的支持!