CT系统参数标定及成像.pptx

A题：CT系统参数标定及成像 Parameter calibration and imaging of CT system CONTENT 问题重述 RESTATEMENT 01 CT系统介绍 一种典型的二维CT系统如图1所示，平行入射的X射线垂直于探测器平面，每个探测器单元看成一个接收点，且等距排列。X射线的发射器和探测器相对位置固定不变，整个发射-接收系统绕某固定的旋转中心逆时针旋转180次。对每一个X射线方向，在具有512个等距单元的探测器上测量经位置固定不动的二维待检测介质吸收衰减后的射线能量，并经过增益等处理后得到180组接收信息。 CT系统安装时往往存在误差，从而影响成像质量，因此需要对安装好的CT系统进行参数标定，即借助于已知结构的样品（称为模板）标定CT系统的参数，并据此对未知结构的样品进行成像。 问题重述 RESTATEMENT 01 建立模型，解决以下问题 (1) 在正方形托盘上放置两个均匀固体介质组成的标定模板，模板的几何信息如图2所示，相应的数据文件见附件1，其中每一点的数值反映了该点的吸收强度，这里称为“吸收率”。对应于该模板的接收信息见附件2。请根据这一模板及其接收信息，确定CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离以及该CT系统使用的X射线的180个方向。 (2) 附件3是利用上述CT系统得到的某未知介质的接收信息。利用(1)中得到的标定参数，确定该未知介质在正方形托盘中的位置、几何形状和吸收率等信息。另外，请具体给出图3所给的10个位置处的吸收率，相应的数据文件见附件4。 (3) 附件5是利用上述CT系统得到的另一个未知介质的接收信息。利用(1)中得到的标定参数，给出该未知介质的相关信息。另外，请具体给出图3所给的10个位置处的吸收率。 (4) 分析(1)中参数标定的精度和稳定性。在此基础上自行设计新模板、建立对应的标定模型，以改进标定精度和稳定性，并说明理由。 图2：几何模版 图3：10个位置点 第一问 FIRST QUESTION 02 Mathematica Mathematica是一款科学计算软件，很好地结合了数值和符号计算引擎、图形系统、编程语言、文本系统、和与其他应用程序的高级连接。很多功能在相应领域内处于世界领先地位，它也是使用最广泛的数学软件之一。Mathematica的发布标志着现代科技计算的开始。Mathematica是世界上通用计算系统中最强大的系统。自从1988发布以来，它已经对如何在科技和其它领域运用计算机产生了深刻的影响。 第一问 FIRST QUESTION 02 探测器单元间的距离 直径为8mm的圆遮挡了29个探测单元，于是可以计算出探测单元之间的距离为0.2759mm，我们还根据当探测器处于模版正上方和正左方时，椭圆的长轴和短轴的长度以及其投影在探测器单元上的个数，分别计算得到的距离为0.2759mm和0.2727mm。 旋转中心在正方形托盘中的位置 旋转中心在探测器上的投影位于第256个探测单元上。以探测器在正上方时最左边的点的下方512个单元的地方建立坐标系，记为坐标系1，那么旋转中心的坐标为(256，256)。 第一问 FIRST QUESTION 02 旋转中心在正方形托盘中的位置 射线穿过的长度越长，衰减的能量越多，我们用数据画出了CT探测器接收到的衰减能量的分布图如下。我们假设探测器在正上方时左边为第1个单元，右边为第512个单元。 当探测器位于模版正左方时（如右上图所示），圆中心与椭圆中心的连线和CT探测器垂直，所以检测到的数值最大的传感器对应的横坐标是235，即表示椭圆中心对应检测器单元序号为235。 当探测器处于正下方，即两个物体成像的中心相距最远时（如右下图所示），小圆中心与椭圆中心的连线与CT射线垂直，所以得到了四组椭圆中心对应的探测器单元序数222、223、224、225，取平均值为223.5。 第一问 FIRST QUESTION 02 旋转中心在正方形托盘中的位置 所以椭圆中心在坐标系1中坐标为（288.5,235），经过坐标变换将旋转中心坐标转换到坐标系2内，得出旋转中心坐标为（-32.5，21）。结合上一部分求出的探测器单元之间的距离，可以计算出旋转中心在距正方形托盘左边41.0345，距上边44.2069处。 由图二中给出的各个长度数据我们可以看出旋转中心在椭圆的内部，右图一为任一角度时的旋转中心与椭圆的位置关系，右图二为旋转中心在坐标系2中的位置。 第一问 FIRST QUESTION 02 CT系统使用的X射线的180个方向 建立几何模型，在选定的直角三角形中利用距离计算，运用反三角函数计算出探测器与水平线之间的角度，然后再转换成射线与水平线的夹角，也就是系统中射线的180个方向。要找出两个圆心