帮我写一段边缘检测算法求边缘点集的MATLAB代码

边缘检测是图像处理中的基础步骤，MATLAB提供了一些内置函数来实现这一过程，比如Sobel算子、Canny边缘检测等。下面是一个简单的使用Sobel算子进行边缘检测的例子：

% 加载图片
img = imread('image.jpg'); % 替换为你的图像文件名

% 转化为灰度图像
gray_img = rgb2gray(img);

% 应用Sobel算子进行水平和垂直方向的梯度计算
sobelx = imfilter(gray_img, fspecial('sobel', 'x'));
sobely = imfilter(gray_img, fspecial('sobel', 'y'));

% 计算梯度幅值和角度
grad_mag = sqrt(sobelx.^2 + sobely.^2);
grad_angle = atan2(sobely, sobelx);

% 设置阈值进行非极大值抑制
[edges, ~] = edge(grad_mag, 'canny', [50 150]); % 可调整阈值

% 找到边缘点集
edge_points = bwmorph(edges, 'thin', Inf); % 使用连通组件分析去除噪声并获取点集

figure;
subplot(2, 1, 1), imshow(gray_img), title('Original Image');
subplot(2, 1, 2), imshow(edges), title('Detected Edges');

注意：这个代码片段假设你已经安装了MATLAB，并且你的系统中包含了`imread`、`rgb2gray`、`imfilter`、`fspecial`、`edge`和`bwmorph`等相关函数。上述例子使用了Canny边缘检测，你可以根据需要选择其他边缘检测算法。