CT扫描成像论文.docx

PAGE 11 CT系统参数标定及成像 对于问题二，首先根据问题一中的标定参数，将附件3中的数据变换为旋转中心在正方形托盘正中心的数据。然后将数据进行反拉东变换得到反投影图像，发现其存在较多伪影，成像效果不佳。因此结合伪影形成的原因，考虑利用傅里叶中心切片定理，设计滤波反投影算法（FBP），先将投影数据进行傅立叶变换，滤波后逆傅里叶变换，将所得的值在反投影平面累加，实现吸收率图像重构。最终得到未知介质的位置、形状以及要求十点的吸收率。 对于问题三，题目要求与问题二基本一致，但附件五数据较附件三复杂的多，干扰因素加多，需对滤波反投影重建模型优化。该干扰因素主要为噪音影响，所以在问题二模型的基础上进行降噪处理，去噪后得到重建图像以及对应十点的吸收率。  一、问题的提出 1.1问题的背景 CT(Computed Tomography)可以在不破坏样品的情况下，利用样品对射线能量的吸收特性对生物组织和工程材料的样品进行断层成像，由此获取样品内部的结构信息。一种典型的二维CT系统如图1所示，平行入射的X射线垂直于探测器平面，每个探测器单元看成一个接收点，且等距排列。X射线的发射器和探测器相对位置固定不变，整个发射-接收系统绕某固定的旋转中心逆时针旋转180次。对每一个X射线方向，在具有512个等距单元的探测器上测量经位置固定不动的二维待检测介质吸收衰减后的射线能量，并经过增益等处理后得到180组接收信息。 CT系统安装时往往存在误差，从而影响成像质量，因此需要对安装好的CT系统进行参数标定，即借助于已知结构的样品（称为模板）标定CT系统的参数，并据此对未知结构的样品进行成像。 1.2目标任务 根据以上背景，以及给出的5个附件，需要解决以下四个问题： (1) 在正方形托盘上放置两个均匀固体介质组成的标定模板，模板的几何信息如图2所示，相应的数据文件见附件1，其中每一点的数值反映了该点的吸收强度，这里称为“吸收率”。对应于该模板的接收信息见附件2。请根据这一模板及其接收信息，确定CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离以及该CT系统使用的X射线的180个方向。 (2) 附件3是利用上述CT系统得到的某未知介质的接收信息。利用(1)中得到的标定参数，确定该未知介质在正方形托盘中的位置、几何形状和吸收率等信息。另外，请具体给出图3所给的10个位置处的吸收率，相应的数据文件见附件4。 (3) 附件5是利用上述CT系统得到的另一个未知介质的接收信息。利用(1)中得到的标定参数，给出该未知介质的相关信息。另外，请具体给出图3所给的10个位置处的吸收率。 (4) 分析(1)中参数标定的精度和稳定性。在此基础上自行设计新模板、建立对应的标定模型，以改进标定精度和稳定性，并说明理由。 (1)-(4)中的所有数值结果均保留4位小数。同时提供(2)和(3)重建得到的介质吸收率的数据文件（大小为256×256，格式同附件1，文件名分别为problem2.xls和problem3.xls） 二、问题的分析 2.1问题一的分析 对于问题一，要求我们根据所给均匀固体介质组成的标定模板的几何信息及其接收信息，确定CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离以及该CT系统使用的X射线的180个方向。本题的难点在于寻找特殊条件确定相关变量。 首先根据图像的几何信息建立合适的直角坐标系，确定已知变量的相关位置；其次对附件2中的数据利用Matlab进行图像还原,通过观察分析图像的特点，并结合原始数据，确定入射方向的近似范围；然后寻找特殊情况，通过对特殊情况的特征分析确定相关数据，经分析，以垂直于椭圆的长轴和短轴的方向为两个特殊情况；最后，根据所得数据求解出三个未知变量。 2.2问题二的分析 题目要求在问题一的标定参数下，确定未知介质的位置、形状及吸收率。通过分析附件三信息可知，所给信息为X射线在该照射方向吸收强度的投影。由此确定我们需要由投影数据来对物体吸收率图像还原，所以需要建立反投影重建模型来还原相关数据。而通过拉东变换可以实现这一点，但重建图像效果有待验证，若效果不理想，需要对重建过程进行改进。 2.3问题三的分析 对于问题三，同样要求确定未知介质的位置、形状等信息，题目要求大致与问题二一致，但由于附件五所给数据比附件三复杂得多，而且分为更多的区域，所以应该需要在问题二模型的基础上对模型进行优化，然后才能更好的运算出结果。 2.4问题四的分析 三、模型的假设 1、假设其中两次入射方向垂直与椭圆的长轴和短轴，无偏差； 2、假设探测器接收的从左向右依次排列数据对应于附件中每列从上至下的数据； 3、假设该CT系统发射平行光，每次的旋转角度一定； 四、符号说明 符号 解释说明 CT系统的旋转中心 探测器单元之间的距离 五、模型的建立与求