CN104939828B - 磁共振成像方法 - Google Patents
磁共振成像方法 Download PDFInfo
- Publication number
- CN104939828B CN104939828B CN201510195370.XA CN201510195370A CN104939828B CN 104939828 B CN104939828 B CN 104939828B CN 201510195370 A CN201510195370 A CN 201510195370A CN 104939828 B CN104939828 B CN 104939828B
- Authority
- CN
- China
- Prior art keywords
- space
- plane
- phase
- subset
- referred
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/482—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
- G01R33/4822—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory in three dimensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
- G01R33/4826—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory in three dimensions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及一种磁共振成像方法。一种用于加速磁共振成像的方法,包括:在3D MRI中,将相位编码平面中的k空间划分为两个对称部分和三个不对称部分。在不同的部分中应用不同的采样密度。当随机采样应用在每个部分中时,通过对代价函数迭代地进行最小化来重建图像。通过利用k空间的共轭对称性,将相位约束项加入到该代价函数中以改进重建的质量。
Description
技术领域
本发明涉及用于三维数据集采集的磁共振成像方法,其中执行根据三个相互正交的磁场梯度的空间编码使得在一个空间方向k1上在读取梯度下读出信号,而通过在信号采集之前在另外两个空间方向k2和k3上应用相位编码梯度来执行所述另外两个空间方向上的空间编码,并且以顺序方式执行数据采集使得在每个采集步骤中在所述读出梯度下但是利用两个相位编码梯度的不同组合来采集信号。
这种类型的方法被称为是在参考文献[0]M.Bomans,K.Hohne,G.Laub,A.Pommert,U.Tiede,Improvent of 3D acquisition and visualization inMRI.Magn.Reson.Imaging 9,597-609(1991)中描述的3DFT(或体)成像。
背景技术
本发明的一般背景
本发明一般涉及磁共振成像(=MRI)技术。它具体涉及数据采集和图像重建方法以及用于MRI的空间编码。
与计算机断层扫描(=CT)相比,磁共振成像是相对新的技术,而第一张MR图像由P.C.Lauterbur于1973年发表在“Image Formation by Induced Local Interactions:Examples of Employing Nuclear Magnetic Resonance”,Nature 242,190491中。它最初是医学成像技术,最常用在放射学中以使身体的结构和功能可视化。它可以提供身体在任何平面中的细节图像。与CT相比,MRI提供身体的不同软组织之间的大得多的对比度,使MRI在神经学、心血管和肿瘤学成像中尤其有用。它使用强大的磁场来排列身体中水的氢原子的核磁化。射频场用来系统地改变这种磁化的排列,从而使氢原子核产生可被扫描仪检测的旋转磁场。该信号可以由另外的磁场操作以建立足够的信息来重建身体的图像。
通常,MRI系统一般沿着经受MRI过程的对象的中心轴建立均匀磁场。该均匀主磁场通过排列形成身体组织的原子和分子中的核自旋来影响要成像的对象的磁性能。如果核自旋的取向混乱不齐,则原子核试图通过场来重新排列它们的自旋。核自旋取向的混乱通常由调谐至感兴趣材料的拉莫尔频率的射频(RF)脉冲的施加引起。在重新排列处理期间,原子核围绕主磁场的方向进动(precess)并且发射电磁信号,这些电磁信号可以由放置在对象上或在对象周围的一个或多个RF检测器线圈检测。
磁共振成像采用时间上和空间上可变的磁场,以通过影响自旋的局部拉莫尔频率来对位置编码。通常用于该目的的梯度线圈产生叠加在主磁场上的空间编码磁场(=SEM)。这使得可以选择图像切片的局部化,并且还可以提供相位编码和频率编码。这种编码允许在图像重建期间识别共振信号的来源。图像质量和分辨率显著地依赖于所施加的编码场的强度和可被控制的方式。一般根据预建立的协议或事件的顺序(被称为脉冲顺序)来执行梯度线圈的控制,从而使不同类型的对比度机制得以成像。
具体背景技术
在半傅里叶成像中,留下k空间数据的一部分不采集以便缩短MRI中的测量时间,而缺少的数据通过在重建期间利用k空间的共轭对称性来合成(也称为部分傅里叶成像)。矩形采样图案常常用于多维半傅里叶成像(参见参考文献[1]、[2])。
偏共振效应和技术缺陷可能导致所采集的MR图像中的空间相位变化,这破坏k空间的共轭对称性。开发了两种类型的相位校正方法来合成未采集数据(参见参考文献[3]、[4])。在这些方法中,在数据合成期间,中间结果中的相位图像完全由估计值替代。
近来开发的k空间随机采样技术(也被称为“压缩感测”,参见参考文献[5])与半傅里叶采集组合以加速MRI(参见参考文献[2]、[6])。在参考文献[2]中,使用参考文献[3]中所提出的零差检测(Homodyne detection)方法来在单独的步骤中重建半傅里叶成像中的缺少的数据。在参考文献[6]中,使用参考文献[4]中所提出的迭代POCS方法来在单独的步骤中重建半傅里叶成像中的缺少的数据。在所有以上的方法中,重建的图像相位直接由估计相位替代,估计相位在具有快速相位变化的一些区域中可能不准确。在利用“直接”相位替代的这些方法中,在具有不良相位估计的区域中存在显著的重建误差。
发明内容
发明目的
本发明提出了一种基本上克服以上所讨论的现有方法的一个或多个缺点的方式。
本发明的一个目的是提出一种用于3D MRI中的半傅里叶成像的数据采集方法,其中针对第一空间维度(被称为k1)中的空间编码在读出梯度下采集数据。通过使用两个相互正交的相位编码梯度(被称为k2和k3)来对第二空间维度和第三空间维度中的数据进行空间编码。在k2-k3平面中执行k空间欠采样(undersampling)。
所述k2-k3平面中的k空间包括五部分
-部分1:中心k空间中的对称采集区域。
-部分2:具有比部分1中更高的空间频率的对称采集区域。
-部分3:具有比部分2中更高的空间频率的下半k空间中的采集区域。
-部分4:具有比部分2中更高的空间频率的上半k空间中的采集区域。
-部分5:具有比部分4中更高的空间频率的上半k空间中的未采集区域。
对所述k空间中的部分1进行完全采样。
用均匀采样密度对所述k空间中的部分2进行欠采样。
用比部分2中更低的均匀采样密度对所述k空间中的部分3进行欠采样。
用比部分3中更低的均匀采样密度对所述k空间中的部分4进行欠采样。
完全不采集所述k空间中的部分5。
所述部分2和部分4的形状的一种选择是椭圆。
所述部分2和部分4的形状的一种选择是矩形。
数据采集的一种选择是在所述部分2、部分3和部分4中应用矩形欠采样(参见参考文献[7])。
对数据采集的一种选择是在所述部分2、部分3和部分4中应用随机欠采样。
本发明的另一目的是提出一种用于通过所述k空间采样图案来采集的数据集的图像重建方法。
发明概要
这些目的通过对以上讨论的方法进行如下修改来实现:对被称为部分1的k2-k3平面的第一子集进行完全采样,部分1关于k空间中心对称;用均匀密度对被称为部分2的k2-k3平面的第二子集进行欠采样,部分2具有比部分1中更高的空间频率且关于k空间中心对称;用均匀密度对被称为部分3的下半k2-k3平面中的第三子集进行欠采样,部分3具有比部分2中更高的空间频率;用均匀密度对被称为部分4的上半k2-k3平面中的第四子集进行欠采样,部分4具有比部分2中更高的空间频率;完全不采集被称为部分5的上半k2-k3平面中的第五子集,部分5具有比部分4中更高的空间频率;
部分2中的采样密度比部分1中低,部分3中的采样密度比部分2中低,部分4中的采样密度比部分3中低,部分5中不采集数据;以及
通过像在非线性共轭梯度方法中那样用下降算法对代价函数迭代地进行最小化来重建图像,该代价函数是多个正则项(regularization term)的加权和,相位约束项被引入该代价函数,被称为Rpc的该相位约束项为:其中x表示迭代最小化中的中间解,||g||1=∑k|gk|,gk是矩阵g的第k个元素;AοB表示矩阵A和矩阵B的Hadamard乘积;|x|表示x的幅值;PR表示要被重建的图像的相位的估计值;W表示权重图(weighting map)。
本发明提出用于以上所讨论的图像重建的下列步骤:
通过仅使用所述部分1和部分2中的采集数据来重建图像(被称为IR)。在所述IR的重建中,将所述部分3、部分4和部分5中的数据设置为零。提取所述IR中的相位图像(被称为PR)。
通过使用下降算法对代价函数fc迭代地进行最小化来获得最终图像(被称为IF),下降算法例如参考文献[5]本身提出的非线性共轭梯度方法。
所述代价函数fc可以表示为下列正则项的加权和,但不限于列出的项:
x表示所述迭代最小化处理中IF的中间解;y是所采集的k空间数据;F是傅里叶变换运算符;U表示具有二进制元(binary entry)的k空间采样掩码;Ψ是稀疏变换运算符(例如小波变换);TV是总变化运算符;Rpc是相位约束项。λl1、λTV是两个实值非负比例因子;λpc是实值正比例因子。
所述代价函数fc中的Rpc表示为:
所述Rpc中的||g||1可以表示为:||g||1=∑k|gk|。其中是gk是矩阵g的第k个元素。
所述Rpc中的AοB表示矩阵A和矩阵B的Hadamard乘积。
所述Rpc中的|x|表示x的幅值。
所述Rpc中的W表示权重图。
一种选择是:所述W=|x|。所述W在迭代最小化的处理中变化。为了找到代价函数fc的所述迭代最小化中的下降方向,计算所述相位约束项Rpc的梯度:
Rpc的梯度为:
其中
其中
是具有元的对角矩阵;
是具有元的对角矩阵。
ε是正的小参数以避免“被零除”。
表示x的共轭。
(A)*表示矩阵A的共轭转置。
一种选择是:所述W在迭代最小化处理中是恒定的。W与通过将λpc设置为零而使用全部或部分采集的k空间来重建的幅值图像成比例。为了找到代价函数fc的所述迭代最小化中的下降方向,计算所述相位约束项Rpc的梯度:
Rpc的梯度为:
其中
其中
是具有元的对角矩阵;
是具有元的对角矩阵。
ε是正的小参数以避免“被零除”。
表示x的共轭。
(A)*表示矩阵A的共轭转置。
所述重建方法优选应用到通过所述数据采集方法采集的数据集,但不限于所述数据采集方法采集的数据集。
所述重建方法不限于3D MRI。
附图说明
通过结合附图来仔细研究本发明目前优选的示例性实施例的下面的详细说明,可以更好地理解和明白本发明的这些以及其它目的和优点。
为了使本发明的前述和其它特征和优点对本领域技术人员更清楚,下面将参考附图来详细描述本发明的优选实施例,附图中相同的标号表示相同的部件。
图1示出k空间中的采样密度的分布的示图:a)部分2和部分4是椭圆形;b)部分2和部分4是矩形;
图2示出通过网格k空间中的随机采样从图1a)中的采样密度的分布产生的采样图案的示图;
图3示出不同重建方法的比较:
a)参考图像;
b)采用本发明方法的重建误差图像;
c)采用参考文献[4]中的POCS方法的重建误差图像;
d)采用参考文献[3]中的零差检测方法的重建误差图像。
具体实施方式
在本发明的优选实施例中,k2-k3平面的子集的形状是矩形。这与k空间常规欠采样兼容,并且通过使用参考文献[7]中的算法有效地重建缺少的数据。
在本发明的替代实施例中,k2-k3平面的子集的形状是椭圆形。这里,通过有限量的采集的k空间数据来有效地估计图像相位。
在本发明的另一实施例中,其中通过仅使用在k2-k3平面中的部分1和部分2中对称采集的数据来从重建的图像获得相位PR的估计值。根据本发明的方法的这种变形简单且直接,并且在具有慢相位变化的区域中提供充分的准确性。
在本发明的另一类实施例中,权重图W设置为中间解x的幅值。图像相位估计的准确性依赖于信噪比(=SNR)。SNR与图像像素的强度成正比。通过应用权重图,对相位约束的依赖性在具有不良相位估计的区域中减小以减小重建误差。通过将中间解x的幅值用作权重图,当迭代重建收敛时权重图更精确地与SNR成比例,这能够充分地抑制在具有不良SNR的区域中的假象(artifact)。
在本发明的替代类实施例中,权重图W设置为不使用Rpc的重建的图像的幅值。在当权重图W为不使用Rpc重建的恒定幅值图像时的所有情况下Rpc为凸项。下降算法变得更鲁降以利用凸正则项找到最优解。
本发明公开了一种用于加速磁共振成像的方法,包括:在3D MRI中,将相位编码平面中的k空间划分为两个对称部分和三个不对称部分。在不同的部分中应用不同的采样密度。当随机采样应用在每个部分中时,通过对代价函数迭代地进行最小化来重建图像。通过利用k空间的共轭对称性,将相位约束项加入到该代价函数中以改进重建的质量。
虽然本文中仅图示和描述了本发明的某些特征,但是本领域技术人员将想到很多变型和更改。因此,应当理解,所附权利要求书意图覆盖落入本发明真正精神的所有这些变型和更改。
参考文献
[0]M.Bomans et al,Magn.Reson.Imaging 9,597-609(1991);
[1]Li F.et al,ISMRM 2011:4368;
[2]Liu F.et al 2012 Phys.Med.Biol.57 N391;
[3]Nell,D.C.,IEEE Trans.Med.Imaging,10,154-163,1991;
[4]EM Hacke,J Magn Reson 1991;92:126-145;
[5]Lustig M et al,MRM 2007,1182-1195;
[6]Doneva M et al,ISMRM 2010:485;
[7]Griswold,M.A.,et al,Magn.Reson.Med.,47:1202-1210.
Claims (8)
1.一种用于采集三维数据集的磁共振成像方法,其中执行根据三个相互正交的磁场梯度的空间编码使得在一个空间方向k1上在读取梯度下读出信号,而通过在信号采集之前在另外两个空间方向k2、k3上应用相位编码梯度来执行所述另外两个空间方向上的空间编码,并且以顺序方式执行数据采集使得在每个采集步骤中在所述读取梯度下但是利用两个相位编码梯度的不同组合来采集信号,其特征在于,
对被称为部分1的k2-k3平面的第一子集进行完全采样,部分1关于k空间中心对称;用均匀密度对被称为部分2的k2-k3平面的第二子集进行欠采样,部分2具有比部分1中更高的空间频率且关于k空间中心对称;用均匀密度对被称为部分3的下半k2-k3平面中的第三子集进行欠采样,部分3具有比部分2中更高的空间频率;用均匀密度对被称为部分4的上半k2-k3平面中的第四子集进行欠采样,部分4具有比部分2中更高的空间频率;完全不采集被称为部分5的上半k2-k3平面中的第五子集,部分5具有比部分4中更高的空间频率;
部分2中的采样密度比部分1中低,部分3中的采样密度比部分2中低,部分4中的采样密度比部分3中低,部分5中不采集数据;以及
通过用下降算法对代价函数迭代地进行最小化来重建图像,该代价函数是多个正则项的加权和,相位约束项被引入该代价函数,被称为Rpc的该相位约束项为:
其中x表示迭代最小化中的中间解;||g||1=∑k|gk|,gk是矩阵g的第k个元素;表示矩阵A和矩阵B的Hadamard乘积;|x|表示x的幅值;PR表示要重建的图像的相位的估计值;W表示权重图。
2.根据权利要求1所述的方法,其中所述下降算法包括非线性共轭梯度方法。
3.根据权利要求1所述的方法,其中k2-k3平面的子集的形状是矩形。
4.根据权利要求1所述的方法,其中k2-k3平面的子集的形状是椭圆形。
5.根据权利要求3所述的方法,其中通过仅使用在k2-k3平面中的部分1和部分2中对称采集的数据来从重建的图像获得相位PR的估计值。
6.根据权利要求4所述的方法,其中通过仅使用在k2-k3平面中的部分1和部分2中对称采集的数据来从重建的图像获得相位PR的估计值。
7.根据权利要求1至6中任一项所述的方法,其中权重图W设置为中间解x的幅值。
8.根据权利要求1至6中任一项所述的方法,其中权重图W设置为不使用Rpc的重建的图像的幅值。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP14162404.9 | 2014-03-28 | ||
EP14162404.9A EP2924457B1 (en) | 2014-03-28 | 2014-03-28 | Half Fourier MRI with iterative reconstruction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104939828A CN104939828A (zh) | 2015-09-30 |
CN104939828B true CN104939828B (zh) | 2018-10-12 |
Family
ID=50382375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510195370.XA Active CN104939828B (zh) | 2014-03-28 | 2015-03-27 | 磁共振成像方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9478051B2 (zh) |
EP (1) | EP2924457B1 (zh) |
CN (1) | CN104939828B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016183572A1 (en) | 2015-05-14 | 2016-11-17 | Ohio State Innovation Foundation | Systems and methods for estimating complex b1+ fields of transmit coils of a magnetic resonance imaging (mri) system |
WO2017206048A1 (en) * | 2016-05-31 | 2017-12-07 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for removing gibbs artifact in medical imaging system |
JP6636676B1 (ja) * | 2016-11-17 | 2020-01-29 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 強度補正された磁気共鳴画像 |
US10890631B2 (en) | 2017-01-19 | 2021-01-12 | Ohio State Innovation Foundation | Estimating absolute phase of radio frequency fields of transmit and receive coils in a magnetic resonance |
US10551458B2 (en) | 2017-06-29 | 2020-02-04 | General Electric Company | Method and systems for iteratively reconstructing multi-shot, multi-acquisition MRI data |
WO2019000348A1 (en) * | 2017-06-29 | 2019-01-03 | Shanghai United Imaging Healthcare Co., Ltd. | SYSTEM AND METHOD FOR ACCELERATING MAGNETIC RESONANCE IMAGING |
CN107622474B (zh) * | 2017-09-26 | 2021-03-30 | 北京大学深圳研究生院 | 基于主视点的全景视频映射方法 |
US20190180481A1 (en) * | 2017-12-13 | 2019-06-13 | General Electric Company | Tomographic reconstruction with weights |
CN108324276B (zh) * | 2018-01-11 | 2021-07-30 | 上海东软医疗科技有限公司 | 磁共振成像方法和装置 |
JP7164320B2 (ja) * | 2018-05-11 | 2022-11-01 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置、医用画像処理装置、及び画像再構成方法 |
EP3598162A1 (en) * | 2018-07-19 | 2020-01-22 | Siemens Healthcare GmbH | Magnetic resonance method and apparatus for suppressing metal artifacts |
US10950014B2 (en) * | 2019-03-01 | 2021-03-16 | Canon Medical Systems Corporation | Method and apparatus for adaptive compressed sensing (CS) to correct motion artifacts in magnetic resonance imaging (MRI) |
US11313933B2 (en) * | 2019-03-20 | 2022-04-26 | The Board Of Trustees Of The Leland Stanford Junior University | Diffusion-weighted MRI with magnitude-based locally low-rank regularization |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7545966B2 (en) * | 2003-05-05 | 2009-06-09 | Case Western Reserve University | Efficient methods for reconstruction and deblurring of magnetic resonance images |
CN100366218C (zh) * | 2003-05-09 | 2008-02-06 | 西门子(中国)有限公司 | 磁共振成像方法 |
WO2005003805A1 (en) * | 2003-06-27 | 2005-01-13 | Case Western University | Efficient method for mr image reconstruction using coil sensitivity encoding |
CN101308202B (zh) * | 2007-05-17 | 2011-04-06 | 西门子公司 | 并行采集图像重建的方法和装置 |
US8089278B1 (en) * | 2008-04-11 | 2012-01-03 | The Regents Of The University Of California | Time-resolved contrast-enhanced magnetic resonance (MR) angiography |
US7592808B1 (en) * | 2008-05-06 | 2009-09-22 | General Electric Company | System and method for reducing MR scan time using partial fourier acquisition and compressed sensing |
EP2317333A1 (de) * | 2009-10-30 | 2011-05-04 | Max-Delbrück-Centrum für Molekulare Medizin (MDC) | MRT-Betriebsverfahren |
CN102232830B (zh) * | 2010-04-30 | 2014-09-03 | 西门子(深圳)磁共振有限公司 | 一种磁共振成像水脂分离方法 |
CN101975936A (zh) * | 2010-09-03 | 2011-02-16 | 杭州电子科技大学 | 一种基于cs压缩感知技术的快速磁共振成像方法 |
CN101975935B (zh) * | 2010-09-03 | 2012-12-05 | 杭州电子科技大学 | 一种基于部分回波压缩感知的快速磁共振成像方法 |
DE102011005445B4 (de) * | 2011-03-11 | 2014-10-09 | Siemens Aktiengesellschaft | Normalisierung von Magnetresonanzbilddaten bei bewegtem Tisch |
US8886283B1 (en) * | 2011-06-21 | 2014-11-11 | Stc.Unm | 3D and 4D magnetic susceptibility tomography based on complex MR images |
-
2014
- 2014-03-28 EP EP14162404.9A patent/EP2924457B1/en active Active
-
2015
- 2015-03-11 US US14/644,221 patent/US9478051B2/en active Active
- 2015-03-27 CN CN201510195370.XA patent/CN104939828B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
US9478051B2 (en) | 2016-10-25 |
US20150279065A1 (en) | 2015-10-01 |
EP2924457A1 (en) | 2015-09-30 |
CN104939828A (zh) | 2015-09-30 |
EP2924457B1 (en) | 2016-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104939828B (zh) | 磁共振成像方法 | |
Wang et al. | Model‐based T 1 mapping with sparsity constraints using single‐shot inversion‐recovery radial FLASH | |
Huang et al. | T2 mapping from highly undersampled data by reconstruction of principal component coefficient maps using compressed sensing | |
JP6513398B2 (ja) | 事前情報に制約される正則化を使用するmr画像再構成 | |
US9983284B2 (en) | MRI with dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field | |
US9766313B2 (en) | MR imaging using apt contrast enhancement and sampling at multiple echo times | |
US10234523B2 (en) | MRI with dixon-type water/fat separation with estimation of the main magnetic field variations | |
US8379951B2 (en) | Auto calibration parallel imaging reconstruction method from arbitrary k-space sampling | |
JP5848713B2 (ja) | 磁気共鳴イメージング装置及びコントラスト強調画像取得方法 | |
CN106574954A (zh) | 针对epi的具有奈奎斯特伪影校正的并行mr成像 | |
US20090238430A1 (en) | Mri acquisition using sense and highly undersampled fourier space sampling | |
CN105814449B (zh) | 具有水/脂肪分离的零回波时间mr成像 | |
WO2014115043A1 (en) | Metal resistant mr imaging | |
Levine et al. | Accelerated three‐dimensional multispectral MRI with robust principal component analysis for separation of on‐and off‐resonance signals | |
US9274197B2 (en) | Magnetic resonance imaging data sampling methods and systems | |
US10324154B2 (en) | Generalized spherical deconvolution in diffusion magnetic resonance imaging | |
EP3931588A1 (en) | Epi mr imaging with distortion correction | |
US20160334487A1 (en) | Apparatus and method for improving balanced steady-state free precision in magnetic resonance imaging | |
FERNANDES | MULTICHANNEL DENOISING STRATEGIES FOR ULTRAFAST MAGNETIC RESONANCE IMAGING | |
CN116953582A (zh) | 磁共振成像设备的多个线圈灵敏度图和sense重建 | |
She | Signal processing techniques for improving image reconstruction of parallel Magnetic Resonance Imaging and dynamic Magnetic Resonance Imaging | |
Mayer et al. | Synthetic aperture MRI | |
Berman | Accelerated Radial Magnetic Resonance Imaging: New Applications and Methods | |
Makhijani | New methods for carotid MRI | |
Bones et al. | Sparse sampling in MRI |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |