CN107533120B - 磁共振成像图像翘曲之校正系统与方法 - Google Patents
磁共振成像图像翘曲之校正系统与方法 Download PDFInfo
- Publication number
- CN107533120B CN107533120B CN201580079293.XA CN201580079293A CN107533120B CN 107533120 B CN107533120 B CN 107533120B CN 201580079293 A CN201580079293 A CN 201580079293A CN 107533120 B CN107533120 B CN 107533120B
- Authority
- CN
- China
- Prior art keywords
- space
- image
- magnetic field
- gradient
- field
- 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
Images
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/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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56572—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field
-
- 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
-
- 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
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Nonlinear Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供了一种校正图像翘曲的方法,所述图像在MRI系统中获取,而所述翘曲由梯度线圈的梯度场分布中之非线性引起,此方法包括a)构建代表所述MRI系统中每个梯度线圈的导电路径的计算机模型;b)从所述模型中,在每所述梯度线圈每个空间点处计算预测磁场;c)在所述MRI系统中每个所述梯度线圈的每个空间点处测量真正实际的磁场;d)在每个所述空间点将所述预测磁场与所述真正实际磁场作比较,而验证所述模型的准确度,然、当所述模型不准确时重就重复a)至d),并且当所述模型不准确时则;根据所述预测磁场的线性偏离构建畸变图表,以便映射真实空间坐标在所述获取图像的翘曲空间之上;以及利用所述所畸变图表来返正所述所获图像的翘曲。
Description
发明背景:
本发明一般与磁共振成像有关。进一步说,本发明与由梯度线圈非线性引起的图像翘曲有关。
发明背景:
磁共振图像(MRI)是一医学上主要的造影技术。MRI能产生各软组织的详细图像,如脑脏、肌肉、及肾脏。利用组织中含有的化合物,如水分或脂肪,其特有的性质来创造图像。举例,在一个强力磁场影响之下,大量拥有核自旋角动量的原子,如水分或脂肪中所含之氢原子,此氢原子的总核磁矩矢量和会产生一个与上述磁场成一直线对准的净磁矩。这净磁矩会进一步进动于一个与磁场正比的确定频率。经由无线射频冲脉激励之后,这净磁化能便可提供一个可检测之信号。
此产生的信号可使用各种机构来进行编码,以允许转换空间图像。例,梯度线圈(高电压电磁体)利用编码空间信息。通过使梯度线圈产生线性变化的磁场来实现位置性编码,该磁场具有放置要被扫描的图像的成像体积中的位置。在现实世界实现中,由梯度线圈产生的场分布与纯线性有偏差。此偏差量取决于图像内的空间位置。与纯线性场的偏差导致获取图像的空间翘曲。因此,在获取成像数据的处理期间需要撤销图像翘曲或「返正」获取图像。
发明内容:
本发明提供了一个对MRI扫描系统和方法的创新系统和方法去消除或减缓至少一个以上确认之现有存在的现有技术缺点。
根据本说明书的一个方面,提出了一种校正MRI系统中的获取图像的翘曲的方法,该方法由梯度线圈的梯度场分布中的非线性引起,包括:
a)构建表示所述MRI系统中每个梯度线圈的导电路径的计算机模型;
b)在所述模型中计算每个所述梯度线圈的每个空间点处的预测磁场;
c)在所述MRI系统中针对每个所述梯度线圈测量空间中的每个点处的真正实际磁场;
d)通过将所述预测磁场与空间中的每个所述点处的所述真正实际磁场进行比较并且在所述模型不准确然后重复a)至d)的情况下验证所述模型的精确度,并且在所述模型是准确的情况下;
e)根据所述预测磁场与线性的偏差,构建用于将真实空间中的坐标映射到所述获取图像的翘曲空间中的坐标的畸变图;以及
f)使用所述畸变图表返回正取得获取图像的翘曲。
这些跟其他方域及其优点后来之可以彰显乃在于本发明之构造和操作详情,而此详情会随后完全地描述和声称。附图之引用至引用相应附图之说明部分亦构成此详情之一部分,其中相同的编号表示其相同的部件。
附图说明
图一示出一个根据实践形式的磁共振成像(MRI)系统的功能子系统之框图。
图2示出根据一种实施方式由图1MRI系统扫描而成的成像体积和其对应切片。
图3示出根据一实施方式示例的一个脉冲序列。
图4示出根据一种实施方式中包含一条接收线的k空间的示意图。
图5示出根据一实施方式的一个脉冲序列示例。
图6示出在成像体内扫描的物体。
图7示出可以由梯度线圈产生的磁场中的非线性产生的示例图像翘曲。
图8示出根据实施例的使用图1MRI系统校正空间翘曲的方法。
图9示出根据实施例的图像混叠的说明性示例。
图10示出根据实施例的体素塌陷的效果。
具体实施例
传统的磁共振图像(MRI)系统是一种成像模式,其主要是利用物体中如氢原子内质子发出的磁共振(MR)信号来构成图像。在医学MRI上,一般是关注MR信号来自组织中主要含氢的成份,如水份和脂肪。
现在参照图一,示于(100)是一个根据实践形式例子的磁共振成像(MRI)系统之框图。此MRI系统实践例子只插图之用,可以容许合各单元的变动,包括增加、减小、或改换其部件部分。
从图1示出,此显示的MRI系统(100)包括一个数据处理系统(105)。此数据处理系统(105)通常可以包括一个或多个输出仪器,如显示器,一个或多个输入设备例如键盘和鼠标,及连接到有易失和持久性组件存储器的一个或多个处理器。此数据处理系统(105)还可包括一个或多个用于执行扫描的介面,其适应与MRI系统(100)的硬件组件作通信及数据交换。
接续图一,此范例MRI系统(100)亦包括一个主磁力场(110)。这主场磁铁(110)可由一个如永久的、超导体的、或电阻性的磁铁去落实。其他类型的磁铁,包括适合这MRI系统(100)使用的混合型磁铁,现这可从一个技术人员实施与考虑。此主场磁铁(110)可被操作而产生一个实质均匀主磁场,而此磁场具有磁场强度B0以及沿主轴的方向。此主力磁场为用作产生一个图像体积。其中所要求的物体中的原子核心如处于水分和脂肪中的氢质子处于对齐磁场排列作扫描准备。在某实践过程中,如这个实践范例中,一个与数据处理系统(105)沟通的主场磁铁操纵单元(115)可作控制主场磁铁(110)之用。
此MRI系统(100)还包括梯度线圈(120)对沿着主磁场的空间信息进行编码,如三垂直梯度轴心。梯度线圈(120)的大小和配置可以令其产生一个可受控制而且均匀的线性梯度。例如,可以设计在主场磁铁(110)内的三对载运正交电流的主线圈去产生所期望的线性梯度磁场。
在某实施方案中,此梯度线圈(120)可被屏蔽并且可包括外层屏蔽线圈,其可产生一个反向磁场而抵消由主梯度线圈产生的梯度磁场,而作成一个主屏蔽线圈对。在这种线圈对中,其「主」线圈可负责产生梯度场,而其「屏蔽」可负责减小在某一体积外的主线圈杂散场,如成像体积。在此梯度线圈(120)的的主及屏蔽线圈可串联连接。而且,在任何一起成屏蔽梯度线圈的梯度轴中,也可有于两层线圈。此屏蔽梯度线圈(120)可减少能导致扫描图像中伪影之涡流和其它干扰。由于通常在此MRI系统(100)的导电部件中流动涡流,是由成像体积外的磁场(又称「边缘场」)而引起的,所以,降低这个从梯度线圈(120)产生的边缘场能一同减低干扰。因此,此主屏蔽线圈对之形状与大小、导线之图案与大小、及电流之幅度与图案可以作为选择,以促成一个在于梯度线圈(120)外部,而强度接近零的净磁场。如对于圆柱形磁体,其两个线圈能以由共同轴心定位,而对于垂直磁场磁体,其两个线圈能以由共同轴盘定位。
此梯度线圈(120)的导电组件,不论屏蔽或非屏蔽,并包括主线圈及屏蔽线圈,可以由一个电导体(如铜,铝等)而组成。当此梯度线圈(120)之端子受一个电压差施加时,其内部电连接可使电流在预想的路径上流动。此在三个梯度轴,用于主梯度线圈和梯度屏蔽线圈的导电部件可从物理分离与或非导电屏障而隔离。
由此梯度线圈(120)产生的磁场,不论齐发与或顺序,可以叠加在主磁场上,以得在成像体积内之物体受选择性空间激发。除了允许空间激发之外,此梯度线圈(120)可将空间特定的频率与周期资料附置在成像体积内的原子核,从而允由所得MR信号重建有用图像。一个梯度线圈操纵单元(125)与其数据处理系统(105)互相通信,以用来操作此梯度线圈(120)。
在某此MRI系统(100)之实施方案中,有可能带有额外电磁线圈(无在图上表出),如垫片线圈(通常包括但不限产生二阶或更高的球谐函数的磁场分布)、或一个均匀场偏移线圈、或其他校正电磁铁。为执行有效的匀场过程(即校正当不同物体放置在系统之内或周围时引入的磁场失真),如此校正电磁铁,如垫片线圈,具有用来提供磁场的一种电流,而令主磁场更均匀。举例,这由此线圈产生的磁场可助于校正主磁场中的不均匀性,其由于主场磁铁(110)中的不完善点、或外部铁磁物体、或由于成像基体内物质的磁化率差、或任何其他静态或时变的现象而引起。
此MRI系统(100)还包括射频(RF)线圈(130)。此RF线圈(130)用于建立具有强度B1的RF磁场,以其激发原子核或「旋转」。此RF线圈(130)亦可检测从被成像对象内旋转「弛豫」时发出的信号。因此,此RF线圈(130)可有分开发射和接收线圈之形式,或可有共合发射和接收线圈形式,其有发射和接收模式可用一个切换机构切换。
此RF线圈(130)可实现为表面线圈,其通常当为纯粹接收线圈,与或体积线圈,其可当为接收和发送线圈。此RF线圈(130)能与主场磁铁(110)之孔綜合。此外,此RF线圈(130)或可以设于靠近扫描对象的位置,例如头部,其并且可具有类似这对象的形状,例如一个紧密头盔。一个射频线圈操纵单元(135)与其数据处理系统(105)互相通信,以用来操作此RF线圈(130)。
有很多不同技术可以用此MRI系统(100)来获取图像,包括T1与T2加权图像。以下一个非限示例简述收获质子密度加权图像的进程,为提供一个本MRI系统(100)功能之简述。根据前述的示例简述来创建图像,当一个物体受一个比较大的磁场影响时,此MRI系统(100)检测一个物体中所存而含有旋转角动量核子之存在,如组织内水分和脂肪所有的氢质子。在该实施例中,此主磁场有强度B0并含有旋转角动量、类似氢质子的核子。该主磁场将在主场磁铁(110)成像中体积里的物体中的氢质子部分极化。然后,此质子会通过适当调整的RF辐射激发,形成如B1强度的RF磁场。最后,当质子从磁相互作用下「弛豫」时,来自此激发质子的弱RF辐射信号会检测为MR信号。而此检测到的MR信号频率(拉莫尔频率)与经受磁场成正比。
从其获得信号物体上可以通过在物体上产生磁场梯度来选择横截面,以使得主磁场的磁场值可沿物体中不同位置转变。因信号频率与此产生磁场的变化成正比,所以这变化允许将一个特定信号频率分配给相位分配给一个在此物体中的位置。因此,由此获得的MR信号中可以足够的信息,以用质子所存为物本体构成一张解剖图,此是传统MRI图像的基础。举例,由于质子密度随组织类型变化,此组织变化可以在所获得的信号被处理之后映射为图像对比度变化。
现在参考图2,为进一步说明此MRI系统(100)的信号获取过程示例,假设一个放置在一个主磁铁110内一个成像体积250之中的物体,其磁铁具有一个强度B0而指向z轴(指示在240)的磁场。该物体随后具有净磁化矢量。在该说明性示例中,一个成像会在一个沿x和y轴(指示在205)平面切片取获。应当注意,在该示例中,兆切片之z轴方向有一个有限的厚度,从而产生体积切片(205)。
为以传统的方式从MRI系统(100)获得图像,一组或多组RF脉冲和梯度波形(统称为「脉冲序列」)会在数据处理系统105处选择。此数据处理系统(105)将所选择的脉冲序列信息传递给其RF控制单元(135)与其梯度控制单元(125),其共同为提供执行扫描发出相关的波形和时恰。
现参考图3,其说明可以用MRI系统(100)来获取图像的脉冲序列(300)。这具体展示了该脉冲序列例子的时序图。时间图示出以时间为函数的脉冲或信号幅度,以于发射(RFt)信号,磁场梯度Gx,Gy和Gz,接收RFx信号和归一化移位信号(RFt)。为了说明,一个简化的理想脉冲序列可以包藏着在RFt处射频脉冲(310)的选择切片,在Gz处梯度脉冲(320)的选择切片,在Gy处梯度脉冲(330)的相位编码,在Gx处脉冲(340)的频率编码以及在RFx处检测到的MR信号(350)。三个梯度脉冲Gx,Gy和Gz表现可由线圈(120)所产生的磁场梯度的幅度和持续时间。脉冲(310)之切片选择可由RF线圈(130)在发射层面上而产生。其检测到的MR信号(350)可以通过RF线圈(130)在接收层面而得到检测.在该说明示例中,其假设RF线圈(130)的发射层面和接收层面各由不同的线圈形成。
第一在脉冲序列(300)中可产生是打开切片选择梯度脉冲(320).在同一时间可应用切片选择RF脉冲(310)。在这说明示例中,切片选择RF脉冲(310)可以是呈正弦函数形状的RF能量爆破。在其他实现中,可以使用其他RF脉冲形状和持续时间。当一旦关闭切片选择RF脉冲(310),切片选择梯度脉冲(320)亦可被关闭,而相位编码梯度脉冲(330)则被为打开。在一些实施方案中,场相转移信号(360)也可以在这一点被接通以改变主磁场强度。一旦相位编码梯度脉冲(330)关闭,频率编码梯度脉冲(340)就可以被接通,并且可以记录检测到的MR信号(350)。应该注意的是,如图4中所示的脉冲和信号的形状,幅度,顺序和持续时间只作说明之用,而在其实现过程中,可通过改变一个或多个其因素或其它以实现要求的扫描结果。
在某变化中,脉冲序列(400)可以作一定次数的重复或迭代,例如256次,以收集产生一个图像所需的所有数据。每个重复通常涉及脉冲序列的变化,以允许接收对应于图像的不同分量的信号。相隔其脉冲序列(400)每次重复之时间可被称为重复时间(TR)。此外,切片选择脉冲(410)的中心点与其检测到MR信号(450)的峰值之间的持续时间可被称为回波时间(TE)。TR和TE两者都可为所需扫描作变化。
为了进一步说明MRI系统(100)的信号采集过程,参照图2并图3。为了选择切片,可以沿z轴施加切片选择梯度脉冲(320),满足位于切片(205)中的质子的共振条件。实际上,切片沿z轴的位置可以部分基于切片选择梯度脉冲(320)而决定。因此,在该示例中,与切片选择梯度脉冲(320)同时生成的切片选择脉冲(310)可以激发位于切片(205)内的质子。位于切片(205)上方和下方的质子通常是不受切片选择脉冲(310)的影响。
继续说明性示例,根据脉冲序列(300),可以在切片选择梯度脉冲(320)之后施加相位编码梯度脉冲(330)。假设这是沿着y轴施加,沿着y轴的不同位置的自旋可以按不同的拉莫尔频率开始进行。当相位编码梯度脉冲(320)关断时,处于不同位置的净磁化矢量可以按相同但具不同的相位的速率进行动作。相位可以由相位编码梯度脉冲(330)的持续时间和幅度而确定。
一旦相位编码梯度脉冲(330)关闭,频率编码梯度脉冲(340)就可以被打通。在这示例中,频率编码梯度是处于X方向上的。频率编码梯度可以使所选切片中质子跟据其在x中位置之速率作处理。因此,现切片内的不同空间位置是由其独特的相位角和进动频率来作表征。当频率编码梯度脉冲(340)接通时,RF接收线圈(130)可用于接收由正在被扫描物体中质子所产生的检测信号(350)。
当MRI系统(100)执行脉冲序列(300)时,所获取的信号可被存储在被称为K空间的临时矩阵中,如图4中(400)所示。通常,K空间是来收集扫描中所测量的信号,并且处于空间频域。可以通过沿x轴(420)(Kx)的频率编码数据和沿着y轴(430)(Ky)的相位编码数据来覆盖K空间。当接收到全部切片的K空间矩阵的行线时(例如,在单个切片的扫描结束时)数据可用数学处理去产生最终图象,例如通过二维傅立叶变换(傅立叶变换)。因此,在重建图像成空间领域之前,K空间是可保存着原始数据。通常,K空间是具有与最终图像相同的行数和列数,并且是填充着在扫描期中的原始数据,通常每个脉冲序列(300)是一行。例如,在(410)处所表示中,K空间(400)的第一行是在完成其扫描片段所生成的脉冲序列之第一次迭代后被填入,并包含用于该脉冲序列迭代的检测信号。在脉冲序列多次迭代后,K空间可被充填。脉冲序列的每个迭代可以作稍微改变,从而获得用于K空间适当部分的信号。应当注意的是,基于不同的脉冲序列,其它充填K空间的方法是可行的是被为考虑的,例如以螺旋方式。
通常由梯度线圈(120)产生的梯度场分布包括与线性场轮廓的偏差。作为说明性的示例,理想地,如图5(a)所示,当通过梯度线圈(120)沿X轴施加梯度场时,由B表示的主磁场的强度沿X轴成像在空间中是线性的并且具有Gx的斜率,Gx是梯度线圈(120)沿着X轴产生的期望的梯度。如图5(b)所示,然而,在(520)处,磁场强度的真正实际轮廓偏离期望的梯度Gx。在图5(b)中,在(530)处用虚线表示理想梯度Gx作为参考。
梯度场分布与线性轮廓的偏差通常导致获得空间翘曲的图像。参考图6,再次示出图2的成像体积(250)。然而,在图6中,物体O已放置在成像体积(250)中,并且根据例如如上所述的脉冲序列,使用MRI系统(100)获取切片(605)的图像。物体O推定为球形壳体,其被选择用于说明性目的。因此,在该说明性示例中,当沿着切片(605)扫描时,物体O应该导致包含圆圈(705)的图像,如图7(a)所示,对应于包含在球形壳体物体O中的圆周切片(605)。然而,如图7(b)中(710)所示,梯度线圈(120)产生的梯度场中的非线性导致收取图像的翘曲。图7(b)所示的翘曲的程度和形状是简单的选择,仅是清楚地说明了翘曲效果。
根据本发明的一方面,此获取图像可以处理以校正由梯度场分布中的非线性引起的空间翘曲。现在参见图8,一种校正空间翘曲的方法一般在(800)表示。为了帮助解释该方法,将假定方法(800)部分地使用如图1所示的MRI系统(100)来操作。此外,对方法(800)的以下讨论导致对系统(100)的进一步理解。然而,应当理解,系统(100)和方法(800)可以变化,并且不需要如本文中所讨论的那样精确地工作,并且这些变化在所附权利要求书的范围之内。
为了校正图像的空间翘曲,一个梯度线圈(120)的计算机模型在(810)创建。该模型表示梯度线圈(120)中存在的单元,例如存在于MRI系统(100)中的线圈的线模式,如包括电线连接。例如,该模型可以包括主线圈的导电单元,并且在变化中可以包括护罩线圈。该模型通常构造为尽可能准确地基于计算机建模技术的梯度线圈(120)的表示。一个示例模型基于单元阵列,即电磁体被建模为一组小电流单元,当连接在一起时,形成电磁体的线模式。此单元阵列模型可以包括梯度线圈(120)的连接路径以及存在于梯度线圈(120)中的任何其它载流线。此基于此单元阵列模型可以构建梯度线圈(120)的完整的空间表示。此单元阵列模型可以包括例如表示梯度线圈(120)的线模式的导电单元的阵列。
继续方法(800),用于梯度线圈(120)的每个所需空间点的预测磁场是基于该模型在(820)中计算。预测磁场的计算可以基于用于计算来自磁体模型的预测磁场的各种方法。例如,基于单元阵列模型,可以使用适当形式的毕奥-萨伐尔方程来计算成像体积(250)中每个期望点处的预测磁场。根据一个实现,毕奥-萨伐尔方程可以采取以下形式:
该公式用于计算由稳定电流I(例如由于单元阵列模型的一个导电单元)产生的位置处的合成磁场强度BP。是从导电单元到计算场的点的全位移向量,r是其大小。dl是传统电流方向的幅度是导电单元的微分单元的长度的向量,μ0是自由空间磁导率的磁常数。
基于毕奥-萨伐尔方程的计算可以适当地进行,以计算成像体积中每个点处的预测磁场强度。在变化中可以使用毕奥-萨伐尔方程的不同公式或其他等效方程式来计算基于梯度场线圈(120)模型的预测磁场。例如,当电流可以近似为通过无限窄的电线运行时,此上面公式可良好地运行。如果导线具有一定的厚度,则可以相应地改变毕奥-萨伐尔方程公式。尽管在该说明性示例中使用毕奥-萨伐尔方程,但在变化中可以使用等效等式或公式来计算空间中每个点处的预测磁场。
继续图8,在(830),通过测量由真正实际构造的梯度线圈产生的磁场来验证预测磁场的准确度。例如,此测量可以使用3轴霍尔探测器或其他类似的现场测量装置进行。将收取的测量值与预测的场值进行比较。当被测点和空间中的被测点的预测场值被发现是一致的,例如在预定阈值内时,该方法继续到(840)。然而,当发现空间中被测试点的测量值和预测场值不一致时,模型在(810)处被适当地更新以提高模型的准确度。在变型中,还可以更新所使用的配方和方程,以提高预测磁场的准确度。
在(840)中,畸变图表(映射函数)是基于预测的磁场构建,通过创建可以在需要之间插值或通过计算各个所须数值的映射。例如,可以使用以下等式来将真实空间中的坐标,例如,如图6所示,物体O空间映射到基于梯度场的获取图像的翘曲空间中的点从线性的偏差:
因此,物体的空间中的点(x,y,z)被映射到图像的点(x',y',z')。B(x,y,z)值是从每个梯度线圈的毕奥-萨伐尔方程中获得的。Bx是沿X轴产生的梯度线圈(120)的预测磁场强度的z分量(或与主静态共轴的分量),By是梯度线圈的预测磁场强度的z分量(120),Bz是沿Z轴产生的梯度线圈(120)的预测磁场强度的z分量。例如,参考图7,有效的畸变图表示图7(a)的与物体空间对应的圆衫如何映射到图7(b)的椭圆形。因此,使用图中的点(x,y)等式如图7(b)所示,例如,图7(a)是决定来映射到点(x',y')。
一旦畸变图表被收取在(840),在畸变空间出现扫描物体的获取图像可以通过基于畸变图表的插值在(850)返正或畸变。
munwarped(x,y)=w(x′,y′)mwarped(x′,y′)
为了确定在真实空间x'和y'中的位置(x,y),munwarped(x,y)处的返正图像的值,使用上述用于构造畸变图表。场畸变图表可以与切片选择脉冲的描述一起使用,以按照上述等式的评估来确定每个(x,y)位置处的激发z位置。此翘曲图像在(x',y')的数值,mwarped(x′,y′)是从对应于mwarped的抽样的相邻(x',y')位置插值来计算,如从傅里叶变换获得k空间数据。w(x^',y^')是在每个(x',y')位置处施加的加权值,以校正由于体素体积在成像体积上的变化而引起的信号堆积。一种计算w(x′,y′)的方式是从翘曲空间体素体积在(x',y')与翘曲空间体素体积在(x,y)的比例之中。
实际上,通过返正,收缩图像被拉伸或以其他方式变换,使得图像中的物体被转换成与原始物体空间更紧密对应的空间。例如,通过内插,图7(b)的翘曲图像可以变换成更圆形的形状,其大小更接近于球形物体O的真正实际圆形横截面,如图7(a)所示。
失真函数可以在图像空间中实现,如使用最近邻,线性或立方插值方法(或其他相等方式)。除了观察此返正问题之外,还可以从翘曲图像空间中的已知值的均匀抽样网格内插值,该问题也可以直接从k空间数据看作(x',y')的计算值:
其中Mwarped(k′x,k′y)代表k空间的数据抽样集合。这难题可以分别为是一个反网格化问题,或,当抽样从均匀到不均位置时,一个不均匀的傅里叶变换。因此,可以使用逆网格化方法来有效地求解,例如通过在k空间中执行变迹修正,随后进行反傅立叶变换,然后使用使用凯泽-贝塞尔核在图像空间内卷积插值。此外,在此反网格校正翘曲法之下,傅里叶变换的大小是由k空间矩阵大小和超抽样因子决定。比蒂(Beatty)的方法指出,高准确度结果可以用很小超抽样因子收集,如1.25至1.5。
如上所述,MR图像通常被重建为比获取矩阵大小更大的矩阵大小。例如,一个320×192采集矩阵可能被重建为一个512×512图像。这通常在k空间矩阵中通过零填充到512×512来执行傅里叶变换来,随后通过图像空间插值进行翘曲校正。然而,这种操作可以在1.25倍以上的抽样网格上更有效地执行,使得能够进行大小为400×240的逆傅立叶变换,随后进行插值(例如,凯泽-贝塞尔核)在一个512×512图像矩阵。由于400×240矩阵比512×512矩阵至小2.7倍,因此可提高计算机和计算机存储器的利用率,特别是对于较大的零填充因子。另外,认识到对于涉及多个接收机通道的重构,可以在400×240矩阵上的信道上组合数据,并且仅对512×512矩阵进行插值仅需要对每个重建图像执行一次。这将进一步提高效率。
在一些变型中,混叠伪影可能发生在翘曲空间中物体的部分映射到图像的视野(FOV)边界之外。例如,图像中的物体可以翘曲超过图像的边界,并返以一定方式回到图像中,因此不能有实际地进行返正。当被成像的物体接近图像的视野(FOV)的大小取决于与当存在大的梯度场非线性时,混叠伪影变得重要。如果翘曲是使物体的一部分翘曲地驻留在FOV之外,则这可能导致不能容易地去除的混叠伪影。参考图9,(900)示出了如何在存在大梯度场非线性的情况下物体O的扫描式切片(605)可以导致混叠图像的简化说明性示例。
图像的FOV可以放大以避免影像混叠。因此,在一个变型中可以基于畸变图表计算给定FOV的预期失真的最大量。然后可以使用计算的失真来创建放大的FOV。例如,在图像的x值从0变到255的情况下,可以基于畸变图来计算x=0和x=255到其翘曲位置的变换。这些转换可以指示视野之外(即低于0及以上255),X轴的变形可以延伸。因此,可以增加FOV以适应变换。类似地,物体可以被扭曲成小于预期的FOV。这将降低图像的分辨率。因此,可以减小FOV以提高图像分辨率。一般来说,畸变图可以用于在真实空间中找到包含规定FOV的翘曲空间中的最小FOV。因此,翘曲空间中产生的FOV可用于创建脉冲序列。在许多情况下,有利的是使用户能够在真实空间中开出非矩形FOV(例如椭圆或圆角角矩形),以便在翘曲空间中实现较小的包围FOV。
基于放大的FOV的确定,物体的图像获得,使获得图像具有较大的FOV。例如,当获取真正实际物体的图像时,k空间可以更精细地抽取。例如,物体的k空间收取可以在512x512的分辨率代替256×256的程度(即相同的图像分辨率)。在返正过程中,图像可以返回给使用者确定FOV。例如,一旦返正图像被收取,它可以根据期望的FOV进行裁剪。
人工增加的图像分辨率也可以帮助返回正梯度不均匀性的能力,而不会由于许多体素在翘曲图像中的单个体素上的崩溃而导致伪影。例如,获得的翘曲图像可以增加到更高的抽样密度,例如通过零填充翘曲图像的k空间。参考图10,体素折叠在图10(a)中示出。可以看出,两个体素的翘曲和返正导致两种体素之间的强度差异的丧失。另一方面,如图10(b)所示,分辨率提高可以减少或避免损失。当在空间频域中执行返回函数时,该过程可以是特别有效的。这是因为高抽样图像的形成可以一起避免,节省计算时间和存储器;在卷积之前,只能向获取k空间的外部添加零。此外,在这个阶段,直接在零填充的k空间中附加一个变迹函数,以便在去除吉布斯振荡伪影时加以应用。
上述实施例旨在作为示例,所属技术领域的技术人员可以对其进行变更和修改,而不脱离仅由所附权利要求书定义的范围。例如,所讨论的方法,系统和实施例可以全部或部分地变动和组合。
Claims (16)
1.一种校正图像翘曲的方法,所述图像在MRI系统中获取,而所述翘曲由梯度线圈的梯度场分布中之非线性引起,此方法包括:
a)构建代表所述MRI系统中每个梯度线圈的导电路径的计算机模型;
b)从所述模型中,在每个所述梯度线圈每个空间点处计算预测磁场;
c)在所述MRI系统中每个所述梯度线圈的每个空间点处测量真正实际的磁场;
d)在每个所述空间点将所述预测磁场与所述真正实际磁场作比较,而验证所述模型的准确度,并且当所述模型不准确时就重复a)至d),并且当所述模型准确时则;
e)根据所述预测磁场的线性偏离构建畸变图表,以便映射真实空间中的坐标到所述获取图像的翘曲空间中的坐标;以及
f)利用所述畸变图表来返正所述获取图像的翘曲,
其中所述预测磁场是从毕奥-萨伐尔或相等方程中之一计算,
通过为每个所述梯度线圈从毕奥-萨伐尔或相等方程取得B(x,y,z)值来构造所述畸变图表,所述畸变图表用于映射在真实空间中每一点(x,y,z)到在图像翘曲空间中的图像的一点(x',y',z'),
并且其中(x',y')的数值是由k空间数据使用从均匀到非均匀抽样位置的反网格来直接计算,所述反网格是通过在k空间中执行变迹校正,随后进行反傅里叶变换,然后用凯泽-贝塞尔核在图像空间中进行卷积插值来实现。
2.根据权利要求1所述的方法,其中所述测量是使用场测量设备进行。
3.根据权利要求2所述的方法,其中所述场测量设备包括3轴霍尔探测器。
4.根据权利要求1所述的方法,其中所述计算机模型是梯度线圈电磁体单元阵列模型,其包括一组连接一起的单元,以为每个所述梯度线圈形成线模式的空间表示。
5.根据权利要求1所述的方法,其中所述毕奥-萨伐尔或相等方程包括:
并且其中Bx代表与预测磁场主静态场共轴的向量,其沿着所述MRI系统x轴产生,其中By代表与预测磁场主静态场共轴的向量,其沿着所述MRI系统y 轴产生,并且其中Bz代表与预测磁场主静态场共轴的向量,其沿着所述MRI系统z轴产生,并且其中Gx代表梯度线圈沿着x轴产生的所需要的梯度,Gy代表梯度线圈沿着y轴产生的所需要的梯度,Gz代表梯度线圈沿着z轴产生的所需要的梯度,并且fx(x,y,z)是真实空间映射到翘曲空间的x坐标畸变图表,fy(x,y,z)是真实空间映射到翘曲空间的y坐标畸变图表,fz(x,y,z)是真实空间映射到翘曲空间的z坐标畸变图表。
6.根据权利要求1所述的方法,其中所述返正是利用最近邻,线性或立方插值其中之一在图像空间中执行。
7.根据权利要求1所述的方法,其中所述傅里叶变换的大小由所述k空间矩阵大小和过抽样因子决定。
8.根据权利要求7所述的方法,其中所述凯泽-贝塞尔核是在插值期间被预先存储和重新抽样。
9.根据权利要求1所述的方法,还包括调整所述获取图像的视野,其通过或者放大或者缩小此视场,以或者消除混叠伪影或者提高图像分辨率,所述混叠伪影在翘曲空间中部分图像映射到真实空间图像视野边界之外而发生。
10.根据权利要求9所述的方法,其中所述放大包括:
使用所述畸变图表计算所述视场的期待最大畸变量;以及
使用所计的期待最大畸变量来创建放大的视野。
11.根据权利要求9所述的方法,还包括基于所调整的视场来产生用于所述MRI系统的脉冲序列。
12.根据权利要求10所述的方法,其中所述图像是根据所述放大视场通过对所述k空间进行过抽样来获得。
13.根据权利要求1所述的方法,还包括提高所述获取图像的抽样密度,以避免图像伪影,所述伪影是从多个体素折叠到所述翘曲空间中的单个体素而导致。
14.根据权利要求13所述的方法,其中所述获取图像的采样密度通过对所述k空间进行零填充而增加。
15.根据权利要求13所述的方法,进一步包括
在所述卷积之前将零添加到所述k空间外部之一部分。
16.根据权利要求15所述的方法,还包括应用变迹函数於零填充的k空间,以便消除吉布斯振荡伪影。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/IB2015/053059 WO2016174497A1 (en) | 2015-04-27 | 2015-04-27 | System and method for image warp correction for magnetic resonance imaging |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107533120A CN107533120A (zh) | 2018-01-02 |
CN107533120B true CN107533120B (zh) | 2020-03-24 |
Family
ID=57148602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580079293.XA Active CN107533120B (zh) | 2015-04-27 | 2015-04-27 | 磁共振成像图像翘曲之校正系统与方法 |
Country Status (7)
Country | Link |
---|---|
US (2) | US9529068B2 (zh) |
JP (1) | JP6499774B2 (zh) |
CN (1) | CN107533120B (zh) |
CA (1) | CA2982047C (zh) |
DE (1) | DE112015006495T5 (zh) |
GB (1) | GB2568540B (zh) |
WO (1) | WO2016174497A1 (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016174497A1 (en) * | 2015-04-27 | 2016-11-03 | Synaptive Medical (Barbados) Inc. | System and method for image warp correction for magnetic resonance imaging |
US20170192072A1 (en) * | 2016-01-06 | 2017-07-06 | General Electric Company | Methods and systems for correcting k-space trajectories |
DE102016218536A1 (de) * | 2016-09-27 | 2018-03-29 | Siemens Healthcare Gmbh | Verbesserte Verzeichniskorrektur von mit Hilfe von Magnetresonanztechnik aufgenommenen Messdaten eines Untersuchungsobjektes |
US10545211B2 (en) | 2017-06-28 | 2020-01-28 | Synaptive Medical (Barbados) Inc. | Method of correcting gradient nonuniformity in gradient motion sensitive imaging applications |
CN110095741B (zh) * | 2019-04-30 | 2021-09-28 | 上海东软医疗科技有限公司 | 磁共振线圈位置的确定方法、装置及磁共振成像系统 |
DE102019214669A1 (de) * | 2019-09-25 | 2021-03-25 | Robert Bosch Gmbh | Verfahren und Vorrichtung zum Betreiben eines mobilen Agenten mit einem Magnetometer |
DE102019129541A1 (de) * | 2019-10-31 | 2021-05-06 | Otto-Von-Guericke-Universität Magdeburg | Verfahren zum Betrieb eines bildgebenden Untersuchungsgeräts, Computerprogramm und Untersuchungsgerät |
CN112365760A (zh) * | 2020-10-29 | 2021-02-12 | 中南民族大学 | 一种磁共振成像设备和人体扫描的虚拟仿真系统及方法 |
EP4102244B1 (de) | 2021-06-09 | 2024-01-17 | Siemens Healthcare GmbH | Computerimplementiertes verfahren zur einstellung eines sichtfelds für eine magnetresonanzaufnahme, magnetresonanzeinrichtung, computerprogramm und elektronisch lesbarer datenträger |
CN114236439A (zh) * | 2021-11-09 | 2022-03-25 | 深圳市联影高端医疗装备创新研究院 | 线圈定位方法、装置、磁共振设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498774A (zh) * | 2008-01-31 | 2009-08-05 | 西门子公司 | 医学磁共振成像回波平面成像测量中的动态失真校正方法 |
CN101939150A (zh) * | 2008-02-06 | 2011-01-05 | 株式会社村田制作所 | 注塑成形件以及磁传感器 |
CN103376433A (zh) * | 2012-04-27 | 2013-10-30 | 西门子(深圳)磁共振有限公司 | 图像畸变校正方法及系统、磁共振成像设备 |
CN103885017A (zh) * | 2014-02-20 | 2014-06-25 | 厦门大学 | 基于单扫描正交时空编码磁共振成像的图像畸变校正方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3187092B2 (ja) * | 1990-09-20 | 2001-07-11 | シーメンス アクチエンゲゼルシヤフト | 核スピントモグラフ |
US6278275B1 (en) * | 1999-10-18 | 2001-08-21 | Picker International, Inc. | Gradient coil set with non-zero first gradient field vector derivative |
US6707300B2 (en) * | 2002-05-17 | 2004-03-16 | Ge Medical Systems Global Technology Co., Llc | Gradient non-linearity compensation in moving table MRI |
US6969991B2 (en) | 2002-12-11 | 2005-11-29 | The Board Of Trustees Of The Leland Stanford Junior University | Correction of the effect of spatial gradient field distortions in diffusion-weighted imaging |
EP1582886B1 (en) * | 2004-04-02 | 2012-08-15 | Universität Zürich | Magnetic resonance apparatus with magnetic field detectors |
DE102005051021A1 (de) * | 2005-10-25 | 2007-04-26 | Universitätsklinikum Freiburg | Apparaturen und Verfahren zur kernspintomographischen Aufnahme mit lokalen Magnetfeldgradienten in Verbindung mit lokalen Empfangsspulen |
US7408345B2 (en) * | 2006-02-06 | 2008-08-05 | The Board Of Trustees Of The Leland Stanford Junior University | Generalized MRI reconstruction with correction for multiple image distortion |
US7348776B1 (en) * | 2006-09-01 | 2008-03-25 | The Board Of Trustees Of The Leland Stanford Junior University | Motion corrected magnetic resonance imaging |
JP5683984B2 (ja) * | 2011-02-03 | 2015-03-11 | 株式会社日立メディコ | 磁気共鳴イメージング装置および非線形性歪み補正方法 |
JP6109508B2 (ja) * | 2012-08-29 | 2017-04-05 | 東芝メディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
WO2016174497A1 (en) * | 2015-04-27 | 2016-11-03 | Synaptive Medical (Barbados) Inc. | System and method for image warp correction for magnetic resonance imaging |
-
2015
- 2015-04-27 WO PCT/IB2015/053059 patent/WO2016174497A1/en active Application Filing
- 2015-04-27 JP JP2017554559A patent/JP6499774B2/ja not_active Expired - Fee Related
- 2015-04-27 GB GB1719260.0A patent/GB2568540B/en active Active
- 2015-04-27 CN CN201580079293.XA patent/CN107533120B/zh active Active
- 2015-04-27 CA CA2982047A patent/CA2982047C/en active Active
- 2015-04-27 US US14/903,636 patent/US9529068B2/en active Active
- 2015-04-27 DE DE112015006495.8T patent/DE112015006495T5/de active Pending
-
2016
- 2016-11-17 US US15/354,287 patent/US9989615B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498774A (zh) * | 2008-01-31 | 2009-08-05 | 西门子公司 | 医学磁共振成像回波平面成像测量中的动态失真校正方法 |
CN101939150A (zh) * | 2008-02-06 | 2011-01-05 | 株式会社村田制作所 | 注塑成形件以及磁传感器 |
CN103376433A (zh) * | 2012-04-27 | 2013-10-30 | 西门子(深圳)磁共振有限公司 | 图像畸变校正方法及系统、磁共振成像设备 |
CN103885017A (zh) * | 2014-02-20 | 2014-06-25 | 厦门大学 | 基于单扫描正交时空编码磁共振成像的图像畸变校正方法 |
Non-Patent Citations (1)
Title |
---|
A cradle-shaped gradient coil to expand the clear-bore width of an animal MRI scanner;K M Gilbert;《PHYSICS IN MEDICINE AND BIOLOGY》;20091221;第502-503页 * |
Also Published As
Publication number | Publication date |
---|---|
JP2018514277A (ja) | 2018-06-07 |
US20170067979A1 (en) | 2017-03-09 |
WO2016174497A1 (en) | 2016-11-03 |
CA2982047C (en) | 2019-04-02 |
CA2982047A1 (en) | 2016-11-03 |
US20160313434A1 (en) | 2016-10-27 |
CN107533120A (zh) | 2018-01-02 |
DE112015006495T5 (de) | 2018-01-18 |
GB201719260D0 (en) | 2018-01-03 |
GB2568540B (en) | 2021-10-13 |
US9529068B2 (en) | 2016-12-27 |
GB2568540A (en) | 2019-05-22 |
US9989615B2 (en) | 2018-06-05 |
JP6499774B2 (ja) | 2019-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107533120B (zh) | 磁共振成像图像翘曲之校正系统与方法 | |
US10401456B2 (en) | Parallel MR imaging with Nyquist ghost correction for EPI | |
CN107110922B (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 | |
EP3322997A1 (en) | Mr imaging with motion detection | |
JP2015128585A (ja) | 磁気共鳴生データの再構成方法および装置 | |
CN109143135B (zh) | 在梯度运动传感成像用途中校正梯度不均匀性的方法 | |
US10338180B2 (en) | System and method for gradient measurement using single-point imaging | |
CN106796274B (zh) | 具有伪迹抑制的propeller-mr成像 | |
CN110720047B (zh) | 双回波Dixon型水/脂肪分离MR成像 | |
US10996305B2 (en) | Propeller MR imaging with relaxation weightings varying between k-space blades and k-space data lines | |
US11914016B2 (en) | Optimized k-space profile ordering for 3D radial MR imaging | |
CN113544526B (zh) | 使用波编码的并行mr成像 | |
US11815577B2 (en) | Parallel MR imaging using wave-encoding | |
JP2006512992A (ja) | 非線形磁場勾配を用いる磁気共鳴法 | |
Ullisch | A navigator-based rigid body motion correction for magnetic resonance imaging | |
US20230400545A1 (en) | Spin echo mr imaging with spiral acquisition | |
Sumpf | Model-based T2 Relaxometry using Undersampled Magnetic Resonance Imaging |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201228 Address after: Toronto, Ontario, Canada Patentee after: SANAP medical Co. Address before: Bridgeton, Barbados Patentee before: SYNAPTIVE MEDICAL (BARBADOS) Inc. |