matlab基于分水岭算法的matlab程序

以下是基于分水岭算法的MATLAB程序示例：

clc;clear all;close all;
I=imread('rice.png'); %读取图像
subplot(2,2,1);imshow(I);title('原图');
I=im2double(I);
%计算梯度幅值
Gmag=imgradient(I,'sobel');
subplot(2,2,2);imshow(Gmag,[]);title('梯度幅值');
%计算梯度的方向
Gdir=imgradient(I,'sobel');
subplot(2,2,3);imshow(Gdir,[]);title('梯度方向');
%计算分水岭
L=watershed(-Gmag);
subplot(2,2,4);imshow(label2rgb(L));title('分水岭');

该程序读取一个Rice.png图像，使用Sobel算子计算梯度幅值和方向，然后通过计算分水岭来分割图像。程序的输出结果包括原始图像、梯度幅值、梯度方向和分水岭图像。