CT系统参数标定和成像.PDFVIP

CT 系统参数标定及成像 摘要 CT 系统（Computed Tomography，电子计算机断层扫描），是利用样品对X 射线、γ射 线、超声波等波的能量吸收特性，再由探测器接收透过样品的波信号，进而分析可得样品的 内部结构信息。本题通过研究典型平行束二维CT 系统的成像原理，通过数学建模的思想方 法，解决相应问题。 对于问题（1），首先，通过确定椭圆长短轴、圆直径的投影宽度对应的接收器单元数量， 根据平行投影等宽的原理，可以确定接收器单元间距。其次，椭圆在某个方向上的投影长度 与该方向的角度存在着函数关系，因此由每条数据求出其相应投影长度进而求出其对应角度。 最后，旋转中心经过探测器平面中垂线，则通过计算椭圆中心在长短轴正投影是与探测器的 中垂线的偏移量确定旋转中心。但该方法所得的旋转中心精度为 0.2mm 左右，接下来通过 遗传算法对旋转中心坐标的求解进行优化。 对于问题 （2 ）、（3 ），在问题（1）中得到标定参数 （旋转中心、探测器单元间距、射线 方向）后，根据Radon 变换及反演原理，编写本题参数下反演的Matlab 程序。再分别对附 件3、附件5 提供的未知介质的吸收信息，反演出问题（2 ）、（3 ）中未知介质的吸收率矩阵， 通过矩阵缩放变换得到题目要求的256*256 的吸收率矩阵。再根据吸收率矩阵，进行成像， 即可得到未知介质的位置、几何形状，而吸收率信息可直接由矩阵反映。 对于问题（4 ），分析了影响参数标定的精度的因素。根据误差传播定律计算了模型误差， 并在分析不同情况下数据变化下的模型稳定性，最后设计了新的标定模板及方法改进精度和 稳定性。 关键词：CT 成像、旋转中心、Radon 变换与反演、遗传算法、误差传播定律 1 一、问题重述 如图1 所示典型的二维CT 系统，平行入射的X 射线垂直于探测器平面，每个探测器单 元看成一个接收点，且等距排列。X 射线的发射器和探测器相对位置固定不变，整个发射- 接收系统绕某固定的旋转中心逆时针旋转 180 次。对每一个X 射线方向，在具有512 个等 距单元的探测器上测量经位置固定不动的二维待检测介质吸收衰减后的射线能量，并经过增 益等处理后得到180 组接收信息。解决以下问题： (1) 在正方形托盘上放置两个均匀固体介质组成的标定模板，请根据这模板及其接收信 息，确定CT 系统旋转中心在正方形托盘中的位置、探测器单元之间的距离以及该CT 系统 使用的X 射线的180 个方向。 (2)利用CT 系统得到的某未知介质的接收信息。利用(1) 中得到的标定参数，确定该未知 介质在正方形托盘中的位置、几何形状和吸收率等信息，并具体给出图2 所给的 10 个位置 处的吸收率。 (3)利用CT 系统得到的另一个未知介质的接收信息。利用(1) 中得到的标定参数，给出该 未知介质的相关信息，并具体给出图2 所给的10 个位置处的吸收率。 (4) 分析(1)中参数标定的精度和稳定性。在此基础上自行设计新模板、建立对应的标定 模型，以改进标定精度和稳定性，并说明理由。 图1. CT 系统示意图 图2. 10 个位置示意图 二、模型假设 (1). 假设X 射线波长一定； (2). 假设探测器的效率为100％； (3). 忽略探测器所探测到的一次或多次散射光子。此外，当射线能量高于1.02 MeV 时，忽 略探测器所探测到的( 由电子对效应产生的)正电子与电子湮没而产生的X 射线； (4). 忽略X 射线源焦点的物理尺度； (5). 忽略探测器单元的物理尺寸。 2 三、符号说明 符号 意义 旋转中心横坐标(mm)