作者：康 宁 赵汗青 李霖龙
楼钟>金牌班长
来源：《现代经济信息》 2018年第14期
引言
二维CT 系统由发射器- 探测器系统组成，发射器发射一组平行的x 射线，经过一段固定的距离，被相应的模板吸收能量后，垂直地被等距排列的探测器接受。探测器上接受的值为经位置固定不动的二维待检测介质吸收衰减后的一个表征射线能量的值，再经过增益等处理后得到一个最终值。整个系统围绕旋转中心旋转180次，对每一个方向，512 个探测器都会返回一个最终值，形成512个数据，再经过180 次旋转，得到180 组信息。 由于CT 系统在安装时存在一些误差，可借助已知模板对系统进行参数标定，并根据参数标定结果对未知样品进行图像重建。本文主要研究了CT 系统成像误差以及信号吸收率的问题，做出了CT系统参数标定和图像重建。
一、CT 系统参数标定
1. 探测器间距
选取第一次旋转的数据，显然不为 0 的探测值是由于对应的X 射线穿过了圆形介质或者椭圆形介质导致的。选取某一合适方向，利用Matlab，根据经过小圆的射线，我们可得出经过小圆的29 个探测器的间隔距离，因为小圆的直径被这29 个探测器所扫描，且其直径与x 射线垂直，所以探测器元件的间隔为：8mm/28mm=0.2857mm 2. 探测器旋转方向
通过对每一方向512 个探测器所接受到的信息，绘出180 幅根据接受信息所拟出的图像，为了进一步研究，利用线性回归、最小二乘法的办法，将本来经过圆的射线的接受信息，变成了按原路径发射但是不受圆介质的吸收效应的影响的接受信息，使得图像只保留X 射线经过椭圆的对应接受信息的拟合图像。先建立一个合理的参考系，以椭圆的中心为原点，将穿过椭圆中心和小圆中心的直线定义为x 轴，将穿过椭圆长轴的直线定义为y 轴。辛普森悖论
接受信息与相应路程d 有如下关系。
J=1.7722×d( 式1.1)
运用穷尽算法的思想[1]，拟合出椭圆大量的中心半径与其相应的斜率之间的关系，并应用matlab 绘
毛细管电谱图。通过对每个方向的曲线分析得到该方向上的极大值，求得相应的斜率，并用反三角函数即得到相应的角度，即方向。乌克兰危机来龙去脉
3. 探测器旋转中心
利用反雷登变换[2] 来实现图像的重建。在Matlab 中运用irandon 函数实现逆变换。程序运行之后得到了下图：
通过确定扫描图像对角线交点与椭圆中心的相对位置来确定原图中的旋转中心。首先画出图像对角线交点以及椭圆长轴与短轴的交点，然后通过测量与计算我们可以得出在变换以后图像中心相对于椭圆中心的相对位置，此相对位置变换到原图就是旋转中心。
在图中作出四条辅助线，其中红线为图像的对角线，黄点为对角线交点，蓝线为以椭圆的两个轴延伸的直线。其中通过辅助线我们可以得出黄点相对于蓝线交点的相对位置为 (-10.3697,6.6461)，所以旋转中心的位置为(-10.3697,6.6461)。
二、CT 系统图像重建
CT 成像的本质是衰减系数成像，通过对接收器各个方向的能量值进行滤波反投影成像，可以求出该介质各体素的衰减系数值，获得衰减系数的二维分布，即衰减系数矩阵。根据CT 值的定义以及CT 值与衰减系数的关系，将衰减系数矩阵转化为CT 值矩阵，并利用MATLAB 将各像素的CT 值转化为灰度，生成CT 影像。
1. 求解衰减系数分布函数
根据滤波反投影重建算法[3-8] 推导出衰减系数分布函数的求解过程：
假设物体在xy 平面内都不均匀，衰减系数f=f(x,y) 为被测物体某断面的线性衰减系数的分布函数。当初始强度为I0 的x 射线穿过物体时，探测器检测到的I 服从Beer 定理：
施韦泽
三、结语
根据上述方法，可以通过已知模板信息对安装好的CT 系统进行参数标定，并基于滤波反投影算法，利用Matlab 软件，进行CT系统的图像重建。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议