matlab实现实现分水岭算法

分水岭算法是一种图像分割算法，可以将图像分割成多个具有明显边界的区域。下面是用 MATLAB 实现分水岭算法的步骤：

1. 读取图像并转换为灰度图像。

img = imread('example.jpg');
gray = rgb2gray(img);

2. 对灰度图像进行梯度计算。

hy = fspecial('sobel'); % 竖直方向的Sobel算子
hx = hy';
Iy = imfilter(double(gray), hy, 'replicate');
Ix = imfilter(double(gray), hx, 'replicate');
gradmag = sqrt(Ix.^2 + Iy.^2);

3. 对梯度图像进行阈值分割，得到二值图像。

level = graythresh(gradmag);
bw = imbinarize(gradmag, level);

4. 对二值图像进行形态学操作，填充空洞和平滑边界。

se = strel('disk', 20);
closedbw = imclose(bw, se);

5. 计算距离变换，得到每个像素点到最近背景像素点的距离。

D = bwdist(~closedbw);

6. 对距离变换进行分水岭变换。

D = -D;
D(~closedbw) = -Inf;
L = watershed(D);

7. 将分割结果可视化。

rgb = label2rgb(L, 'jet', 'w', 'shuffle');
imshow(rgb)

以上就是用 MATLAB 实现分水岭算法的步骤，其中包括了图像预处理、阈值分割、形态学操作、距离变换和分水岭变换等步骤。