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

梯度分水岭分割是一种基于图像梯度的边缘检测算法，在MATLAB中可以这样实现：

1. **加载图像**：
首先，你需要使用`imread`函数读取MRI.png图像：

   img = imread('MRI.png');

2. **预处理**：
对图像进行灰度化并计算梯度，这通常通过`imgradient`函数完成：

   [Ix, Iy] = imgradient(img);
   gradMag = sqrt(Ix.^2 + Iy.^2);

3. **分水岭变换**：
使用`imregionalmin`找出局部最小值作为种子点，然后用`watershed`函数进行分割：

   seeds = bwareaopen(imregionalmin(gradMag), 50); % 可调整阈值
   labels = watershed(-gradMag, seeds);

4. **展示结果**：
- 显示原始图像：`imshow(img)`
- 梯度图像：`imshow(gradMag, 'XData', [-1, 1], 'YData', [-1, 1])`
- 分水岭分割图像：`imshow(labels)`
- 脊线可以用`bwboundaries`函数获取，但MATLAB的这个功能需要额外安装某些工具箱：

     lines = bwboundaries(labels);
     plot(lines)

5. **闭-开操作**：
对于平滑效果，你可以使用`imerode`和`imdilate`组合进行闭-开（Closing）操作：

   structuringElement = strel('disk', 2); % 创建结构元素
   closedLabels = imclose(imbwlabel(labels), structuringElement);

最后，再次显示处理过的分割图像。