梯度分水岭分割法 MATLAB编程，对附件的核磁图像 MRI.png 进行分水岭法分割.要求显示原图、梯 度图像、梯度分水岭分割图像、分水岭脊线、经闭-开运算平滑后的梯度分水岭分 割图像。

梯度分水岭分割是一种基于图像梯度的边缘检测方法，常用于医学图像分析，如MRI扫描中区分组织边界。在MATLAB中，你可以按照以下步骤操作：

1. **加载图像**：

   img = imread('MRI.png'); % 加载MRI图像
   figure, imshow(img), title('原始图像');

2. **灰度处理**：
如果MRI是彩色的，可以转换为灰度：

   img_gray = rgb2gray(img);
   figure, imshow(img_gray), title('灰度图像');

3. **计算梯度**：

   [gradient, edge] = imgradientmag(img_gray); % 计算梯度幅值和方向

4. **构造种子点**：
分水岭算法通常从局部最大像素开始，创建一个标记数组（种子矩阵）。

5. **分水岭变换**：

   markers = bwdist(edge, 'inf'); % 使用距离变换找到种子点
   ws = watershed(-gradient, markers); % 进行分水岭分割
   figure, imshow(ws), title('梯度分水岭分割图像');

6. **绘制分水岭脊线**：

   watershed_lines = bwboundaries(ws); % 提取分水岭轮廓
   plot(watershed_lines), hold on;

7. **闭-开运算平滑**：
这是为了去除小的断裂和噪声：

   ws_smoothed = imclose(ws, strel('disk', 5)); % 使用5x5圆盘结构元素进行膨胀然后腐蚀
   figure, imshow(ws_smoothed), title('平滑后的梯度分水岭分割图像');

8. **显示结果**：
最后，展示所有步骤的结果。

记得运行以上代码之前，确保你的MATLAB环境安装了`imgradientmag`, `bwmorph`等函数，如果没有，需要先安装它们。