[发明专利]一种X射线图像线性重建方法

摘要

本发明公开了一种X射线图像线性重建方法，首先获取X射线投影图像，引入基于全变差先验的正则项，构造重建目标函数；引入超先验参数，构造分层贝叶斯模型；运用变量分裂方法，引入分裂变量，分离数据保真项和TV正则项，得到分裂形式下的联合概率密度分布；基于Jefferys先验定义超先验变量，得到各变量的条件分布；迭代更新超先验参数，对含TV正则项的分裂变量的条件分布进行求解；利用正向矩阵的低秩性质，近似待求参数的全条件概率密度分布，计算低秩近似的目标分布得到关于待求参数的闭合解；计算采样样本的均值，对待求参数进行估计。本发明能够有效地解决求解大规模线性逆问题中存在的计算开销大等问题。