example6_2.m

《matlab6.5辅助图像处理》源代码

%首先，提取图像的局部边缘，见图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