一种基于拟合精确度加权步长的弥散加权成像运动矫正优化方法
【专利摘要】本发明公开了一种基于拟合精确度加权步长的弥散加权成像运动矫正优化方法，依次包括：(1)利用IVIM双指数模型对采集的多b值图像进行量化参数估计，获得各个量化参数值(D,D*,f)及拟合精确度(R2)；(2)再利用获得的拟合精确度于对应位置上优化步长进行加权；(3)根据不同位置上拟合精确度加权的步长来对我们构建的目标函数进行优化，获得最终的运动矫正变形场结果；本发明对多b值的弥散加权成像运动矫正效果更好。
【专利说明】
一种基于拟合精确度加权步长的弥散加权成像运动矫正优化 方法
技术领域
[0001] 本发明属于磁共振图像处理领域，具体涉及不同模态磁共振图像间的运动矫正优 化方法。
【背景技术】
[0002] 磁共振成像(Magnetic Resonance Imaging，MRI)由于其无电离福射、能够获得丰 富的组织对比度信息以及具有非入侵式检测等优点，已广泛应用于临床医学影像检查。弥 散加权成像(Diffusion Weighted Imaging，DWI)因其能够在分子学水平上观测到组织内 水分子的运动而获得十分广泛的应用。定性和定量的分析水分子的扩散运动越来越受到广 泛的关注。然而要获得定量的结果需要不同程度加权的弥散图像，然后对不同程度加权的 弥散图像进行参数估计。现如今存在许多对水分子进行定量的参数估计模型，然而这些模 型都要求不同程度加权的弥散图像间要完全对齐。由于磁共振数据采集较为缓慢，不同程 度加权的弥散图像间由于心跳，呼吸等影响难免会存在运动。运动矫正是获得准确的量化 参数估计的前提。呼吸触发的图像采集方式能够一定限度上抑制运动，但是呼吸触发的方 式需要更多的采集时间，并且呼吸不规律时，运动会同样存在。图像配准就成为了不可或缺 的预处理步骤，然而，弥散加权图像的信号会随着b值的增大而衰减，这就导致不同b值图像 间信号强度存在较大的差异。现存的各种配准的模型方法，不能很好地矫正弥散加权图像 的运动。
【发明内容】
[0003] 本发明的目的在于针对现有技术采集到的弥散加权图像不同b值间存在运动，导 致量化参数估计的不准确性。用本发明优化方法来矫正运动，从而获取更加准确的量化参 数。
[0004] 本发明的上述目的通过如下技术方案实现：
[0005] -种基于拟合精确度加权步长的弥散加权成像运动矫正方法，依次包括如下步 骤：
[0006] (1)利用IVIM双指数模型对采集的多弥散权重系数值b图像进行量化参数估计，获 得水分子扩散系数D、微循环血流灌注系数D*、灌注分数f及拟合精确度R 2;
[0007] (2)再利用获得的拟合精确度计算步长权重，于对应位置上对步长进行加权操作；
[0008] (3)根据不同位置上拟合精确度加权的步长来对我们构建的目标函数进行优化， 获得最终的运动矫正变形场结果。
[0009] 作为一种优选的技术方案，所述步骤(1)中，根据采集数据点和IV頂双指数模型的 关系Sb/S〇 = f · e-b'D*+(l-f) · e-b'D拟合获得各个参数值(0,0*彳），其中沁为与弥散加权程 度相关的值，S b为对应b值观测到的信号强度，So为b = 0时的信号强度;根据采集数据点和拟 合数据点的关系来求得拟合精确度
4&为采集数据点和拟合数据点的离差平 方和，SST为采集数据点的平均值。
[0010]作为一种优选的技术方案，步骤（2)中，根据拟合精确度与步长权重的关系 ?(i?2) =e-μ2计算获得权重矩阵ω (R2)，其中R2为拟合精确度，h为用于放大拟合精确度差异 的因子。
[0011]作为一种优选的技术方案，步骤(3)中，采用加权的权重 来更新变形场，其中T1为上一步获得的变形场，T1+1为用所求权重更新的变形场，μ为优化的 步长，目标函数的梯度。
[0012] 作为一种优选的技术方案，步骤（3)中，采用目标函数Ε(Τ)=_ΕΜ?(Α，Τ(Β))+λ · ESMClth(T)来求解运动矫正的变形场，其中ΕΜΙ(Α，Τ(Β))为图像Α与矫正后的图像Β的互信息 量，E SMClth(T)是对变形场进行约束项，λ为权重因子，用于平衡图像间互信息量和变形场的 平滑性。
[0013] 采用了上述技术方案，本发明的有益效果为:相比于现存的确定一个全局的步长 的优化方法，本发明所提出的基于拟合精确度加权的优化方法，利用有关变形场的先验信 息来对变形场进行优化，可以跳出某些局部极值，找到更优的变形场。
【附图说明】
[0014] 利用附图对本发明作进一步的说明，但附图中的内容不构成对本发明的任何限 制。
[0015] 图1为本发明一种基于拟合精确度加权步长的弥散加权成像运动矫正优化方法的 步长权重图；
[0016] 图2为肝脏结果比较图。
【具体实施方式】
[0017] 下面结合具体的实施例对本发明进行详细描述。
[0018] 如图1所示，一种基于拟合精确度加权步长的弥散加权成像运动矫正方法，依次包 括如下步骤：
[0019] 步骤一、利用IVIM双指数模型对采集的多弥散权重系数值b图像进行量化参数估 计，获得水分子扩散系数D、微循环血流灌注系数D*、灌注分数f及拟合精确度R 2。
[0020] 根据采集数据点和IVIM双指数模型的关系Sb/So = f · eTb'w+(l-f) · eTb'D拟合获 得各个参数值(D，D*，f)，其中，b为与弥散加权程度相关的值，S b为对应b值观测到的信号强 度，So为b = 0时的信号强度；根据采集数据点和拟合数据点的关系来求得拟合精确度
6&为采集数据点和拟合数据点的离差平方和，SST为采集数据点的平均值。
[0021 ]步骤二、再利用获得的拟合精确度于对应位置上通过求解目标函数来优化步长进 行加权。
[0022] 根据拟合精确度与步长权重的关系β/炉）=e^计算获得权重矩阵co(R2)，其中 R 2为拟合精确度，h为用于放大拟合精确度差异的因子。
[0023] 步骤三、根据不同位置上拟合精确度加权的步长来对我们构建的目标函数进行优 化，获得最终的运动矫正变形场结果。
[0024]
?来更新变形场，其中T1为上一步获得的变 形场，τ1+1为用所求权重更新的变形场，μ为优化的步长，目标函数的梯度。
[0025] 采用目标函数E(T)=-EMI(A，T(B))+λ·Esm。。th(T)来求解运动矫正的变形场，其中 ΕΜΙ(Α，Τ(Β))为图像Α与矫正后的图像Β的互信息量，Ε_〇*(Τ)是对变形场进行约束项，λ为权 重因子，用于平衡图像间互信息量和变形场的平滑性。
[0026] 由于我们所得到的变形场不是每一个体素的变形场，我们需要对变形场进行插 值。变形场处理后的图像未必每个点都能落在其整数位置上，我们还需要对变形后的图像 进行插值，才能获取最终的输出图像。为了获得较好的输出图像，我们采用Β样条进行数据 插值。如图2所示，为肝脏结果比较图，其中（a)为运动矫正的参考图像，（b)为黄色区域放大 的效果图，红色边界为从参考图像中提取的边界，（c)，（f)，（i)为运动矫正前b值分别为 150,400和1000的图，（d)，（g)，（ j)为全局的步长优化b值分别为150,400和1000的运动矫正 结果，（e)，（h)，（k)为采用本发明的方法获得的运功矫正结果图
[0027] 与现有技术中的优化方法相比较。本发明中利用变形场的先验信息，使用加权的 步长，能够跳出目标函数的某些局部极值点，获得更好的运动矫正效果。
[0028]经过实验证明，本发明的方法能够很好的矫正弥散加权成像中的运功。
【主权项】
1. 一种基于拟合精确度加权步长的弥散加权成像运动矫正方法，其特征在于：依次包 括如下步骤： (1) 利用IVIM双指数模型对采集的多弥散权重系数值b图像进行量化参数估计，获得水 分子扩散系数D、微循环血流灌注系数D*、灌注分数f及拟合精确度R2; (2) 再利用获得的拟合精确度计算步长权重，于对应位置上对步长进行加权操作； (3) 根据不同位置上拟合精确度加权的步长来对我们构建的目标函数进行优化，获得 最终的运动矫正变形场结果。2. 根据权利要求1所述的一种基于拟合精确度加权步长的弥散加权成像运动矫正优化 方法，其特征在于：所述步骤（1)中，根据采集数据点和IVIM双指数模型的关系拟合获得各个参数值(D，D*，f)，其中，b为与弥散加权程度相 关的值，Sb为对应b值观测到的信号强度，So为b = 0时的信号强度； 根据采集数据点和拟合数据点的关系来求得拟合精确度，SSe为采集数据 点和拟合数据点的离差平方和，SSt为采集数据点的平均值。3. 根据权利要求1所述的一种基于拟合精确度加权步长的弥散加权成像运动矫正优化 方法，其特征在于:步骤(2)中，根据拟合精确度与步长权重的关系如知）十算获得权 重矩阵ω(Κ2)，其中R2为拟合精确度，h为用于放大拟合精确度差异的因子。4. 根据权利要求1所述的一种基于拟合精确度加权步长的弥散加权成像运动矫正优化 方法，其特征在于:步骤(3)中，采用加权的权重杉更新变形场，其 中了1为上一步获得的变形场，TW为用所求权重更新的变形场，μ为优化的步长，巧恣目标函 数的梯度。5. 根据权利要求1所述的一种基于拟合精确度加权步长的弥散加权成像运动矫正优化 方法，其特征在于:步骤(3)中，采用目标函数Ε(Τ) = -Εμι(A，Τ(Β)) +λ · EsmDDth(Τ)来求解运动 矫正的变形场，其中Emi(A，T(B))为图像A与矫正后的图像B的互信息量，EsmDDth(T)是对变形 场进行约束项，λ为权重因子，用于平衡图像间互信息量和变形场的平滑性。
【文档编号】G06T5/00GK106097262SQ201610391030
【公开日】2016年11月9日
【申请日】2016年6月2日 公开号201610391030.9, CN 106097262 A, CN 106097262A, CN 201610391030, CN-A-106097262, CN106097262 A, CN106097262A, CN201610391030, CN201610391030.9
【发明人】后桂荣, 郭义昊, 卢振泰, 张静, 冯前进, 陈武凡, 许乙凯, 冯衍秋
【申请人】南方医科大学