(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211123054.8 (22)申请日 2022.09.15 (71)申请人 京心禾 （北京） 医疗科技有限公司 地址 102299 北京市昌平区未来科 学城南 区国家电投集团科学技术研究院有限 公司院内1号楼 216室 (72)发明人 陈思　李琨　 (74)专利代理 机构 佛山市禾才知识产权代理有 限公司 4 4379 专利代理师 刘羽波 (51)Int.Cl. G06T 11/00(2006.01) G06T 5/00(2006.01) G06T 3/40(2006.01) G06F 17/16(2006.01) (54)发明名称 一种断层图像迭代重建方法 (57)摘要 本发明提供一种断层图像迭代重建方法， 包 括： 获取目标图像的成像目标三维空间分布数据 和投影数据， 对投影数据进行降噪处理并获得降 噪后投影数据； 对成像目标三维空间分布数据到 降噪后投影数据的过程进行建模并产生系统传 输矩阵； 通过系统传输矩阵进行断层图像迭代重 建； 断层图像迭代重建的目标为结合了降噪后投 影数据匹配度与 目标图像的先验估计的代价函 数最小化问题， 在利用梯度下降法进行迭代求解 的过程中， 根据每一投影数据单元与当前目标图 像的像素的关联性强弱与匹配度强弱进行权重 函数的计算与调整。 与相关技术相比， 本发明的 断层图像迭代重建方法的重建图像的分辨率高 且信噪比高。 权利要求书2页 说明书8页 附图6页 CN 115423890 A 2022.12.02 CN 115423890 A 1.一种断层图像迭代重建方法， 其特 征在于， 该 方法包括如下步骤： 步骤S1、 获取目标图像的成像目标三维空间分布数据和投影数据， 对所述投影数据进 行降噪处 理并获得降噪后投影数据； 步骤S2、 对所述成像目标三维空间分布数据到所述降噪后投影数据的过程进行建模并 产生系统传输矩阵； 步骤S3、 通过所述系统传输矩阵进行断层图像迭代重建； 所述断层图像迭代重建的目 标为结合了降噪后投影数据匹配度与所述目标图像的先验估计的代价函数最小化问题， 在 利用梯度下降法进行迭代求解的过程中， 根据每一所述投影数据单元与当前所述目标图像 的像素的关联性强弱与匹配度强弱进行预设的权重函数计算， 其中， 所述权重函数根据所 述迭代求 解的次数进行相应的更新。 2.根据权利要求1所述的断层图像迭代重建方法， 其特征在于， 所述步骤S1中， 所述降 噪处理为二维数据降噪。 3.根据权利要求2所述的断层图像迭代重建方法， 其特征在于， 所述降噪处理 的降噪方 法采用高斯滤波、 Butter worth滤波以及应用基于数据训练的滤波函数中的任意一种； 如果 使用的是基于规则的滤波方法， 则在所述步骤S2中， 对所述降噪处理所采用的所述滤波函 数进行建模。 4.根据权利要求3所述的断层图像迭代重建方法， 其特征在于， 所述降噪后投影数据y ( θ,x,z)为角度 θ、 向量x以及向量z三个方向的函数； 所述降噪处理为对不同的角度 θ， 在(x, z)方向分别独立应用的二维降噪处 理。 5.根据权利要求4所述的断层图像迭代重建方法， 其特征在于， 所述步骤S2中， 所述系 统传输矩阵为A， 并满足 公式： Ax＝y(1)； 其中， 向量y为所述降噪后投影数据的向量， 向量x为所述目标图像的向量。 6.根据权利要求5所述的断层图像迭代重建方法， 其特征在于， 所述系统传输矩阵A 的 第i行、 第j列的元素为Aij， 所述元素Aij为图像向量的第j个体素的信号对投影数据第i个像 素的贡献概率； 所述建模根据元素Aij对衰减、 散射以及导致系统分辨率降低的因素进行数 学建模。 7.根据权利要求6所述的断层图像迭代重建方法， 其特征在于， 所述步骤S3中， 所述断 层图像迭代重建将公式(1)方程 求解转化为一个代价 函数最小化的问题， 并满足如下公式： 其中， 为所述目标图像的向量x的最优化估计， argminx为寻找最小化上述公式的x， L (y， Ax)为实测投影数据与基于当前所述目标图像x估计的投影数据的匹配程度； fp(x)为对 所述目标图像x的图像估计的先验概 率或者惩罚项， 用于约束图像中的噪声和伪影。 8.根据权利要求7所述的断层图像迭代重建方法， 其特征在于， 所述步骤S3中， 将公式 (2)的最优化问题通过梯度下降类的数值方法求 解， 再经过多次迭代， 并满足如下公式： 其中， 为第n次迭代的图像估计， 为对x求梯度函数， 为沿梯度方向搜索的权　利　要　求　书 1/2 页 2 CN 115423890 A 2步长。 9.根据权利要求8所述的断层图像迭代重建方法， 其特征在于， 对于每一个体素j， 公式 (2)等效如下公式： 定义 为采样角度θ方向的投影数据中与xj相关的单元集合， 定义 为图像体素 中 与yi相关的单 元集合， 定义 为： 其分量形式为： 公式(2a)进一 步表述为： 10.根据权利要求9所述的断层图像迭代重建方法， 其特征在于， 所述步骤S3 中， 将公式 (6)引入权 重函数 以提升图像分辨 率， 如以下公式： 权　利　要　求　书 2/2 页 3 CN 115423890 A 3