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,dt