CN102438518A - 磁共振成像装置及磁共振成像方法 - Google Patents

磁共振成像装置及磁共振成像方法 Download PDF

Info

Publication number
CN102438518A
CN102438518A CN2011800017350A CN201180001735A CN102438518A CN 102438518 A CN102438518 A CN 102438518A CN 2011800017350 A CN2011800017350 A CN 2011800017350A CN 201180001735 A CN201180001735 A CN 201180001735A CN 102438518 A CN102438518 A CN 102438518A
Authority
CN
China
Prior art keywords
echo
data
signal
collection
template data
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
Application number
CN2011800017350A
Other languages
English (en)
Other versions
CN102438518B (zh
Inventor
草原博志
葛西由守
待井丰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of CN102438518A publication Critical patent/CN102438518A/zh
Application granted granted Critical
Publication of CN102438518B publication Critical patent/CN102438518B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615Echo 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/5616Echo 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56554Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI

Landscapes

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

Abstract

本发明提供一种磁共振成像装置及磁共振成像方法。在一个实施方式中,MRI装置(20)具有第1采集部、第2采集部和校正部(100)。第1采集部及第2采集部通过执行含有相位编码方向梯度磁场的施加的EPI的回波信号采集序列,产生多个回波信号,将多个回波信号作为第1、第2模板数据分别进行采集。第2采集部在第1模板数据的采集之后,按照读取方向梯度磁场的施加开始时刻与第1模板数据的采集时错开的序列,采集第2模板数据。校正部(100)用第1及第2模板数据对包含在回波信号中的相位误差进行校正。

Description

磁共振成像装置及磁共振成像方法
技术领域
本发明的实施方式涉及磁共振成像装置及磁共振成像方法。
背景技术
MRI是以拉莫尔频率的RF脉冲对置于静磁场中的受检体的原子核自旋进行磁激励,从伴随该激励而产生的MR信号重建图像的摄像法。再者,上述MRI是磁共振成像(Magnetic Resonance Imaging)的意思,RF脉冲是高频脉冲(radio frequency pulse)的意思,MR信号是核磁共振信号(nuclear magnetic resonance signal)的意思。
在MRI中,有一种被称为EPI(echo planar imaging:回波平面成像)的摄像法。EPI是高速摄像法的一种,它是针对1次的核磁激励使梯度磁场(倾斜磁場)高速连续翻转(反転),连续生成回波而采集MR信号的方法。
具体而言,在EPI中,在施加激励脉冲后、xy平面内的磁化经过横向弛豫(T2弛豫)过程衰减而丧失之前,利用相位编码的阶跃(ステツプ)发生连续的梯度回波(グラジエントエコ—),采集重建图像所必需的所有数据。
EPI中,有采集在激励脉冲以及重聚脉冲之后产生的自旋回波信号的利用了自旋回波法(SE:spin echo)的SE EPI、采集发生在施加激励脉冲之后产生的回波信号的利用了场回波法(FE:field echo)的FE EPI、利用了FFE(Fast FE)法的FFE EPI。
此外,把经多次施加激励脉冲而得到的回波链的数据合并而生成1份图像数据的EPI被称为多次激发(シヨツト)EPI,而仅以1次施加的激励脉冲重建图像的EPI被称为单次激发EPI。
由于在EPI序列中一边使梯度磁场高速翻转一边摄像,所以导致了在所采集的回波数据中含有相位误差。从而在基于这样的回波数据重建的图像中发生畸变。EPI中的相位误差的主要原因可以考虑为静磁场不均匀性和因梯度磁场的切换(スイツチング)而产生的涡流磁场这两个原因(参照,例如,非专利文献1及非专利文献2)。
上述由于静磁场不均匀性引起的相位误差因摄像对象的信号强度分布和摄像区域内的静磁场不均匀性的空间分布而异。另一方面,上述由于涡流磁场引起的相位误差主要在实际空间的引出(リ—ドアウト)方向上具有1次梯度,由于回波数据采集中的引出方向梯度磁场的极性,其相位梯度的方向发生翻转。
作为减少由于静磁场不均匀性引起的相位误差的现有技术,有专利文献1中所记载的静磁场不均匀性的校正方法等。
另一方面,作为减少由于静磁场不均匀性以外的原因引起的相位误差的方法,专利文献2中所记载的技术为已知。
具体而言,在专利文献2中,在正式摄像前采集使引出方向梯度磁场的极性翻转后的2次模板激发(テンプレ—トシヨツト)A、B的回波数据。并且,由于通过模板激发A、B分别采集的一对回波数据具有相同的回波时间,所以静磁场不均匀性引起的相位错开的值相等,基于这一点,消除了由于静磁场不均匀性引起的相位误差成分。由此,仅抽取静磁场不均匀性以外的原因引起的相位误差成分,以此作为相位校正数据。
(专利文献)
专利文献1:日本特开2006-255046号公报
专利文献2:日本特开平9-276243号公报
(非专利文献)
非专利文献1:Self-Correcting EPI Reconstruction Algorithm;A.Jesmanowicz.et.al.;Proceedings of SMR,No.619,1995
非专利文献2:Phase Correction for EPI Using Internal ReferenceLine;
A.Jesmanowicz.et.al.;Proceedings of SMR,No.1239,1995
发明内容
对回波数据所包含的相位误差中起因于静磁场不均匀性的相位误差和由其他原因引起的相位误差进行校正的方法如上所述有多种多样。虽然利用现有技术的EPI在实用上已经得到了十分良好的图像,但是还是希望起因于相位误差的图像畸变尽可能地小。
本发明所要解决的问题是,提供用于进一步减少起因于EPI中的相位误差的图像畸变的、与现有技术不同的技术。
在一个实施方式中,MRI装置能够执行EPI,即:向静磁场内的受检体发射用于引发核磁共振的激励脉冲,采集反复进行读取方向梯度磁场的极性翻转所产生的多个回波信号,基于多个所述回波信号重建图像数据,该MRI装置具有第1采集部、第2采集部和校正部。
第1采集部通过执行包含相位编码方向梯度磁场的施加的EPI的回波信号采集序列,从而将所产生的多个回波信号作为第1模板数据进行采集。
该第2采集部在第1模板数据的采集后,通过执行EPI的回波信号采集序列,将所产生的多个回波信号作为第2模板数据进行采集,其中,EPI的回波信号采集序列包含相位编码方向梯度磁场的施加,并且读取方向梯度磁场的施加开始时刻与第1模板数据的采集时相错开。
校正部,用所述第1及第2模板数据,至少进行包含在回波信号中的相位误差的校正。
根据上述构成的MRI装置,利用与现有技术不同的技术,能够进一步减少起因于EPI中的相位误差的图像畸变。
在一个实施方式中,MRI方法为包括EPI的方法,具有以下的步骤。
一个是通过执行包含相位编码方向梯度磁场的施加的EPI的回波信号采集序列,从而将所产生的多个回波信号作为第1模板数据进行采集的步骤。
一个是在第1模板数据的采集后,通过执行EPI的回波信号采集序列,从而将所产生的多个回波信号作为第2模板数据进行采集的步骤,其中,EPI的回波信号采集序列包含相位编码方向梯度磁场的施加,并且读取方向梯度磁场的施加开始时刻与第1模板数据的采集时相错开。
一个是用第1及第2模板数据,至少进行包含在回波信号中的相位误差的校正的步骤。
根据上述构成的MRI方法,利用与现有技术不同的技术,能够进一步减少起因于EPI中的相位误差的图像畸变。
附图说明
图1是表示一个实施方式中的MRI装置的整体构成的框图。
图2是图1所示计算机58的功能框图。
图3是表示在第1实施方式中,作为自旋回波系单次激发EPI的主扫描、模板激发1、2的脉冲序列之一例的时间图。
图4是表示相位编码及频率编码的矩阵元素数为256×256的场合,在k空间上的MR信号的排列顺序之一例的示意图。
图5是表示从模板数据1得到的第1相位图像之一例的示意图。
图6是表示从模板数据2得到的第2相位图像之一例的示意图。
图7是表示用第1相位图像和第2相位图像,基于公式(3)得到的差分图像(磁场校正映射图)之一例的示意图。
图8是不在主扫描之前进行利用磁场校正映射图的静磁场不均匀性的校正,对均匀的模型(フアントム)进行摄像而得到的图像之一例。
图9是在主扫描之前进行利用磁场校正映射图的静磁场不均匀性的校正后,对与图8相同的模型进行摄像而得到的图像之一例。
图10是通过抽取出图8中的4处畸变区域,表示不进行基于磁场校正映射图的校正的场合(图8)与进行基于磁场校正映射图的校正的场合(图9)之不同的示意图。
图11是表示第1实施方式中的MRI装置的动作流程的流程图。
图12是表示在第2实施方式中,作为自旋回波系单次激发EPI的主扫描、模板激发1、2’的脉冲序列之一例的时间图。
图13是表示第2实施方式中的MRI装置的动作流程的流程图。
图14是表示在第3实施方式中,作为自旋回波系单次激发EPI的主扫描、模板激发1、2’、3的脉冲序列之一例的时间图。
图15是表示第3实施方式中的MRI装置的动作流程的流程图。
图16是表示作为FID型的单次激发EPI的主扫描、模板激发1”、2”的脉冲序列之一例的时间图。
具体实施方式
以下,基于附图说明MRI装置以及MRI方法的实施方式。再者,各图中,相同的要素赋予相同的附图标记,省略重复的说明。
(第1实施方式)
图1是表示第1实施方式中的MRI装置20的整体构成的框图。如图1所示,MRI装置20具有用于形成静磁场的筒状静磁场用磁铁22、在静磁场用磁铁22的内侧同轴设置的筒状匀场线圈24、梯度磁场线圈26、RF线圈28、控制装置30和载放受检体P的台板32。
在此作为一个例子,对装置坐标系的彼此垂直交叉的X轴、Y轴、Z轴以如下方式进行定义。首先,静磁场用磁铁22及匀场线圈24是以其轴向与铅垂方向垂直交叉的方式配置的,将静磁场用磁铁22及匀场线圈24的轴向设为Z轴方向。此外,设铅垂方向为Y轴方向,台板32是以其顶板的载放用的面的法线方向成为Y轴方向的方式配置的。
控制装置30具有静磁场电源40、匀场线圈电源42、梯度磁场电源44、RF发射器46、RF接收器48、序列控制器56和计算机58。
梯度磁场电源44具有X轴梯度磁场电源44x、Y轴梯度磁场电源44y和Z轴梯度磁场电源44z。此外,计算机58具有运算装置60、输入装置62、显示装置64和存储装置66。
静磁场用磁铁22与静磁场电源40连接,利用从静磁场电源40供给的电流在摄像空间中形成静磁场。
上述所谓摄像空间是指,例如,放置受检体P、被施加了静磁场的门架(ガントリ)内的空间。所谓门架,是以含有静磁场磁铁22、匀场线圈24、梯度磁场线圈26和RF线圈28的方式形成为例如圆筒状的构造体。门架及台板32以载有受检体P的台板32能够在门架内部移动的方式构成。再者,图1中由于烦杂,所以没有对门架自身进行图示,而是将门架内的静磁场磁铁22、匀场线圈24、梯度磁场线圈26和FR线圈28作为构成要素进行了图示。
摄像区域是指,例如,生成“一个图像”或者“一套图像”用的MR信号的采集范围,是作为摄像空间的一部分被设定的区域。这里的所谓“一个图像”及“一套图像”既可以是二维图像,也可以是三维图像。这里的所谓“一套图像”是指,例如,像多薄片摄像等那样,在一个脉冲序列内“多个图像”的MR信号被一并采集场合的“多个图像”。摄像区域由,例如,装置坐标系限定为三维区域。这里作为一个例子,摄像区域如果是厚度较薄的区域则称之为薄片(スライス),如果是有一定程度的厚度的区域则称之为厚片(スラブ)。
匀场线圈24与匀场线圈电源42连接,利用从匀场线圈电源42供给的电流使该静磁场均匀化。静磁场用磁铁22大多数场合用超导线圈构成,一般情况下在励磁时连接到静磁场电源40而被供给电流,一经励磁后则被置于非连接状态。再者,不设置静磁场电源40,用永磁铁构成静磁场用磁铁22也是可以的。
梯度磁场线圈26具有X轴梯度磁场线圈26x、Y轴梯度磁场线圈26y和Z轴梯度磁场线圈26z,其在静磁场用磁铁22的内侧形成为筒状。X轴梯度磁场线圈26x、Y轴梯度磁场线圈26y和Z轴梯度磁场线圈26z分别与X轴梯度磁场电源44x、Y轴梯度磁场电源44y和Z轴梯度磁场电源44z连接。
利用分别从X轴梯度磁场电源44x、Y轴梯度磁场电源44y和Z轴梯度磁场电源44z供给到X轴梯度磁场线圈26x、Y轴梯度磁场线圈26y和Z轴梯度磁场线圈26z的电流,在摄像区域中分别形成X轴方向梯度磁场Gx、Y轴方向梯度磁场Gy和Z轴方向梯度磁场Gz。
即,将装置坐标系的X、Y、Z轴方向的梯度磁场Gx、Gy、Gz进行合成,能够任意设定作为理论轴的薄片选择方向梯度磁场Gss、相位编码方向梯度磁场Gpe以及读取方向(频率编码方向)梯度磁场Gro。薄片选择方向、相位编码方向以及读取方向的各梯度磁场Gss、Gpe、Gro与静磁场重叠。
RF发射器46基于从序列控制器56输入的控制信息,生成用于引发核磁共振的拉莫尔频率的RF脉冲(RF电流脉冲),将该RF脉冲发射到发射用的RF线圈28。RF线圈28有内置于门架的发射、接收RF脉冲的用的全身用线圈和设在台板32或受检体P附近的接收RF脉冲用的局部线圈等。
发射用的RF线圈28从RF发射器46接收RF脉冲后将其发射到受检体P。接收用的RF线圈28接收因受检体P内部的原子核自旋受到RF脉冲激励而产生的MR信号(高频信号),该MR信号由RF接收器48进行检测。
RF接收器48在对所检测到的MR信号进行了前置增幅、中间频率变换、相位检波、低频增幅、过滤等各种信号处理之后,通过进行A/D(模拟(analog)-数字(digital))变换,生成数字化的复变数据的原始数据(raw data)。RF接收器48将所生成的MR信号的原始数据输入到序列控制器56。
运算装置60是进行MRI装置20整体系统控制的装置。
序列控制器56根据运算装置60的指令存储为驱动梯度磁场电源44、RF发射器46及RF接收器48所必需的控制信息。这里的所谓控制信息,是记述了,例如,应该施加在梯度磁场电源44上的脉冲电流的强度、施加时间、施加时刻等动作控制信息的序列信息。
序列控制器56根据所存储的规定序列对梯度磁场电源44、RF发射器46及RF接收器48进行驱动,由此发生X轴梯度磁场Gx、Y轴梯度磁场Gy、Z轴梯度磁场Gz以及RF脉冲。此外,序列控制器56接收从RF接收器48输入的MR信号的原始数据(raw data)并将其输入到运算装置60。
图2是图1所示计算机58的功能模块图。计算机的运算装置60具有MPU(Micro Processor Unit)86、系统总线88、图像重建部90、图像数据库94、图像处理部96、显示控制部98和校正部100。
MPU86在摄像条件的设定、摄像动作以及摄像后的图像显示中,通过系统总线88等的配线进行MRI装置20整体的系统控制。
所谓“摄像条件”,是指,例如,由自旋回波或EPI等中的哪一种脉冲序列,在什么样的条件下发射RF信号等,在什么样的条件下从受检体采集MR信号的意思。
作为“摄像条件”的各个项目,可以列举出,例如,作为摄像空间内的位置性信息的摄像区域、翻转角、重复时间(TR:RepetitionTime)、薄片数、相位编码方向以及频率编码方向的阶跃数、EPI或自旋回波法等的脉冲序列的种类等。
此外,MPU86也发挥着作为摄像条件设定部的功能,基于来自输入装置62的指令信息设定包含脉冲序列在内的摄像条件,并将所设定的摄像条件输入到序列控制器56。因此,MPU86对显示控制部98进行控制,使摄像条件的设定用画面信息显示在显示装置64中。
输入装置62向用户提供设定摄像条件和图像处理条件的功能。
图像重建部90内部具有k空间数据库92。图像重建部90将从序列控制器56输入的MR信号的原始数据作为k空间数据配置在k空间数据库92中形成的k空间里。图像重建部90对k空间数据施以图像重建处理,生成受检体P的各薄片的图像数据,并将所生成的图像数据保存在图像数据库94中。
图像处理部96从图像数据库94获取图像数据,并对其施以规定的图像处理,将图像处理后的图像数据作为显示用图像数据存储在存储装置66中。
存储装置66使上述显示用图像数据的生成所使用的摄像条件和受检体P的信息(患者信息)等作为附带信息附属于该显示用图像数据进行存储。
显示控制部98根据MPU86的控制,将摄像条件的设定用画面和摄像生成的图像数据所表示的图像显示在显示装置64上。
校正部100基于在作为目的图像用MR信号采集(主扫描)之前作为预扫描进行的模板激发1、2的采集数据,进行以下二项工作。
第1、校正部100生成用于对静磁场的不均匀性进行校正(均匀化)的磁场校正映射图。
第2、校正部100生成从通过主扫描采集的MR信号重建图像数据时用于消除MR信号中所含有的相位误差的影响的相位校正数据。以下将就这二项工作的详细情况进行说明。
图3是表示主扫描及模板激发1、2的各脉冲序列之一例的时间图。图3的上段、中段和下段分别表示主扫描(MAIN SCAN)、模板激发1(TEMPLATE SHOT 1)和模板激发2(TEMPLATE SHOT 2)的各脉冲序列,各横轴表示经过时间t。在图3的上段、中段和下段中,RF表示RF脉冲,Gss表示薄片选择方向的梯度磁场,Gpe表示相位编码方向梯度磁场,Gro表示读取方向梯度磁场,Signal表示回波(MR信号)。再者,作为本实施方式的特征之一举出的是在模板激发1、2中施加有相位编码方向梯度磁场Gpe这一点。
在图3的主扫描中作为一个例子表示了自旋回波系的单次激发EPI。即,在施加了翻转角为90°的激励脉冲(RF脉冲)后,在经过了有效编码时间的大约一半的时刻施加180°的激励脉冲,采集MR信号。这里作为一个例子,模板激发1用与主扫描相同的脉冲序列。
模板激发1、2在脉冲序列上的不同之处仅为下面的2点。
作为第1个不同点,在模板激发2中,使相位编码阶跃脉冲的“施加开始时刻”以及读取方向梯度磁场Gro的施加开始时刻比模板激发1延迟回波间隔2倍的时间差。
这里的所谓“施加开始时刻”,是指以,例如,180°RF脉冲的施加开始时刻为基准的时刻。再者,这里作为一个例子,关于主扫描以及各模板激发1、2、90°脉冲、180°脉冲的施加开始时刻是共通的(下文中的第2实施方式的模板激发2’、第3实施方式的模板激发3亦同)。
此外,上述所谓“回波间隔”,是指从读取方向梯度磁场Gro的极性翻转开始到再次翻转为止的时间范围,以下记为“Gro翻转间隔”。由此,模板激发2中的、以施加90°激励脉冲时为基准的回波(MR信号)的发生开始时刻,比模板激发1的场合延迟Gro翻转间隔2倍的时间差。
此外,上述所谓“相位编码阶跃脉冲”是在图3的相位编码方向梯度磁场Gpe的例子中,除了最初被施加了的预脉冲之外的其它的相位编码方向梯度磁场脉冲,其以Gro翻转间隔施加。预脉冲是在图3中横轴(时间轴)下侧以梯形形状表示的部分。关于相位编码方向梯度磁场Gpe,只有最初施加的预脉冲与其它的不同,第2个以后以规定的时间间隔施加的相位编码阶跃脉冲是彼此相同的。
通过使这些相位编码阶跃脉冲与读取方向梯度磁场的翻转同步施加,相位编码方向的位置性信息被赋予了伴随梯度磁场的翻转连续产生的各MR信号中。
作为第2个不同点,相位编码方向梯度磁场Gpe的预脉冲的施加强度的时间积分值的绝对值在模板激发1、2中不同(二者的施加时间相同)。
具体而言,k空间数据的中心线上排列有在相位编码阶跃为零的时刻接收的MR信号。所谓相位编码阶跃为零是指在单次激发的EPI场合,以规定间隔依次施加的相位编码阶跃脉冲的总计面积等于相位编码方向梯度磁场Gpe的预脉冲的面积的时刻。这里的所谓面积是指信号强度的时间积分值的绝对值。此外,这里的所谓相位编码阶跃为零,与利用通常的自旋回波法变换相位编码阶跃每次采集MR信号的场合,不施加相位编码方向梯度磁场Gpe进行采集的阶跃相对应。
在模板激发1中,以下的2个时刻被定为相同时刻。其中一方是,由于根据自旋回波法施加90°及180°RF脉冲,选择薄片中的所有氢原子核的自旋相位(各自旋的横磁化矢量的方向)一致,MR信号强度应当成为最大的时刻。该时刻作为有效回波时间(EFFECTIVE ECHOTIME)在图3中以纵向的虚线表示。另外一方是,信号强度由于相位编码方向梯度磁场而应当成为最大的时刻,即,相位编码阶跃为零的时刻。
即,在模板激发1中,在相位编码方向梯度磁场Gpe中,最初的预脉冲的面积等于紧接其后的5个相位编码阶跃脉冲的总计面积。在这种场合,在图3中以虚线所示的有效回波时间里,信号强度最大的MR信号被采集,该数据被配置在k空间的中心线上。
另一方面,在模板激发2的相位编码方向梯度磁场Gpe中,最初的预脉冲的面积等于紧接其后的4个相位编码阶跃脉冲的总计面积。还有,在模板激发2中,MR信号的产生开始时间也比模板激发1延迟使读取方向梯度磁场的施加开始时刻延迟的部分——延迟2个回波的量。因此,在模板激发2中,紧跟在施加(从时间上早的一方算起)第4个相位编码阶跃脉冲后面采集的MR信号,与相位编码阶跃为零相当,被配置在k空间的中心线上(该MR信号成为最大强度。)
即,在模板激发2中,k空间中心线的MR信号(相位编码阶跃为零时刻的采集信号)与模板激发1相比,采集时刻延迟1个回波的量。为此,从模板激发2得到的相位图像与从模板激发1得到的相位图像相比,由于静磁场的不均匀性等原因,相位提前了。
因此,通过利用上述模板激发1、2的采集数据的差分,能够生成用于校正相位误差的相位校正数据以及磁场校正映像图。
这样,作为第1实施方式的原理,可以说有以下的2个条件。
作为第1个条件,出于生成磁场校正映射图的观点,在模板激发1、2中,最好使与相位编码阶跃为零相对应的MR信号的产生时间错开(例如Gro翻转间隔的自然数倍)。这里的所谓“错开”,是指以,例如,在相位编码方向梯度磁场Gpe的预脉冲的施加开始时,或者180°RF脉冲的施加开始时为基准错开。
作为第2个条件,出于生成相位校正数据的观点,使与相位编码阶跃为零的时刻的读取方向梯度磁场Gro的极性,在模板激发1、2中彼此相反。这是因为,由于使在上述时刻的读取方向梯度磁场Gro极性翻转后的2个模板激发的数据差分,静磁场的不均匀性所引起的相位误差成分被消除,能够仅抽取出相位误差成分。
因此,在不生成磁场校正映射图而生成相位校正数据的场合,无需满足上述第1条件,使模板激发1、2中与相位编码阶跃为零对应的MR信号的产生时刻相同也是可以的(参照下文的第2实施方式及图12)。
同样,在不考虑生成相位校正数据而仅着眼于生成磁场校正映射图的场合,不必满足上述第2条件。即,相位编码阶跃为零的时刻的读取方向梯度磁场Gro的极性在2个模板激发中相同也是可以的。这一点将以下文的第3实施方式的模板激发1、3和图4进行说明。
在第1实施方式中,为了生成磁场校正映射图及相位校正数据,模板激发1、2满足上述第1及第2条件。
具体而言,在图3的模板激发1的例子中,以虚线表示的有效回波时间中,相位编码阶跃变为零,这一时刻的读取方向梯度磁场Gro的极性为负。
另一方面,在模板激发2中,在仅比模板激发1的场合延迟了Gro翻转间隔(仅1个回波的量)的时刻,相位编码阶跃变为零,这一时刻的读取方向梯度磁场Gro的极性为正。
此外,虽然在图3的例子中在模板激发1、2间错开的MR信号的产生开始时刻,以Gro翻转间隔为基准,设成了它的2倍,但是使其以其它的自然数倍的时间差延迟也是可以的。这可以通过将相位编码阶跃脉冲以及读取方向梯度磁场Gro的施加开始时刻适宜地错开来实现。
在本实施方式中,用上述模板激发1、2的采集数据生成磁场校正映像图以及相位校正数据,先就相位校正数据的生成方法进行说明。
首先,进行从通过模板激发1采集的各MR信号(以下称为模板数据1)减去输送频率的余弦函数等处理,得到k空间数据的实数部分。此外,进行从模板数据1减去输送频率的正弦函数等处理,得到k空间数据的虚数部分。以下,就k空间数据的制作方法之一例用图4进行具体说明。
图4是表示相位编码及频率编码的矩阵元素数为256×256的场合,k空间中的MR信号的排列顺序之一例的示意图。图4中,TR为重复时间(Repetition Time),横向的Ts为采样时间(Sampling Time),纵向为相位编码阶跃(Phase Encode Step)。
通常的自旋回波法那样的脉冲序列,是使相位编码变换256次,从所采集的256行的MR信号分别减去输送频率的余弦函数或正弦函数。将进行了该处理后的256个MR信号,按照图4那样的相位编码阶跃顺序,自下而上排列为-127、-126、......-1、0、1、......127、128。由此得到由256×256的矩阵元素构成的矩阵数据,即k空间数据的实数部分或虚数部分。
另一方面,像本实施方式这样的单次激发的EPI,如果在相位编码阶跃为零的时刻之前只能采集4行份的MR信号,则采集数为(256/2)+4=132行。在这种场合,没被采集的124行在k空间上例如作为数据填上零。
在模板激发1的场合,在相位编码阶跃为-127~-6的各行中,作为MR信号的数据,填上例如零(参照图4)。并且,通过模板激发1采集的MR信号,如图3所示,相位编码阶跃成为按照时间顺序的-5、-4、-3、-2、-1、0、1、......127、128的顺序。即,在模板激发1中,在最早时间被采集的MR信号在k空间中被配置在相位编码阶跃为-5的行上。此外,在模板激发1中,早在第6时间被采集的MR信号在k空间中被配置在相位编码阶跃为零的行(k空间的中心线)上。
这里,在图4的横向,例如,把各MR信号的采样时间Ts以256按相等间隔分割成ΔTs,把在每个ΔTs的MR信号的强度作为矩阵值。由此,分别关于实数和虚数求256行256列的矩阵数据,以此作为k空间数据。
接下来,通过对从模板激发1得到的k空间数据的实数部分的中心线的数据施以1维傅立叶逆变换,从而得到以横轴为频率、以纵轴为光谱强度的数据。接下来,使该1维傅立叶逆变换后的数据从-∞变成+∞(单位为赫兹或弧度)为止,将频率积分值作为Real1。
同样,对从模板激发1得到的k空间数据的虚数部分的中心线的数据进行1维傅立叶逆变换后,从-∞到+∞为止,将频率积分值作为Imag1。接下来,由利用了反正切(arctangent)的如下公式,计算出模板激发1的相位角度Ph1。
Ph1=arctan(Imag1/Real1)    ......(1)
同样,进行从通过模板激发2采集的各MR信号(以下称为模板数据2)减去输送频率的余弦函数或正弦函数等处理,分别获得k空间数据的实数部分及虚数部分。接下来,对从模板激发2得到的k空间数据的实数部分的中心线的数据进行1维傅立叶逆变换后,从-∞到+∞为止,将频率积分值作为Real2。
同样,对从模板激发2得到的k空间数据的虚数部分的中心线的数据进行1维傅立叶逆变换后,从-∞到+∞为止,将频率积分值作为Imag2。然后,由如下公式(2),计算出模板激发2的相位角度Ph2。
Ph2=arctan(Imag2/Real2)    ......(2)
在以上述相位角度Ph1、Ph2为相位校正数据从通过主扫描采集的MR信号重建图像数据时,基于相位校正数据进行相位校正。
接下来,就磁场校正映射图的生成方法,以相位编码阶跃数以及频率编码数均为256的情况为例进行说明。不过,相位编码阶跃数以及频率编码阶跃数为256以外的数字也是可以的。
首先,作为第1相位图像,用从模板数据1得到的k空间数据的实数部分以及虚数部分生成纵横像素数为256×256的相位图像。具体而言,通过k空间数据的实数部分256×256的各矩阵元素和k空间数据虚数部分256×256的各矩阵元素,基于相同位置上的元素之比,计算出反正切。
例如,设k空间数据实数部分第1行第1列的矩阵元素的值为a,设k空间数据的虚数部分的第1行第1列的矩阵元素的值为b。
然后,计算出用θ=arctan(b/a)给出的角度θ,以此作为第1相位图像中第1行第1列的矩阵元素值。在以同样方法计算出所有256×256的矩阵元素值后,通过对所计算出的矩阵元素值的数据施以2维傅立叶变换,得到像素数为256×256的第1相位图像。
而且,模板数据2的相位图像2也用同样方法生成,以此作为第2相位图像。
接下来,(用处于相同位置的像素彼此)计算出第1相位图像的各像素值与第2相位图像的各像素值的差分后,乘以规定的系数生成差分图像(像素数256×256)。具体而言,例如,设第1相位图像的第1行第1列的像素值为θ1,设第2相位图像的第1行第1列的像素值为θ2。然后,设用以下的公式(3)给出的θsub为差分图像的第1行第1列的像素值。
θsub=γ×(θ1-θ2)/DT    ......(3)
公式(3)中,γ为氢原子的回转磁比率(42.6MHz),DT为在模板激发1、2中错开了的数据采集开始时间的时间差(秒)。
在图3的例子中,Gro翻转间隔的2倍成为DT。其它的像素也利用公式(3)同样地计算出像素值,藉此生成差分图像。将如此生成的第1相位图像与第2相位图像的差分图像作为磁场校正映射图。
说的更详细一点,在图3的模板激发1中,以使基于自旋回波法的有效回波时间中的MR信号数据被配置在k空间的中心线上的方式施加相位编码方向梯度磁场Gpe。基于自旋回波法的有效回波时间中的MR信号数据,由于施加了90°及180°RF脉冲,其选择薄片内的全氢原子核的自旋相位理想地相一致,只要将其配置在k空间的中心,即可生成基本上对称的第1相位图像作为磁场映射图。
另一方面,在模板激发2中,如前面所述,因自旋回波法的90°及180°RF脉冲的施加而使MR信号应当成为最大强度的时刻、与相位编码阶跃为零的时刻,错开一个Gro翻转间隔。
为此,在模板激发2中,作为与相位的错开量相应地非对称的磁场映射图,生成第2相位图像。即,在模板激发2中,作为与第1模板激发1的场合相比磁场不均匀的磁场映射图,生成第2相位图像。因此,通过如上所述生成第1相位图像与第2相位图像的差分图像,能够将其作为磁场校正映射图有效地利用。
图5是表示从模板数据1得到的第1相位图像之一例的示意图。图5中,为了简化,像素数粗取为20×20(下文的图6、图7亦同)。此外,图5中,像素值越为大的像素则越白(提高辉度等级),像素值越为小的像素则越黑(下文的图6、图7、图8、图9亦同)。如图5所示,从模板激发1能够得到基本上对称的作为磁场映射图的第1相位图像。
图6是表示从模板数据2得到的第2相位图像之一例的示意图。如图6所示,从模板激发2得到对称性比模板激发1的场合差的第2相位图像。
图7是表示用第1相位图像和第2相位图像基于公式(3)得到的差分图像(磁场校正映射图)之一例的示意图。如图7所示,通过获取分别从模板激发1、2得到的相位图像的差分,作为像素值的高低的梯度(勾配),相位差显现出来。
图8表示不在主扫描之前进行利用磁场校正映射图的静磁场不均匀性的校正而对均匀的模型进行摄像而得到的图像之一例。
图9表示在主扫描之前进行了利用磁场校正映射图的静磁场不均匀性的校正之后对与图8相同的模型进行摄像而得到的图像之一例。
图10是用于说明不进行基于磁场校正映射图的校正的情况(图8)与进行的情况(图9)之不同的示意图。图10抽取出存在于图8的图像中的4处畸变区域,将其作为畸变区域122、124、126、128。
不进行基于磁场校正映射图的校正的情况,如图8与图10的对比所示,出现了4处畸变区域122、124、126、128。进行基于磁场校正映射图的校正的场合,如图9所示,畸变区域122移动到了模型的左上侧外缘,畸变区域124移动到了模型的右上侧外缘,畸变区域126移动到了模型的右下侧外缘,畸变区域128移动到了模型的左下侧外缘。由此,畸变得到了大的改善,本来就应该均匀的模型被改善得比图8的场合均匀了。
图11是表示第1实施方式中的MRI装置20的动作流程的流程图。以下,适宜地参照前述各图,按照图11所示步骤序号,说明MRI装置20的动作。
(步骤S1)MPU86(参照图2)基于通过输入装置62对MRI装置20输入的摄像条件,进行MRI装置20的初始设定。在该初始设定中,设定定位图像摄像时的FR脉冲的临时中心频率等。
(步骤S2)序列控制器56驱动未图示的台板控制系,以使受检体P的摄像部位位于成为摄像空间的门架内的磁场中心的方式移动台板32(参照图1)。摄像部位是指将头部、胸部、腹部、腰部、腿部等受检体P的哪一部分作为摄像区域进行图像化的意思。
接下来,MRI装置20发射定位图像的数据采集用的RF脉冲等,MR信号由RF接收器48进行检测。RF接收器48对MR信号施以规定的信号处理,生成数字化的MR信号的原始数据,并将其输入到序列控制器56。
序列控制器56将MR信号的原始数据输入到图像重建部90,图像重建部90对该原始数据施以规定的处理后生成定位图像的图像数据,并将其输入到图像数据库94。
图像处理部96对被输入的图像数据施以规定的图像处理,存储装置66存储图像处理后的定位图像的图像数据。此后,显示控制部98按照MPU86的指令使定位图像显示在显示装置64上,基于此设定感兴趣的区域等摄像条件。
(步骤S3)MRI装置20按照以下的<1>~<4>的流程,执行用于生成相位校正数据以及磁场校正数据的2个种类的模板激发作为预扫描,之后,执行作为目的的图像用的数据采集(主扫描)。
再者,以下的模板激发1、2针对空间上彼此相同的区域,并且与主扫描的摄像区域相同的区域分别进行。
<1>MPU86,例如,把与基于到步骤S2为止被输入了的摄像条件所设定的主扫描的脉冲序列相同的脉冲序列设定为模板激发1的脉冲序列(参照图3)。接下来,MPU86将模板激发1的脉冲序列输入到序列控制器56,指示采集开始。
由此,MRI装置20按照模板激发1的脉冲序列,发射数据采集用的RF脉冲等,由RF接收器48检测MR信号。
序列控制器56将由RF接收器48检测并生成的MR信号的原始数据输入到图像重建部90。
图像重建部90在形成在k空间数据库92中的k空间里将该原始数据作为模板激发1的k空间数据进行配置。
<2>MPU86,把使相位编码阶跃脉冲以及读取方向梯度磁场Gro的施加开始时刻比模板激发1延迟了的序列设定为模板激发2的脉冲序列。此外,使相位编码阶跃为零的时刻的、模板激发1、2的各读取方向梯度磁场Gro的极性相反。关于模板激发1、2的区别的详细情况,由于已经用图3说明过了,所以略去说明。
接下来,MPU86将模板激发2的脉冲序列输入到序列控制器56,指示采集开始。由此,与上述<1>同样,在形成在k空间数据库92中的k空间里配置模板激发2的k空间数据。
<3>校正部100从图像重建部90的k空间数据库92分别获取模板激发1、2的k空间数据,生成第1相位图像,即,模板激发1的相位图像和第2相位图像,即,模板激发2的相位图像。
接下来,校正部100用第1相位图像和第2相位图像,利用公式(3)生成差分图像,将该差分图像作为磁场校正映射图输入到MPU86。第1及第2相位图像和差分图像的生成方法如上文所言相同。
<4>在基于磁场校正映射图进行了以使静磁场均匀化为目的的校正(匀场)之后,按照所设定的摄像条件进行主扫描的MR信号采集。
因此,MPU86计算用于抵消磁场校正映射图所表示的磁场的不均匀性使静磁场均匀化的补偿磁场。并且,MPU86以使补偿磁场重叠在静磁场上进行施加的方式通过序列控制器56控制各部。再者,一般情况下,静磁场的不均匀成分中的2次以上的成分能够通过,例如,向匀场线圈24供给电流的调节来匀场。此外,静磁场的不均匀成分中的1次成分能够通过,例如,各向X轴、Y轴、Z轴梯度磁场线圈26x、26y、26z供给电流的调节来匀场。
作为具体的动作,首先,利用被静磁场电源40励磁后的静磁场用磁铁22,在摄像空间里形成静磁场。然后,当从输入装置62向MPU86输入摄像开始指令时,MPU86把包括主扫描的脉冲序列等在内的摄像条件输入到序列控制器56。
序列控制器56按照所输入的脉冲序列驱动梯度磁场电源44、RF发射器46以及RF接收器48,藉此在摄像区域中形成薄片选择方向梯度磁场Gss、相位编码方向梯度磁场Gpe和读取方向梯度磁场Gro,并且从RF线圈28发生RF信号。
此时,在从梯度磁场电源44向X轴、Y轴、Z轴梯度磁场线圈26x、26y、26z的供给电流上叠加用于生成上述补偿磁场的量,相当于静磁场基本上被均匀化。
然后,由受检体P内的核磁共振所产生的MR信号被RF线圈28接收而被接收器48检测。RF接收器48对所检测到的MR信号进行上文所述的规定的信号处理,生成MR信号的原始数据,并将其输入到序列控制器56。
序列控制器56将原始数据输入到图像重建部90,图像重建部90将原始数据作为k空间数据配置在k空间数据库中形成的k空间里。
MRI装置20将以上<1>~<4>的处理作为1个薄片的数据采集,进行所有薄片的数据采集。
(步骤S4)校正部100针对所有薄片分别计算出相位校正数据。校正部100将各个相位校正数据与成为各相位校正数据的计算对象的各薄片分别进行关联后输入到图像重建部90。
具体而言,校正部100在计算出对模板数据1的k空间实数部分以及虚数部分的各中心线进行1维傅立叶逆变换后的各数据的各频率积分值之后,如上文所述,用反正切计算出模板数据1的相位角度Ph1。校正部100以同样的方法计算出模板数据2的相位角度Ph2,并将相位角度Ph1、Ph2作为与该模板数据1、2相对应的薄片的相位校正数据。
(步骤S5)图像重建部90从k空间数据库92读入k空间数据,对其施以包括傅立叶变换在内的图像重建处理,对各薄片的每个薄片重建图像数据。
在该图像重建处理中,图像重建部90用针对每一个薄片个别计算出的相位校正数据一并进行消除包含在MR信号中的相位误差影响的校正。再者,到生成相位校正数据为止的处理与以往的不同,但是生成相位校正数据之后的相位误差的校正处理可以与现有技术相同(参照专利文献2等),故省略详细的说明。
图像重建部90如上述那样生成所有薄片的图像数据,将这些图像数据保存在图像数据库94中。图像处理部96从图像数据库94读入图像数据,通过对其施以规定的图像处理生成2维的显示用图像数据,并将该显示用图像数据保存在存储装置66中。
(步骤S6)显示控制部98按照MPU86的控制从存储装置66获取显示用图像数据,并使显示用图像数据所表示的图像显示在显示装置64中。以上是第1实施方式的MRI装置20的动作说明。
这样,在第1实施方式中,在主扫描之前执行包括与各回波相对应的相位编码方向梯度磁场的施加并且回波的发生开始时刻彼此错开了的模板激发1、2,得到模板数据1、2。
并且,由于基于模板数据1、2生成磁场校正映射图和相位校正数据,所以能够一并进行主扫描中的使静磁场均匀化的校正和主扫描的MR信号采集后的相位误差校正。
此时,先基于磁场校正映射图进行使静磁场均匀化的校正,在对主扫描的MR信号进行了采集后,再进行基于相位校正数据的相位误差校正。即,主扫描的MR信号是从施加了基于磁场校正映射图进行了良好的均匀化之后的静磁场的摄像区域采集的。因此,能够提高相位误差的校正精度。
而且,由于基于从摄像薄片的同一断面采集的数据生成磁场校正数据而对各薄片中的每个薄片进行静磁场不均匀性的校正,所以静磁场的校正效果变高。
因此,比起以往,能够减轻起因于EPI中的相位误差的图像畸变(参照图8、图9)。
再者,对包含在各回波中的相位误差中起因于静磁场不均匀性的相位误差和其他原因引起的相位误差进行校正的方法多种多样,以往是分别单独进行的。
(第2实施方式)
第2及第3实施方式的MRI装置,作为装置构成,与第1实施方式的MRI装置20相同。第2实施方式中,不进行基于磁场校正映射图的使静磁场均匀化的校正,只进行相位误差的校正。以下,以与第1实施方式的不同之处为中心,对第2实施方式进行说明。
图12是表示第2实施方式中作为自旋回波系的单次激发EPI的主扫描、模板激发1、2’的脉冲序列之一例的时间图。图12中,横轴等的格式与图3相同,图12的上段的主扫描以及中段的模板激发1分别与第1实施方式的主扫描以及模板激发1相同。
图12中,只有下段的模板激发2’(TEMPLATE SHOT 2’)与第1实施方式的模板激发2不同。
例如,以相位编码梯度磁场Gpe的预脉冲的施加开始时为基准,在模板激发2’中,对应于相位编码阶跃脉冲为零的MR信号的产生时刻与模板激发1相同。这是因为,为了生成磁场校正映射图,在模板激发1、2中使对应于相位编码阶跃脉冲为零的MR信号的产生时刻错开,而第2实施方式中不生成磁场校正映射图。
在模板激发2’中,相位编码梯度磁场Gpe的预脉冲的面积(强度的时间积分值的绝对值)等于最初的4个相位编码阶跃脉冲的总计面积。即,在模板激发2’中,在时间上第5个被采集的MR信号与相位编码阶跃为零相对应,该MR信号的强度成为最大。
这里,注意与相位编码阶跃为零相对应的MR信号的产生时刻(图12中以纵虚线表示的有效回波时间)。该时刻的读取方向梯度磁场Gro的极性在模板激发1中为负,在模板激发2’中为正。这是因为,为了得到相位校正数据,最好使相位编码阶跃为零时刻的读取方向梯度磁场Gro的极性在模板激发1、2’中彼此相反。
通过使用以上的模板激发1、2’,能够以与第1实施方式同样的原理生成相位校正数据。
图13是表示第2实施方式中的MRI装置20的动作流程的流程图。以下,按照图13中所示步骤序号说明第2实施方式的MRI装置20的动作。
(步骤S21)与第1实施方式的步骤S1同样,进行MRI装置20的初始设定。
(步骤S22)与第1实施方式的步骤S2同样,生成定位图像的图像数据,将其在显示装置64中进行图像显示,设定感兴趣的区域等摄像条件。
(步骤S23)MRI装置20,按照以下的<1>~<3>的流程,执行用于生成相位校正数据的2个种类的模板激发作为预扫描,之后,执行主扫描。再者,模板激发1、2’针对空间上彼此相同的区域,并且与主扫描的摄像区域相同的区域分别进行。
<1>MRI装置20,与第1实施方式的步骤S3的<1>同样,执行模板激发1的脉冲序列。由此,将模板激发1的k空间数据配置(存储)在k空间数据库92中形成的k空间里。
<2>MPU86,在设定了上文所述的模板激发2’的脉冲序列后,将该脉冲序列输入到序列控制器56,指示采集开始。由此,在k空间数据库92中形成的k空间里配置模板激发2’的k空间数据。
<3>不进行基于磁场校正映射图的静磁场均匀化的校正,根据所设定的摄像条件,进行主扫描的MR信号的采集。关于MR信号采集的具体动作,除了不使用磁场校正映射图之外,与第1实施方式的步骤S3的<4>是同样的。
MRI装置20将以上<1>~<3>的处理作为一个薄片的数据采集进行所有薄片的数据采集。
(步骤S24)校正部100通过使用模板激发1、2’的k空间数据,与第1实施方式的步骤S4同样,针对所有薄片分别计算出相位校正数据。校正部100将计算出的各相位校正数据与计算对象的各薄片进行关联后输入到图像重建部90。
(步骤S25)与第1实施方式的步骤S5同样,图像重建部90一边基于各薄片的每一个的相位校正数据对包含在MR信号中的相位误差进行校正,一边重建各薄片的图像数据。此后,与第1实施方式的步骤S5同样,生成显示用图像数据,并将其保存在存储装置66中。
(步骤S26)与第1实施方式的步骤S6同样,进行显示用图像数据的图像显示。以上是对第2实施方式的MRI装置20的动作说明。
这样,在第2实施方式中,也能够以与第1实施方式同样的原理减少EPI中的相位误差。再者,在现有技术中,使相位误差减少的模板激发中没有施加相位编码方向梯度磁场Gpe。
(第3实施方式)
第3实施方式中,与第1实施方式同样,进行基于磁场校正映射图的使静磁场均匀化的校正和相位误差的校正。在第3实施方式中,进行3个模板激发。以下,以与第1实施方式的不同之处为中心,对第3实施方式进行说明。
图14是表示第3实施方式中作为自旋回波系的单次激发EPI的模板激发1、2’、3的脉冲序列之一例的时间图。在图14中,横轴等的格式与图3相同。
第3实施方式的主扫描的脉冲序列由于与第1实施方式的主扫描以及第3实施方式的模板激发1相同,所以没有进行图示。
图14的上段中表示的模板激发1(TEMPLATE SHOT 1)与第1实施方式的模板激发1相同。
图14的中段中表示的模板激发2’(TEMPLATE SHOT 2’)与第2实施方式的模板激发2’相同。在第3实施方式中,通过使用模板激发1、2’,与第2实施方式同样,生成相位校正数据,进行相位误差的校正。
图14的下段中表示的模板激发3(TEMPLATE SHOT 3)是用来生成磁场校正映射图的。在第3实施方式中,通过使用模板激发1、3,与第1实施方式同样,生成磁场校正映射图,进行静磁场均匀性的校正。为了生成磁场校正映射图,最好使模板激发1、3中与相位编码阶跃为零相对应的MR信号的产生时刻错开。
因此,在模板激发3中,与相位编码阶跃为零相对应的MR信号的产生时刻比模板激发1延迟Gro翻转间隔的1倍。这里的所谓“延迟”是指,例如,以相位编码方向梯度磁场Gpe的预脉冲的施加开始时为基准。
在模板激发3中,相位编码梯度磁场Gpe的预脉冲的面积(强度的时间积分值的绝对值)等于最初的5个相位编码阶跃脉冲的总计面积。即,在模板激发3中,按时间在第6个被采集的MR信号与相位编码阶跃为零相对应地被配置在k空间的中心线上,该MR信号的强度变为最大。
这里,注意与相位编码阶跃为零相对应的MR信号的产生时刻。该时刻在模板激发1中与图14的以纵向的虚线表示的有效回波时间(EFFECTIVE ECHO TIME)相对应,在模板激发3中,比模板激发1延迟一个回波的量。
这样做的话,与相位编码阶跃为零相对应的MR信号发生时刻的读取方向梯度磁场Gro的极性,在模板激发1、3双方中为负(参照图14)。即,在第1实施方式中,模板激发2由于用来生成相位校正数据以及磁场校正映射图的双方,所以使该时刻的读取方向梯度磁场Gro的极性相反。这是因为,其为生成相位校正数据所希望的。
然而,第3实施方式的模板激发3不用于相位校正数据的生成。在仅着眼于磁场校正映射图的生成的场合,无需使相位编码阶跃为零时刻的读取方向梯度磁场Gro的极性在2个模板激发中彼此相反。
图15是表示第3实施方式中的MRI装置20的动作流程的流程图。以下,按照图15所示的步骤序号说明第3实施方式的MRI装置20的动作。
(步骤S31)与第1实施方式的步骤S1同样,进行MRI装置20的初始设定。
(步骤S32)与第1实施方式的步骤S2同样,生成定位图像的图像数据,将其在显示装置64中进行图像显示,设定关心的区域等摄像条件。
(步骤S33)MRI装置20按照以下的<1>~<5>的流程,执行用于生成相位校正数据以及磁场校正映射图的3个种类的模板激发作为预扫描,之后,执行主扫描。再者,模板激发1、2’、3针对空间上彼此相同的区域,并且与主扫描的摄像区域相同的区域分别进行。
<1>MRI装置20与第1实施方式的步骤S3的<1>同样,执行模板激发1的脉冲序列。由此,将模板激发1的k空间数据配置(存储)在k空间数据库92中形成的k空间里。
<2>MRI装置20与第2实施方式的步骤S23的<2>同样,执行模板激发2’的脉冲序列。由此,将模板激发2’的k空间数据配置(存储)在k空间数据库92中形成的k空间里。
<3>MPU86在设定了上文所述的模板激发3的脉冲序列后,将该脉冲序列输入到序列控制器56,指示采集开始。由此,在k空间数据库92中形成的k空间里配置模板激发3的k空间数据。
<4>校正部100从图像重建部90的k空间数据库92分别获取模板激发1、3的k空间数据。校正部100生成第1相位图像,即,模板激发1的相位图像和第2相位图像,即,模板激发3的相位图像。
接下来,校正部100用第1相位图像和第2相位图像,与上文所述同样,利用公式(3)生成差分图像,将该差分图像作为磁场校正映射图输入到MPU86。
<5>与第1实施方式的步骤S3的<4>同样,在基于磁场校正映射图进行了以使静磁场均匀化为目的的校正之后,按照所设定的摄像条件进行主扫描的MR信号的采集。
MRI装置20将以上<1>~<5>的处理作为一个薄片的数据采集进行所有薄片的数据采集。
(步骤S34)校正部100通过使用模板激发1、2’的k空间数据,与第1实施方式的步骤S4同样,针对所有薄片分别计算出相位校正数据。校正部100将所计算出的相位校正数据与计算对象的各薄片进行关联后输入到图像重建部90。
(步骤S35)与第1实施方式的步骤S5同样,图像重建部90一边基于各薄片的每一个的相位校正数据对包含在MR信号中的相位误差进行校正,一边重建各薄片的图像数据。此后,与第1实施方式的步骤S5同样,生成显示用图像数据,并将其保存在存储装置66中。
(步骤S36)与第1实施方式的步骤S6同样,进行显示用图像数据的图像显示。以上是对第3实施方式的MRI装置20的动作说明。
这样,在第3实施方式中,也能够得到与第1实施方式同样的效果。
(实施方式的补充事项)
(1)如图3、图12、图14所示,在第1~第3实施方式中描述了在主扫描和模板激发1中使脉冲序列完全相同的例子。本发明的实施方式并非为这样的形态所限定。在主扫描和模板激发1中脉冲序列多少加些变更也是可以的。
但是,除了相位编码阶跃脉冲以及读取方向梯度磁场Gro的施加开始时间之外,最好使模板激发1、2、2’、3的脉冲序列相同。这是因为改变条件时仅抽取相位的差分可能会变得困难的缘故。关于以上这点,对以下的补充事项(2)、(3)也是同样的。
(2)像第1及第3实施方式那样,在进行了基于磁场校正映射图的静磁场均匀性的校正的基础上进行主扫描的MR信号采集的场合,最好各模板激发1、2、2’、3在主扫描之前进行。这是因为,在这种场合,如果没有在主扫描之前生成磁场校正映射图,则不能执行基于磁场校正映射图的静磁场均匀性的校正。在这种场合,各模板激发1、2、2’、3,先执行哪一方都是可以的。
但是,在像例如以下的变型例那样,在主扫描的MR信号采集时不执行基于磁场校正映射图的静磁场均匀性校正的场合,各模板激发(2、2’等)在主扫描之后进行也是可以的。
具体而言,在第1~第3实施方式中,不进行模板激发1,基于从通过主扫描采集的回波信号群(MR信号群)得到的数据,进行相位误差校正也是可以的。在这种场合,将,例如,主扫描的MR信号作为模板数据1,利用公式(1)计算出相位校正数据。由此,能够省略模板激发1,所以能够缩短摄像时间。
这里,在第1及第3实施方式中进行了基于磁场校正映射图用匀场线圈24等使静磁场均匀化的基础上执行主扫描的说明。但是,考虑以下那样的第1及第3实施方式的变型例。具体而言,在第1实施方式中省略模板激发1而使用通过主扫描得到的数据的变型例中,代替使用了匀场线圈24等的静磁场均匀化,只要进行,例如,以下那样的处理即可。即,将通过主扫描采集的MR信号群作为模板数据1,基于该模板数据1和模板数据2,在主扫描后的步骤S4中也进行磁场校正映射图的生成。然后,基于磁场校正映射图计算出畸变量,能够针对在对MR信号施以相位校正后重建的图像数据,进行畸变校正(像素的移动)。
在第3实施方式中省略模板激发1、3而使用由主扫描得到的数据的变型例,只要与上述第1实施方式的变型例同样进行即可。
(3)在第1~3实施方式中描述了自旋回波系的单次激发EPI的例子。上述各实施方式的原理也可以应用在场回波系等其他EPI中,并且不限于单次激发EPI,也可以应用在多次激发EPI中。
图16是以与图3同样的形式表示FID型(自由感应衰减信号:Free Induction Decay)的单次激发EPI中的主扫描、模板激发1”、2”的脉冲序列之一例的时间图。
图16中,上段、中段、下段分别表示主扫描、模板激发1”、2”的各脉冲序列。在图16中,将模板激发1’、2’中读取方向梯度磁场Gro的施加开始错开了回波间隔的2倍,但这不过是一个例子,错开2倍之外的自然数倍也是可以的。
(4)在上述各实施方式中,描述了将薄片设定为摄像区域,生成2维图像数据的例子。本发明的实施方式并非为这样的形态所限定。摄像区域为厚片也是可以的。此外,所生成的图像数据可以是2维图像用的图像数据,也可以是3维的体数据。
(5)作为MRI装置20,描述了在包含有静磁场用磁铁22、匀场线圈24、梯度磁场线圈单元26、RF线圈28的门架的外面存在RF接收器48的例子(参照图1)。本发明的实施方式并非为这样的形态所限定。RF接收器48包含在门架内的形态也是可以的。
具体而言,例如,将相当于RF接收器48的电子回路基板配设在内门架内。并且,由接收用RF线圈把从电磁波变换成模拟电信号的MR信号,通过该电子回路基板内的前置放大器增幅后,作为数字信号输出到门架外面,输入到序列控制器56也是可以的。在向门架外面输出时,由于只要使用,例如,光缆作为光数字信号发射,即可减轻外部噪音的影响,所以最为理想。
(6)对技术方案的用语和实施方式的对应关系进行说明。再者,以下所表示的对应关系不过是为了参考而表示的一种解释而已,并不是用来对本发明进行限定的。
静磁场用磁铁22、匀场线圈24、梯度磁场线圈26、RF线圈28、控制装置30整体(参照图1)在基于EPI序列的主扫描中采集MR信号的功能,是技术方案记载的摄像数据采集部的一个例子。
静磁场用磁铁22、匀场线圈24、梯度磁场线圈26、RF线圈28、控制装置30整体执行模板激发1后采集MR信号的功能,是技术方案记载的第1采集部的一个例子。但是,不执行模板激发1,而以主扫描的MR信号作为模板数据1利用公式(1)计算相位校正数据的场合,静磁场用磁铁22、匀场线圈24、梯度磁场线圈26、RF线圈28、控制装置30整体采集主扫描的MR信号的功能,是技术方案记载的第1采集部的一个例子。
静磁场用磁铁22、匀场线圈24、梯度磁场线圈26、RF线圈28、控制装置30整体执行模板激发2(或2’)后采集MR信号的功能,是技术方案记载的第2采集部的一个例子。
静磁场用磁铁22、匀场线圈24、梯度磁场线圈26、RF线圈28、控制装置30整体执行模板激发3后采集MR信号的功能,是技术方案记载的第3采集部的一个例子。
用于生成相位校正数据以及磁场校正映射图的校正部100的功能、基于相位校正数据进行相位误差的校正的图像重建部90的功能、基于磁场校正映射图在主扫描之前进行以使静磁场均匀化为目的的磁场校正的梯度磁场电源44、梯度磁场线圈26、控制装置30整体的功能,是技术方案记载的校正部的一个例子。
(6)虽然对本发明的几个实施方式进行了说明,这些实施方式是作为例子而提出的,并无限定发明的范围之意。这些实施方式可以以其他各种形态实施,在不脱离本发明的宗旨的范围内,可以进行各种省略、置换、变更。这些实施方式或其变形包含在发明的范围或宗旨内,同样,其包含在专利的技术方案书所记载的发明及其等同范围内。

Claims (20)

1.一种磁共振成像装置,所述磁共振成像装置能够执行EPI,即平面回波成像,所述EPI向静磁场内的受检体发射用于引发核磁共振的激励脉冲,对反复进行读取方向梯度磁场的极性翻转所产生的多个回波信号进行采集,基于多个所述回波信号重建所述受检体的图像数据,其特征在于,具有:
第1采集部,该第1采集部通过执行所述EPI的回波信号采集序列,从而将所产生的多个所述回波信号作为第1模板数据进行采集,该第1采集部所执行的所述EPI的回波信号采集序列包含相位编码方向梯度磁场的施加;
第2采集部,该第2采集部在所述第1模板数据的采集后通过执行所述EPI的回波信号采集序列,从而将所产生的多个所述回波信号作为第2模板数据进行采集,该第2采集部所执行的所述EPI的回波信号采集序列包含所述相位编码方向梯度磁场的施加,并且所述读取方向梯度磁场的施加开始时刻与所述第1模板数据的采集时相错开;以及
校正部,该校正部用所述第1及第2模板数据,至少进行包含在所述回波信号中的相位误差的校正。
2.如权利要求1所记载的磁共振成像装置,其特征在于:
所述校正部用所述第1及第2模板数据,进行所述相位误差的校正及用于使所述静磁场均匀化的磁场校正的双方。
3.如权利要求2所记载的磁共振成像装置,其特征在于:
还具有摄像数据采集部,该摄像数据采集部在所述第1及第2模板数据的采集后,通过执行所述EPI的回波信号序列,来采集来自所述受检体的多个所述回波信号;
所述校正部在由所述摄像数据采集部进行的多个所述回波信号的采集之前进行所述磁场校正,并且,在从由所述摄像数据采集部采集的多个所述回波信号重建所述图像数据时进行所述相位误差的校正。
4.如权利要求3所记载的磁共振成像装置,其特征在于:
所述第1及所述第2采集部针对与成为所述摄像数据采集部进行的多个所述回波信号的采集对象的摄像区域相同的区域,分别执行所述EPI的回波信号采集序列,籍此从所述相同的区域采集所述第1及第2模板数据。
5.如权利要求4所记载的磁共振成像装置,其特征在于:
所述第2采集部使所述读取方向梯度磁场的施加开始时刻与所述第1模板数据的采集时错开时间差,该时间差为从所述读取方向梯度磁场的极性翻转开始到再次翻转为止的时间差的自然数倍。
6.如权利要求4所记载的磁共振成像装置,其特征在于:
所述第2采集部,以使与相位编码阶跃为零相对应的所述回波信号发生时的所述读取方向梯度磁场的极性与所述第1采集部进行的所述第1模板数据的采集时相反的方式,进行所述第2模板数据的采集。
7.如权利要求6所记载的磁共振成像装置,其特征在于:
所述校正部,针对从所述第1及第2模板数据得到的各k空间数据,分别计算出基于对虚数数据的中心线进行了1维傅立叶逆变换后的数据的频率积分值与对实数数据的中心线进行了1维傅立叶逆变换后的数据的频率积分值之比的值的反正切,基于所述反正切进行所述相位误差的校正。
8.如权利要求7所记载的磁共振成像装置,其特征在于:
所述第2采集部,以使与相位编码阶跃为零相对应的所述回波信号的产生时刻与所述第1采集部进行的所述第1模板数据的采集时错开的方式,进行所述第2模板数据的采集;
所述校正部,生成基于所述第1模板数据的相位图像和基于所述第2模板数据的相位图像的差分图像作为所述磁场校正映射图,基于所述磁场校正映射图进行所述磁场校正。
9.如权利要求4所记载的磁共振成像装置,其特征在于:
所述第2采集部,以使与相位编码阶跃为零相对应的所述回波信号的产生时刻与所述第1采集部进行的所述第1模板数据的采集时错开的方式,进行所述第2模板数据的采集。
10.如权利要求9所记载的磁共振成像装置,其特征在于:
所述校正部生成与所述第1模板数据和所述第2模板数据的相位差相当的磁场校正映射图,基于所述磁场校正映射图进行所述磁场校正。
11.如权利要求10所记载的磁共振成像装置,其特征在于:
所述校正部生成基于所述第1模板数据的相位图像和基于所述第2模板数据的相位图像的差分图像作为所述磁场校正映射图。
12.如权利要求2所记载的磁共振成像装置,其特征在于:
所述校正部构成为,进行用于使所述静磁场均匀化的磁场校正,此后,进行所述相位误差的校正。
13.如权利要求1所记载的磁共振成像装置,其特征在于:
所述第2采集部,以使与相位编码阶跃为零相对应的所述回波信号发生时的所述读取方向梯度磁场的极性与所述第1采集部进行的所述第1模板数据的采集时相反的方式,进行所述第2模板数据的采集。
14.如权利要求2所记载的磁共振成像装置,其特征在于:
所述第2采集部,以使与相位编码阶跃为零相对应的所述回波信号的产生时刻与所述第1采集部进行的所述第1模板数据的采集时错开的方式,进行所述第2模板数据的采集。
15.如权利要求1所记载的磁共振成像装置,其特征在于:
还具有第3采集部,该第3采集部在所述第2模板数据的采集后通过执行所述EPI的回波信号采集序列,从而将所产生的多个所述回波信号作为第3模板数据进行采集,该第3采集部所执行的所述EPI的回波信号采集序列包含所述相位编码方向梯度磁场的施加,并且所述读取方向梯度磁场的施加开始时刻与所述第1模板数据的采集时相错开;
所述第2采集部,以使与相位编码阶跃为零相对应的所述回波信号发生时的所述读取方向梯度磁场的极性与所述第1采集部进行的所述第1模板数据的采集时相反的方式,进行所述第2模板数据的采集;
所述第3采集部,以使与相位编码阶跃为零相对应的所述回波信号的产生时刻与所述第1采集部进行的所述第1模板数据的采集时错开的方式,进行所述第3模板数据的采集;
所述校正部,用所述第1及第2模板数据进行所述相位误差的校正,用所述第1及第3模板数据进行用于使所述静磁场均匀化的磁场校正。
16.如权利要求15所记载的磁共振成像装置,其特征在于:
还具有摄像数据采集部,该摄像数据采集部在所述第1及第3模板数据的采集后,通过执行所述EPI的回波信号序列,来采集来自所述受检体的多个所述回波信号;
所述校正部在所述摄像数据采集部进行的多个所述回波信号的采集之前进行所述磁场校正,并且,在从由所述摄像数据采集部采集的多个所述回波信号重建所述图像数据时进行所述相位误差的校正。
17.如权利要求16所记载的磁共振成像装置,其特征在于:
所述第1、第2、第3采集部针对与成为所述摄像数据采集部进行的多个所述回波信号的采集对象的摄像区域相同的区域,分别执行所述EPI的回波信号采集序列,籍此从所述相同的区域采集所述第1、第2、第3模板数据;
所述校正部根据基于所述第1模板数据的相位图像和基于所述第3模板数据的相位图像的差分图像进行所述磁场校正,针对从所述第1及第2模板数据得到的各k空间数据,分别计算出基于对虚数数据的中心线进行了1维傅立叶逆变换后的数据的频率积分值与对实数数据的中心线进行了1维傅立叶逆变换后的数据的频率积分值之比的值的反正切,基于所述反正切进行所述相位误差的校正。
18.一种磁共振成像装置,所述磁共振成像装置能够执行EPI,即平面回波成像,所述EPI向静磁场内的受检体发射用于引发核磁共振的激励脉冲,对反复进行读取方向梯度磁场的极性翻转所产生的多个回波信号进行采集,基于多个所述回波信号重建所述受检体的图像数据,其特征在于,具有:
第1采集部,该第1采集部通过按照包含相位编码方向梯度磁场的施加的所述EPI的回波信号采集序列执行主扫描,来采集第1群的所述回波信号群;
第2采集部,该第2采集部通过执行所述EPI的回波信号采集序列,来采集第2群的所述回波信号群,该第2采集部所执行的所述EPI的回波信号采集序列包含所述相位编码方向梯度磁场的施加,并且所述读取方向梯度磁场的施加开始时刻与所述主扫描的所述EPI的回波信号采集时相错开;以及
校正部,该校正部用所述第1及第2群的回波信号群,至少进行包含在所述主扫描的回波信号中的相位误差的校正。
19.一种磁共振成像方法,所述磁共振成像方法包含EPI,即平面回波成像,所述EPI向静磁场内的受检体发射用于引发核磁共振的激励脉冲,对反复进行读取方向梯度磁场的极性翻转所产生的多个回波信号进行采集,基于多个所述回波信号重建所述受检体的图像数据,其特征在于,具有:
通过执行包含相位编码方向梯度磁场的施加的所述EPI的回波信号采集序列,从而将所产生的多个所述回波信号作为第1模板数据进行采集的步骤;
在所述第1模板数据的采集后,通过执行所述EPI的回波信号采集序列,来将所产生的多个所述回波信号作为第2模板数据进行采集的步骤,该步骤所执行的所述EPI的回波信号采集序列包含所述相位编码方向梯度磁场的施加,并且所述读取方向梯度磁场的施加开始时刻与所述第1模板数据的采集时相错开;以及
用所述第1及第2模板数据,至少进行包含在所述回波信号中的相位误差的校正的步骤。
20.一种磁共振成像方法,所述磁共振成像方法包含EPI,即平面回波成像,所述EPI向静磁场内的受检体发射用于引发核磁共振的激励脉冲,对反复进行读取方向梯度磁场的极性翻转所产生的多个回波信号进行采集,基于多个所述回波信号重建所述受检体的图像数据,其特征在于,具有:
通过按照包含相位编码方向梯度磁场的施加的所述EPI的回波信号采集序列执行主扫描,来采集第1群的所述回波信号群的步骤;
通过执行所述EPI的回波信号采集序列,来采集第2群的所述回波信号群的步骤,该步骤所执行的所述EPI的回波信号采集序列包含所述相位编码方向梯度磁场的施加,并且所述读取方向梯度磁场的施加开始时刻与所述主扫描的所述EPI的回波信号采集时相错开;以及
用所述第1及第2群的回波信号群,至少进行包含在所述主扫描的回波信号中的相位误差的校正的步骤。
CN201180001735.0A 2010-08-11 2011-08-02 磁共振成像装置及磁共振成像方法 Active CN102438518B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2010-180389 2010-08-11
JP2010180389 2010-08-11
PCT/JP2011/068014 WO2012020720A1 (ja) 2010-08-11 2011-08-02 磁気共鳴イメージング装置および磁気共鳴イメージング方法

Publications (2)

Publication Number Publication Date
CN102438518A true CN102438518A (zh) 2012-05-02
CN102438518B CN102438518B (zh) 2015-06-10

Family

ID=45567687

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180001735.0A Active CN102438518B (zh) 2010-08-11 2011-08-02 磁共振成像装置及磁共振成像方法

Country Status (5)

Country Link
US (1) US8890525B2 (zh)
EP (1) EP2457503B1 (zh)
JP (1) JP5835989B2 (zh)
CN (1) CN102438518B (zh)
WO (1) WO2012020720A1 (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102727203A (zh) * 2011-03-31 2012-10-17 株式会社东芝 磁共振成像装置
CN102768347A (zh) * 2012-08-09 2012-11-07 厦门大学 核磁共振波谱仪上缩减相位编码数快速三维梯度匀场方法
CN102890255A (zh) * 2012-10-19 2013-01-23 浙江大学 一种用于缩短回波时间的由内而外平面回波成像方法
CN103424726A (zh) * 2012-05-21 2013-12-04 西门子公司 连续校正磁共振测量序列的相位误差的方法和磁共振设备
CN104035059A (zh) * 2013-03-06 2014-09-10 上海联影医疗科技有限公司 平面回波成像序列图像的重建方法
WO2017113178A1 (zh) * 2015-12-30 2017-07-06 深圳先进技术研究院 一种聚焦超声位移成像方法及装置
CN108226831A (zh) * 2016-12-12 2018-06-29 西门子保健有限责任公司 加速的磁共振成像
CN109741409A (zh) * 2018-11-30 2019-05-10 厦门大学 回波平面成像涡流伪影的无参考扫描校正方法
CN109917315A (zh) * 2019-04-30 2019-06-21 上海联影医疗科技有限公司 磁共振成像扫描方法、装置、计算机设备和存储介质

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102011005728B4 (de) * 2011-03-17 2012-10-31 Siemens Aktiengesellschaft Verzeichnungsfreies Bestimmen von Magnetresonanzdaten
DE102011006230B4 (de) * 2011-03-28 2013-01-24 Siemens Aktiengesellschaft Pixelweise Korrektur einer Phaseninformation in MR-Bildern unter Verwendung eines Navigatorsignals
WO2013118117A1 (en) * 2012-02-08 2013-08-15 Anatech Advanced Nmr Algorithms Technologies Ltd Method and system for inspection of composite material components
JP6058477B2 (ja) * 2013-06-14 2017-01-11 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
JP6495057B2 (ja) * 2015-03-16 2019-04-03 キヤノンメディカルシステムズ株式会社 Mri装置及び撮像時間短縮方法
JP2016198392A (ja) 2015-04-13 2016-12-01 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
KR101819030B1 (ko) 2016-02-17 2018-01-16 삼성전자주식회사 자기 공명 영상 장치 및 자기 공명 영상 장치를 위한 방법
CA2964945A1 (en) * 2016-04-22 2017-10-22 Sunnybrook Research Institute System and method for producing distortion free magnetic resonance images using dual-echo echo-planar imaging
JP6283134B2 (ja) * 2017-03-07 2018-02-21 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置及び画像処理装置
EP3441781A1 (de) * 2017-08-11 2019-02-13 Siemens Healthcare GmbH Beschleunigte magnetresonanz-messung
US20210361167A1 (en) * 2018-09-17 2021-11-25 Tatiana WOLFE Signal isolation magnetic resonance image (simri) and methods thereof
WO2020238818A1 (en) * 2019-05-24 2020-12-03 Shanghai United Imaging Healthcare Co., Ltd. Imaging systems and methods
EP3796023A1 (de) * 2019-09-18 2021-03-24 Siemens Healthcare GmbH Verbessertes magnetresonanz-dixon-verfahren
EP3865892A1 (en) 2020-02-17 2021-08-18 Koninklijke Philips N.V. Iterative reconstruction of gradient echo magnetic resonance images
CN112731235B (zh) * 2020-12-09 2024-04-19 中国科学院深圳先进技术研究院 一种磁共振化学交换饱和转移成像方法以及相关设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09276243A (ja) * 1996-04-16 1997-10-28 Toshiba Corp 磁気共鳴診断装置
US5810727A (en) * 1996-04-24 1998-09-22 U.S. Philips Corporation Method of and device for the imaging of objects by means of magnetic resonance
JP2001238866A (ja) * 2000-03-01 2001-09-04 Hitachi Ltd 磁気共鳴イメージング装置
JP2005319074A (ja) * 2004-05-10 2005-11-17 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2006255046A (ja) * 2005-03-16 2006-09-28 Hamano Life Science Research Foundation 磁気共鳴映像法および画像処理装置
JP2007159718A (ja) * 2005-12-12 2007-06-28 Toshiba Corp 磁気共鳴イメージング装置および画像補正評価方法

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4885542A (en) * 1988-04-14 1989-12-05 The Regents Of The University Of California MRI compensated for spurious NMR frequency/phase shifts caused by spurious changes in magnetic fields during NMR data measurement processes
US4937526A (en) * 1988-11-23 1990-06-26 Mayo Foundation For Medical Education And Research Adaptive method for reducing motion and flow artifacts in NMR images
US5270654A (en) * 1991-07-05 1993-12-14 Feinberg David A Ultra-fast multi-section MRI using gradient and spin echo (grase) imaging
JP3152690B2 (ja) 1991-09-06 2001-04-03 株式会社東芝 磁気共鳴映像装置
JPH09253070A (ja) 1996-03-26 1997-09-30 Hitachi Ltd 核磁気共鳴を用いた検査装置
US5672969A (en) * 1996-04-25 1997-09-30 General Electric Company Reduction of Nyquist ghost artifacts in oblique echo planar imaging
US6249595B1 (en) * 1998-01-22 2001-06-19 General Electric Company Iterative reconstruction for EPI
JP3018076B2 (ja) * 1998-04-13 2000-03-13 技術研究組合医療福祉機器研究所 核磁気共鳴を用いた検査装置
US6160397A (en) * 1998-12-30 2000-12-12 General Electric Company Fast spin echo prescan for magnetic resonance imaging systems
US6728568B1 (en) * 1999-05-31 2004-04-27 Hitachi Medical Corporation Magnetic resonance imaging method and device
US6700374B1 (en) * 2000-03-29 2004-03-02 Koninklijke Philips Electronics, N.V. EPI calibration method to minimize ghosting in reconstructed images
JP4619674B2 (ja) * 2004-03-24 2011-01-26 株式会社東芝 磁気共鳴イメージング装置
US7098662B2 (en) * 2004-09-30 2006-08-29 General Electric Company Method and system of MR imaging with variable readout gradient filtering
US7102352B2 (en) * 2004-09-30 2006-09-05 General Electric Company Method and system for image artifact reduction using nearest-neighbor phase correction for echo planar imaging
US7259557B2 (en) * 2005-10-19 2007-08-21 General Electric Company Method and apparatus of echo planar imaging with real-time determination of phase correction coefficients
US7622924B2 (en) * 2007-06-12 2009-11-24 General Electric Company Method and apparatus for k-space and hybrid-space based image reconstruction for parallel imaging and artifact correction

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09276243A (ja) * 1996-04-16 1997-10-28 Toshiba Corp 磁気共鳴診断装置
US5810727A (en) * 1996-04-24 1998-09-22 U.S. Philips Corporation Method of and device for the imaging of objects by means of magnetic resonance
JP2001238866A (ja) * 2000-03-01 2001-09-04 Hitachi Ltd 磁気共鳴イメージング装置
JP2005319074A (ja) * 2004-05-10 2005-11-17 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2006255046A (ja) * 2005-03-16 2006-09-28 Hamano Life Science Research Foundation 磁気共鳴映像法および画像処理装置
JP2007159718A (ja) * 2005-12-12 2007-06-28 Toshiba Corp 磁気共鳴イメージング装置および画像補正評価方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NAN-KUEI CHEN, ET AL.: "Spectral-Domain fMRl based on Multiple Asymmetric Spin-Echo EPI", 《PROC. INTL. SOT. MAG. RESON. MED.》 *
XIAOPING HU,ET AL.: "Artifact Reduction in EPI with Phase-Encoded Reference Scan", 《MAGNETIC RESONANCE IN MEDICINE》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102727203A (zh) * 2011-03-31 2012-10-17 株式会社东芝 磁共振成像装置
CN103424726A (zh) * 2012-05-21 2013-12-04 西门子公司 连续校正磁共振测量序列的相位误差的方法和磁共振设备
CN103424726B (zh) * 2012-05-21 2016-06-22 西门子公司 连续校正磁共振测量序列的相位误差的方法和磁共振设备
CN102768347A (zh) * 2012-08-09 2012-11-07 厦门大学 核磁共振波谱仪上缩减相位编码数快速三维梯度匀场方法
CN102768347B (zh) * 2012-08-09 2014-07-09 厦门大学 核磁共振波谱仪上缩减相位编码数快速三维梯度匀场方法
CN102890255A (zh) * 2012-10-19 2013-01-23 浙江大学 一种用于缩短回波时间的由内而外平面回波成像方法
CN104035059B (zh) * 2013-03-06 2015-05-13 上海联影医疗科技有限公司 平面回波成像序列图像的重建方法
CN104035059A (zh) * 2013-03-06 2014-09-10 上海联影医疗科技有限公司 平面回波成像序列图像的重建方法
WO2017113178A1 (zh) * 2015-12-30 2017-07-06 深圳先进技术研究院 一种聚焦超声位移成像方法及装置
CN107205718A (zh) * 2015-12-30 2017-09-26 深圳先进技术研究院 一种聚焦超声位移成像方法及装置
CN108226831A (zh) * 2016-12-12 2018-06-29 西门子保健有限责任公司 加速的磁共振成像
CN108226831B (zh) * 2016-12-12 2021-03-19 西门子保健有限责任公司 磁共振成像方法和磁共振成像设备
CN109741409A (zh) * 2018-11-30 2019-05-10 厦门大学 回波平面成像涡流伪影的无参考扫描校正方法
CN109917315A (zh) * 2019-04-30 2019-06-21 上海联影医疗科技有限公司 磁共振成像扫描方法、装置、计算机设备和存储介质
CN109917315B (zh) * 2019-04-30 2021-09-28 上海联影医疗科技股份有限公司 磁共振成像扫描方法、装置、计算机设备和存储介质

Also Published As

Publication number Publication date
EP2457503B1 (en) 2019-07-17
JP2012055684A (ja) 2012-03-22
CN102438518B (zh) 2015-06-10
US20120146640A1 (en) 2012-06-14
US8890525B2 (en) 2014-11-18
JP5835989B2 (ja) 2015-12-24
WO2012020720A1 (ja) 2012-02-16
EP2457503A1 (en) 2012-05-30
EP2457503A4 (en) 2014-12-31

Similar Documents

Publication Publication Date Title
CN102438518B (zh) 磁共振成像装置及磁共振成像方法
Stockmann et al. In vivo B0 field shimming methods for MRI at 7 T
JP6058477B2 (ja) 磁気共鳴イメージング装置
US9335393B2 (en) MR parallel imaging system reducing imaging time
JP5063279B2 (ja) 磁気共鳴装置
Liao et al. High‐fidelity, high‐isotropic‐resolution diffusion imaging through gSlider acquisition with and T1 corrections and integrated ΔB0/Rx shim array
JP4619674B2 (ja) 磁気共鳴イメージング装置
WO2014129516A1 (ja) 磁気共鳴イメージング装置
CN102342832B (zh) 磁共振成像方法、磁共振成像装置及其控制装置
CN104224179A (zh) 一种磁共振成像系统的磁场稳定方法和装置
JP2002095646A (ja) 磁気共鳴イメージング装置
JP6285669B2 (ja) 磁気共鳴イメージング装置
CN103238082A (zh) 使用多点Dixon技术和低分辨率校准的MR成像
RU2015116879A (ru) Опорное сканирование при устойчивой к металлам мр визуализации
WO2013042466A1 (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP5666779B2 (ja) 磁気共鳴イメージング装置
CN107669272A (zh) 用于磁共振成像的连续波恒定振幅产生共振和偏共振自旋锁的系统和方法
JP5971683B2 (ja) 磁気共鳴イメージング装置
Cao et al. DTI‐MR fingerprinting for rapid high‐resolution whole‐brain T1, T2, proton density, ADC, and fractional anisotropy mapping
JP2001070283A (ja) Mrイメージングのデータ収集方法及びmri装置
JP2009172360A (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング装置の制御プログラム
Zhang et al. Validation of transmit SENSE with reciprocity
JP2002301041A (ja) 磁気共鳴イメージング装置
JP5395619B2 (ja) 磁気共鳴イメージング装置
JP2006087499A (ja) Mriデータの形成方法およびそれを用いたイメージング装置

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
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160713

Address after: Japan Tochigi

Patentee after: Toshiba Medical System Co., Ltd.

Address before: Tokyo, Japan, Japan

Patentee before: Toshiba Corp

Patentee before: Toshiba Medical System Co., Ltd.