2017数学建模A题二等奖论文.pdfVIP

CT 系统参数标定及成像 摘要 目前，如何通过CT 机的探测器已知的投影信息还原出物体几何形状以及内 部结构，探测已知模块得到的数据推导CT 机的参数以及校正，是研究精密仪器 的重要方面。本文通过建立数学模型，给出了通过分析表格数据来解决上述所有 问题的方法。 对于问题一，以正方形左下角为原点建立平面直角坐标系对方向角进行几何 抽象，通过介质上任意一点在转动过程中在光屏上的投影点的运动轨迹都是在旋 转中心的投影点的上下波动这一几何特性，引入三角函数、直流分量和交流分量 并与几何方法结合，建立了已知模板接收信息、吸收率、CT 参数、余弦轨迹的 傅里叶展开之间关系的综合模型。 y a  a cosn w b sinn w ， n 0  1 1  yk tangxk (b a tan) 求解模型时，将附件的探测器接收数据和带入，利用光路方向角与椭圆中心 和圆心之间的距离的几何关系求出180 个方向角度，旋转中心坐标 （50.2775,56.1050）,相邻探测器间距离为0.2775mm。 对于问题二和问题三，本文根据中心切片定理，建立起投影物体的密度与其 二维傅里叶变换后的函数的关系模型。在每一个确定角度下，分别使用直接反投 影法和滤波反投影法 f (x ,y ) | f (R cos,R sin) | G  (R ) G () F(,) F [ f ( x, y) | ]   将变换后的函数作傅立叶逆变换就可以得到原物体的密度函数，并进行插值 填充由MATLAB 生成吸收率矩阵。问题二的介质几何中心为（43.1641,47.4609）， 形状大致是椭圆形的嵌套 （见图5-9）。问题二的10 个点吸收率 （附件4 坐标 点自上到下）分别为：0、0、2.0053、1.9714、1.9613、1.9644、0、0、0、0。 图片见本文4.2.2 模型求解部分。问题三的10 个点吸收率（附件4 坐标点自上 到下）分别为：0、0、0.6479、0、0、4.7242、、0、10.5331、0、0。图片轮廓 形状类似于组织细胞图，图片见本文4.2.3 模型求解部分。 对于问题四，问题1 中的相邻探测器间距离参数，是影响成像精度的重要参 数，距离越小精确度越高。旋转中心在托盘上的位置和探测器的个数，是影响成 像稳定性的重要参数，增加单位长度的探测器个数，并且调整旋转中心和探测器 中心相对位置，使得探测器每转一个角度，接收到信息探测器尽可能在中部（避 免下次转过一个角度时，介质的投影信息溢出光屏）。探测器绕旋转中心扫过的 方向角度尽量稍微大于180°，使信息采集的更加完备，从而增加了CT 工作的 稳定性。 关键词：余弦轨迹 滤波反投影法 中心切片定理,傅里叶变换 傅里叶级数分 析 1 一、问题重述 1.1 CT (Computed Tomography)可以在不破坏样品的情况下，利用样品对射线能 量的吸收特性对生物组织和工程材料的样品进行断层成像，由此获取样品内部的 结构信息。一种典型的二维CT 系统如图1 所示，平行入射的X 射线垂直于探测 器平面，每个探测器单元看成一个接收点，且等距排列。X 射线的发射器和探测 器相对位置固定不变，整个发射-接收系统绕某固定的旋转中心逆时针旋转180 次。对每一个X 射线方向，在具有512 个等距单元的探测器上测量经位置固定不 动的二维待检测介质吸收衰减后的射线能量，并经过增益等处理后得到180 组接 收信息。 CT 系统安装时往往存在误差，从而影响成像质量，因此需要对安装好的CT 系统进行参数标定，即借助于已知结构的样品（称为模板）标定CT 系统的参数， 并据此对未知结构的样品进行成像。 1.2 问题的提出 1.在正方形托盘上放置两个均匀固体介质组成的标定模板，其中每一点的数 值反映了该点的吸收强度，这里称为“吸收率”。请根据这一模板及其接收信息， 确定CT 系统旋转中心在正方形托盘中的位置、探测器单元之间的距离以及该CT 系统使用的X 射线的180 个方向。