ct系统与matlab成像,基于MATLAB的CT系统参数标定及成像研究

92|电子制作2017年12月 实验研究 CT(Computed Tomography)，即电子计算机断层扫描， 它是利用精确准直的X射线、超声波等，与灵敏度极高的 探测器一同围绕人体的某一部位作一个接一个的断面扫描， 具有扫描时间快，图像清晰等特点。CT用X射线照射人体， 由于人体内不同的组织或器官拥有不同的密度与厚度，故其 对X射线产生不同程度的衰减作用，从而形成不同组织或 器官的灰阶影像对比分布图，进而以病灶的相对位置、形状 和大小等改变来判断病情。本文主要借助赛题中的样品的数 据(称为模板)建立相应的数学模型和算法，标定CT系数 的参数，并据此对未知结构的样品进行成像。 1.模型假设 (1)假设CT转动时，发射-接收系统，始终垂直于 正方形托盘； (2)假设光束波动不明显，即忽略光的衍射现象， CT 的射线源发出的是能量相等的射线束，经过对射线束处理里 面的光线是相互平行的； (3)依据附表2中第1列和第180列的数据在 MATLAB程序中的成的图形可知，两个图形是近似对称，假 设两列数据对应的X射线相差 π ，且发射器是等角度逆时 针旋转，且每次旋转180 π ； (4)假设 = 2 πθ 时，X射线恰好与椭圆长轴平行，且射 线过椭圆的中心，接收信息即投影取最大值。 2.模型的建立与求解 ■■2.1■问题一 以正方形托盘中心(即椭圆中心)为原点建立笛卡尔直 角坐标系，以横轴为x轴，纵轴为y轴，角 θ 为X射线与x 轴正向的夹角，如图1所示。发射器与具有512个等距单位 的探测器的位置固定不变，整个系统的旋转中心为( 00 ,xy)。 当射线穿过均匀物质时，其吸收强度 是按指数规律衰 减，则 图1 0 l I I e µ −= ， 其中l为射线穿过物质的直线长度； µ 为被测物质的吸 收强度，即 “吸收率”。 若多物质分段均匀，各段物体的“吸收率”为分别为 ( )1,2,i in µ =  ，线性长度分别为： i l ( )1,2,in=  ，其吸收 强度I： 1 1 2 2( + ) nn l l lI I e µ µ µ −++ ， 一般地，如果物体在xy− 平面内都不均匀，“吸收率” 为 ( ),f x y，在某一方向沿某一路径L的射线吸收强度I为： ( ), 0 L f x y dl I I e −∫= ， 射线吸收强度I经过取负对数后，计为p， 0 ln I p I =   称为射线穿透物体后的投影，增益后即为附表2，附表 3，附表5中的数据。 通过上述内容知：附表2中的数值是CT成像中的投影 p增益后的数据，当 θ = 2 π 时，依据图1可知射线穿过物 体达到最大长度为80mm,其投影增益值达到最大值，根据 MATLAB程序对附表2数据分析知，元(223,151)最大值为 141.7794，即发射器第151次发射，通过对151列数据分 析知,1~45束的数据都为0，46~74束X射线经过图1中的 圆形图案，75~168束的数据都为0，第169~277束X射线 基于MATLAB的CT系统参数标定及成像研究 作者/孔建霞、曹莉、薛康、杨素青、谢桃枫、贾荣，内蒙古医科大学计算机信息学院 基金项目：内蒙古医科大学教学改革研究课题(NYJGB2201507),内蒙古医科大学实验室开放基金开放性项目(2017KF15)。摘要：本文基于2017年全国大学生数学建模大赛A题所给出的条件，在模型假设的基础上运用MATLAB中plot命令建立各点吸收强度的空间 图形、水平及垂直入射的X射线吸收强度图形，根据附件2中投影的变化，参