代码16-10.txt

matlab宝典书籍源码 全部代码 希望大家喜欢

      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的逆变换结果