CT系统标定和成像.pdfVIP

CT 系统参数标定及成像问题的研究 摘要 本文通过建立数学模型，给出了由大量数据标定CT 系统参数和数据转化为 图像的方法。 针对问题一，主要是对附件一、附件二中数据的分析和处理，标定CT 系统 的各个参数。根据附件一中模板的吸收率，对样品做连续化处理得到标定模板的 方程，建立以椭圆中心为原点的坐标系，引入旋转中心等参数。用Matlab 绘制 出模板接收信息的图像，将其与标定模板示意图进行对比，又因为介质为均匀介 质，分析得到：表格数据大小反映射线扫描到的介质长度，每一列不为零的数据 个数反映了射线照射到有介质区域的单元探测器个数。通过模板在探测器上的投 影占据的探测器单元个数和单个数据对应的探测器单元的相对位置，确定了 CT 系统的探测器单元间距为 0.2767mm,旋转中心坐标为 （-9.2233,6,0182）。每个 方向每个单元探测器接收到信息的数值与表格数据相对应，列出512 个方程，得 到180 个方程组，用最小二乘法算得误差最小时所对应的X 射线旋转的180 个方 向。 针对问题二，根据问题一得到的标定参数，利用上述 CT 系统和某未知介质 的接收信息，模板接收信息矩阵作radon 逆变换、平移、像素转换等操作后，得 到模板原始图像的 “吸收率”为模板接收信息的2.0751 倍，由未知介质的接收 信息、标定参数及真实吸收率与接收信息间的关系，确定该未知介质外部轮廓为 一个椭圆，分布在正方形区域中心位置，其中第三和第四象限分别被挖空了一个 小椭圆，一、二象限间填充了两个小椭圆。 针对问题三，处理方法与问题二类似，由于图像受数字噪声影响，所以先降 噪，再运用滤波反投影法 （FBP），先修正、后投影重建图像的做法，得到原始图 像的吸收率信息。 针对问题四，首先对问题一中射线旋转方向、探测器单元间距和旋转中心坐 标分别进行精度分析和稳定性分析，又根据附件中所提供的数据，分别算出三个 参数的精度，得到精度为0.01 左右。在此基础上建立新的模板，给出相关参量， 最后对新模板的参数进行分析并作出评价。 关键词：CT 系统、最小二乘原理、Radon 逆变换、滤波反投影重建法 一、问题重述 CT (Computed Tomography)可以在不破坏样品的情况下，利用样品对射线能 量的吸收特性对生物组织和工程材料的样品进行断层成像，由此获取样品内部的 结构信息。一种典型的二维CT 系统如图1 所示，平行入射的X 射线垂直于探测 器平面，每个探测器单元看成一个接收点，且等距排列。X 射线的发射器和探测 器相对位置固定不变，整个发射-接收系统绕某固定的旋转中心逆时针旋转 180 次。对每一个X 射线方向，在具有512 个等距单元的探测器上测量经位置固定不 动的二维待检测介质吸收衰减后的射线能量，并经过增益等处理后得到180 组接 收信息。而CT 系统安装时往往存在误差，从而影响成像质量，因此需要对安装 好的CT 系统进行参数标定，即借助于已知结构的样品（称为模板）标定CT 系统 的参数，并据此对未知结构的样品进行成像。 请建立相应的数学模型和算法，解决以下问题： (1) 在正方形托盘上放置两个均匀固体介质组成的标定模板，模板的几何信 息如图2 所示，相应的数据文件见附件1，其中每一点的数值反映了该点的吸收 强度，这里称为“吸收率”。对应于该模板的接收信息见附件2。请根据这一模 板及其接收信息，确定CT 系统旋转中心在正方形托盘中的位置、探测器单元之 间的距离以及该CT 系统使用的X 射线的180 个方向。 (2) 附件3 是利用上述CT 系统得到的某未知介质的接收信息。利用(1)中得 到的标定参数，确定该未知介质在正方形托盘中的位置、几何形状和吸收率等信 息。另外，请具体给出图3 所给的 10 个位置处的吸收率，相应的数据文件见附 件4。 (3) 附件5 是利用上述CT 系统得到的另一个未知介质的接收信息。利用(1) 中得到的标定参数，给出该未知介质的相关信息。另外，请具体给出图3 所给的 10 个位置处的吸收率。 (4) 分析(1)中参数标定的精度和稳定性。在此基础上自行设计新模板、建 立对应的标定模型，以改进标定精度和稳定性，并说明理由。 图1. CT 系统示意图 图2. 模板示意图（单位：mm） 图3. 1