用matlab做图像分割的分水岭分割，我要代码

以下是MATLAB实现分水岭分割的代码示例：

% 读取图像
I = imread('image.png');

% 转换为灰度图像
Igray = rgb2gray(I);

% 通过 Sobel 滤波器计算水平和垂直梯度
hy = fspecial('sobel');
hx = hy';

Iy = imfilter(double(Igray), hy, 'replicate');
Ix = imfilter(double(Igray), hx, 'replicate');

% 计算梯度幅值和方向
gradmag = sqrt(Ix.^2 + Iy.^2);
gradmag = gradmag / max(gradmag(:));

% 设置分水岭算法的参数
se = strel('disk', 20);
Io = imopen(Igray, se);
Ie = imerode(Io, se);
Iobr = imreconstruct(Ie, Igray);
Iobrd = imdilate(Iobr, se);
Iobrcbr = imreconstruct(imcomplement(Iobrd), imcomplement(Iobr));
Iobrcbr = imcomplement(Iobrcbr);
fgm = imregionalmax(Iobrcbr);
D = bwdist(bw);
DL = watershed(D);
bgm = DL == 0;

% 计算分割结果
gradmag2 = imimposemin(gradmag, bgm | fgm);
L = watershed(gradmag2);

% 显示结果
figure;
imshow(I);
hold on;
h = imshow(label2rgb(L, 'jet', 'w'));
set(h, 'AlphaData', 0.5);

这段代码实现了分水岭分割，可以用于图像分割。需要注意的是，这段代码需要一些参数调整才能适用于不同的图像，比如 `strel` 函数的参数和 `imregionalmax` 函数的阈值等。