[发明专利]基于CT投影数据的直接迭代多材料分解图像重建方法

权利要求书

1.一种基于CT投影数据的直接迭代多材料分解图像重建方法，其特征在于，包含以下步骤：

步骤1)，针对具有n个能段的CT成像系统，建立基于CT投影数据的m种基材料直接分解图象重建模型：

其中，P表示CT成像系统采集的投影向量，n表示CT成像系统的探测能段数，n≥1，m表示基材料数，m≥2且m≥n，P(E_i)表示第i个能段的投影向量，E_i表示第i个能段的平均能量，i＝1,2,...,n，μ_j(E_i)分别为第j种基材料在能量E_i下的线性衰减系数，j＝1,2,...,m，A表示投影算子，(x₁,x₂,…,x_m)表示m种基材料的分解图像，表示联合投影算子；

步骤2)，构建基于CT投影数据的m种基材料分解直接迭代重建算法，运用该算法对m种基材料直接分解图象重建模型进行求解获得m种基材料分解图像，所述迭代重建算法按照如下公式进行迭代计算：

其中，表示m种基材料第k+1次迭代时的分解图像，表示m种基材料第k次迭代时的分解图像，k表示迭代次数，S表示矩阵B中各列元素和的向量，R^k表示正则项；

当n＝m时，其中，β₁、…、β_m为预设的正则化参数，U(·)表示能量函数；

当nm时，其中，γ₁、…、γ_m为预设的惩罚参数，表示第l种基材料的L1范数。