example6_2.m

来自「《MATLAB 6.5辅助图像处理》 本书是基于MATLAB 6.5的图像处」· M 代码 · 共 14 行

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