计算方框图像从0°~180°每隔1°增加多个radon变换.m

i=zeros(100,100);
i(25:85,25:85)=1;
theta=0:180;
[r,xp]=radon(i,theta);
imagesc(theta,xp,r);
title('r_{\theta}(x\prime)');
xlabel('\theta(degrees)');
ylabel('x\prime');
set(gca,'xtick',0:20:180);
colormap(hot);
colorbar