计算方框图像在0°到45°方向上的radon变换.m

%计算方框图像在0°到45°方向上的radon变换
i=zeros(100,100);
i(25:75,25:75)=1;%产生一个正方形的黑框
imshow(i);
[r,xp]=radon(i,[0 45]);
figure;
%显示图像在0°方向上的radon变换
plot(xp,r(:,1));
title('r_{0^o}(x\prime)');
figure;
%显示图像在45°方向上的radon变换
plot(xp,r(:,2));
title('r_{45^o}(x\prime)');