CN102415882A - 磁共振成像装置 - Google Patents

磁共振成像装置 Download PDF

Info

Publication number
CN102415882A
CN102415882A CN2011102957822A CN201110295782A CN102415882A CN 102415882 A CN102415882 A CN 102415882A CN 2011102957822 A CN2011102957822 A CN 2011102957822A CN 201110295782 A CN201110295782 A CN 201110295782A CN 102415882 A CN102415882 A CN 102415882A
Authority
CN
China
Prior art keywords
mentioned
data
processing
phase place
afi
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
CN2011102957822A
Other languages
English (en)
Other versions
CN102415882B (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 CN102415882A publication Critical patent/CN102415882A/zh
Application granted granted Critical
Publication of CN102415882B publication Critical patent/CN102415882B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • 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
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/0037—Performing a preliminary scan, e.g. a prescan for identifying a region of interest

Landscapes

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

Abstract

根据本发明的实施方式的磁共振成像装置包括成像单元和数据处理条件设定单元。成像单元从受检体收集与在k空间中在波数方向上为非对称的取样区域对应的磁共振数据，通过包含相位修正和取得复共轭的滤波处理的数据处理，基于上述磁共振数据生成图像数据。数据处理条件设定单元根据上述相位修正中使用的相位分布或对上述相位分布有影响的拍摄条件，设定上述数据处理的条件。

Description

磁共振成像装置
技术领域
本实施方式涉及磁共振成像(MRI)装置。
背景技术
MRI是用拉莫尔(Larmor)频率的高频(RF)信号磁激励放置在静磁场中的受检体的原子核自旋，根据伴随该激励产生的核磁共振(NMR：Nuclear Magnetic Resonance)信号重建图像的的拍摄法。
作为MRI中的成像法的一种，已知有AFI(非对称傅里叶成像)法。AFI法是，以在k空间中(k-space)在波数方向上为非对称的方式对数据取样，用基于被取样的自数据推定的相位分布进行相位修正后重建图像数据的方法。如果用AFI法，则可以生成与根据在k空间中被对称地取样的数据生成的图像数据同等的图像数据。
关于AFI法，提出了Margosian法、FIR(finite impulse response，有限脉冲响应)、MoFIR(改良FIR)法、POCS(projection on to convexsets，凸面集合投影)法或混合法等各种方法。另外，还已知虽然不是AIF法，但以在K空间成为非对称的方式对数据取样，在没有数据的部分填入0(0-filling，填零)，然后进行傅里叶变换(FT：Fouriertransform)的非对称数据取样法。
在Margosian法中，对被非对称取样的k空间数据施加作为窗函数的零差滤波(homodyne filter)，然后通过FT生成与非对称k空间数据对应的实空间(r-space)数据。另一方面，基于被非对称取样的k空间数据中的、被对称取样的k空间中心附近的低频区域中的k空间数据推定相位分布。然后，用被推定的相位分布进行与非对称k空间数据对应的实空间数据的相位修正。
POCS法是Margosian法的改良，是在Margosian法之后进行POCS循环处理的方法。POCS循环处理是通过反复进行实化处理、合成处理以及相位修正处理而使虚部的变化收敛到阈值以下的处理，其中实化处理指只保留相位修正后的实空间数据的实部、使虚部为0的处理，合成处理指把在使实化了的实空间数据的相位复原后进行逆傅里叶变换(IFT：inverse Fourier transform)而得到的k空间数据的非取样部与原数据的取样部合成的处理，相位修正处理指对通过合成处理后的k空间数据的FT得到的实空间数据进行的相位修正处理。如果相位修正彻底，则POCS法基于使实空间数据的虚部为0的原理。如果用POCS法，则通过反复进行数次POCS循环处理，可以降低在Margosian法中因零差滤波处理生成的相位修正的错误。
另一方面，在FIR法中，在对被非对称取样的k空间数据施加零差滤波之前进行相位修正。即，在FIR法中，在通过非对称k空间数据的FT生成的实空间数据的相位修正后，通过IFT把相位修正后的实空间数据变换成k空间数据。然后，对相位修正后的k空间数据施加零差滤波。在该FIR法中，由于数据处理时间比Margosian法长出两次FT的大小，但在零差滤波处理前进行相位修正，所以可以降低因零差滤波处理生成的相位修正的错误。
MoFIR法是FIR法的改良法，以使得基于不仅仅是k空间数据被对称取样的低频区域、还包含被非对称取样的部分的全部k空间数据推定相位修正用的相位分布。即，在FIR法中，仅仅根据k空间中心附近的低频区域中的k空间数据推定相位修正用的低频区域的相位分布，而在MoFIR法中，基于被非对称取样的全部k空间数据推定相位修正用的相位分布。因此，在MoFIR法中，与FIR法相比，虽然与本来的相位分布不同，但可以推定更高频区域的相位分布。因此，在MoFIR法中，可以降低在Margosian法、FIR法中的因零差滤波处理生成的相位修正的错误。
另一方面，在作为最简单的针对非对称取样数据的重建法的填零中，图像产生模糊。但是，在像被对称取样的数据的70％以上被取样那样的、非对称度比较小的非对称取样的情况下，图像的模糊在允许范围内。另外，在进行填零时，不需要特别的处理，而且AFI中的相位的过修正造成的伪影(artifact)也不会出现。因此，在非对称度比较小的情况下，仍然大量使用填零。
混合法是把填零和AFI法组合起来的方法。即，混合法是把由填零得到的填零图像与由AFI得到的AFI图像加权并相加合成的方法。更具体地说，调整权重以使得在填零图像与AFI图像之间，相位差或振幅差大的部分是填零图像，相位差或振幅差小的部分是AFI图像。
在现有的AFI中，如果相位推定不好，相位误差大，则存在有时与进行填零时相比，伪影显著出现的问题。因此，希望把为了图像重建所需的数据处理时间的增加抑制到最小限度，且以与根据在k空间中被对称取样的数据生成的图像同等的精度，根据被非对称取样的数据生成图像数据。
发明内容
本发明的目的在于提供可以抑制数据处理时间的增加、且基于在k空间中被非对称取样的MR数据生成更高精度的图像数据的磁共振成像装置。
根据本发明的实施方式的磁共振成像装置，包括成像单元和数据处理条件设定单元。成像单元从受检体收集与在k空间中在波数方向上为非对称的取样区域对应的磁共振数据，通过包含相位修正和取得复共轭的滤波处理的数据处理，生成基于上述磁共振数据的图像数据。数据处理条件设定单元根据上述相位修正中使用的相位分布或对上述相位分布有影响的拍摄条件，设定上述数据处理的条件。
附图说明
图1是示出根据本发明的实施方式的磁共振成像装置的构成图。
图2是示出图1所示的RF线圈的详细构成的一例的图。
图3是图1所示的计算机的功能框图。
图4是说明通过预扫描确定AFI中的数据取样区域的边界的方法的图。
图5是示出图3所示的数据处理部中执行的第一AFI处理的例子的图。
图6是示出图3所示的数据处理部中执行的第二AFI处理的例子的图。
图7是示出图3所示的数据处理部中执行的第三AFI处理的例子的图。
图8是示出通过一维模拟由第一AFI处理生成的实空间数据I(x)的例子的图。
图9是示出通过一维模拟由第二AFI处理生成的实空间数据I(x)的例子的图。
图10是示出通过一维模拟由第三AFI处理生成的实空间数据I(x)的例子的图。
图11是示出通过一维模拟对被非对称取样的k空间进行填零生成的实空间数据I0fill(x)的例子的图。
图12是示出在图8～图12的模拟中使用了的部分的原k空间数据Sorig(k)的图。
图13是示出在图12所示的部分的原k空间数据Sorig(k)的制作中使用的全部频率区域中的全部k空间数据Sfull(k)的图。
图14是示出把图13所示的全部k空间数据Sfull(k)变换得到的实空间数据I(x)和相位分布Φ(x)的图。
图15是示出与非对称取样区域对应的零差滤波的增益Hhomo.max＝2的第三AFI处理的模拟结果的图。
图16是示出与非对称取样区域对应的零差滤波的增益Hhomo.max＝4的第三AFI处理的模拟结果的图。
图17是示出在第三AFI处理中改变了与非对称取样区域对应的零差滤波的最大增益Hhomo.max时RMSE的变化的图。
图18是说明确定零差滤波的强度分布形状的参数的图。
图19是说明确定零差滤波的最大增益Hhomo.max以使得基于通过填零生成的数据的相位修正错误的指标值最小的方法的图。
图20是示出相位修正后的k空间数据SPcor(k)的信号的绝对值的图。
图21是示出图20所示的相位修正后的k空间数据SPcor(k)的生成中使用的原始数据的图。
图22是示出把零差滤波的与非对称取样部分对应的增益设定成向着高频部分衰减至1的例子的图。
图23是示出AFI处理中使用的滤波的具体形状的例子的图。
图24是示出根据AFI处理的种类可变地设定AFI处理中使用的零差滤波的最大增益Hhomo.max(k)的例子的图。
图25是示出利用图1所示的磁共振成像装置20通过AFI和PI收集诊断图像时的流程的流程图。
具体实施方式
参照附图，说明根据本发明的实施方式的磁共振成像装置。
图1是示出根据本发明的实施方式的磁共振成像装置的构成图。
磁共振成像装置20包括：形成静磁场的筒状的静磁场用磁体21、在该静磁场用磁体21的内部设置的匀磁(Shim)线圈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作为拍摄区域，受检体P置于机台37上。在RF线圈24中具有：内置在机架中的接收发送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信号造成的激励产生的NMR信号并赋予接收器30的功能。
图2是示出图1所示的RF线圈24的详细构成的一例的图。
像图2所示的那样，RF线圈24包括筒状的全身用(WB)线圈24a和相位阵列(phased array)线圈24b。相位阵列线圈24b包括多个线圈要素24c，在被检测体P的体表侧和背面侧分别配置多个线圈要素24c。
另一方面，接收器30包括：收发转换器(duplexer)30a、放大器30b、切换合成器30c和接收系统电路30d。收发转换器30a与发送器29、WB线圈24a和WB线圈24a用的放大器30b连接。放大器30b按各线圈要素24c和WB线圈24a的数目设置，分别单独与各线圈要素24c和WB线圈24a连接。设置一个或多个切换合成器30c，切换合成器30c的输入侧经由多个放大器30b与多个线圈要素24c或WB线圈24a连接。接收系统电路30d按所希望的数目设置，以使其数目为各线圈要素24c和WB线圈24a的数目以下，且设置在切换合成器30c的输出侧。
WB线圈24a可以用作RF信号的发送用的线圈。另外，作为NMR信号的接收用的线圈可以使用各线圈要素24c。而且，WB线圈24a也可以用作接收用的线圈。
因此，收发转换器30a构成为，把从发送器29输出的发送用的RF信号赋予WB线圈24a，把在WB线圈24a中接收的NMR信号经由接收器30内的放大器24d赋予切换合成器30c。另外，还构成为，在各线圈要素24c中接收的NMR信号也分别经由对应的放大器24d输出到切换合成器30c。
构成为，切换合成器30c进行从线圈要素24c、WB线圈24a接收的NMR信号的合成处理和切换，向对应的接收系统电路30d输出。换言之，构成为，可以按照接收系统电路30d的数目在切换合成器30c中进行从线圈要素24c、WB线圈24a接收的NMR信号的合成处理和切换，用所希望的多个线圈要素24c形成与拍摄部位对应的灵敏度分布，接收来自各拍摄部位的NMR信号。
但是，也可以不设置线圈要素24c，只用WB线圈24a接收NMR信号。另外，也可以不设置切换合成器30c，直接把在线圈要素24c、WB线圈24a中接收的NMR信号输出到接收系统电路30d。而且，也可以在整个宽范围上配置更多的线圈要素24c。
另一方面，控制系统25的序列控制器31与倾斜磁场电源27、发送器29和接收器30连接。序列控制器31具有存储记述了为了驱动倾斜磁场电源27、发送器29和接收器30所需的控制信息，例如应向倾斜磁场电源27施加的脉冲电流的强度、施加时间、施加定时等的动作控制信息的序列信息的功能；以及通过根据所存储的预定的序列驱动倾斜磁场电源27、发送器29和接收器30产生X轴倾斜磁场Gx、Y轴倾斜磁场Gy、Z轴倾斜磁场Gz和RF信号的功能。
另外，序列控制器31构成为，接收通过接收器30中的NMR信号的检波和A/D(模/数)变换得到的复数数据即原始数据(raw data)并赋予计算机32。
因此，发送器29具有基于从序列控制器31接收的控制信息向RF线圈24赋予RF信号的功能，而接收器30具有通过对从RF线圈24接收的NMR信号进行检波、实施所需的信号处理并进行A/D变换，生成数字化了的复数数据即原始数据的功能和把生成的原始数据赋予序列控制器31的功能。
而且，磁共振成像装置20中，具有取得受检体P的心电图(ECG：electro cardiogram)信号的ECG单元38。构成为，由ECG单元38取得的ECG信号经由序列控制器31输出到计算机32。另外，也可以取代把心跳表示为心搏信息的ECG信号，而取得把心跳表示为脉搏信息的脉搏同步(PPG：peripheral pulse gating)信号。
另外，通过用运算装置35执行在计算机32的存储装置36中保存的程序，计算机32具有各种功能。但是，也可以取代程序的至少一部分，而在磁共振成像装置20中设置具有各种功能的特定电路。
图3是图1所示的计算机32的功能框图。
计算机32的运算装置35通过执行在存储装置36中保存的程序，而用作拍摄条件设定部40、数据处理部41和AFI处理条件设定部42。另外，存储装置36用作k空间数据库43、图像数据库44、相位数据库45和AFI条件数据库46。数据处理部41具有PI(平行成像)处理部41A和AFI处理部41B。
拍摄条件设定部40具有基于来自输入装置33的指示信息设定包含脉冲序列的拍摄条件，通过把设定的拍摄条件输出到序列控制器31而驱动控制序列控制器31的功能。
特别地，拍摄条件设定部40构成为可以设定在k空间中在波数方向上为非对称的方式对MR数据取样的AFI的拍摄条件和用多个线圈要素24c接收数据的PI的拍摄条件。AFI和PI可以在包含磁共振血管摄影法(MRA：magnetic resonance angiography)、扩散加强成像(DWI：diffusion weighted imaging)的各种成像方式中。然后，在AFI中，从受检体P收集与k空间中在波数方向上非对称的取样区域对应的MR数据。
另外，在PI中有SENSE(Sensitivity Encoding，灵敏度编码)、SMASH(Shimadzu minimum angle shot，岛津最小角度拍摄)、GRAPPA(generalized autocalibrating partially parallel acquisition，统一自动标准部分平行取得)等种类。在SENSE中在实空间中执行数据处理，而在SMASH中在k空间中处理数据。GRAPPA是SMASH的发展，所以是在针对每个线圈要素24c的k空间数据执行作为PI特有的后处理的展开处理和图像重建处理之后，把展开处理后的每个线圈要素24c的图像数据合成的方法。
作为用于进行AFI的拍摄条件，如果是二维(2D)取样，则有k空间数据的读出方向和相位编码方向中的一个方向k的波数方向上的取样区域(-Kc≤k≤Kmax)。取样区域的边界Kc，例如，可以根据来自AFI处理条件设定部42的设定信息可变地设定。非取样区域可以是k方向的正侧和负侧中的任一侧，但在此以把负侧作为非取样区域的情形为例进行说明。
或者，也可以基于通过预扫描收集的数据确定取样区域的边界Kc。此时，在拍摄条件设定部40中设定预扫描用的数据收集条件。具体地，把使取样区域的边界Kc缓慢变化而进行MR数据收集的条件设定为预扫描用的数据收集条件。
图4是说明通过预扫描确定AFI中的数据取样区域的边界Kc的方法的图。
像图4(A)所示的那样，将使取样区域的边界Kc缓慢变化成Kc1、Kc2、Kc3......而进行数据收集的条件设定为用于预扫描。如果把用于求该边界Kc的预扫描称为Kc-PREP扫描，则Kc-PREP扫描的数据收集条件尽可能与成像扫描的数据收集条件一致有助于求出更合适的边界Kc。另一方面，也可以为了缩短数据收集时间而使相位编码为0，即作为一维(1D)的投影数据改变数据取样区域的边界Kc。
如果针对通过Kc-PREP扫描收集的非对称取样数据进行后述的AFI处理，则像图4(B)所示的那样，在彼此不同的多个取样区域的边界Kc1、Kc2、Kc3......分别生成对应的多个2D图像数据I(Kc1)、I(Kc2)、I(Kc3)......。这样的话，使用者通过进行通过输入装置33输入的图像选择信息或阈值处理等的图像处理，可以选择合适图像质量的图像I(Kc_opt)。然后，可以基于来自这多个图像数据的图像数据的选择信息，把与合适图像质量的图像I(Kc_opt)对应的数据取样区域的边界Kc_opt确定为成像扫描用的拍摄条件和AFI数据处理条件。
然后，像图4(C)所示的那样，把被选择的数据取样区域的边界Kc_opt作为拍摄条件执行成像扫描。
进而，拍摄条件设定部40具有设定预扫描用的数据收集条件的功能，预扫描用的数据收集条件用来设定其它的成像扫描用的拍摄条件、取得用成像扫描收集的数据的处理所必需的数据。例如，在进行静磁场的匀磁时，设定匀磁用的预扫描的数据收集条件。根据通过匀磁用的预扫描收集的MR数据测定用来确定应从匀磁线圈电源28向匀磁线圈22供给的电流的磁场分布图。作为另一个例子，有用来测定针对通过PI收集的数据的展开处理所必需的线圈要素24c的灵敏度图的预扫描。
数据处理部41具有：从序列控制器31接收利用AFI和PI用的条件收集的原始数据并作为k空间数据配置到在k空间数据库43中形成的k空间中的功能；以及通过执行针对k空间数据的AFI处理和由PI造成的PI处理而生成诊断用的图像数据的功能。AFI处理是，基于在k空间中在波数方向上被非对称取样的k空间数据，伴随着相位修正处理生成与根据被对称取样的数据生成的图像同等的图像数据的处理。即，AFI处理是通过包含相位修正和取得复共轭的滤波处理的数据处理，生成基于被对称取样的MR数据的图像数据的处理。通过PI处理和AFI处理生成的诊断用的图像数据可以显示在显示装置34上。
数据处理部41的PI处理部41A具有执行PI处理的功能，AFI处理部41B具有执行AFI处理的功能。在PI处理和AFI处理中，有时k空间数据与实空间数据之间的变换是必需的。于是在PI处理部41A和AFI处理部41B中具有设置空间数据与实空间数据之间的变换处理功能，并且进行从k空间数据库43读入k空间数据、向k空间数据库43写入k空间数据、从图像数据库44读入实空间数据和向图像数据库44写入实空间数据的功能。
因此，在k空间数据库43中保存k空间数据，在图像数据库44中保存实空间数据。
在PI处理中，有通过把发生折回的图像数据展开而除去折回的展开处理和把每个线圈要素24c的图像数据合成的合成处理。另一方面，在AFI处理中，可举出利用了复共轭对称性的数据填补用的k空间数据的滤波处理、实空间数据的相位修正处理、相位修正错误的减少处理等的处理。
图5是示出图3所示的数据处理部41中执行的第一AFI处理的例子的图。
图5示出对在k空间的kx方向上非对称的k空间数据进行取样的例子，但在对在ky方向上非对称的k空间数据进行取样时也是一样的。另外，也可以在2D方向上对三维(3D)k空间数据非对称性地取样，进行AFI处理。下面，为了简化说明，以在k方向进行1D的非对称取样的情形为例进行说明。
另外，如果把对数据非对称取样的方向作为读出方向，则可以缩短回波时间(TE：Echo time)。另一方面，如果把对数据非对称取样的方向作为编码方向，则可以缩短拍摄时间。
像图5所示的那样，在AFI法中，在k空间中以关于中心呈非对称的方式收集MR数据，作为k空间数据配置在k空间中。因此，在k空间中，形成关于k＝0非对称的取样区域，把一部分k空间数据填充到非对称取样区域中。图5示出以1D方式把除去了一方的高频侧的数据的位于横轴方向的-Kc≤k≤Kmax范围内的k空间数据作为部分的原k空间数据Sorig(k)而配置的例子。
在第一AFI处理中，首先，像式(1)所示的那样，对部分的原k空间数据Sorig(k)(-Kc≤k≤Kmax)施加零差滤波Hhomo(k)后，通过FT生成原实空间数据Ihomo(x)。
零差滤波是与向k空间的无数据部分填充复共轭数据等价的滤波。即，通过利用零差滤波进行的部分的原k空间数据Sorig(k)的窗化，执行利用了复共轭对称的k空间数据的填补处理。零差滤波是k空间数据的对称取样部分与非对称取样部分权重不同的窗函数。
Ihomo(x)＝FT{Hhomo(k)Sorig(k)}  (1)
另一方面，像式(2)所示的那样，在部分的原k空间数据Sorig(k)上施加LPF(低通滤波)Hlow(k)，抽出在k空间中对称的低频区域的原k空间数据Slow(k)(-Kc≤k≤Kc)。然后，通过被抽出的低频区域的原k空间数据的FT，生成与低频数据对应的低通滤波了的实空间数据Ilow(x)。
Ilow(x)＝FT{Hlow(k)Sorig(k)}  (2)
然后，像式(3-1)和式(3-2)所示的那样，根据低通滤波了的实空间数据Ilow(x)求低频区域中的相位分布Φlow(x)，用低频区域中的相位分布Φlow(x)执行原实空间数据Ihomo(x)的相位修正后，进行除去虚数部分的实化处理。由此，生成AFI处理后的实空间数据Icor(x)。
exp{-iΦlow(x)}＝Ilow *(x)/abs{Ilow(x)}        (3-1)
Icor(x)＝Re[Ihomo(x)exp{-iΦlow(x)}]          (3-2)
其中，式(3-2)中Re[]是输出复数的实部的函数。
即，AFI是利用复共轭的对称性埋入非取样部分的k空间数据的成像法，但实际上不能无视相位分布，由于虚部不为0，所以对称性不成立。于是，可以通过基于相位分布的相位修正确保复共轭的对称性。但是，即使利用该相位修正有时也产生修正错误。
于是，针对相位修正后的实空间数据Icor(x)执行预定次数的用来减少相位修正的错误的循环处理。该循环处理基于如果相位修正中没有错误(误差)则为了使k空间中的复共轭为对称而使相位修正后的实空间数据Icor(x)的虚部为0的原理。因此，循环处理是使实空间数据Icor(x)的虚部收敛到0的收敛计算。然后，通过反复进行循环处理的收敛计算，减少零差滤波处理造成的相位修正的错误。
具体地，像式(4)所示的那样，把相位修正后的只有实部的实空间数据Icor(x)作为循环处理中的初始的实空间数据Ij(x)，执行使只有实部的实空间数据Ij(x)的相位复原到相位修正前的相位的处理。由此生成相位复原了的实空间数据I′j(x)。
I’j(x)＝Ij(x)exp{iΦlow(x)}   (4)
然后，像式(5)所示的那样，通过IFT把相位复原处理后的相位复原了的实空间数据I′j(x)变换成k空间数据Sj(k)。
Sj(k)＝IFT{I’j(x)}   (5)
然后，像式(6)所示的那样，把k空间数据Sj(k)的一部分即-Kc≤k≤Kmax的范围置换成利用数据切出滤波Hmerge(k)从部分的原k空间数据Sorig(k)切出的k空间数据。换言之，把在式(5)中得到的k空间数据Sj(k)中的k＜-Kc范围的k空间数据与部分的原k空间数据Sorig(k)中-Kc≤k≤Kmax范围的k空间数据合成。
Sj+1(k)＝Hmerge(k)Sorig(k)+{1-Hmerge(k)}Sj(k)   (6)
然后，针对通过在式(6)中得到的合成后的k空间数据Sj+1(k)的FT得到的实空间数据进行相位修正后，进行除去虚数部分的实化处理。即，通过把置换和合成后的k空间数据Sj+1(k)变换成实空间数据并进行相位修正的处理、以及取得相位修正后的实空间数据的实部的处理，像式(7)所示的那样，生成再次相位修正后的只有实部的实空间数据Ij+1(x)。
Ij+1(x)＝Re[FT{Sj+1(k)}exp{-iΦlow(x)}]  (7)
使j＝0，1，2，3......反复进行式(4)～式(7)的循环处理，直到满足预定的收敛条件。即，执行使作为取得实部的处理的对象的实空间数据的虚部收敛到0的收敛计算。收敛条件可以是例如，像式(8)所示的那样，第j个实空间数据Ij(x)与第j+1个实空间数据Ij+1(x)的差小于阈值Eth的情形。
|Ij+1(x)-Ij(x)|＜Eth   (8)
然后，通过反复进行循环处理直到满足收敛条件，生成作为诊断图像数据的实空间数据I(x)。因此，在第一AFI处理中，如果反复进行循环处理的次数为j次，则FT的次数为2+j次，IFT的次数为j次。因此，第一AFI处理中的数据处理时间包含2+j次的FT和j次的IFT所需的时间。
图6是示出图3所示的数据处理部41中执行的第二AFI处理的例子的图。
在图6所示的第二AFI处理中，与第一AFI处理的不同之处主要在于，在对部分的原k空间数据Sorig(k)施加零差滤波Hhomo(k)之前进行相位修正处理。因此，针对与第一AFI处理同样的处理省略说明。
在第二AFI处理中也是，与第一AFI处理同样地，利用式(2)根据部分的原k空间数据Sorig(k)生成与低频数据对应的低通滤波了的实空间数据Ilow(x)。
另一方面，像式(9)所示的那样，利用切出全部频率区域的数据的切出滤波Hwhole，从部分的原k空间数据Sorig(k)切出非对称的原k空间数据，通过FT生成原实空间数据Iwhole(x)。
Iwhole(x)＝FT{HwholeSorig(k)}  (9)
然后，像式(10)所示的那样，用从低通滤波了的实空间数据Ilow(x)求得的低频区域中的相位分布Φlow(x)执行原实空间数据Iwhole(x)的相位修正。由此，生成相位修正后的实空间数据IPcor(x)。
Ipcor(x)＝Iwhole(x)exp{-iΦlow(x)}]   (10)
然后，像式(11)所示的那样，通过IFT把相位修正后的实空间数据IPcor(x)变换成相位修正后的k空间数据SPcor(k)。
SPcor(k)＝IFF{Ipcor(x)}  (11)
然后，像式(12)所示的那样，对相位修正后的k空间数据SPcor(k)施加零差滤波Hhomo(k)后，通过FT变换成实空间数据。进一步进行除去实空间数据的虚数部分的实化处理。由此，生成零差滤波处理和相位修正后的实空间数据Icor(x)。
Icor(x)＝Re[FT{Hhomo(k)SPcor(k)}]    (12)
对这样生成的相位修正后的实空间数据Icor(x)以预定次数进行与第一AFI处理同样的循环处理。因此，在第二AFI处理中，如果反复进行循环处理的次数为j次，则FT的次数为3+j次，IFT的次数为1+j次。因此，第二AFI处理中的数据处理时间比第一AFI处理长相当于FT和IFT各自一次的时间。但是，在对部分的原k空间数据Sorig(k)施加零差滤波Hhomo(k)之前进行相位修正处理。因此，在第二AFI处理中，与第一AFI处理相比可以更加减少零差滤波处理造成的相位修正的错误。
图7是示出图3所示的数据处理部41中执行的第三AFI处理的例子的图。
在图7所示的第三AFI处理中，与第二AFI处理的不同之处主要在于，基于被非对称取样的部分的原k空间数据Sorig(k)的全部推定全部频率区域的相位分布Φwhole(x)，用全部频率区域中的相位分布Φwhole(x)执行相位修正。即，通过不仅使用被对称取样的部分，还使用被非对称取样的部分的部分的原k空间数据Sorig(k)，推定全部频率区域的相位分布Φwhole(x)用于相位修正。关于其它部分由于与第二AFI处理相同，所以省略详细说明。
在第三AFI处理中，在相位修正中使用全部频率区域中的相位分布Φwhole(x)。因此，在第三AFI处理中，与第二AFI处理相比可以更加减少相位修正的错误。另外，在第三AFI处理中，FT执行3+j次，IFT执行1+j次。
在上述的AFI处理方法中，从作为修正对象的k空间数据自身推定相位修正用的相位分布，但也可以事先根据除了在诊断图像数据的生成中使用的MR数据以外的其它MR数据求得相位修正用的相位分布。
特别地，如果不用AFI法进行数据收集，则在诊断图像数据的生成中得到比被非对称取样的k空间数据更宽的频率区域中的k空间数据。例如，如果收集全部频率区域的k空间数据，则可以求全部频率区域的相位分布Φwhole(x)。因此，可以用根据与利用AFI法生成图像数据时使用的MR数据不同、且在比MR数据更宽的频率区域中收集的其它MR数据求得的相位分布进行相位修正。
在不从自数据而从其它数据求相位分布时，根据同一受检体的同一拍摄部位收集的数据求相位分布在精度上是实用的。即，可以推定更接近本来的相位分布的相位分布。此时，重要的是，在至少比自数据高的频率区域中对其它数据取样。
例如，通过用来测定静磁场的匀磁用的磁场分布图的匀磁序列，得到被大致完全取样的k空间数据。另外，为了测定PI处理用的线圈要素24c的灵敏度分布，也通过预扫描收集比自数据宽的k空间(频率区域)的数据。因此，可以利用通过匀磁序列、灵敏度分布图取得用序列收集到的数据，求全部频率区域中的相位修正用的相位分布。即，可以利用根据通过用来测定静磁场的匀磁用的磁场分布图的预扫描或用来测定多个线圈要素24c的灵敏度分布的预扫描收集的MR数据求得的相位分布，进行相位修正。
于是，AFI处理部41B中具有基于通过预扫描等的所希望的扫描收集的其它数据求相位修正用的相位分布的功能、把求得的相位分布写入相位数据库45的功能、和从相位数据库45读入相位分布的功能。因此，在相位数据库45中保存事先取得的AFI处理中的相位修正用的相位分布。
根据其它数据求相位分布时，也可以用其它数据和自数据这两者。例如，如果表示与序列对应的其它数据的相位偏移与作为相位修正对象的数据的相位偏移的类似度的指标为阈值以上，则把根据其它数据求得的全部频率区域中的相位分布作为相位修正用的相位分布。而如果表示相位偏移间的类似度的指标小于阈值，则通过k空间中的窗化抽出被大致完全取样的其它k空间数据中的、作为相位修正对象的自数据中未被取样的部分的数据。然后，把抽出的k空间数据变换成实空间数据，把基于得到的实空间数据算出的相位分布与基于自数据算出的相位分布合成。
由此，可以求基于其它数据和自数据这两者的全部频率区域中的相位分布图。通过这样的相位分布的计算方法，虽然严格地讲在相位修正中使用与本来的相位分布不同的相位分布，但可以减少相位修正错误的发生。
另外，在上述的例子中，对k空间数据执行相位修正，但也可以对实空间数据执行相位修正。
另外，PI处理和AFI处理可以通过优化处理顺序缩短处理时间。为此，PI处理部41A和AFI处理部41B构成为分别把处理中或处理后的中间数据赋予另一方。例如，如果对与多个线圈要素24c对应的各k空间数据执行了AFI处理中的零差滤波处理后，进行PI处理中的展开处理和图像数据合成处理，则用来减少相位修正错误的循环处理的对象是单一的实空间数据。换言之，如果在PI处理中的展开处理和图像数据合成处理后进行循环处理，则循环处理的对象只是在线圈要素24c间合成了的图像数据。因此，可以减少循环处理的次数，减少全部的数据处理量和数据处理时间。
AFI处理条件设定部42具有根据相位修正处理中使用的相位分布或对相位分布有影响的拍摄条件把AFI处理的条件设定成合适的条件的功能。作为根据相位修正中使用的相位分布或对相位分布有影响的拍摄条件设定的AFI处理的条件，像上述那样，有AFI处理的种类、循环处理的次数j、取样区域的边界Kc、以及零差滤波的强度和形状等。零差滤波的强度是被非对称取样的MR信号的增益。
另外，通过使取样区域的边界Kc缓慢变化而进行MR数据收集的Kc-PREP扫描生成的多个2D图像数据，是具有与相位修正中使用的相位分布或对相位分布有影响的拍摄条件对应的图像质量的图像数据。因此，通过选择具有合适的图像质量的2D图像数据来设定取样区域的边界Kc，相当于设定与相位修正中使用的相位分布或对相位分布有影响的拍摄条件对应的合适的取样区域的边界Kc。
通过例如把三种变更方法参数化，可以确定AFI处理的种类。第一种变更方法是在零差滤波处理之前进行相位修正处理还是在零差滤波处理之后进行相位修正处理的变更方法。第二种变更方法是在相位修正中使用与数据处理范围对应的全部频率区域的相位分布Φwhole(x)、还是在相位修正中使用与非对称取样区域中的对称部分对应的低频区域的相位分布Φlow(x)的变更方法。第三种变更方法是仅从作为修正对象的k空间数据自身即自数据推定相位修正中使用的相位分布Φ(x)、还是仅从其它数据推定，还是从其它数据和自数据这两者推定的变更方法。
通过利用输入装置33的操作向AFI处理条件设定部42输入设定信息，可以根据用于相位修正的相位分布或对相位分布产生影响的拍摄条件可变地设定作为AFI用的数据处理条件的各种参数。另外，也可以预先准备与相位分布或拍摄条件对应的多个AFI处理的条件作为候补，从而可以利用输入装置33的操作进行选择。而且，也可以把与各种拍摄条件对应的多个拍摄规程分别择一地与合适的AFI处理条件的参数相关联而保存起来。
在AFI条件数据库46中保存预先通过模拟或凭经验确定的、与多个拍摄条件对应的AFI用的数据处理条件或与拍摄条件对应的AFI用的数据处理条件的候补。然后，AFI处理条件设定部42构成为，在从拍摄条件设定部40取得拍摄规程的选择信息等的拍摄条件信息时，从作为存储装置的AFI条件数据库46取得与对相位分布有影响的拍摄条件对应的AFI用的数据处理条件或AFI用的数据处理条件的候补。此时，由于如果利用例如输入装置33的操作选择拍摄规程，则使用者无意识地自动设定对应的AFI处理条件的参数，所以有助于提高使用者的便利性。
但是，关于取样区域的边界Kc，也可以像上述那样，通过执行Kc-PREP扫描来确定。
作为对相位分布有影响的拍摄条件，可以举出成像序列的种类、TE、拍摄部位等的条件。例如，在用FASE(fast advanced spin echo或fast asymmetric spin echo，快速高级自旋回波或快速非对称自旋回波)等的FE(field echo，场回波)系的序列、TE长的GE(gradientecho，梯度回波)系的序列、EPI(echo planar imaging，回波平面成像)序列、稳态自由进动(SSFP：steady state free precession)序列进行成像时，有时因相位修正的误差而使伪影显著出现。
另外，在进行拍摄部位是血管的MRA、DWI时、磁共振成像装置20是强磁场装置时，要求高精度的相位修正。在MRA中有FS-BB(flow-sensitive black-blood，流动敏感性黑血)法、TOF(time offlight，飞行时间)法等。FS-BB法是在把血液作为低信号区域而描绘成黑色的BB(black blood，黑血)法中，通过作为去相(dephase)倾斜磁场脉冲施加MPG(motion probing gradient，移动探测梯度)脉冲，选择性地降低来自关注区域中的血液的信号的方法。因此，在FS-BB法中，血液的相位积极地分散，即被去相化。
因此，在这些要求高精度的相位修正的条件下进行成像时，如果相位修正存在误差就会使伪影显著出现。由此，通过设定针对每个拍摄条件具体化了的合适的AFI处理条件，可以提高减少伪影的效果。
下面，示出通过模拟确定合适的AFI处理条件的参数的例子。首先，示出确定AFI处理的种类的例子。
图8是通过一维模拟示出由第一AFI处理生成的实空间数据I(x)的例子的图；图9是通过一维模拟示出由第二AFI处理生成的实空间数据I(x)的例子的图；图10是通过一维模拟示出由第三AFI处理生成的实空间数据I(x)的例子的图；图11是通过一维模拟示出对被非对称取样的k空间进行填零生成的实空间数据I0fill(x)的例子的图；图12是示出在图8～图12的模拟中使用了的部分的原k空间数据Sorig(k)的图；图13是示出在图12所示的部分的原k空间数据Sorig(k)的制作中使用的全部频率区域中的全部k空间数据Sfull(k)的图；图14是示出把图13所示的全部k空间数据Sfull(k)变换得到的实空间数据I(x)和相位分布Φ(x)的图。
在图8～图11中，各横轴表示一维的像素位置x，各纵轴表示像素位置x处的实空间数据I(x)。另外，在图8～图11中，(A)示出不进行循环处理时即j＝0时的实空间数据I0(x)，(B)示出进行了4次循环处理时即j＝4时的实空间数据I4(x)。
另外，在图12中，横轴表示k空间的位置k，纵轴表示位置k处的部分的原k空间数据Sorig(k)。像图12所示的那样，通过图5、图6和图7分别示出的第一、第二和第三AFI处理，根据Kc＝16的部分的原k空间数据Sorig(k)(-16≤k≤128)生成了像图8、图9和图10分别所示的那样的实空间数据I0(x)、I4(x)。另外，为了比较参考，示出进行图11的填零生成的实空间数据I0fill(x)。
另外，图12所示的部分的原k空间数据Sorig(k)是从图13所示的全部频率区域中的全部k空间数据Sfull(k)切出非对称取样区域-16≤k≤128的数据并移动5个相位得到的数据。图14(A)示出把图13所示的全部频率区域中的全部k空间数据Sfull(k)FT得到的实空间数据I(x)，图14(B)示出从图14(A)所示的实空间数据I(x)得到的相位分布Φ(x)。
根据图8～图11所示的模拟结果可知，通过第二AFI处理(即在零差滤波处理之前进行使用了低频区域的相位分布Φlow(x)的相位修正处理)的方法和第三AFI处理(即在零差滤波处理之前进行使用了全部频率区域的相位分布Φwhole(x)的相位修正处理)的方法，可以减少相位修正错误。
图5、图6和图7所示的减少相位修正错误的循环处理中，循环处理前的相位修正错误越少则收敛越快。因此，从AFI处理的高速化的角度可以看出，第二AFI处理或第三AFI处理是合适的。尤其是，在部分的原k空间数据Sorig(k)的峰值从k空间的中心偏移时等相位误差大时，优选执行在零差滤波处理之前进行相位修正处理的第二AFI处理或第三AFI处理。
另外，可以看出，不进行循环处理时，使用全部频率区域的相位分布Φwhole(x)的第三AFI处理的相位修正错误最少。因此，在既要缩短数据处理时间又要确保图像质量时，第三AFI处理是合适的。
另外，可以看出，在第三AFI处理中，如果执行循环处理则相位修正错误改善。因此，在重视进一步提高图像质量时，在第三AFI处理中执行循环处理是合适的处理条件。但是，通过进行模拟确认了，虽然用循环处理改善了一次以下的低次的相位修正错误，但即使反复进行4次以上的循环处理也没有显著的改善。
另一方面，可以看出，即使在使用低频区域的相位分布Φlow(x)的第二AFI处理中执行循环处理，也几乎没有改善相位修正错误。因此，在得不到全部频率区域的相位分布Φwhole(x)时，不执行循环处理从而不会增加数据处理时间的第二AFI处理是合适的。
相反，可以看出，在零差滤波处理之后进行相位修正处理的第一AFI处理的情况下，虽然如果有一次的峰移动等较大的相位变化则利用相位修正也不能充分地修正，但通过反复进行循环处理相位修正错误得到改善。但是也确认了，在第一AFI处理的情况下即使反复进行4次以上的循环处理也没有显著的改善。
下面，示出通过模拟确定合适的零差滤波的增益的例子。
如果假定相位修正后的相位修正错误可以忽视，则与k空间数据的非对称取样区域对应的零差滤波的最大增益Hhomo.max的理论值为2≤Hhomo.max≤4。另外，现有技术中，与非对称取样区域对应的零差滤波的最大增益Hhomo.max为固定值2。
但是，如果在相位分布的推定错误残存、或者推定错误由于相位修正加强了时，把零差滤波的最大增益Hhomo.max设定成2～4的值，有时伪影反而会增加。因此，有时零差滤波的与非对称取样区域对应的最大增益Hhomo.max是比理论值小的值最合适。
设定零差滤波的最大增益Hhomo.max＝1相当于执行填零。因此，如果设定零差滤波的最大增益Hhomo.max＞1，则至少空间分辨率比填零改善。于是，只要零差滤波的最大增益Hhomo.max设定成例如Hhomo.max＝1.5等的1＜Hhomo.max＜2，则有希望相位修正错误不被加强，而空间分辨率比填零改善。
零差滤波的与非对称取样区域对应的最大增益Hhomo.max在第一、第二和第三AFI处理中可以分别确定成合适的值。
在零差滤波处理之前进行相位修正的第二或第三AFI处理中，如果把相位修正后的实空间数据变换成k空间数据，则信号向非取样部分漏失。即，取样部分中的信号成分成为非取样部分中的信号成分。
尤其是在第三AFI处理中，通过使用了全部频率区域的相位分布的相位修正使虚部信号成为近0的值。因此，k空间数据的实部关于k＝0几乎对称。因此，非取样部分中的k空间数据的信号强度为1/2左右。因此，在进行第三AFI处理时，如果相位修正错误可以忽视，则为了通过使非取样部分的信号强度为0而改善图像模糊，把零差滤波的最大增益Hhomo.max设定为2×2＝4是妥当的。
在第二AFI处理中也是，由于信号向非取样部分漏失，如果相位修正错误可以忽视，则把零差滤波的最大增益Hhomo.max设定为比2大的值是合适的。
由以上可知，在进行第一AFI处理时，即，在零差滤波处理之后进行相位修正处理时，为了避免相位修正错误被加强，优选地，把零差滤波的最大增益Hhomo.max设定成1＜Hhomo.max≤2或1＜Hhomo.max＜2。另外，在进行第二或第三AFI处理时，即，在零差滤波处理之前进行相位修正处理时，如果考虑相位修正错误，则优选地，把零差滤波的最大增益Hhomo.max设定成2≤Hhomo.max≤4或2＜Hhomo.max≤4。这样，可以把零差滤波的最大增益Hhomo.max可变地设定成2以外的值。
下面，说明着眼于相位修正错误的程度确定零差滤波的最大增益Hhomo.max的例子。
图15是示出与非对称取样区域对应的零差滤波的增益Hhomo.max＝2的第三AFI处理的模拟结果的图。图16是示出与非对称取样区域对应的零差滤波的增益Hhomo.max＝4的第三AFI处理的模拟结果的图。
图15(A)示出与非对称取样区域对应的最大增益Hhomo.max＝2的零差滤波的k方向上的强度(增益)Hhomo的分布。如果针对相位修正后的k空间数据适用图15(A)所示的零差滤波，不伴随循环处理地(j＝0)执行第三AFI处理，则生成图15(B)所示的实空间数据I0(x)。另外，如果适用图15(A)所示的零差滤波，伴随4次循环处理地(j＝4)执行第三AFI处理，则生成图15(C)所示的实空间数据I4(x)。
另一方面，图16(A)示出与非对称取样区域对应的最大增益Hhomo.max＝4的零差滤波的k方向上的强度(增益)Hhomo的分布。如果针对相位修正后的k空间数据适用图16(A)所示的零差滤波，不伴随循环处理地(j＝0)执行第三AFI处理，则生成图16(B)所示的实空间数据I0(x)。另外，如果适用图16(A)所示的零差滤波，伴随4次循环处理地(j＝4)执行第三AFI处理，则生成图16(C)所示的实空间数据I4(x)。
根据图15和图16可以看出，不进行循环处理时，如果把零差滤波的最大增益Hhomo.max设定为4，则与把零差滤波的最大增益Hhomo.max设定为2时相比，高频部分的信号强度增加。另一方面，如果进行循环处理，不管是把零差滤波的最大增益Hhomo.max设定为2还是设定为4都生成同样的信号强度分布的实空间数据I(x)。
考虑这是因为，在循环处理中把实空间数据变换成k空间数据后，k空间数据的一部分被替换成原数据。即，因为通过循环处理使相位修正后的信号的虚部为0，相位修正错误得以修正，所以与零差滤波的最大增益Hhomo.max无关地，生成同等信号强度分布的实空间数据I(x)。
另外，分别计算了改变零差滤波的最大增益Hhomo.max和循环处理的次数j而生成的各实空间数据I(x)的针对全部k空间数据Sfull(k)的平均二次方差(RMSE)。其结果为，Hhomo.max＝2，j＝0时，RMSE＝1.44；Hhomo.max＝2，j＝4时，RMSE＝1.55；Hhomo.max＝4，j＝0时，RMSE＝1.98；Hhomo.max＝4，j＝4时，RMSE＝1.54。即，在不执行循环处理时，设与非对称取样区域对应的零差滤波的最大增益Hhomo.max＝2可以减小RMSE。
因此，基于图15和图16所示的模拟，只要满足空间分辨率、相位修正错误的条件，从减小RMSE的观点看，与非对称取样区域对应的零差滤波的最大增益设成为Hhomo.max＝2是妥当的。
图17是示出在第三AFI处理中改变与非对称取样区域对应的零差滤波的最大增益Hhomo.max时RMSE的变化的图；图18是说明确定零差滤波的强度分布形状的参数的图。
图17中，横轴表示零差滤波的最大增益Hhomo.max，纵轴表示-128≤x≤128这256个象素的RMSE的总和。另外，图18中，横轴表示k空间的位置k，纵轴表示零差滤波的强度Hhomo。
使最大增益Hhomo.max以0.25为刻度从0.5变化到4.0，对第三AFI处理(即使用了全部频率区域的相位分布Φwhole(x)的相位修正后的k空间数据)施加了零差滤波。然后，在对取样区域的边界Kc＝16的部分的原k空间数据Sorig(k)(-Kc＝-16≤k≤Kmax＝128)通过第三AFI处理生成的-128≤x≤128整个区域中的实空间数据Ix的谱线(profile)与对应于完全取样数据Sfull(k)(-Kmax-128≤k≤Kmax＝128)的-128≤x≤128整个区域中的实空间数据的谱线之间，计算了RMSE。
在图17的例子中，零差滤波的最大增益Hhomo.max＝1.75时RMSE最小。RMSE最小的零差滤波的最大增益Hhomo.max的理论值为4。因此，在图17的例子中，RMSE最小的零差滤波的最大增益Hhomo.max比理论值大幅度减小。其理由，考虑是因为在相位的空间分布中高频成分多。
由以上可知，在第二和第三AFI处理中，如果着眼于信号向非取样部分漏失，则把零差滤波的最大增益设定成2≤Hhomo.max≤4是妥当的。
另一方面，在相位的空间分布中高频成分多，相位修正错误不能忽视时，有时为了使RMSE最小化把零差滤波的最大增益Hhomo.max设定成1＜Hhomo.max＜2是妥当的。因此，也可以根据相位修正用的相位的频率分布把零差滤波的最大增益Hhomo.max的能够取值范围设定为1＜Hhomo.max≤2或2≤Hhomo.max≤4。
另外，在图17中，在RMSE的计算中使用的确定零差滤波的强度分布形状Hhomo的参数K1＝16。即，像图18所示的那样，确定零差滤波的强度分布形状Hhomo的参数除了取样区域的边界Kc以外，还有K1和K2。
取样区域的边界Kc是零差滤波的强度Hhomo的上升沿的起始点。另外，在从边界Kc前进k空间上的距离K1大小的位置处成为零差滤波的强度Hhomo一定的最大值Hhomo.max。另外，零差滤波的强度Hhomo为最大值Hhomo.max的1/2时的k空间上的位置与零差滤波的强度Hhomo为最大值Hhomo.max时的k空间上的位置之间的距离用K2表示。另外，K1的值越小，则滤波的强度分布形状越接近矩形。
在图17的例子中，基于实空间数据I(x)的整个谱线的差计算了RMSE，但也可以只基于相位错误大的区域的谱线的差计算RMSE。
与非对称取样区域对应的零差滤波的最大增益Hhomo.max的最合适值与作为零差滤波的对象的数据有关，严格上讲，除非是被对称取样的k空间数据否则不能计算。于是，除了上述的基于AFI处理的种类、RMSE的方法以外，可以用各种方式确定零差滤波的最大增益Hhomo.max。
例如，也可以与取样区域的边界Kc同样地，使零差滤波的最大增益Hhomo.max缓慢地变化，通过AFI处理生成与多个不同的最大增益Hhomo.max对应的图像数据，把生成的多个图像数据显示到显示装置34上。此时，通过从输入装置33向AFI处理条件设定部42输入与相位分布对应的图像质量更好的图像数据的选择信息，AFI处理条件设定部42可以设定与图像质量良好的图像数据对应的零差滤波的最大增益Hhomo.max作为AFI处理条件。图像质量良好的图像数据的选择可以通过使用者的目视进行，也可以通过图像处理自动地进行。
或者，也可以用预定的数据收集条件与AFI中用的序列一致的相同种类的序列，基于以AFI用的MR数据收集条件与预定的条件一致的数据收集条件从同一对象脏器或同一拍摄部位收集的其它MR数据，事先用实验方法确定零差滤波的最大增益Hhomo.max。由此，可以确定与对相位分布有影响的拍摄条件对应的合适的零差滤波的最大增益Hhomo.max。此时，对相位分布有影响的TE等的数据收集条件的一致度为所要的比例，这是很重要的。另外，在以后的对应的AFI中可以把事先确定的零差滤波的最大增益Hhomo.max作为恒定值使用。
作为再一种方法，还可以为了简化数据处理，把根据AFI处理的种类取得的零差滤波的最大增益Hhomo.max的范围的平均值确定为AFI处理条件。例如，在进行第一AFI处理时，可以把零差滤波的最大增益Hhomo.max固定在1＜Hhomo.max≤2的中间值即1.5，在进行第二或第三AFI处理时，可以把零差滤波的最大增益Hhomo.max固定在2≤Hhomo.max≤4的中间值即3。
另一方面，还可以通过计算确定合适的零差滤波的最大增益Hhomo.max。在此，说明通过使相位修正错误的指标值为最小的优化计算来计算零差滤波的最大增益Hhomo.max的例子和通过使用了表示相位修正造成的信号向k空间数据的非取样部漏失的程度的指标值的计算来计算零差滤波的最大增益Hhomo.max的例子。相位修正错误的指标值和表示相位修正造成的信号向k空间数据的非取样部漏失的程度的指标值都是与相位修正中使用的相位分布以及对相位分布有影响的拍摄条件对应的值。
相位修正错误的指标值例如可以用部分的原k空间数据Sorig(k)的通过填零生成的实空间数据I0fill(x)和通过AFI处理生成的实空间数据IAFI(x)的差值来定义。然后，在取得零差滤波的最大增益Hhomo.max的范围内例如1＜Hhomo.max≤2或2≤Hhomo.max≤4中，可以通过使相位修正错误的指标值为最小的优化计算来确定零差滤波的最大增益Hhomo.max。
图19是说明确定零差滤波的最大增益Hhomo.max以使得基于通过填零生成的数据的相位修正错误的指标值最小的方法的图。
图19中，横轴表示象素位置x，纵轴表示通过AFI处理生成的实空间数据IAFI(x)和通过填零生成的实空间数据I0fill(x)在象素位置x处的差的绝对值Isub(x)。像图19所示的那样，对AFI实空间数据IAFI(x)和填零实空间数据I0fill(x)的差的绝对值Isub(x)设定阈值Ith，可以把超过阈值Ith的绝对值Isub(x)的总和作为相位修正错误的指标值ErrorSum。
即，可以用式(13)所示的算法计算相位修正错误的指标值ErrorSum。
Figure BDA0000095555280000281
然后，通过进行在可取得零差滤波的最大增益Hhomo.max的范围内变化的收敛计算，以使相位修正错误的指标值ErrorSum为最小，可以求出更合适的零差滤波的最大增益Hhomo.max。
另外，通过针对填零实空间数据I0fill(x)与AFI实空间数据IAFI(x)的差的阈值处理，可以选择性地抽出差大的范围用于相位修正错误的指标值ErrorSum的计算。由此可以避免在用填零实空间数据I0fill(x)与AFI实空间数据IAFI(x)的全部位置处的差计算相位修正错误的指标值ErrorSum时产生的模糊。
另外，在式(13)所示的算法中，把AFI实空间数据IAFI(x)与填零实空间数据I0fill(x)的振幅差的绝对值Isub(x)加起来得到值作为相位修正错误的指标值ErrorSum，但也可以基于AFI实空间数据IAFI(x)与填零实空间数据I0fill(x)的相位差定义相位修正错误的指标值ErrorSum。此时也可以通过对相位差设定阈值减少模糊。
作为另一个例子，也可以把掩蔽了与填零实空间数据I0fill(x)的相位分布图超过预定阈值的范围对应的区域的、AFI实空间数据IAFI(x)以及填零实空间数据I0fill(x)的二次方平均值定义为相位修正错误的指标值ErrorSum，进行收敛计算以使该指标值ErrorSum为最小，从而计算最合适的零差滤波的最大增益Hhomo.max。
为了提高图像质量，也可以每次执行AFI时都执行上述的使相位修正错误的指标值ErrorSum为最小的收敛计算。但是，该收敛计算伴随着相位修正处理和零差滤波处理的反复进行，所以导致数据处理时间的增加。于是，也可以事先用与AFI中用的序列相同种类的序列，从同一对象脏器或同一拍摄部位收集数据，通过基于用相同种类的序列收集的数据的收敛计算求出相位修正错误的指标值ErrorSum为最小的零差滤波的最大增益Hhomo.max。由此，可以抑制数据处理时间的增加。
另外，在AFI或先于AFI进行的成像中进行多体层数据收集时，严格地说，每个体层的收敛计算的结果不同。于是，也可以把与多个体层对应的零差滤波的最大增益Hhomo.max的平均值作为AFI处理用的零差滤波的最大增益Hhomo.max，或者针对一个或少数个代表性的体层求零差滤波的最大增益Hhomo.max。由此，可以抑制数据处理时间的增加。
另外，在进行用来减少相位修正错误的循环处理时，循环处理的反复进行等价于优化零差滤波的最大增益Hhomo.max以使相位修正错误的指标值ErrorSum为最小。因此，进行循环处理时，不需要确定零差滤波的最大增益Hhomo.max的最优值以使相位修正错误的指标值ErrorSum最小的过程。
但是，即使在进行循环处理时，零差滤波的最大增益Hhomo.max越接近最优值，循环处理的收敛时间越短。因此，为了缩短循环处理的时间，也可以进行使相位修正错误的指标值ErrorSum为最小的零差滤波的最大增益Hhomo.max的优化。
下面，说明用在第二和第三AFI处理中生成的针对k空间数据的向非取样部漏失的程度的指标值确定合适的零差滤波的最大增益Hhomo.max的方法。与使用相位修正错误的指标值ErrorSum的方法相比，如果用该方法，则可以更简单地确定合适的零差滤波的最大增益Hhomo.max。
针对相位修正造成的k空间中的信号漏失程度的指标，例如，可以定义为相位修正后的k空间数据SPcor(k)的信号的绝对值的对称度。
图20是示出相位修正后的k空间数据SPcor(k)的信号的绝对值的例图；图21是示出图20所示的相位修正后的k空间数据SPcor(k)的生成中使用的原始数据的图。
图20和图21中，各横轴表示k空间的1D方向的位置k，纵轴表示k空间数据S(k)的信号值的绝对值。
图20(A)是第二AFI处理中的即用低频区域的相位进行相位修正后的k空间数据SPcor_low(k)的信号的绝对值|SPcor_low(k)|的谱线的放大图，图20(B)是第三AFI处理中的即用全部频率区域的相位进行相位修正后的k空间数据SPcor_whole(k)的信号的绝对值|SPcor_whole(k)|的谱线的放大图。
另外，从图21(A)所示的全部取样k空间数据Sfull(k)的绝对值|Sfull(k)|，切出Kc＝16的-Kc≤k≤Kmax的范围后的部分的原k空间数据Sorig(k)的绝对值|Sorig(k)|是图21(B)。然后，通过图21(B)所示的部分的原k空间数据Sorig(k)的相位修正生成图21(A)、(B)所示的相位修正后的各k空间数据SPco_rlow(k)、SPcor- whole(k)。
像图20(A)所示的那样，可以看出，如果进行使用了低频区域的相位的相位修正，则取样区域的信号成分向k＜-Kc的非取样区域有若干漏失。另一方面，像图20(B)所示的那样，可以看出，如果进行使用了全部频率区域的相位的相位修正，则取样区域的信号成分向k＜-Kc的非取样区域有漏失从而对于k＝0成为大致对称。
像图20(A)、(B)所示的那样，从取样区域向非取样区域漏失的信号成分强度越大，取样区域中的信号成分的强度越低。于是，只要把零差滤波的最大增益Hhomo.max设定成，向非取样区域漏失的信号成分强度越大，在取得零差滤波的最大增益Hhomo.max的范围内越大即可。即，可以根据向非取样区域漏失的信号成分的强度确定合适的零差滤波的最大增益Hhomo.max。
在此，作为表示信号向非取样区域漏失的程度的指标，示出两个定义的例子。第一指标SymIndex1可以基于相位修正前的部分的原k空间数据Sorig(k)和相位修正后的k空间数据SPcor(k)定义。即，像式(14)所示的那样，第一指标SymIndex1可以定义为，部分的原k空间数据Sorig(k)的取样区域的非对称部分Kc≤k≤Kmax中的信号强度的绝对值|Sorig(k)|的积分值与相位修正后的k空间数据SPcor(k)的取样区域的非对称部分-Kc≤k≤Kmax中的信号强度的绝对值|SPcor(k)|的积分值的比。
SymIndexl = ∫ Kc K max | S Pcor ( k ) | dk / ∫ Kc K max | S orig ( k ) | dk - - - ( 14 )
如果像式(14)那样定义第一指标SymIndex1，只要没有信号向非取样区域漏失，则由于取样区域的非对称部分Kc≤k≤Kmax中的信号强度在相位修正前后无变化，所以SymIndex1＝1。相反，如果信号向非取样区域有漏失从而使得对于k＝0成为对称，则由于相位修正后的k空间数据SPcor(k)的取样区域的非对称部分Kc≤k≤Kmax中的信号强度为1/2，所以SymIndex1＝0.5。
即，0.5≤SymIndex1≤1。因此，可以用例如式(15)确定零差滤波的最大增益Hhomo.max。
Hhomo.max＝2a1/SymIndex1
∴2a1≤Hhomo.max≤4a1                     (15)
其中，式(15)中，a1是0＜a1≤1的系数。可以根据错误的程度等条件调整确定系数a1。
另外，第二指标SymIndex2只基于相位修正后的k空间数据SPcor(k)定义。即，像式(16)所示的那样，第二指标SymIndex2可以定义为，相位修正后的k空间数据SPcor(k)中非取样区域-Kmax≤k≤-Kc中的信号强度的绝对值|SPcor(k)|的积分值与取样区域的非对称部分Kc≤k≤Kmax中的信号强度的绝对值|SPcor(k)|的积分值的比。
SymIndex 2 = ∫ - Kc - K max | S Pcor ( k ) | dk / ∫ Kc K max | S Pcor ( k ) | dk - - - ( 16 )
如果像式(16)那样定义第二指标SymIndex2，只要没有信号向非取样区域漏失，则SymIndex2＝0。相反，如果信号向非取样区域有漏失从而使得对于k＝0成为对称，则由于非取样区域-Kmax≤k≤-Kc和取样区域的非对称部分Kc≤k≤Kmax中相位修正后的k空间数据SPcor(k)对称，所以SymIndex2＝1。
即，0≤SymIndex2≤1。因此，可以用例如式(17)确定零差滤波的最大增益Hhomo.max。
Hhomo.max＝2a2(1+SymIndex2)
∴2a2≤Hhomo.max≤4a2                   (17)
其中，式(17)中，a2是0＜a2≤1的系数。可以根据错误的程度等条件调整确定系数a2。
到此为止说明了把针对零差滤波的非对称取样部分的增益Hhomo.max设定成固定值的情况，但也可以设定成不恒定的值。
图22是示出把零差滤波的与非对称取样部分对应的增益设定成向着高频部分衰减至1的例子的图。
图22中，横轴表示k空间中的1D方向的位置k，纵轴表示位置k处的零差滤波的增益Hhomo(k)。相位修正后的k空间数据SPcor(k)的相位修正错误向着高频区域变大的可能性高。于是，通过向零差滤波的与非对称取样部分对应的增益Hhomo(k)赋予向着高频部分衰减至1的特性，可以相对地减小相位修正错误大的区域的信号强度。
另外，图22示出零差滤波的最大增益Hhomo.max(k)的取得范围为1＜Hhomo.max≤2的例子，但在零差滤波的最大增益Hhomo.max(k)的取得范围为2≤Hhomo.max≤4时，把零差滤波的增益Hhomo设定成在2≤Hhomo.max≤4范围内衰减至2即可。
像以上那样，在AFI处理条件设定部42中设定零差滤波的包含最大增益Hhomo.max(k)的形状。此外，AFI处理中使用的其它滤波的强度(增益)和形状也在AFI处理条件设定部42中设定。
图23是示出AFI处理中使用的滤波的具体形状的例子的图；图24是示出根据AFI处理的种类可变地设定AFI处理中使用的零差滤波的最大增益Hhomo.max(k)的例子的图。
图23(A)、(B)、(C)和图24(A)、(B)中，横轴表示k空间中的1D方向的位置k。
图23(A)示出用来从非对称取样的部分的原k空间数据Sorig(k)抽出对称的低频部分的k空间数据Slow(k)(-Kc≤k≤Kc)的LPF Hlow(k)的形状的例子；图23(B)示出从部分的原k空间数据Sorig(k)切出全部频率区域的全部数据切出的滤波Hwhole(k)的形状的例子；图23(C)示出为了循环处理内的合成处理而从部分的原k空间数据Sorig(k)切出-Kc≤k≤Kmax范围的数据切出滤波Hmerge(k)的形状的例子。
可以分别像式(18-1)、式(18-2)、式(18-3)和式(18-4)所示的那样，计算LPF Hlow(k)、全部频率区域的数据切出滤波Hwhole(k)、合成处理用的数据切出滤波Hmerge(k)和零差滤波Hhomo(k)。
Hlow(k)＝1：|k|≤Kc-K1
       ＝exp[(-ln2){(k-(Kc-K1))/K2}2]：Kc-K1＜|k|≤Kmax    (18-1)
Hwhole(k)＝Hlow(k)：k≤0
         ＝1：0≤k≤Kmax  (18-2)
Hmerge(k)＝Hlow(k)：k≤0
         ＝1：0≤k≤Kmax(18-3)
Hhomo(k)＝Hlow(k)：k≤0
        ＝Hhomo.max-(Hhomo.max-1)Hlow(k)：0≤k≤Kmax   (18-4)
由此，生成图23(A)、(B)和(C)所示的那样的形状的LPFHlow(k)、全部频率区域的数据切出滤波Hwhole(k)和合成处理用的数据切出滤波Hmerge(k)。另外，为了使增益的迁移部分平滑地变化，使用高斯函数。
在式(18-4)中，可以根据AFI处理的种类可变地设定零差滤波的最大增益Hhomo.max(k)。例如，在使用低频区域的相位分布Φlow(x)进行相位修正的第一和第二AFI处理的情况下考虑到相位修正错误可以使零差滤波的最大增益Hhomo.max(k)＝2，在使用全部频率区域的相位分布Φwhole(x)进行相位修正的第三AFI处理的情况下可以使零差滤波的最大增益Hhomo.max(k)＝4。
图24(A)示出零差滤波的最大增益Hhomo.max(k)＝2的第一和第二AFI处理用的零差滤波的形状Hhomo(k)。而图24(B)示出零差滤波的最大增益Hhomo.max(k)＝4的第三AFI处理用的零差滤波的形状Hhomo(k)。
式(18-1)、式(18-2)、式(18-3)和式(18-4)以及图23和图24是在k空间的负侧区域设定非取样区域时的例子。在k空间的正侧区域设定非取样区域时，只要把k空间的负侧与正侧中的滤波强度的值反转即可。即，在图23和图24中只要左右反转即可。
另外，在上述的说明中，示出了在AFI处理中即对利用零差滤波Hhomo(k)或全部频率区域的数据切出滤波Hwhole(k)进行滤波处理后的实空间数据Ihomo(x)或Iwhole(x)进行相位修正的例子，但也可以确定数据处理条件以使得作为AFI处理的前处理在k空间中进行0次相位修正和1次相位修正。即，可以在零差滤波处理和在零差滤波处理后执行的相位修正之前，进行k空间中的MR数据的0次和1次中的至少一方的相位修正。
针对k空间数据S(k)的1次相位修正相当于使k空间数据的回波峰的绝对值|Speak(k)|在k空间中的位置移动到与中心吻合的位置。因此，如果作为AFI处理的前处理进行0次相位修正和1次相位修正或至少进行0次相位修正，则可以有助于减少AFI处理中的相位修正错误。
尤其是在零差滤波处理之后进行相位修正处理的第一AFI处理的情况下，如果非对称的部分的原k空间数据Sorig(k)的相位误差大，则因AFI处理产生相位修正错误。因此，在进行第一AFI处理的情况下，有时作为AFI处理的前处理进行低次的相位修正是合适的。
另外，也可以像式(19-1)所示的那样，把对通过AFI处理生成的AFI实空间数据IAFI(x)和通过填零生成的填零实空间数据I0fill(x)用权重系数w加权并相加得到的实空间数据Iw(x)作为诊断图像数据。此时，权重系数w，作为简单的例子，可以像式(19-2)所示的那样是AFI实空间数据IAFI(x)与填零实空间数据I0fill(x)的差。
Iw(x)＝w*I0fill(x)+(1-w)*IAFI(x)  (19-1)
w＝|I0fill(x)-IAFI(x)|     (19-2)
关于是否作为AFI处理的前处理进行上述的0次和1次相位修正、以及是否进行式(19)所示的加权相加运算，也可以在AFI处理条件设定部42中作为数据处理条件设定。
下面，说明磁共振成像装置20的动作和作用。
图25是示出利用图1所示的磁共振成像装置20通过AFI和PI收集诊断图像时的流程的流程图。在此，描述通过预扫描确定非对称取样区域的边界Kc时的情形。
首先，在步骤S1中，在成像扫描之前，执行包含用来确定AFI中的数据取样区域的边界Kc的Kc-PREP扫描的预扫描。作为预扫描，除了Kc-PREP扫描以外，还有匀磁用的预扫描、用来测定PI处理用的多个线圈要素24c的灵敏度分布图的预扫描。
即，拍摄条件设定部40设定包含Kc-PREP扫描的预扫描用的拍摄条件，把设定了的拍摄条件输出到序列控制器31。由此，以与后述的成像扫描同样的流程收集必需的MR数据。也可以根据通过匀磁用的预扫描或用来测定线圈要素24c的灵敏度分布图的预扫描收集的数据，在整个宽范围的频率区域上求AFI处理的相位修正用的相位分布。
根据通过Kc-PREP扫描收集的MR数据，像图4(B)所示的那样，通过AFI处理生成图像数据。因此，如果在Kc-PREP扫描之前执行匀磁用的预扫描等的其它预扫描而求相位分布，则可以在针对通过Kc-PREP扫描收集的数据的相位修正中利用。然后，通过根据通过Kc-PREP扫描收集的图像数据选择图像质量良好的图像数据，可以确定成像扫描用的数据取样区域的边界Kc。
然后，在步骤S2中，执行与k空间中的数据取样区域-Kc≤k≤Kmax对应的MR数据的收集。即，拍摄条件设定部40设定包含用来从非对称数据取样区域-Kc≤k≤Kmax收集MR数据的希望的脉冲序列的拍摄条件。可以通过用输入装置33的操作从在显示装置34上显示的多个拍摄规程的候补选择所希望的拍摄规程，进行拍摄条件的设定。
然后，拍摄条件设定部40向序列控制器31输出包含脉冲序列的拍摄条件。这样一来，序列控制器31、静磁场用磁体21等的用来执行扫描的磁共振成像装置20的构成要素根据设定了的拍摄条件从受检体P收集与k空间的非对称取样区域-Kc≤k≤Kmax对应的成像数据。
为此，预先把受检体P置于机台37上，在用静磁场电源26激励的静磁场用磁体21(超导磁体)的拍摄区域上形成静磁场。另外，基于通过匀磁用的预扫描收集的数据，从匀磁线圈电源28向匀磁线圈22供给规定的电流，把拍摄区域上形成的静磁场均匀化。
然后，序列控制器31通过根据从拍摄条件设定部40接收的脉冲序列驱动倾斜磁场电源27、发送器29和接收器30，在放置受检体P的拍摄区域上形成倾斜磁场，并且从RF线圈24产生RF信号。
因此，受检体P的内部的因核磁共振产生的NMR信号由RF线圈24接收并赋予接收器30。接收器30从RF线圈24接收NMR信号，生成数字数据的NMR信号即原始数据。接收器30把生成的原始数据赋予序列控制器31。序列控制器31向数据处理部41赋予原始数据，数据处理部41把原始数据作为非对称取样区域-Kc≤k≤Kmax的部分的原k空间数据Sorig(k)，配置在形成在k空间数据库43中的k空间中。
另一方面，在步骤S3中，AFI处理条件设定部42设定与拍摄条件对应的合适的AFI处理的种类、用来减少相位修正错误的循环处理的次数j、零差滤波的强度和形状等的AFI处理的数据处理条件。即，确定如下事项：是在零差滤波处理之前进行相位修正处理还是在零差滤波处理之后进行相位修正处理；是在相位修正中使用全部频率区域的相位分布Φwhole(x)还是在相位修正中使用低频区域的相位分布Φlow(x)；是仅根据作为修正对象的k空间数据自身推定相位修正处理中使用的相位分布Φ(x)，还是仅根据其它数据推定，还是根据自数据和其它数据这两者推定等。
例如，AFI处理条件设定部42可以通过从拍摄条件设定部40取得拍摄规程的选择信息，从AFI条件数据库46读入预先与被选择的拍摄规程相关联的AFI处理条件，来设定AFI处理条件。另外，在通过基于部分的原k空间数据Sorig(k)的收敛计算等的计算求出与零差滤波的非对称取样部对应的最大增益Hhomo.max(k)时，在AFI处理条件设定部42中计算最大增益Hhomo.max(k)。
另外，如果设定AFI处理的种类，则构成AFI处理的相位修正、滤波处理等的各种处理和PI处理的顺序也由AFI处理条件设定部42确定。例如，为了缩短处理时间，把数据处理顺序确定成，在希望在与多个线圈要素24c对应的图像数据的合成处理之前执行的零差滤波处理等的AFI处理中的前级处理之后执行PI展开处理和图像数据的合成处理，在合成处理之后执行AFI处理的循环处理等的剩余的后段的处理。
而且，根据拍摄条件如果有必要，在AFI处理条件设定部42中把作为AFI处理的前处理进行0次和1次相位修正处理设定为数据处理条件。
然后，AFI处理条件设定部42把包含前处理信息的AFI处理条件和数据处理顺序赋予AFI处理部41B。
然后，在步骤S4中，AFI处理部41B根据AFI处理条件执行AFI处理的前处理和前级处理。
然后，在步骤S5中，PI处理部41A执行PI展开处理和图像数据的合成处理。把处理后的数据从PI处理部41A赋予AFI处理部41B。
然后，在步骤S6中，AFI处理部41B根据AFI处理条件执行循环处理等的AFI处理的后级处理。由此，生成与通过对称取样的k空间数据的图像重建处理生成的图像数据同等图像质量的诊断图像数据。
然后，在步骤S7中，数据处理部41把通过PI处理和AFI处理生成的诊断用的图像数据显示到显示装置34上。由此，使用者可以观察根据拍摄条件以合适的条件AFI处理了的高图像质量的诊断图像。
即，以上那样的磁共振成像装置20在进行AFI时，可以根据拍摄条件、相位分布适当地、可变地设定数据处理条件。
因此，根据磁共振成像装置20，在AFI中可以比现有技术提高图像的精度。即，可以提供与根据在k空间中被对称取样的MR数据生成的图像更接近的图像。
另外，如果在零差滤波处理前进行相位修正处理，则由于可以加快用来减少错误的循环处理的收敛，所以可以把为了提高图像精度所需的数据处理的时间的增加抑制到最小限度。
以上，描述了特定的实施方式，但在此描述的实施方式都不过是一个例子，并非用来限定发明的范围。在此描述的新的方法和装置可以以其它的各种方式实现。另外，在此描述的方法和装置的方式，在不脱离发明的主要构思的范围内，可以进行各种省略、替换、变更。这样的各种方式及其变形例都包含在发明的范围和主要构思内，且包含在权利要求书记载的发明及其等价的范围内。

Claims (18)

1.一种磁共振成像装置，其特征在于包括：
成像单元，从受检体收集与在k空间中在波数方向上为非对称的取样区域对应的磁共振数据，通过包含相位修正和取得复共轭的滤波处理的数据处理，生成基于上述磁共振数据的图像数据；以及
数据处理条件设定单元，根据上述相位修正中使用的相位分布或对上述相位分布有影响的拍摄条件，设定上述数据处理的条件。
2.如权利要求1所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为设定以下条件中的至少一种：上述取样区域的边界、上述滤波处理的最大增益、上述滤波处理中使用的滤波的形状、用来减少错误的循环处理的次数、是在上述滤波处理之前进行上述相位修正还是在上述滤波处理之后进行上述相位修正、是在上述相位修正中使用与数据处理范围对应的频率区域的相位分布还是在上述相位修正中使用与上述取样区域中的对称部分对应的频率区域的相位分布、是只根据上述磁共振数据求得上述相位修正中使用的相位分布还是只根据其它磁共振数据求得还是根据上述磁共振数据和上述其它磁共振数据这两者求得；
上述循环处理是通过反复进行取得上述相位修正和上述滤波处理后的实空间数据的实部的处理、使上述实部的相位复原以成为相位修正前的相位的处理、把上述相位复原的处理后的实空间数据变换得到的k空间数据的一部分置换成上述磁共振数据的处理、把上述置换后的k空间数据变换成实空间数据进行上述相位修正的处理，而把作为上述取得实部的处理的对象的实空间数据的虚部收敛到0的收敛计算。
3.如权利要求1所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：能够把上述滤波处理的最大增益可变地设定成2以外的值。
4.如权利要求3所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：在上述相位修正中使用与数据处理范围对应的频率区域的相位分布，且把上述滤波处理的最大增益设定成2以上且4以下。
5.如权利要求3所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：在上述相位修正中使用与上述取样区域中的对称部分对应的频率区域的相位分布，且把上述滤波处理的最大增益设定成1以上且2以下。
6.如权利要求3所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：把上述滤波处理的最大增益设定成1.5或3。
7.如权利要求3所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：根据上述相位分布把上述滤波处理的最大增益的能够取值范围设定为1以上且2以下，或2以上且4以下。
8.如权利要求3所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：通过在上述滤波处理的最大增益的能够取值范围内使表示相位修正错误的程度的指标值最小的优化计算来设定上述滤波处理的最大增益。
9.如权利要求3所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：根据表示信号向上述磁共振数据的非取样部漏失的程度的指标设定上述滤波处理的最大增益。
10.如权利要求3所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：基于以预定的条件与上述磁共振数据的收集条件一致的数据收集条件，从同一对象脏器或同一拍摄部位收集的其它磁共振数据，设定上述滤波处理的最大增益。
11.如权利要求1所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：根据上述相位分布把上述滤波处理的增益的特性设定成向着高频区域衰减的特性。
12.如权利要求3所述的磁共振成像装置，其特征在于：
上述成像单元构成为：用根据与上述图像数据的生成中使用的上述磁共振数据不同且在比上述磁共振数据更宽的频率区域中收集的其它磁共振数据求得的相位分布，进行上述相位修正。
13.如权利要求12所述的磁共振成像装置，其特征在于：
上述成像单元构成为：用根据通过用来测定静磁场的匀磁用的磁场分布图的预扫描或用来测定多个线圈要素的灵敏度分布的预扫描收集的磁共振数据求得的相位分布，进行上述相位修正。
14.如权利要求1所述的磁共振成像装置，其特征在于：
上述数据处理条件设定单元构成为：把上述数据处理的条件设定成执行包含收敛计算的循环处理，且在上述滤波处理之前进行上述相位修正，该收敛计算通过反复进行取得上述相位修正和上述滤波处理后的实空间数据的实部的处理、使上述实部的相位复原以成为相位修正前的相位的处理、把上述相位复原处理后的实空间数据变换得到的k空间数据的一部分置换成上述磁共振数据的处理、把上述置换后的k空间数据变换成实空间数据进行上述相位修正的处理，而把作为上述取得实部的处理的对象的实空间数据的虚部收敛到0。
15.如权利要求1所述的磁共振成像装置，其特征在于：
上述成像单元构成为：对使用了多个线圈要素的多个磁共振数据进行收集而生成与上述多个线圈要素对应的多个图像数据，在包含收敛计算的循环处理之前执行上述多个图像数据的合成处理，该收敛计算通过反复进行取得上述相位修正和上述滤波处理后的实空间数据的实部的处理、使上述实部的相位复原以成为相位修正前的相位的处理、把上述相位复原处理后的实空间数据变换得到的k空间数据的一部分置换成上述磁共振数据的处理、把上述置换后的k空间数据变换成实空间数据进行上述相位修正的处理，而把作为上述取得实部的处理的对象的实空间数据的虚部收敛到0。
16.如权利要求1所述的磁共振成像装置，其特征在于：
还包括改变上述取样区域的边界而收集磁共振数据，生成与多个不同的取样区域的边界对应的多个图像数据的预扫描单元；
上述数据处理条件设定单元构成为：基于根据上述多个图像数据选择图像数据的选择信息，把上述取样区域的边界设定为上述数据处理的条件。
17.如权利要求1所述的磁共振成像装置，其特征在于：
还包括存储与多个拍摄条件对应的多个数据处理条件的存储单元，
上述数据处理条件设定单元构成为：从上述存储单元取得与对上述相位分布有影响的拍摄条件对应的数据处理条件。
18.如权利要求1所述的磁共振成像装置，其特征在于：
上述成像单元构成为：在上述滤波处理和在上述滤波处理后执行的上述相位修正之前，进行k空间中的上述磁共振数据的0次和1次中的至少一种的相位修正。
CN201110295782.2A 2010-09-28 2011-09-28 磁共振成像装置 Active CN102415882B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2010-217981 2010-09-28
JP2010217981A JP5502682B2 (ja) 2010-09-28 2010-09-28 磁気共鳴イメージング装置

Publications (2)

Publication Number Publication Date
CN102415882A true CN102415882A (zh) 2012-04-18
CN102415882B CN102415882B (zh) 2014-03-12

Family

ID=45869996

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110295782.2A Active CN102415882B (zh) 2010-09-28 2011-09-28 磁共振成像装置

Country Status (3)

Country Link
US (1) US8648599B2 (zh)
JP (1) JP5502682B2 (zh)
CN (1) CN102415882B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103536360A (zh) * 2012-07-09 2014-01-29 西门子公司 由医学图像数据组提取数据组的方法及医学图像拍摄装置
CN103747728A (zh) * 2012-08-04 2014-04-23 株式会社东芝 磁共振成像装置以及磁共振成像方法
CN104644172A (zh) * 2013-11-25 2015-05-27 深圳先进技术研究院 时间反转稳态进动快速成像序列中黑血的控制方法和系统
CN105334479A (zh) * 2014-06-10 2016-02-17 西门子(深圳)磁共振有限公司 一种磁共振成像方法和装置
CN106999093A (zh) * 2015-02-06 2017-08-01 株式会社日立制作所 磁共振成像装置以及磁共振成像方法
CN109581254A (zh) * 2017-09-29 2019-04-05 西门子(深圳)磁共振有限公司 相位偏差获取方法及系统、相位校准方法及系统
CN116744844A (zh) * 2021-03-24 2023-09-12 科隆有限责任公司 用于mri系统的rf接收线圈电路

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8427147B2 (en) * 2009-04-10 2013-04-23 Wisconsin Alumni Research Foundation Magnetic resonance imaging with fat suppression by combining phase rotating data with phase shifted data in K-space
JP5670159B2 (ja) 2009-11-26 2015-02-18 株式会社東芝 磁気共鳴イメージング装置
DE102010012948B4 (de) * 2010-03-26 2012-04-26 Siemens Aktiengesellschaft Verfahren zum Ermitteln von Phasenkorrekturparametern und Magnetresonanzvorrichtung
JP5502682B2 (ja) * 2010-09-28 2014-05-28 株式会社東芝 磁気共鳴イメージング装置
WO2013118117A1 (en) * 2012-02-08 2013-08-15 Anatech Advanced Nmr Algorithms Technologies Ltd Method and system for inspection of composite material components
WO2014042026A1 (ja) 2012-09-11 2014-03-20 株式会社 東芝 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
WO2014120953A1 (en) * 2013-01-31 2014-08-07 The Regents Of The University Of California Method for accurate and robust cardiac motion self-gating in magnetic resonance imaging
KR102038628B1 (ko) * 2013-02-19 2019-10-30 삼성전자주식회사 Ｒｆ트랜스시브 코일 및 rf 리시브 코일을 이용한 ｂ１ 자기장의 쉬밍 방법, 장치 및 자기 공명 영상 시스템.
US10534060B2 (en) * 2013-06-10 2020-01-14 Toshiba Medical Systems Corporation Parallel MRI with spatially misregistered signal
DE102013216529B4 (de) * 2013-08-21 2019-05-23 Siemens Healthcare Gmbh Verfahren insbesondere zur patientenadaptiven B0-Homogenisierung von MR-Systemen unter Verwendung unterschiedlicher Typen von Shim-Spulen
KR101517606B1 (ko) * 2013-10-29 2015-05-04 한국과학기술원 K-공간의 저주파 정보를 기반으로 설계된 필터뱅크를 이용한 병렬영상기법
DE102014206917B3 (de) * 2014-04-10 2015-09-03 Siemens Aktiengesellschaft Reduzierung von Artefakten in MR-Bildern aufgrund von Suszeptibilitätssprüngen durch Ersetzen von Rohdaten aus Messungen mit einer zusätzlichen Dephasierung der Spins
JP6495057B2 (ja) * 2015-03-16 2019-04-03 キヤノンメディカルシステムズ株式会社 Ｍｒｉ装置及び撮像時間短縮方法
JP7455508B2 (ja) * 2018-12-26 2024-03-26 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置および医用複素数画像処理装置
EP3851866A1 (en) * 2020-01-16 2021-07-21 Siemens Healthcare GmbH In-phase and out-of-phase magnetic resonance imaging
DE102020211107A1 (de) * 2020-09-03 2022-03-03 Siemens Healthcare Gmbh Erzeugung von kombinierten Bilddaten basierend auf MR-Daten
US20220248972A1 (en) * 2021-02-10 2022-08-11 GE Precision Healthcare LLC System and method for reconstruction using a high-resolution phase in magnetic resonance images

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2307746A (en) * 1995-11-28 1997-06-04 Univ California A programmable digital low pass filter for magnetic resonance imaging
JP2004321791A (ja) * 2003-04-10 2004-11-18 Toshiba Corp Ｍｒｉ装置及びｍｒイメージング方法
JP2006061235A (ja) * 2004-08-25 2006-03-09 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2008093070A (ja) * 2006-10-10 2008-04-24 Hitachi Medical Corp 磁気共鳴イメージング装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3322943B2 (ja) 1993-07-14 2002-09-09 株式会社日立メディコ Ｍｒｉ装置
JPH09108197A (ja) * 1995-10-19 1997-04-28 Hitachi Medical Corp 磁気共鳴測定方法およびその装置
JP5366370B2 (ja) * 2006-09-06 2013-12-11 株式会社東芝 磁気共鳴イメージング装置
JP5666779B2 (ja) 2009-01-27 2015-02-12 株式会社東芝 磁気共鳴イメージング装置
JP5377220B2 (ja) 2009-10-23 2013-12-25 株式会社東芝 磁気共鳴イメージング装置
JP5502682B2 (ja) * 2010-09-28 2014-05-28 株式会社東芝 磁気共鳴イメージング装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2307746A (en) * 1995-11-28 1997-06-04 Univ California A programmable digital low pass filter for magnetic resonance imaging
JP2004321791A (ja) * 2003-04-10 2004-11-18 Toshiba Corp Ｍｒｉ装置及びｍｒイメージング方法
JP2006061235A (ja) * 2004-08-25 2006-03-09 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2008093070A (ja) * 2006-10-10 2008-04-24 Hitachi Medical Corp 磁気共鳴イメージング装置

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103536360A (zh) * 2012-07-09 2014-01-29 西门子公司 由医学图像数据组提取数据组的方法及医学图像拍摄装置
CN103747728A (zh) * 2012-08-04 2014-04-23 株式会社东芝 磁共振成像装置以及磁共振成像方法
CN103747728B (zh) * 2012-08-04 2016-06-29 株式会社东芝 磁共振成像装置以及磁共振成像方法
US10916007B2 (en) 2012-08-04 2021-02-09 Toshiba Medical Systems Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN104644172A (zh) * 2013-11-25 2015-05-27 深圳先进技术研究院 时间反转稳态进动快速成像序列中黑血的控制方法和系统
CN105334479A (zh) * 2014-06-10 2016-02-17 西门子(深圳)磁共振有限公司 一种磁共振成像方法和装置
CN105334479B (zh) * 2014-06-10 2019-02-22 西门子(深圳)磁共振有限公司 一种磁共振成像方法和装置
CN106999093A (zh) * 2015-02-06 2017-08-01 株式会社日立制作所 磁共振成像装置以及磁共振成像方法
CN109581254A (zh) * 2017-09-29 2019-04-05 西门子(深圳)磁共振有限公司 相位偏差获取方法及系统、相位校准方法及系统
US10823810B2 (en) 2017-09-29 2020-11-03 Siemens Healthcare Gmbh Phase deviation acquisition method and system, and phase calibration method and system
CN116744844A (zh) * 2021-03-24 2023-09-12 科隆有限责任公司 用于mri系统的rf接收线圈电路

Also Published As

Publication number Publication date
US8648599B2 (en) 2014-02-11
US20120074940A1 (en) 2012-03-29
JP2012070939A (ja) 2012-04-12
JP5502682B2 (ja) 2014-05-28
CN102415882B (zh) 2014-03-12

Similar Documents

Publication Publication Date Title
CN102415882B (zh) 磁共振成像装置
Poser et al. Pulse sequences and parallel imaging for high spatiotemporal resolution MRI at ultra-high field
JP6188597B2 (ja) 磁気共鳴システムの作動方法および制御装置
CN102870000B (zh) 使用单个和多个通道接收器线圈的用于同时多切片磁共振成象的方法
US5909119A (en) Method and apparatus for providing separate fat and water MRI images in a single acquisition scan
US9770186B2 (en) Systems and methods for magnetic resonance imaging
DE102010039693B4 (de) Phasensensitive Magnetresonanz-Bildgebung
US10288703B2 (en) MRI method of hybrid acquisition in 3D TSE
US20150177353A1 (en) Method for Magnetic Resonance Imaging with Controlled Aliasing
EP0529527B1 (en) Method and apparatus for high speed magnetic resonance imaging with improved image quality
Pipe Pulse sequences for diffusion-weighted MRI
US10552955B2 (en) Imaging acceleration methods for MRI parameter mapping
CA2530298A1 (en) Methods and apparatuses for fast chemical shift magnetic resonance imaging
US9645211B2 (en) Magnetic resonance imaging using steering-propeller
CN102621510A (zh) 用于抑制mr成像中伪影的系统
US9971007B2 (en) Method and apparatus for accelerated magnetic resonance imaging
CN105785298A (zh) 一种高精度三维化学位移成像方法
JP4225620B2 (ja) Ｍｒｉ装置
US10775468B2 (en) Magnetic resonance method and apparatus
CN107076819A (zh) 具有对流伪影的抑制的Dixon MR成像
JP2000201903A5 (zh)
Sprenger et al. NeuroMix—A single‐scan brain exam
Sbrizzi et al. Dictionary-free MR Fingerprinting reconstruction of balanced-GRE sequences
Kim et al. Improved cine displacement‐encoded MRI using balanced steady‐state free precession and time‐adaptive sensitivity encoding parallel imaging at 3 T
Chang et al. Correction of geometric distortion in Propeller echo planar imaging using a modified reversed gradient approach

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160726

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.