CT实验一.docxVIP

实验课程：现代医学影像设备 生物医学工程学院实验技术管理中心 实验题目： 1. 傅立叶切片定理证明 实验目的： 加深对傅立叶切片定理的理解以及强化matlab编程能力。 实验原理： 傅立叶切片定理 实验内容： 用matlab证明在x线束平行于x轴时傅立叶切片定理的正确性 实验步骤 1.运行matlab程序，在命令窗口输入以下exit moxing 输入以下代码： I=phantom(512); //生成shepp-logan体模 H=imrotate(I,270); // figure,imshow(I); figure,imshow(H); 显示图像如下： 图1 shepp-ladon体模 512X512 图2体模旋转180度 矩阵512X512 2、X线以与Y轴平行的方向进行投影，根据傅立叶切片定理知道=180度， 1）获取在=180度角度下的体模的投影数据 projection=sum(H,1); figure,plot(1:512,projection); 图如下： 2）再把此角度下的投影数据进行傅立叶变换，图如下 f=fft(projection); %进行傅立叶变换 fab=f.*conj(f); %由于变换是复数，取模 f3=fftshift(fab); %对所得的傅立叶变换进行频谱转移 figure,plot(