技术特征：
1.一种基于行与列同时低秩约束的磁共振波谱重建方法，包括以下步骤：(1)输入一个二维磁共振波谱信号的时空编码数据y和欠采样矩阵u，对于待重建的具有n个行向量和m个列向量的二维磁共振波谱信号x，其第n行中第m列的元素表示为x
n，m
；(2)针对x中第m列的向量x
m
，构建其对应的汉克尔矩阵可表示为：其中，表示将列向量转换为汉克尔矩阵的线性算符，l表示尺寸参数；(3)针对x中第n行的向量x
n
，构建其对应的汉克尔矩阵可表示为：其中，表示将行向量转换为汉克尔矩阵的线性算符，k表示尺寸参数；(4)采用非凸的h
ε
范数约束构建的汉克尔矩阵的低秩特性，表示任意矩阵l的h
ε
范数，其定义式为：其中，σ
r
为l的第r个奇异值，ε＞0是一个小常数；(5)在利用h
ε
范数同时约束x中行与列构建的汉克尔矩阵低秩特性的基础上，建立关于二维磁共振波谱的重建模型：其中，f
m
表示m
×
m大小的傅里叶变换矩阵，i
n
表示n
×
n大小的单位矩阵，表示矩阵的克罗内克积算符，vec(
·
)表示将矩阵的各列依次堆叠成向量的算符，表示向量的二范数的平方，λ表示时空编码数据保真项的正则化参数；(6)在重建模型中引入中间变量和并根据交替方向乘子法求解重建模型中各优化变量的子问题：(6a)固定x和r
n
，求解关于l
m
的子问题：其中，β表示惩罚参数，t
m
表示针对l
m
引入的拉格朗日乘子，表示矩阵的frobenius范数的平方；(6b)固定x和l
m
，求解关于r
n
的子问题：
其中，z
n
表示针对r
n
引入的拉格朗日乘子；(6c)固定l
m
和r
n
，求解关于x的子问题：(7)在得到各优化变量的解后，更新拉格朗日乘子t
m
和z
n
，以及惩罚参数β：，以及惩罚参数β：β＝τβ其中τ＞1表示增长因子，重复步骤(6)～(7)，直到重建的磁共振波谱满足条件或迭代次数达到预设上限。2.根据权利要求1所述的一种基于行与列同时低秩约束的磁共振波谱重建方法，其特征在于，步骤(6a)和(6b)中求解关于l
m
和r
n
的子问题，可以按照以下步骤进行：(6a1)对矩阵进行奇异值分解，得到：其中，p表示左奇异向量矩阵，v表示右奇异向量矩阵，δ表示奇异值矩阵；(6a2)计算关于矩阵奇异值的阈值参数μ：(6a3)对奇异值矩阵δ中所有对角元素进行阈值处理：其中，th(
·
)表示阈值函数，δ表示δ中任意对角元素；(6a4)对δ中所有对角元素阈值处理后的结果可表示为th(δ)，则关于l
m
的子问题的最优解表达式为：l
m
＝pth(δ)v
h
(6b1)对矩阵进行奇异值分解，得到：其中，e表示左奇异向量矩阵，q表示右奇异向量矩阵，σ表示奇异值矩阵；(6b2)通过阈值函数th(
·
)对奇异值矩阵σ中所有对角元素进行阈值处理，结果可表示为th(σ)，则关于r
n
的子问题的最优解表达式为：r
n
＝eth(σ) q
h
。3.根据权利要求1所述的一种基于行与列同时低秩约束的磁共振波谱重建方法，其特征在于，步骤(6c)中求解关于x的子问题，可以按照以下步骤进行：(6c1)该子问题是一个大规模的最小二乘问题，为快速计算其近似解，需定义与
的共轭运算，其定义式为：的共轭运算，其定义式为：其中，表示的共轭算符，a表示任意l
×
(ν-l+1)大小的矩阵，a
k，i
表示a中第k行中第i列的元素，表示的共轭算符，b表示任意k
×
(m-k+1)大小的矩阵，b
k，i
表示b中第k行中第i列的元素，min(
·
，
·
)表示最小值函数，max(
·
，
·
)表示最大值函数，(
·
)
t
表示矩阵的转置算符；(6c2)利用共轭算符与计算辅助矩阵j1与j2：：另外，计算ν
×
ν大小的对角权重矩阵与m
×
m大小的对角权重矩阵其中的第i个对角元素与的第i个对角元素的表达式为：的表达式为：(6c3)引入中间变量s，关于中间变量s的计算式为：其中，γ＞0是一个小常数，i
mn
表示mn
×
mn大小的单位矩阵，i
m
表示m
×
m大小的单位矩阵，i
n
表示n
×
n大小的单位矩阵，表示矩阵的克罗内克积算符，(
·
)-1
表示矩阵的逆，x
last
表示上一次迭代中该子问题的解，vec(
·
)表示将矩阵的各列依次堆叠成向量的算符；(6c4)在得到中间变量s后，关于x的子问题的近似解表达式为：其中，是m
×
m大小的傅里叶变换矩阵f
m
的共轭转置矩阵，u
h
是欠采样矩阵u的共轭转置矩阵，将得到的vec(x)转换为矩阵x即为获得的近似解。
技术总结
本发明公开了一种基于行与列同时低秩约束的磁共振波谱重建方法。属于数字图像处理技术领域。它是一种利用二维磁共振波谱的二元指数结构对其行与列分别构成的汉克尔矩阵进行低秩约束的磁共振波谱重建方法。首先将二维磁共振波谱的各行与各列构建对应汉克尔矩阵，然后采用一种非凸范数约束汉克尔矩阵的低秩特性，最后通过交替方向乘子法优化求解建立的磁共振波谱重建模型。本发明对二维磁共振波谱施加非凸的低秩约束，有效提取了其二元指数结构，同时分别对重建模型的各子问题进行精确求解，使得从欠采样的时空编码数据中重建的磁共振波谱有效抑制了伪影现象，而且重建波谱中的谱峰更接近于真实的谱峰，因此可用于磁共振波谱的重建。谱的重建。谱的重建。
技术研发人员：刘书君 曹建鑫 田新雨 张奎 张新征
受保护的技术使用者：重庆大学
技术研发日：2022.07.20
技术公布日：2022/9/23