CN105676155B - 磁共振成像系统的磁场不均匀性值获取方法和装置和失真校正方法和装置 - Google Patents
磁共振成像系统的磁场不均匀性值获取方法和装置和失真校正方法和装置 Download PDFInfo
- Publication number
- CN105676155B CN105676155B CN201410663071.XA CN201410663071A CN105676155B CN 105676155 B CN105676155 B CN 105676155B CN 201410663071 A CN201410663071 A CN 201410663071A CN 105676155 B CN105676155 B CN 105676155B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- magnetic
- pixel
- imaging system
- resonance imaging
- 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
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- 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/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
-
- 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/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/243—Spatial mapping of the polarizing 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/387—Compensation of inhomogeneities
- G01R33/3875—Compensation of inhomogeneities using correction coil assemblies, e.g. active shimming
-
- 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/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Medical Informatics (AREA)
- Quality & Reliability (AREA)
Abstract
本发明公开了一种磁共振成像系统的磁场不均匀性值获取方法和失真校正方法。所述磁场不均匀性值获取方法,包括如下步骤:为了将磁共振成像系统的原始磁场均匀化为一目标磁场,针对磁共振成像系统通过动态匀场方法提供一磁场补偿量,其中动态匀场方法包括进行一三维低分辨率双回波梯度回波序列;利用一第三通式获取磁场不均匀性值,第三通式是:ΔB=ΔBoriginal+ΔBcompensatins其中,ΔB是磁场不均匀性值,ΔBariginal是原始磁场与目标磁场之间的差值,ΔBcompensating是磁场补偿量。根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值方法节省了大量用于再次映射磁场的时间,进而缩短了磁共振成像的时间,提高了磁共振成像的效率,所以这项技术对于医院大量繁重的扫描工作而言非常宝贵。
Description
技术领域
本发明涉及磁共振成像技术领域,特别是磁共振成像系统的磁场不均匀性值获取方法和失真校正方法。
背景技术
磁共振成像(Magnetic Resonance Imaging,MRI)是利用磁共振现象进行成像的一种技术。磁共振现象的原理主要包括:包含单数质子的原子核,例如人体内广泛存在的氢原子核,其质子具有自旋运动,犹如一个小磁体,并且这些小磁体的自旋轴没有一定的规律,如果施加外在磁场,这些小磁体将按外在磁场的磁力线重新排列,具体为在平行于或反平行于外在磁场磁力线的两个方向排列,将上述平行于外在磁场磁力线的方向称为正纵向轴,将上述反平行于外在磁场磁力线的方向称为负纵向轴;原子核只具有纵向磁化分量,该纵向磁化分量既具有方向又具有幅度。用特定频率的射频(RadioFrequency,RF)脉冲激发处于外在磁场中的原子核,使这些原子核的自旋轴偏离正纵向轴或负纵向轴,产生共振,这就是磁共振现象。上述被激发的原子核的自旋轴偏离正纵向轴或负纵向轴之后,该原子核就具有了横向磁化分量。
停止发射射频脉冲后,被激发的原子核发射回波信号,将吸收的能量逐步以电磁波的形式释放出来,其相位和能级都恢复到激发前的状态,将原子核发射的回波信号经过空间编码等进一步处理即可重建图像。
回波平面成像(Echo Planar Imaging,EPI)是一种存在几何失真的快速磁共振(MR)成像方法。为了对这种几何失真进行校正,现有技术的磁共振成像系统已经提出若干失真校正方法。这些方法并且已经广泛用于(具体来说)在EPI以及血氧水平依赖功能磁共振成像(Blood Oxygen Level Dependent Functional Magnetic Resonance Imaging,BOLDfMRI)中的失真校正;然而,这些方法依赖于预先获得的磁场图或点扩散函数。具体而言,在例如BOLD fMRI或弥散张量成像(Diffusion Tensor Imaging,DTI)等具有长测量时间的应用中,用于测量磁场图或点扩散函数的时间成本是微不足道的;但是,对于例如弥散加权成像(Diffusion Weighted Imaging,DWI)等具有短测量时间的应用,额外的时间不能忽略,具体而言,弥散加权成像基于回波平面成像(Echo Planar Imaging,EPI)方法(一种快速磁共振成像方法)。
发明内容
有鉴于此,本发明提供一种磁共振成像系统的磁场不均匀性值获取方法,包括如下步骤:
为了将磁共振成像系统的原始磁场均匀化为一目标磁场,针对所述磁共振成像系统通过动态匀场方法提供一磁场补偿量,其中所述动态匀场方法包括进行一三维低分辨率双回波梯度回波序列;
利用一第三通式获取所述磁场不均匀性值,所述第三通式是:
ΔB=ΔBoriginal+ΔBcompensating,
其中,ΔB是所述磁场不均匀性值,ΔBoriginal是所述原始磁场与所述目标磁场之间的差值,ΔBcompensating是所述磁场补偿量。
优选地,利用一第四通式获取所述差值,所述第四通式是:
ΔBoriginal=Δφ/(γ·ΔTE),
其中,ΔTE是所述三维低分辨率双回波梯度回波序列的双回波的回波时间的差值,Δφ是所述三维低分辨率双回波梯度回波序列生成的两个梯度回波图像的相位差,γ是旋磁比。
本发明还提供一种磁共振成像系统的失真校正方法,包括如下步骤:
根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移;
利用各个所述像素偏移对各个所述像素进行失真校正。
优选地,利用一第一通式根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移的步骤,所述第一通式是:
ΔnPE=α·ΔB+β,
其中,ΔnPE是各个所述像素偏移,ΔB是所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值,α是一转换参数,β是一调整参数。
优选地,所述磁共振图像是由回波平面成像方法生成的,其中,根据一第二通式得出所述转换参数,所述第二通式是:
α=γ·Tesp·NPE,
其中,α是所述转换参数,γ是旋磁比,Tesp是所述回波平面成像方法所得到的磁共振回波间距,NPE是所述磁共振图像在相位编码方向上的步进数。
本发明还提供一种磁共振成像方法,包括如上任一所述的磁场不均匀性值获取方法和/或如上任一所述的失真校正方法。
本发明还提供一种磁共振成像系统的磁场不均匀性值获取装置,包括:
一匀场单元,用于为了将磁共振成像系统的原始磁场均匀化为一目标磁场,针对所述磁共振成像系统通过动态匀场方法提供一磁场补偿量,其中所述动态匀场方法包括进行一三维低分辨率双回波梯度回波序列;
一获取单元,用于利用一第三通式获取所述磁场不均匀性值,所述第三通式是:
ΔB=ΔBoriginal+ΔBcompensating,
其中,ΔB是所述磁场不均匀性值,ΔBoriginal是所述原始磁场与所述目标磁场之间的差值,ΔBcompensating是所述磁场补偿量。
优选地,所述获取单元,进一步用于利用一第四通式获取所述差值,所述第四通式是:
ΔBoriginal=Δφ/(γ·ΔTE),
其中,ΔTE是所述三维低分辨率双回波梯度回波序列的双回波的回波时间的差值,Δφ是所述三维低分辨率双回波梯度回波序列生成的两个梯度回波图像的相位差,γ是旋磁比。
本发明还提供一种磁共振成像系统的失真校正装置,包括如下步骤:
一像素偏移计算单元,用于根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移;
一校正单元,用于利用各个所述像素偏移对各个所述像素进行失真校正。
优选地,所述像素偏移计算单元,用于利用一第一通式根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移的步骤,所述第一通式是:
ΔnPE=α·ΔB+β,
其中,ΔnPE是各个所述像素偏移,ΔB是所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值,α是一转换参数,β是一调整参数。
优选地,所述磁共振图像是由回波平面成像方法生成的,其中,根据一第二通式得出所述转换参数,所述第二通式是:
α=γ·Tesp·NPE,
其中,α是所述转换参数,γ是旋磁比,Tesp是所述回波平面成像方法所得到的磁共振回波间距,NPE是所述磁共振图像在相位编码方向上的步进数。
本发明还提供一种磁共振成像系统,包括如上任一所述的磁场不均匀性值获取装置和/或如上任一所述的失真校正装置。
从上述方案中可以看出,根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值方法可以直接利用动态匀场方法中产生的数据计算磁场不均匀性,因此节省了大量用于再次映射磁场的时间,进而缩短了磁共振成像的时间,提高了磁共振成像的效率,所以这项技术对于医院大量繁重的扫描工作而言非常宝贵。
根据本发明的具体实施例的磁共振成像方法的失真校正方法可以直接利用磁场不均匀性值计算得出像素偏移而无需进行测量,所以对于某些磁共振成像方法(例如回波平面成像方法和迪克逊水脂成像方法)而言,能够节省大量时间成本。
附图说明
下面将通过参照附图详细描述本发明的优选实施例,使本领域的普通技术人员更清楚本发明的上述及其它特征和优点,附图中:
图1是根据本发明的具体实施例的磁共振成像系统的失真校正方法的步骤图。
图2是根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值获取方法的步骤图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下举具体实施例对本发明进一步详细说明。
根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值方法的核心在于,使用磁场不均匀性值计算得到磁共振图像的各个像素在相位编码方向上的像素偏移,利用像素偏移对各个像素进行失真校正(即,图像配准)。
图1是根据本发明的具体实施例的磁共振成像方法的失真校正方法的步骤图。如图1所示,根据本发明的具体实施例的磁共振成像方法的失真校正方法100包括:步骤101,根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移;和步骤102,利用各个所述像素的所述像素偏移对各个所述像素进行失真校正。
具体而言,利用一第一通式根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移的步骤,所述第一通式是:
ΔnPE=α·ΔB+β,
其中,ΔnPE是各个所述像素在相位编码方向上的像素偏移,ΔB是所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值,α是一转换参数,β是一调整参数。可见,所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值和各个所述像素在相位编码方向上的像素偏移之间存在线性映射关系。
所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值和各个所述像素在相位编码方向上的像素偏移之间,根据不同磁共振成像方法分别存在一定的映射关系,线性或非线性。
以基于回波平面成像(Echo Planar Imaging,EPI)方法的磁共振成像方法为例,所述磁共振图像是由回波平面成像方法生成的,其中,根据一第二通式得出所述转换参数,所述第二通式是:
α=γ·Tesp·NPE,
其中,α是所述转换参数,γ是旋磁比,Tesp是所述回波平面成像方法所得到的磁共振回波间距,NPE是所述磁共振图像在相位编码方向上的步进数,所述调整参数是一用户设定值或经过线性拟合得出的计算值。具体而言,在本发明的具体实施例的磁共振成像方法的失真校正方法中,所述调整参数β是0。
使用如上步骤得出的各个像素在相位编码方向上的像素偏移,在后处理中对磁共振图像中的各个像素进行失真校正,即图像配准,也就是将各个像素在相位编码方向上的原始坐标与在相位编码方向上的像素偏移取相加,从而得出各个像素在相位编码方向上的校正坐标,进而根据校正坐标得出各个像素构成的磁共振图像。
根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值获取方法的核心在于:首先进行动态匀场,然后利用动态匀场的过程中产生的数据获取磁场不均匀性值。
图2是根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值获取方法的步骤图。如图2所示,根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值获取方法200,包括如下步骤:步骤101,为了将磁共振成像系统的原始磁场均匀化为一目标磁场,针对所述磁共振成像系统通过动态匀场方法提供一磁场补偿量,其中所述动态匀场方法包括进行一三维低分辨率双回波梯度回波序列;步骤102,利用一第三通式获取所述磁场不均匀性值,所述第三通式是:
ΔB=ΔBoriginal+ΔBcompensating
其中,ΔB是所述磁场不均匀性值,ΔBoriginal是所述原始磁场与所述目标磁场之间的差值,ΔBcompensating是所述磁场补偿量。
具体而言,在磁共振成像系统中包括采用动态匀场方法来改进磁场的均匀性,也就是将磁共振成像系统的原始磁场均匀化为一目标磁场。在动态匀场方法中,将原始磁场叠加一磁场补偿量,从而构成匀场后的一实际磁场。而磁场不均匀性值就是目标磁场和实际磁场之间的差值,同时磁场不均匀性值也就是原始磁场与目标磁场之间的差与磁场补偿量的和。
具体而言,利用一第四通式获取所述差值,所述第四通式是:
ΔBoriginal=Δφ/(γ·ΔTE)
其中,ΔTE是所述三维低分辨率双回波梯度回波序列的双回波的回波时间的差值,Δφ是所述三维低分辨率双回波梯度回波序列生成的两个梯度回波图像的相位差,γ是旋磁比。同时,ΔBoriginal是所述原始磁场与所述目标磁场之间的差值还可以通过多种其它方式获得,比如测量等等。
具体而言,利用一第五通式获取所述差值,所述第五通式是:
第五通式描述的是在球坐标系中的情况,其中,(r,θ,φ)是球坐标系的坐标;由于存在多项匀场线圈,因此利用(m,n)来区分各项匀场线圈;R0表示匀场区域的半径;和表示第(m,n)项匀场线圈的灵敏度;表示通入第(m,n)项线圈的电流的大小;是勒让德多项式。其中,R0,和都可以从磁共振成像系统中获取。通常来说R0,和跟系统硬件相关,不会变化;是动态匀场技术根据ΔBoriginal计算得到的,会随着扫描对象的变化而变化。
在动态匀场方法中,磁场测量序列,即三维低分辨率双回波梯度回波序列,获取将要进行匀场的区域的三维体数据,用以评估在每个匀场线圈中所需的匀场电流从而优化磁场均匀性。
图3是根据本发明的具体实施例的磁共振成像方法的步骤图。如图3所示并且如上所述,根据本发明的具体实施例的磁共振成像方法300包括:步骤301,进行根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值获取方法;步骤302,进行一回波平面成像方法(或一迪克逊水脂成像方法);步骤303,进行根据本发明的具体实施例的磁共振成像方法的失真校正方法。
根据本发明的具体实施例的磁共振成像系统的磁场不均匀性值方法可以直接利用动态匀场方法中产生的数据计算磁场不均匀性,因此节省了大量用于再次映射磁场的时间,进而缩短了磁共振成像的时间,提高了磁共振成像的效率,所以这项技术对于医院大量繁重的扫描工作而言非常宝贵。
根据本发明的具体实施例的磁共振成像方法的失真校正方法可以直接利用磁场不均匀性值计算得出像素偏移而无需进行测量,所以对于某些磁共振成像方法(例如回波平面成像方法和迪克逊水脂成像方法)而言,能够节省大量时间成本。
具体而言,以回波平面成像方法和迪克逊水脂成像方法为例:其采用基于梯度回波的双回波序列来映射原始磁场;随后基于原始磁场进行匀场,也就是为了使磁共振系统的主磁场成为一匀质磁场(即目标磁场),通过将原始磁场叠加补偿磁场从而生成一实际磁场(但是实际磁场与目标磁场之间仍存在差距),例如,匀场操作包括多种方法,例如,静态匀场方法和动态匀场方法,其中动态匀场方法即利用匀场电流生成补偿磁场。由此,原始磁场和补偿磁场的向量和构成实际磁场,进而目标磁场和实际磁场之间的差构成磁场不均匀性值。在根据本发明的具体实施例的磁共振成像方法的失真校正方法中,利用磁场不均匀性值计算像素偏移,因此节省了大量用于再次映射磁场的时间,进而缩短了磁共振成像的时间,提高了磁共振成像的效率,所以这项技术对于医院大量繁重的扫描工作而言非常宝贵。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种磁共振成像系统的磁场不均匀性值获取方法,包括如下步骤:
为了将磁共振成像系统的原始磁场均匀化为一目标磁场,针对所述磁共振成像系统通过动态匀场方法提供一磁场补偿量,其中所述动态匀场方法包括进行一三维低分辨率双回波梯度回波序列;
利用一第三通式获取所述磁场不均匀性值,所述第三通式是:
ΔB=ΔBoriginal+ΔBcompensating,
其中,ΔB是所述磁场不均匀性值,ΔBoriginal是所述原始磁场与所述目标磁场之间的差值,ΔBcompensating是所述磁场补偿量,
利用一第四通式获取所述差值,所述第四通式是:
ΔBoriginal=Δφ/(γ·ΔTE),
其中,ΔTE是所述三维低分辨率双回波梯度回波序列的双回波的回波时间的差值,Δφ是所述三维低分辨率双回波梯度回波序列生成的两个梯度回波图像的相位差,γ是旋磁比。
2.一种磁共振成像系统的失真校正方法,包括如下步骤:
根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移;
利用各个所述像素偏移对各个所述像素进行失真校正,
利用一第一通式根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移的步骤,所述第一通式是:
ΔnPE=α·ΔB+β,
其中,ΔnPE是各个所述像素偏移,ΔB是所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值,α是一转换参数,β是一调整参数。
3.如权利要求2所述的失真校正方法,其特征在于,所述磁共振图像是由回波平面成像方法生成的,其中,根据一第二通式得出所述转换参数,所述第二通式是:
α=γ·Tesp·NPE,
其中,α是所述转换参数,γ是旋磁比,Tesp是所述回波平面成像方法所得到的磁共振回波间距,NPE是所述磁共振图像在相位编码方向上的步进数。
4.一种磁共振成像方法,包括如权利要求1所述的磁场不均匀性值获取方法和/或如权利要求2-3任一所述的失真校正方法。
5.一种磁共振成像系统的磁场不均匀性值获取装置,包括:
一匀场单元,用于为了将磁共振成像系统的原始磁场均匀化为一目标磁场,针对所述磁共振成像系统通过动态匀场方法提供一磁场补偿量,其中所述动态匀场方法包括进行一三维低分辨率双回波梯度回波序列;
一获取单元,用于利用一第三通式获取所述磁场不均匀性值,所述第三通式是:
ΔB=ΔBoriginal+ΔBcompensating,
其中,ΔB是所述磁场不均匀性值,ΔBoriginal是所述原始磁场与所述目标磁场之间的差值,ΔBcompensating是所述磁场补偿量,
所述获取单元,进一步用于利用一第四通式获取所述差值,所述第四通式是:
ΔBoriginal=Δφ/(γ·ΔTE),
其中,ΔTE是所述三维低分辨率双回波梯度回波序列的双回波的回波时间的差值,Δφ是所述三维低分辨率双回波梯度回波序列生成的两个梯度回波图像的相位差,γ是旋磁比。
6.一种磁共振成像系统的失真校正装置,包括如下步骤:
一像素偏移计算单元,用于根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移;
一校正单元,用于利用各个所述像素偏移对各个所述像素进行失真校正,
所述像素偏移计算单元,用于利用一第一通式根据所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值得出各个所述像素在相位编码方向上的像素偏移的步骤,所述第一通式是:
ΔnPE=α·ΔB+β,
其中,ΔnPE是各个所述像素偏移,ΔB是所述磁共振成像系统的主磁场在所述磁共振成像系统得出的磁共振图像的各个像素上的磁场不均匀性值,α是一转换参数,β是一调整参数。
7.如权利要求6所述的失真校正装置,其特征在于,所述磁共振图像是由回波平面成像方法生成的,其中,根据一第二通式得出所述转换参数,所述第二通式是:
α=γ·Tesp·NPE,
其中,α是所述转换参数,γ是旋磁比,Tesp是所述回波平面成像方法所得到的磁共振回波间距,NPE是所述磁共振图像在相位编码方向上的步进数。
8.一种磁共振成像系统,包括如权利要求5所述的磁场不均匀性值获取装置和/或如权利要求6-7任一所述的失真校正装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410663071.XA CN105676155B (zh) | 2014-11-19 | 2014-11-19 | 磁共振成像系统的磁场不均匀性值获取方法和装置和失真校正方法和装置 |
US14/945,761 US9824448B2 (en) | 2014-11-19 | 2015-11-19 | Method for acquiring a magnetic field inhomogeneity value and distortion correction method for magnetic resonance imaging system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410663071.XA CN105676155B (zh) | 2014-11-19 | 2014-11-19 | 磁共振成像系统的磁场不均匀性值获取方法和装置和失真校正方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105676155A CN105676155A (zh) | 2016-06-15 |
CN105676155B true CN105676155B (zh) | 2019-02-22 |
Family
ID=55961465
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410663071.XA Active CN105676155B (zh) | 2014-11-19 | 2014-11-19 | 磁共振成像系统的磁场不均匀性值获取方法和装置和失真校正方法和装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US9824448B2 (zh) |
CN (1) | CN105676155B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102013216529B4 (de) * | 2013-08-21 | 2019-05-23 | Siemens Healthcare Gmbh | Verfahren insbesondere zur patientenadaptiven B0-Homogenisierung von MR-Systemen unter Verwendung unterschiedlicher Typen von Shim-Spulen |
CN108272452A (zh) * | 2017-01-05 | 2018-07-13 | 上海康达卡勒幅医疗科技有限公司 | 一种磁共振偏中心成像二阶匀场方法 |
CN108387857B (zh) * | 2017-12-25 | 2020-11-10 | 深圳先进技术研究院 | 一种用于磁共振成像的局部匀场系统及匀场方法 |
CN108761364B (zh) * | 2018-05-21 | 2020-06-05 | 上海健康医学院 | 一种核磁共振弛豫分析仪锁场系统及其应用方法 |
CN110074786B (zh) * | 2019-04-30 | 2022-12-06 | 上海东软医疗科技有限公司 | 核磁共振匀场方法、装置、计算设备及核磁共振成像系统 |
CN110610529B (zh) * | 2019-09-20 | 2023-01-03 | 上海联影医疗科技股份有限公司 | 磁共振扫描仪的主磁场均匀化方法、设备和介质 |
CN112986881B (zh) * | 2019-12-17 | 2022-12-20 | 上海联影医疗科技股份有限公司 | 磁共振系统校准方法、成像方法以及磁共振系统 |
CN113625212A (zh) * | 2021-08-18 | 2021-11-09 | 中国人民解放军东部战区总医院 | 一种磁共振系统b0场图的校正方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1910471A (zh) * | 2004-01-14 | 2007-02-07 | 皇家飞利浦电子股份有限公司 | 具有实时磁场绘制的磁共振成像 |
CN101470179A (zh) * | 2007-12-29 | 2009-07-01 | 西门子(中国)有限公司 | 磁共振成像中失真校准的方法和装置 |
CN101498774A (zh) * | 2008-01-31 | 2009-08-05 | 西门子公司 | 医学磁共振成像回波平面成像测量中的动态失真校正方法 |
CN101604008A (zh) * | 2009-07-09 | 2009-12-16 | 北京海思威科技有限公司 | 磁体一阶匀场方法 |
CN101966081A (zh) * | 2010-08-05 | 2011-02-09 | 华东师范大学 | 校正呼吸信息影响静息态功能磁共振脑成像的方法 |
CN102508182A (zh) * | 2011-11-30 | 2012-06-20 | 苏州安科医疗系统有限公司 | 一种用于磁共振成像的多阶谐波动态匀场方法 |
CN102597795A (zh) * | 2009-11-05 | 2012-07-18 | 皇家飞利浦电子股份有限公司 | 使用导航器的mr成像 |
CN103376433A (zh) * | 2012-04-27 | 2013-10-30 | 西门子(深圳)磁共振有限公司 | 图像畸变校正方法及系统、磁共振成像设备 |
CN103901373A (zh) * | 2012-12-28 | 2014-07-02 | 上海联影医疗科技有限公司 | 一种磁共振成像匀场方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5485085A (en) * | 1994-11-28 | 1996-01-16 | Board Of Trustees Of The Leland Stanford Junior University | Method for correcting magnetic field inhomogeneity distortion in MRI images having disconnected regions |
DE19954925C2 (de) * | 1999-11-16 | 2001-10-04 | Bruker Medical Gmbh | Verfahren zum Korrigieren von Feldinhomogenitäten höherer Ordnung in einer Apparatur der magnetischen Resonanz |
US8085041B2 (en) * | 2008-04-10 | 2011-12-27 | General Electric Company | Three-point method and system for fast and robust field mapping for EPI geometric distortion correction |
US9678184B2 (en) * | 2012-01-11 | 2017-06-13 | Hitachi, Ltd. | Method for increment of RF-phase based on static magnetic field inhomogeneity in MRI apparatus |
-
2014
- 2014-11-19 CN CN201410663071.XA patent/CN105676155B/zh active Active
-
2015
- 2015-11-19 US US14/945,761 patent/US9824448B2/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1910471A (zh) * | 2004-01-14 | 2007-02-07 | 皇家飞利浦电子股份有限公司 | 具有实时磁场绘制的磁共振成像 |
CN101470179A (zh) * | 2007-12-29 | 2009-07-01 | 西门子(中国)有限公司 | 磁共振成像中失真校准的方法和装置 |
CN101498774A (zh) * | 2008-01-31 | 2009-08-05 | 西门子公司 | 医学磁共振成像回波平面成像测量中的动态失真校正方法 |
CN101604008A (zh) * | 2009-07-09 | 2009-12-16 | 北京海思威科技有限公司 | 磁体一阶匀场方法 |
CN102597795A (zh) * | 2009-11-05 | 2012-07-18 | 皇家飞利浦电子股份有限公司 | 使用导航器的mr成像 |
CN101966081A (zh) * | 2010-08-05 | 2011-02-09 | 华东师范大学 | 校正呼吸信息影响静息态功能磁共振脑成像的方法 |
CN102508182A (zh) * | 2011-11-30 | 2012-06-20 | 苏州安科医疗系统有限公司 | 一种用于磁共振成像的多阶谐波动态匀场方法 |
CN103376433A (zh) * | 2012-04-27 | 2013-10-30 | 西门子(深圳)磁共振有限公司 | 图像畸变校正方法及系统、磁共振成像设备 |
CN103901373A (zh) * | 2012-12-28 | 2014-07-02 | 上海联影医疗科技有限公司 | 一种磁共振成像匀场方法 |
Also Published As
Publication number | Publication date |
---|---|
US9824448B2 (en) | 2017-11-21 |
US20160139220A1 (en) | 2016-05-19 |
CN105676155A (zh) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105676155B (zh) | 磁共振成像系统的磁场不均匀性值获取方法和装置和失真校正方法和装置 | |
JP6462306B2 (ja) | 医用画像処理装置及び磁気共鳴イメージング装置 | |
Jezzard | Correction of geometric distortion in fMRI data | |
JP6037652B2 (ja) | 拡散強調磁気共鳴データの生成方法、磁気共鳴システムおよびコンピュータ読み取り可能な記憶媒体 | |
US10488479B2 (en) | MRI apparatus and image processing apparatus | |
Chung et al. | Rapid B1+ mapping using a preconditioning RF pulse with TurboFLASH readout | |
US9513359B2 (en) | Systems and methods for shim current calculation | |
US20140210471A1 (en) | Method and control device for operating a magnetic resonance system | |
US10031201B2 (en) | Method and apparatus for magnetic resonance imaging | |
RU2015116879A (ru) | Опорное сканирование при устойчивой к металлам мр визуализации | |
JP2009050615A (ja) | 磁気共鳴イメージング装置および磁気共鳴画像表示方法 | |
Xu et al. | Prospective and retrospective high order eddy current mitigation for diffusion weighted echo planar imaging | |
JPWO2012077543A1 (ja) | 磁気共鳴イメージング装置及びコントラスト強調画像取得方法 | |
Haskell et al. | Off‐resonance artifact correction for MRI: A review | |
Bastin et al. | On the application of a non‐CPMG single‐shot fast spin‐echo sequence to diffusion tensor MRI of the human brain | |
US9494662B2 (en) | Magnetic resonance method and apparatus for automatic calculation of a maximum pulse-length of an excitation pulse | |
JP2007159718A (ja) | 磁気共鳴イメージング装置および画像補正評価方法 | |
US9945924B2 (en) | Method and magnetic resonance system for slice-selective magnetic resonance imaging | |
KR101617937B1 (ko) | Mr 데이터를 획득하고 b1 자기장을 판단하는 방법 및 이에 대응하여 설계된 자기 공명 시스템 | |
CN111615637B (zh) | 用于执行扩散加权磁共振测量的方法 | |
van de Looij et al. | Diffusion tensor echo planar imaging using surface coil transceiver with a semiadiabatic RF pulse sequence at 14.1 T | |
KR20130135776A (ko) | 자기 공명 기법에서 측정 볼륨 내의 검사 대상의 대상-특정 b1 분포를 결정하는 방법, 자기 공명 시스템, 컴퓨터 프로그램, 및 전자적으로 판독가능한 데이터 매체 | |
US9772390B2 (en) | Magnetic resonance imaging device and method for generating image using same | |
US20140086468A1 (en) | Method and apparatus for correction of artifacts in magnetic resonance images | |
US9417304B2 (en) | Method and magnetic resonance system to acquire MR data with diffusion information |
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 |