图像分割分水岭算法MATLAB源代码

f=imread（‘e:\pic\man.bmp‘）;                                   % 读入图像
bw=im2bw（fgraythresh（f））;                                 % 转换为黑白二值图像
bwc=~bw;                                                % 图像反色
dst=bwdist（bwc）;                                          % 图像距离
ws=watershed（-dst）;
w=ws==0;
rf=bw&~w;
figureimshow（f）;                                        % 显示原始图像
title（‘Original Image‘）;                                    % 设置图像标题
figureimshow（bw）;                                      % 显示处理后的图像
title（‘Negative Image‘）;                                   % 设置图像标题
figureimshow（ws）;                                       % 显示处理后的图像
title（‘Watershed - Distance Transform‘）;                     % 设置图像标题
figureimshow（rf）;                                        % 显示处理后的图像
title（‘Superimposed - Watershed and original image‘）;         % 设置图像标题
h=fspecial（‘sobel‘）;                                       % 设置滤波器 
fd=im2double（f）;                                         % 数据类型转换
sq=sqrt（imfilter（fdh‘replicate‘）.^2+imfilter（fdh‘