radon变换检测图像中的直线

rgb=imread（‘gantrycrane.png‘）;
 i=rgb2gray（rgb）;
 bw=edge（i‘sobel‘）;
 figure（1）imshow（rgb）;figure（2）imshow（i）;figure（3）imshow（bw）;
 theta_step=1;
 theta=0:theta_step:360;
 [rxp]=radon（bwtheta）;
 figure（4）imagesc（thetaxpr）;colormap（hot）;
 xlabel（‘theta（degrees）‘）;ylabel（‘x\prime‘）;
title（‘r_{theta}（x\prime）‘）;
colorbar;
 max_r=max（max（r））;
threshold=0.75;
[IIJJ]=find（r>=（max_r*threshold））;
[nd]=size（II）;
figure（5