%首先，提取圖像的局部邊緣，見圖6-24。
I = imread('ic.tif');
BW = edge(I);
imshow(I)
figure, imshow(BW)           
%然后，對邊緣圖像做Radon變換，見圖6-25。
theta = 0:179;
[R,xp] = radon(BW,theta);
imagesc(theta, xp, R); 
colormap(hot);
xlabel('\theta (degrees)'); 
ylabel('x');
colorbar