程序后面两句有问题,可这样改:
%显示坐标刻度
iptsetpref('ImshowAxesVisible','on');
%创建图像并显示
I=zeros(100,100);
I(25:75,25:75)=1;
figure(1);
imshow(I);title('原始图像');
%定义方向并计算Radon变换
theta=0:180;
[R,xp]=radon(I,theta);
%结果显示
figure(2);
%imshow(theta,xp,R,[],'InitialMagnification','fit');
imshow(R,[],'Xdata',theta,'Ydata',xp,'InitialMagnification','fit')
title('Radon变换结果');
xlabel('\theta(degree)');
ylabel('x''');
colormap(gca,hot),colorbar