MATLAB图像处理之图像重构

P = phantom(256);%生成头骨幻影图像
imshow§%显示头骨幻影图像
theta1 = 0:10:170;%投影的角度，步长为10
[R1,xp] = radon(P,theta1);%radon变换
num_angles_R1 = size(R1,2);%角度的个数
theta2 = 0:5:175;%投影角度，步长为5
[R2,xp] = radon(P,theta2);%radon变换
num_angles_R2 = size(R2,2);%角度的个数
theta3 = 0:2:178;%投影角度，步长为2
[R3,xp] = radon(P,theta3);%radon变换
num_angles_R3 = size(R3,2);%角度的个数
N_R1 = size(R1,1);%角度步长为10时对角线的长度
N_R2 = size(R2,1);%角度步长为5时对角线的长度
N_R3 = size(R3,1);%角度步长为2时对角线的长度
figure,
imagesc(theta3,xp,R3);%显示角度步长为2时的Radon变换
colormap(hot); colorbar
xlabel(‘旋转角度’); ylabel(‘感知器位置’);
output_size = max(size§);%确定变换后图像的大小
dtheta1 = theta1(2) - theta1(1);%步长
I1 = iradon(R1,dtheta1,output_size);%radon逆变换
figure,
subplot(1,3,1), imshow(I1)%显示逆变换的图像
dtheta2 = theta2(2) - theta2(1);%步长
I2 = iradon(R2,dtheta2,output_size);%radon逆变换
subplot(1,3,2), imshow(I2)%显示逆变换的图像
dtheta3 = theta3(2) - theta3(1);%步长
I3 = iradon(R3,dtheta3,output_size);%radon逆变换
subplot(1,3,3), imshow(I3)%显示逆变换的图像
D = 250; dsensor1 = 2;
F1 = fanbeam(P,D,‘FanSensorSpacing’,dsensor1);%fanbeam变换
dsensor2 = 1;
F2 = fanbeam(P,D,‘FanSensorSpacing’,dsensor2);%fanbeam变换
dsensor3 = 0.25;
[F3, sensor_pos3, fan_rot_angles3] = fanbeam(P,D,…
‘FanSensorSpacing’,dsensor3);%fanbeam变换
figure,
imagesc(fan_rot_angles3, sensor_pos3, F3)%显示fanbeam变换的图像
colormap(hot); colorbar
xlabel(‘旋转角度’); ylabel(‘感知器位置’)
Ifan1 = ifanbeam(F1,D,‘FanSensorSpacing’,…
dsensor1,‘OutputSize’,output_size);%fambeam逆变换
figure,
subplot(1,3,1), imshow(Ifan1)%显示重构图像
Ifan2 = ifanbeam(F2,D,‘FanSensorSpacing’,…
dsensor2,‘OutputSize’,output_size);%fanbeam逆变换
subplot(1,3,2), imshow(Ifan2)%显示重构图像
Ifan3 = ifanbeam(F3,D,‘FanSensorSpacing’,…
dsensor3,‘OutputSize’,output_size);%fanbeam逆变换
subplot(1,3,3), imshow(Ifan3)%显示重构图像