CT系统参数标定及反投影重建成像 (1).docVIP

1

CT系统参数标定及反投影重建成像

摘要

CT技术在当代社会已广泛应用在临床医学、工业工程等领域。在本题中首先通过建立离散模型并对其简化从而对CT仪器参数进行了标定；建立平行束滤波反投影重建模型，通过Radon变换及R-L滤波器解决了成像问题，并通过对平行束滤波反投影重建模型进行优化，使模型具有降噪能力，从而得到了更精确符合实际的图形与吸收率值。

对问题一，由于附件2中的数据存在系统增益，通过对附件2数据作消除增益(增益系数为1.776974)处理后，建立离散化模型，由于方程组模型难于求解，故建立一种简化离散模型，且通过MATLAB、Mathmatica等软件对方程进行求解，解得180个方向的精确角度(结果见表5.1)。通过Excel对附件2数据使用条件约束进行色阶处理，可以一定程度上将小圆与椭圆的扫描情况清晰地显示。我们利用附件2中的数据与模板存在的几何关系，求得探测器单元间距范围为(0.2759,0.2796),以椭圆几何中心为原点建立的坐标系xoy中，解得旋转中心坐标为(-9.2996mm,5.5520mm)。

对问题二，由于传统反投影重建算法会引入星状伪迹，我们决定使用基于吸收率的Radon变化及R-L滤波器构成的平行束滤波反投影重建模型。由于系统旋转中心和几何中心不重合，使用Radon变换时会造成数据缺失，故我们对正方形托盘进行了“镶边”处理。在对数据除去增益(此处增益系数为2.0033)处理并且在滤波结束后对“镶边”进行去除后，我们对被检测物体实现了重建并且得到题目所需十个点的吸收率分别为：0.000、0.9722、0.0036、1.1761、1.0426、1.4652、1.2849、0.0007、0.0000、0.0175,其

对应的坐标见表5.2。

对问题三，我们对平行束滤波反投影重建模型进行了优化，加入了降噪处理函数。使用问题二中模型对附件5进行重建后发现所得重建图像存在边界模糊不清等问题。通过对模板进行数据反演后与标准值对比，发现吸收率数据会受非系统因素影响，存在一定波动。所以本文中通过多项对比选取自适应滤波函数这一最优降噪方法对数据进行降噪处理后，减小了波动范围，并拟合出了边界更为清晰的图形(见图5.17),实现了重建并且得到题目所需十个点的吸收率分别为：0.0126、2.2902、5.9159、0.0163、0.0823、3.1336、6.0333、0.0000、7.7184、0.0861,其对应的坐标见表5.3。。

对问题四，首先我们基于标定CT的各项参数，建立归一化均方差评价模型从而对精度和稳定性进行分析。此外基于原模板存在的精度与稳定性的问题，我们建立了新的模板(见图5.20、5.21),新模板具有更高的精度与稳定性，能够对该CT进行更精确地

参数标定。

本题中我们所建立的一系列模型均能够满足题目要求，且层层递进、环环相扣，具有较高的精度与稳定性。其中本文建立的平行束滤波反投影重建模型还可以由二维推广至三维，可以用于太空垃圾的形状确定与分类，航空器在太空中对未知天体的形状确定等领域。

关键词：离散模型、二维平行波反投影重建、Radon变化、R-L滤波器、降噪处理

2

一、问题重述

CT可以在不破坏样品的情况下，利用样品对射线能量的吸收特性对生物组织和工程材料的样品进行断层成像，由此获取样品内部的结构信息。本题X射线的发射器和探测器的相对位置固定不变，整个发射-接收系统绕位于正方形托盘下方某处旋转中心逆时针旋转180次。对每一个X射线方向，发射接收装置装有512个等距单元探测器，用于测量位置固定不动的二维待测介质吸收衰减后的射线能量，并且通过增益等处理方式得到180组接收信息。然而由于存在系统误差，所以需要对安装好的CT系统进行参数标定，通过已知模板对CT系统的参数进行标定，并根据标定的参数对未知结构的样品进行成像。

具体问题重述为下：

(1)在正方形托盘上放置两个均匀固体介质组成标定模板，模板的几何信息如图2给出，相应的数据文件见附件1,其中每一点的数值反映了该点的吸收强度“吸收率”。对应用于模板的接收信息见附件2。问题一要求根据模板及其接收信息对CT系统进行参数标定，确定出此CT系统实际的旋转中心在正方形托盘中的位置、探测器单元间距以及该CT系统使用的X射线的180个方向。

(2)利用问题一所标定的CT系统相关参数以及所建立模型，使用附件3所给未知介质的接收信息，确定出该未知介质在正方形托盘中的位置、几何形状以及吸收率等信息。另外利用附件4中数据给出图3中所给10个位置处的吸收率。

(3)附件5为利用该CT系统得到了另一未知介质的