CN107750338A - 使用星形叠层采集的mr成像 - Google Patents
使用星形叠层采集的mr成像 Download PDFInfo
- Publication number
- CN107750338A CN107750338A CN201680035043.0A CN201680035043A CN107750338A CN 107750338 A CN107750338 A CN 107750338A CN 201680035043 A CN201680035043 A CN 201680035043A CN 107750338 A CN107750338 A CN 107750338A
- Authority
- CN
- China
- Prior art keywords
- space
- signals
- radial
- images
- opening position
- 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
- 238000003475 lamination Methods 0.000 title claims abstract description 23
- 238000009826 distribution Methods 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000003384 imaging method Methods 0.000 claims abstract description 26
- 238000007689 inspection Methods 0.000 claims abstract description 6
- 230000002093 peripheral effect Effects 0.000 claims description 16
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 9
- 239000010931 gold Substances 0.000 claims description 9
- 229910052737 gold Inorganic materials 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 230000002829 reductive effect Effects 0.000 abstract description 6
- 238000005070 sampling Methods 0.000 description 41
- 230000005415 magnetization Effects 0.000 description 8
- 239000002872 contrast media Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000018199 S phase Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000005520 cutting process Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 230000001965 increasing effect Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000010030 laminating Methods 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- 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
- 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/5601—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
-
- 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/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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56308—Characterization of motion or flow; Dynamic imaging
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Vascular Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及对放置在MR设备(1)的检查体积中的至少一个对象(10)的MR成像的方法。本发明的一个目的是使用具有降低的条纹伪影水平的星形叠层采集方案来实现快速MR成像。本发明的方法包括:‑使所述对象(10)经受至少一个RF脉冲和切换的磁场梯度的成像序列,‑根据星形叠层方案采集MR信号,其中,所述MR信号从沿切片方向布置在不同位置处的多个平行切片(21‑27)被采集为径向k空间分布(S1‑S12),其中,所述k空间分布(S1‑S12)的径向密度根据所述切片位置而变化,其中,所述径向密度在更靠中心的k空间位置处更高,并且在更靠***的k空间位置处更低,并且其中,相比于从更靠***的k空间位置处的切片采集k空间分布,以更高的时间密度来从更靠中心的位置处的切片采集k空间分布,并且‑根据所述MR信号重建MR图像。
Description
技术领域
本发明涉及磁共振(MR)成像领域。本发明涉及一种对对象的MR成像的方法。本发明还涉及MR设备和在MR设备上运行的计算机程序。
背景技术
利用磁场与核自旋之间的相互作用以便形成二维或三维图像的图像形成MR方法现今被广泛使用,特别是在医学诊断领域,因为对于软组织的成像,它们在许多方面优于其他成像方法,不需要电离辐射并且通常不是有创的。
根据一般的MR方法,将例如待检查的患者的身体的对象布置在强而均匀的磁场中,该磁场的方向同时限定了测量所依据的坐标系的轴(通常为z轴)。根据磁场强度,磁场为个体核自旋产生不同的能级,其能够通过应用限定频率(所谓的拉莫尔频率,或MR频率)的电磁交变场(RF场)而被激励(自旋共振)。从宏观角度来看,个体核自旋的分布产生整体磁化,通过施加适当频率的电磁脉冲(RF脉冲)整体磁化能够偏转出平衡状态,从而磁化执行关于z轴的进动运动。进动运动描述了其孔径角被称为翻转角的锥体表面。翻转角的大小取决于施加的电磁脉冲的强度和持续时间。在所谓的90°脉冲的情况下,自旋从z轴偏转到横向平面(翻转角90°)。
在RF脉冲终止之后,磁化驰豫回到原始平衡状态,其中,在z方向上的磁化以第一时间常数T1(自旋晶格或纵向弛豫时间)再次建立,并且在垂直于z方向的方向上的磁化以第二时间常数T2(自旋-自旋或横向弛豫时间)弛豫。磁化的变化能够借助于在MR设备的检查体积内布置和取向的接收RF线圈以在垂直于z轴的方向上测量磁化的变化的方式来检测。在施加例如90°脉冲之后,横向磁化的衰减伴随着核自旋(由局部磁场不均匀性引起)从具有相同相位的有序状态过渡到在其中全部相位角均匀分布(移相)的状态的过渡。移相能够借助重新聚焦脉冲(例如180°脉冲)来补偿。这在接收线圈中产生回波信号(自旋回波)。
为了实现身体内的空间分辨率,沿着三个主轴延伸的恒定磁场梯度被叠加在均匀磁场上,导致自旋共振频率的线性空间依赖性。接收线圈中拾取的信号则包含不同频率的分量,这些分量能够与身体中的不同位置相关联。经由接收线圈获得的信号数据对应于空间频率域并被称为k空间数据。k空间数据通常包括用不同相位编码采集的多条线。每条线都通过收集多个样本进行数字化。借助图像重建算法将一组k空间数据转换成MR图像。
在已知的所谓的三维(3D)星形叠层(stack-of-stars)采集方案(参见例如WO2013/159044A1)中,应用了多个空间非选择性或板选择性RF激励,每个空间非选择性或板选择性RF激励之后是一个或多个MR信号(例如梯度回波信号)的采集,其中,每个MR信号表示k空间分布。MR信号被从多个平行切片采集为径向k空间分布。切片沿切片方向布置在不同的位置。在切片方向(例如kz方向)上执行标准笛卡尔相位编码,同时在沿着围绕中心(kx=ky=0)旋转的径向“辐条”的每个单个切片内采集MR信号。这导致由叠置的圆盘(“星形叠层”)组成的圆柱形k空间覆盖。从技术上讲,这是通过同时在切片的平面内方向上生成磁场梯度并调制其幅度来实现的。能够使用不同的方案来选择k空间分布采集步骤的时间顺序。例如,能够在采集不同角位置处的k空间分布之前顺序地采集沿切片方向的所有相位编码步骤。这确保了笛卡尔采样的周期保持较短,这导致在切片叠层中的高度的数据一致性,并且保持了对星形叠层方法的径向采样的一般运动鲁棒性。笛卡尔相位编码步骤可以从中心切片到k空间***(中心向外)执行,或者以从-kz,最大到+kz,最大的线性顺序执行。对于角排序,成像序列能够使用具有多个交错的等距角采样或者所谓的黄金角度方案。在等距方案中,角距离根据ΔΦ=180°/n总计算,其中,n总是辐条的总数。使用多个交错(或“旋转”)来采集辐条可能是有益的,因为交错减少了k空间中的时间相干性。因此,运动不一致性在k空间中分散并且伪影被衰减。在黄金角度方案中,k空间分布的角度每次递增ΔΦ=111.25°,这对应于180°乘以黄金比例。因此,随后采样的辐条总是增加补充信息,同时填补先前采样的辐条集中最大的间隙。因此,采集的辐条的任何顺序集大致均匀地覆盖了k空间,这实现了时间子帧的重建,并且使得黄金角度方案非常适合于动态成像研究。
上述的3D径向星形叠层方案为临床MR成像提供了几个有前景的优点,如高运动鲁棒性,良性混叠伪影和k空间中心的连续更新。然而,尽管有这些益处,条纹伪影往往是一个问题。结果是所需的径向采样密度更高,并且因此扫描时间更长。
D.C.Peters等人的ISMRM-摘要2007第207页的“Radial undersampling that isvariable in kz”公开了一种在k空间中的星形叠层采样方法。使用在kz维上变化的径向欠采样模式,采集用于中心kz空间的更多数量的投影。
发明内容
从上述内容容易理解,需要一种改进的MR成像技术。因此,本发明的目的是使用具有降低的条纹伪影水平的星形叠层采集方案来实现快速MR成像。
根据本发明,公开了一种对放置在MR设备的检查体积中的对象的MR成像的方法。该方法包括:
-使对象经受至少一个RF脉冲和切换的磁场梯度的成像序列,
-根据星形叠层方案采集MR信号,其中,所述MR信号在多个不同的角位置处从沿切片方向布置在不同位置处的多个平行切片被采集为径向k空间分布,其中,所述k空间分布的径向密度根据切片位置而变化,其中,所述径向密度在更靠中心的k空间位置处更高,并且在更靠***的k空间位置处更低,其中,相比于从更靠***的k空间位置处的切片采集k空间分布,以更高的时间密度(即,每单位时间采集的k空间分布的数量)来从更靠中心的k空间位置处的切片采集k空间分布,并且
-根据所述MR信号重建MR图像。
本发明提出应用修改的星形叠层采集策略。修改是径向k空间采样密度(即,每个切片的不同取向的径向k空间分布的数量)在切片之间变化。k空间分布的径向密度在更靠近k空间中心(kz=0)定位的切片中更高,而k空间分布的径向密度在更远离k空间中心定位的切片中更低。
k空间中心附近的径向k空间采样的增加的密度导致降低的条纹伪影水平,而总扫描时间能够保持在最小。
由于k空间的中心区域的规则更新,所以在对比增强的MR成像的背景下,径向k空间采样通常被已知是感兴趣的,因为其允许跟随感兴趣区域中的造影剂的摄取。尽管如此,即使与k空间加权图像对比方法(参见Song等人,Magnetic Resonance in Medicine,第44卷,2004年,第825-832页)结合以用于对比度操控,从常规径向的星形叠层采集重建的MR图像的时间序列显示模糊的对比度增强分布,使得造影剂到达特定图像体素的时间点不能被清楚地区分。这能够通过本发明的方法来改进。
在优选的实施例中,本发明提出了一种采集顺序,其中,相比于从更靠***的k空间位置处的切片采集径向k空间分布,以更高的时间密度来从更靠中心的k空间位置处的切片采集径向k空间分布。换句话说,中央奈奎斯特椭圆体内的径向k空间分布是在接近的时间点获得的,而***径向k空间分布(即,在奈奎斯特椭圆体之外)在更远离的时间点被采样。注意,随着来自k空间的更多采样被采集,在其内部采样的k空间分布的(径向)密度等于或者高于针对视场的奈奎斯特准则的中央奈奎斯特椭圆体的尺寸将在k空间中径向向外逐渐增加,而***k空间中采样的k空间分布的(径向)密度等于或小于针对视场的奈奎斯特准则。奈奎斯特椭圆体的实际尺寸确定了重建的MR图像的空间分辨率。
这种方法改进了对比增强的磁共振成像。中心k空间椭圆体处增加的时间密度确保造影剂到达特定图像体素的对比时间点被捕获。
然后使用奈奎斯特(Nyquist)椭圆体内的仅在各个最接近的时间点处采集的k空间分布与从***k空间区域采集的k空间分布(其在更远离的时间点处采样)相结合来重建每个MR图像。这意味着在中央完全采样或过采样的椭圆体k空间区域中选择时间上彼此最接近地采集的MR信号部分用于重建MR图像。对于***k空间MR椭圆体,仅选择k空间中的从向外的径向k空间位置延伸的部分径向辐条。随着这些部分辐条向外移动到***k空间MR椭圆体,这些部分辐条具有更高的角密度。
通过这种方法,能够更清楚地区分造影剂到达特定图像体素的时间点,因为仅在时间上紧密采集的MR信号被用于k空间中心的MR图像重建。因此降低所测量的对比增强分布的模糊。
在重建MR图像之前,可以移除在椭圆体之外的k空间位置处的MR信号。这构成了k空间中的椭圆“快门”。这有效地消除了条纹伪影。在其他有利的实施方式中,k空间的中心区域可以被过采样(以比乃奎斯特准则更高的采样密度),这使得MR信号的采集对于运动伪影更加鲁棒。
在用于重建MR图像的更精细的实施方式中,来自***k空间的MR信号沿着k空间中的从向外的径向k空间位置延伸的部分径向辐条进行选择,并且这些部分辐条在k空间中具有较高的角密度,因为它们向外的径向k空间位置在k空间中更为径向向外。部分辐条从它们向外的径向k空间位置延伸,向外的径向k空间位置在每个kz平面中与星形叠层的kz轴径向分开。这可以在各体kz平面中对中心区域进行完全采样或甚至过采样来实施,例如通过首先沿着在k空间中以黄金角度角偏移的(三个)径向辐条采样。在k空间中的一个或若干围绕中心区域的中间环形区域中,除了沿起源于kz轴并且还采样中心k空间区域的径向辐条的选择之外,还沿着从偏离kz轴的向外的径向位置延伸的若干部分径向辐条来选择用于重建的MR信号。在k空间的径向超出一个或多个中间环形区域的***,选择沿k空间中从超过中间环形区域的向外的径向位置延伸的另外的部分辐条被添加。在k空间中部分径向辐条的角密度沿径向向外增加。
对k空间的中心区域进行采样,使得以比在向外(在k空间轴向以及径向)区域中添加样本更高的时间率添加样本。在信号采集中,优选沿着从星形叠层的其切片的中心径向向外延伸的径向分布采集k空间分布。为了重建,只选择来自部分辐条的MR信号。必须指出的是,在已知的PROPELLER成像技术中采集的k空间分布也被认为是本发明含义内的径向k空间分布。在PROPELLER概念中,MR信号在N个条带中在k空间中被采集,每个条带包括L个并行的k空间分布,对应于基于笛卡尔的k空间采样方案中的L个最低频率相位编码线。也被称为K空间叶片的每个条带在K空间中旋转例如180°/N的角度,使得MR信号组的总集合在K空间中跨越圆形。在本发明的含义内,星形叠层因此可以同样是“叶片叠层(stack-of-blades)”。
根据本发明的优选实施例,k空间分布的径向密度发生变化,使得根据给定视场(FOV)的奈奎斯特准则在围绕k空间的中心的椭圆体内被满足。椭圆体的尺寸能够根据个体成像任务的采集速度和图像质量之间的所需折衷来选择。本发明使总扫描时间最小化,同时在图像能量主导的中央k空间中提供过采样(即,比奈奎斯特准则所要求的更高的径向采样密度)。径向采样密度可以逐渐减小到***k空间区域中的较低的径向采样密度。在***区域中,即在椭圆体外部,径向k空间采样密度甚至可以低于奈奎斯特阈值,而对图像质量没有显著的影响。如上所述,根据本发明的方法在奈奎斯特椭球内发生k空间的过采样。这允许仅选择在过采样的中央k空间区域中时间上彼此最接近地采集的MR信号,以重建每个MR图像而不损害图像质量。
根据本发明的优选实施例,连续采集的k空间分布的角度增量是黄金角度。这对应于上述的黄金角度方案,其中,k空间分布的角度每次递增ΔΦ=111.25°。这对于动态成像是特别有利的。
优选地,在采集不同角位置处的另外的k空间分布之前,在相同角位置处从至少两个不同切片采集k空间分布。换句话说,在采集不同黄金角度位置处的k空间分布之前,顺序地采集沿切片方向的相位编码步骤。这确保了笛卡尔采样的周期保持较短,这导致在切片叠层内的高度数据一致性并且保持了一般的运动鲁棒性。
根据本发明,以可变密度对k空间进行采样。这可以被采用以便通过(运动优化的)可变密度压缩感测来减少扫描时间。k空间分布甚至可以不规则地分布在***k空间区域中,即在k空间被完全采样的椭圆体之外。已知CS的理论对于显著的信号数据减少具有巨大的潜力。在CS理论中,通过应用合适的正则化算法,能够从欠采样测量恢复在变换域中具有稀疏表示的信号数据集。欠采样的可能性导致显著减少的采集时间。作为信号采样和重建的数学框架,CS规定了以下条件:在该条件下,信号数据集能够被精确重建或至少以高图像质量重建,甚至在k空间采样密度远低于奈奎斯特准则的情况下,并且其也为这种重建提供了方法。在大多数现有的基于CS的MR采集和重建方案中,使用基本的CS公式,其仅利用以下先决条件:MR信号数据在变换域中是稀疏的。还已知的是,由于在某些应用中可能有关于未知MR图像的附加先验信息可用,将该先验信息并入CS重建是有利的。
本发明的技术实现了3D可变密度k空间采样方案,该方案非常适用于具有完全采样的小型中央k空间椭圆体和在***k空间区域中欠采样的CS。
本发明的方法可以与传统视图共享、键孔或k-t采样技术以直接方式组合。可以使用诸如半扫描的加速技术,或者诸如SENSE或迭代SENSE的并行成像技术来采集k空间分布并将其重建为MR图像。
例如,MR图像的时间序列可以从连续采集的MR信号重建,其中,从***k空间位置采集的MR信号在该序列的两个或更多个连续MR图像之间共享(视图共享)以便提高时间分辨率。
迄今为止描述的本发明的方法能够借助MR设备来执行,所述MR设备包括用于在检查体积内生成均匀的静态磁场的至少一个主磁体线圈,用于在在检查体积内在不同空间方向上生成切换的磁场梯度的多个梯度线圈,用于在检查体积内生成RF脉冲和/或用于从位于检查体积中的对象接收MR信号的至少一个RF线圈,用于控制RF脉冲和切换的磁场梯度的时间相继性的控制单元,以及用于根据所接收的MR信号重建MR图像的重建单元。本发明的方法能够例如通过MR设备的重建单元和/或控制单元的对应编程来实施。
目前在临床使用的大多数MR设备中能够有利地执行本发明的方法。为此,仅需要利用计算机程序来控制MR设备,使得其执行本发明的上述方法步骤。该计算机程序可以存在于数据载体上或存在于数据网络中,以便下载以安装在MR设备的控制单元中。
附图说明
附图公开了本发明的优选实施例。然而,应该理解的是,附图的设计仅仅是为了说明的目的,而不是作为本发明的限制的定义。在附图中:
图1示出了用于执行本发明的方法的MR设备;
图2示出了k空间的示意图,其示意性示出本发明的采样方案的实施例的;
图3示出了与通常以接近奈奎斯特极限的较低的径向密度采集的MR图像(图3a)相比,根据本发明在中央k空间切片中以更高的径向密度采集并重建的MR图像(图3b);
图4示出了k空间的示意图,其示意性示出了本发明用于动态成像的应用。
具体实施方式
参考图1,示出了MR设备1。该设备包括超导或电阻性主磁体线圈2,使得沿着通过检查体积的z轴创建基本均匀的时间恒定的主磁场。
磁共振生成和操控***应用一系列RF脉冲和切换的磁场梯度来反转或激励核磁自旋,诱导磁共振,重新聚焦磁共振,操控磁共振,空间地且以其他方式编码磁共振,饱和自旋等来执行MR成像。
更具体地,梯度脉冲放大器3沿着检查体积的x轴,y轴和z轴将电流脉冲施加到全身梯度线圈4、5和6中的选定全身梯度线圈。数字RF频率发射器7经由发送/接收开关8将RF脉冲或脉冲包发送到全身体积RF线圈9,以将RF脉冲发送到检查体积中。典型的MR成像序列由彼此一起取得的短持续时间的RF脉冲段的包组成,并且任何施加的磁场梯度实现对核磁共振的选定操控。RF脉冲用于饱和、激励共振,反转磁化,重新聚焦共振或操控共振,并选择位于检查体积中的身体10的部分。MR信号也被全身体积RF线圈9拾取。
为了生成身体10的有限区域的MR图像,将一组局部阵列RF线圈11、12、13放置成与被选择用于成像的区域邻接。阵列线圈11、12、13能够用来接收由身体线圈RF传输诱导的MR信号。
得到的MR信号由全身体积RF线圈9拾取和/或由阵列RF线圈11、12、13拾取,并且由优选包括前置放大器(未示出)的接收器14解调。接收器14经由发送/接收开关8连接到RF线圈9、11、12和13。
主计算机15控制梯度脉冲放大器3和发射器7以生成诸如回波平面成像(EPI)、回波体积成像、梯度和自旋回波成像,快速自旋回波(TSE)成像等多个MR成像序列中的任意。对于所选择的序列,接收器14在每个RF激励脉冲之后快速连续地接收单个或多个MR数据线。数据采集***16执行接收信号的模数转换,并将每个MR数据线转换成适合于进一步处理的数字格式。在现代MR设备中,数据采集***16是专门用于采集原始图像数据的独立计算机。
最终,数字原始图像数据通过应用傅立叶变换或其他合适的重建算法的重建处理器17被重建为图像表示。MR图像可以表示通过患者的平面切片、平行平面切片的阵列、三维体积等。然后将该图像存储在图像存储器中,在图像存储器中可以访问该图像,以将图像表示的切片、投影或其他部分转换为用于可视化的适当格式,例如经由视频监视器18,其提供得到的MR图像的人类可读的显示。
继续参考图1并进一步参考图2-4,解释本发明的成像方法的实施例。
如图2a所示,根据3D星形叠层采集方案采集MR信号(例如使用传统的梯度回波成像序列)。在多个空间非选择性或板选择性RF激励中的每一个之后,采集一个或多个MR信号,其中,每个MR信号表示k空间分布。MR信号是从多个(在图2a的实施例中为7个)平行切片21、22、23、24、25、26、27被采集为径向k空间分布。切片沿着切片方向kz布置在不同的位置。在kz方向上执行笛卡尔相位编码,同时在沿着围绕中心(kx=ky=0)旋转的径向“辐条”S1、S2、S3、S4、S5、S6、S7、S8的每个单个切片内采集MR信号。这导致由堆叠的盘组成的圆柱形k空间覆盖。在所描绘的实施例中,在采集不同角位置处的k空间分布之前,顺序地采集沿切片方向kz的相位编码步骤。对于辐条S1、S2、S3、S4、S5、S6、S7、S8的角排序,采用上述黄金角度方案。辐条的角度每次递增ΔΦ=111.25°。根据本发明,每个切片的k空间分布的径向密度,即采集的辐条S1、S2、S3、S4、S5、S6、S7、S8的数量根据切片位置而变化,其中,径向密度在更靠中心的k空间位置处更高,并且在更靠***k空间位置处更低。例如,这在图2a的实施例中如下实现:在第一步骤中,从内部的三个切片23、24和25采集辐条S1。对辐条S2和S3相继进行该操作。在对不同黄金角位置处的k空间分布进行采样之前,执行沿着切片方向的相位编码步骤的连续采集,这对确保高度数据一致性和一般的运动鲁棒性至关重要。在第二步骤中,仅从内部的五个切片22、23、24、25和26采集辐条S4和S5。在第三步骤中,从所有七个切片21、22、23、24、25、26和27采集辐条S6、S7和S8,如黑色箭头所示。这样,k空间中心(kz=0附近)比k空间***更密集地采样。以这种方式改变k空间分布(辐条)的径向密度,使得根据给定的FOV的奈奎斯特准则在围绕k空间的中心(kx=ky=kz=0)的椭圆体29(见图2b)内被满足。椭圆体沿着椭圆29与切片21相交,如图2a所示。在椭圆29内满足奈奎斯特准则,而在椭圆29之外k空间被欠采样。切片22、23、24中的各个椭圆的半径对应地朝向kz=0增加。任选地,由于椭圆体28覆盖足够大的k空间的中心区域,该足够大的椭圆体可以充当k空间中的椭圆体“快门”,其通过在根据所采集的MR信号重建MR图像之前移除椭圆体29之外的k空间位置处的MR信号而被采用。
本发明使符合奈奎斯特准则的k空间采样的总扫描时间最小化。然而,其也能够在图像能量主导的中心k空间中提供过采样。径向采样密度从中心切片23、24、25逐渐减小到***切片21、27中的较低径向采样密度。在椭圆体29外部,径向k空间密度甚至可以低于奈奎斯特阈值,而不会显著影响图像质量。
由于k空间中心周围的较高的径向采样密度,根据本发明在最小的扫描时间内有效地减少了条纹伪影。
图3显示了更高的径向密度的影响。与用较低的径向密度采集的图3a的MR图像相比,在根据本发明的利用较高的径向密度进行采样的图3b的MR图像中有效地减少了条纹伪影。
由于k空间的中心区域的定期更新,本发明的k空间采样技术在动态对比增强MR成像的背景下是感兴趣的。根据本发明,根据连续采集的MR信号重建MR图像的时间序列,其中,相比于从更靠***的k空间位置处的切换采集径向k空间分布,以更高的时间密度来从更靠中心的k空间位置处的切片采集径向k空间分布。
这在图4a中示出,图4a显示作为时间t的函数的k空间采样的进展。采集从在kz=0的中央k空间切片采集三个k空间辐条(黄金角采集)开始(左示意图)。奈奎斯特准则在所描绘的实心圆内被满足。在下一步中,从三个内部切片中的每一个采集另外三个辐条。在接下来的步骤中,从所有五个切片沿着另外三个径向方向采集k空间分布。在所有三个示意图中,满足奈奎斯特准则的区域被表示为实心圆。该分布顺序方案得到了中心3D k空间奈奎斯特椭圆体28,其中,在接近的时间点采集MR信号,而在k空间的椭圆体k空间快门29外侧的***在更远离的时间点进行采样,如图4b所示。
k空间的过采样发生在奈奎斯特椭圆体28内。图4c示出了具有连续采集的辐条S1-S12的图4a的切片叠层的kz=0切片。在k空间的最中心的圆形区域中,仅选择辐条S1、S2、S3的MR信号用于图像重建。这个MR信号数据补充有来自k空间的中间环形区域中的辐条S1、S2、S3、S4、S5、S6的MR信号。MR信号数据由来自k空间的其余***区域中的所有辐条S1-S12的MR信号完成。以这种方式,仅使用分别在最接近的时间点处采集的奈奎斯特椭圆体内的k空间分布与在更远离的时间点采样的从***k空间区域采集的k空间分布相结合来重建时间序列的每个MR图像。在图4a所示的示例中,随着时间的推移,首先对更靠中心的k空间位置进行采样,并且然后对更靠***的k空间位置进行采样。这个顺序可以颠倒。或者,采样位置的kz值可以例如随着时间的推移而线性增加或减小。
该策略实现了改进的对比增强MR成像。相比于现有技术,能够以更高的准确性来跟随造影剂在感兴趣区域中的摄取。能够更清楚地区分造影剂到达特定图像体素的时间点,因为仅在时间上紧密采集的MR信号被用于k空间中心的MR图像重建。因此,降低了所测量的对比度增强分布的模糊。
Claims (13)
1.一种对放置在MR设备(1)的检查体积中的对象(10)的MR成像的方法,所述方法包括:
-使所述对象(10)经受至少一个RF脉冲和切换的磁场梯度的成像序列,
-根据星形叠层方案采集MR信号,其中,所述MR信号在多个不同的角位置处并且从沿切片方向布置在不同位置处的多个平行切片(21-27)被采集为径向k空间分布(S1-S12),其中,所述k空间分布(S1-S12)的径向密度根据切片位置而变化,其中,所述径向密度在更靠中心的k空间位置处更高,并且在更靠***的k空间区域中的位置处更低,其中,
-相比于从更靠***的k空间位置处的切片采集k空间分布,以更高的时间密度来从更靠中心的k空间位置处的切片采集k空间分布(S1-S12),并且
-根据所述MR信号重建一幅或多幅MR图像。
2.根据权利要求1所述的方法,其中,为了重建所述MR图像,沿着k空间中的径向辐条采集的并且在时间上彼此最接近的MR信号在围绕k空间的中心的椭圆体(28)中被选择,在所述椭圆体中,所述k空间分布的径向密度发生变化使得根据给定视场的奈奎斯特准则在围绕k空间的中心的所述椭圆体(28)内被满足。
3.根据权利要求1-2中的任一项所述的方法,其中,为了重建所述MR图像,仅选择来自超过中心椭圆体的***k空间区域中的k空间中的部分径向辐条的MR信号,并且k空间中的这些部分径向辐条从向外的径向k空间位置延伸,并且随着这些部分辐条向外移动至***的k空间MR椭圆体,这些部分辐条具有较高的角密度。
4.根据权利要求1-3中的任一项所述的方法,其中,在重建所述MR图像之前,移除所述椭圆体(29)外部的k空间位置处的MR信号。
5.根据权利要求1-4中的任一项所述的方法,其中,在采集不同的角位置处的另外的k空间分布(S1-S12)之前,在相同的角位置处从至少两个不同的切片(21-27)采集k空间分布(S1-S12)。
6.根据权利要求1-5中的任一项所述的方法,其中,所述k空间分布的径向密度发生变化,使得相比于椭圆体(28)外侧,在所述椭圆体中,使用比根据给定视场的奈奎斯特准则更高的径向密度,并且所述椭圆体(28)外侧的所述k空间分布的径向密度等于或小于奈奎斯特准则。
7.根据权利要求1-6中的任一项所述的方法,其中,连续采集的k空间分布(S1-S12)的角度增量(ΔΦ)是黄金角度。
8.如权利要求1-7中的任一项所述的方法,其中,所述k空间分布(S1-S12)不规则地分布在k空间的***区域中。
9.根据权利要求1-8中的任一项所述的方法,其中,使用压缩感测来重建所述MR图像。
10.根据权利要求1-9中的任一项所述的方法,其中,从连续采集的MR信号重建MR图像的时间序列,其中,从超出所述椭圆体的***k空间位置采集的MR信号在所述序列的两个或更多个连续MR图像之间共享。
11.根据权利要求1-10中的任一项所述的方法,其中,所述MR信号经由多个RF接收天线(11、12、13)并行地被采集,其中,所述MR图像是使用SENSE重建的。
12.一种MR设备,包括:至少一个主磁体线圈(2),其用于在检查体积内生成均匀的静态磁场;多个梯度线圈(4、5、6),其用于在所述检查体积内的不同的空间方向上生成切换的磁场梯度;至少一个RF线圈(9),其用于在所述检查体积内生成RF脉冲和/或用于接收来自位于所述检查体积中的对象(10)的MR信号;控制单元(15),其用于控制RF脉冲和切换的磁场梯度的时间相继性;以及重建单元(17),其用于根据所接收的MR信号重建MR图像,其中,所述MR设备(1)被布置为执行以下步骤:
-使所述对象(10)经受至少一个RF脉冲和切换的磁场梯度的成像序列,
-根据星形叠层方案采集MR信号,其中,所述MR信号从沿着切片方向布置在不同位置处的多个平行切片(21-27)被采集为径向k空间分布(S1-S12),其中,所述k空间分布(S1-S12)的径向密度根据切片位置而变化,其中,所述径向密度在更靠中心的k空间位置处更高并且在更靠***的k空间区域中更低,其中,
-相比于从更靠***的k空间位置处的切片采集k空间分布,以更高的时间密度来从更靠中心的k空间位置处的切片采集k空间分布(S1-S12),并且
-根据所述MR信号重建一幅或多幅MR图像。
13.一种在MR设备上运行的计算机程序,所述计算机程序包括用于以下的指令:
-生成至少一个RF脉冲和切换的磁场梯度的成像序列,
-根据星形叠层方案采集MR信号,其中,所述MR信号从沿着切片方向布置在不同位置处的多个平行切片(21-27)被采集为径向k空间分布(S1-S12),其中,所述k空间分布(S1-S12)的径向密度根据切片位置而变化,其中,所述径向密度在更靠中心的k空间位置处更高并且在更靠***的k空间区域中更低,其中,
-相比于从更靠***的k空间位置处的切片采集k空间分布,以更高的时间密度来从更靠中心的k空间位置处的切片采集k空间分布(S1-S12),并且
-根据所述MR信号重建一幅或多幅MR图像。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP15172087.7 | 2015-06-15 | ||
EP15172087 | 2015-06-15 | ||
PCT/EP2016/063413 WO2016202707A1 (en) | 2015-06-15 | 2016-06-12 | Mr imaging using a stack-of-stars acquisition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107750338A true CN107750338A (zh) | 2018-03-02 |
CN107750338B CN107750338B (zh) | 2020-12-22 |
Family
ID=53397943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201680035043.0A Expired - Fee Related CN107750338B (zh) | 2015-06-15 | 2016-06-12 | 使用星形叠层采集的mr成像方法和设备 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10670678B2 (zh) |
EP (1) | EP3308186B1 (zh) |
JP (2) | JP6865697B2 (zh) |
CN (1) | CN107750338B (zh) |
WO (1) | WO2016202707A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109696647A (zh) * | 2019-02-21 | 2019-04-30 | 奥泰医疗***有限责任公司 | 三维多次激发扩散加权成像的k空间采集方法及重建方法 |
CN111257809A (zh) * | 2020-01-23 | 2020-06-09 | 上海东软医疗科技有限公司 | 磁共振成像方法、装置、存储介质及医疗设备 |
CN112557978A (zh) * | 2020-11-26 | 2021-03-26 | 上海交通大学 | 一种基于三维非笛卡尔轨迹的多回波采样方法 |
CN113614558A (zh) * | 2019-03-14 | 2021-11-05 | 皇家飞利浦有限公司 | 使用具有软运动门控的3d径向或螺旋采集的mr成像 |
US20220146615A1 (en) * | 2020-11-12 | 2022-05-12 | Oxford University Innovation Limited | Methods and systems for optimizing sampling efficiency of dynamically sampled magnetic resonance imaging (mri) |
CN115430828A (zh) * | 2022-09-22 | 2022-12-06 | 济南海圣机电科技有限公司 | 一种浇注机铁水定量定速浇注控制方法 |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10353034B2 (en) * | 2016-05-05 | 2019-07-16 | Siemens Healthcare Gmbh | Multi-echo pseudo-golden angle stack of stars thermometry with high spatial and temporal resolution using k-space weighted image contrast |
EP3413071A1 (en) * | 2017-06-09 | 2018-12-12 | Koninklijke Philips N.V. | Mr imaging using a stack-of-stars acquisition with variable contrast |
US20190054194A1 (en) * | 2017-08-17 | 2019-02-21 | Case Western Reserve University | System and Method For Evaluation of Subjects Using Magnetic Resonance Imaging and Oxygen-17 |
EP3454070A1 (en) | 2017-09-06 | 2019-03-13 | Commissariat à l'Energie Atomique et aux Energies Alternatives | Method and apparatus for accelerated magnetic resonance imaging |
EP3467533A1 (en) * | 2017-10-09 | 2019-04-10 | ETH Zurich | Magnetic resonance imaging method with hybrid filling of k-space |
EP3470869A1 (en) | 2017-10-10 | 2019-04-17 | Koninklijke Philips N.V. | Mr imaging using motion-dependent radial or spiral k-space sampling |
EP3517988A1 (en) | 2018-01-29 | 2019-07-31 | Koninklijke Philips N.V. | Mr imaging using a stack-of-stars acquisition with intrinsic motion correction |
EP3579009A1 (en) * | 2018-06-05 | 2019-12-11 | Koninklijke Philips N.V. | Zero echo time mr imaging with water-fat separation |
EP3598162A1 (en) * | 2018-07-19 | 2020-01-22 | Siemens Healthcare GmbH | Magnetic resonance method and apparatus for suppressing metal artifacts |
US10746830B2 (en) * | 2018-08-28 | 2020-08-18 | General Electric Company | Systems and methods for hybrid slice encoding in three-dimensional magnetic resonance imaging |
EP3709040A1 (de) * | 2019-03-13 | 2020-09-16 | Siemens Healthcare GmbH | Passive magnetfeldkamera und verfahren zum betrieb der passiven magnetfeldkamera |
JP7292930B2 (ja) * | 2019-04-03 | 2023-06-19 | キヤノンメディカルシステムズ株式会社 | 画像再構成方法及び再構成装置 |
EP3748385A1 (en) | 2019-06-04 | 2020-12-09 | Koninklijke Philips N.V. | Optimized k-space profile ordering for 3d radial or spiral mr imaging |
EP3842819A1 (en) | 2019-12-23 | 2021-06-30 | Koninklijke Philips N.V. | 3d mr imaging with intrinsic motion detection |
JP7509546B2 (ja) * | 2020-02-12 | 2024-07-02 | キヤノンメディカルシステムズ株式会社 | 画像再構成装置 |
GB2608404B (en) * | 2021-06-30 | 2023-07-05 | Siemens Healthcare Gmbh | Method of reconstructing a dynamic series of motion-compensated magnetic resonance images |
US11940516B2 (en) | 2022-01-29 | 2024-03-26 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for magnetic resonance imaging |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19816961A1 (de) * | 1998-04-17 | 1999-06-17 | Siemens Ag | Verfahren zur Gewinnung eines dreidimensionalen Rohdatensatzes für die MR-Bildgebung und Vorrichtung zur Durchführung des Verfahrens |
US20090105582A1 (en) * | 2006-04-21 | 2009-04-23 | Lawrence Dougherty | Rapid 3-diamensional bilateral breast MR imaging |
US7853060B2 (en) * | 2006-04-13 | 2010-12-14 | Siemens Aktiengesellschaft | Method and MR system for generating MR images |
CN103229069A (zh) * | 2010-12-02 | 2013-07-31 | 皇家飞利浦电子股份有限公司 | 使用多点狄克逊技术的mr 成像 |
JP2013192578A (ja) * | 2012-03-16 | 2013-09-30 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置及びプログラム |
CN103477238A (zh) * | 2011-03-24 | 2013-12-25 | 皇家飞利浦有限公司 | 采用来自先验采集的约束的压缩感测mr图像重建 |
CN104068859A (zh) * | 2013-03-25 | 2014-10-01 | 西门子公司 | 用于确定多个磁共振图像的方法和磁共振设备 |
CN104395772A (zh) * | 2012-03-19 | 2015-03-04 | 皇家飞利浦有限公司 | 具有在对中心和外周k空间区的采样期间的呼吸运动检测的磁共振图像重建方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6404194B1 (en) | 1994-08-05 | 2002-06-11 | The Board Of Trustees Of The Leland Stanford Junior University | Three dimensional magnetic resonance imaging using spiral trajectories in k-space |
US5713358A (en) | 1996-03-26 | 1998-02-03 | Wisconsin Alumni Research Foundation | Method for producing a time-resolved series of 3D magnetic resonance angiograms during the first passage of contrast agent |
US5881728A (en) | 1996-07-26 | 1999-03-16 | Wisconsin Alumni Research Foundation | Digital subtraction magnetic resonance angiography with image artifact suppression |
US5873825A (en) | 1997-04-11 | 1999-02-23 | Wisconsin Alumni Research Foundation | Three dimensional digital subtraction magnetic resonance angiography with limited k-space mask |
JPH11276456A (ja) * | 1998-03-30 | 1999-10-12 | Shimadzu Corp | Mrイメージング装置 |
WO2006120584A2 (en) * | 2005-02-11 | 2006-11-16 | Koninklijke Philips Electronics N.V. | Magnetic resonance imaging with short echo times |
US8120358B2 (en) | 2006-04-13 | 2012-02-21 | The Regents Of The University Of California | Magnetic resonance imaging with high spatial and temporal resolution |
JP5643789B2 (ja) * | 2006-09-06 | 2014-12-17 | 株式会社東芝 | 磁気共鳴イメージング装置 |
WO2009082531A1 (en) * | 2007-12-20 | 2009-07-02 | Wisconsin Alumni Research Foundation | Method for prior image constrained image reconstruction |
EP2350698B1 (en) * | 2008-10-14 | 2020-04-01 | Mayo Foundation For Medical Education And Research | System and method for moving table mri |
WO2012102338A1 (ja) | 2011-01-28 | 2012-08-02 | 株式会社 日立メディコ | 磁気共鳴イメージング装置及びそれを用いた血管撮像方法 |
DE102011078273B4 (de) * | 2011-06-29 | 2014-11-20 | Siemens Aktiengesellschaft | MR-Angiographie mit nicht-kartesischer Signalakquisition |
WO2013159044A1 (en) | 2012-04-19 | 2013-10-24 | New York University | System, method and computer-accessible medium for highly-accelerated dynamic magnetic resonance imaging using golden-angle radial samplng and compressed sensing |
US8744551B2 (en) | 2012-07-12 | 2014-06-03 | Northshore University Healthsystem | Method for non-contrast enhanced magnetic resonance angiography |
WO2014085288A1 (en) | 2012-11-29 | 2014-06-05 | The Regents Of The University Of California | Noninvasive 4-d time-resolved dynamic magnetic resonance angiography |
US10052033B2 (en) * | 2015-03-12 | 2018-08-21 | Siemens Healthcare Gmbh | System and method for motion-robust 3D magnetic resonance imaging of vessel walls |
-
2016
- 2016-06-12 US US15/580,402 patent/US10670678B2/en not_active Expired - Fee Related
- 2016-06-12 CN CN201680035043.0A patent/CN107750338B/zh not_active Expired - Fee Related
- 2016-06-12 JP JP2017564875A patent/JP6865697B2/ja active Active
- 2016-06-12 WO PCT/EP2016/063413 patent/WO2016202707A1/en active Application Filing
- 2016-06-12 EP EP16730779.2A patent/EP3308186B1/en active Active
-
2021
- 2021-01-07 JP JP2021001435A patent/JP2021090760A/ja active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19816961A1 (de) * | 1998-04-17 | 1999-06-17 | Siemens Ag | Verfahren zur Gewinnung eines dreidimensionalen Rohdatensatzes für die MR-Bildgebung und Vorrichtung zur Durchführung des Verfahrens |
US7853060B2 (en) * | 2006-04-13 | 2010-12-14 | Siemens Aktiengesellschaft | Method and MR system for generating MR images |
US20090105582A1 (en) * | 2006-04-21 | 2009-04-23 | Lawrence Dougherty | Rapid 3-diamensional bilateral breast MR imaging |
CN103229069A (zh) * | 2010-12-02 | 2013-07-31 | 皇家飞利浦电子股份有限公司 | 使用多点狄克逊技术的mr 成像 |
CN103477238A (zh) * | 2011-03-24 | 2013-12-25 | 皇家飞利浦有限公司 | 采用来自先验采集的约束的压缩感测mr图像重建 |
JP2013192578A (ja) * | 2012-03-16 | 2013-09-30 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置及びプログラム |
CN104395772A (zh) * | 2012-03-19 | 2015-03-04 | 皇家飞利浦有限公司 | 具有在对中心和外周k空间区的采样期间的呼吸运动检测的磁共振图像重建方法 |
CN104068859A (zh) * | 2013-03-25 | 2014-10-01 | 西门子公司 | 用于确定多个磁共振图像的方法和磁共振设备 |
Non-Patent Citations (3)
Title |
---|
HIMANSHU BHAT 等: "Contrast-Enhanced Whole-Heart Coronary Magnetic Resonance Angiography at 3T with Radial EPI", 《MAGNETIC RESONANCE IN MEDICINE》 * |
LI FENG 等: "K-t Radial SPARSE-SENSE: Combination of Compressed Sensing and Parallel Imaging with Golden Angle Radial Sampling for Highly Accelerated Volumetric Dynamic MRI", 《PROCEEEDINGS OF TJE INTERNATIAL SOCIETY FOR MAGNETIC RESONANCE IN MEDICINE,20TH ANNUAL MEETING & EXHIBITION 》 * |
ZHITAO LI 等: "Highly Accelerated 3D Dynamic Imaging with Variable Density Golden Angle Stack-of-Stars Sampling", 《PROCEEDINGS OF THE INTERNATIONAL SOCIETY FOR MAGNETIC RESONANCE IN MEDICINE》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109696647A (zh) * | 2019-02-21 | 2019-04-30 | 奥泰医疗***有限责任公司 | 三维多次激发扩散加权成像的k空间采集方法及重建方法 |
CN109696647B (zh) * | 2019-02-21 | 2021-05-28 | 奥泰医疗***有限责任公司 | 三维多次激发扩散加权成像的k空间采集方法及重建方法 |
CN113614558A (zh) * | 2019-03-14 | 2021-11-05 | 皇家飞利浦有限公司 | 使用具有软运动门控的3d径向或螺旋采集的mr成像 |
CN111257809A (zh) * | 2020-01-23 | 2020-06-09 | 上海东软医疗科技有限公司 | 磁共振成像方法、装置、存储介质及医疗设备 |
US20220146615A1 (en) * | 2020-11-12 | 2022-05-12 | Oxford University Innovation Limited | Methods and systems for optimizing sampling efficiency of dynamically sampled magnetic resonance imaging (mri) |
US11802926B2 (en) * | 2020-11-12 | 2023-10-31 | Oxford University Innovation Limited | Methods and systems for optimizing sampling efficiency of dynamically sampled magnetic resonance imaging (MRI) |
CN112557978A (zh) * | 2020-11-26 | 2021-03-26 | 上海交通大学 | 一种基于三维非笛卡尔轨迹的多回波采样方法 |
CN112557978B (zh) * | 2020-11-26 | 2022-03-25 | 上海交通大学 | 一种基于三维非笛卡尔轨迹的多回波采样方法 |
CN115430828A (zh) * | 2022-09-22 | 2022-12-06 | 济南海圣机电科技有限公司 | 一种浇注机铁水定量定速浇注控制方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2016202707A1 (en) | 2016-12-22 |
EP3308186B1 (en) | 2022-02-02 |
CN107750338B (zh) | 2020-12-22 |
JP2018519050A (ja) | 2018-07-19 |
US20180149721A1 (en) | 2018-05-31 |
EP3308186A1 (en) | 2018-04-18 |
JP2021090760A (ja) | 2021-06-17 |
JP6865697B2 (ja) | 2021-04-28 |
US10670678B2 (en) | 2020-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107750338A (zh) | 使用星形叠层采集的mr成像 | |
CN103477238B (zh) | 采用来自先验采集的约束的压缩感测mr图像重建 | |
CN106030329B (zh) | 使用多回波分段k空间采集的MR成像方法 | |
JP6998218B2 (ja) | 動き検出を用いるmr撮像 | |
US6954067B2 (en) | Three-dimensional phase contrast imaging using interleaved projection data | |
CN105074491B (zh) | 具有使用压缩感知的图像重建的动态mri | |
CN105143906B (zh) | 金属抗mr成像 | |
CN104204839B (zh) | 使用apt对比增强和多回波时间采样的mr成像 | |
CN104685368B (zh) | 用于耐金属mr成像的方法和设备 | |
CN104714199B (zh) | 一种磁共振成像方法和装置 | |
US11067653B2 (en) | Magnetic resonance imaging method including acquiring echo signals from each k-space slice with different relaxation time weightings after echo signals from different k-space splices are acquired | |
CN106574954A (zh) | 针对epi的具有奈奎斯特伪影校正的并行mr成像 | |
CN105051563B (zh) | 使用相位调制rf脉冲的并行多切片mr成像 | |
US20120262172A1 (en) | Motion correction in accelerated t1-weighted magnetic resonance imaging | |
CN105005011A (zh) | 在三维快速自旋回波中的混合采集的磁共振成像方法 | |
CN107209238A (zh) | 具有对边带伪迹的抑制的并行多切片mr成像 | |
CN106796274B (zh) | 具有伪迹抑制的propeller-mr成像 | |
US20120046539A1 (en) | Dual-contrast mr imaging using fluid-attenuation inversion recovery (flair) | |
CN107548464A (zh) | 自旋回波mr成像 | |
CN107076818A (zh) | 零回波时间mr成像 | |
JP2021512671A (ja) | 固有運動補正を有するスタック・オブ・スター収集を用いたmr撮像 | |
WO2020078907A1 (en) | Magnetic resonance imaging using motion-compensated image reconstruction | |
US20200300952A1 (en) | Mr imaging using motion-dependent radial or spiral k-space sampling | |
JP2021505296A (ja) | 脂肪抑制を用いる拡散mrイメージング | |
EP3185029A1 (en) | Mr imaging using propeller acquisition with t2 decay correction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201222 |