CN101470179A - 磁共振成像中失真校准的方法和装置 - Google Patents
磁共振成像中失真校准的方法和装置 Download PDFInfo
- Publication number
- CN101470179A CN101470179A CNA2007103016660A CN200710301666A CN101470179A CN 101470179 A CN101470179 A CN 101470179A CN A2007103016660 A CNA2007103016660 A CN A2007103016660A CN 200710301666 A CN200710301666 A CN 200710301666A CN 101470179 A CN101470179 A CN 101470179A
- Authority
- CN
- China
- Prior art keywords
- image
- distortion
- zone
- pixel
- calibration
- 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
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/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
-
- 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
-
- 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/5617—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 RF refocusing, e.g. RARE
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (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
本发明公开了一种磁共振成像中失真校准方法,包括:在复数个数据读出方向上采集K空间数据,并转换为复数个图像,为所述每个图像分别生成对应的像素偏移图,根据所述像素偏移图对对应的图像进行几何失真校准,再对所有经过几何失真校准的图像进行结合。本发明还公开了一种磁共振成像中失真校准的装置。由于运动失真通常存在于数据读出方向,因此,从复数个数据读出方向上采集K空间数据可以有效地对运动失真进行校准。而且,根据像素偏移图对从复数个数据读出方向上采集的数据转换后的图像进行几何失真校准,可以减小结合后生成的最终图像的几何失真。使用了本发明的方案,不仅可以对MRI图像进行运动失真的校准,也可以实现MRI图像的几何失真校准。
Description
技术领域
本发明涉及磁共振成像(MRI,Magnetic Resonance Imaging)技术,特别是涉及一种MRI中失真校准的方法和装置。
背景技术
在MRI过程中,由于各种因素的影响,会产生失真。从失真的产生来源划分,MRI中普遍存在的失真主要包括:磁场相关的失真和与人体相关的失真。
磁场相关的失真主要是由于磁场不均匀,尤其是主磁场的不均匀性导致的几何失真。当存在较大的磁场不均匀性,很多信号可能互相重叠,从而出现几何失真。很多因素可以导致磁场不均匀,例如病人体存在内置金属设备、金属固定物等等。
与人体相关的失真主要是运动失真,其产生的原因是,在长时间的MR扫描过程中,病人有可能会失去耐心并移动身体,对于年龄非常小的或躁动的病人,其运动伪影尤其明显。如果运动伪影比较严重的话,图像质量会很差,对诊断造成不良影响。
目前普遍采用螺旋桨成像(PROPELLER,Periodically Rotated Overlapping ParallELLines Enhanced Reconstruction)技术进行运动校准,即采用一种独特的K空间填充模式,辐射状的叶片(blade)以螺旋的方式采集数据,指导整个K空间数据填满。这一技术可以在最复杂困难的情况下消除运动伪影,获得高质量图像。K空间也称为频率空间,其空间的纵横坐标分别代表在相位编码方向和频率编码方向的频率变化。
虽然PROPELLER技术可以有效地对运动失真进行校准,但是到目前为止,PROPELLER技术还不能对磁场不均匀性或磁化率的几何失真进行校准。
发明内容
本发明的一个目的在于提出一种MRI中失真校准的方法,能够有效地对MRI中的运动失真和几何失真进行校准。
本发明的目的在于提供一种与上述MRI中失真校准的方法相应的装置。
为了实现上述目的,本发明提出一种MRI中失真校准的方法,在复数个数据读出方向上采集K空间数据,并转换为复数个图像;为所述每个图像分别生成对应的像素偏移图;根据所述像素偏移图对对应的图像进行几何失真校准;对所有经过几何失真校准的图像进行结合。
其中,所述根据像素偏移图对对应的图像进行几何失真校准包括:根据像素偏移图计算对应图像中每个像素的偏移位置;根据该偏移位置以及该偏移位置上像素的值计算几何失真校准后的图像的像素值。
所述方法还包括:在对所有图像进行几何失真校准之后,确定所有经过几何失真校准的图像中是否存在严重失真的区域;利用其他图像中与所述严重失真区域相对应的区域的像素值替换所述严重失真区域的像素值。
其中,所述确定所有经过几何失真校准的图像中是否存在严重失真的区域包括:对于所有经过几何失真校准的图像中的每个图像,根据像素偏移图计算图像保真图,所述图像保真图反映均匀体模的像素偏移并重叠的程度;当所述图像保真图中存在像素值超过阈值的区域时,确定为所述经过几何失真校准的图像中存在严重失真区域。
其中,所述利用其他图像中与所述严重失真区域对应的区域的像素值替换所述严重失真区域的像素值包括:搜索其他图像中与所述严重失真区域对应的区域中是否存在严重失真的像素;当与所述严重失真区域对应的区域中不存在严重失真像素的图像的数量为一时,利用该图像中与所述严重失真区域对应的区域的像素值替换所述严重失真区域中像素的值;当与所述严重失真区域对应的区域中不存在严重失真像素的图像的数量大于一时,将所述一个以上图像中与所述严重失真区域对应的区域相结合,将结合后的值替换所述严重失真区域中像素的值。
其中,所述对所有经过几何失真校准的图像进行结合包括:在图像域对所有经过几何失真校准的图像相结合。
其中,所述对所有经过几何失真校准的图像进行结合包括:将所有经过几何失真校准的图像转换到K空间中;将对应于所有经过几何失真校准的图像的K空间数据相结合;将结合后的K空间数据转换为最终图像。
为了实现上述目的,本发明提出一种MRI中失真校准的装置,包括:采集单元用于在复数个数据读出方向上采集K空间数据,并转换为复数个图像;几何失真校准单元,用于为每个图像生成对应的像素偏移图,并根据所述像素偏移图对所有图像分别进行几何失真校准;图像结合单元,用于对所有经过几何失真校准的图像进行结合。
其中,所述几何失真校准单元包括:像素偏移图生成模块,用于为每个图像生成对应的像素偏移图;几何失真校准模块,用于根据像素偏移图计算对应图像中每个像素的偏移位置,并根据该偏移位置以及该偏移位置上像素的值计算几何失真校准后的图像的像素值。
所述装置还包括:严重失真区域替换单元,位于所述几何失真校准单元和图像结合单元之间,用于确定所有经过几何失真校准的图像中是否存在严重失真的区域,并利用其他图像中与所述严重失真区域相对应的区域的像素值替换所述严重失真区域的像素值。
其中,所述严重失真区域替换单元对于所有经过几何失真校准的图像中的每个图像,根据对应的像素偏移图计算图像保真图,所述图像保真图反映均匀体模的像素偏移并重叠的程度;当所述图像保真图中存在像素值超过阈值的区域时,确定为所述经过几何失真校准的图像中存在严重失真区域。
其中,所述严重失真区域替换单元搜索其他图像中与所述严重失真区域对应的区域中是否存在严重失真的像素;当与所述严重失真区域对应的区域中不存在严重失真像素的图像的数量为一时,利用该图像中与所述严重失真区域对应的区域的像素值替换所述严重失真区域中像素的值;当与所述严重失真区域对应的区域中不存在严重失真像素的图像的数量大于一时,将所述一个以上图像中与所述严重失真区域对应的区域相结合,将结合后的值替换所述严重失真区域中像素的值。
其中,所述图像结合单元在图像域将所有经过几何失真校准的图像相结合。
其中,所述图像结合单元将所有经过几何失真校准的图像转换到K空间中;将对应于所有经过几何失真校准的图像的复数个K空间数据相结合;再将结合后的K空间数据转换为最终图像。
由上述技术方案可以看出,在本发明中,在复数个数据读出方向上采集K空间数据,并转换为复数个图像,生成像素偏移图,根据像素偏移图对复数个图像分别进行几何失真校准,再对复数个经过几何失真校准的图像进行结合。由于运动失真通常存在于数据读出方向,因此,从复数个数据读出方向上采集K空间数据可以有效地对运动失真进行校准。而且,根据像素偏移图对从复数个数据读出方向上采集的数据转换后的图像进行几何失真校准,可以减小结合后生成的最终图像的几何失真。综上所述,使用了本发明的方案,不仅可以对MRI图像进行运动失真的校准,也可以实现MRI图像的几何失真校准。
附图说明
下面将通过参照附图详细描述本发明的优选实施例,使本领域的普通技术人员更清楚本发明的上述及其它特征和优点,相同的标号表示相同的部件,附图中:
图1是本发明提出的MRI中失真校准的方法的流程图;
图2是根据本发明实施例一的MRI中失真校准的方法的流程图;
图3是PROPELLER技术采用的一个回波链中基于TSE的序列示意图;
图4是PROPELLER技术中一次采集K空间的示意图;
图5是PROPELLER技术中采集时叶片旋转的示意图;
图6是由于磁场不均匀导致的像素位置偏移和失真的示意图;
图7是根据本发明实施例二的MRI中失真校准的方法的流程图;
图8是根据本发明实施例三的MRI中失真校准的装置结构图;
图9是根据本发明实施例四的MRI中失真校准的装置结构图;
图10示出了不同的RO方向导致不同重叠模式的示意图;
图11是对本发明进行仿真的流程图;
图12是创建的体模的原始图像;
图13是用于失真模拟的磁场不均匀分布图;
图14是对基于TSE序列的模拟原始数据直接逆傅立叶变换后的所有叶片的失真图像;
图15是使用磁场图对所有叶片进行失真校准后的图像;
图16是根据保真图对所有叶片中严重失真的区域进行修复的图像;
图17是原始无失真图像、采用PROPELLER采集并经校准后的图像、采用PROPELLER采集但未经校准的图像的对比示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
图1是本发明提出的MRI中失真校准的方法的流程图。如图1所示,在步骤S101,在复数个数据读出方向上采集K空间数据,并转换为复数个图像;在步骤S102,生成分别对应于复数个图像的复数个像素偏移图;在步骤S103,根据复数个像素偏移图对复数个图像分别进行几何失真校准;在步骤S104,对复数个经过几何失真校准的图像进行结合。
下面通过几个具体实施例对本发明进行详细阐述。
图2是根据本发明实施例一的MRI中失真校准的方法的流程图。如图2所示,在MRI中,进行失真校准主要包括如下步骤:
在步骤S201,使用PROPELLER技术对用于成像的K空间数据进行采集,并转换为复数个图像。
对于快速自旋回波(TSE,Turbo Spin Echo)序列,失真仅出现在数据读出(RO,ReadOut)方向,因此,如果一种成像方法仅具有单一的RO方向,这些重叠的信号较难被恢复出来。而本发明使用PROPELLER采集技术,可以在复数个RO方向上采集数据。
图3是PROPELLER技术采用的一个回波链中基于TSE的序列示意图。如图3所示,在本实施例中,通过将传统的读出梯度和相位编码梯度相组合,实现采集不同叶片时,其实际编码方向的旋转。图4是PROPELLER技术中一次采集K空间的示意图。如图4所示,在一次采集中,在K空间中心采集L条相位编码线。这些线组成叶片,叶片的宽度是一次采集的相位数,也就是回波链的长度。叶片的宽度L的值可以任意设置。
在扫描并采集完一个叶片的数据之后,进行叶片旋转。图5是本实施例的PROPELLER采集中叶片旋转的示意图。如图5所示,初始的RO方向与横坐标轴重合,在对一个叶片进行采集之后,通过调整传统读出梯度和相位编码梯度的幅度,将有效的RO方向旋转角度θ,至一个新的方向a(n+1),采集另一个叶片。这里,θ可以任意设置。通常,磁场不均匀性越严重,θ的值越小。按照相同的方式,连续旋转RO的方向,直至所有的叶片基本上覆盖了整个K空间,即K空间被填满。整个K空间的大小由理想的视野(FOV)和空间分辨率确定。
在采集完成之后,将所有叶片的数据转换为复数个图像,以进行后处理。
在步骤S202,生成复数个像素偏移图IΔx。
像素偏移图IΔx表示了由磁场不均匀性导致的信号偏移的距离。像素偏移图IΔx与磁场不均匀性图ΔB(x,y)成正比。例如,采用梯度回波(GRE,gradient echo)序列,计算两幅具有不同回波时间(TE,echo time)的图像之间的相位差,该相位差图可反映磁场图ΔB(x,y),根据公式(1)计算像素偏移图IΔx。
其中,Δx是在RO方向上的像素偏移距离;Nro表示RO方向上的图像矩阵大小;Gro表示RO梯度的幅度;Lro表示RO方向上的FOV。当然,也可以采用其他计算像素偏移图测量的方法。
在步骤S203,使用像素偏移图校准几何失真。
图6是由于磁场不均匀导致的像素位置偏移和失真的示意图。如图6所示,在RO方向上,如果主磁场是均匀的,位置Xn+1和Xn+2之间的信号形成傅立叶变换后的重建图像中的像素Pn+1。由于磁场不均匀,位置Xn+1偏移到位置Xn+1’,位置Xn+2偏移到位置Xn+2’,因此实际上位于Pn+1的部分信号偏移到像素Pm+1,其他Pn+1的信号偏移到傅立叶变换后的重建图像中的Pm+2。这里,位置Xn+1’和位置Xn+2’可以从像素偏移图中计算出来。
在此介绍一种可能的使用像素偏移图校准失真的方式。为了恢复像素Pn+1,首先,从像素偏移图IΔx计算Xn+1’,Xn+2’,然后根据公式(2)计算:
其中,Pm+1和Pm+2是直接逆傅立叶变换后的图像中的像素。
以相同的方式,对所有像素进行校准。在下一次处理之前,对所有叶片进行校准。
在步骤S204,将所有叶片的数据相结合,得到最终图像。
具体的结合方法可以是在图像空间中直接结合,得到最终图像。也可以在K空间结合,然后转换到图像空间得到最终图像。当采用在K空间结合时,将所有叶片的数据转换到K空间,将这些数据进行栅格重置(regrid),至笛卡儿坐标系中,形成一个完整的K空间,然后通过逆傅立叶变换获得最终图像。
图7是根据本发明实施例二的MRI中失真校准的方法的流程图。本实施例与实施例一的差别在于,本实施例在数据后处理阶段,增加恢复严重失真区域的处理,从而进一步加强几何失真校准的性能,得到更佳的图像质量。如图7所示,在本实施例中,进行失真校准主要包括如下步骤:
在步骤S701,使用PROPELLER类技术对用于成像的K空间数据进行采样,并转换为复数个图像。
在步骤S702,生成复数个像素偏移图IΔx。
在步骤S703,使用像素偏移图校准复数个图像的几何失真。
以上步骤S701至S703采用与实施例一的步骤S201至S203相同或类似的处理方式。
在步骤S704,确定每个叶片数据中的严重失真区域。
根据像素偏移图IΔx,计算图像保真图(image fidelity map)If。具体的方法可以是:
首先,假设设置一个均匀的体模,即所有像素具有相同的值的体模,该体模处在同一磁场不均匀度BΔ之下,然后计算所有叶片的RO方向上的失真图像。在这些失真图像中,像素值越大,重叠的信号越多。
第二步,计算这些重叠图像的来源,在这些位置应用失真的像素的值,以形成所有叶片的保真图。例如,假设均匀体模中的三个像素Pl=1,Pl+1=1,Pl+2=1一起偏移并重叠到失真图像的像素Pk=3,然后保真图If中的这三个位置具有失真图像中Pk的值。这样,可以计算出所有叶片的一系列保真图If1,If2,...。因此,图像保真图反映了均匀体模的像素偏移并重叠的程度。
第三步,设置一个阈值,对于各保真图If中大于该阈值的值,标记为各叶片的严重失真区域。
在步骤S705,恢复每个叶片的严重失真区域。
对于一个叶片,如果存在严重失真区域的话,使用其他叶片的相同区域内的对应像素替代该严重失真区域。
具体方法可以是:搜索其他叶片的图像中与严重失真区域对应的区域中是否存在严重失真的像素;当与严重失真区域对应的区域中不存在严重失真像素的图像的数量为一时,利用该图像中与严重失真区域对应的区域的像素值替换严重失真区域中像素的值;当与严重失真区域对应的区域中不存在严重失真像素的图像的数量大于一时,将一个以上图像中与严重失真区域对应的区域相结合,将结合后的值替换严重失真区域中像素的值。
在步骤S706,将所有的叶片的数据相结合,得到最终图像。本步骤可采用与实施例的步骤S204相同或类似的处理方式。
本发明还提出了与上述MRI中失真校准的方法对应的装置。
图8是本发明实施例三的MRI中失真校准的装置结构图。如图8所示,该装置包括:采集单元801、几何失真校准单元802、图像结合单元803。其中,采集单元801在复数个数据读出方向上采集K空间数据,并转换为复数个图像;几何失真校准单元802生成分别对应于复数个图像的复数个像素偏移图,并根据复数个像素偏移图对复数个图像分别进行几何失真校准;图像结合单元803对复数个经过几何失真校准的图像进行结合。
几何失真校准单元802包括像素偏移图生成模块8021和几何失真校准模块8022。其中,像素偏移图生成模块8021生成分别对应于复数个图像的复数个像素偏移图;几何失真校准模块8022对复数个图像中的每个图像,根据对应的像素偏移图计算图像中每个像素的偏移位置,并根据该偏移位置以及该偏移位置上像素的值计算几何失真校准后的图像的像素值。
图像结合单元803在图像域将复数个经过几何失真校准的图像相结合。或者,图像结合单元803将复数个经过几何失真校准的图像转换到K空间中;在对应于复数个经过几何失真校准的图像的复数个K空间数据相结合;再将结合后的K空间数据转换为最终图像。
图9是本发明实施例四的MRI中失真校准的装置结构图。如图9所示,该装置包括:采集单元901、几何失真校准单元902、图像结合单元903和严重失真区域替换单元904。其中采集单元901、几何失真校准单元902和图像结合单元903与采集单元801、几何失真校准单元802、图像结合单元803的功能相同。
严重失真区域替换单元904位于几何失真校准单元902和图像结合单元903之间,确定复数个经过几何失真校准的图像中是否存在严重失真的区域,并利用复数个图像中的其他图像的与严重失真区域相对应的区域的像素值替换严重失真区域的像素值。
严重失真区域替换单元904对于复数个经过几何失真校准的图像中的每个图像,根据对应的像素偏移图计算图像保真图,图像保真图反映了均匀体模的像素偏移并重叠程度;如果图像保真图中存在像素值超过阈值的区域,确定为经过几何失真校准的图像中存在严重失真区域。
具体来说,严重失真区域替换单元904搜索其他图像中的相同区域中是否存在严重失真的像素;当与严重失真区域对应的区域中不存在严重失真像素的图像的数量为一时,利用该图像中与严重失真区域对应的区域的像素值替换严重失真区域中像素的值;当与严重失真区域对应的区域中不存在严重失真像素的图像的数量大于一时,将一个以上图像中与严重失真区域对应的区域相结合,将结合后的值替换严重失真区域中像素的值。
采用上述方法和装置,可以将运动校准和几何失真校准相结合。在长时间扫描过程中,病人有可能非常不耐烦并可能挪动身体,基于PROPELLER的技术可以校准这类伪影,提供更高的图像质量,这是因为,对于RO梯度方向单一的方案,如果在失真图像中存在严重的重叠,恢复图像的保真可大大降低,然而,在PROPELLER的采集方案中,从不同的RO方向分别采集数据,由于不同的RO方向导致不同的失真模式,这样在一个叶片的严重失真位置有可能在其他叶片上是不存在的,因此从所有叶片中进行信号选择和结合叶片将提供更佳的恢复图像。而且,本发明的方法和装置中,对各叶片的几何失真进行校准,从而在磁场不均匀等导致几何失真的因素存在的情况下,可进一步提高图像的质量。在本发明的方法中,运动失真的校准和几何失真的校准是互相独立的,相互之间没有干扰,因此,几何失真校准的效果不影响运动失真校准的效果。
图10示出了不同的RO方向导致不同重叠模式的示意图。如图10所示,假设阴影区域为感兴趣区,对于阴影区域,在方向ROm上的磁场变化比ROn大,如果RO方向处于ROm并且其极性未被正确设置,在直接逆傅立叶变换后的重建图像上将出现严重的重叠。然而,这种重叠在ROn方向上不会出现。因此,我们可以利用叶片ROn而不是ROm来初步恢复阴影区域的信息。
以下是使用MATLAB进行的仿真,以示出本发明的优点。
图11是对本发明进行仿真的流程图。如图11所示,按照以下方式进行仿真:
在步骤S1101,创建体模和磁场不均匀分布图。
图12是创建的体模的原始图像,其大小是256×256。图13是用于失真模拟的磁场不均匀分布图。图13的磁场不均匀分布图被正则化,使一个单位对应一个像素的偏移,根据场不均匀分布图可计算像素偏移图。
在步骤1102,确定RO极性,对PROPELLER采集进行模拟。
模拟的采集方案为:9个叶片,每叶片有64个PE线,RO方向的角度旋转步进为π/9。图14是对基于TSE序列的模拟原始数据直接逆傅立叶变换后的所有叶片的失真图像。
在步骤1103,使用磁场图进行失真校准。
图15是使用磁场图对所有叶片进行失真校准后的图像。
在步骤1104,恢复严重失真区域。
图16是根据保真图对所有叶片中严重失真的区域进行修复的图像。
在步骤1105,结合所有叶片数据,形成最终图像。
图17是原始无失真图像、采用PROPELLER采集并经校准后的图像、采用PROPELLER采集但未经校准的图像的对比示意图。图17的左图是原始无失真图像;中图为采用PROPELLER采集并经几何失真校准后的图像;右图为采用PROPELLER采集但未经几何失真校准的图像。从图16可以看出,采用PROPELLER并经几何失真校准后,图像质量较未经几何失真校准的图像的质量大幅提高。说明采用了本发明的方法和装置,可以有效地改善PROPELLER采集的图像中存在的几何失真问题。从而,应用了本发明的方法和装置,既可以对MRI进行有效的运动失真校准,也可以进行有效的几何失真校准。
需要说明的是,本发明不仅可以应用在上述实施例举出的PROPELLER类的技术中,还可以应用在其他在复数个RO向上采集数据的技术中。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (14)
1、一种磁共振成像中失真校准的方法,其特征在于,包括:
在复数个数据读出方向上采集K空间数据,并转换为复数个图像;
为所述每个图像分别生成对应的像素偏移图;
根据所述像素偏移图对对应的图像进行几何失真校准;
对所有经过几何失真校准的图像进行结合。
2、根据权利要求1所述的方法,其特征在于,所述根据像素偏移图对对应的图像进行几何失真校准包括:
根据像素偏移图计算对应图像中每个像素的偏移位置;
根据该偏移位置以及该偏移位置上像素的值计算几何失真校准后的图像的像素值。
3、根据权利要求1所述的方法,其特征在于,还包括:
在对所有图像进行几何失真校准之后,确定所有经过几何失真校准的图像中是否存在严重失真的区域;
利用其他图像中与所述严重失真区域相对应的区域的像素值替换所述严重失真区域的像素值。
4、根据权利要求3所述的方法,其特征在于,所述确定所有经过几何失真校准的图像中是否存在严重失真的区域包括:
对于所有经过几何失真校准的图像中的每个图像,
根据像素偏移图计算图像保真图,所述图像保真图反映均匀体模的像素偏移并重叠的程度;
当所述图像保真图中存在像素值超过阈值的区域时,确定为所述经过几何失真校准的图像中存在严重失真区域。
5、根据权利要求3所述的方法,其特征在于,所述利用其他图像中与所述严重失真区域对应的区域的像素值替换所述严重失真区域的像素值包括:
搜索其他图像中与所述严重失真区域对应的区域中是否存在严重失真的像素;
当与所述严重失真区域对应的区域中不存在严重失真像素的图像的数量为一时,利用该图像中与所述严重失真区域对应的区域的像素值替换所述严重失真区域中像素的值;
当与所述严重失真区域对应的区域中不存在严重失真像素的图像的数量大于一时,将所述一个以上图像中与所述严重失真区域对应的区域相结合,将结合后的值替换所述严重失真区域中像素的值。
6、根据权利要求1所述的方法,其特征在于,所述对所有经过几何失真校准的图像进行结合包括:
在图像域对所有经过几何失真校准的图像相结合。
7、根据权利要求1所述的方法,其特征在于,所述对所有经过几何失真校准的图像进行结合包括:
将所有经过几何失真校准的图像转换到K空间中;
将对应于所有经过几何失真校准的图像的K空间数据相结合;
将结合后的K空间数据转换为最终图像。
8、一种磁共振成像中失真校准的装置,其特征在于,包括:
采集单元(801,901),用于在复数个数据读出方向上采集K空间数据,并转换为复数个图像;
几何失真校准单元(802,902),用于为每个图像生成对应的像素偏移图,并根据所述像素偏移图对所有图像分别进行几何失真校准;
图像结合单元(803,903),用于对所有经过几何失真校准的图像进行结合。
9、根据权利要求8述的装置,其特征在于,所述几何失真校准单元(802,902)包括:
像素偏移图生成模块(8021,9021),用于为每个图像生成对应的像素偏移图;
几何失真校准模块(8022,9022),用于根据像素偏移图计算对应图像中每个像素的偏移位置,并根据该偏移位置以及该偏移位置上像素的值计算几何失真校准后的图像的像素值。
10、根据权利要求8所述的装置,其特征在于,还包括:
严重失真区域替换单元(904),位于所述几何失真校准单元(902)和图像结合单元(903)之间,用于确定所有经过几何失真校准的图像中是否存在严重失真的区域,并利用其他图像中与所述严重失真区域相对应的区域的像素值替换所述严重失真区域的像素值。
11、根据权利要10所述的装置,其特征在于,所述严重失真区域替换单元(904)对于所有经过几何失真校准的图像中的每个图像,根据对应的像素偏移图计算图像保真图,所述图像保真图反映均匀体模的像素偏移并重叠的程度;当所述图像保真图中存在像素值超过阈值的区域时,确定为所述经过几何失真校准的图像中存在严重失真区域。
12、根据权利要求11所述的装置,其特征在于,所述严重失真区域替换单元(904)搜索其他图像中与所述严重失真区域对应的区域中是否存在严重失真的像素;当与所述严重失真区域对应的区域中不存在严重失真像素的图像的数量为一时,利用该图像中与所述严重失真区域对应的区域的像素值替换所述严重失真区域中像素的值;当与所述严重失真区域对应的区域中不存在严重失真像素的图像的数量大于一时,将所述一个以上图像中与所述严重失真区域对应的区域相结合,将结合后的值替换所述严重失真区域中像素的值。
13、根据权利要求8所述的装置,其特征在于,所述图像结合单元(803,903)在图像域将所有经过几何失真校准的图像相结合。
14、根据权利要求8所述的装置,其特征在于,所述图像结合单元(803,903)将所有经过几何失真校准的图像转换到K空间中;将对应于所有经过几何失真校准的图像的复数个K空间数据相结合;再将结合后的K空间数据转换为最终图像。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007103016660A CN101470179B (zh) | 2007-12-29 | 2007-12-29 | 磁共振成像中失真校准的方法和装置 |
US12/344,900 US8131047B2 (en) | 2007-12-29 | 2008-12-29 | Method and device for distortion correction in magnetic resonance imaging |
JP2009000374A JP5486811B2 (ja) | 2007-12-29 | 2009-01-05 | 磁気共鳴イメージングにおける歪み補正方法及び装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007103016660A CN101470179B (zh) | 2007-12-29 | 2007-12-29 | 磁共振成像中失真校准的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101470179A true CN101470179A (zh) | 2009-07-01 |
CN101470179B CN101470179B (zh) | 2012-06-27 |
Family
ID=40798521
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007103016660A Active CN101470179B (zh) | 2007-12-29 | 2007-12-29 | 磁共振成像中失真校准的方法和装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8131047B2 (zh) |
JP (1) | JP5486811B2 (zh) |
CN (1) | CN101470179B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102144923A (zh) * | 2010-02-04 | 2011-08-10 | 西门子公司 | 用于减小扩散成像中的失真的方法和磁共振设备 |
CN102279375A (zh) * | 2010-03-31 | 2011-12-14 | 西门子公司 | 用于降低在扩散成像中的失真的方法和磁共振设备 |
CN102998642A (zh) * | 2011-09-07 | 2013-03-27 | 西门子公司 | 磁共振设备内检查对象的部分区域的成像方法 |
CN103181763A (zh) * | 2011-12-29 | 2013-07-03 | 上海联影医疗科技有限公司 | 磁共振成像变形的自适应矫正方法 |
CN103576116A (zh) * | 2012-08-07 | 2014-02-12 | 通用电气公司 | 磁共振成像中的运动校正的多激励叶片获取 |
CN105676155A (zh) * | 2014-11-19 | 2016-06-15 | 西门子(深圳)磁共振有限公司 | 磁共振成像系统的磁场不均匀性值获取方法和失真校正方法 |
CN106574954A (zh) * | 2014-08-22 | 2017-04-19 | 皇家飞利浦有限公司 | 针对epi的具有奈奎斯特伪影校正的并行mr成像 |
CN106680748A (zh) * | 2016-12-30 | 2017-05-17 | 上海联影医疗科技有限公司 | 磁共振图像获取方法及磁共振系统 |
CN106842084A (zh) * | 2016-12-30 | 2017-06-13 | 上海联影医疗科技有限公司 | 一种磁共振成像方法及装置 |
CN106999091A (zh) * | 2014-12-04 | 2017-08-01 | 通用电气公司 | 用于组成材料的改进分类的方法和系统 |
CN109073721A (zh) * | 2016-04-28 | 2018-12-21 | 皇家飞利浦有限公司 | 处置计划评价工具 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7541808B2 (en) * | 2007-04-11 | 2009-06-02 | Allegheny-Singer Research Institute | Rapid MRI dynamic imaging using MACH |
US8405394B2 (en) * | 2009-10-20 | 2013-03-26 | Allegheny-Singer Research Institute | Targeted acquisition using holistic ordering (TACHO) approach for high signal to noise imaging |
CN102841329B (zh) * | 2011-06-24 | 2016-03-30 | 西门子(深圳)磁共振有限公司 | 磁共振信号处理方法及装置 |
US9581671B2 (en) * | 2014-02-27 | 2017-02-28 | Toshiba Medical Systems Corporation | Magnetic resonance imaging with consistent geometries |
EP3114495B1 (en) * | 2014-03-05 | 2020-10-21 | Koninklijke Philips N.V. | Mri propeller with motion correction, water-fat separation and estimation of magnetic field inhomogeneity information |
US10180483B2 (en) | 2014-04-24 | 2019-01-15 | David W Holdsworth | 3D printed physical phantom utilized in measuring 3D geometric distortion occurring in MRI and CT images |
WO2017009391A1 (en) | 2015-07-15 | 2017-01-19 | Koninklijke Philips N.V. | Mr imaging with motion detection |
CN113962876A (zh) * | 2019-07-17 | 2022-01-21 | 中国电子科技集团公司第十三研究所 | 一种像素畸变的校正方法、校正装置及终端 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62179449A (ja) * | 1986-01-31 | 1987-08-06 | 横河メディカルシステム株式会社 | 核磁気共鳴画像の歪補正装置 |
JPS6321048A (ja) * | 1986-07-16 | 1988-01-28 | 横河メディカルシステム株式会社 | 核磁気共鳴画像の歪み補正方法 |
DE69326202T2 (de) * | 1992-06-29 | 2000-04-13 | Koninkl Philips Electronics Nv | Verfahren und Gerät zur Kernresonanzabbildung |
DE19540837B4 (de) * | 1995-10-30 | 2004-09-23 | Siemens Ag | Verfahren zur Verzeichnungskorrektur für Gradienten-Nichtlinearitäten bei Kernspintomographiegeräten |
US5943433A (en) * | 1996-12-30 | 1999-08-24 | General Electric Company | Method for correcting inhomogeneity of spatial intensity in an aquired MR image |
US6201985B1 (en) * | 1998-08-14 | 2001-03-13 | General Electric Company | Segmented k-space method for three-dimensional MR imaging |
US6794869B2 (en) * | 2001-03-30 | 2004-09-21 | General Electric Company | Moving table MRI with frequency-encoding in the z-direction |
EP1554595B1 (en) * | 2001-12-17 | 2008-05-14 | Koninklijke Philips Electronics N.V. | Magnetic resonance method |
JP2005152175A (ja) * | 2003-11-25 | 2005-06-16 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴撮像装置および方法 |
JP4193690B2 (ja) * | 2003-12-25 | 2008-12-10 | 株式会社日立製作所 | Mri装置 |
JP4612000B2 (ja) * | 2005-01-20 | 2011-01-12 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
US7023207B1 (en) * | 2005-02-16 | 2006-04-04 | General Electric Company | Method and system of MR imaging with reduced radial ripple artifacts |
JP2006304955A (ja) * | 2005-04-27 | 2006-11-09 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置 |
CN101232845B (zh) * | 2005-07-27 | 2010-08-04 | 株式会社日立医药 | 磁共振成像装置 |
CN100586181C (zh) * | 2006-01-18 | 2010-01-27 | 华中科技大学 | 一种图像数据水印信息的嵌入、检测方法及系统 |
US7535222B2 (en) * | 2007-01-02 | 2009-05-19 | The Board Of Trustees Of The Leland Stanford Junior University | MRI data acquisition using propeller k-space data acquisition |
DE102008007048B4 (de) * | 2008-01-31 | 2010-06-17 | Siemens Aktiengesellschaft | Dynamische Verzeichnungskorrektur bei EPI-Messungen in der medizinischen Magnet-Resonanz-Bildgebung |
-
2007
- 2007-12-29 CN CN2007103016660A patent/CN101470179B/zh active Active
-
2008
- 2008-12-29 US US12/344,900 patent/US8131047B2/en active Active
-
2009
- 2009-01-05 JP JP2009000374A patent/JP5486811B2/ja active Active
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102144923B (zh) * | 2010-02-04 | 2014-12-17 | 西门子公司 | 用于减小扩散成像中的失真的方法和磁共振设备 |
CN102144923A (zh) * | 2010-02-04 | 2011-08-10 | 西门子公司 | 用于减小扩散成像中的失真的方法和磁共振设备 |
CN102279375A (zh) * | 2010-03-31 | 2011-12-14 | 西门子公司 | 用于降低在扩散成像中的失真的方法和磁共振设备 |
CN102279375B (zh) * | 2010-03-31 | 2015-08-19 | 西门子公司 | 用于降低在扩散成像中的失真的方法和磁共振设备 |
CN102998642A (zh) * | 2011-09-07 | 2013-03-27 | 西门子公司 | 磁共振设备内检查对象的部分区域的成像方法 |
CN102998642B (zh) * | 2011-09-07 | 2016-12-21 | 西门子公司 | 磁共振设备内检查对象的部分区域的成像方法 |
CN103181763A (zh) * | 2011-12-29 | 2013-07-03 | 上海联影医疗科技有限公司 | 磁共振成像变形的自适应矫正方法 |
CN103576116B (zh) * | 2012-08-07 | 2018-01-09 | 通用电气公司 | 磁共振成像中的运动校正的多激励叶片获取 |
CN103576116A (zh) * | 2012-08-07 | 2014-02-12 | 通用电气公司 | 磁共振成像中的运动校正的多激励叶片获取 |
CN106574954B (zh) * | 2014-08-22 | 2020-04-07 | 皇家飞利浦有限公司 | 针对epi的具有奈奎斯特伪影校正的并行mr成像 |
CN106574954A (zh) * | 2014-08-22 | 2017-04-19 | 皇家飞利浦有限公司 | 针对epi的具有奈奎斯特伪影校正的并行mr成像 |
CN105676155B (zh) * | 2014-11-19 | 2019-02-22 | 西门子(深圳)磁共振有限公司 | 磁共振成像系统的磁场不均匀性值获取方法和装置和失真校正方法和装置 |
CN105676155A (zh) * | 2014-11-19 | 2016-06-15 | 西门子(深圳)磁共振有限公司 | 磁共振成像系统的磁场不均匀性值获取方法和失真校正方法 |
CN106999091A (zh) * | 2014-12-04 | 2017-08-01 | 通用电气公司 | 用于组成材料的改进分类的方法和系统 |
CN106999091B (zh) * | 2014-12-04 | 2020-12-01 | 通用电气公司 | 用于组成材料的改进分类的方法和系统 |
CN109073721A (zh) * | 2016-04-28 | 2018-12-21 | 皇家飞利浦有限公司 | 处置计划评价工具 |
CN109073721B (zh) * | 2016-04-28 | 2022-02-08 | 皇家飞利浦有限公司 | 处置计划评价工具 |
CN106842084A (zh) * | 2016-12-30 | 2017-06-13 | 上海联影医疗科技有限公司 | 一种磁共振成像方法及装置 |
CN106680748A (zh) * | 2016-12-30 | 2017-05-17 | 上海联影医疗科技有限公司 | 磁共振图像获取方法及磁共振系统 |
CN106680748B (zh) * | 2016-12-30 | 2019-10-08 | 上海联影医疗科技有限公司 | 磁共振图像获取方法及磁共振系统 |
CN106842084B (zh) * | 2016-12-30 | 2019-11-12 | 上海联影医疗科技有限公司 | 一种磁共振成像方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
US8131047B2 (en) | 2012-03-06 |
JP5486811B2 (ja) | 2014-05-07 |
CN101470179B (zh) | 2012-06-27 |
JP2009160409A (ja) | 2009-07-23 |
US20090169083A1 (en) | 2009-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101470179B (zh) | 磁共振成像中失真校准的方法和装置 | |
CN101711126B (zh) | 磁共振成像装置及由倾斜磁场引起的误差修正方法 | |
US8344729B2 (en) | Magnetic resonance imaging apparatus | |
JP6691797B2 (ja) | 磁気共鳴イメージング装置 | |
US6448771B1 (en) | Magnetic resonance method for forming a fast dynamic image | |
CN101843487B (zh) | 借助并行采集技术建立图像的方法 | |
CN101680937B (zh) | 磁共振设备和方法 | |
JP3929047B2 (ja) | 磁気共鳴イメージング装置 | |
CN101435858B (zh) | 基于部分并行获取进行磁共振成像的方法和装置 | |
CN101669038A (zh) | 用于propeller mri的磁共振设备和方法 | |
Anderson III et al. | Adaptive retrospective correction of motion artifacts in cranial MRI with multicoil three‐dimensional radial acquisitions | |
JP2000229073A (ja) | 磁気共鳴画像装置 | |
CN101672905A (zh) | 在患者卧榻连续移动情况下的失真校正方法和装置 | |
CN101498774A (zh) | 医学磁共振成像回波平面成像测量中的动态失真校正方法 | |
CN101105523A (zh) | 在磁共振断层造影中加速自旋编码成像的方法和设备 | |
WO2022236851A1 (zh) | 卷褶视野磁共振图像的重建方法、计算机设备及存储介质 | |
CN101082659A (zh) | 用接收器线圈阵列改进灵敏编码磁共振成像的方法和设备 | |
CN100443917C (zh) | 磁共振方法和设备 | |
CN101236240B (zh) | 磁共振断层造影成像中改进的三维层选择的多层激励方法 | |
Mickevicius et al. | Simultaneous acquisition of orthogonal plane cine imaging and isotropic 4D-MRI using super-resolution | |
CN102841329B (zh) | 磁共振信号处理方法及装置 | |
JP5722212B2 (ja) | 磁気共鳴イメージング装置及び方法 | |
JP2020175194A (ja) | 磁気共鳴イメージング装置及び磁気共鳴イメージング方法 | |
US20050253580A1 (en) | Method of pseudopolar acquisition and reconstruction for dynamic MRI | |
Bao et al. | Unsupervised deep learning model for correcting Nyquist ghosts of single‐shot spatiotemporal encoding |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210114 Address after: Room 516, 5th floor, 38 Yinglun Road, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai 200031 Patentee after: SIEMENS HEALTINEERS Ltd. Address before: 100102 No. 7 South Central Road, Chaoyang District, Beijing, China Patentee before: Siemens (China) Co.,Ltd. |