dfr重建图像matlab,CT系统的参数标定及成像建模

CT(Computed Tomoraphy)可在不破坏样品的情况下，利用样品对射线能量的吸收特性对生物组织和工程材料的样品进行断层成像，由此获取样品的结构信息[1]。CT系统安装时往往存在误差，从而影响成像质量，因此，需要对已安装的CT系统进行参数标定，即借助于已知结构的模板标定CT系统的参数，并据此对未知结构的样品进行成像。本文拟解决文献[1]提出的相关问题。

1　模型一的建立和求解

问题一[1]要求，根据已知结构模板及其X射线接收信息，确定CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离以及该CT系统使用的X射线的180个方向。

3.5 【指南建议】 建议对先前的胎儿、婴儿或其他相关家庭成员的不确定或缺乏信息的测序结果进行回顾性分析，因为这些结果有可能增加将来妊娠的风险。若距离上次报告分析的时间太久，而患者处于备孕或妊娠状态，应在知情同意前提下进行重新分析。

CT图像重建的原理[2-3]就是，X射线照射在被检测物体上，一部分射线能量被物体吸收而使射线强度产生衰减。X射线在穿过被检测物体时，衰减呈指数变化，遵循Lambert-Beer吸收定律。对于均匀同质物体，X射线衰减是线性的，若其厚度为x，线性衰减系数(又称吸收率)为μ(问题一中吸收率为1)，X射线穿过物体前的初始强度为I0，在经过距离 x 后其强度变为 I，则有 I=I0e-μx，显然

附件2中的数据是经过增益处理后得到的接收信息。增益处理中的增益关系为

即M=CL，其中M为增益处理后X射线总衰减值，C为增益系数，L为X射线穿透的物体厚度，在问题一中L即为X射线与椭圆和小圆相交的截距和。

以模板示意图中的椭圆中心为原点，水平方向为x轴，垂直方向为y轴，单位长度为mm，易得到椭圆和圆的方程分别为

和

设512个探测器的对称轴所在直线方程为y=kx+b，X射线与已知图形的位置情况如图1所示。

 

图1　X射线与已知图形位置关系

观察用DFR方法直接重建的图像发现存在伪影。为消除伪影用先卷积后反投影的方法(CBP)重建图像，该方法利用了傅里叶变换的时域卷积定理。时域卷积定理内容如下：若f1(t)↔F1(ω)，f2(t)↔F2(ω)，则f1(t)·f2(t)↔F1(ω)F2(ω)。此定理说明，频域中两者的积是时域中两者的卷积。

机器算法中知识库的作用是存放着搜集来的知识，通过系统地表达或模块化，使得计算机能够推论和解决问题。知识库中的知识通常包含两种。一是知识本身，即物质及概念的实体分析、彼此关系；另一种是人类专家特有的经验法则、判断与直觉[7]。

先尝试直接利用傅里叶变换(DFR)方法进行图像重建，该方法以中心切片定理为基础。中心切片定理是指某图像f(x，y)在视角为φ时的投影pφ(r)的一维傅里叶变换，是f(x，y)的二维