restrictions.txt

在matlab环境中利用产生的大脑图实现图像的重建。

P=phantom(256);   %用phantom函数产生Sheep-Logan的大脑图，n为图像p中的行列数，默认为256

imshow(P)

title('原始图像')

%以下为三种不同角度的投影模式

theta1=0:10:170;[R1,xp]=radon(P,theta1);   %存在18个角度投影

theta2=0:5:175;[R2,xp]=radon(P,theta2);    %存在36个角度投影

theta3=0:2:178;[R3,xp]=radon(P,theta3);    %存在90个角度投影

figure,imagesc(theta3,xp,R3);colormap(hot);colorbar; 

%显示图像Sheep-Logan的radon变换

title('经radon变换后的图像')

xlabel('\theta');ylabel('x\prime');     %定义坐标轴

%用三种情况的逆radon变换来重建图像

I1=iradon(R1,10);

I2=iradon(R2,5);

I3=iradon(R3,2);

figure,imshow(I1)

title('角度增值为10时的iradon变换图像')

figure,imshow(I2)

title('角度增值为5时的iradon变换图像')

figure,imshow(I3)

title('角度增值为2时的iradon变换图像')