CN102342832A - 磁共振成像方法、磁共振成像装置及其控制装置 - Google Patents

磁共振成像方法、磁共振成像装置及其控制装置 Download PDF

Info

Publication number
CN102342832A
CN102342832A CN2011102056065A CN201110205606A CN102342832A CN 102342832 A CN102342832 A CN 102342832A CN 2011102056065 A CN2011102056065 A CN 2011102056065A CN 201110205606 A CN201110205606 A CN 201110205606A CN 102342832 A CN102342832 A CN 102342832A
Authority
CN
China
Prior art keywords
mentioned
magnetic field
vortex flow
shooting area
leaning
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
CN2011102056065A
Other languages
English (en)
Other versions
CN102342832B (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 CN102342832A publication Critical patent/CN102342832A/zh
Application granted granted Critical
Publication of CN102342832B publication Critical patent/CN102342832B/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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56518Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
    • 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/483NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • G01R33/4836NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices using an RF pulse being spatially selective in more than one spatial dimension, e.g. a 2D pencil-beam excitation pulse
    • 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/5617Echo 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

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装置,它包括:信号收集部、图像生成部、位置取得部和修正部。信号收集部通过在放置受检体的拍摄空间内施加倾斜磁场和RF脉冲,收集从受检体产生的磁共振信号。图像生成部基于磁共振信号对受检体的图像数据进行重建。位置取得部取得拍摄区域作为拍摄空间内的位置信息。修正部通过以使得因伴随着倾斜磁场的施加产生的涡电流所产生的磁场即涡电流磁场相抵销的方式基于位置信息和涡电流磁场的时间常数使倾斜磁场的波形变形,而使产生的倾斜磁场分布接近目标分布。

Description

磁共振成像方法、磁共振成像装置及其控制装置
本申请基于巴黎公约将申请日为2010年7月23日的日本专利申请2010-166271和申请日为2011年6月29日的日本专利申请2011-144834作为优先权基础,享受它们的优先权。本申请通过援引这些申请而包含这些申请的全部内容。
技术领域
本发明的实施方式涉及磁共振成像方法、磁共振成像装置及其控制装置。
背景技术
MRI是用拉莫尔(Larmor)频率的RF脉冲磁激励放置在静磁场中的受检体的原子核自旋,根据伴随该激励产生的MR信号重建图像的拍摄法。另外,上述MRI是指磁共振成像(Magnetic ResonanceImaging),RF脉冲是指高频脉冲(Radio frequency pulse),MR信号是指核磁共振信号(Nuclear magnetic resonance signal)。
已知倾斜磁场分布的畸曲是MRI的图像质量劣化的主要原因之一。理想地,体层选择方向、相位编码方向、频率编码方向的各倾斜磁场的分布是例如磁场强度与沿施加方向的位置对应而线性地变化。但是实际情况是,一旦向倾斜磁场线圈供给脉冲电流,就会产生涡电流,在倾斜磁场上叠加因涡电流所导致的磁场,发生倾斜磁场分布的畸曲。
通常的涡电流补偿,由于以因涡电流所导致的磁场的各成分中的、在空间上呈线形变化的一次成分为对象,修正倾斜磁场波形的输入,所以不能对因涡电流所导致的磁场的二次以上的成分进行补偿。于是,有通过使电流流过高次的调整线圈而对因涡电流所导致的磁场的二次以上的成分进行补偿的技术。
另外,在日本特开平3-195539号公报中记载了,通过从预先设置的多个涡电流补偿电路中选择一个或多个来对因涡电流所导致的磁场进行补偿的构成。
发明内容
一般地,离磁场中心越远,倾斜磁场分布的畸曲越大。因此,在只对因涡电流所导致的磁场的一次成分进行补偿的现有技术中,尤其在以偏离中心区(off-center)为拍摄区域时,由于存在因涡电流所导致的磁场的二次以上的成分,伴随着倾斜磁场分布的畸曲产生图像质量劣化。另外,上述偏离中心区是指离开磁场中心的位置。
利用高次调整线圈对因涡电流所导致的磁场的二次以上的成分进行补偿的现有技术有为了避免与其它线圈耦合等的影响而设置例如几秒钟的等待时间的缺点。另外,有时高次调整成分与因涡电流所导致的磁场的高次成分的信道不对应。具体地说,即使调整线圈仅仅具有补偿XZ、YY等的各二次成分的各线圈,在因涡电流所导致的磁场的三次成分、四次成分大时,也不能充分地补偿因涡电流所导致的磁场。
本发明要解决的问题是,具有用来通过简单地补偿因涡电流所导致的磁场成分来改善图像质量的与以往不同的技术的MRI装置及其控制装置、MRI方法。
在一个实施方式中,MRI装置包括:信号收集部、图像生成部、位置取得部和修正部。信号收集部通过在放置受检体的拍摄空间内施加倾斜磁场和RF脉冲,收集从受检体产生的磁共振信号。图像生成部基于磁共振信号对受检体的图像数据进行重建。位置取得部取得拍摄区域,作为拍摄空间内的位置信息。修正部通过以使得因伴随着倾斜磁场的施加产生的涡电流所导致的磁场即涡电流磁场相抵销的方式,根据位置信息和涡电流磁场的时间常数使倾斜磁场的波形变形,从而使产生的倾斜磁场分布接近目标分布。如果采用上述构成的MRI装置,则可以通过简单地补偿因涡电流所导致的磁场成分来改善图像质量。
在一个实施方式中,控制装置是MRI装置的控制装置,该MRI装置通过向倾斜磁场线圈供给电流而在放置受检体的拍摄空间内施加倾斜磁场,伴随倾斜磁场执行受检体的磁共振成像,该控制装置包括上述位置取得部和上述修正部。如果采用上述构成的MRI装置的控制装置,则可以通过简单地补偿因涡电流所导致的磁场成分来改善图像质量。
在一个实施方式中,MRI方法包括:取得步骤、修正步骤、收集步骤和重建步骤。在取得步骤中,取得MRI的拍摄区域,作为在放置受检体的拍摄空间内的位置信息。在修正步骤中,通过以使得因伴随着倾斜磁场的施加产生的涡电流所导致的磁场即涡电流磁场相抵销的方式,基于上述位置信息和涡电流磁场的时间常数使倾斜磁场的波形变形,从而使产生的倾斜磁场的分布接近目标分布。在收集步骤中,通过在拍摄空间内施加通过修正步骤变形了的倾斜磁场和RF脉冲,收集从受检体产生的磁共振信号。在重建步骤中,基于磁共振信号对受检体的图像数据进行重建。如果采用上述构成的MRI方法,则可以通过简单地补偿因涡电流所导致的磁场成分来改善图像质量。
附图说明
图1是示出根据本实施方式的MRI装置的整体构成的框图。
图2是图1所示的计算机58的功能框图。
图3是示出在Y轴正方向上离磁场中心5cm的区域中的因涡电流所导致的磁场的一成分的X轴方向强度分布的一例的示意图。
图4是示出X偏离中心区处的各涡磁场成分的各参数的值的一例的表。
图5是示出Y偏离中心区处的各涡磁场成分的各参数的值的一例的表。
图6是示出Z偏离中心区处的各涡磁场成分的各参数的值的一例的表。
图7是对比受到涡电流影响的倾斜磁场波形的一例和与其对应的理想倾斜磁场波形的一例的示意图。
图8是示出在Y轴正方向上离磁场中心5cm的区域中的、因涡电流所导致的磁场的另一成分的X轴方向强度分布的一例的示意图。
图9是示出多体层拍摄中的拍摄体层的一例的示意图。
图10是示出以自旋回波法为例的多体层拍摄中的脉冲序列的一例的示意图。
图11是示出根据本实施方式的MRI装置的动作流程的流程图。
图12是示出不进行考虑了涡磁场成分的倾斜磁场的修正而拍摄体模得到的图像的一例的图。
图13是示出根据本实施方式的方法进行倾斜磁场分布的修正,拍摄与图12相同的体模得到的图像的一例的图。
具体实施方式
下面,基于附图说明MRI装置及其控制装置、MRI方法的实施方式。另外,在各图中对相同要素赋予相同的附图标记,重复说明省略。
(本实施方式的构成)
图1是示出根据本实施方式的MRI装置20的整体构成的框图。像图1所示的那样,MRI装置20包括:形成静磁场的筒状的静磁场磁体22;以及在该静磁场磁体22的内侧同轴地设置的筒状的调整线圈24、倾斜磁场线圈26、RF线圈28、控制装置30以及供受检体H躺在其上的床板32。
在此,作为一例,像以下那样定义装置坐标系的相互正交的X轴、Y轴、Z轴。首先,将静磁场磁体22和调整线圈24配置成使它们的轴方向与铅直方向正交,以静磁场磁体22和调整线圈24的轴方向为Z轴方向。另外,以铅直方向为Y轴方向,将床板32配置成其顶板的用来载置的面的法线方向为Y轴方向。
MRI装置20的控制装置30具有:例如,静磁场电源40、调整线圈电源42、倾斜磁场电源44、RF发送器46、RF接收器48、床板驱动装置50、序列控制器56和计算机58。
倾斜磁场电源44由X轴倾斜磁场电源44x、Y轴倾斜磁场电源44y和Z轴倾斜磁场电源44z构成。另外,计算机58由运算装置60、输入装置62、显示装置64和存储装置66构成。
静磁场磁体22与静磁场电源40连接,利用从静磁场电源40供给的电流在拍摄空间中形成静磁场。
上述拍摄空间是指,例如,放置受检体H,被施加静磁场的机架内的空间。所谓机架,是以包含静磁场磁体22、调整线圈24、倾斜磁场线圈26、RF线圈28的方式形成为例如圆筒状的结构体。机架和床板32构成为,受检体H所躺的床板32可以在机架的内部移动。另外,在图1中,由于太复杂,所以作为构成要素只示出机架内的静磁场磁体22、调整线圈24、倾斜磁场线圈26、RF线圈28,未示出机架自身。
拍摄区域是指,例如,用于生成“一幅图像”或“一组图像”的MR信号的收集范围,是作为拍摄空间的一部分设定的区域。此处的“一幅图像”和“一组图像”有时是二维图像,也有时是三维图像。此处的“一组图像”是例如像多体层拍摄等那样,在一个脉冲序列内一并收集“多幅图像”的MR信号时的“多幅图像”。拍摄区域例如用装置坐标系三维地规定。在此,作为一例,拍摄区域如果是厚度薄的区域则称为拍摄体层,如果是有一定程度的厚度的区域则称为拍摄板区。
调整线圈24与调整线圈电源42连接,利用从调整线圈电源42供给的电流,将其静磁场均匀化。静磁场磁体22多数情况下用超导线圈构成,在激励时与静磁场电源40连接而被供给电流,但一般情况下一旦被激励后就变成非连接状态。另外,也可以不设置静磁场电源40,用永磁体构成静磁场磁体22。
倾斜磁场线圈26包括X轴倾斜磁场线圈26x、Y轴倾斜磁场线圈26y和Z轴倾斜磁场线圈26z,在静磁场磁体22的内侧筒状地形成。X轴倾斜磁场线圈26x、Y轴倾斜磁场线圈26y和Z轴倾斜磁场线圈26z分别与倾斜磁场电源44的X轴倾斜磁场电源44x、Y轴倾斜磁场电源44y和Z轴倾斜磁场电源44z连接。
然后,利用分别从X轴倾斜磁场电源44x、Y轴倾斜磁场电源44y和Z轴倾斜磁场电源44z向X轴倾斜磁场线圈26x、Y轴倾斜磁场线圈26y和Z轴倾斜磁场线圈26z供给的电流,在拍摄空间中分别形成X轴方向的倾斜磁场Gx、Y轴方向的倾斜磁场Gy和Z轴方向的倾斜磁场Gz。
即,可以将装置坐标系的三轴方向的各倾斜磁场Gx、Gy、Gz合成,任意地设定作为逻辑轴的体层选择方向倾斜磁场Gss、相位编码方向倾斜磁场Gpe和读出方向(频率编码方向)倾斜磁场Gro的各方向。体层选择方向、相位编码方向和读出方向的各倾斜磁场Gss、Gpe、Gro与静磁场重叠。
RF发送器46基于从序列控制器56输入的控制信息,生成用来产生核磁共振的拉莫尔频率的RF脉冲(RF电流脉冲),将它发送到发送用的RF线圈28。在RF线圈28中具有:内置在机架中的接收发送RF脉冲用的全身用线圈、设置在床板32或受检体H附近的接收RF脉冲用的局部线圈等。发送用的RF线圈28从RF发送器46接收RF脉冲并向受检体H发送。接收用的RF线圈28接收通过用RF脉冲激励受检体H内部的原子核自旋而产生的MR信号,并用RF接收器48检测该MR信号。
RF接收器48通过对检测到的MR信号实施前置放大、中间频率变换、相位检波、低频放大、滤波等的各种信号处理,然后进行A/D变换,生成数字化了的复数数据即原始数据。RF接收器48将生成的MR信号的原始数据输入给序列控制器56。
运算装置60进行MRI装置20整体的系统控制。
序列控制器56根据运算装置60的指令,存储为了驱动倾斜磁场电源44、RF发送器46和RF接收器48所需的控制信息。此处的“控制信息”是指,例如,记述了应向倾斜磁场电源44施加的脉冲电流的强度、施加时间、施加定时等的动作控制信息的序列信息。
序列控制器56通过根据所存储的预定的序列驱动倾斜磁场电源44、RF发送器46和RF接收器48,产生倾斜磁场Gx、Gy、Gz和RF脉冲。另外,序列控制器56接收从RF接收器48输入的MR信号的原始数据,将它输入运算装置60。
床板驱动装置50经由序列控制器56与运算装置60连接。序列控制器56通过根据运算装置60的指令控制床板驱动装置50,使床板32的顶板移动。
图2是图1所示的计算机58的功能框图。计算机58的运算装置60包括:MPU(微处理单元)86、系统总线88、图像重建部90、图像数据库94、图像处理部96、显示控制部98、位置取得部100和修正部102。
MPU 86在拍摄条件的设定、拍摄动作和拍摄后的图像显示中,通过系统总线88等的布线进行MRI装置20整体的系统控制。另外,MPU 86作为拍摄条件设定部起作用,基于来自输入装置62的指示信息设定包含脉冲序列的拍摄条件,将设定的拍摄条件输入序列控制器56。因此,MPU86控制显示控制部98,在显示装置64上显示拍摄条件的设定用画面信息。
输入装置62向用户提供设定拍摄条件、图像处理条件的功能。
图像重建部90在内部具有k空间数据库92。图像重建部90在k空间数据库92中形成的k空间中配置从序列控制器56输入的MR信号的原始数据作为k空间数据。图像重建部90对k空间数据实施包含二维傅里叶变换等的图像重建处理,生成受检体H的各体层的图像数据。图像重建部90将生成的图像数据存储到图像数据库94中。
图像处理部96从图像数据库94取入图像数据,对它实施预定的图像处理,将图像处理后的图像数据作为显示用图像数据存储到存储装置66中。
存储装置66对上述显示用图像数据附加该显示用图像数据的生成中使用的拍摄条件、受检体H的信息(患者信息)等,作为附带信息,并存储起来。
显示控制部98根据MPU 86的控制,在显示装置64上显示拍摄条件的设定用画面、表示通过拍摄生成的图像数据的图像。
位置取得部100取得拍摄区域,作为在拍摄空间内的基于装置坐标系的位置信息,将它输入到修正部102。
修正部102算出通过向倾斜磁场线圈26供给脉冲电流而生成的涡电流所导致的磁场的大小。
下面,说明因涡电流所导致的磁场的大小的计算方法。首先,在此作为一例,使磁场中心位于装置坐标系的原点。另外,“X偏离中心区”的含义是沿X轴方向从磁场中心偏离了的位置。同样地,“Y偏离中心区”、“Z偏离中心区”的含义分别是沿Y轴方向、Z轴方向从磁场中心偏离了的位置。
图3是示出在Y轴正方向上离磁场中心5cm的区域中的因涡电流所导致的磁场的一个成分的X轴方向强度分布的一例。图3的横轴表示偏离中心区的轴,即X轴方向的位置。图3的纵轴是观测轴,表示因涡电流所导致的磁场的Y成分的二次成分中的某时间常数τ的磁场强度。在Y轴负方向上离磁场中心5cm的区域中的同一磁场成分的磁场强度的X轴方向分布,是将图3的纵轴方向正负逆转得到的。
像上述的一例那样,因涡电流产生的磁场强度,不管哪个成分一般情况下都是离磁场中心越远越大,因涡电流所导致的磁场的影响度因拍摄区域的位置不同而不同。于是,修正部102基于因涡电流所导致的磁场的大小和拍摄区域的位置,算出向X轴、Y轴、Z轴的各倾斜磁场线圈26x、26y、26z供给的电流的修正值,以使各倾斜磁场Gss、Gpe、Gro的分布成为目标分布。在此,因涡电流所产生的磁场的各成分像以下那样地定义。
在X偏离中心区中,施加X轴方向的倾斜磁场Gx时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是XX_X、XY_X、XZ_X。
在X偏离中心区中,施加Y轴方向的倾斜磁场Gy时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是YX_X、YY_X、YZ_X。
在X偏离中心区中,施加Z轴方向的倾斜磁场Gz时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是ZX_X、ZY_X、ZZ_X。
在Y偏离中心区中,施加X轴方向的倾斜磁场Gx时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是XX_Y、XY_Y、XZ_Y。
在Y偏离中心区中,施加Y轴方向的倾斜磁场Gy时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是YX_Y、YY_Y、YZ_Y。
在Y偏离中心区中,施加Z轴方向的倾斜磁场Gz时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是ZX_Y、ZY_Y、ZZ_Y。
在Z偏离中心区中,施加X轴方向的倾斜磁场Gx时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是XX_Z、XY_Z、XZ_Z。
在Z偏离中心区中,施加Y轴方向的倾斜磁场Gy时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是YX_Z、YY_Z、YZ_Z。
在Z偏离中心区中,施加Z轴方向的倾斜磁场Gz时因涡电流所产生的磁场的X成分、Y成分、Z成分分别是ZX_Z、ZY_Z、ZZ_Z。
即,由包含下横线的4个字符表示的各符号的最前面的字符表示对应于施加倾斜磁场Gx、Gy、Gz中的哪一个时的情形。另外,各符号的第二个字符表示着眼于(作为因涡电流所产生的磁场的观测轴)X成分、Y成分、Z成分中的哪一个。另外,各符号的末尾的字符表示对应于X偏离中心区、Y偏离中心区、Z偏离中心区中的哪一个。
下面,将因涡电流所产生的磁场的各成分称为涡磁场成分,像涡磁场成分XX_X那样地来标示。各涡磁场成分,虽然是振幅和按照时间常数τ呈指数函数关系衰减的衰减项的积,但有时也具有不同的时间常数τ的成分。于是,对于上述27个涡磁场成分,针对每个时间常数分别算出振幅。
另外,涡电流磁场的时间常数τ由涡电流的时间常数确定,例如由以下因素确定。第一是X轴、Y轴、Z轴倾斜磁场线圈26x、26y、26z与静磁场磁体22的相对位置关系。第二是倾斜磁场线圈26的形状、材质等。例如,通过在安装调整MRI装置20时测定涡磁场成分而算出这些时间常数τ,将算出数据做成像表格数据那样,存储到修正部102中即可。或者,也可以基于倾斜磁场线圈26与静磁场磁体22的相对位置关系、倾斜磁场线圈26的形状、材质等,通过模拟算出涡电流磁场的时间常数τ,将算出结果存储到修正部102中。
下面,说明对于上述27个涡磁场成分针对每个时间常数算出振幅的方法。
第一,对于上述27个涡磁场成分内的偏离中心区轴与观测轴相同的涡磁场成分(符号的末尾的字符与第二个字符相同),作为一例用下式计算振幅。具体地说,对XX_X、YX_X、ZX_X、XY_Y、YY_Y、ZY_Y、XZ_Z、YZ_Z、ZZ_Z这9个涡磁场成分用下式(1)计算振幅L。
L=A1+A2×{pol(p)}×|p|
+A3×|p|2+A4×{pol(p)}×|p|3   …(1)
在式(1)中,|p|表示p的绝对值,pol(p)表示正或负的p的极性,根据p的极性代入+1或-1。
第二,对于上述27个涡磁场成分内的偏离中心区轴与观测轴不同的涡磁场成分,作为一例用下式(2)计算振幅。具体地说,对XY_X、XZ_X、YY_X、YZ_X、ZY_X、ZZ_X、XX_Y、XZ_Y、YX_Y、YZ_Y、ZX_Y、ZZ_Y、XX_Z、XY_Z、YX_Z、YY_Z、ZX_Z、ZY_Z这18个涡磁场成分用下式(2)计算振幅L。
L=A1×p+A2×p2+A3×p3+A4×p4  …(2)
式(1)和式(2)中的A1、A2、A3、A4针对每个涡磁场成分不同。关于A1~A4,通过例如在安装调整MRI装置20时测定各涡磁场成分而取得,将取得的数据做成像查询表(表格数据)那样,输入设定并存储到修正部102中即可。
图4是示出X偏离中心区处的各涡磁场成分的各参数的值的一例的表。
图5是示出Y偏离中心区处的各涡磁场成分的各参数的值的一例的表。
图6是示出Z偏离中心区处的各涡磁场成分的各参数的值的一例的表。
另外,图4~图6中的时间常数τ的单位是毫秒。像上述那样,27种涡磁场成分存在根据时间常数进一步分成多项的情况。在图4的例子中,涡磁场成分XX_X、YY_X、ZZ_X分成多项。另外,在图4~图6中,A1、A3、A4一律为零,但这不过是为了使以下的说明简化而采用的一个例子而已,在实际中基于测定值等使用适当的值。
在本实施方式中,通过使用像式(1)、式(2)、图4~图6那样预先存储的多项式模型化了的数据,进行针对X偏离中心区、Y偏离中心区、Z偏离中心区处的因涡电流所产生的磁场的一次和二次以上的成分的修正。由此,使倾斜磁场的分布接近目标分布。
具体地说,按照装置坐标系将拍摄空间内的(作为拍摄体层或拍摄板区的)拍摄区域的中心的坐标位置规定为(x1,y1,z1)。此时,用X坐标值的x1计算X偏离中心区成分(上述27个涡磁场成分内符号末尾字符为X的成分)。同样地,用Y坐标值的y1计算Y偏离中心区成分,用Z坐标值的z1计算Z偏离中心区成分。
由此,针对全部27个涡磁场成分的各时间常数的项计算其振幅。作为一例,具体地说明将基于装置坐标系的拍摄体层的中心坐标设为(-10,5,2),计算各涡磁场成分的振幅的例子。
对于涡磁场成分XX_X,由于(因符号的第二个字符与末尾的字符相同)观测轴与偏离中心区轴相同,所以用上述式(1)根据图4的各参数值计算。另外,由于涡磁场成分XX_X(因为符号的末尾字符)是X偏离中心区成分,所以式(1)中的p值用X坐标值-10。这样一来,像下式(3)那样计算出涡磁场成分XX_X的第1项(对应于时间常数τ=40毫秒)的振幅LXX_X1
LXX_X1=A1+A2×{pol(p)}×|p|
+A3×|p|2+A4×{pol(p)}×|p|3
=0+30×(-1)×|-10|+0×|-10|2+0×(-1)×|-10|3
=-300…(3)
即,如果倾斜磁场Gx的施加开始时刻为时间t=0,在X坐标值为-10的X偏离中心区中,施加X轴方向的倾斜磁场Gx时因涡电流所产生的磁场的X成分的第一项,则大致用下式(4)表示。
LXX_X1×{exp(-t/τ)}
=-300×{exp(-t/0.04)}…(4)
另外,由式(4)给出的磁场强度的单位为例如特斯拉。同样地,像下式(5)那样计算出涡磁场成分XX_X的第2项(对应于时间常数τ=3毫秒)的振幅LXX_X2
LXX_X2=A1+A2×{pol(p)}×|p|
+A3×|p|2+A4×{pol(p)}×|p|3
=-32…(5)
对于涡磁场成分YX_Z,由于偏离中心区轴与观测轴不同,所以用上述式(2)根据图6的各参数值计算。另外,由于涡磁场成分YX_Z是Z偏离中心区的成分,所以用Z坐标值2作为式(2)中的p值。在图6的例子中,涡磁场成分YX_Z只有对应于时间常数τ=40毫秒的项,所以像下式(6)那样计算出其振幅LYX_Z1
LYX_Z=A1×p+A2×p2+A3×p3+A4×p4
=0×2+1.5×22+0×23+0×24=6  …(6)
这样,通过使用像式(1)、式(2)、图4~图6那样的表格数据,可以对于27个涡磁场成分,针对每个时间常数全部计算其振幅。由此,可以算出因涡电流所产生的磁场分布的随时间的变化。
在修正倾斜磁场分布时,通过修正向倾斜磁场线圈26供给的电流以使全部涡磁场成分相抵销,使倾斜磁场的时间波形变形即可。例如,在不考虑涡电流时,算出向施加作为目标的倾斜磁场分布(目标分布)的X轴、Y轴、Z轴的各倾斜磁场线圈26x、26y、26z供给的电流值,作为“暂定值”。
然后,算出向施加与涡磁场成分符号相反的磁场分布的各倾斜磁场线圈26x、26y、26z供给的电流值,作为“调整值”。然后,通过将“暂定值”和“调整值”相加来确定向各倾斜磁场线圈26x、26y、26z最终供给的电流值即可。
或者,也可以修正向各倾斜磁场线圈26x、26y、26z供给的电流值,以施加从作为目标的倾斜磁场分布减去全部涡磁场成分得到的倾斜磁场分布。
在此,涡磁场成分像式(4)那样按照时间常数衰减。因此,作为考虑了该时间常数的修正,使针对向倾斜磁场线圈26供给的电流值的“调整值”与从各倾斜磁场脉冲的施加开始时刻(与后述的图7的t1相当)算起经过的时间对应地变化。
另外,所谓“各倾斜磁场脉冲”在后述的图10所示的多体层拍摄的脉冲序列的情况下,是指例如体层选择脉冲200、202、相位编码脉冲204、读出脉冲206、208的各个。即,“各倾斜磁场脉冲”一般是在体层选择方向、相位编码方向、读出方向的各逻辑轴上,在时间轴上离散的脉冲。
例如,考虑上述的通过将“暂定值”和“调整值”相加来确定向各倾斜磁场线圈26x、26y、26z最终供给的电流值的修正。此时,由于“调整值”是向施加与涡磁场成分符号相反的磁场分布的倾斜磁场线圈26供给的电流值,所以随着时间的经过,按照上述时间常数减少。
这是因为,如果涡磁场成分按照时间常数衰减,向倾斜磁场线圈26供给的电流值内的作为修正成分的“调整值”的绝对值也应当随着时间经过而减少。即,是随着从各倾斜磁场脉冲的开始施加算起的时间经过而与修正量相当的“调整值”减少的那样的修正。
在像以上那样使倾斜磁场波形变形的修正中,如果不考虑涡电流,则得到从作为目标的倾斜磁场分布(目标分布)减去全部涡磁场成分得到的倾斜磁场分布,但在实际中涡磁场成分相重叠,先减去的全部涡磁场成分相抵销。即,得到大致为目标的倾斜磁场分布。
图7是对比受到涡电流影响的倾斜磁场波形的一例和与其对应的理想倾斜磁场波形的一例的示意图。图7(A)中,用粗线表示受到涡电流影响的倾斜磁场波形的一例,用虚线表示与该波形对应的理想倾斜磁场波形。图7(B)用粗线表示理想倾斜磁场波形,与图7(A)的虚线部分相同。另外,图7(A)、(B)中,纵轴表示倾斜磁场的强度,横轴表示经过的时间t。
图7中,时刻t1为倾斜磁场脉冲的施加开始时刻(与向倾斜磁场线圈26供给的电流供给开始时刻大致相同)。在不进行针对涡电流磁场的补偿时,像图7(A)所示的那样,倾斜磁场波形称为曲线状(钝化)。这是因为,与向倾斜磁场线圈26开始供给电流同步,还产生涡电流磁场,涡电流磁场重叠在倾斜磁场上。但是,涡电流磁场随着时间经过,按照上述时间常数τ减少。因此,在图7中的梯形的理想倾斜磁场波形的平坦部分的后半部分几乎没有涡电流磁场,所以理想倾斜磁场波形与受到涡电流影响的倾斜磁场波形基本上一致。
但是,在本实施方式中,由于使倾斜磁场的波形变形以使涡电流磁场相抵销,所以得到像图7(B)那样的理想倾斜磁场波形。这是因为,在本实施方式中,像上述那样,将向倾斜磁场线圈26供给的电流值修正成基于拍摄区域的位置和涡电流磁场的时间常数τ随时间变化。
图3中的拍摄区域Img1是中心坐标为(像上述的(-10、5、2)的例子那样)在Y轴负方向上离磁场中心10cm的位置的一例。例如,在拍摄肩、手腕等时,像拍摄区域Img1那样位于X轴方向的一侧。
此时,通过修正倾斜磁场Gx、Gy、Gz以使得以拍摄区域Img1的中心位置为基准的涡磁场成分相互抵销,可以在拍摄区域Img1的中心处基本上完全除去涡磁场成分的影响。在拍摄区域Img1的外缘区域,由于各涡磁场成分的振幅与拍摄区域Img1的中心不同,虽然不能完全地修正,但与只对一次成分进行修正的现有技术相比可以改善图像质量。
图8是与图3同样地,示出在Y轴正方向上离磁场中心5cm的区域中的、因涡电流所产生的磁场的另一成分的X轴方向强度分布的一例的示意图。图8的横轴表示X轴方向的位置。图8的纵轴是因涡电流所产生的磁场的Y成分的二次成分中的某时间常数τ的磁场强度。
图8中用斜线区域表示的拍摄区域Img2,作为一例,是在X轴方向涵盖-10cm~10cm的范围、Y轴方向涵盖-10cm~10cm的范围、Z轴方向涵盖-0.5cm~0.5cm的范围的长方体状的拍摄板区。即,拍摄区域Img2的中心与装置坐标系的原点(在该例中是磁场中心)一致。
在此,对于上述涡磁场成分的振幅,也可以用拍摄区域内的平均值求出。以拍摄区域Img2为例,说明针对X轴方向、Y轴方向,使涡磁场成分的偏离平均化的一例。
例如,涡磁场成分XY_X,在图4的例子中只有时间常数50毫秒的一项,通过将拍摄区域Img2的X坐标值代入而作为式(2)中的p,可以计算其振幅LXY_X。在此,作为X坐标值,像-10,-9,……0,1,……9,10那样以1cm刻度抽出21个计测点,以涵盖拍摄区域Img2的X轴方向的一端到另一端。然后在各计测点计算LXY_X,计算用计测点数将它平均化而得到的平均振幅AVLXY_X。在此,由于在图4的例子中A1=0、A3=0、A4=0,所以可以像下式(7)那样计算平均振幅AVLXY_X
AVL XY _ X = 1 21 × Σ p = - 10 10 { A 1 × p + A 2 × p 2 + A 3 × p 3 + A 4 × p 4 }
= 1 21 × Σ p = - 10 10 { A 2 × p 2 }
= A 2 21 × { ( - 10 ) 2 + ( - 9 ) 2 · · · + 10 2 }
≈ A 2 × 36.7 · · · ( 7 )
另外,例如,涡磁场成分XY_Y在图5的例子中只有时间常数100毫秒的一项,通过将拍摄体层的Y坐标值作为式(1)中的p代入而可以计算其振幅LXY_Y。与上述同样地,作为Y坐标值,像-10,-9,……,10那样以1cm刻度抽出21个计测点,以涵盖拍摄区域Img2的Y轴方向的一端到另一端,在各计测点计算LXY_Y,计算将它取平均化而得到的平均振幅AVLXY_Y。在此,由于在图5的例子中A1=0、A3=0、A4=0,所以可以像下式(8)那样计算平均振幅AVLXY_Y
AVL XY _ Y = 1 21 × Σ p = - 10 10 [ A 1 + A 2 × { pol ( p ) } × | p | + A 3 × | p | 2 + A 4 × { pol ( p ) } × | p | 3 ]
= 1 21 × Σ p = - 10 10 [ A 2 × { pol ( p ) } × | p | ]
= A 2 21 × { ( - 1 ) × | 10 | + ( - 1 ) × | 9 | + · · · + 1 × | 9 | + 1 × | 10 | }
= A 2 × 0 · · · ( 8 )
对于其它涡磁场成分,也可以同样地计算平均振幅。由此,可以作为图像整体将倾斜磁场分布的因涡电流产生的偏离量平均化。
另外,进行多体层拍摄时,优选地,针对各拍摄体层的每一个进行上述的倾斜磁场分布的修正。
图9是示出头部的多体层拍摄中的拍摄体层的一例的示意图。在多体层拍摄中,在选择某一拍摄区域的同时从该拍摄区域收集MR信号,然后在该拍摄区域的经过重复时间(TR)为止的等待时间中,发送改变了频率的RF脉冲,选择另一拍摄区域,并且向该另一拍摄区域施加激励脉冲,收集MR信号。这样的话,在一个MR信号收集序列中收集不同的多个拍摄体层的MR信号。
图10是作为针对图9的拍摄体层SL1~SL7的多体层拍摄,示出自旋回波法的脉冲序列的一例的示意图。在图10中,Gss表示体层选择方向倾斜磁场,Gpe表示相位编码方向倾斜磁场,Gro表示读出方向倾斜磁场,Sig表示检测到的MR信号,RF表示90°激励脉冲和180°再会聚脉冲。在图9和图10的例子中,来自拍摄体层SL1~SL7的各MR信号用一个MR信号收集序列收集。
(本实施方式的动作说明)
图11是示出MRI装置20的动作流程的流程图。下面,一边适当参照上述的图1~图10,一边按照图11所示的步骤编号,以利用自旋回波法的多体层拍摄为例,说明包含倾斜磁场分布的修正的MRI装置20的动作。
〔步骤S1〕MPU 86(参照图2)基于通过输入装置62对运算装置60输入的输入信息等进行MRI装置20的初始设定。此处的“输入信息”包含规定“拍摄条件”的至少一部分的输入。作为“拍摄条件”,可举出例如倒(flip?)角、重复时间、体层数、(作为位置信息的)拍摄区域、相位编码方向和频率编码方向的体层数、EPI(回波平面成像)、自旋回波法等的脉冲脉列的种类等。
在此,作为一例,在该初始设定中设定多体层拍摄的拍摄条件。
〔步骤S2〕MRI装置20通过进行预扫描,计算是在什么样的条件下拍摄。例如,进行作为使原子核自旋的纵磁化成分旋转90°的RF脉冲所必需的功率(90°条件)、RF脉冲的中心频率的计算等。
〔步骤S3〕MPU 86基于步骤S1中的输入信息等,确定轴向断面、冠状断面、弧矢断面等的各定位图像的拍摄区域。位置取得部100从MPU86取得各定位图像的拍摄区域,作为在拍摄空间内的位置信息,将它们输入修正部102。
此处的“位置信息”是按照装置坐标系规定拍摄区域的中心位置和三维范围(拍摄区域的三维外缘)的信息。例如,如果作为拍摄板区是有一定厚度的长方体形状的拍摄区域,则也可以将拍摄区域的8个顶点的装置坐标系中的坐标位置作为“位置信息”。或者,是作为拍摄体层的拍摄区域时,也可以将设置了其中心坐标、体层厚度、断面方向和体层的长宽尺寸的信息作为“位置信息”。
在此,作为一例,修正部102以定位图像的拍摄区域的中心坐标为基准,对于27种涡磁场成分,针对每个时间常数全部计算出涡磁场成分的振幅。通过使用在修正部102中预先存储的式(1)、式(2)和这两式的A1~A4等的各参数值(参照图4~图6),用上述的方法计算这些振幅。另外,关于涡磁场成分的振幅,也可以用式(7)、式(8)像前面描述过的那样,作为在拍摄区域的一端侧到另一侧的多个坐标位置平均化了的值算出。
然后,修正部102用上述的方法算出为了得到作为目标的倾斜磁场分布(目标分布)而向倾斜磁场线圈26供给的电流值,将它输入MPU 86。此处的向倾斜磁场线圈26供给的电流值是指例如对向倾斜磁场线圈26供给的电流值进行了修正以施加从作为目标的倾斜磁场分布减去全部涡磁场成分得到的倾斜磁场分布的电流值。
〔步骤S4〕MPU 86向序列控制器56输入包含像上述那样算出的向倾斜磁场线圈26供给的电流值的、定位图像的MR信号收集用的序列。然后,生成定位图像。
具体地说,MPU 86向序列控制器56输入包含脉冲序列的拍摄条件。由用静磁场电源40激励的静磁场磁体22在机架内的拍摄空间中形成静磁场。另外,从调整线圈电源42向调整线圈24供给电流,将拍摄空间中形成的静磁场均匀化。
序列控制器56通过根据被输入的脉冲序列驱动倾斜磁场电源44、RF发送器46和RF接收器48,在包含受检体H的拍摄对象部位的拍摄区域上形成倾斜磁场,并且从RF线圈28产生RF信号。此时,向X轴倾斜磁场线圈26x、Y轴倾斜磁场线圈26y、Z轴倾斜磁场线圈26z供给的电流值是像上述那样进行了修正的值。因此,像上述那样,根据拍摄区域的位置和涡电流磁场的时间常数,使倾斜磁场的波形变形以使涡电流磁场相抵销。即,在拍摄区域中形成大致成为目标的倾斜磁场分布。
然后,受检体H的内部的因核磁共振产生的MR信号由RF线圈28接收,由RF接收器48检测。RF接收器48通过对检测到的MR信号实施预定的信号处理,然后通过对它进行A/D变换,生成数字化了的MR信号即原始数据。RF接收器48将生成的原始数据输入给序列控制器56。序列控制器56将原始数据输入给图像重建部90,图像重建部90在形成在k空间数据库92中的k空间中配置原始数据,作为k空间数据。
图像重建部90从k空间数据库92取入k空间数据,通过对它实施包含傅里叶变换等的图像重建处理,重建图像数据,将得到的图像数据存储到图像数据库94中。
图像处理部96从图像数据库94取入图像数据,通过对它实施预定的处理,生成二维的显示用图像数据,将该显示用图像数据存储到存储装置66中。
〔步骤S5〕显示控制部98根据MPU 86的指令在显示装置64的监视器上显示定位图像。基于该显示图像,(由用户)通过输入装置62对运算装置60输入设定拍摄区域等的信息。
具体地说,例如,在定位图像上用长方形的框输入设定拍摄区域的二维范围。另外,通过输入装置62输入设定例如拍摄区域的厚度。
MPU 86通过取得在此输入的条件,基于例如装置坐标系三维地算出多个拍摄体层的范围和中间位置,作为拍摄空间内的“位置信息”存储起来。另外,这不过是一例,关于“位置信息”,也可以二维地算出。像上述那样,是作为拍摄体层的拍摄区域时,也可以将设置了其中心坐标、体层厚度、断面方向和体层的长宽尺寸的信息作为“位置信息”。
在此作为一例,将图9所示的体层SL1~SL7等作为拍摄区域而设定。
〔步骤S6〕位置取得部100从MPU 86取得各拍摄体层的位置和规定三维范围的“位置信息”,将它们输入修正部102。
修正部102与步骤S3同样地,针对每个拍摄体层,计算出27种涡磁场成分的针对每个时间常数的振幅。此处的振幅可以以各拍摄体层的中心坐标为基准算出,也可以用式(7)和式(8)像前面描述过的那样,作为在拍摄体层内平均了的值算出。然后,修正部102与步骤S3同样地,算出为了得到作为目标的倾斜磁场分布而向倾斜磁场线圈26供给的电流值,输入MPU 86。
〔步骤S7〕在到步骤S6为止设定的多体层拍摄的序列下,从被设定的拍摄区域(拍摄体层SL1~SL7等)收集MR信号。
具体地说,MPU 86向序列控制器56输入包含多体层拍摄的脉冲序列的拍摄条件。与拍摄定位图像时同样地,形成静磁场并均匀化。序列控制器56,通过根据被输入的脉冲序列驱动倾斜磁场电源44、RF发送器46和RF接收器48,在拍摄区域上形成倾斜磁场,并且从RF线圈28产生RF信号。
具体地说,例如,在第一相位编码的MR信号收集中,在拍摄体层SL1的重复时间TR期间收集拍摄体层SL1~SL7这7个体层的MR信号(参照图10)。此时,由于在步骤S6中已经针对各拍摄体层SL1~SL7中的每一个算出了向倾斜磁场线圈26供给的电流的修正值,所以针对各拍摄体层SL1~SL7中的每一个进行上述的倾斜磁场的修正,不管哪个拍摄体层的MR信号都在作为目标的倾斜磁场分布下被收集。此处的倾斜磁场是体层选择方向倾斜磁场Gss、相位编码方向倾斜磁场Gpe和读出方向倾斜磁场Gro。
如果第一相位编码的MR信号的收集结束,则同样地从体层SL1到体层SL7依次进行第二相位编码的MR信号的收集。这样,以重复时间TR的时间间隔每一步骤增加相位编码,同时以相位编码方向的步骤数目反复进行MR信号的收集。
然后,与定位图像的情况同样地,根据由RF线圈28接收的MR信号生成原始数据,在形成于k空间数据库92中的k空间中配置原始数据,作为k空间数据。
〔步骤S8〕图像重建部90从k空间数据库92取入k空间数据,通过对它实施包含傅里叶变换等的图像重建处理,重建所有拍摄体层的图像数据,将得到的图像数据存储到图像数据库94中。
图像处理部96从图像数据库94取入图像数据,通过对它实施预定的图像处理,生成二维的显示用图像数据,将该显示用图像数据存储到存储装置66中。显示控制部98根据MPU 86的指令,在显示装置64的监视器上将显示用图像数据作为MR图像显示。
以上是本实施方式的MRI装置20的动作说明。
(本实施方式的效果)
这样,在本实施方式中,根据拍摄区域(的位置),不仅加入一次成分而且加入二次以上的成分地针对每个时间常数算出涡磁场成分的振幅。然后,通过修正向倾斜磁场线圈26供给的电流以使涡磁场磁场相抵销,从而使倾斜磁场波形变形,因此可以得到大致成为目标的倾斜磁场分布(目标分布)。因此,可以改善倾斜磁场分布畸曲造成的图像质量劣化。
在拍摄区域的中心与磁场中心吻合时,一般情况下倾斜磁场分布畸曲造成的图像质量劣化相对较小,但在本实施方式中,根据拍摄体层的位置,针对每个拍摄体层修正向倾斜磁场线圈26供给的电流。因此,即使在拍摄体层从磁场中心偏离时,也可以通过得到得到大致成为目标的倾斜磁场分布,改善倾斜磁场分布畸曲造成的图像质量劣化。
另外,在多体层拍摄中,基于各拍摄体层的位置,算出每个时间常数的涡磁场成分的振幅,基于该算出结果,针对每个拍摄体层进行向倾斜磁场线圈26供给的电流的修正。因此,即使在多体层拍摄中,由于可以在作为目标的倾斜磁场分布下收集各拍摄体层的MR信号,所以可以改善图像质量劣化。
即,根据以上说明过的实施方式,通过在MRI中简单地补偿因涡电流所产生的磁场成分,可以改善图像质量。
图12是示出不进行基于涡磁场成分的倾斜磁场的修正而拍摄体模(phantom)得到的弧矢断面图像的一例的图。该幻影大致为圆形,其中央具有黑的环状的低信号区域(几乎没有氢原子的区域)。
图13是示出根据本实施方式的方法进行倾斜磁场分布的修正,拍摄与图12相同的体模得到的弧矢断面图像的一例的图。图12、图13都是用FSE(快速自旋回波)法在X偏离中心区拍摄的图像。图像的横方向为Y轴方向。在图12中,从图像中心向左侧跨过黑的环状区域而呈稍黑的部分是图像的劣化部分。它是因FSE法的CPMG(Carr-PurcellMeiboom-Gill序列)条件的偏离造成的灵敏度不均。在图13中,如果着眼于从图像中心向左侧的区域,则可以确认上述的图像的劣化部分被清楚地呈现。
(本实施方式的补充事项)
〔1〕在本实施方式中描述了基于式(1)、式(2)、以及图4~图6,将涡磁场成分的振幅考虑到四次而算出的例子。本发明的实施方式不限于这样的方式。式(1)、式(2)及图4~图6仅是涡磁场成分的振幅的算出方法的一例而已,也可以采用其它的计算方法。另外,也可以是考虑五次以上的成分的计算方法。
另外,虽然设为磁场中心位于装置坐标系的原点,但这仅是一例。即使是磁场中心与装置坐标系的原点不一致的构成,也可以适用本实施方式的方法。
〔2〕在图9~图11中,以利用自旋回波法的多体层拍摄为例说明了说明包含倾斜磁场分布的修正的MRI装置20的动作,但本发明的实施方式不限于这样的方式。关于MR信号收集序列,也可以是利用场回波法等的其它方法的多体层拍摄。或者,也可以不是多体层拍摄,而是倾斜(oblique)拍摄等的其它拍摄序列。
〔3〕作为MRI装置20,描述了在包含静磁场磁体22、调整线圈24、倾斜磁场线圈单元26、RF线圈28的机架外存在RF接收器48的例子(参照图1)。本发明的实施方式不限于这样的方式。也可以是在机架内包含RF接收器48的方式。
具体地说,例如,在机架内配设与RF接收器48相当的电子电路基板。然后,利用该电子电路基板内的预放大器将由接收RF线圈从电磁波变换成模拟电气信号的MR信号放大,作为数字信号输出到机架外,输入序列控制器56。在输出到机架外时,如果用例如光通信缆线作为光数字信号发送,则减轻了外部噪声的影响,所以是优选的。
〔4〕说明一下权利要求书的用语与实施方式的对应关系。另外,以下所示的对应关系仅仅是为了参考而做的一种解释,不对本发明构成限制。
静磁场磁体22、调整线圈24、倾斜磁场线圈26、RF线圈28、控制装置30整体(参照图1)通过伴随倾斜磁场和RF脉冲的施加的拍摄从受检体H收集MR信号的功能是权利要求书中记载的信号收集部的一例。
基于收集到的MR信号生成受检体H的图像数据的图像重建部90、图像数据库94、图像重建部90(参照图2)的功能是权利要求书中记载的图像生成部的一例。
〔5〕虽然说明了本发明的几个实施方式,但这些实施方式都是作为例子提出的,并非用来限定本发明的范围。这些实施方式可以以其它的各种方式实施,在不脱离发明的主要构思的范围内,可以进行各种省略、置换、变更。这些实施方式及其变形都包含在发明的范围和主要构思内,且包含在权利要求书记载的发明及其等价的范围内。

Claims (19)

1.一种磁共振成像装置,包括:
信号收集部,通过在放置受检体的拍摄空间内施加倾斜磁场和RF脉冲,收集从上述受检体产生的磁共振信号;和
图像生成部,基于上述磁共振信号对上述受检体的图像数据进行重建,
其特征在于,还包括:
位置取得部,取得拍摄区域作为在上述拍摄空间内的位置信息;以及
修正部,通过以使得因伴随着上述倾斜磁场的施加产生的涡电流所产生的磁场即涡电流磁场相抵销的方式基于上述位置信息和上述涡电流磁场的时间常数使上述倾斜磁场的波形变形,而使产生的上述倾斜磁场的分布接近目标分布。
2.如权利要求1所述的磁共振成像装置,其特征在于,
还包括产生静磁场的静磁场磁体和产生与被供给的电流对应的上述倾斜磁场的倾斜磁场线圈,
上述修正部预先存储上述涡电流磁场的时间常数,作为基于上述倾斜磁场线圈与上述静磁场磁体的相对位置关系、上述倾斜磁场线圈的形状和材质的值,根据该时间常数使上述倾斜磁场的波形变形。
3.如权利要求1所述的磁共振成像装置,其特征在于:
还包括产生与被供给的电流对应的上述倾斜磁场的倾斜磁场线圈,
上述修正部通过对上述信号收集部向上述倾斜磁场线圈供给的电流进行修正,使上述倾斜磁场的波形变形。
4.如权利要求3所述的磁共振成像装置,其特征在于:
上述修正部根据从上述位置信息得到的上述拍摄区域与磁场中心的距离,进行针对每个上述倾斜磁场的施加方向分别算出上述涡电流磁场的装置坐标系的X轴方向成分的振幅、上述涡电流磁场的装置坐标系的Y轴方向成分的振幅、上述涡电流磁场的装置坐标系的Z轴方向成分的振幅的运算处理,基于上述运算处理的结果修正向上述倾斜磁场线圈供给的电流。
5.如权利要求4所述的磁共振成像装置,其特征在于:
上述修正部通过使上述信号收集部向上述倾斜磁场线圈供给的电流值根据上述涡电流磁场的时间常数随时间变化,使上述倾斜磁场的分布接近上述目标分布。
6.如权利要求4所述的磁共振成像装置,其特征在于:
上述修正部基于上述位置信息取得上述拍摄区域的中心位置,根据上述中心位置与上述磁场中心的距离进行上述运算处理。
7.如权利要求6所述的磁共振成像装置,其特征在于:
上述修正部通过使上述信号收集部向上述倾斜磁场线圈供给的电流值根据上述涡电流磁场的时间常数随时间变化,使上述倾斜磁场的分布接近上述目标分布。
8.如权利要求6所述的磁共振成像装置,其特征在于:
上述信号收集部构成为从与多个图像分别对应的多个上述拍摄区域收集上述磁共振信号;
上述位置取得部构成为取得多个上述拍摄区域作为多组上述位置信息;
上述修正部构成为针对各个上述拍摄区域的每一个修正向上述倾斜磁场线圈供给的电流。
9.如权利要求4所述的磁共振成像装置,其特征在于:
上述修正部通过使用上述拍摄区域内的多个位置处的与上述磁场中心的各距离,进行上述运算处理,以使得上述涡电流磁场的装置坐标系的X轴方向成分、装置坐标系的Y轴方向成分、装置坐标系的Z轴方向成分的各振幅的算出结果成为在整个上述拍摄区域的宽度上平均了的值。
10.如权利要求9所述的磁共振成像装置,其特征在于:
上述修正部通过使上述信号收集部向上述倾斜磁场线圈供给的电流值根据上述涡电流磁场的时间常数随时间变化,使上述倾斜磁场的分布接近上述目标分布。
11.如权利要求9所述的磁共振成像装置,其特征在于:
上述信号收集部构成为从与多个图像分别对应的多个上述拍摄区域收集上述磁共振信号;
上述位置取得部构成为取得多个上述拍摄区域作为多组上述位置信息;
上述修正部构成为针对各个上述拍摄区域的每一个修正向上述倾斜磁场线圈供给的电流。
12.如权利要求3所述的磁共振成像装置,其特征在于:
上述信号收集部构成为从与多个图像分别对应的多个上述拍摄区域收集上述磁共振信号;
上述位置取得部构成为取得多个上述拍摄区域作为多组上述位置信息;
上述修正部构成为针对各个上述拍摄区域的每一个修正向上述倾斜磁场线圈供给的电流。
13.如权利要求3所述的磁共振成像装置,其特征在于:
作为多体层拍摄,上述信号收集部通过施加选择激励某一拍摄区域的上述倾斜磁场和上述RF脉冲,从该拍摄区域收集上述磁共振信号,并且在经过该拍摄区域的重复时间为止的期间,通过施加选择激励另一拍摄区域的上述倾斜磁场和上述RF脉冲,从上述另一拍摄区域收集上述磁共振信号;
上述位置取得部将由上述多体层拍摄选择激励的各个拍摄区域作为上述位置信息;
上述修正部针对各个拍摄区域的每一个算出上述涡电流磁场的装置坐标系的X轴方向成分的振幅、上述涡电流磁场的装置坐标系的Y轴方向成分的振幅、上述涡电流磁场的装置坐标系的Z轴方向成分的振幅,基于上述算出结果针对各个拍摄区域的每一个修正向上述倾斜磁场线圈供给的电流。
14.如权利要求3所述的磁共振成像装置,其特征在于:
上述修正部通过使上述信号收集部向上述倾斜磁场线圈供给的电流值根据上述涡电流磁场的时间常数随时间变化,使上述倾斜磁场的分布接近上述目标分布。
15.如权利要求1所述的磁共振成像装置,其特征在于:
上述位置信息是规定上述拍摄区域的三维范围的信息。
16.如权利要求1所述的磁共振成像装置,其特征在于:
上述位置信息是上述拍摄区域的外缘的各顶点的在装置坐标系中的坐标位置。
17.如权利要求1所述的磁共振成像装置,其特征在于:
上述位置信息是设置了上述拍摄区域的中心坐标、体层厚度、断面方向和体层的长宽尺寸的信息。
18.一种磁共振成像装置的控制装置,该磁共振成像装置通过向倾斜磁场线圈供给电流而在放置受检体的拍摄空间内施加倾斜磁场,伴随上述倾斜磁场执行上述受检体的磁共振成像,
其特征在于包括:
位置取得部,取得上述磁共振成像的拍摄区域,作为在上述拍摄空间内的位置信息;以及
修正部,通过以使得因伴随着上述倾斜磁场的施加产生的涡电流所产生的磁场即涡电流磁场相抵销的方式基于上述位置信息和上述涡电流磁场的时间常数使上述倾斜磁场的波形变形,而使产生的上述倾斜磁场的分布接近目标分布。
19.一种磁共振成像方法,其特征在于包括以下步骤:
取得步骤,取得磁共振成像的拍摄区域,作为在放置受检体的拍摄空间内的位置信息;
修正步骤,通过以使得因伴随着上述倾斜磁场的施加产生的涡电流所产生的磁场即涡电流磁场相抵销的方式基于上述位置信息和上述涡电流磁场的时间常数使上述倾斜磁场的波形变形,而使产生的上述倾斜磁场的分布接近目标分布;
收集步骤,通过在上述拍摄空间内施加通过上述修正步骤变形了的上述倾斜磁场和RF脉冲,收集从上述受检体产生的磁共振信号;以及
重建步骤,基于上述磁共振信号对上述受检体的图像数据进行重建。
CN201110205606.5A 2010-07-23 2011-07-22 磁共振成像方法、磁共振成像装置及其控制装置 Active CN102342832B (zh)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2010166271 2010-07-23
JP2010-166271 2010-07-23
JP2011144834A JP2012040362A (ja) 2010-07-23 2011-06-29 磁気共鳴イメージング方法、磁気共鳴イメージング装置およびその制御装置
JP2011-144834 2011-06-29

Publications (2)

Publication Number Publication Date
CN102342832A true CN102342832A (zh) 2012-02-08
CN102342832B CN102342832B (zh) 2015-07-22

Family

ID=45493097

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110205606.5A Active CN102342832B (zh) 2010-07-23 2011-07-22 磁共振成像方法、磁共振成像装置及其控制装置

Country Status (3)

Country Link
US (1) US8890530B2 (zh)
JP (1) JP2012040362A (zh)
CN (1) CN102342832B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107615080A (zh) * 2015-05-29 2018-01-19 梅约医学教育与研究基金会 用于在具有非对称梯度的磁共振成像中进行伴随场校正的系统和方法
CN110850344A (zh) * 2018-08-21 2020-02-28 西门子医疗有限公司 操作mri装置的方法

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012063654A1 (ja) * 2010-11-09 2012-05-18 株式会社 日立メディコ 磁気共鳴イメージング装置、計測空間座標補正方法、及び、画像再構成方法
JP2012183233A (ja) * 2011-03-07 2012-09-27 Toshiba Corp 磁気共鳴イメージング装置
KR101473872B1 (ko) * 2013-02-05 2014-12-18 삼성전자 주식회사 자기공명영상장치 및 그 제어방법
JP6697261B2 (ja) 2014-12-26 2020-05-20 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置、拡散強調画像の生成方法及び画像処理装置
US10234530B2 (en) 2015-02-23 2019-03-19 Toshiba Medical Systems Corporation Magnetic resonance imaging apparatus
DE102015222114A1 (de) * 2015-11-10 2017-05-11 Siemens Healthcare Gmbh Verfahren zum Ansteuern einer Shimeinheit, Steuereinheit und Magnetresonanzgerät
KR101797674B1 (ko) * 2016-02-12 2017-11-15 삼성전자주식회사 자기 공명 영상 촬영 장치 및 그에 따른 자기 공명 영상 촬영 방법
JP6817775B2 (ja) 2016-10-11 2021-01-20 株式会社東芝 補正装置、補正方法及び磁気共鳴画像装置
TR201914968A1 (tr) * 2017-04-06 2022-02-21 Ihsan Dogramaci Bilkent Ueniversitesi Optimum faz kayması darbe genişlik modülasyonu modelinin uygulanması vasıtasıyla gradyan dizisi sisteminde minimum akım dalgalanmasının algoritması ve uygulanması.

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1521515A (zh) * 2003-02-13 2004-08-18 确定对涡流场进行补偿设置的方法
US20100148774A1 (en) * 2008-12-17 2010-06-17 Mitsukazu Kamata Magnetic resonance imaging apparatus and magnetic resonance imaging method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4703275A (en) * 1985-07-25 1987-10-27 Picker International, Inc. Method and apparatus to compensate for eddy currents in magnetic resonance imaging
US4647858A (en) * 1985-07-29 1987-03-03 General Electric Company Methods for overcoming transient magnetic field inhomogeneity in nuclear magnetic resonance imaging
JPH03195539A (ja) 1989-12-26 1991-08-27 Toshiba Corp 磁気共鳴イメージング装置
US7112964B2 (en) * 2004-08-02 2006-09-26 General Electric Company Eddy current measurement and correction in magnetic resonance imaging systems with a static phantom
US8008915B2 (en) 2007-12-28 2011-08-30 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1521515A (zh) * 2003-02-13 2004-08-18 确定对涡流场进行补偿设置的方法
US20100148774A1 (en) * 2008-12-17 2010-06-17 Mitsukazu Kamata Magnetic resonance imaging apparatus and magnetic resonance imaging method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
P. JEHENSON 等: "Analytical Method for the Compensation of Eddy-Current Effects Induced by Pulsed Magnetic Field Gradients in NMR Systems", 《JOURNAL OF MAGNETIC RESONANCE》 *
马超 等: "磁共振成像中涡流效应的仿真", 《清华大学学报(自然科学版)》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107615080A (zh) * 2015-05-29 2018-01-19 梅约医学教育与研究基金会 用于在具有非对称梯度的磁共振成像中进行伴随场校正的系统和方法
US10712420B2 (en) 2015-05-29 2020-07-14 Mayo Foundation For Medical Education And Research Systems and methods for concomitant field correction in magnetic resonance imaging with asymmetric gradients
CN110850344A (zh) * 2018-08-21 2020-02-28 西门子医疗有限公司 操作mri装置的方法

Also Published As

Publication number Publication date
US8890530B2 (en) 2014-11-18
JP2012040362A (ja) 2012-03-01
CN102342832B (zh) 2015-07-22
US20120019251A1 (en) 2012-01-26

Similar Documents

Publication Publication Date Title
CN102342832A (zh) 磁共振成像方法、磁共振成像装置及其控制装置
JP4114989B2 (ja) 磁気共鳴システムの磁場を補償する方法
US20120146640A1 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN102680930B (zh) 用于设置匀场电流和射频中心频率的方法以及磁共振设备
CN105517495B (zh) 用于校正与磁共振成像无线射频线圈中的发送和接收时钟之间的异步定时偏移相关联的误差的方法和装置
US8942945B2 (en) System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging
US8482279B2 (en) System and method of parallel imaging for magnetic resonance imaging near metallic implants
US7141970B2 (en) Magnetic resonance imaging device
US7847545B2 (en) System and method for correcting flow velocity measurements in phase contrast imaging using magnetic field monitoring
DE19750637A9 (de) Verfahren zur Messung und Kompensation von durch Wirbelströme induzierten sich örtlich und zeitlich ändernden Magnetfeldern
WO2013042466A1 (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JPH0449949A (ja) 磁気共鳴サーモグラフィー方法
Scheffler et al. Spread‐spectrum magnetic resonance imaging
CN103957786A (zh) 磁共振成像装置及磁共振成像中的线圈选择支援方法
CN102670202B (zh) 磁共振成像装置
US11294016B1 (en) System and method for calibration of asymmetric gradient concomitant field correction parameters
CN100589756C (zh) 磁共振摄影装置
EP0135143B1 (en) Imaging method and apparatus using nuclear magnetic resonance
Cho et al. A new nuclear magnetic resonance imaging technique for unambiguous unidirectional measurement of flow velocity
JP2013017811A (ja) 磁気共鳴イメージング装置
US20230358838A1 (en) Method for determining a radio-frequency transmission pulse for a magnetic resonance scan, a magnetic resonance apparatus, and a computer program product
Zhao et al. A fast B1 mapping method for transmit/receive coils for parallel transmit (PTX) applications
JP3317552B2 (ja) Mri装置
CN102309323B (zh) Epi中使用自导航实时相位校正的图像伪影减少系统和方法
JP2011098031A (ja) 磁気共鳴イメージング装置及び画像補正方法

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: 20160804

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.