CT实验平行束投影与反投影重建报告.docxVIP

专业课程实验报告------“CT成像及原理”姓名：郭炳映指导老师：傅健学 号： 院 系： 机械工程及自动化 2012年11月15日专业课程实验任务书CT成像原理及算法一、教学目标结合计算机层析成像原理（CT）与数字图像软件平台，让学生在层析成像理论和软件调试等方面得到综合训练，提高学生对所学知识的综合应用能力，为毕业设计及今后从事相关工作培养分析和解决工程实际问题的能力。二、设备及原理1) 已有的数字图像平台2) 二维平行束计算机层析成像原理3) 二维扇束计算机层析成像原理三、内容与步骤1．学习二维平行束计算机层析成像原理（2周）2．编制二维平行束投影，滤波与重建程序（2周）3．将二维平行束CT程序融于数字图像平台软件（1周）4．学习二维扇束计算机层析成像原理（2周）5．编制二维扇束投影，滤波与重建程序（1周）6．将二维扇束CT程序融于数字图像平台软件。（1周）7．写出课程设计报告（1周）四、参考书目：1.《CT原理与算法》庄天戈编著上海交通大学出版社2.《由投影重建图像》赫尔曼编著科学出版社3.《计算机断层成像技术》 Jiang Hsieh 编著科学出版社4．《C语言程序设计》谭浩强编著清华大学出版社真实投影数据的模拟投影的原理：实际的投影数据及真实的投影数据是由检测器得到的，在计算机模拟时，则取自仿真模型。仿真模型常由一系列大小不同、位置各异，具有给定密度的椭圆组成，以模拟脑部不同的组织（图A.1）。每一椭圆的位置，用其中心坐标（G、H），长轴与y轴的交角α来表示（图A.3）。其大小自然由长轴2a与短轴2b确定。先讨椭圆中心在坐标原点，长轴与y轴重叠时，穿过椭圆的某一射线投影的求法。投影的计算机算法实现：for(i=0;i180;i++) //角度，从0到180 {for(d=-128, j=0;d128;j++,d++) //距离，从-128到128 { r=p.a*p.a*cos(pi*i/180-pi*p.angle/180)*cos(pi*i/180-pi*p.angle/180)+p.b*p.b*sin(pi*i/180-pi*p.angle/180)*sin(pi*i/180-pi*p.angle/180); d1=d-p.x*cos(pi*i/180)-p.y*sin(pi*i/180);if((r-d1*d1) 0) pa[i*256+j]=2*p.a*p.b*sqrt(r-d1*d1)/r; //代入公式求出投影数据else pa[i*256+j] = 0; } }得到的投影效果图：以长轴100，短轴50，中心坐标（0，0），偏转角45度为例：滤波投影2.1 滤波的原理及滤波函数的选取：2.2 滤波投影的计算机实现：for(i=0;i180;i++)//补全投影数据{for(j=0;j256;j++)q[j]=(pa[i*256]+pa[i*256+1])/2; //前255个点取投影0和1的平均值for(j=256;j511;j++)q[j]=pa[i*256+j-256]; //中间255个点不变for(j=511;j766;j++)q[j]=(pa[i*256+254]+pa[i*256+255])/2; //后255个点取投影254和255的平均值for(int l=0;l256;l++) //用S-L滤波函数滤波{sum=0;for(int k=0;k511;k++)sum+=q[k+l]*(-2)/(pi*pi*(4*(k-255)*(k-255)-1));pa[i*256+l]=sum;}}2.3 滤波投影效果图：以长轴100，短轴50，中心坐标（0，0），偏转角45度为例：反投影重建3.1 反投影重建基础——中心切片定理：3.2 卷积反投影算法：3.3 反投影的计算机实现：for(x=0;x256;x++)for(y=0;y256;y++){sum1=0;for(i=0;i180;i++){int n;float m;m=128*(1-cos(pi*i/180)-sin(pi*i/180))+(x-1)*cos(pi*i/180)+(y-1)*sin(pi*i/180);m=m+1;if(m0||m255)continue;else {n=floor(m);m=m-n;sum1+=((1-m)*pa[i*256+n]+m*pa[i*256+n+1]); }}a[x*256+y]=sum1 *pi/180;}3.4 反投影效果图：对本课程的感想CT的投影与重建是在完成数字图像处理平台后进行的，在经过了对位图