[发明专利]一种CT图像重建方法及系统、终端及可读存储介质

权利要求书

1.一种CT图像重建方法，其特征在于，所述方法包括：

原始重建图像优化：对原始重建图像进行去伪去噪优化，获得优化的原始重建图像；

差值数据运算：对所述原始重建图像和所述优化的原始重建图像进行求差运算，获得差值数据；对CT扫描数据进行原始的CT图像重建计算，获得原始重建图像；

特征算子提取：对所述优化的原始重建图像进行图像特征提取，获得特征算子；所述特征算子Fⁿ⁺¹通过如下方式获得：其中，M表示所述以像素点i为中心的尺寸为的图像块的宽度，D表示所述像素点i到所述的图像块的边缘的距离；以像素点i为中心的尺寸为的图像块所对应的子特征算子fⁿ⁺¹通过如下方式获得：其中，C表示常量；

CT图像重建：对所述差值数据与所述特征算子进行卷积，并将得到的卷积数据添加至所述优化的原始重建图像中，获得CT重建图像；所述CT图像重建步骤具体根据下述公式获得CT重建图像：其中，表示所述CT重建图像，μⁿ表示所述原始重建图像，μⁿ⁺¹表示所述优化的原始重建图像，Fⁿ⁺¹表示所述特征算子，表示点乘法，A表示非负控制符号，A满足

按照预置的图像标准，对所述CT重建图像进行达标判断；所述预置的图像标准由SSIM算法确定；

若确定所述CT重建图像达标，则确定当前的所述CT重建图像为最终的CT重建图像；

若确定所述CT重建图像未达标，则将所述CT重建图像作为所述原始重建图像，返回执行所述原始重建图像优化步骤。

2.如权利要求1所述的方法，其特征在于，所述特征算子提取步骤具体包括如下步骤：

对所述优化的原始重建图像做线性高斯滤波处理，获得模糊图像；

分别从所述优化的原始重建图像和所述模糊图像中，提取以像素点i为中心的尺寸为的图像块，分别获得优化的原始重建图像的子图像块和模糊图像的子图像块；其中，N表示像素数量，为正整数，i的初始值为i＝[(N-1)/2]+1；

分别对所述优化的原始重建图像的子图像块和所述模糊图像的子图像块进行均方差计算，获得所述优化的原始重建图像的子图像块的均方差和所述模糊图像的子图像块的均方差，并求得所述优化的原始重建图像的子图像块与所述模糊图像的子图像块之间的均方差；

利用所述优化的原始重建图像的子图像块的均方差、所述模糊图像的子图像块的均方差，以及所述优化的原始重建图像的子图像块与所述模糊图像的子图像块之间的均方差，获得以像素点i为中心的尺寸为的图像块所对应的子特征算子；

令i＝i+1，返回执行所述分别从所述优化的原始重建图像和所述模糊图像中，提取以像素点i为中心的尺寸为的图像块，获得优化的原始重建图像的子图像块和模糊图像的子图像块操作，依次迭代循环，直至i＝W-[(N-1)/2]退出，获得若干个所述子特征算子；其中，W表示所述优化的原始重建图像或所述模糊图像的宽度；

对所述若干个所述子特征算子求和，获得所述特征算子。

3.如权利要求2所述的方法，其特征在于：

所述优化的原始重建图像的子图像块的均方差σ_p(i)通过如下方式获得：其中，P(i)表示所述优化的原始重建图像的子图像块的像素的均值，p_i表示所述优化的原始重建图像的子图像块，μⁿ⁺¹表示所述优化的原始重建图像；

所述模糊图像的子图像块的均方差σ_q(i)通过如下方式获得：其中，Q(i)表示所述模糊图像的子图像块的像素的均值，q_i表示所述模糊图像的子图像块，μ_d表示所述模糊图像；

所述优化的原始重建图像的子图像块与所述模糊图像的子图像块之间的均方差σ_qp(i)通过如下方式获得：