CN102551723A - 一种多支撑向量机模型的磁共振并行成像方法 - Google Patents
一种多支撑向量机模型的磁共振并行成像方法 Download PDFInfo
- Publication number
- CN102551723A CN102551723A CN2012100120610A CN201210012061A CN102551723A CN 102551723 A CN102551723 A CN 102551723A CN 2012100120610 A CN2012100120610 A CN 2012100120610A CN 201210012061 A CN201210012061 A CN 201210012061A CN 102551723 A CN102551723 A CN 102551723A
- Authority
- CN
- China
- Prior art keywords
- alpha
- training
- parallel imaging
- imaging method
- kernel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 26
- 238000010606 normalization Methods 0.000 claims abstract description 20
- 238000005070 sampling Methods 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 10
- 238000007689 inspection Methods 0.000 claims description 21
- 238000012706 support-vector machine Methods 0.000 claims description 14
- 238000005457 optimization Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 238000012886 linear function Methods 0.000 claims description 4
- 239000004744 fabric Substances 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000009977 dual effect Effects 0.000 claims description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 10
- 230000035945 sensitivity Effects 0.000 description 4
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000004223 radioprotective effect Effects 0.000 description 1
- 238000009394 selective breeding Methods 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种多支撑向量机模型的磁共振并行成像方法,属于磁共振并行成像领域,它的步骤如下:(1)对K空间中间区域进行全采样,将其划分为训练集和检验集,其它区域加速采样后作为预测集,并对各集合内的数据进行归一化处理;(2)将训练集划分为多组训练子集,利用支撑向量机,选择不同的参数对各训练子集进行训练,得到不同的联合权重函数模型;(3)在检验集上,对各联合权重函数进行检验,选择最佳的几个子模型;(4)分别使用最佳的几个子模型对预测集进行预测,取其平均值作为未采集点的值,反归一化处理后,将K空间数据转换为图像。上述并行成像方法,利用支撑向量机拟合的权重函数具有良好的泛化能力,整体重建误差较小。
Description
技术领域
本发明属于磁共振并行成像领域,具体涉及一种多支撑向量机模型的磁共振并行成像方法。
背景技术
磁共振成像(Magnetic resonance imaging,MRI)由于其无核辐射,分辨率高,能多方位和多参数成像等优点,已成为临床医学影像检查的重要手段之一。然而受到傅里叶编码方式和奈奎斯特采样定理限制,磁共振成像速度较慢,这不但给患者带来一定的不适,而且容易产生运动伪影。同时,长的扫描时间限制了MRI对运动物体的成像,如婴儿,血流,心脏等。经过近几十年的发展,依靠提高硬件性能来加速采集的方式以达到人体承受极限。
并行成像技术使用多个线圈同时采集信号,利用各个线圈的空间敏感度来代替部分傅里叶编码,从而达到减少扫描时间。磁共振并行成像方法有很多种,其中灵敏度编码和广义自校准并行采集(Generalized Auto-calibrating Partially Parallel Acquisitions,GRAPPA)是临床上较为常用的两种。GRAPPA算法是假设线圈敏感度的线性联合能产生空间谐波,达到傅里叶编码的效果。首先,通过中间区域全采样的数据,作为校准数据,通过构造线性方程组,求解线圈联合权重函数。再对其它区域未采集到的点进行联合权重填充。对于第l个线圈坐标位置在(kx,ky+mΔky)的点,可以通过公式:
其中kx,ky为K空间频率编码和相位编码方向上的坐标值,m为相位编码方向上的偏移量,Δky为相位编码方向上单位间隔,s代表K空间信号值,Nb and Na为联合的相邻的相位编码方向上所用的联合权重行数(分别上和下),Ha和Hb为频率编码方向所用的联合联合权重列数(分别为左和右),j代表线圈编号,h为重建的频率编码方向所用联合权重点序号,b为相位编码方向所用的联合权重点序号,w为联合权重函数,R为采样加速因子。为了得到权重函数,GRAPPA算法在中间区域全采样,然后假设一些采集到的点为需要拟合的点,组成一个知道输入和输出的线性系统方程:
AW=B
通过最小二乘法求解得:W=argmin||B-AW||2
实际上,由于联合权重函数的维数未知,只能通过人为的选择,因此GRAPPA算法具有一定的随意性。如何选择最佳的一个相邻子集以及这个集合的大小一直是GRAPPA算法的难点,一些方法通过迭代或者交叉验证的方法,去选择误差最小的相邻子集,但计算时间太长,不能满足实时性要求。同时,通过最小二乘法求解,只能使校准区域内的匹配误差最小,而没考虑联合权重函数的复杂度,容易出现过匹配,造成对未采集点的预测误差增大。
因此,如何提供一种具有最小全局误差,使得重建伪影较少的重建算法,同时计算速度较快的重建算法已成为业界急需解决的问题。
发明内容
本发明的目的在于提供一种多支撑向量机模型的磁共振并行成像方法,其具有最小全局误差,使得重建伪影较少、计算速度快。
为了实现上述目的,本发明的技术方案为:一种多支撑向量机模型的磁共振并行成像方法,其特征在于,包括如下步骤:
(1)用多通道线圈对K空间中间区域进行全采样后,将其划分为训练集和检验集,其它区域加速采样后作为预测集,并对各集合内的数据进行归一化处理;
(2)将训练集划分为多组训练子集,利用支撑向量机,选择不同的核函数和拟合参数(C,v)对各训练子集进行训练,得到不同的联合权重函数模型;
(3)利用检验集,对不同的联合权重函数进行检验,选择最佳的几个模型,
(4)使用最佳的几个联合权重函数对预测集进行预测,取其平均值作为未采集点的值,反归一化处理后,将K空间数据转换为图像。
进一步地,所述训练子集为Si={(x,y)i,(x,y)∈T,i=1,...,Nt},其是由训练集T={(x,y)}等分或者随机的划分为Nt个子集,其中y为全采样中一些要拟合的值,x为所有线圈在坐标上相应的邻域内采集到的值,将所有的y和对应的x列在一起组成训练集,
所述检验集为A={(x,y),(x,y)∈T},其是从全采样的数据中随机选取一部分,
所述预测集为P={(x,y0)},其中y0为未采集到的点,设为0。
进一步地,在步骤(1)中,对数据进行归一化处理包括:将K空间数据映射到[-1,1]上,其采取的公式为:
其中x为参考点的原始值,y为目标点的原始值,x′为归一化后参考点的值,y′为目标点归一化后的值。
进一步地,步骤(2)中,选择一种改进的支撑向量机(v-SVR),对不同的训练子集选择不同的核函数和拟合参数(C,v)进行训练,解其约束优化问题:
subjectto((w·xi)+b)-yi≤ε+ξi
其中,设求解的线性函数为y=wx+b,x为训练集中参考点的原始值,w为各参考点的联合权重值,b为线性函数在Y轴上的截距,y为训练集中目标点的原始值,C是平衡因子,调节目标函数的平滑度和误差,ε为误差精度控制项,可以被另一个变量v自动的调节,v∈(0,1),代表误差点所占比例的上限,l为训练样本总数,ξ和ξ*为松驰变量,通过引入核函数并解此优化问题的对偶形式:
subject to
αi, i=1,...l
进一步地,对于不同的训练子集,核函数分别选取线性核、径向基函数核和正则化的傅里叶核,线性核的形式为:k(xi,x)=xi Hx,
径向基函数核形式为:
傅里叶核形式为: 0<q<1,
其中γ和q分别为高斯核函数和傅里叶核函数的宽度参数,控制函数的径向作用范围。
进一步地,对于不同的训练子集,选取不同的平衡因子C和变量v,平衡因子C的范围从0.005到500,其变化步长为10,变量v从0.1到0.9,其变化步长为0.1。
进一步地,在步骤(3)中将步骤(2)中训练得到的各个子集的权重函数应用到检验集,对于每一个权重函数模型,其在检验集上的均方根误差定义为:
进一步地,在步骤(4)中利用步骤(3)选出的最佳模型分别对预测集内的点进行预测,最后通过线性联合平均其中为每一个子模型预测结果,然后对数据进行反归一化处理;并将K空间数据经过二维的快速离散傅里叶变换转换为图像,最后,通过平方和的方法将各个线圈的图像联合成一幅最终图像。
进一步地,步骤(4)中的数据反归一化公式为: 且将数据从[-1,1]映射到原始空间。
本发明与现有技术相比具有如下优点:选择一种改进的支撑向量机(v-SVR),因引入变量v,能自动调节不敏感损失函数中参数ε的大小,使用起来比支撑向量回归算法得方便和准确;利用支撑向量机求解得到的线圈联合权重函数,具有良好的泛化能力,对K空间噪声具有良好的抑制能力,同时,可以选择非线性的训练核函数进行非线性的重建,使得当采样加速因子过大时能较好去地除去混叠伪影,并用多个模型同时拟合,减少运算量,选择最佳的几个模型,使重建误差和伪影更少。
附图说明
图1为并行成像图像重建算法的流程图;
图2为MRI加速扫描数据填充方式和GRAPPA重建的示意图。
具体实施方式
下面结合具体的实施例及附图对本发明进行详细描述
如图1所示,一种多支撑向量机模型的磁共振并行成像方法,包括以下步骤:
对于自校准类的并行成像算法,勿需另外单独扫描得到线圈敏感度,只需在扫描时进行混合采样,对图2所示,其在K空间中,一部分相位编码按奈奎斯特采样速率进行采样,其它部分则进行加速采样,对于R倍的加速,则每采集一条相位编码后,隔R-1条相位编码步后再采集一条相位编码线。
对混合采样后的K空间数据划分为不同的集合,把全采样的数据作为训练集,其中一些点作为要拟合的值即y,而其所有线圈在坐标上相对应的邻域内采集到的值作为x(如图2),将所有的y和对应的x列在一起组成训练集T={(x,y)},把训练集等分或者随机的划分为Nt个子集Si={(x,y)i,(x,y)∈T,i=1,...,Nt},同时从全采样的数据中随机选取一部分作为检验集A={(x,y),(x,y)∈T},而把全部未采集到的数据作为预测集。
为了加快训练速度,需对训练集,检验集和预测集分别进行归一化处理,其归一化采取公式:
利用v-SVR,对不同的训练子集选择不同的核和参数进行训练,求解凸二次优化问题:
αi, i=1,...,l
其中,αi,为拉格朗日变量,C是平衡因子,调节目标函数的平滑度和误差,ε为误差精度控制项,可以被另一个变量v自动的调节,v∈(0,1),代表误差点所占比例的上限,l为训练样本总数。
核函数k(xi,xj)的选取不仅可以选择线性的核,还可以选取正则化的傅里叶核,以及径向基函数核。在本实例中,对于不同的训练子集,分别选取线性核,径向基函数核和正则化的傅里叶核,其形式为:
线性核:k(xi,x)=xi Hx;
径向基函数核:k(xi,x)=exp(-||x-xi||2/γ2);
傅里叶核: 0<q<1。
其中γ和q分别为高斯核函数和傅里叶核函数的宽度参数,控制函数的径向作用范围。在优选的实施方式中,γ和q分别设置为0.05和0.3。
在优选的实施方式中,对于不同的训练子集,还有选取不同的参数C和v,C的范围从0.005到500,其变化步长可以10,v从0.1到0.9。其变化步长可为0.1。
通过上述核函数拟合得到的联合权重函数为非线性的,还可以选择如多项式,B样条等其它非线性核。
对于不同的子模型可以得到不同的f(x),然后将这些函数模型对检验集上的数据进行预测,得到
将上述训练得到的各个联合权重函数模型应用到检验集,对于每一个权重函数模型,其在检验集上的均方根误差定义为:n为检验集的大小,为通过权重函数计算得到的估计值,y是检验集里本身采集到的值。最后选取均方根误差最小的Nm个子模型,经验显示,一般选择最佳子模型数达到10到15个时,其重建误差减小的趋势随着子模型数的增加将变得缓慢。
并将K空间数据经过二维的快速离散傅里叶变换转换为图像。最后,通过平方和的方法将各个线圈的图像联合成一幅最终图像。
Claims (10)
1.一种多支撑向量机模型的磁共振并行成像方法,其特征在于,包括如下步骤:
(1)用多通道线圈对K空间中间区域进行全采样后,将其划分为训练集和检验集,其它区域加速采样后作为预测集,并对各集合内的数据进行归一化处理;
(2)将训练集划分为多组训练子集,利用支撑向量机,选择不同的核函数和拟合参数对各训练子集进行训练,得到不同的联合权重函数模型;
(3)在检验集上,对各联合权重函数模型进行检验,选择最佳的几个子模型;
(4)分别使用最佳的几个子模型对预测集进行预测,取其平均值作为未采集点的值,反归一化处理后,将K空间数据转换为图像。
2.根据权利要求1所述的一种多支撑向量机模型的磁共振并行成像方法,其特征在于:所述训练子集为Si={(x,y)i,(x,y)∈T,i=1,...,Nt},其是由训练集T={(x,y)}等分或者随机的划分为Nt个子集,其中y为全采样中一些要拟合的值,x为所有线圈在坐标上相应的邻域内采集到的值,将所有的y和对应的x列在一起组成训练集,
所述检验集为A={(x,y),(x,y)∈T},其是从全采样的数据中随机选取一部分,
所述预测集为P={(x,y0)},其中y0为未采集到的点,设为0。
3.根据权利要求1所述的一种多支撑向量机模型的磁共振并行成像方法,其特征在于,在步骤(1)中,对数据进行归一化处理包括:将K空间数据映射到[-1,1]上,其采取的公式为:
其中x为参考点的原始值,y为目标点的原始值,x′为归一化后参考点的值,y′为目标点归一化后的值。
4.根据权利要求1或3所述的一种多支撑向量机模型的磁共振并行成像方法,其特征在于:步骤(2)中,选择一种改进的支撑向量机v-SVR,其目标优化问题为:
subjectto((w·xi)+b)-yi≤ε+ξi
其中,设求解的线性函数为y=wx+b,x为训练集中参考点的原始值,w为各参考点的联合权重值,b为线性函数在Y轴上的截距,y为训练集中目标点的原始值,C是平衡因子,调节目标函数y的平滑度和误差,ε为误差精度控制项,可以被另一个变量v自动的调节,v∈(0,1),代表误差点所占比例的上限,l为训练样本总数,ξ和ξ*为松驰变量,通过引入核函数并解此优化问题的对偶形式:
αi, i=1,...l
5.根据权利要求4所述的一种多支撑向量机模型的磁共振并行成像方法,其特征在于:对于不同的训练子集,核函数分别选取线性核、径向基函数核和正则化的傅里叶核,线性核的形式为:k(xi,x)=xi Hx,
径向基函数核形式为:
傅里叶核形式为: 0<q<1,
其中γ和q分别为高斯核函数和傅里叶核函数的宽度参数,控制函数的径向作用范围。
7.根据权利要求4所述的一种多支撑向量机模型的磁共振并行成像方法,其特征在于:对于不同的训练子集,选取不同的平衡因子C和变量v,平衡因子C的范围从0.005到500,其变化步长为10,变量v从0.1到0.9,其变化步长为0.1。
8.根据权利要求1所述的一种多支撑向量机模型的磁共振并行成像方法,其特征在于:在步骤(3)中将步骤(2)中训练得到的各个子集的权重函数应用到检验集,对于每一个权重函数模型,其在检验集上的均方根误差定义为:
其中,n为检验集的大小,为通过权重函数计算得到的估计值,y是检验集里本身采集到的值,最后选取均方根误差最小的Nm个子模型。
10.根据权利要求9所述的一种多支撑向量机模型的磁共振并行成像方法,其特征在于:步骤(4)中的数据反归一化公式为: 且将数据从[-1,1]映射到原始空间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210012061.0A CN102551723B (zh) | 2012-01-16 | 2012-01-16 | 一种多支撑向量机模型的磁共振并行成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210012061.0A CN102551723B (zh) | 2012-01-16 | 2012-01-16 | 一种多支撑向量机模型的磁共振并行成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102551723A true CN102551723A (zh) | 2012-07-11 |
CN102551723B CN102551723B (zh) | 2014-01-15 |
Family
ID=46399298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210012061.0A Expired - Fee Related CN102551723B (zh) | 2012-01-16 | 2012-01-16 | 一种多支撑向量机模型的磁共振并行成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102551723B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930567A (zh) * | 2012-09-25 | 2013-02-13 | 电子科技大学 | 多核加权最小二乘支撑向量机的磁共振并行成像重建方法 |
CN103064047A (zh) * | 2012-12-25 | 2013-04-24 | 深圳先进技术研究院 | 磁共振弹性成像的弹性重建修正方法和系统 |
CN106842089A (zh) * | 2017-01-23 | 2017-06-13 | 上海联影医疗科技有限公司 | 一种磁共振成像方法及系统 |
CN109069059A (zh) * | 2016-04-22 | 2018-12-21 | 通用电气公司 | 用于对移动的主体成像的系统和方法 |
CN109155012A (zh) * | 2016-12-30 | 2019-01-04 | 谷歌有限责任公司 | 评估机器学习模型的准确度 |
CN110940943A (zh) * | 2019-12-06 | 2020-03-31 | 上海联影医疗科技有限公司 | 搏动伪影校正模型的训练方法和搏动伪影校正方法 |
CN111175681A (zh) * | 2018-11-13 | 2020-05-19 | 西门子(深圳)磁共振有限公司 | 基于刀锋序列的磁共振成像方法、装置及其存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1557792A2 (en) * | 2004-01-19 | 2005-07-27 | Konica Minolta Medical & Graphic, Inc. | Medical image processing apparatus |
CN101308202A (zh) * | 2007-05-17 | 2008-11-19 | 西门子(中国)有限公司 | 并行采集图像重建的方法和装置 |
CN102043137A (zh) * | 2009-10-21 | 2011-05-04 | 凯斯西储大学 | 全程放射grappa校准 |
US20110257505A1 (en) * | 2010-04-20 | 2011-10-20 | Suri Jasjit S | Atheromatic?: imaging based symptomatic classification and cardiovascular stroke index estimation |
-
2012
- 2012-01-16 CN CN201210012061.0A patent/CN102551723B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1557792A2 (en) * | 2004-01-19 | 2005-07-27 | Konica Minolta Medical & Graphic, Inc. | Medical image processing apparatus |
CN101308202A (zh) * | 2007-05-17 | 2008-11-19 | 西门子(中国)有限公司 | 并行采集图像重建的方法和装置 |
CN102043137A (zh) * | 2009-10-21 | 2011-05-04 | 凯斯西储大学 | 全程放射grappa校准 |
US20110257505A1 (en) * | 2010-04-20 | 2011-10-20 | Suri Jasjit S | Atheromatic?: imaging based symptomatic classification and cardiovascular stroke index estimation |
Non-Patent Citations (4)
Title |
---|
D.A. KARRAS: "Efficient MRI Reconstruction Using a Hybrid Framework for Integrating Stepwise Bayesian Restoration and Neural Network Models in a Memory Based Priors System", 《HAIS》, 31 December 2008 (2008-12-31) * |
E. G. KHOLMOVSKI等: "Spatially Variant GRAPPA", 《PROC. INTL.SOC.MAG.RESON.MED.》, 31 December 2006 (2006-12-31) * |
LESLIE YING等: "Parallel MRI Using Phased Array Coils:Multichannel Sampling Theory Meeting Spin Physics", 《IEEE SIGNAL PROCESSING MAGAZINE》, 31 December 2010 (2010-12-31) * |
黎维娟等: "基于支持向量机的脑部MR图像细分类", 《电路与系统学报》, vol. 15, no. 1, 31 December 2010 (2010-12-31) * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930567A (zh) * | 2012-09-25 | 2013-02-13 | 电子科技大学 | 多核加权最小二乘支撑向量机的磁共振并行成像重建方法 |
CN102930567B (zh) * | 2012-09-25 | 2016-01-20 | 电子科技大学 | 多核加权最小二乘支撑向量机的磁共振并行成像重建方法 |
CN103064047A (zh) * | 2012-12-25 | 2013-04-24 | 深圳先进技术研究院 | 磁共振弹性成像的弹性重建修正方法和系统 |
CN103064047B (zh) * | 2012-12-25 | 2014-12-10 | 深圳先进技术研究院 | 磁共振弹性成像的弹性重建修正方法和系统 |
CN109069059A (zh) * | 2016-04-22 | 2018-12-21 | 通用电气公司 | 用于对移动的主体成像的系统和方法 |
CN109155012A (zh) * | 2016-12-30 | 2019-01-04 | 谷歌有限责任公司 | 评估机器学习模型的准确度 |
CN106842089A (zh) * | 2017-01-23 | 2017-06-13 | 上海联影医疗科技有限公司 | 一种磁共振成像方法及系统 |
CN111175681A (zh) * | 2018-11-13 | 2020-05-19 | 西门子(深圳)磁共振有限公司 | 基于刀锋序列的磁共振成像方法、装置及其存储介质 |
CN111175681B (zh) * | 2018-11-13 | 2022-08-30 | 西门子(深圳)磁共振有限公司 | 基于刀锋序列的磁共振成像方法、装置及其存储介质 |
US11474181B2 (en) * | 2018-11-13 | 2022-10-18 | Siemens Healthcare Gmbh | MRI method and device based on a blade sequence, and storage medium |
CN110940943A (zh) * | 2019-12-06 | 2020-03-31 | 上海联影医疗科技有限公司 | 搏动伪影校正模型的训练方法和搏动伪影校正方法 |
CN110940943B (zh) * | 2019-12-06 | 2021-11-02 | 上海联影医疗科技股份有限公司 | 搏动伪影校正模型的训练方法和搏动伪影校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102551723B (zh) | 2014-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102551723B (zh) | 一种多支撑向量机模型的磁共振并行成像方法 | |
Hoppe et al. | Deep learning for magnetic resonance fingerprinting: a new approach for predicting quantitative parameter values from time series | |
Sandino et al. | Compressed sensing: From research to clinical practice with deep neural networks: Shortening scan times for magnetic resonance imaging | |
Vishnevskiy et al. | Deep variational network for rapid 4D flow MRI reconstruction | |
CN104933683B (zh) | 一种用于磁共振快速成像的非凸低秩重建方法 | |
US9910118B2 (en) | Systems and methods for cartesian dynamic imaging | |
US10996306B2 (en) | MRI system and method using neural network for detection of patient motion | |
US10852376B2 (en) | Magnetic resonance imaging method and device | |
CN102930567B (zh) | 多核加权最小二乘支撑向量机的磁共振并行成像重建方法 | |
US9430854B2 (en) | System and method for model consistency constrained medical image reconstruction | |
CN112150568A (zh) | 基于Transformer模型的磁共振指纹成像重建方法 | |
CN103278784B (zh) | 一种多约束滑动窗的磁共振并行成像方法 | |
US20150276908A1 (en) | Magnetic resonance imaging apparatus and image processing apparatus | |
Korkmaz et al. | Deep MRI reconstruction with generative vision transformers | |
CN108447102A (zh) | 一种低秩与稀疏矩阵分解的动态磁共振成像方法 | |
CN103027681A (zh) | 用于重构并行获取的mri图像的系统 | |
CN103876739B (zh) | 用于并行成像应用的多阶段磁共振重建 | |
CN104778341B (zh) | 磁共振线圈合并系数计算方法、磁共振成像方法及其装置 | |
CN113628296A (zh) | 一种从时频域信号到二维图像的磁粒子成像重建方法 | |
CN110782452A (zh) | 一种t2定量图像成像方法及系统 | |
CN103519816B (zh) | 脑功能磁共振成像方法和系统 | |
CN106991651A (zh) | 基于合成分析反卷积网络的快速成像方法及系统 | |
CN112617798B (zh) | 一种基于Lp范数联合全变分的并行磁共振成像重构方法 | |
Lei et al. | Joint learning of multiple longitudinal prediction models by exploring internal relations | |
CN110123324A (zh) | 一种婴儿大脑t1加权磁共振成像优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140115 |
|
CF01 | Termination of patent right due to non-payment of annual fee |