SVD Compression for Magnetic Resonance Fingerprinting in the Time Domain,本文源自IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 33, NO. 12, DECEMBER 2014,
文章可以从这里下载,这里主要写写本人对这篇文章的翻译与理解。
摘要
磁共振指纹是一种获取和处理磁共振数据的技术,它可以同时通过模式识别算法同时提供不同组织参数的定量图谱。一个预先定义的字典模型可以通过计算观测信号和字典中的每一个信号的内积来实现。这个字典以磁共振参数和模式识别算法的不同组合,使用布洛赫方程方程来仿真可能的信号演变。尽管这一匹配算法已经表明能够准确地预测感兴趣的磁共振参数,我们希望能有一个更有效的方法来获取定量图像。我们提出采用奇异值分解来压缩字典,这将会提供一个低秩逼近方法。通过压缩时间域内的字典大小,我们能够在不牺牲前文描述的原始方案的高信噪比的情况下,加速模式识别算法,加速因子为3.4到4.8。
关键词
降维,磁共振成像,模式识别与分类,奇异值分解
介绍
磁共振指纹(MRF)定量成像
字典可以用
D∈Cx×t
,其中
n
表示参数组合数,
x∗ 表示向量 x 的共轭转置,
奇异值分解背景
每个矩阵
A∈Cp×q
都可以写成奇异值分解形式,可用用下式表示
U∈Cp×p , V∈Cq×q , Σ∈Rp×q 是一个对角矩阵,它包含了不递增的奇异值
σi,i=1,...,min{p,q} 。 U 的列可用
矩阵A的秩k逼近可用矩阵的截取前k个秩对应的矩阵的和来表示,写作下式
文献[26]可以证明,下式是A最好的低秩逼近
其中下确界由所有的 p×q 矩阵 B 决定,它们的秩都小于等于
A
的总能量可以定义为它所有奇异值的平方和
能量比[6]表示秩k逼近保留的能量 A(k) 在总能量 E 中的占比。
对低秩逼近而言,能量比在决定一个适当的截断索引时很有用,它能像期望的那样保留初始矩阵的大量信息。
方法
我们将会使用模式识别算法,类似于由Turk和Pentland在文献[15]中提出的本征脸算法,人脸识别是通过将图像投影到特征向量张成的上低维子空间来执行的, 这些特征向量对应对打的特征值。
我们把奇异值分解应用到磁共振指纹的字典上,写作
这里 U,V 和 Σ 和前文部分描述的一样。在此文中令 r=rank(D) ,我们假定 n>t 。
对于一个给定的索引 k,1≤k≤r ,则(2)式中截断奇异值分解可以写成矩阵形式,产生以下字典的低秩逼近
其中 Uk=[u1,...,uk] 表示包含前k个左奇异向量,同理,适用于 Σk , Vk 。
奇异值分解的另一条性质是前r个右奇异向量
{v1,...,vr}
构成
D
的行的正交基, 也即,每一个字典条目可以写成这些正交向量的线性组合。将字典投影到由前
为了简便, 我们将把这个低维空间称为”奇异值分解空间(SVD space)“。
奇异值分解空间的模板匹配
将观测到的磁共振信号
x
投影到由
到此为止(1)式中的模板匹配就可以在奇异值分解空间中计算了。
我们搜索字典词目就变成了下式
注意到由于 U 是一个酉矩阵,积