帮我写一份MATLAB的利用radon变换检测直线代码

好的，以下是MATLAB利用radon变换检测直线的代码：

1、读取图像并转化为灰度图像

img = imread('image.png'); % 读取图像
gray_img = rgb2gray(img); % 转换为灰度图像

2、应用radon变换

theta = 0:179; % 设定变换角度
[R,xp] = radon(gray_img, theta); % radon变换

3、寻找最大峰值

[maxValue,columnIndex] = max(R(:)); % 找到最大峰值及其所在位置
[rowIndex, columnIndex] = ind2sub(size(R),columnIndex); % 将位置转换为行列坐标

4、绘制检测结果

figure;
imshow(gray_img,[],'InitialMagnification','fit');
hold on;
x = 0:size(gray_img,2);
y = (cosd(theta(rowIndex))/-sind(theta(rowIndex)))*x + (xp(columnIndex)/sind(theta(rowIndex)));
line(y+1,x+1,'Color','r','LineWidth',2);
title(sprintf('Line Detection Using Radon Transform\nAngle = %0.1f degrees',theta(rowIndex)))
axis normal
hold off;

希望这份代码对您有所帮助！