www.gusucode.com > 《matlab图像处理与界面编程宝典》秦襄培 编著,每章的MATLAB源代码程序 > 第16章/代码16-10.txt
P = phantom(256); % 生成头部图像 imshow(P) % 显示图像 theta1 = 0:10:170; [R1,xp] = radon(P,theta1); %计算R1的radon变换 theta2 = 0:5:175; [R2,xp] = radon(P,theta2); %计算R2的radon变换 theta3 = 0:2:178; [R3,xp] = radon(P,theta3); %计算R3的radon变换 figure, imagesc(theta3,xp,R3); % 建立新窗口显示图像 colormap(hot); % 设置色彩索引图 colorbar % 显示色彩索引条 xlabel('\theta'); % x轴标签 ylabel('x\prime'); % y轴标签 I1 = iradon(R1,10); % 计算R1的逆变换 I2 = iradon(R2,5); % 计算R2的逆变换 I3 = iradon(R3,2); % 计算R3的逆变换 imshow(I1) % 显示R1的逆变换结果 figure, imshow(I2) % 显示R2的逆变换结果 figure, imshow(I3) % 显示R3的逆变换结果