shepp-logan大脑图是由图像处理工具箱函数phantom创建的。 首先产生一个256灰度级的大脑图: p=phantom(256); imshow(p); 另外投影和重建的程序我也替你编了吧! randon变换: theta1=0:10:170,theta2=0:5:175;theta3=0:2:178; =radon(P,theta1); =radon(P,theta2); =radon(P,theta3); figure,imagesc(theta3,xp,R3); colormap(hot);colorbar xlabel(' \theta ' );ylabel(' x\prime ' ); %利用不同部分的randon逆变换来重构图像 I1=iradon(R1,10);I2=iradon(R2,5);I3=iradon(R3,2); subplot(1,3,1); imshow(I1);title('用R1重构图像'); subplot(1,3,2); imshow(I2);title('用R2重构图像'); subplot(1,3,3); imshow(I3);title('用R3重构图像'); 你的新问题,关于三维图形如何旋转。我的方案是:
使用rotate3D函数命令。例: surf(peaks(20)); rotate3D 在图形窗口的图形区域内按住鼠标左键不放,来调节视角。可以从不同的角度来观察所绘得的图形。 希望对你有用,有问题继续交流。 |
|