【数学建模大赛】ct系统参数的标定与成像问题.docxVIP

CT系统参数标定及成像 芋 本文主要研究CT系统参数的标定与成像问题，通过建立模型利用结构 样品标定CT系统的参数，并对未知结构样品的图像进行重建。 针对问题一，先对圆分析，无论射线的方向如何，所得到的投影长度总是 等于圆的直径，由此关系可求得相邻探测器单元的距离为。再利用几 何知识求出圆和椭圆在某个投影方向上的投影长度与该方向角度的关系式，并将 其有较大跳跃或异常的点去除，利用傅里叶变换拟合出一个函数，再将附件2 中180个角度对应的数据与拟合函数在同一角度得出的数据做差，用最小二乘法 拟合出一个更加精确的函数，再将附件2中的数据代入进去求出180个角度，其 中前三个以y轴正方向的夹角分别是29.0661°, 30.5702°, 31.5755°。以托盘左 下角为原点，水平方向为横坐标，竖直方向为纵坐标建立直角坐标系，同样利用 最小二乘法得出旋转中心为(40.7543,56.2520)。 针对问题二，首先找到投影坐标与原图像坐标之间的转化关系，再利用滤波 反投影法得出图像的衰减系数分布函数，即所有经过滤波后的射线投影在 0-180度范围的衰减系数的叠加。考虑到CT系统是逆时针旋转180次，每一 次对应一个角度，x射线对物体进行平行束扫描，得到是180个离散投影数据， 于是将衰减系数分布函数离散化，将180个角度得到的投影数据叠加，得到图像 的像素点，从而重建出图像，并算出了附件四给出的10个位置的吸收率分别为: 、、、、、、、、、 0. 0307o 针对问题三，利用问题二中的模型重建出图像并解出10个位置的吸收率分 别为：、、、、、、、、 、0. 0774o 针对问题四，关于参数标定的精度与稳定性问题，利用CT系统的结构以及 成像原理可知，x射线对物体进行平行束扫描获得了 180个离散数据，由于CT系 统的探测器之间具有一定的间隔，所以物体的边界可能有一局部无法被收集到信 息，从而造成投影长度缩小的情况，导致所求探测器的单元间隔的增大，并且造 成旋转中心与x射线旋转角的误差。通过分析椭圆投影到探测器上长度所对应的 误差分析，得到投影长度越长，模板以下图形面积越小，所对应的误差就越小，即 可得出一个精度较高的新模板。 关键词：最小二乘法；滤波反射法；衰变系数；傅里叶变换; 一、问题的重述 问题的背景： 计算机层析成像（Computed Tomography, CT）是通过对物体进行不同角度的 射线投影测量而获取物体横截面信息的成像技术，涉及多个学科领域。CT技术 不但给诊断医学带来革命性的影响.还成功地应用于无损检测、产品反求和材料 组织分析等工业领域，其实质是由扫描所得到的投影数据反求出成像平面上每个 点的衰减系数值。有一种典型的二维CT系统，平行入射的X射线垂直于探测器 平面，每个探测器单元看成一个接收点，且等距排列。X射线的发射器和探测器 相对位置固定不变，整个发射-接收系统绕某固定的旋转中心逆时针旋转180次。 对每一个X射线方向，在具有512个等距单元的探测器上测量经位置固定不动的 二维待检测介质吸收衰减后的射线能量，并经过增益等处理后得到180组接收信 ，息、O 问题的目的： CT系统安装时往往存在误差，从而影响成像质量，因此需要对安装好的CT 系统进行参数标定，即借助于结构的样品（称为模板）标定CT系统的参数, 并据此对未知结构的样品进行成像。 所要解决的问题： 问题一：在正方形托盘上放置两个均匀固体介质组成的标定模板，模板的 几何信息已经给定，相应的数据文件见附件1,其中每一点的数值反映了该点的 吸收强度，这里称为“吸收率”。对应于该模板的接收信息见附件2。请根据这 一模板及其接收信息，确定CT系统旋转中心在正方形托盘中的位置、探测器单 元之间的距离以及该CT系统使用的X射线的180个方向。 问题二：附件3是利用上述CT系统得到的某未知介质的接收信息。利用⑴ 中得到的标定参数，确定该未知介质在正方形托盘中的位置、几何形状和吸收率 等信息。另外，请具体给出所给的10个位置处的吸收率，相应的数据文件见附 件4。 问题三：附件5是利用上述CT系统得到的另一个未知介质的接收信息。利 用（1）中得到的标定参数，给出该未知介质的相关信息。另外，请具体给出所给 出的10个位置处的吸收率。 问题四：分析（1）中参数标定的精度和稳定性。在此根底上自行设计新模板、 建立对应的标定模型，以改良标定精度和稳定性，并说明理由。 二、问题的分析 本文主要研究CT系统两方面的问题，一是利用结构的样品和它的投影 信息得到CT系统的参数，二是用己知的标定CT系统的参数来对未知结构的样品 进行成像。 2. 1问题一的分析： 该问题要求利用附件1和附件2的信息来确定出CT系统的旋转中心、相邻 探测单元之间的距离以及X射线每次旋转的方向。针对这