基于单目标优化模型和图像重建算法的CT系统研究.pdfVIP

基于单目标优化模型和图像重建算法的CT 系统研究 摘要 本文针对CT系统参数标定和未知介质信息的确定问题，基于黄金分割法、 单目标优化模型、搜索算法、图像重建算法，求出了CT系统的旋转中心、探测 器单元距离和射线的180 个方向。根据标定的系统参数，求出了与接收数据对应 的未知介质的位置、几何形状和吸收率。最后设计了新模板，对系统参数的精度 和稳定性进行了改进。 对于探测器单元距离的求解问题，根据几何关系可以求出介质厚度的理论值， 基于此建立单目标优化模型，目标函数为相邻接受信息理论比值与实际比值的最 小误差平方和，遍历射线到介质边缘的距离，决策变量为探测器单元距离。通过 黄金分割算法，逐渐缩小探测器单元距离范围直到满足精度要求，求得探测器单 元距离为0.2768 。 对于CT系统旋转中心位置的确定问题，首先建立以椭圆中心为原点的直角 坐标系，以短轴向右建立轴，长轴向上建立轴。建立单目标优化模型，目标 函数为相邻接受信息理论比值与实际比值的最小误差平方和，决策变量为穿透介 质的第一条射线到介质边缘的距离。通过遍历搜索算法，找到最优距离。然后根 据射线之间的间距建立中心线方程组。分别研究射线平行 ，轴入射时的情况， 得到旋转中心的坐标为(−9.3040,6.2149) 。 对于射线的180 个方向的确定问题，首先通过最小二乘法原理求得接收数据 与厚度之间的比例系数为1.7725。对于一个方向的射线，通过Radon变换求得512 个介质厚度值，通过比例系数得到接收数据的理论值，建立单目标优化模型，目 标函数为接收数据的理论值与实际值的最小误差平方和，决策变量为射线与轴 正方向的夹角，通过遍历搜索算法，求得误差最小的方向夹角。通过迭代法解得 180 个方向与轴正方向的夹角。 对于未知介质的位置、几何形状和吸收率的求解问题，基于中心切片定理， 设计了滤波反投影求解算法。首先根据反投影算法重建图像，然后进行滤波和降 噪两步操作，得到具有理想效果的滤波反投影图像重建算法，带入附件3 和附件 5 的接收数据，重建出未知介质图像。对于给定的 10 个位置，附件 3 的吸收率 为0.0000，1.0004，0.0000，1.1846，1.0488，1.4817，1.2965，0.0000，0.0000， 0.0000，附件5 的吸收率为0.0000，2.8091，7.0353，0.0000，0.0101，3.2351， 6.1054，0.0000，8.1176，0.0000。 我们对CT系统标定的参数分别进行了精度和稳定性分析，设计了由正方形 和圆形均匀固体介质组成的新模板，建立了相应的标定模型，并阐述了新模板的 优点。 关键词 黄金分割法 单目标优化模型 遍历搜索算法 滤波反投影重建 1 一、问题重述 CT可以在不破坏样品的情况下，利用样品对射线能量的吸收特性对生物组 织和工程材料的样品进行断层成像，由此获取样品内部的结构信息。一种典型的 二维CT系统中，平行入射的X 射线垂直于探测器平面，每个探测器单元看成一 个接收点，且等距排列。X 射线的发射器和探测器相对位置固定不变，整个发射 -接收系统绕某固定的旋转中心逆时针旋转 180 次。对每一个X 射线方向，在具 有 512 个等距单元的探测器上测量经位置固定不动的二维待检测介质吸收衰减 后的射线能量，并经过增益等处理后得到180 组接收信息。 CT系统安装时往往存在误差，从而影响成像质量，因此需要对安装好的CT系 统进行参数标定，即借助于模板标定CT系统的参数，并据此对未知结构的样品 进行成像。 建立相应的数学模型和算法，解决以下问题： (1)正方形托盘上放置两个模板，已知每一点的吸收率。根据这一模板及其 接收信息，确定CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距 离以及该CT系统使用的X 射线的180 个方向。 (2)利用(1)中得到的标定参数，根据附件3 和附件5 所给未知介质的接收信 息，分别确定未知介质在正方形托盘中的位置、几何形状和吸收率等信息。另外， 具体求出给定的10 个位置处的吸收率。 (3) 分析(1)中参数标定的精度和稳定性。在此基础上自行设计新模板、建 立对应的标定模型，以改进标定精度和稳定性，并说明理由。