题目CT系统参数标定与成像.PDFVIP

题目：CT系统参数标定及成像 摘要 本文通过模板的几何信息及对应的吸收率，求得模板参数，进而利用得到的 标定参数实现介质图像的重建，根据重建的图像确定介质所在位置、几何形状及 吸收率等介质信息。在对参数标定进行精度与稳定性分析之后，优化标定模板。 针对问题1，首先考虑利用Excel 中色阶分级对附录2 的数据进行初步处理，根 据CT系统工作原理以沿射线投影方向线积分为基础，推测增益处理后结果与线 积分数据有相关性，使用Matlab编程对一定范围内散点进行拟合，得到二者线性 相关的结论。以此为基础，结合附录2数据，利用解析几何在极坐标系与直角坐 标系下不同的方程，解得线性相关系数 ，进一步得到探测器间距 ，  D 0.2768mm 角度 =119.6490 , 298.6356 ，间隔不均匀，稳定在0.9999 。随后根据CT旋 1 180 转中心在探测器中点的垂直线上这一假设，取两特殊射线方向联立，求出旋转中 心在以正方形托盘中心为原点建立的坐标系中位置为（-9.2728,5.9512 ）从而完成 对CT系统参数的标定。 针对问题2 、3，根据问题1中标定的180个投影方向以及探测器单元的距 离，利用matlab 中iradon函数对附件三进行滤波反投影成像，将成像结果进行 radon变换与附件二对比，以比较不同的滤波函数及差值方法效果，选取最好的 方法对附件三、五进行反投影成像。结果是以CT系统旋转中心为原点的分辨率 为502502 图像，利用问题1中标定的旋转中心按比例截取出362362 的正方形 托盘，再将此图像插值缩放为256256 分辨率的标准规格得到几何形状。每个 像素的大小表示相对密度，以质心表示该图像的几何位置。附件二取同样方法 得到附件一的相对密度，由于相对密度和吸收率成正比，再根据附件一的吸收 率按比例可得到问题二、三的吸收率。10个位置点将其转换为相应的栅格位置 坐标即可得到吸收率。 针对问题4 ，投影方向得出1806 的数据，按行分析其稳定性，按列求得的 平均值进行线性回归分析精度。探测器单元的距离求解和最大吸收率的选取有关， 可对比不同最大吸收率得到的单元距离分析精度，利用同一吸收率对应的单元距 离进行稳定性分析。旋转中心的确定需要利用特征值，故采用反投影方法重新确 立旋转中心的坐标，由于栅格像素本身具有宽度，而X射线无宽度，对比问题一 中的求解结果，利用两者误差是否在相应允许范围内以分析精度。最后，由于原 模型为圆形、椭圆形曲线，实际图形与栅格图形不完全符合，产生误差，降低精 度，故考虑采用正方形。利用射线入射方向水平或竖直时的特性计算模板参数， 用模板接收信息的投影长度表示，找出主要误差来源，利用误差传播定律计算出 每个参数的相对误差，用以表示参数的精度，分析其稳定性。 关键词：radon变换 解析几何 滤波反投影成像 线性回归 误差传播定律 一、问题重述 CT(Computed Tomography) 可以利用样品对射线能量的吸收特性对样品断 层成像，由此获取其内部结构信息。一种二维CT系统发射平行且等距的X射线， 垂直于探测器表面且与探测器保持相对位置不变，绕一固定旋转中心逆时针旋转 180次，而使二维待测介质的位置固定不动。探测器具有512个等距单元，可测量 待测介质吸收衰减后的射线能量，经过增益等处理可得到180组接收信息。 由于CT系统安装存在误差，需对安装好的CT系统利用模板进行参数标定， 才可对位置结构的样品成像。 以此为背景，建立相应的数学模型和算法，解决如下4个问题： (1) 标定模板由两个均匀固体介质组成，放置于正方形托盘上，几何信息已知， 每点的吸收强度由附件1 中256 ×256方格内数值表示，接受信息见附件 2 。由此确定CT 系统旋转中心在正方形托盘中的位置，探测器单元之间 的距离以及CT 系统使用的X 射线的180 个方向。 (2) 利用得到的CT 标定参数，附件3 中未知介质的接收信息，确定该未知介