CN102711602A - 磁共振成像装置和磁共振成像方法 - Google Patents

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

Info

Publication number
CN102711602A
CN102711602A CN2011800030340A CN201180003034A CN102711602A CN 102711602 A CN102711602 A CN 102711602A CN 2011800030340 A CN2011800030340 A CN 2011800030340A CN 201180003034 A CN201180003034 A CN 201180003034A CN 102711602 A CN102711602 A CN 102711602A
Authority
CN
China
Prior art keywords
distortion correction
mentioned
view data
value
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
CN2011800030340A
Other languages
English (en)
Other versions
CN102711602B (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 CN102711602A publication Critical patent/CN102711602A/zh
Application granted granted Critical
Publication of CN102711602B publication Critical patent/CN102711602B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341—Diffusion imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • 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/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/56518—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Vascular Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Image Processing (AREA)

Abstract

实施方式的磁共振成像装置具备成像单元和失真校正单元。成像单元通过伴随施加轴相同且与不同的b值对应的MPG脉冲的施加来执行扩散强调成像，从而收集与上述不同的b值对应的多个扩散强调图像数据。失真校正单元基于与上述多个扩散强调图像数据中的成为失真校正的对象的扩散强调图像数据所对应的b值不同的b值所对应的扩散强调图像数据，计算出针对成为上述失真校正的对象的扩散强调图像数据的失真校正系数，通过使用所计算出的上述失真校正系数来执行成为上述失真校正的对象的扩散强调图像数据的失真校正，从而生成被实施了上述失真校正的图像数据。

Description

磁共振成像装置和磁共振成像方法
技术领域
本发明实施方式涉及磁共振成像(MRI：Magnetic Resonance Imaging)装置和磁共振成像方法。
背景技术
MRI是如下这样的摄像法：使置于静磁场中的被检测体的原子核自旋利用拉莫频率(Larmor frequency)的高频(RF：radio frequency(射频))信号进行磁激励，根据伴随该激励所产生的磁共振(MR：magneticresonance)信号来重建图像。
在该磁共振成像的领域中，有称为回波平面成像(EPI：echo planarimaging)的摄像法。EPI是MRI中的高速摄像法的一种，其是对于1次核磁激励使倾斜磁场高速地连续反转，使回波连续产生来进行扫描。更具体地说，在EPI中，在施加了激励脉冲(FLIP PULSE)后，在xy平面内的磁化因横缓和(T2缓和)而衰减并消失之前，利用相位编码(PE：phaseencode)的步骤来产生连续的梯度回波，收集图像重建所需的全部的数据。在EPI中，有使用对在激励脉冲和聚相脉冲(FLOP PULSE)之后产生的自旋回波信号进行收集的自旋回波法(SE：spin echo)的SE EPI、以及使用对在激励脉冲的施加后产生的回波信号进行收集的场回波法(FE：fieldecho)的FE EPI或使用FFE(Fast FE)法的FFE EPI。此外，结合施加历经多次的激励脉冲而获得的回波链(echo train)的数据来做成1张的量的图像数据的EPI，被称为多发EPI，与其相比，仅由1次激励脉冲的施加来重建图像的EPI，被称为单发(SS：single shot)EPI。此外，还有称为混合EPI(Hybrid EPI)的EPI。
此外，作为EPI的应用技术，已知有扩散强调成像(DWI：diffusionweighted imaging(又称弥散加权成像))。DWI是如下这样的摄像法：通过施加称为MPG(motion probing gradient：运动梯度场)脉冲的较强强度的倾斜磁场，从而对由摄像对象的活动引起的相位移动进行强调，获得强调了摄像对象的扩散效果的图像。一般来说，基于改变MPG脉冲的施加方向而收集的多个扩散强调图像(DWI：diffusion weighted image)和未施加MPG脉冲而收集的基准图像，制作用于诊断的扩散系数(ADC：ApparentDiffusion Coefficient(表观扩散系数))图像、各向同性的(isotropic)DWI等的参数图像。
进而，作为DWI的应用，已知有扩散张量成像(DTI：diffusion tensorimaging)。一般来说，DTI是如下这样的摄像法：改变MPG脉冲的施加方向来进行多个DWI的收集，进行基于DWI的张量解析，在数学上对摄像部位进行图像化。
在DWI和DTI中使用的EPI中，由于强的MPG脉冲而感应出涡流。因此，有时会因涡流引起的不均一的磁场的影响而使图像产生失真。该图像的失真依赖于MPG脉冲的施加定时、大小和方向。
因此，不但在改变MPG脉冲的施加方向、大小而收集的DWI之间，有时还在不施加MPG脉冲而收集的基准图像与DWI之间产生由与MPG脉冲的条件相应的图像失真引起的位置偏移。当基于有这样的位置偏移的图像而生成ADC图像时，会产生伪影(artifact)，会导致分辨率的降低。此外，当基于具有失真的图像而生成各向同性DWI(isotropic DWI)时，会成为模糊的图像。由于在DTI的情况下，也会在改变MPG脉冲的施加方向而收集到的图像数据中产生不同的失真，所以有时难以进行正确的张量解析。
根据这样的背景，考虑了在图像空间上对由MPG脉冲引起的图像的失真进行校正的各种各样的技术。
现有技术文献
专利文献
专利文献1：日本特开2007-167374号公报；
专利文献2：日本特开2009-195584号公报。
发明内容
发明要解决的问题
在DWI和DTI中，进一步良好地校正图像的失真是个问题。
本发明的目的是，提供一种能获得更良好地实施了失真校正的DWI或者DTI的磁共振成像装置和磁共振成像方法。
用于解决问题的方案
本发明实施方式的磁共振成像装置具备成像单元和失真校正单元。成像单元通过伴随施加轴相同且与不同的b值对应的MPG脉冲的施加来执行扩散强调成像，从而收集与上述不同的b值对应的多个扩散强调图像数据。失真校正单元基于与上述多个扩散强调图像数据中的成为失真校正的对象的扩散强调图像数据所对应的b值不同的b值所对应的扩散强调图像数据，计算出针对成为上述失真校正的对象的扩散强调图像数据的失真校正系数，通过使用计算出的上述失真校正系数来执行成为上述失真校正的对象的扩散强调图像数据的失真校正，从而生成被实施了上述失真校正的图像数据。
此外，本发明实施方式的磁共振成像方法具有如下步骤：通过伴随施加轴相同且与不同的b值对应的MPG脉冲的施加来执行扩散强调成像，从而收集与上述不同的b值对应的多个扩散强调图像数据的步骤；以及基于与上述多个扩散强调图像数据中的成为失真校正的对象的扩散强调图像数据所对应的b值不同的b值所对应的扩散强调图像数据，计算出针对成为上述失真校正的对象的扩散强调图像数据的失真校正系数，通过使用计算出的上述失真校正系数来执行成为上述失真校正的对象的扩散强调图像数据的失真校正，从而生成被实施了上述失真校正的图像数据的步骤。
附图说明
图1是本发明实施方式的磁共振成像装置的构成图。
图2是图1所示的计算机的功能框图。
图3是表示在图2所示的摄像条件设定部中设定的EPI序列的一个例子的图。
图4是表示在图2所示的摄像条件设定部中设定的2个DWI序列的一个例子的图。
图5是说明在图2所示的失真校正部中求取失真校正系数并进行失真校正的第一方法的图。
图6是说明在图2所示的失真校正部中求取失真校正系数并进行失真校正的第二方法的图。
图7是说明在图2所示的失真校正部中求取失真校正系数并进行失真校正的第三方法的图。
图8是表示由图1所示的磁共振成像装置伴随失真校正系数的计算处理来进行DWI时的流程的流程图。
图9是表示由图1所示的磁共振成像装置伴随使用预先取得的失真校正系数的失真校正处理来进行DWI时的流程的流程图。
具体实施方式
参照附图，对本发明实施方式的磁共振成像装置和磁共振成像方法进行说明。
图1是本发明实施方式的磁共振成像装置的构成图。
磁共振成像装置20具备形成静磁场的筒状的静磁场用磁铁21、设在该静磁场用磁铁21的内部的匀场线圈(shim coil)22、倾斜磁场线圈23和RF线圈24。
此外，在磁共振成像装置20中配备有控制系统25。控制系统25具备静磁场电源26、倾斜磁场电源27、匀场线圈电源28、发射器29、接收器30、序列控制器31和计算机32。控制系统25的倾斜磁场电源27包括X轴倾斜磁场电源27x、Y轴倾斜磁场电源27y和Z轴倾斜磁场电源27z。此外，在计算机32中配备有输入装置33、显示装置34、运算装置35和存储装置36。
静磁场用磁铁21与静磁场电源26连接，具有利用从静磁场电源26供给的电流而使摄像区域形成静磁场的功能。另外，静磁场用磁铁21多是由超导线圈构成，一般是在励磁时与静磁场电源26连接并被供给电流，在一旦励磁之后成为非连接状态。此外，也有时用永磁铁来构成静磁场用磁铁21，不设置静磁场电源26。
此外，在静磁场用磁铁21的内侧在同轴上设有筒状的匀场线圈22。匀场线圈22构成为：与匀场线圈电源28连接，从匀场线圈电源28对匀场线圈22供给电流，使静磁场均一化。
倾斜磁场线圈23包括X轴倾斜磁场线圈23x、Y轴倾斜磁场线圈23y和Z轴倾斜磁场线圈23z，并在静磁场用磁铁21的内部形成为筒状。在倾斜磁场线圈23的内侧设有床37并作为摄像区域，在床37放置被检测体P。在RF线圈24中，有内置于吊架(gantry)中的RF信号收发用的全身用线圈(WBC：whole body coil)和设于床37、被检测体P附近的RF信号接收用的局部线圈等。
此外，倾斜磁场线圈23与倾斜磁场电源27连接。倾斜磁场线圈23的X轴倾斜磁场线圈23x、Y轴倾斜磁场线圈23y和Z轴倾斜磁场线圈23z分别与倾斜磁场电源27的X轴倾斜磁场电源27x、Y轴倾斜磁场电源27y和Z轴倾斜磁场电源27z连接。
而且，构成为：利用从X轴倾斜磁场电源27x、Y轴倾斜磁场电源27y和Z轴倾斜磁场电源27z分别对X轴倾斜磁场线圈23x、Y轴倾斜磁场线圈23y和Z轴倾斜磁场线圈23z供给的电流，能在摄像区域分别形成X轴向的倾斜磁场Gx、Y轴向的倾斜磁场Gy和Z轴向的倾斜磁场Gz。
RF线圈24与发射器29和接收器30的至少一方连接。发射用的RF线圈24具有从发射器29接受RF信号并将其发射给被检测体P的功能，接收用的RF线圈24具有：接收伴随被检测体P内部的原子核自旋的RF信号的激励而产生的MR信号并将其提供给接收器30的功能。
另一方面，控制系统25的序列控制器31与倾斜磁场电源27、发射器29和接收器30连接。序列控制器31具有如下功能：存储为了使倾斜磁场电源27、发射器29和接收器30驱动所需的控制信息，例如记述了应施加于倾斜磁场电源27的脉冲电流的强度、施加时间、施加定时等的动作控制信息的序列信息的功能；以及通过按照所存储的规定的序列使倾斜磁场电源27、发射器29和接收器30驱动，而使X轴倾斜磁场Gx、Y轴倾斜磁场Gy、Z轴倾斜磁场Gz和RF信号产生的功能。
此外，序列控制器31构成为：接受由接收器30中的MR信号的检波和A/D(analog to digital)变换而获得的复数数据即原始数据(raw data)，并将其提供给计算机32。
因此，在发射器29中具备基于从序列控制器31接受的控制信息而将RF信号提供给RF线圈24的功能，另一方面，在接收器30中具备下述功能：通过对从RF线圈24接受的MR信号进行检波并执行所需要的信号处理并且对其进行A/D变换，从而生成作为数字化后的复数数据的原始数据的功能；以及将所生成的原始数据提供给序列控制器31的功能。
此外，通过运算装置35来执行在计算机32的存储装置36中保存的程序，从而在计算机32中具备各种功能。但是，也可以代替程序的至少一部分，而将具有各种功能的特定的电路设置于磁共振成像装置20中。
图2是图1所示的计算机32的功能框图。
计算机32的运算装置35通过执行在存储装置36中保存的程序而作为摄像条件设定部40和数据处理部41发挥功能。数据处理部41具有图像重建部41A、失真校正部41B和图像处理部41C。此外，存储装置36作为k空间数据存储部42、图像数据存储部43和失真系数存储部44发挥功能。
摄像条件设定部40具有如下功能：基于从输入装置33输入的信息，对包含EPI序列的摄像条件进行设定，并将所设定的摄像条件提供给序列控制器31，由此进行驱动控制。
图3是表示在图2所示的摄像条件设定部40中设定的EPI序列的一个例子的图。
在图3中横轴表示时间，RF表示RF发射脉冲，GSS表示切片选择(SS：slice selection)倾斜磁场脉冲，GPE表示相位编码(PE：phase encode)倾斜磁场脉冲，GRO表示读出(RO：readout)倾斜磁场脉冲，ECHO表示MR接收回波信号。
如图3所示那样，EPI序列是如下这样的序列，即，在利用90°RF脉冲和180°RF脉冲的施加激励了由SS倾斜磁场脉冲的施加而选择出的切片后，反复施加PE倾斜磁场脉冲和RO倾斜磁场脉冲，由此根据所选择的切片来连续地收集多个MR接收回波信号的脉冲序列。从90°RF脉冲的施加时刻起至k空间中心的回波信号的收集时刻为止的时间被称为有效回波时间(effective TE：effective echo time)。
此外，摄像条件设定部40具有设定脉冲序列的功能，其中该脉冲序列是用于将作为MPG脉冲的效果的指标的b值(b factor)设定为各种值并进行DWI的序列。即，在摄像条件设定部40中，能设定用于收集与相互不同的b值对应的多个DWI数据的摄像条件。特别是，摄像条件设定部40还具备如下功能：设定伴随施加轴和b值的绝对值相互相同且极性相互不同的MPG脉冲的施加的2个DWI序列的功能。
图4是表示在图2所示的摄像条件设定部40中设定的2个DWI序列的一个例子的图。
在图4(A)、(B)中横轴表示时间，RF表示RF发射脉冲，GSS表示SS倾斜磁场脉冲，GPE表示PE倾斜磁场脉冲，GRO表示RO倾斜磁场脉冲，ECHO表示MR接收回波信号。
图4(A)表示伴随使PE方向为施加轴的极性为正的MPG脉冲的施加的DWI序列的例子。此外，图4(B)表示伴随使PE方向为施加轴的极性为负的MPG脉冲的施加的DWI序列的例子。如图4(A)、(B)所示那样，能在MR接收回波信号的收集、也就是说PE倾斜磁场脉冲的施加之前，将伴随施加轴和面积相同且极性分别为正和负的MPG脉冲(+MPGPULSE，-MPG PULSE)的施加的2个DWI序列作为摄像条件，在摄像条件设定部40中进行设定。
另外，图4(A)、(B)虽然表示MPG脉冲的施加轴为PE方向的例子，但也有时根据诊断目的而将其他方向设为施加轴。
此外，DWI数据能通过取多个图像数据的平均的平均处理来将其收集为AVERAGE图像数据。只要进行平均处理，就能获得使SNR(signal tonoise ratio)提高的DWI数据。因此，作为为了取平均而进行加法运算的图像数据的数的AVERAGE数、与成为加法对象的DWI数据对应的MPG脉冲的施加轴和成为加法对象的DWI数据的b值，被作为摄像条件由摄像条件设定部40进行设定。
数据处理部41具有如下功能：从序列控制器31接受MR信号并将其作为k空间数据配置在形成于k空间数据存储部42的k空间中的功能；利用包含针对在k空间数据存储部42中保存的k空间数据的图像重建处理和失真校正处理的数据处理，生成实施了失真校正的图像数据的功能；将所生成的图像数据写入图像数据存储部43的功能；求取在DWI数据的失真校正中使用的失真校正系数并将其写入失真系数存储部44的功能；以及使显示装置34显示图像数据的功能。
图像重建部41A具有如下功能：通过从k空间数据存储部42取入k空间数据并实施包含傅立叶变换(FT：Fourier transform)的图像重建处理，从而重建图像数据。施加MPG脉冲而收集的图像数据成为DWI数据，各DWI数据成为与b值和MPG脉冲的施加轴对应的图像数据。另外，不是DWI数据的图像数据相当于b值为0的图像数据。
此外，图像重建部41A具有如下功能：根据需要，通过取将MPG脉冲的施加轴和b值设定为相同而收集的多个DWI数据的平均的第一平均处理，生成AVERAGE图像数据。即，能通过对利用相同的DWI的条件收集的DWI数据进行加法运算，从而使与特定的施加轴和b值对应的DWI数据的SNR提高。
失真校正部41B具有如下功能：求取用于对伴随MPG脉冲的施加而收集的DWI数据进行失真校正的失真校正系数的功能；以及通过使用所求取的失真校正系数或者参照失真系数存储部44而取得的失真校正系数，进行针对成为失真校正的对象的DWI数据的失真校正，从而生成被实施了失真校正的图像数据的功能。
特别是，失真校正部41B构成为：基于与由DWI收集的多个DWI数据中的成为失真校正对象的DWI数据所对应的b值不同的b值所对应的DWI数据，计算出针对成为失真校正的对象的DWI数据的失真校正系数。此外，失真校正系数的计算能根据需要不施加MPG脉冲地还基于不是使b值为0而收集的DWI的图像数据来进行。
更具体地说，能基于施加轴和b值的绝对值相互相同且与不同极性的MPG脉冲对应的2个DWI数据，更高精度地求取针对2个DWI数据的一方的失真校正系数。
图5是说明在图2所示的失真校正部41B中求取失真校正系数并进行失真校正的第一方法的图。
在图5(A)、(B)、(C)、(D)中各纵轴表示x轴。例如，当将b值设定为0、+500、-500并分别收集图像数据时，如图5(A)所示那样，使b值为0而收集的基准图像数据Ib0由于没有MPG脉冲的影响，所以理想的是成为没有失真的图像数据。另一方面，b值与+500对应的DWI数据Ib+500和b值与-500对应的DWI数据Ib-500分别受到MPG脉冲的影响而成为具有失真的图像数据。
此外，b值与+500对应的DWI数据Ib+500和b值与-500对应的DWI数据Ib-500由于分别受到相互相反的极性的MPG脉冲的影响，所以若能忽略其他影响，则失真方向的极性关于原点呈相互对称。另外，图5(A)示出为了简化说明而在DWI数据Ib+500、Ib-500产生x轴向的1维失真的情况的例子。不限于b值为±500的情况，与绝对值相同且符号相反的b值对应的2个DWI数据若能忽略其他影响，则失真方向关于原点呈对称，失真量变得相同。
通过对于与这样的0、+500、-500的b值对应的各图像数据Ib0、Ib+500、Ib-500进行使用阈值的二值化处理，从而如图5(B)所示那样，分别生成与0、+500、-500的b值对应的二值化图像数据I′b0、I′b+500、I′b-500。
接着，一边对与+500的b值对应的二值化DWI数据I′b+500的各图像点进行坐标变换使其移动，一边计算与0的b值对应的二值化基准图像数据I′b0的各图像点中的像素值和与+500的b值对应的二值化DWI数据I′b+500的各图像点中的像素值之间的平方差的总和。即，将表示二值化DWI数据I′b+500的坐标变换的各系数作为参数，执行使二值化DWI数据I′b+500相对于二值化基准图像数据I′b0的像素值的平方误差为最小的最优化处理。
坐标变换导致能简化成为以平行移动、旋转移动和伸缩变形的组合表示的仿射变换的处理。根据仿射变换的坐标变换由式(1)来表示。
[数学式1]
f ( x , y , z ) = x ′ y ′ z ′ = p 0 p 1 p 2 p 3 p 4 p 5 p 6 p 7 p 8 x y z + T x T y T z - - - ( 1 )
其中，在式(1)中，f是表示坐标变换的函数，(x，y，z)是坐标变换前的图像点的空间坐标，(x′，y′，z′)是坐标变换后的图像点的空间坐标，p0、p1、p2、...、p8、Tx、Ty、Tz是坐标变换的系数。
当将式(1)的各坐标变换系数作为参数，进行使二值化DWI数据I′b+500的像素值的平方误差为最小的最优化处理时，如图5(C)所示那样生成使相对于二值化基准图像数据I′b0的位置偏移为最小的二值化校正DWI数据I′b+500_cor，并且获得多个坐标变换系数并将其做成矩阵。即，能获得基准图像数据Ib0和与正的b值对应的DWI数据Ib+500间的二值化处理和坐标变换处理后的各像素值之差分值的总和为最小的坐标变换系数的组合K1。另外，最优化处理成为包含坐标变换系数的逐次探索处理和像素值的内插处理的处理。
能以相同的方法，与基准图像数据Ib0和与负的b值对应的DWI数据Ib-500间的二值化处理和坐标变换处理后的各像素值之差分值的总和为最小的坐标变换系数的组合K2一起，还获得二值化校正DWI数据I′b+500_cor。
而且，能选择针对二值化基准图像数据I′b0的像素值的平方差的总和变得较小的二值化校正DWI数据(I′b+500_cor或者I′b-500_cor)的生成中所使用的坐标变换系数的组合(K1或者K2)，并将其用作失真校正系数。
进而，通过使用失真校正系数对针对二值化基准图像数据I′b0的像素值的平方差的总和变得较小的二值化校正DWI数据(I′b+500_cor或者I′b-500_cor)的生成中所使用的DWI数据(Ib+500或者Ib-500)进行坐标变换，从而能如图5(D)所示那样获得失真校正后的DWI数据Ib500_cor。
图6是说明在图2所示的失真校正部41B中求取失真校正系数并进行失真校正的第二方法的图。
在图6(A)、(B)、(C)、(D)中各纵轴表示x轴。通过与图5(A)、(B)同样地，收集如图6(A)所示那样的无失真的基准图像数据Ib0和b值与±500对应的失真方向相互相反的2个DWI数据Ib+500、Ib 500并进行二值化处理，从而能如图6(B)所示那样，分别生成与0、+500、-500的b值对应的二值化图像数据I′b0、I′b+500、I′b-500。
在此，假设当使用使与+500的b值对应的二值化DWI数据I′b+500的相对于二值化基准图像数据I′b0的各像素值的平方差的总和为最小的最优的各坐标变换系数变为2倍的坐标变换系数的组合2*K1，来进行与+500的b值对应的二值化DWI数据I′b+500的坐标变换时，由坐标变换生成的二值化变换DWI数据I′b+500_double在理论上变得和与-500的b值对应的二值化DWI数据I′b+500相等。
因此，将与+500的b值对应的二值化DWI数据I′b+500的各坐标变换系数作为参数，执行使SumError1+SumError2为最小的最优化处理，其中，该SumError1+SumError2是，与+500的b值对应的二值化DWI数据I′b+500的相对于二值化基准图像数据I′b0的各像素值的平方差的总和SumError1、和由2倍的坐标变换系数的组合进行坐标变换的二值化变换DWI数据I′b+500_double和与-500的b值对应的二值化DWI数据I′b+500之间的各像素值的平方差的总和SumError2之和。
该最优化处理的结果是，生成如图6(C)所示的那样的二值化校正DWI数据I′b+500_cor2和二值化变换DWI数据I′b+500_double。此外，能获得坐标变换系数的组合K12，该组合K12是使与+500的b值对应的二值化DWI数据I′b+500和与0的b值对应的二值化基准图像数据I′b0之间的位置偏移、以及与+500的b值对应的二值化DWI数据I′b+500和与-500的b值对应的二值化DWI数据I′b+500之间的实质位置偏移的双方为最小的坐标变换系数的组合。
而且，能将所获得的坐标变换系数的组合K12设定为针对与+500的b值对应的DWI数据Ib+500的失真校正系数。进而，通过使用失真校正系数对与+500的b值对应的DWI数据Ib+500进行坐标变换，从而能获得如图6(D)所示那样的失真校正后的DWI数据Ib500_cor。
另外，虽然在上述的例子中，针对将与+500的b值对应的二值化DWI数据I′b+500的各坐标变换系数作为参数进行最优化处理的情况进行了说明，但是也能同样地通过进行将与-500的b值对应的二值化DWI数据I′b-500的各坐标变换系数作为参数的最优化处理，从而求取失真校正系数。
此外，还能不使用b值为0的图像数据地求取失真校正系数。
图7是说明在图2所示的失真校正部41B中求取失真校正系数并进行失真校正的第三方法的图。
在图7(A)、(B)、(C)、(D)中各纵轴表示x轴。与图5(A)、(B)同样地，通过收集如图7(A)所示那样的b值与±500对应的失真方向相互相反的2个DWI数据Ib+500、Ib-500并进行二值化处理，从而能如图7(B)所示那样，分别生成与+500、-500的b值对应的二值化图像数据I′b+500、I′b-500。
接着，将与-500的b值对应的二值化图像数据I′b-500的各坐标变换系数作为参数，执行使与+500的b值对应的二值化DWI数据I′b+500和与-500的b值对应的二值化图像数据I′b-500之间的各像素值的平方差的总和为最小的最优化处理。
该最优化处理的结果是，生成如图7(C)所示那样的与-500的b值对应的二值化图像数据I′b-500的坐标变换后的二值化变换DWI数据I′b- 500_double。此外，能获得使与+500的b值对应的二值化DWI数据I′b+500和与-500的b值对应的二值化DWI数据I′b+500之间的实质位置偏移为最小的坐标变换系数的组合K21。
于是，在二值化变换DWI数据I′b-500_double的生成中所使用的各坐标变换系数的组合K21，相当于使与+500的b值对应的二值化DWI数据I′b+500的相对于二值化基准图像数据I′b0的各像素值的平方差的总和为最小的最优的各坐标变换系数的2倍。
因此，能将在二值化变换DWI数据I′b-500_double的生成中所使用的各坐标变换系数的1/2的坐标变换系数的组合(K21)/2设定为针对与+500的b值对应的DWI数据Ib+500的失真校正系数。进而，通过使用失真校正系数并对与+500的b值对应的DWI数据Ib+500进行坐标变换，从而能获得如图7(D)所示那样的失真校正后的DWI数据Ib500_cor。
另外，虽然在图5、图6和图7所示的第一、第二和第三方法中，使最优化处理的对象为像素值的平方差的总和，但也可以将表示二值化图像数据间的位置偏移的相关系数等的其他指标作为最优化处理的对象。此外，也可以不进行作为边缘检测处理的二值化处理地将针对DWI数据Ib+500、Ib -500的坐标变换系数作为参数来执行使图像数据间的位置偏移的指标最小化的最优化处理。
或者，也可以通过包含二值化处理或者不伴有二值化处理的任意的边缘检测处理，检测包含DWI数据的多个图像数据的各边缘，进行使表示所检测的边缘间的位置偏移的指标为最小的最优化处理。而且，能通过包含针对2个DWI数据的边缘检测的处理，来计算出失真校正系数。
因此，图5所示的第一方法可说是如下这样的方法：将2个DWI数据或者基于2个DWI数据生成的针对2个图像数据的各坐标变换系数作为参数，通过使2个DWI数据或者2个图像数据的相对于使b值为0而收集的图像数据的位置偏移的指标为最小的最优化处理，来计算各坐标变换系数，使与更小的指标对应的坐标变换系数为失真校正系数。
此外，图6所示的第二方法可说是如下这样的方法：将针对2个DWI数据的一方或者基于2个DWI数据生成的2个图像数据的一方的坐标变换系数作为参数，基于通过最优化处理而计算的坐标变换系数，来决定失真校正系数所述最优化处理是使上述一方的相对于使b值为0而收集的图像数据的位置偏移的指标和相对于另一方的位置偏移的指标为最小的处理。
进而，图7所示的第三方法可说是如下这样的方法：将2个DWI数据的一方或者基于2个DWI数据而生成的针对2个图像数据的一方的坐标变换系数为参数，基于通过使一方相对于另一方的位置偏移的指标为最小的最优化处理而计算出的坐标变换系数，来决定失真校正系数。
至此，虽然对收集与绝对值相同且符号相互相反的b值对应的2个DWI数据，失真校正部41B基于2个DWI数据来计算出针对2个DWI数据的一方的失真校正系数的例子进行了说明，但是也能在失真校正部41B中，与绝对值和与成为校正对象的DWI数据对应的b值不同的b值对应，而且基于MPG脉冲的施加轴相同的DWI数据来求取失真校正系数。以下，对收集绝对值与不同的b值对应的多个DWI数据，且失真校正部41B基于与成为基准的b值对应的失真校正系数来计算与比成为基准的b值大的b值对应的失真校正系数的例子进行说明。
具体地说，能基于b值与各失真校正系数K(b)的关系式K(b)＝f(b)，求取与各b值对应的失真校正系数K(b)。关系式的系数能根据与系数的数对应的成为数的基准的b值(bref)以及与成为基准的b值(bref)对应的失真校正系数Kref来求取。
b值与各失真校正系数K(b)的关系式能利用高次式、指数函数、对数函数等任意的函数来进行近似。在近似使用的适当函数的种类和用高次式近似的情况下能否分别忽略2次以上的各项，依赖于是否通过控制MPG脉冲的振幅而调整b值、是否通过控制MPG脉冲的施加时间而调整b值或者是否通过控制MPG脉冲的振幅和施加时间的双方来调整b值等这样的MPG脉冲的设计方法。
因此，通过使用MPG脉冲的条件的模拟、使用装置的安置时或者MPG脉冲的设计时的幻像的试验摄像或者利用伴随MPG脉冲的施加的预扫描而收集的图像数据的图像处理，能预先决定表示b值与各失真校正系数K(b)的关系的近似式。
在b值与各失真校正系数视为正比的情况下，如式(2)所示那样，能基于与成为1组的基准的b值(bref)对应的失真校正系数Kref来求取与各b值对应的失真校正系数K(b)。
K(b)＝(Kref/bref)b    (2)
成为基准的b值(bref)不是0且成为能设定为DWI的摄像条件的较小的值，这从以与失真小的DWI对应的失真校正系数Kref为基准，精度良好地求取与各b值对应的失真校正系数K(b)的观点来看是优选的。即，能基于与更小的b值对应的失真校正系数，推定与更大的b值对应的失真校正系数。例如，在收集b＝0的图像数据和b＝500、b＝1000的DWI数据的情况下，能基于与作为b＝500而收集的DWI数据对应的失真校正系数，利用式(2)计算与作为b＝1000收集的DWI数据对应的失真校正系数。
与成为基准的b值(bref)对应的失真校正系数Kref，如参照图5、图6和图7如上所述那样，能基于与符号不同的b值(+bref、-bref)对应的2个DWI数据Ib+ref、Ib-ref或者与符号不同的b值(+bref、-bref)对应的2个DWI数据Ib+ref、Ib-ref和b值为0的基准图像数据Ib0来求取。
或者，也可以是基于与正或者负的单一的符号的b值(bref)对应的DWI数据Ibref和b值为0的基准图像数据Ib0，求取与成为基准的b值(bref)对应的失真校正系数Kref。在这种情况下，能将对DWI数据Ibref进行二值化的二值化DWI数据I′bref的坐标变换系数作为参数，通过使二值化DWI数据I′bref的相对于二值化基准图像数据I′b0的各像素值的平方差的总和为最小的最优化计算，求取与成为基准的b值(bref)对应的失真校正系数Kref来作为坐标变换系数。
这样在失真校正部41B中计算出的失真校正系数能存储于失真系数存储部44中。即，能将与MPG脉冲的施加轴和每个b值对应的失真校正系数保存在失真系数存储部44中，失真校正部41B能从失真系数存储部44读入并取得与DWI的摄像条件对应的失真校正系数。此外，也可以是能将根据b值求取失真校正系数的函数本身或者函数的系数保存在失真系数存储部44中，能根据b值来计算失真校正系数。
因此，基于在使用幻像(phantom)的试验摄像、过去的摄像中获得的DWI，求取失真校正系数，只要将求取的失真校正系数保存在失真系数存储部44中，就能仅通过失真校正部41B参照失真系数存储部44来取得失真校正系数。其结果是，能省略在各摄像时利用失真校正部41B的失真校正系数的计算处理。
图像处理部41C具有如下功能：对于在数据处理部41中生成的图像数据或者从图像数据存储部43读入的图像数据，实施用于图像显示所需的图像处理的功能；以及使显示装置34显示图像处理后的显示用的图像数据的功能。
例如，能根据需要对失真校正后的多个DWI数据实施第二平均处理。由此，能获得具有诊断所需的SNR的失真校正后的DWI数据。
在图像重建部41A中执行的第一平均处理对失真校正前的DWI数据实施。因此，第一平均处理的对象是，将失真量同等的、MPG脉冲的施加轴和b值相同的多个DWI数据作为对象为前提条件的。
与此相对，第二平均(averaging)处理对失真校正后的多个DWI数据实施。因此，第二平均处理能对与相互不同的b值对应的多个DWI数据来进行。因此，在通过将MR信号的绝对值作为信号值的图像重建处理来生成DWI数据的情况下，优选将与绝对值相同且极性不同的b值等相互不同的b值对应的多个DWI数据作为对象，来执行第二平均处理。由此，能为了SNR的提高而利用为了求取失真校正系数而追加收集的DWI数据。
即，在为了失真校正系数的计算而收集与诊断用非必需的b值对应的DWI数据的情况下，能为了获得诊断用的DWI数据所需的SNR而使用。其结果是，能减少为了获得需要的SNR而成为加法对象的、与诊断用所需的b值对应的DWI数据的数。
另一方面，在为了失真校正系数的计算而收集与诊断用所需的b值对应的DWI数据的情况下，能保持原样地在诊断用的DWI数据的生成用中使用。换言之，能将为了诊断用而本来收集的DWI数据用作失真校正系数的计算用的DWI数据。
因此，能避免伴随上述那样的失真校正部41B中的失真校正系数的计算的实质的数据收集时间和数据收集量的增加。换言之，能使每个单位数据收集时间的SNR提高率为恒定。
第二平均处理可以是将MPG脉冲的施加轴和b值相同的多个DWI数据作为对象来执行。此外，也可以是不执行第一和第二平均处理的一方或者双方。因此，也可以是进行取失真校正前和上述失真校正后的规定的多个DWI数据的至少一方的平均的平均处理。
而且，通过具有这样的构成的磁共振成像装置20，能伴随施加轴相同且与不同的b值对应的MPG脉冲的施加来执行DWI，能收集与不同的b值对应的多个DWI数据。此外，还能使b值为0来收集图像数据。
接着，对磁共振成像装置20的动作和作用进行说明。首先，对通过磁共振成像装置20伴随失真校正系数的计算处理来进行DWI的情况进行说明。
图8是表示通过图1所示的磁共振成像装置20伴随失真校正系数的计算处理来进行DWI时的流程的流程图。
首先在步骤S1中，摄像条件设定部40基于从输入装置33输入的信息，根据如图3所示那样使b值为0的EPI序列和诊断目的，设定与需要的b值对应的DWI序列。进而，摄像条件设定部40设定绝对值与最小的b值相同且符号相反的b值所对应的DWI序列。即，在摄像条件设定部40中设定摄像条件，该摄像条件至少包含伴随如图4所示那样施加方向相同且极性相互相反的MPG脉冲的施加的2个DWI序列。此外，根据需要，将为了决定DWI数据的收集和加法运算的次数所需的第一和第二平均处理用的各AVERAGE数，设定为摄像条件。
在此，例如，设与b＝500和b＝1000对应的2个DWI数据在诊断时是必要的。在该情况下，设定与b＝0、b＝+500、b＝-500、b＝1000对应的EPI序列和多个DWI序列。
此外，如果将具有相同的b值的DWI数据作为对象的第一平均处理的AVERAGE数为2，则例如按与b＝0对应的EPI序列、与b＝+500的第一AVERAGE对应的DWI序列、与b＝-500的第一AVERAGE对应的DWI序列、与b＝1000的第一AVERAGE对应的DWI序列、与b＝+500的第二AVERAGE对应的DWI序列、与b＝-500的第二AVERAGE对应的DWI序列和与b＝1000的第二AVERAGE对应的DWI序列的顺序来设定序列的顺序。
接着，在步骤S2中，按照在摄像条件设定部40中设定的摄像条件，来执行包含DWI扫描的成像扫描。
因此，预先在床37设置被检测体P，在由静磁场电源26励磁的静磁场用磁铁21(超导磁铁)的摄像区域形成静磁场。此外，从匀场线圈电源28对匀场线圈22供给电流，对形成于摄像区域的静磁场进行均一化。
而且，从摄像条件设定部40对序列控制器31依次提供作为摄像条件而设定的多个EPI序列。于是，序列控制器31通过按照脉冲序列使倾斜磁场电源27、发射器29和接收器30驱动，从而使设置了被检测体P的摄像区域形成倾斜磁场，并且，使RF信号从RF线圈24产生。
因此，通过被检测体P的内部的核磁共振而产生的MR信号，由RF线圈24接收并提供给接收器30。接收器30从RF线圈24接受MR信号，在执行所需要的信号处理后，进行A/D(analog to digital)变换，生成作为数字数据的MR信号的原始数据。而且，接收器30将所生成的MR信号提供给序列控制器31。进而，序列控制器31将MR信号输出到计算机32的数据处理部41。
于是，数据处理部41将MR信号作为k空间数据配置于形成在k空间数据存储部42中的k空间中。其结果是，在k空间数据存储部42中，依次写入与b＝0对应的k空间数据、与b＝+500的第一AVERAGE对应的k空间数据、与b＝-500的第一AVERAGE对应的k空间数据、与b＝1000的第一AVERAGE对应的k空间数据、与b＝+500的第二AVERAGE对应的k空间数据、与b＝-500的第二AVERAGE对应的k空间数据和与b＝1000的第二AVERAGE对应的k空间数据，并保存起来。
接着在步骤S3中，图像重建部41A通过从k空间数据存储部42取入k空间数据并实施图像重建处理，从而对图像数据进行重建。其结果是，分别生成作为与b＝0对应的非DWI数据的基准图像数据、与b＝+500的第一AVERAGE和第二AVERAGE对应的2帧的量的DWI数据、与b＝-500的第一AVERAGE和第二AVERAGE对应的2帧的量的DWI数据、与b＝1000的第一AVERAGE和第二AVERAGE对应的2帧的量的DWI数据。
此外，图像重建部41A通过对将MPG脉冲的施加轴和b值设定为相同而收集的多个DWI数据进行加法运算的第一平均处理，生成AVERAGE图像数据。即，与b＝+500的第一AVERAGE和第二AVERAGE对应的2帧的量的DWI数据彼此的相加、与b＝-500的第一AVERAGE和第二AVERAGE对应的2帧的量的DWI数据彼此的相加、以及与b＝1000的第一AVERAGE和第二AVERAGE对应的2帧的量的DWI数据彼此的相加，在图像重建部41A中进行。其结果是，生成分别与b＝+500、b＝-500和b＝1000对应的3帧的量的DWI AVERAGE数据。
接着，在步骤S4中，失真校正部41B基于b值小并且符号相反的2帧的量的DWI AVERAGE数据，计算出与小的b值对应的失真校正系数。即，失真校正部41B基于与b＝±500对应的2帧的量的DWI数据，计算出与b＝500对应的失真校正系数。
失真校正系数如图5、图6或者图7所示那样，能通过将针对与二值化处理后的b＝+500对应的二值化DWI AVERAGE数据以及与二值化处理后的b＝-500对应的二值化DWI AVERAGE数据的一方或者双方的坐标变换系数作为参数的最优化处理，进行计算。
在由图5或者图6示出的方法中计算出失真校正系数的情况下，还使用与二值化处理后的b＝0对应的二值化基准图像数据。而且，在利用图5所示的方法计算出失真校正系数的情况下，在最优化处理中成为最小的对象被设定为针对二值化基准图像数据的与b＝+500和b＝-500对应的二值化DWI AVERAGE数据的各像素值的平方差的总和。
此外，在利用图6所示的方法计算出失真校正系数的情况下，在最优化处理中成为最小的对象被设定为：与b＝+500对应的二值化DWIAVERAGE数据的相对于二值化基准图像数据和与b＝-500对应的二值化DWI AVERAGE数据的各像素值的平方差的总和。
此外，在利用图7所示的方法计算出失真校正系数的情况下，在最优化处理中成为最小的对象被设定为：与b＝-500对应的二值化DWIAVERAGE数据的相对于与b＝+500对应的二值化DWI AVERAGE数据的各像素值的平方差的总和。
利用这样的最优化处理计算出的失真校正系数由于基于与失真方向相互相反的b＝±500对应的DWI数据来计算出，所以和通过使与b＝+500对应的二值化DWI AVERAGE数据的相对于二值化基准图像数据的各像素值的平方差的总和为最小的最优化处理计算出的失真校正系数相比，成为高精度。
接着，在步骤S5中，失真校正部41B基于与小的b值对应的失真校正系数，求取与大的b值对应的失真校正系数。即，失真校正部41B设为b值与失真校正系数成比例，利用式(2)，根据与b＝500对应的失真校正系数，计算出与b＝1000对应的失真校正系数。
由此，能简易地取得被设定为摄像条件的与全部b值对应的失真校正系数。按所取得的每个b值的失真校正系数能根据需要而写入失真系数存储部44中并保存起来。
接着，在步骤S6中，失真校正部41B使用与各b值对应的失真校正系数，分别进行DWI AVERAGE数据的失真校正。即，使用与b＝500和b＝1000分别对应的各失真校正系数，对与b＝500和b＝1000对应的DWI AVERAGE数据的失真进行校正。
接着，在步骤S7中，图像处理部41C针对与b＝0对应的图像数据以及失真校正后的DWI AVERAGE数据，实施图像显示用所需的图像处理。例如，能在与不同的b值对应的失真校正后的多个DWI AVERAGE数据间实施第二平均处理。在本例中，在诊断中需要与b＝500和b＝1000对应的2个DWI数据，而与b＝-500对应的DWI AVERAGE数据在诊断中不需要。因此，能针对与b＝±500对应的DWI AVERAGE数据执行第二平均处理。由此，能生成使SNR提高了的与b＝500对应的诊断用的DWI AVERAGE数据。
接着，在步骤S8中，图像处理部41C使显示装置34显示图像处理后的显示用的与b＝0对应的图像数据以及DWI AVERAGE数据。由此，用户能确认显示于显示装置34的与b＝0对应的图像、与b＝500对应的DWI和与b＝1000对应的DWI。
此时显示于显示装置34的DWI由于基于精度良好地计算出的失真校正系数来进行失真校正，所以成为良好地校正了失真的图像。此外，显示于显示装置34的DWI通过第一和第二平均处理而成为SNR提高了的图像。
接着，对利用磁共振成像装置20伴随使用了预先取得的失真校正系数的失真校正处理而进行DWI的情况进行说明。
图9是表示利用图1所示的磁共振成像装置20伴随使用预先取得的失真校正系数的失真校正处理而进行DWI时的流程的流程图。另外，在与图8所示的流程图的步骤相同的步骤中，赋予相同符号并省略说明。
首先，在预先使用了幻像的试验摄像、过去的摄像中求取的失真校正系数保存在失真系数存储部44中。其中，失真校正系数如图5、图6或者图7所示那样，通过将基于与绝对值相同且符号相互相反的b值对应的2帧的量的DWI AVERAGE数据的坐标变换系数作为参数的最优化处理，来进行计算。此外，在失真校正系数的计算中，使用与能设定为DWI的摄像条件的小的b值对应的DWI AVERAGE数据。例如，与b＝500对应的失真校正系数保存在失真系数存储部44中。
接着，在步骤S11中，摄像条件设定部40基于从输入装置33输入的信息，根据使b值为0的EPI序列以及诊断目的，来对与需要的b值对应的DWI序列进行设定。此外，根据需要设定DWI的AVERAGE数。
另外，与图8的步骤S1不同，无需对与绝对值和某个b值相同且符号相反的b值对应的DWI序列进行设定。因此，例如对与b＝0、b＝500和b＝1000对应的EPI序列进行设定。
接着，在步骤S12中，按照在摄像条件设定部40中设定的摄像条件，执行包含DWI扫描的成像扫描。即，按照与图8的步骤S1相同的流程，收集与绝对值相互不同的多个b值(b＝0、500、1000)分别对应的k空间数据。
接着，在步骤S13中，图像重建部41A生成与绝对值相互不同的多个b值(b＝0、500、1000)分别对应的图像数据以及DWI数据。因此，除了执行针对k空间数据的图像重建处理之外，还根据被设定为摄像条件的AVERAGE数，执行针对DWI数据的第一平均处理。在执行第一平均处理的情况下，生成DWI AVERAGE数据。
接着，在步骤S14中，失真校正部41B从失真系数存储部44，取得与基于b值小且符号相反的2帧的量的DWI AVERAGE数据或者DWI数据而计算出的小的b值对应的失真校正系数。即，失真校正部41B从失真系数存储部44，取得与b＝500对应的失真校正系数。
由此，在步骤S5中，失真校正部41B能基于与b＝500对应的失真校正系数来计算出与b＝1000对应的失真校正系数。而且，能与图8的步骤S6至步骤S8同样地，使显示装置34显示实施了使用失真校正系数的失真校正后的DWI所需的图像处理的图像。
换句话说，以上这样的磁共振成像装置20通过施加施加轴和b值的绝对值相同且极性相互相反的MPG脉冲来进行扩散成像，从而收集失真极性相反的2个DWI数据，基于失真极性相反的2个DWI数据，精度良好地求取失真校正系数。
因此，根据磁共振成像装置20，能比以往更精度良好地获得实施了失真校正的DWI。特别是，在被检测体P的腹部，没有仅在限定的方向扩散的神经束。因此，在设摄像部位为腹部来收集DWI的情况下，将MPG脉冲的施加轴作为1轴方向，能更精度良好地求取失真校正系数，能获得良好地校正了失真的DWI。另一方面，即使MPG脉冲的施加轴是多个方向，通过按每个施加轴来施加极性相互不同的MPG脉冲，从而也能更精度良好地计算失真校正系数。
此外，根据磁共振成像装置20，能基于与某个b值对应的失真校正系数，来计算与其他b值对应的失真校正系数。因此，能使失真校正系数的计算所需的数据处理量减少。
虽然上面对特定的实施方式进行了记载，但所记载的实施方式只不过是一个例子，并不是对发明的范围进行限定。在此所记载的新的方法和装置能通过各种其他的样式来具体化。此外，在此所记载的方法和装置的样式中，能在不从发明的要旨脱离的范围内，进行各种省略、置换和变更。一同附上的权利要求书以及其等同物，作为包含在发明的范围和要旨中的实施例，还包含这样的各种样式和变形例。

Claims (12)

1.一种磁共振成像装置，具备：
成像单元，通过伴随施加轴相同且与不同的b值对应的MPG脉冲的施加来执行扩散强调成像，从而收集与上述不同的b值对应的多个扩散强调图像数据；以及
失真校正单元，基于与上述多个扩散强调图像数据中的成为失真校正的对象的扩散强调图像数据所对应的b值不同的b值所对应的扩散强调图像数据，计算出针对成为上述失真校正的对象的扩散强调图像数据的失真校正系数，通过使用计算出的上述失真校正系数来执行成为上述失真校正的对象的扩散强调图像数据的失真校正，从而生成被实施了上述失真校正的图像数据。
2.根据权利要求1所述的磁共振成像装置，其中，
上述成像单元构成为：收集与绝对值相同且符号相反的b值对应的2个扩散强调图像数据，
上述失真校正单元构成为：基于上述2个扩散强调图像数据，计算出针对上述2个扩散强调图像数据的一方的失真校正系数。
3.根据权利要求1所述的磁共振成像装置，其中，
上述成像单元构成为：收集与绝对值不同的b值对应的多个扩散强调图像数据，
上述失真校正单元构成为：基于与成为基准的b值对应的失真校正系数，计算出与比上述成为基准的b值大的b值对应的失真校正系数。
4.根据权利要求1所述的磁共振成像装置，其中，
还具备：存储单元，存储由上述失真校正单元计算出的上述失真校正系数，
上述失真校正单元构成为：参照上述存储单元能够取得上述失真校正系数。
5.根据权利要求2所述的磁共振成像装置，其中，
上述成像单元构成为：使b值为0来进一步收集图像数据，
上述失真校正单元构成为：还基于使上述b值为0而收集到的图像数据，计算出上述失真校正系数。
6.根据权利要求2所述的磁共振成像装置，其中，
上述失真校正单元构成为：将针对上述2个扩散强调图像数据或者基于上述2个扩散强调图像数据生成的2个图像数据的各坐标变换系数作为参数，通过使上述2个扩散强调图像数据或者上述2个图像数据相对于使b值为0而收集到的图像数据的位置偏移的指标为最小的最优化处理，来计算上述各坐标变换系数，将与较小的指标对应的坐标变换系数作为上述失真校正系数。
7.根据权利要求2所述的磁共振成像装置，其中，
上述失真校正单元构成为：将针对上述2个扩散强调图像数据的一方或者基于上述2个扩散强调图像数据生成的2个图像数据的一方的坐标变换系数作为参数，基于通过最优化处理计算出的坐标变换系数，决定上述失真校正系数，其中，上述最优化处理是使上述一方相对于使b值为0而收集到的图像数据的位置偏移的指标和相对于另一方的位置偏移的指标为最小的处理。
8.根据权利要求2所述的磁共振成像装置，其中，
上述失真校正单元构成为：将针对上述2个扩散强调图像数据的一方或者基于上述2个扩散强调图像数据生成的2个图像数据的一方的坐标变换系数作为参数，基于通过使上述一方相对于另一方的位置偏移的指标为最小的最优化处理而计算出的坐标变换系数，决定上述失真校正系数。
9.根据权利要求2所述的磁共振成像装置，其中，
上述失真校正单元构成为：通过包含针对上述2个扩散强调图像数据的边缘检测在内的处理，计算出上述失真校正系数。
10.根据权利要求1所述的磁共振成像装置，其中，
还具备：图像处理单元，进行取上述失真校正前和上述失真校正后的规定的多个扩散强调图像数据的至少一方的平均的平均处理。
11.根据权利要求10所述的磁共振成像装置，其中，
上述图像处理单元构成为：在上述失真校正后对与相互不同的b值对应的多个扩散强调图像数据进行上述平均处理。
12.一种磁共振成像方法，具有如下步骤：
通过伴随施加轴相同且与不同的b值对应的MPG脉冲的施加来执行扩散强调成像，从而收集与上述不同的b值对应的多个扩散强调图像数据的步骤；以及
基于与上述多个扩散强调图像数据中的成为失真校正的对象的扩散强调图像数据所对应的b值不同的b值所对应的扩散强调图像数据，计算出针对成为上述失真校正的对象的扩散强调图像数据的失真校正系数，通过使用计算出的上述失真校正系数来执行成为上述失真校正的对象的扩散强调图像数据的失真校正，从而生成被实施了上述失真校正的图像数据的步骤。
CN201180003034.0A 2011-01-11 2011-12-19 磁共振成像装置和磁共振成像方法 Active CN102711602B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2011003509 2011-01-11
JP2011-003509 2011-01-11
PCT/JP2011/079413 WO2012096105A1 (ja) 2011-01-11 2011-12-19 磁気共鳴イメージング装置及び磁気共鳴イメージング方法

Publications (2)

Publication Number Publication Date
CN102711602A true CN102711602A (zh) 2012-10-03
CN102711602B CN102711602B (zh) 2015-07-22

Family

ID=46507017

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180003034.0A Active CN102711602B (zh) 2011-01-11 2011-12-19 磁共振成像装置和磁共振成像方法

Country Status (4)

Country Link
US (1) US8581581B2 (zh)
JP (1) JP5868160B2 (zh)
CN (1) CN102711602B (zh)
WO (1) WO2012096105A1 (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104688226A (zh) * 2013-12-06 2015-06-10 株式会社东芝 医用图像处理装置以及磁共振成像装置
CN104799859A (zh) * 2014-01-27 2015-07-29 西门子(深圳)磁共振有限公司 一种磁共振弥散加权成像图像的校正方法和装置
CN105989596A (zh) * 2015-02-12 2016-10-05 东芝医疗系统株式会社 医学图像处理装置和方法以及磁共振成像设备
CN108492307A (zh) * 2018-03-26 2018-09-04 苏州朗润医疗系统有限公司 一种磁共振adc图像分割方法及使用该方法的磁共振系统
CN112704484A (zh) * 2019-10-25 2021-04-27 通用电气精准医疗有限责任公司 磁共振成像方法及系统，非暂态计算机可读存储介质
CN113490859A (zh) * 2019-02-25 2021-10-08 皇家飞利浦有限公司 使用多幅扩散加权mr图像来计算b0图像
WO2022095142A1 (zh) * 2020-11-09 2022-05-12 无锡鸣石峻致医疗科技有限公司 一种adc-t2二维图谱的测量方法、装置、计算机设备及非均匀场磁共振系统

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6317574B2 (ja) * 2013-11-29 2018-04-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置
CN104599245B (zh) * 2014-12-23 2017-06-16 中国科学院深圳先进技术研究院 磁共振弥散张量去噪方法和系统
JP6579908B2 (ja) * 2015-10-20 2019-09-25 株式会社日立製作所 磁気共鳴イメージング装置及び拡散強調画像計算方法
KR101747029B1 (ko) 2016-04-01 2017-06-27 경희대학교 산학협력단 전류를 주입 없이 mri를 이용한 저주파 전도도 영상 복원 장치 및 그 방법
JP6752064B2 (ja) * 2016-06-23 2020-09-09 株式会社日立製作所 磁気共鳴イメージング装置、画像処理装置、及び拡散強調画像計算方法
JP6843706B2 (ja) * 2017-05-24 2021-03-17 株式会社日立製作所 磁気共鳴イメージング装置及び拡散強調画像の補正方法
JP6824132B2 (ja) * 2017-09-13 2021-02-03 ゼネラル・エレクトリック・カンパニイ 画像処理装置及び磁気共鳴イメージング装置
EP3495832A1 (en) * 2017-12-11 2019-06-12 Koninklijke Philips N.V. Diffusion mr imaging with fat suppression
DE102020212036A1 (de) * 2020-09-24 2022-03-24 Siemens Healthcare Gmbh Verfahren zur Aufnahme diffusionsgewichteter Messdaten mittels einer Magnetresonanzanlage
CN114983389B (zh) 2022-06-15 2023-01-10 浙江大学 基于磁共振扩散张量成像的人脑轴突密度定量评估方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007159718A (ja) * 2005-12-12 2007-06-28 Toshiba Corp 磁気共鳴イメージング装置および画像補正評価方法
US20100171498A1 (en) * 2009-01-02 2010-07-08 Karin Auslender Magnetic resonance method and apparatus to reduce distortions in diffusion imaging

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6969991B2 (en) * 2002-12-11 2005-11-29 The Board Of Trustees Of The Leland Stanford Junior University Correction of the effect of spatial gradient field distortions in diffusion-weighted imaging
US6815952B1 (en) * 2003-05-12 2004-11-09 The University Of Queensland Magnetic resonance diffusion imaging with eddy-current compensation
JP5181235B2 (ja) * 2008-02-25 2013-04-10 株式会社日立メディコ 画像処理装置、医用画像撮影装置、プログラムおよび画像処理方法
DE102010029932B4 (de) * 2010-06-10 2012-10-04 Siemens Aktiengesellschaft Korrektur von Bildverzerrungen in diffusionsgewichteten Magnetresonanz-Bildern

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007159718A (ja) * 2005-12-12 2007-06-28 Toshiba Corp 磁気共鳴イメージング装置および画像補正評価方法
US20100171498A1 (en) * 2009-01-02 2010-07-08 Karin Auslender Magnetic resonance method and apparatus to reduce distortions in diffusion imaging

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
H HANSSON ET AL: "A novel robust algorithm to correct for eddy current distortions in high b-value diffusion MRI", 《PROC. INT. SOC. MAGN. RESON. MED》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104688226A (zh) * 2013-12-06 2015-06-10 株式会社东芝 医用图像处理装置以及磁共振成像装置
CN104799859A (zh) * 2014-01-27 2015-07-29 西门子(深圳)磁共振有限公司 一种磁共振弥散加权成像图像的校正方法和装置
CN104799859B (zh) * 2014-01-27 2019-01-08 西门子(深圳)磁共振有限公司 一种磁共振弥散加权成像图像的校正方法和装置
US10620289B2 (en) 2014-01-27 2020-04-14 Siemens Aktiengesellschaft Correction method and apparatus for magnetic resonance diffusion weighted imaging image
CN105989596A (zh) * 2015-02-12 2016-10-05 东芝医疗系统株式会社 医学图像处理装置和方法以及磁共振成像设备
CN105989596B (zh) * 2015-02-12 2019-09-10 东芝医疗系统株式会社 医学图像处理装置和方法以及磁共振成像设备
CN108492307A (zh) * 2018-03-26 2018-09-04 苏州朗润医疗系统有限公司 一种磁共振adc图像分割方法及使用该方法的磁共振系统
CN108492307B (zh) * 2018-03-26 2022-07-05 苏州朗润医疗系统有限公司 一种磁共振adc图像分割方法及使用该方法的磁共振系统
CN113490859A (zh) * 2019-02-25 2021-10-08 皇家飞利浦有限公司 使用多幅扩散加权mr图像来计算b0图像
CN113490859B (zh) * 2019-02-25 2024-03-26 皇家飞利浦有限公司 使用多幅扩散加权mr图像来计算b0图像
CN112704484A (zh) * 2019-10-25 2021-04-27 通用电气精准医疗有限责任公司 磁共振成像方法及系统，非暂态计算机可读存储介质
WO2022095142A1 (zh) * 2020-11-09 2022-05-12 无锡鸣石峻致医疗科技有限公司 一种adc-t2二维图谱的测量方法、装置、计算机设备及非均匀场磁共振系统

Also Published As

Publication number Publication date
US20130015855A1 (en) 2013-01-17
US8581581B2 (en) 2013-11-12
WO2012096105A1 (ja) 2012-07-19
CN102711602B (zh) 2015-07-22
JP2012157687A (ja) 2012-08-23
JP5868160B2 (ja) 2016-02-24

Similar Documents

Publication Publication Date Title
CN102711602A (zh) 磁共振成像装置和磁共振成像方法
US9778336B2 (en) System and method for rapid, multi-shot segmented magnetic resonance imaging
US8648599B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN102697501B (zh) 校正磁共振图像中的相位信息的方法和相应的磁共振设备
US8032201B2 (en) Magnetic resonance imaging device
JP5366484B2 (ja) 磁気共鳴イメージング装置およびこの磁気共鳴イメージング装置における脂肪抑制効果の分析方法
US9317917B2 (en) Method, reconstruction device, and magnetic resonance apparatus for reconstructing magnetic resonance raw data
JP4619674B2 (ja) 磁気共鳴イメージング装置
US8085041B2 (en) Three-point method and system for fast and robust field mapping for EPI geometric distortion correction
WO2012020720A1 (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング方法
US10605881B2 (en) Magnetic resonance imaging apparatus and image processing method
US20120271583A1 (en) System and method for retrospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging
JP4812420B2 (ja) 磁気共鳴イメージング装置および画像補正評価方法
CN102713658B (zh) 磁化率梯度绘图
US20130207653A1 (en) Magnetic resonanace imaging apparatus and method for measuring irradiation magnetic field
JP2012040362A (ja) 磁気共鳴イメージング方法、磁気共鳴イメージング装置およびその制御装置
US7689015B2 (en) Magnetic resonance imaging apparatus and image correction estimating method
JP5666779B2 (ja) 磁気共鳴イメージング装置
US9846216B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN111093495A (zh) 磁共振成像装置、奈奎斯特重影校正方法以及奈奎斯特重影校正用程序
US11965945B2 (en) Magnetic resonance system and shimming method and imaging method thereof
Wang et al. Simultaneous pure T2 and varying T2’-weighted BOLD fMRI using Echo Planar Time-resolved Imaging (EPTI) for mapping laminar fMRI responses
US20220057467A1 (en) Epi mr imaging with distortion correction
Lee et al. MaxGIRF: image reconstruction incorporating Maxwell fields and gradient impulse response function distortion
JP4625677B2 (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: 20160725

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.