CT系统标定及成像-(1).docxVIP

CT系统参数标定及成像研究 摘要 本文通过建立数学模型，给出了通过分析大量数据标定CT参数的方法，得到了通过数据借助模型和相关的算法及函数生成图像的方法；以及其逆向运算方法。 针对问题一，主要是对附件一、附件二中数据的分析和处理，标定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)-(4)中的所有数值结果均保留4位小数。同时提供(2)和(3)重建得到的介质吸收率的数据文件（大小为256×256，格式同附件1，文件名分别为problem2.xls和problem3.xls） 图1. CT系统示意图 图2. 模板示意图（单位：mm） 图3. 10个位置示意图 二、问题分析 1、问题一 该问题主要为数据分析处理类问题。首先对附件一进行分析，通过对模板每一点吸收率连续化处理分别得到椭圆和圆的方程。 分析附件二中数据，并用MATLAB画出对应矩阵的图像，每一列不为0的数据个数为Mi，反映第i次扫描有介质区域的长度，对比模板，反映探测器与模板的相对位置，Mi最大时探测器平行于y轴，Mi最时探测器平行于x轴；模板由均匀固体介质组成，则经过增益后的数据与实际数据成比例，附件中每个单元格中数据即为增益后的截模板弦长 选取表格中数据量较大处的数据与相应的实际数据用近似最小二乘定理算得误差最小时的比例系数μ 三、模型假设 1、假设CT系统的旋转中心位于探测器中垂线上 2、假设探测器上每一点发出光强相同，且均匀分布在整个探测器上 3、忽略光的衍射 四、符号系统 Mi：第i列不为0的数据个数(i=0,1,2,… x轴：与椭圆短半轴重合的数轴 y轴：过椭圆中心且与长半轴重合的数轴 x0 y0 d：探测器单元之间的距离 Bij：第i次旋转第j个 Pij：表格中第i列第j行的数据(i=1,2…，180，j=1,2,… Lij：探测器第i次扫描时第j个单元探测器扫描到的介质弦长和(i=1,2…，180，j=1,2,… μ：问题一中弦长与衰减后射线能量的比例系数 ki：射线斜率(i=1,2,… θi：旋转第i次与X正半轴的夹角(i=0,1,2,… nj：第n个探测器单元（j=1,2,3，… 五、模型建立 1、问题一 该问题主要为数据分析处理类问题。首先对附件一进行分析，通过对模板每一点吸收率连续化处理分别得到椭圆和圆的方程与理论一致分别为： x 分析附件二中数据，并用MATLAB画出对应矩阵的图像（如图一），根据图像，探测器逆时针旋转，则探测器处于平行于y轴时，沿x负半轴方向分别为第1个到第512个探测单元。 Mi反映第i次扫描有介质区域的长度，反映探测器与模板的相对位置，对比模板实际位置，Mi最大时探测器平行于y轴，Mi最时探测器平行于x