CN101351721A - 床连续移动的灵敏度编码磁共振并行成像 - Google Patents

床连续移动的灵敏度编码磁共振并行成像 Download PDF

Info

Publication number
CN101351721A
CN101351721A CNA2006800112918A CN200680011291A CN101351721A CN 101351721 A CN101351721 A CN 101351721A CN A2006800112918 A CNA2006800112918 A CN A2006800112918A CN 200680011291 A CN200680011291 A CN 200680011291A CN 101351721 A CN101351721 A CN 101351721A
Authority
CN
China
Prior art keywords
space
data
coil
coil sensitivity
formation method
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.)
Pending
Application number
CNA2006800112918A
Other languages
English (en)
Inventor
J·柯普
P·博尔纳特
B·阿尔德菲尔德
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN101351721A publication Critical patent/CN101351721A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image 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/56375Intentional motion of the sample during MR, e.g. moving table imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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/5611Parallel 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

Landscapes

  • 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)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

在连续移动成像对象(12)使其通过扫描仪视场(20)期间,采用多个视频线圈(24,26)来采集k-空间数据。该采集包括在至少一个欠采样方向上欠采样k-空间。对于该至少一个欠采样方向限定从k-空间至真实空间的加权变换(62)。该加权变换组合了与患者位置相关的线圈灵敏度加权因子和傅立叶变换。对所采集的k-空间数据沿连续移动方向进行混合变换,以限定具有连续移动变换方向上的真实空间尺寸和横过该连续移动方向的横向方向上的k-空间尺寸的混合空间数据。沿该横向方向变换混合空间数据以产生重建图像。混合变换和该变换应用视欠采样的相应方向而定的已限定加权变换(62)。

Description

床连续移动的灵敏度编码磁共振并行成像
本发明涉及磁共振成像领域。其特别应用于癌症筛查、血管造影术和其它对全身或者身体大部分有利地进行的诊断成像方法,并将特别参考进行描述。但是,其通常还应用于连续移动工作台的磁共振成像。
磁共振扫描仪具有通常比平均人体尺寸小得多的有限视场。因此,磁共振成像不易用于进行“全身”或者其它扩大的感兴趣区域的成像,其中所感兴趣区域大于扫描仪视场。这限制了用于诊断任务例如癌症筛查、血管造影术等等的磁共振成像的有效性。
一种对大于扫描仪视场的扩大的感兴趣区域成像的方法为多站方法,其中成像以不连续的步骤进行。在每一个步骤成像一个扫描仪视场。随后,成像对象被推进通过扫描仪腔等于扫描仪视场轴向长度的距离,并且成像另一个扫描仪视场。将连续成像的视场拼接在一起以形成扩大的感兴趣区域的图像。多站方法具有某些缺陷。起动-停止工作台的动作会烦扰患者、增加患者运动伪像并增加了总的成像持续时间。还在连续扫描仪视场之间边界的组合图像上引入不连续点。
在连续移动工作台磁共振成像中,通常以恒定的工作台速度连续轴向推进患者通过扫描仪腔,在连续的工作台推进期间采集磁共振成像数据。和多站方法相比,因为工作台连续移动而不是以“起动-停止”模式移动,所以采用连续移动工作台的成像通常更快、减少或者消除了图像不连续点、并且通常不太烦扰患者。
并行成像是另一种进行较快扫描的技术。在并行成像中,采用多个线圈接收磁共振信号。例如,在灵敏度编码(SENSE)成像中,稀疏采样k-空间,并将该多个线圈的数据用于补偿任一个线圈采样的稀疏性。例如采用四个线圈时,仅仅可获得每条第四相位编码线。从每个线圈采集的数据重建分散的图像,同时产生一组和该四个线圈对应的四个“折叠”重建图像。考虑到线圈灵敏度因子组合、或者展开这些折叠图像以产生最终的图像,其特征基本上和通过每条相位编码线的单个线圈采样所获得的图像数据相同。在该实例中,因为仅仅采集每个第四相位编码线,所以可获得成像速度提高四个因子。基于线圈数和SENSE因子,还可另外提高成像速度。其它并行成像技术例如SMASH也有各种优势,例如提高成像速度、图像分辨率较高,等等。
因为当患者移动通过扫描仪视场时工作台移动引入与时间相关的线圈灵敏度因子,所以难以进行采用连续移动工作台和固定接收线圈的并行成像。可替换地采用可移动线圈(例如设置在移动工作台上或者移动工作台内,或直接置于患者上面)。但是,因为工作台平移将移动的线圈移位至扫描仪视场外,所以患者的整个身体都必须覆盖该种线圈。对某些应用,需要混合设置,其中并行成像阵列的一些线圈相对于扫描仪固定,而其它线圈随着患者移动。在这些混合设置中,固定线圈具有与时间相关的线圈灵敏度因子,而可移动线圈具有与时间无关的灵敏度因子。对于采用连续移动工作台和固定或者混合线圈设置的并行成像,标准重建方法(SENSE、SMASH)不适用。
本发明考虑了克服上述限制和其它问题的改进装置和方法。
根据一个方面,提供一种成像方法。连续移动成像对象通过扫描仪视场。在连续移动期间,采用多个射频线圈采集k-空间数据。该采集包括在至少一个欠采样方向上欠采样k-空间。对采集的k-空间数据沿连续移动方向进行傅立叶变换,以限定具有连续移动变换方向上的真实空间尺寸和横过连续移动方向的横向方向上的k-空间尺寸的混合空间数据。沿该横向方向变换该混合数据以产生重建图像。对于该至少一个欠采样方向限定从射频线圈的大量k-空间数据至真实空间的加权变换。该加权变换组合了与患者位置相关的线圈灵敏度权重因子。
根据另一方面,公开一种磁共振成像扫描仪,其执行上段所述的成像方法。
根据另一方面,公开一种处理器,用于对在成像对象连续移动期间使用多个射频线圈所采集的k-空间数据执行图像数据处理方法并包括在至少一个欠采样方向上的k-空间欠采样。通过该处理器执行的图像数据处理方法包括:沿连续移动方向混合变换所采集的k-空间数据,以限定具有连续移动变换方向上的真实空间尺寸和在横过连续移动方向的横向方向上的k-空间尺寸的混合空间数据;对于该至少一个欠采样方向限定从所有射频线圈的k-空间至真实空间的加权变换,该加权变换组合了与患者位置相关的线圈灵敏度权重因子并包括沿该横向方向的混合空间数据的傅立叶变换以产生重建图像;该方法沿欠采样方向应用所限定的加权变换。
根据另一方面,公开一种存储介质,其通过相关数字处理器编码可执行指令以对在成像对象连续移动期间使用多个射频线圈所采集的k-空间数据进行图像数据处理方法并包括在至少一个欠采样方向上的k-空间欠采样。该图像数据处理方法包括:沿连续移动方向混合变换所采集的k-空间数据,以限定具有连续移动变换方向上的真实空间尺寸和在横过连续移动方向的横向方向上的k-空间尺寸的混合空间数据;对于该至少一个欠采样方向限定从所有射频线圈的k-空间至真实空间的加权变换,该加权变换组合了与患者位置相关的线圈灵敏度权重因子并包括沿该横向方向的混合空间数据的傅立叶变换以产生重建图像;该方法沿欠采样方向应用所限定的加权变换。
一个优点在于通过组合连续移动工作台磁共振成像和并行成像而较快成像。
另一个优点在于通过采用结合并行成像进行的连续移动工作台磁共振成像而并行成像的视场比扫描仪视场大。
另一个优点在于以固定线圈促进连续移动工作台磁共振成像和并行成像的组合。
另一个优点在于以固定线圈和移动线圈促进连续移动工作台磁共振成像和并行成像的组合。
另一个优点在于将线圈负载效应可选择地合并到与患者位置相关的、并且结合到重建的线圈灵敏度因子。
本领域技术人员在阅读下面优选实施例的详细描述时可清楚多种其它优势和益处。
本发明具体化为各种部件以及部件的设置,以及各种过程操作以及过程操作的设置。附图仅仅为描述优选实施例而不应解释为限制本发明。
图1示意性地示出了连续移动工作台磁共振成像系统的实例,其包括采用多个磁共振接收线圈的并行成像。
图2示意性地示出了重叠在扫描仪视场两次连续扫描同时成像对象连续移动时采集的有效视场。图2另外还示意性地示出了在重建每个有效视场中采用的与患者位置相关的线圈灵敏度权重因子。
图3示意性地示出了从图2扫描所示出的图像,并添加有效视场和重叠区域。
参考图1,磁共振成像扫描仪10对位于连续移动工作台或者床14上的成像对象12进行成像。该连续移动工作台以图1中被标记的块箭头所示的速度v进行移动。该磁共振成像扫描仪10通常包括各种包含在外罩16中的部件并由此在图1中不可见,例如至少在其中设置有一部分成像对象12的扫描仪视场20(在图1的透视图中扫描仪视场20被扫描仪罩16从直接视野遮挡,因此在图1中以虚线表示)内产生基本上时间和空间恒定的B0磁场的主磁场线圈,和至少在扫描仪视场20中选择性地产生磁场梯度的磁场梯度线圈。
另外,多个射频线圈选择性地与扫描仪视场20中的部分成像对象耦合。在图1的实例中,这些线圈包括设置在外罩16上面或者里面的正交鸟笼式线圈22;该正交鸟笼式线圈22的端环在图1中以虚线表示。图1的示例线圈还包括设置在扫描仪外罩16上的固定线圈24。图1的示例线圈还包括嵌入在连续移动工作台14内的可移动线圈26(在图1的透视图中可以看到两个固定线圈24和一个可移动线圈26)。
成像时,扫描仪10在扫描仪视场20内产生基本上空间和时间恒定的B0磁场,正交鸟笼式线圈22在1H氢拉莫尔频率或者另一合适的磁共振频率下注入射频激励以激励成像对象12的磁共振信号选择部分。在并行成像模式下将固定线圈24用于采集磁共振成像数据。可选地,可移动线圈26也用于并行成像。在另一种考虑的并行成像模式中,仅仅采用可移动线圈26采集磁共振成像数据。成像时,扫描仪10产生沿选择的k-空间轨迹编码磁共振信号的磁场梯度。笛卡儿编码通常包括读出方向以及一个或两个相位编码方向(分别对于两维和三维成像)。可替换地,使用多个两维切片采集来进行三维成像。同时在成像期间,工作台14以速度v连续移动。可选地,对于患者解剖结构的不同部分可改变速度值。通过工作台位置监视器28来监视工作台14的位置。
为了利用通过采用多个通常具有不同线圈灵敏度的线圈24,26来采集k-空间数据而进行的并行成像,该采集应当包括在至少一个欠采样方向上k-空间的欠采样。例如,如果采用笛卡儿k-空间采集轨迹,则通过从k空间采集省略一些相位编码线可适合地获得相位编码方向上的欠采样。通过欠采样可改进成像速度。通过在后续图像重建中包括采用多个线圈采集的k-空间数据来补偿由欠采样所引起k-空间数据的稀疏性。将所采集的k-空间磁共振成像样本存储在k-空间存储器30中,而每个存储的k-空间样本以由工作台位置监视器28所提供的、表示在采集k-空间样本时工作台14位置的值来表示。可替换地或者另外地,因为通过速度v(或者如果对不同解剖区域采用不同速度则通过所应用的多个速度值)使时间和工作台位置相关,所以可通过采集时间来表示每个k-空间值。
为进行并行成像,灵敏度处理器32对合适的线圈灵敏度校准数据采集进行灵敏度分析以确定线圈24,26的线圈灵敏度sγ(r)34(这里r表示空间位置)。在一个合适的方法中,线圈灵敏度校准数据采集包括通过使用鸟笼式线圈22和每个线圈24,26获得的图像扫描。通过使用每个局部线圈24,26采集的比值复杂强度值以及通过使用鸟笼式线圈22采集的相应复杂强度值提供了线圈灵敏度34。该线圈灵敏度采集数据采集为低分辨率,并且例如可通过采集对于工作台14的多个固定位置的数据或者通过采用连续工作台移动的校准扫描在连续工作台移动扫描之前采集。可替换地,可在连续工作台移动扫描期间将低分辨率线圈灵敏度校准数据采集插入在k-空间数据的采集中。
继续参考图1并简单参考图2,当以速度v连续移动成像对象12时,在扫描仪视场20内进行一次或多次扫描。但是,因为在每次扫描的采集期间以速度v连续移动成像对象12,所以每次扫描相对于成像对象12的有效视场较大。通常,如果扫描仪视场20的每次扫描花费时间Tscan,则相对于成像对象12每次扫描的有效视场以大约相当于|v|·Tscan的因子加长或者增加,这里|v|为工作台速度值。图2示出两个这种具有重叠区域Roverlap的有效视场FOVcff1和FOVcff2。通过放射线学者或者其它用户经用户界面36选择扫描参数例如扫描时间Tscan和速度v。
再参考图1,重建处理器40对通过工作台14连续移动所采集的并行成像数据进行重建以产生重建图像,并存储在重建图像存储器42中。通过在重建中包括以多个射频线圈采集的k-空间数据来补偿至少一个欠采样方向上的k-空间欠采样。但是,因为连续移动,线圈灵敏度通常与时间相关,因此不能准确进行图像空间中的常规SENSE-型展开。相反地,如这里所述,将与患者位置相关的线圈灵敏度权重因子直接组合到k-空间至真实空间的变换中,同时允许适当解释线圈灵敏度的时间相关性。可将重建图像适当显示在用户界面36或者另一个更高分辨率的显示设备上,或者打印、通过因特网或者局域网通讯、存储在非易失性存储介质上,或者以其它方式使用。
重建处理器40对说明了线圈灵敏度34和工作台14连续移动的并行成像数据进行重建。一维混合空间变换处理器50沿连续移动方向(即平行于工作台速度v)对所采集的k-空间数据进行傅立叶变换,以限定具有在连续移动变换方向上的真实空间尺寸并保持横过连续移动方向的横向方向上的k-空间尺寸的混合空间数据。在所描述的成像中,k-空间采集使用笛卡儿k-空间轨迹,而读出方向对应于连续移动方向,被标记的相位编码(p.e.)方向对应于横向方向。对于三维成像,可采用横过读出方向和被标记的p.e.方向的第二(未标记的)相位编码方向。和所述处理被标记的p.e.方向相类似,在重建期间适当地处理第三方向。使读取方向与连续移动方向对准的好处在于,容易在读取方向上过采样该数据,这有利于重建无过渡图像。但是,可采用其它k-空间轨迹,例如沿连续移动方向采用相位编码的k-空间轨迹。
将混合空间数据存储在混合空间数据存储器52中。横向变换处理器56沿横向方向变换混合空间数据以产生重建图像,并存储在重建图像存储器42中。如果采用两个相位编码方向进行三维成像(其一个或者两个可以次采样),则沿每个相位编码方向进行一维变换。
在一个合适的实施例中,一维混合空间变换处理器50沿连续移动方向对所采集的k-空间数据进行傅立叶变换,以产生未校正的混合空间数据,并对于连续移动调整未校正的混合空间数据以产生混合空间数据。对于连续移动的校正调整基于通过工作台位置监视器28注释为所采集k-空间数据的工作台位置数值,或者等效地基于所注释的采集时间,连同工作台速度v的已知速度。例如在Kruger等人的“ContinuouslyMoving Table Data Acquisition Method for Long FOV Contrast-Enchanced MRA and Whole Body MRI”,Magnetic Resonance inMedicine(MRM)volume 47,pages 224-31(2002)中描述了合适的校正调整技术。经过对于连续移动进行校正后,混合空间数据对应于由成像对象12所限定的参考帧。因此,应用横向变换处理器56将混合数据变换为重建图像。
但是,因为k-空间数据在至少一个相位编码方向上欠采样,所以变换处理器50,56应当考虑线圈灵敏度以通过组合来自多个射频线圈24,26的k-空间数据来补偿欠采样。通过灵敏度处理器32测量线圈灵敏度sγ(r)34。为解释连续移动,加权变换定义处理器60构造与患者位置相关的线圈灵敏度权重因子,这些因子在功能上与由以速度v在时间t连续移动所产生空间位移v·t偏置的基准位置值r(通过灵敏度处理器32对于和多个工作台位置相应的参考位置r检测的Sγ(r)34)相关。将与患者位置相关的线圈灵敏度权重因子直接组合到k-空间至真实空间的变换以对至少一个欠采样方向限定加权变换62。
将所描述的实施例用作实例。在所描述的实施例中,读取方向与工作台速度v反平行,而相位编码方向横过工作台速度v。对相位编码方向进行欠采样以提供SENSE-型数据,但另外的连续移动叠加至所采集的k-空间数据。一维混合空间变换处理器50沿读取方向对所采集的k-空间数据进行傅立叶变换而不利用线圈灵敏度因子,并对连续移动进行校正调整。这产生了具有转换成真实空间的读取方向但欠采样的相位编码方向保持在k-空间坐标中的混合空间数据。对不进行连续移动的固定线圈24而言,适合应用于相位编码方向的一维加权变换62对应于解决矩阵方程:
m i , γ , κ , λ hyb = Σ λ E i , γ , κ , λ · ρ i , λ - - - ( 1 )
对于真实空间信号密度ρi,λ,这里miγ,k,λ hyb表示混合空间数据,i表示混合空间数据中的相位编码线,λ表示沿混合空间数据相位编码线的横向位置,ρi,λ表示在横向位置λ的线i上的信号密度,γ表示多个射频线圈,k表示相位编码线集合,而Ei,γ,k,λ为对应下式的编码矩阵:
Ei,γ,κ,λ=sγ(ri,λ-v·tκ)·exp(ik(tκ)·ri,λ)  (2)
这里v表示连续移动的速度,k(tk)表示相位编码步长,tk表示以k表示的相位编码线集合的采集时间,ri,λ表示对应于线i和横向位置λ的空间位置,sγ表示线圈密度。因为时间(或者等效地,患者位置)相关性,编码矩阵Ei,γ,k,λ的灵敏度编码和傅立叶编码方面不可分离;相反,必须完全反转Ei,γ,k,λ。例如可采样LU-分解来计算伪逆矩阵,并对真实空间信号密度ρi,λ逐行解答方程(1)。
根据线圈类型来处理线圈灵敏度数据。对于不连续移动的固定线圈24,与患者位置相关的线圈灵敏度权重因子根据相位编码步骤的如式子sγ(ri,λ-v·tk)所示的记录位置来改变灵敏度。对于连续移动的移动线圈而言,选择当前处于扫描仪视场20内的灵敏度模式部分,从而限定与患者位置无关的线圈灵敏度权重因子sγ(ri,λ)。
通常通过在连续移动期间完成扫描仪视场20的多次扫描来实现成像。将转换处理器50,56应用于扫描仪视场20的每次扫描以为每次扫描产生相应的重建图像。如果对于并行成像只采用固定线圈24,则可在每次扫描的重建中采用相同的加权变换62。即,因为对每次扫描都相同,因此不必对每次扫描重新限定加权变换62。另一方面,如果在并行成像中既采用固定线圈24又采用移动线圈26,则移动线圈26平移通过(也可能离开)扫描仪视场20将通常表示对每次扫描重新限定加权变换62。
可选地,加权变换62还考虑由连续移动通常空间不均匀的成像对象12通过扫描仪视场20所造成的与时间相关或者与患者位置相关的线圈负载效应。例如,如果灵敏度处理器32处理以处于适当位置的成像对象12所采集的线圈灵敏度对准数据集合,则得到的线圈灵敏度34将包括以成像对象12装载线圈的效应。这种情况下,即使仅仅采用固定线圈24,当成像对象12移动通过扫描仪视场20时对后来每次的扫描重新限定加权变换62以解释线圈负载的变化是有利的。
通过适当选择相对于在连续移动方向上扫描仪视场20的长度的工作台速度v,其可设置用于从扫描仪视场20的随后扫描重建的空间相邻图像以精确邻接,从而产生包括随后扫描有效视场的放大视场组合扫描。但是,在邻接点可产生图像的不连续。在连续移动方向(在实例k-空间轨迹中的读取方向)上的过采样可使随后的扫描在邻接界面上平滑。
因此,参考图2,选择工作台速度v以部分重叠空间相邻的重建图像是有利的。在图2所示的实例中,两次随后扫描的FOVeff1和FOVeff2在图像空间上重叠一个量Roverlap。可将该冗余用于平滑相邻图像之间的不连续。图2还对FOVeff1的扫描示出了以[sγ(r+vt)]1表示的与患者位置相关的线圈灵敏度因子,对FOVeff2的扫描示出了以[sγ(r+vt)]2表示的与患者位置相关的线圈灵敏度因子。在重叠区域Roverlap通常期望这些灵敏度函数是不同的。因此,采用分别应用于每次扫描k-空间数据的重建处理器40来分别重建两个视场FOVeff1和FOVeff2,以产生与两个视场FOVeff1和FOVeff2对应的两个部分重叠的真实空间图像是有利的。然后在图像空间内将该两个图像的重叠部分组合以例如通过平方和方法来平滑重叠区域Roverlap的不连续。为促进平滑重叠区域Roverlap,在读取方向上过采样是有利的。
图3示意性地示出了得到的组合图像的选择冠状切片,具有两个视场FOVeff1和FOVeff2以及重叠区域Roverlap,随同示出了相位编码(p.e.)和读取方向为该特别实例成像配置所选择的工作台移动速度v。
应理解可以以多种方法物理地表现重建处理器40。例如,重建处理器40可以是可编程数字计算机或者处理器,例如光盘、磁盘、网络服务器非易失性存储器等等(未示出)的存储介质编码可由可编程计算机或者处理器执行的指令以进行重建处理。
如果采用两个相位编码方向以采集三维图像数据,则可对选择的相位编码方向进行欠采样或者对两个相位编码方向进行欠采样。在这种情况下,加权变换定义处理器60限定两个加权变换,对于每个欠采样的相位编码方向各一个。可采用其它k-空间轨迹设置,例如使编码方向与连续移动方向平行或者反平行,使读取方向横过连续移动方向。而且,在连续移动方向上的k-空间采样以及横过连续移动方向上的k-空间采样均可欠采样。这种情况下,加权变换定义处理器60限定两个加权变换,对于每个欠采样方向各一个,并且该两个变换处理器50,56的每一个使用适当的加权变换。
已经参考优选实施例描述了本发明。显然,对于阅读和理解上述详细描述的其它人来说可进行更改和变更。期望将本发明理解为包括所有落入附加权利要求或者其等同物范围内的那些更改和变更。

Claims (18)

1.一种成像方法,包括:
使成像对象(12)连续移动通过扫描仪视场(20);
在连续移动期间,使用多个射频线圈(24,26)采集k-空间数据,该采集包括在至少一个欠采样方向上对k-空间进行欠采样;
对于该至少一个欠采样方向限定从k-空间至真实空间的加权变换(62),该加权变换(62)组合与患者位置相关的线圈灵敏度权重因子和傅立叶变换;
对采集的k-空间数据沿连续移动方向进行变换,以限定具有连续移动变换方向上的真实空间尺寸和横过连续移动方向的横向方向上的k-空间尺寸的混合空间数据;以及
沿该横向方向变换该混合数据以产生重建图像;
该变换包括沿欠采样方向应用所限定的加权变换(62)。
2.根据权利要求1的成像方法,其中沿连续移动方向的变换包括:
对采集的k-空间数据沿连续移动方向进行傅立叶变换,以产生未校正的混合空间数据;以及
调整连续移动的未校正混合空间数据。
3.根据权利要求1的成像方法,其中该加权变换(62)包括:
与患者位置相关的线圈灵敏度权重因子,这些因子在功能上取决于通过由连续移动速度v所产生的空间位移偏置的基准位置值。
4.根据权利要求3的成像方法,其中该加权变换(62)还包括:
由与患者位置相关的线圈灵敏度因子之一来确定的傅立叶变换。
5.根据权利要求1的成像方法,其中所述线圈(24,26)包括(i)一个或多个不连续移动的固定线圈(24),和(ii)一个或多个连续移动的可动线圈(26),并且所述加权变换(62)包括:
与一个或多个固定线圈(24)关联的一个或多个与患者位置相关的线圈灵敏度权重因子;
与一个或多个可动线圈(26)关联的一个或多个与患者位置无关的线圈灵敏度权重因子。
6.根据权利要求1的成像方法,还包括:
对于成像对象(12)相对于扫描仪视场(20)的多个不同位置采集多个线圈灵敏度校准数据集合;
基于线圈灵敏度校准数据集合,计算在所述多个不同位置的每一个上多个射频线圈(24,26)的线圈灵敏度;以及
从在多个不同位置的每一个上已计算的线圈灵敏度获得与患者位置相关的线圈灵敏度权重因子。
7.根据权利要求6的成像方法,其中采集多个线圈灵敏度数据集合包括:
在连续移动期间将线圈灵敏度校准数据集合插入在k-空间数据的采集中。
8.根据权利要求1的成像方法,其中在连续移动期间采集k-空间数据包括:
在连续移动期间完成扫描仪视场(20)的多次扫描,所述变换用于每次扫描以限定相应于每次扫描的重建图像。
9.根据权利要求8的成像方法,其中所述多个射频线圈为不连续移动的固定线圈(24),并且在混合变换和每次扫描的变换中使用相同的加权变换。
10.根据权利要求8的成像方法,其中选择连续移动的速度以部分重叠从多次扫描重建的空间相邻的图像,并且所述方法还包括:
在图像空间中组合空间相邻的重建图像的重叠图像部分。
11.根据权利要求10的成像方法,其中所述至少一个欠采样方向为横向方向,并且使用多个射频线圈在连续移动期间采集k-空间数据,包括:
在连续移动方向上对k-空间过采样。
12.根据权利要求1的成像方法,其中(i)在连续移动期间采集k-空间数据包括沿笛卡儿坐标采集k-空间数据,同时在相位编码方向上对k-空间欠采样,以及(ii)所述采用所限定的加权变换的变换在相位编码方向上产生未展开的真实空间尺寸。
13.根据权利要求12的成像方法,其中所述相位编码方向为横向方向,并且所述加权变换对应于解答矩阵方程:
m i , γ , κ , λ hyb = Σ λ E i , γ , κ , λ · ρ i , λ
对于ρi,λ,这里mi,γ,k,λ hyb表示混合空间数据,i表示混合空间数据中的相位编码线,λ表示沿混合空间数据的相位编码线的横向位置,ρi,λ表示在横向位置λ的线i上的信号密度,γ表示多个射频线圈(24,26),k表示相位编码线集合,而Ei,γ,k,λ为对应于下式的编码矩阵:
Ei,γ,k,λ=sγ(ri,λ-v·tκ)·exp(ik(tκ)·ri,λ)
这里v表示连续移动速度,k(tk)表示相位编码步长,tk表示以k表示的相位编码线集合的采集时间,ri,λ表示对应于线i和横向位置λ的空间位置,sγ表示线圈灵敏度。
14.根据权利要求12的成像方法,其中所述欠采样方向为横向方向,并且沿笛卡儿坐标采集k-空间数据的读取方向与连续移动的方向平行。
15.一种执行权利要求1所述成像方法的磁共振成像扫描仪(10,40)。
16.一种处理器(40),被编程来对在成像对象(12)连续移动期间使用多个射频线圈(24,26)所采集的k-空间数据执行图像数据处理方法并包括在至少一个欠采样方向上的k-空间欠采样,所述图像数据处理方法包括:(i)对于该至少一个欠采样方向限定从k-空间至真实空间的加权变换(62),所述加权变换组合与患者位置相关的线圈灵敏度权重因子和傅立叶变换;(ii)沿连续移动方向变换所采集的k-空间数据,以限定具有连续移动变换方向上的真实空间尺寸和横过连续移动方向的横向方向上的k-空间尺寸的混合空间数据;以及(iii)沿横向方向变换混合空间数据以产生重建图像;所述方法沿欠采样方向应用所限定的加权变换(62)。
17.一种磁共振成像扫描仪(10,40),包括:
患者支撑(14),用于连续移动成像对象通过扫描仪视场;
多个射频线圈(24,26),设置为采集沿至少一个欠采样方向被欠采样的k-空间数据;以及
如权利要求16所述的处理器。
18.一种存储介质,其通过相关数字处理器(40)编码可执行指令以对在成像对象(12)连续移动期间使用多个射频线圈(24,26)所采集的k-空间数据执行图像数据处理方法并包括在至少一个欠采样方向上的k-空间欠采样,所述图像数据处理方法包括:沿连续移动方向混合变换所采集的k-空间数据,以限定具有连续移动变换方向上的真实空间尺寸和横过连续移动方向的横向方向上的k-空间尺寸的混合空间数据;对于该至少一个欠采样方向限定从所有射频线圈的k-空间至真实空间的加权变换(62),所述加权变换组合与患者位置相关的线圈灵敏度权重因子并包括沿横向方向的混合空间数据的傅立叶变换以产生重建图像;所述方法沿欠采样方向应用所限定的加权变换(62)。
CNA2006800112918A 2005-04-06 2006-03-27 床连续移动的灵敏度编码磁共振并行成像 Pending CN101351721A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US66889005P 2005-04-06 2005-04-06
US60/668,890 2005-04-06

Publications (1)

Publication Number Publication Date
CN101351721A true CN101351721A (zh) 2009-01-21

Family

ID=36613397

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2006800112918A Pending CN101351721A (zh) 2005-04-06 2006-03-27 床连续移动的灵敏度编码磁共振并行成像

Country Status (4)

Country Link
US (1) US20090003674A1 (zh)
EP (1) EP1869487A2 (zh)
CN (1) CN101351721A (zh)
WO (1) WO2006106448A2 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102551978A (zh) * 2012-01-19 2012-07-11 江苏美时医疗技术有限公司 一种磁共振兼容婴儿培养箱
CN102670201A (zh) * 2011-03-11 2012-09-19 西门子公司 用于产生磁共振图像数据的方法和磁共振断层造影系统
CN103764025A (zh) * 2011-10-12 2014-04-30 株式会社日立制作所 磁共振成像装置以及磁化率强调图像生成方法
CN102043137B (zh) * 2009-10-21 2015-03-04 凯斯西储大学 全程放射grappa校准
CN107209242A (zh) * 2014-11-07 2017-09-26 皇家飞利浦有限公司 用于生成在其环境中的移动的对象的mr图像的方法和系统
CN109069059A (zh) * 2016-04-22 2018-12-21 通用电气公司 用于对移动的主体成像的系统和方法
CN110471014A (zh) * 2018-05-10 2019-11-19 西门子医疗有限公司 磁共振成像中的条纹伪影减少
CN113466765A (zh) * 2020-03-31 2021-10-01 通用电气精准医疗有限责任公司 磁共振扫描方法及系统、计算机可读存储介质

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005047920A1 (en) * 2003-11-13 2005-05-26 Koninklijke Philips Electronics, N.V. Readout ordering in collection of radial magnetic resonance imaging data
DE102008047218B4 (de) * 2008-09-15 2010-08-19 Siemens Aktiengesellschaft Verfahren zur Erstellung von MR-Bildern mittels einer Magnetresonanzanlage sowie entsprechend ausgestaltete Magnetresonanzanlage
CN107422361A (zh) * 2010-01-22 2017-12-01 登克特有限公司 用于多照相机x射线平板检测器的方法和装置
JP6042069B2 (ja) * 2012-01-25 2016-12-14 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
DE102012217262B4 (de) * 2012-09-25 2023-10-26 Siemens Healthcare Gmbh Verfahren zur automatischen Magnetresonanz-Bildgebung und Magnetresonanzanlage
CN103777162B (zh) * 2012-10-19 2018-06-05 上海联影医疗科技有限公司 磁共振成像k空间运动伪影矫正并行采集重建方法
US20160029970A1 (en) * 2014-07-31 2016-02-04 Lg Cns Co., Ltd. Apparatus and method for determining validity of biological data
WO2020161204A1 (en) * 2019-02-06 2020-08-13 Koninklijke Philips N.V. Reconstruction and reuse of phase oversampled data in magnetic resonance spine imaging

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6289232B1 (en) * 1998-03-30 2001-09-11 Beth Israel Deaconess Medical Center, Inc. Coil array autocalibration MR imaging
GB2345139A (en) * 1998-12-24 2000-06-28 Marconi Electronic Syst Ltd MRI apparatus with continuous movement of patient
EP1218768A2 (en) * 1999-12-15 2002-07-03 Koninklijke Philips Electronics N.V. Magnetic resonance method for forming a fast dynamic image
US6946836B2 (en) * 2000-04-25 2005-09-20 Kabushiki Kaisha Toshiba Magnetic resonance imaging involving movement of patient's couch
US6564082B2 (en) * 2001-04-19 2003-05-13 General Electric Company Method for incremental field-of-view-MR imaging
US7009396B2 (en) * 2002-09-12 2006-03-07 General Electric Company Method and system for extended volume imaging using MRI with parallel reception
CN1977182A (zh) * 2004-04-28 2007-06-06 皇家飞利浦电子股份有限公司 涉及对比度处理和扫描参数更新的连续移动平台的磁共振成像

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102043137B (zh) * 2009-10-21 2015-03-04 凯斯西储大学 全程放射grappa校准
CN102670201B (zh) * 2011-03-11 2016-08-03 西门子公司 用于产生磁共振图像数据的方法和磁共振断层造影系统
CN102670201A (zh) * 2011-03-11 2012-09-19 西门子公司 用于产生磁共振图像数据的方法和磁共振断层造影系统
CN103764025A (zh) * 2011-10-12 2014-04-30 株式会社日立制作所 磁共振成像装置以及磁化率强调图像生成方法
CN103764025B (zh) * 2011-10-12 2015-12-09 株式会社日立制作所 磁共振成像装置以及磁化率强调图像生成方法
CN102551978B (zh) * 2012-01-19 2015-03-25 美时医疗技术(上海)有限公司 一种磁共振兼容婴儿培养箱
CN102551978A (zh) * 2012-01-19 2012-07-11 江苏美时医疗技术有限公司 一种磁共振兼容婴儿培养箱
CN107209242A (zh) * 2014-11-07 2017-09-26 皇家飞利浦有限公司 用于生成在其环境中的移动的对象的mr图像的方法和系统
CN107209242B (zh) * 2014-11-07 2020-07-03 皇家飞利浦有限公司 用于生成在其环境中的移动的对象的mr图像的方法和系统
CN109069059A (zh) * 2016-04-22 2018-12-21 通用电气公司 用于对移动的主体成像的系统和方法
CN110471014A (zh) * 2018-05-10 2019-11-19 西门子医疗有限公司 磁共振成像中的条纹伪影减少
US11435419B2 (en) 2018-05-10 2022-09-06 Siemens Healthcare Gmbh Streak artifact reduction in magnetic resonance imaging
CN113466765A (zh) * 2020-03-31 2021-10-01 通用电气精准医疗有限责任公司 磁共振扫描方法及系统、计算机可读存储介质

Also Published As

Publication number Publication date
US20090003674A1 (en) 2009-01-01
WO2006106448A3 (en) 2007-10-11
WO2006106448A2 (en) 2006-10-12
EP1869487A2 (en) 2007-12-26

Similar Documents

Publication Publication Date Title
CN101351721A (zh) 床连续移动的灵敏度编码磁共振并行成像
US7768263B2 (en) Magnetic resonance imaging apparatus and method
US6949928B2 (en) Parallel MR imaging using high-precision coil sensitivity map
US7860291B2 (en) Method and apparatus for correcting motion in multi-shot diffusion-weighted magnetic resonance imaging
JP5100374B2 (ja) 磁気共鳴イメージング装置および方法
US9018951B2 (en) Accelerated multispectral data magnetic resonance imaging system and method
CN101305908B (zh) 提高并行采集图像重建的信噪比损失的准确度的方法
JP3514320B2 (ja) 磁気共鳴イメ−ジング法
CN101171527A (zh) 磁共振成像系统的各个信号通道中的独立运动校正
CN101248458A (zh) 高约束图像重构方法
JP2006021049A (ja) 連続的テーブル移動中の可変視野からmriデータを取得する方法
CN101484823A (zh) 磁共振成像设备和方法
US6564082B2 (en) Method for incremental field-of-view-MR imaging
JP3556176B2 (ja) Mr投影画像の自動修正
EP1156343A1 (en) Method and apparatus for nuclear magnetic resonance imaging
CN105051563A (zh) 使用相位调制rf脉冲的并行多切片mr成像
JP5385499B2 (ja) 連続的テーブル移動により取得された磁気共鳴画像におけるアーチファクトの除去方法
US6915152B2 (en) Method for MR imaging with an array of RF coils
US20040150400A1 (en) Magnetic resonance method for forming a fast dynamic imaging
JP2006507071A (ja) 磁気共鳴方法
US7417429B2 (en) Moving table MRI with subsampling in parallel
US7436178B2 (en) Fast continuous moving bed magnetic resonance imaging with multiple stationary receive coils for sense reconstruction
JP4675936B2 (ja) 核磁気共鳴撮影装置
JP3983792B2 (ja) 核磁気共鳴撮影装置
JP3952310B2 (ja) 核磁気共鳴撮影装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20090121