具体技术细节
[0004] 本发明针对背景技术中存在的问题，提出了一种基于共轭梯度的多约束模型求解方法，设计增加两个对角子带的变分信息来减少阶梯效应的双全变分(DTV)方法。此外，申请人设计了不同的加速策略来优化共轭梯度，包括快速共轭梯度(FCG)和改进的自适应矩估计共轭梯度(IADAMCG)。
[0005] 技术方案：
[0006] 一种基于多方向全变分正则化和改进自适应矩估计的快速磁共振图像重建方法，该方法基于以下磁共振图像重建数学模型：
[0007]
[0008]
[0009] 其中，编码矩阵A＝ΦΨT，Φ代表k‑空间采样模式，Ψ是正交矩阵；θ是投影系数矩阵；Y是通过欠采样测量获得的k‑空间数据；λ是正则化参数；ω1表示TV(θ)项权重，ω2表示TVdiag(θ)项权重；TV(θ)表示水平和垂直子带方向的变分，TVdiag(θ)表示对角子带方向的变分；
[0010] ε代表测得数据的噪声水平；
[0011] 具体包括以下步骤：
[0012] S1、输入编码矩阵A，在k‑空间中的欠采样MR图像数据Y，是正交矩阵Ψ；
[0013] S2、数学模型更新，引入一阶矩mt和二阶矩vt：
[0014] mt＝β1·mt‑1+(1‑β1)·gt
[0015]
[0016] 一阶矩mt表示梯度的均值，用于控制模型更新的方向；二阶矩vt表示梯度的方差，用于控制迭代步长；β1和β2分别是一阶和二阶矩的衰减系数；gt表示梯度在指向目标函数后的值；
[0017] S3、对一阶矩mt和二阶矩vt校正偏差获得校正后一阶矩 和校正后二阶矩[0018] S4、更新投影系数矩阵θt+1和gt；
[0019] S5、t＝t+1，并判断t是否等于设定阈值，否则返回S4，是则转入S6；
[0020] S6、基于更新后的数学模型，输出重建磁共振图像X＝ΨTθ。
[0021] 优选的，S2中，二阶矩vt的表达式改写如下：
[0022] 优选的，m1＝0,v1＝0。
[0023] 优选的，β1＝0.5,β2＝0.5。
[0024] 优选的，更新投影系数矩阵θt+1：
[0025]
[0026] 式中，γ是学习速率，Δ是平滑项。
[0027] 优选的，γ1＝0.05,Δ＝10‑8。
[0028] 优选的，γt通过下式进行指数衰减：
[0029] γt＝δtγt‑1
[0030] 其中，δ表示衰减因子。
[0031] 优选的，δ＝0.99。
[0032] 优选的，S4中，gt的更新公式如下：
[0033] gt＝gt‑ηtgt‑1
[0034] 式中， 是一个标量，表示方向调节参数。
[0035] 具体的：
[0036] TV(θ)表示为：
[0037]
[0038] 其中， 表示水平方向变分， 表示垂直方向变分；
[0039] TVdiag(θ)表示为：
[0040]
[0041] 其中， 表示45度方向变分， 表示135度方向变分。
[0042] 本发明的有益效果
[0043] 本发明提出的一种基于多方向全变分正则化和改进自适应矩估计的快速磁共振图像重建方法IADAMCG_DTV。通过比较重建图像、PSNR、TEI、NMI和迭代时间等方面的CG系列、FCG系列和IADAMCG系列算法的性能，IADAMCG_DTV算法保留了更多重建图像的全局和局部信息。
[0044] 本发明的基于多方向全变分正则化和改进自适应矩估计的快速磁共振图像重建方法通过接入磁共振成像系统中，不仅可以提升磁共振成像系统的吞吐量和检查效率，减少患者等待时间，还显著提高了医院的经济效益和诊断准确性。此外，它还促进了相关科研发展，推动医疗影像技术的进步，对医疗行业的发展具有重要意义。
法律保护范围
涉及权利要求数量10：其中独权1项，从权-1项
1.一种基于多方向全变分正则化和改进自适应矩估计的快速磁共振图像重建方法，其特征在于该方法基于以下磁共振图像重建数学模型：
T
其中，编码矩阵A＝ΦΨ ，Φ代表k‑空间采样模式，Ψ是正交矩阵；θ是投影系数矩阵；Y是通过欠采样测量获得的k‑空间数据；λ是正则化参数；ω1表示TV(θ)项权重，ω2表示TVdiag(θ)项权重；TV(θ)表示水平和垂直子带方向的变分，TVdiag(θ)表示对角子带方向的变分；
ε代表测得数据的噪声水平；
具体包括以下步骤：
S1、输入编码矩阵A，在k‑空间中的欠采样MR图像数据Y，是正交矩阵Ψ；
S2、数学模型更新，引入一阶矩mt和二阶矩vt：
mt＝β1·mt；1+(1‑β1)·gt
一阶矩mt表示梯度的均值，用于控制模型更新的方向；二阶矩vt表示梯度的方差，用于控制迭代步长；β1和β2分别是一阶和二阶矩的衰减系数；gt表示梯度在指向目标函数后的值；
S3、对一阶矩mt和二阶矩vt校正偏差获得校正后一阶矩 和校正后二阶矩
S4、更新投影系数矩阵θt:1和gt；
S5、t＝t+1，并判断t是否等于设定阈值，否则返回S4，是则转入S6；
T
S6、基于更新后的数学模型，输出重建磁共振图像X＝Ψ θ。
2.根据权利要求1所述的方法，其特征在于S2中，二阶矩vt的表达式改写如下：
3.根据权利要求1所述的方法，其特征在于：m1＝0,v1＝0。
4.根据权利要求1所述的方法，其特征在于：β1＝0.5,β2＝0.5。
5.根据权利要求1所述的方法，其特征在于S4中，更新投影系数矩阵θt:1：
式中，γ是学习速率，Δ是平滑项。
；8
6.根据权利要求5所述的方法，其特征在于γ1＝0.05,Δ＝10 。
7.根据权利要求5所述的方法，其特征在于γt通过下式进行指数衰减：
t
γt＝δ γt；1
其中，δ表示衰减因子。
8.根据权利要求7所述的方法，其特征在于δ＝0.99。
9.根据权利要求1所述的方法，其特征在于S4中，gt的更新公式如下：
gt＝gt‑ηtgt；1
式中， 是一个标量，表示方向调节参数。
10.根据权利要求1所述的方法，其特征在于：
TV(θ)表示为：
其中， 表示水平方向变分， 表示垂直方向变分；
TVdiag(θ)表示为：
其中， 表示45度方向变分， 表示135度方向变分。
