CN110988765B - 一种磁共振相位校正方法和装置 - Google Patents
一种磁共振相位校正方法和装置 Download PDFInfo
- Publication number
- CN110988765B CN110988765B CN201911357124.4A CN201911357124A CN110988765B CN 110988765 B CN110988765 B CN 110988765B CN 201911357124 A CN201911357124 A CN 201911357124A CN 110988765 B CN110988765 B CN 110988765B
- Authority
- CN
- China
- Prior art keywords
- phase error
- data
- image data
- order linear
- echo image
- 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/56545—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- Medical Informatics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请提供了一种磁共振相位校正方法和装置,先获得原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差,再根据所述两个回波图像数据确定读出方向的一阶线性相位误差;并在初始相位误差中去除一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差,最后根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。如此可利用采集数据自身信息去得到初始相位误差,且通过计算并去除一阶线性相位误差的方式避免相位缠绕,能较好的对磁共振进行相位校正。
Description
技术领域
本说明书涉及图像技术领域,尤其涉及一种磁共振相位校正方法和装置。
背景技术
核磁共振是指质子在外磁场作用下自旋方向分布满足玻尔兹曼分布,在外加特定频率的射频磁场作用下吸收能量,撤去射频磁场后发生弛豫释放能量的现象。磁共振成像(Magnetic Resonance Imaging,MRI)主要利用该原理,结合空间编码和傅里叶变换等技术,利用检测到的核磁共振信号还原出成像物体内部结构信息。早期研究表明特定MR参数,如质子密度(PD)、纵向弛豫时间(T1)、横向弛豫时间(T2/T2*),对于人体病例组织的检测具有高敏感度和组织特异性,定量磁共振技术应运而生。随着生活条件的进步,肥胖也成为重大和普遍疾病的一种。如果能准确的测量人体关键组织内的脂肪成分,对于早期预防很多疾病非常有用。
腹部定量磁共振技术可以定量测绘组织的R2*以及脂肪分数,通常用来辅助诊断病人是否有脂肪肝以及铁沉积。此方法一般采用多回波梯度回波序列采集数据,由于呼吸的原因,通常需要病人配合屏气,因此对采集时间要求非常严格,采用双极性梯度回波序列采集能够缩短回波间隔,进而缩短扫描时长,而且获得的图像信噪比更高,相比于同极性采集,也可以在相同时间内采集更多的回波数据。
然而,双极性梯度采集也会带来一些问题,快速切换梯度磁场会造成梯度延时和涡流,进而引入空间变化(XYZ)的相位误差。
发明内容
针对上述技术问题,本说明书实施例提供一种磁共振相位校正方法和装置,技术方案如下:
根据本说明书实施例的第一方面,提供一种磁共振相位校正方法,所述方法包括:
获取原始数据,所述原始数据为根据水脂信号模型建立的,具有多个回波的双极性梯度回波图像数据;
获得所述原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差值,所述初始相位误差包括线性相位误差与非线性相位误差;
获取所述原始数据中两个相邻回波的回波图像数据,根据所述两个回波图像数据确定读出方向的一阶线性相位误差;
在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差;
根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。
根据本说明书实施例的第二方面,提供一种磁共振相位校正装置,所述装置包括:
原始数据获取模块,被配置为获取原始数据,所述原始数据为根据水脂信号模型建立的,具有多个回波的双极性梯度回波图像数据;
初始相位误差确定模块,被配置为获得所述原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差,所述初始相位误差包括线性相位误差与非线性相位误差;
线性相位误差确定模块,被配置获取所述原始数据中两个相邻回波的回波图像数据,根据所述两个回波图像数据确定读出方向的一阶线性相位误差;
剩余相位误差确定模块,被配置为在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差;
误差校正模块。被配置为根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。
根据本说明书实施例的第三方面,提供一种智能计算机设备,其包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其中,处理器执行所述程序时实现如第一方面所述的磁共振相位校正方法。
根据本说明书实施例的第四方面,提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如第一方面所述的磁共振相位校正方法。
本说明书实施例所提供的技术方案,提供了一种磁共振相位校正方法,先获得原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差,再根据所述两个回波图像数据确定读出方向的一阶线性相位误差;并在初始相位误差中去除一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差,最后根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。如此可利用采集数据自身信息去得到初始相位误差,且通过计算并去除一阶线性相位误差的方式避免相位缠绕,能较好的对磁共振进行相位校正。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本说明书实施例。
此外,本说明书实施例中的任一实施例并不需要达到上述的全部效果。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书实施例中记载的一些实施例,对于本领域普通技术人员来讲,还可以根据这些附图获得其他的附图。
图1是本说明书一示例性实施例示出的磁共振相位校正设备的一种示意图;
图2是本说明书一示例性实施例示出的磁共振相位校正方法的一种流程图;
图3是本说明书一示例性实施例示出的磁共振相位校正方法的另一种流程图;
图4是本说明书一示例性实施例示出的磁共振相位校正方法的另一种流程图;
图5是本说明书一示例性实施例示出的磁共振相位校正方法的另一种流程图;
图6a-6e是本说明书一示例性实施例示出的多种水脂分离结果的示意图;
图7是本说明书一示例性实施例示出的磁共振相位校正装置的一种示意图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下所描述的示例性实施方式并不代表与本说明书相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本说明书的一些方面相一致的装置和方法的例子。
在本说明书使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本说明书。在本说明书和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本说明书可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本说明书范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
核磁共振是指质子在外磁场作用下自旋方向分布满足玻尔兹曼分布,在外加特定频率的射频磁场作用下吸收能量,撤去射频磁场后发生弛豫释放能量的现象。磁共振成像(Magnetic Resonance Imaging,MRI)主要利用该原理,结合空间编码和傅里叶变换等技术,利用检测到的核磁共振信号还原出成像物体内部结构信息。早期研究表明特定MR参数,如质子密度(PD)、纵向弛豫时间(T1)、横向弛豫时间(T2/T2*),对于人体病例组织的检测具有高敏感度和组织特异性,定量磁共振技术应运而生。随着生活条件的进步,肥胖也成为重大和普遍疾病的一种。如果能准确的测量人体关键组织内的脂肪成分,对于早期预防很多疾病非常有用。
腹部定量磁共振技术可以定量测绘组织的R2*以及脂肪分数,通常用来辅助诊断病人是否有脂肪肝以及铁沉积。此方法一般采用多回波梯度回波序列采集数据,由于呼吸的原因,通常需要病人配合屏气,因此对采集时间要求非常严格,采用双极性梯度回波序列采集能够缩短回波间隔,进而缩短扫描时长,而且获得的图像信噪比更高,相比于同极性采集,也可以在相同时间内采集更多的回波数据。
然而,双极性梯度采集也会带来一些问题,快速切换梯度磁场会造成梯度延时和涡流,进而引入空间变化(XYZ)的相位误差。
为了更好的理解本申请,首先介绍MRI系统。如图1所示,简单示意了MRI系统的组成,主要包括检查床110、磁体120、梯度线圈131-133、射频线圈140、主计算机150、梯度放大器160、射频控制器170和控制台180。
磁体120是产生主磁场的装置。梯度线圈131-133、梯度放大器160等组成梯度系统,主要用于产生梯度磁场以能够进行磁共振信号的空间定位编码。其中,梯度线圈由三组独立的线圈构成,分别为X轴梯度线圈131,Y轴梯度线圈132和Z轴梯度线圈133。一般来讲,X方向的空间编码为频率编码,Y方向的空间编码为相位编码,Z方向的空间编码为选层。射频线圈140、射频控制器170等组成射频系统,主要用于向被检体发射射频信号,再从被检体接收射频信号,因此射频线圈140有射频发射线圈和射频接收线圈之分,射频发射线圈和射频接收线圈可以为同一个线圈,也可以分为不同的线圈。为了简单起见,图1中只画了一个线圈140。主计算机150负责磁共振成像序列的发送、采集的射频数据的运算、磁共振图像的重建和显示等。
图2是根据一示例性实施例示出的一种磁共振相位校正方法的流程图,该磁共振相位校正方法可以用于能够执行磁共振相位校正的设备上,如图2所示,该磁共振相位校正方法包括以下步骤S21-S25:
在步骤S21中,获取原始数据,所述原始数据为根据水脂信号模型建立的,具有多个回波的双极性梯度回波图像数据;
在一实施例中,在执行步骤S21时,可以采用但不限于以下实现方式:
(1-1)采用MRI梯度回波序列进行信号采集,得到多回波模值图像数据和相位图像数据;
(1-2)将多回波模值图像数据和相位图像数据变换为多回波复值图像数据,将所述数据确定为原始数据。
在一实施例中,步骤S21中使用的水脂信号模型如下所示:
其中,S(n)为回波时间t对应的磁共振复数信号,W为水幅度值,F为脂肪幅度值,tn为第n个波的回波时间,相对模值am(m=1,2,...,M)衡量每个脂肪峰对总脂肪信号的相对贡献,M为脂肪峰数目;对于全体脂肪峰,有对于不同部位脂肪组织,M选取不同值,Δfm为每个脂肪峰所相对于水的化学位移。分别表示奇数波和偶数波的相位误差。
在步骤S22中,获得所述原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差值,所述初始相位误差包括线性相位误差与非线性相位误差;
在一实施例中,在执行步骤S22时,可以采用但不限于以下实现方式:
(2-1)获得所述原始数据中相邻三个回波的回波图像数据,对所述三组回波图像数据分别进行多倍欠采,将三组回波图像数据分别变换为对应的低分辨率回波图像数据;
(2-2)基于所述三组低分辨率的回波数据确定低分辨率的相位误差,将所述低分辨率的相位误差插值回原始分辨率,将所述原始分辨率的相位误差确定为初始相位误差。
在一实施例中,在执行(2-1)时,该相邻三个回波图像数据为原始数据的第一个回波图像数据、第二个回波图像数据和第三个回波图像数据。
在一实施例中,在执行(2-1)时,对所述三组回波图像数据分别进行多倍欠采为分别进行四倍欠采。
另外,上述(2-1)至(2-2)的具体实现过程,可见图3所示实施例。
在步骤S23中,获取所述原始数据中两个相邻回波的回波图像数据,根据所述两个回波图像数据确定读出方向的一阶线性相位误差;
在一实施例中,在执行步骤S23时,可以采用但不限于以下实现方式:
(3-1)获得所述两个回波数据的两条中心K线,将所述两条中心k线分别进行一维傅立叶变换,得到所述两条中心K线对应的图像域数据;
(3-1)将得到的两个图像域数据相除并确定计算结果的相位值,得到一阶线性相位误差。
在步骤S24中,在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差;
在一实施例中,在执行步骤S24时,可以采用但不限于以下实现方式:
(4-1)确定所述一阶线性相位误差的共轭值;
(4-2)将所述初始相位误差与所述共轭值相乘,以在所述初始相位误差中去除所述一阶线性相位误差;
(4-3)获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差。
具体而言,初始误差为包括线性误差和非线性误差的总体误差,其中包含大量的相位缠绕,而由于初始误差不是一阶数值,需要将初始误差取弧度值后用于校正,而存在相位缠绕时,计算初始误差的弧度值会由于相位缠绕的存而使计算不准确。
为了解决这个问题,本实施例先计算一阶线性相位误差,再在初始误差中先去除该一阶线性相位误差,以去除相位缠绕的主要产生原因,再确定剩余相位误差的弧度值,如此该剩余相位误差和一阶相位误差都可直接用于校正,基于不会有相位缠绕。
另外,上述(4-1)至(4-3)的具体实现过程,可见图4所示实施例。
在步骤S25中,根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。
在一实施例中,在执行步骤S25时,可以采用但不限于以下实现方式:
(5-1)对所述剩余相位误差进行平滑处理,得到平滑后的剩余相位误差;
(5-2)将所述平滑后的剩余相位误差与所述一阶线性相位误差相加,得到总相位误差;
(5-3)利用所述总相位误差校正所述原始数据。
另外,上述(5-1)至(5-3)的具体实现过程,可见图5所示实施例。
由上述实例可见,本实例先获得原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差,再根据所述两个回波图像数据确定读出方向的一阶线性相位误差;并在初始相位误差中去除一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差,最后根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。如此可利用采集数据自身信息去得到初始相位误差,且通过计算并去除一阶线性相位误差的方式避免相位缠绕,能较好的对磁共振进行相位校正。
图3是根据一示例性实施例示出的另一种磁共振相位校正方法的流程图,该磁共振相位校正方法可以用于能够执行磁共振相位校正的设备上,并建立在图2所示方法的基础上,如图3所示,在执行步骤S22时,可以包括以下步骤S31-S32:
在步骤S31中,获得所述原始数据中相邻三个回波的回波图像数据,对所述三组回波图像数据分别进行多倍欠采,将三组回波图像数据分别变换为对应的低分辨率回波图像数据;
在步骤S32中,基于所述三组低分辨率的回波数据确定低分辨率的相位误差,将所述低分辨率的相位误差插值回原始分辨率,将所述原始分辨率的相位误差确定为初始相位误差。
具体而言,根据所述三个回波图像数据得到初始相位误差的求解公式如下所示:
图4是根据一示例性实施例示出的另一种磁共振相位校正方法的流程图,该磁共振相位校正方法可以用于能够执行磁共振相位校正的设备上,并建立在图4所示方法的基础上,如图4所示,在执行步骤S24时,可以包括以下步骤S41-S43:
在步骤S41中,确定所述一阶线性相位误差的共轭值;
在步骤S42中,将所述初始相位误差与所述共轭值相乘,以在所述初始相位误差中去除所述一阶线性相位误差;
在步骤S43中,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差。
很多时候光学成像得到的信号是复数形式,包含了幅度值和相位值。但从复数信号中提取真实相位时,相位值会被限制在[-π,π]的区间内,位于该区间外的真实相位被缠绕到这一区间内。此现象称为相位缠绕,得到的相位称为缠绕相位。从缠绕相位中得到真实的相位称为相位解缠绕。
在水脂分离的误差计算过程中,一阶线性相位误差是多为缠绕相位,尤其是在一阶系数比较大的情况下,而初始相位误差是包含一阶线性相位误差在内的误差值,会有较严重的缠绕情况,对后续的校正分析造成困难,相位解缠绕算法在这种情况下也容易不稳定。
因此,本步骤在所述初始相位误差中去除所述一阶线性相位误差,如此可避免初始相位误差中大部分相位缠绕。
由于初始相位误差与一阶线性相位误差均为复数值,因此需要确定所述一阶线性相位误差的共轭值,将初始相位误差与得到的共轭值相乘,才能在初始相位误差中去除一阶线性相位误差。
在本申请一实施例中,可在得到剩余相位误差以常规解缠绕算法再进行一次相位解缠绕处理,去除可能存在的少许缠绕,此时应用相位解缠绕算法会更加稳定,错误较少。
图5是根据一示例性实施例示出的另一种磁共振相位校正方法的流程图,该磁共振相位校正方法可以用于能够执行磁共振相位校正的设备上,并建立在图5所示方法的基础上,如图5所示,在执行步骤S25时,可以包括以下步骤S51-S53:
在步骤S51中,对所述剩余相位误差进行平滑处理,得到平滑后的剩余相位误差;
在步骤S52中,将所述平滑后的剩余相位误差与所述一阶线性相位误差相加,得到总相位误差;
在步骤S53中,利用所述总相位误差校正所述原始数据。
本申请一实施例中,根据所述剩余相位误差与所述一阶线性相位误差得到总相位误差,基于所述总相位误差校正所述原始数据。该过程可称为粗校正过程,在这一过程执行完毕后,还可进一步执行精校正过程,精校正过程包括:根据校正后的原始数据得到场图,根据所述场图代入所述水脂信号模型的展开式,通过迭代计算得到进一步校正的数据。
其中,可将场图作为ψ的初始值,0作为θ的初始值,代入下列公式迭代计算,可以进一步修正相位以及模图误差:
S=E·D·A·P
之后,可根据泰勒公式,计算出Δψ,Δθ,ΔW,ΔF,直到满足迭代停止条件(例如Δψ<1)。
具体而言,将上文粗校正后的场图数据代入所述水脂信号模型的展开式进一步校正,相比将原始场图数据直接代入所述水脂信号模型进行校正来说,校正结果更加准确,可去除数据中的大部分相位误差和幅值误差。
在本申请一实施例中,获得进一步校正的数据后,还可利用T2*-ideal方法处理所述进一步校正的数据,得到最终的水图、脂图、R2*和/或脂肪分数。
具体校正效果可参考图6a-图6e,其中,图6a为只校正读出方向线性误差后的水、脂分离结果,可见只校正读出方向的相位误差是远远不够的,水、脂分离得不够彻底,脂图中残留大量水信号。
图6b通过中心k线和三个相邻回波进行上文所述粗校正后的水、脂分离结果;与图6a相比,结果有所改善,但还有一些水、脂分离不够彻底的地方,例如皮下脂肪,肝脏(箭头处)。
图6c为上文所述精校正后的水、脂分离结果。水脂分离效果最佳,不会影响临床诊断。
图6d和图6e为将校正后的数据带入T2*-ideal算法得到的定量信息,图6d为定量R2*/s-1图。图6e为脂肪定量图,最大值为1,代表全是脂肪,最小值为0,代表全是水。
可以看出,本申请利用采集数据自身信息去得到初始相位误差,避免相位缠绕,且分两步求解的相位校正方法,先粗略校正线性相位差,得到合适的初始值后再进行精准校正;同时校正线性,非线性相位误差以及幅值误差的精确校正方法针对腹部的应用,最后可提供多种诊断信息。通过对脂肪含量、R2*提供的铁沉积信息等综合分析可应用于肝硬化的分级和精准定量。
相应于上述方法实施例,本说明书实施例还提供一种磁共振相位校正装置,参见图7所示,所述装置包括:原始数据获取模块710,初始相位误差确定模块720,线性相位误差确定模块730,剩余相位误差确定模块740和误差校正模块750。
原始数据获取模块710,被配置为获取原始数据,所述原始数据为根据水脂信号模型建立的,具有多个回波的双极性梯度回波图像数据;
初始相位误差确定模块720,被配置为获得所述原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差,所述初始相位误差包括线性相位误差与非线性相位误差;
线性相位误差确定模块730,被配置获取所述原始数据中两个相邻回波的回波图像数据,根据所述两个回波图像数据确定读出方向的一阶线性相位误差;
剩余相位误差确定模块740,被配置为在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差;
误差校正模块750。被配置为根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。
本说明书实施例还提供一种磁共振相位校正设备,所述磁共振相位校正设备至少包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其中,处理器执行所述程序时实现如下磁共振相位校正方法:
获取原始数据,所述原始数据为根据水脂信号模型建立的,具有多个回波的双极性梯度回波图像数据;
获得所述原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差值,所述初始相位误差包括线性相位误差与非线性相位误差;
获取所述原始数据中两个相邻回波的回波图像数据,根据所述两个回波图像数据确定读出方向的一阶线性相位误差;
在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差;
根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。
本说明书实施例还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如下磁共振相位校正方法:
获取原始数据,所述原始数据为根据水脂信号模型建立的,具有多个回波的双极性梯度回波图像数据;
获得所述原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差值,所述初始相位误差包括线性相位误差与非线性相位误差;
获取所述原始数据中两个相邻回波的回波图像数据,根据所述两个回波图像数据确定读出方向的一阶线性相位误差;
在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差;
根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。
对于装置实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本说明书方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机,计算机的具体形式可以是个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件收发设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任意几种设备的组合。
以上所述仅是本说明书的具体实施例,应当指出,对于本技术领域的普通技术人员来说,在不脱离本说明书实施例原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本说明书实施例的保护范围。
Claims (16)
1.一种磁共振相位校正方法,其特征在于,所述方法包括:
获取原始数据,所述原始数据为根据水脂信号模型建立的,具有多个回波的双极性梯度回波图像数据;
获得所述原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差,所述初始相位误差包括线性相位误差与非线性相位误差;
获取所述原始数据中两个相邻回波的回波图像数据,根据所述两个回波图像数据确定读出方向的一阶线性相位误差;
在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差;
根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。
2.如权利要求1所述的方法,其特征在于,所述获得所述原始数据中相邻三个回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差值,包括:
获得所述原始数据中相邻三个回波的回波图像数据,对所述三个回波图像数据分别进行多倍欠采,将三个回波图像数据分别由原始分辨率的回波图像数据变换为对应的低分辨率回波图像数据;
基于所述三个低分辨率的回波数据确定低分辨率的相位误差,将所述低分辨率的相位误差插值回原始分辨率,将所述原始分辨率的相位误差确定为初始相位误差。
3.如权利要求1所述的方法,其特征在于,所述在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差,包括:
确定所述一阶线性相位误差的共轭值;
将所述初始相位误差与所述共轭值相乘,以在所述初始相位误差中去除所述一阶线性相位误差;
获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差。
4.如权利要求1所述的方法,其特征在于,所述根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据,包括:
对所述剩余相位误差进行平滑处理,得到平滑后的剩余相位误差;
将所述平滑后的剩余相位误差与所述一阶线性相位误差相加,得到总相位误差;
利用所述总相位误差校正所述原始数据。
5.如权利要求1所述的方法,其特征在于,所述根据所述两个回波数据确定读出方向的一阶线性相位误差,包括:
获得所述两个回波数据的两条中心K线,将所述两条中心k线分别进行一维傅立叶变换,得到所述两条中心K线对应的图像域数据;
将得到的两个图像域数据相除并确定计算结果的相位值,得到一阶线性相位误差。
6.如权利要求1所述的方法,其特征在于,所述根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据后,还包括:
根据校正后的原始数据得到场图,根据所述场图代入所述水脂信号模型的展开式,通过迭代计算得到进一步校正的数据。
7.如权利要求6所述的方法,其特征在于,所述获得进一步校正的数据后,还包括:
利用T2*-ideal方法处理所述进一步校正的数据,得到最终的水图、脂图、R2*和/或脂肪分数。
8.一种磁共振相位校正装置,其特征在于,所述装置包括:
原始数据获取模块,被配置为获取原始数据,所述原始数据为根据水脂信号模型建立的,具有多个回波的双极性梯度回波图像数据;
初始相位误差确定模块,被配置为获得所述原始数据中三个相邻回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差,所述初始相位误差包括线性相位误差与非线性相位误差;
线性相位误差确定模块,被配置获取所述原始数据中两个相邻回波的回波图像数据,根据所述两个回波图像数据确定读出方向的一阶线性相位误差;
剩余相位误差确定模块,被配置为在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差;
误差校正模块,被配置为根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据。
9.如权利要求8所述的装置,其特征在于,所述获得所述原始数据中相邻三个回波的回波图像数据,根据所述三个回波图像数据得到初始相位误差值,包括:
获得所述原始数据中相邻三个回波的回波图像数据,对所述三个回波图像数据分别进行多倍欠采,将三个回波图像数据分别变换为对应的低分辨率回波图像数据;
基于所述三个低分辨率的回波数据确定低分辨率的相位误差,将所述低分辨率的相位误差插值回原始分辨率,将所述原始分辨率的相位误差确定为初始相位误差。
10.如权利要求8所述的装置,其特征在于,所述在所述初始相位误差中去除所述一阶线性相位误差,获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差,包括:
确定所述一阶线性相位误差的共轭值;
将所述初始相位误差与所述共轭值相乘,以在所述初始相位误差中去除所述一阶线性相位误差;
获取去除结果的弧度值,确定为无相位缠绕的剩余相位误差。
11.如权利要求8所述的装置,其特征在于,所述根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据,包括:
对所述剩余相位误差进行平滑处理,得到平滑后的剩余相位误差;
将所述平滑后的剩余相位误差与所述一阶线性相位误差相加,得到总相位误差;
利用所述总相位误差校正所述原始数据。
12.如权利要求8所述的装置,其特征在于,所述根据所述两个回波数据确定读出方向的一阶线性相位误差,包括:
获得所述两个回波数据的两条中心K线,将所述两条中心k线分别进行一维傅立叶变换,得到所述两条中心K线对应的图像域数据;
将得到的两个图像域数据相除并确定计算结果的相位值,得到一阶线性相位误差。
13.如权利要求8所述的装置,其特征在于,所述根据所述剩余相位误差与所述一阶线性相位误差校正所述原始数据后,还包括:
根据校正后的原始数据得到场图,根据所述场图代入所述水脂信号模型的展开式,通过迭代计算得到进一步校正的数据。
14.如权利要求13所述的装置,其特征在于,所述获得进一步校正的数据后,还包括:
利用T2*-ideal方法处理所述进一步校正的数据,得到最终的水图、脂图、R2*和/或脂肪分数。
15.一种磁共振相位校正设备,其特征在于,其包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其中,处理器执行所述程序时实现如权利要求1-7任一项所述的磁共振相位校正方法。
16.一种计算机可读存储介质,其特征在于,其上存储有计算机程序,该程序被处理器执行时实现如权利要求1-7任一项所述的磁共振相位校正方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911357124.4A CN110988765B (zh) | 2019-12-25 | 2019-12-25 | 一种磁共振相位校正方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911357124.4A CN110988765B (zh) | 2019-12-25 | 2019-12-25 | 一种磁共振相位校正方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110988765A CN110988765A (zh) | 2020-04-10 |
CN110988765B true CN110988765B (zh) | 2022-03-04 |
Family
ID=70076539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911357124.4A Active CN110988765B (zh) | 2019-12-25 | 2019-12-25 | 一种磁共振相位校正方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110988765B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113281690B (zh) * | 2021-05-18 | 2022-08-12 | 上海联影医疗科技股份有限公司 | 一种磁共振成像方法和系统 |
CN115656904B (zh) * | 2022-12-13 | 2023-04-07 | 中国科学技术大学先进技术研究院 | 磁共振梯度延时自动校正方法、装置、设备及存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4321619A (en) * | 1979-10-24 | 1982-03-23 | Ampex Corporation | Second order velocity error correction for time base correctors |
CA1137177A (fr) * | 1979-04-27 | 1982-12-07 | Michel Levy | Procede de reduction des bruits de phase a la reception d'une transmission de donnees |
EP1199578A2 (en) * | 2000-10-20 | 2002-04-24 | GE Medical Systems Global Technology Company LLC | Phase error measuring method and apparatus, phase error correcting method and apparatus, recording medium and magnetic resonance imaging apparatus |
CN103901428A (zh) * | 2014-03-28 | 2014-07-02 | 西安电子科技大学 | 弹载sar子孔径前斜视高阶非线性调频变标成像方法 |
CN104569882A (zh) * | 2014-12-29 | 2015-04-29 | 苏州朗润医疗系统有限公司 | 磁共振成像系统快速自旋回波的相位校正方法及其应用 |
CN106054188A (zh) * | 2016-06-24 | 2016-10-26 | 西安电子科技大学 | 无人机合成孔径雷达成像的图像偏移自聚焦方法 |
CN106772342A (zh) * | 2017-01-11 | 2017-05-31 | 西南石油大学 | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 |
CN110118950A (zh) * | 2019-06-19 | 2019-08-13 | 华东师范大学 | 一种腹部定量磁化率成像中双极读出梯度的相位校正方法 |
-
2019
- 2019-12-25 CN CN201911357124.4A patent/CN110988765B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1137177A (fr) * | 1979-04-27 | 1982-12-07 | Michel Levy | Procede de reduction des bruits de phase a la reception d'une transmission de donnees |
US4321619A (en) * | 1979-10-24 | 1982-03-23 | Ampex Corporation | Second order velocity error correction for time base correctors |
EP1199578A2 (en) * | 2000-10-20 | 2002-04-24 | GE Medical Systems Global Technology Company LLC | Phase error measuring method and apparatus, phase error correcting method and apparatus, recording medium and magnetic resonance imaging apparatus |
CN103901428A (zh) * | 2014-03-28 | 2014-07-02 | 西安电子科技大学 | 弹载sar子孔径前斜视高阶非线性调频变标成像方法 |
CN104569882A (zh) * | 2014-12-29 | 2015-04-29 | 苏州朗润医疗系统有限公司 | 磁共振成像系统快速自旋回波的相位校正方法及其应用 |
CN106054188A (zh) * | 2016-06-24 | 2016-10-26 | 西安电子科技大学 | 无人机合成孔径雷达成像的图像偏移自聚焦方法 |
CN106772342A (zh) * | 2017-01-11 | 2017-05-31 | 西南石油大学 | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 |
CN110118950A (zh) * | 2019-06-19 | 2019-08-13 | 华东师范大学 | 一种腹部定量磁化率成像中双极读出梯度的相位校正方法 |
Non-Patent Citations (2)
Title |
---|
Higher order ionospheric effects in precise GNSS positioning;M. Mainul Hoque 等;《JOURNAL OF GEODESY》;20070430;第81卷;第259-268页 * |
基于区域识别和区域扩展的相位解缠算法;毕海霞 等;《电波科学学报》;20150430;第30卷(第2期);第244-251页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110988765A (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5717839B2 (ja) | 正規化非線形逆再構成法を用いて磁気共鳴画像シーケンスを再構成するための方法及び装置 | |
US10634753B2 (en) | MR imaging with motion detection | |
US8228063B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
JP2013521955A5 (zh) | ||
RU2672151C2 (ru) | Коррекция посторонних эхосигналов epi | |
US11002815B2 (en) | System and method for reducing artifacts in echo planar magnetic resonance imaging | |
WO2003092497A1 (fr) | Dispositif d'imagerie par resonance magnetique | |
JPWO2010035569A1 (ja) | 磁気共鳴イメージング装置 | |
CN108459289A (zh) | 一种基于数据一致性的多次激发弥散加权磁共振成像方法 | |
CN110988765B (zh) | 一种磁共振相位校正方法和装置 | |
JP6986659B2 (ja) | 磁気共鳴イメージング装置、ナイキストゴースト補正方法、及び、ナイキストゴースト補正用プログラム | |
JP4679158B2 (ja) | 磁気共鳴イメージング装置 | |
JP2017529960A (ja) | アーチファクト抑制を有するプロペラmrイメージング | |
US9535148B2 (en) | Dynamic contrast enhanced magnetic resonance imaging with high spatial-temporal resolution | |
JP2000157507A (ja) | 核磁気共鳴撮影装置 | |
JP4975614B2 (ja) | 磁気共鳴イメージング装置及び方法 | |
JP2010233907A (ja) | 磁気共鳴イメージング装置及び感度補正方法 | |
EP3394628A1 (en) | Dti with correction of motion-induced diffusion gradient inconsistency | |
JP4675936B2 (ja) | 核磁気共鳴撮影装置 | |
JP3983792B2 (ja) | 核磁気共鳴撮影装置 | |
JP5156958B2 (ja) | 磁気共鳴イメージング装置 | |
US20190254596A1 (en) | Method for generating a magnetic resonance image dataset, computer program product, data medium, and magnetic resonance system | |
EP3931588A1 (en) | Epi mr imaging with distortion correction | |
JP3952310B2 (ja) | 核磁気共鳴撮影装置 | |
JP4906952B2 (ja) | 核磁気共鳴撮影装置 |
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 |