CN103533887A - 磁共振成像装置、以及重构图像取得方法 - Google Patents

Abstract

本发明的磁共振成像装置为了在非正交测量中不牺牲速度而提高对比度和图像质量，而在非正交系测量中，在将组合由一个拍摄取得多个回波信号的高速摄像时序的摄像中，将配置各拍摄的回波列的叶片形状设为具备由以k空间的原点为中心的圆的半径和圆弧构成的扇形区域和与相邻的叶片重叠的区域。在测量时，进行控制使得将希望的TE的回波信号配置在各叶片的k空间的低空间频率区域中，在图像重构时，使用重叠的区域的数据修正叶片之间的体动。

Description

磁共振成像装置、以及重构图像取得方法

技术领域

本发明涉及一种核磁共振成像（Magnetic Resonance Imaging：以下称为MRI）技术，特别涉及一种非正交系测量技术，其测量来自被检测体中的氢、磷等的核磁共振（Nuclear Magnetic Resonance：以下称为NMR）信号，对原子核的密度分布、缓和时间分布等进行图像化。

背景技术

用于MRI的MRI装置是测量被检测体、特别是构成人体组织的原子核自旋所产生的NMR信号（回波信号），将其头部、腹部、四肢等的形态、功能进行二维或三维的图像化（摄影）的装置。回波信号作为位置信息，被赋予因倾斜磁场而不同的相位编码、频率编码，并与位置信息对应地配置在k空间中。通过对配置在k空间中的回波信号进行二维或三维傅立叶变换来重构图像。

在MRI中，测量回波信号使得沿着k空间的预定扫描轨迹取得数据。K空间的扫描轨迹由所施加的倾斜磁场图案所决定，大致分为在正交坐标系k空间上取得数据的正交系测量的轨迹、在非正交坐标系k空间上取得数据的非正交系测量的轨迹。正交坐标系k空间是通过相互正交的2轴或3轴的坐标轴所规定的二维或三维的数据空间，非正交坐标系k空间是通过大小和偏角所规定的二维或三维的数据空间。非正交系测量一边使偏角变化一边扫描k空间，因此对k空间的中心附近进行重复扫描（例如参照非专利文献1）。因此，是因呼吸等的运动引起的影响被平均化，而在特定方向上不进行成像伪影的稳妥的测量方法。

作为MRI的摄像法，有FSE法，即在施加一次激励脉冲后，到下一次激励脉冲的施加为止的TR期间，多次施加再收敛脉冲，高速地取得多个回波信号。此外，在FSE中，将一次激励脉冲的施加称为拍摄，将在一个拍摄中得到的多个回波信号称为回波列。有一种方法（混合径向法），即将非正交系测量与FSE法组合，高速地得到伪影少的图像。在混合径向法中，在被称为一个叶片的长方形的正交坐标系k空间内对各回波列进行正交系测量，在每个拍摄中使叶片在k空间内旋转。在该情况下，叶片的长轴方向与频率编码对应，短轴方向与相位编码对应。

另外，作为高速地填充k空间的测量法，有将频率编码方向的读出倾斜磁场和相位编码方向的尖峰信号倾斜磁场组合起来进行测量的EPI法。非正交系测量也可以与该EPI法组合起来。在该情况下，将叶片的短轴方向设为频率编码，将长轴方向设为相位编码（例如参照专利文献1）。通过将两者组合来抑制伪影，还可以缩短频率编码倾斜磁场的每一次施加时间缩短，从而降低图像失真。

现有技术文献

专利文献

专利文献1：美国专利第7535222号说明书

非专利文献

非专利文献1：Magnetic Resonance in Medicine42:963-969(1999).MotionCorrection With PROPELLER MRI:Application to Head Motion andFree-Breathing Cardiac Imaging.James G.Pipe.

发明内容

发明要解决的问题

重构图像的对比度由配置在k空间的中心区域（低空间频率区域）中的回波信号来决定。因此，在使用通过一个拍摄取得多个回波信号的测量法的情况下，进行控制使得所希望的对比度的回波信号配置在低空间频率区域（低频率）中。将从激励脉冲的施加到得到该希望的对比度的回波信号的时间称为实效TE。

但是，在非正交系测量中，将全部回波信号配置在k空间低频率附近，因此将在实效TE以外的时间取得的回波信号也配置在k空间的低频率中。因此，所得到的图像与所希望的对比度不同，对比度降低。

另外，如果将k空间分割为相互不重叠的多个叶片而取得数据，则在叶片之间有体动的情况下，无法掌握该体动，因此无法进行适当的修正，进而无法提高图像质量。

本发明就是鉴于上述情况而提出的，其目的在于：在非正交测量中，不牺牲速度地提高对比度，并且修正体动，还提高图像质量。

解决问题的方案

本发明的非正交系测量中，在将通过一个拍摄取得多个回波信号的高速摄像时序组合起来的摄像中，将配置各拍摄的回波列的叶片形状设为具备由以k空间的原点为中心的圆的半径和圆弧构成的扇形区域和与相邻的叶片重叠的区域。在测量时，进行控制使得将希望的TE的回波信号配置在各叶片的k空间的低空间频率区域中，在图像重构时，使用重叠的区域的数据来修正叶片之间的体动。

发明效果

根据本发明，在非正交系测量中，能够高速地修正体动而提高对比度和图像质量。

附图说明

图1是表示第一实施方式的MRI装置100的整体结构的框图。

图2是第一实施方式的信号处理系统70的功能框图。

图3是第一实施方式的摄影处理的流程图。

图4（A）是用于说明现有的摄像时序的说明图，图4（B）是用于说明第一实施方式的摄像时序的说明图。

图5（A）是用于说明现有的长方形叶片的说明图，图5（B）是用于说明第一实施方式的扇形叶片的说明图。

图6（A）是用于说明第一实施方式的扇形叶片整体的形状的说明图，图6（B）是用于说明第一实施方式的体动检测区域的说明图。

图7是第一实施方式的时序生成处理的流程图。

图8是用于说明第一实施方式的分割区域的说明图。

图9（A）和图9（B）是用于说明现有的长方形叶片的回波信号的扫描轨迹的说明图，图9（C）和图9（D）是用于说明第一实施方式的扇形叶片的回波信号的扫描轨迹的说明图。

图10是用于说明第一实施方式的钟摆状轨迹的说明图。

图11（A）是第一实施方式的图像重构处理的流程图，图11（B）是第一实施方式的图像重构处理的变形例子的流程图。

图12是用于说明第一实施方式的体动检测区域内的扫描轨迹的变形例子的说明图。

图13（A）是用于说明第一实施方式的体动检测区域的变形例子的说明图，图13（B）和图13（C）是用于说明变形例子的扫描轨迹的说明图。

图14（A）和图14（B）是用于说明第二实施方式的扇形叶片的说明图。

图15（A）是用于说明第一实施方式的扇形叶片的k空间配置的说明图，图15（B）是用于说明第二实施方式的扇形叶片的k空间配置的说明图。

图16是第二实施方式的时序生成处理的流程图。

具体实施方式

<第一实施方式>

以下，说明应用本发明的第一实施方式。以下，在用于说明本发明的实施方式的全部图中，对具有相同功能的部分附加相同符号，省略其重复的说明。

最初，根据图1说明本实施方式的MRI装置100的概要。图1是表示本实施方式的MRI装置100的整体结构的框图。本实施方式的MRI装置100利用NMR现象得到被检测体的断层图像，具备静磁场产生系统20、倾斜磁场产生系统30、时序产生器40、发送系统50、接收系统60、以及信号处理系统70。

静磁场产生系统20如果是垂直磁场方式，则在被检测体10的周围的空间，在与其体轴垂直的方向上产生均匀的静磁场，如果是水平磁场方式，则在体轴方向上产生均匀的静磁场，在被检测体10的周围配置永久磁铁方式、正常传导方式、或超导方式的静磁场产生源。

倾斜磁场产生系统30具备作为在MRI装置100的坐标系（静止坐标系）的X、Y、Z的3轴方向上卷绕的倾斜磁场线圈31、驱动各个倾斜磁场线圈31的倾斜磁场电源32，根据后述的时序产生器40的指令来驱动各个线圈的倾斜磁场电源32，由此在X、Y、Z的3轴方向施加倾斜磁场Gx、Gy、Gz。在摄像时，在与切片（摄像截面）垂直的方向施加切片方向倾斜磁场脉冲（Gs），设定与被检测体10对应的切片面，在与该切片面垂直、并且相互垂直的剩余2个方向施加相位编码方向倾斜磁场脉冲（Gp）和频率编码方向倾斜磁场脉冲（Gf），对NMR信号（回波信号）将各个方向的位置信息进行编码。

时序产生器40根据后述的信号处理系统70所具备的CPU71的控制信号来控制倾斜磁场产生系统30、发送系统50、以及接收系统60，从而使得重复施加高频磁场脉冲（以下称为“RF脉冲”）和倾斜磁场脉冲。

发送系统50为了使构成被检测体10的生物体组织的原子的原子核自旋产生核磁共振，而向被检测体10照射RF脉冲，具备高频振荡器（合成器）52、调制器53、高频放大器54、发送方的高频线圈（发送线圈）51。在时序产生器40发出的指令的定时通过调制器53对合成器52输出的高频脉冲进行振幅调制，通过高频放大器54对该振幅调制后的高频脉冲进行放大，提供给与被检测体10接近地配置的发送线圈51，由此对被检测体10照射RF脉冲。

接收系统60检测由构成被检测体10的生物体组织的原子核自旋的核磁共振而释放的回波信号（NMR信号），具备接收侧的高频线圈（接收线圈）61、信号放大器62、正交相位检波器63、A/D变频器64。通过与被检测体10接近地配置的接收线圈61检测通过由发送线圈51照射的电磁波而引起的被检测体10的响应回波信号，在信号放大器62被放大后，在时序产生器40发出的指令的定时通过正交相位检波器63被分割为正交的2个系统的信号，分别通过A/D变频器64变换为数字量，发送到信号处理系统70。

信号处理系统70进行各种数据处理和处理结果的显示以及保存等，具备CPU71、存储装置72、外部存储装置73、显示装置74、输入装置75。

本实施方式的信号处理系统70例如根据摄像时序，向时序产生器40提供控制信号，从接收系统60收集用于生成被检测体10的断层图像的数据。然后，使用收集到的数据，对正在摄影的被检测体10的运动进行修正（体动修正），构成被检测体10的断层图像。此外，根据预先保存在存储装置72等中的脉冲时序，使用操作者经由输入装置75进行输入的摄像参数来生成摄像时序。

为了实现该情况，本实施方式的信号处理系统70如图2所示，具备：受理部710，其从操作者受理摄像参数；时序生成部720，其使用所受理的摄像参数决定预先保存的脉冲时序的倾斜磁场形状，生成摄像时序；测量控制部730，其根据所生成的摄像时序，执行回波信号的测量，将测量出的回波信号配置在k空间中；图像重构部740，其使用配置在k空间中的回波信号（数据），一边进行体动修正一边重构图像。

通过在信号处理系统70中，由CPU71将预先保存在存储装置72等中的程序装载到存储器中并执行，从而实现这些功能。

显示装置74显示所重构出的断层图像，并且与输入装置75一起构成用于操作者输入各种控制信息的接口。输入装置75例如由跟踪球或鼠标、以及键盘等构成。存储装置72和外部存储装置存储由操作者输入的信息、在信号处理系统70的处理中途、以及处理时所生成的信息等。

此外，在图1中，在插入被检测体10的静磁场产生系统20的静磁场空间内，如果是垂直磁场方式则与被检测体10相对地设置发送线圈51和倾斜磁场线圈31，如果是水平磁场方式，则将发送线圈51和倾斜磁场线圈31设置成围住被检测体10。另外，接收线圈61被设置成与被检测体10相对、或围住被检测体10。

作为在临床普及的情况，现在的MRI装置的摄像对象原子核种类是被检测体为主的构成物质的氢原子核（质子）。通过对与质子密度的空间分布、激励状态的缓和时间的空间分布有关的信息进行图像化，对人体头部、腹部、四肢等的形态或功能进行二维或三维的摄像。

首先，说明通过信号处理系统70所实现的各功能实现的本实施方式的MRI装置100在摄影时的摄影处理的流程。图3是本实施方式的摄影处理的处理流程。

受理部710受理操作者经由输入装置75输入的摄像参数（步骤S1101）。在输入的摄像参数中，包含频率编码步长数（频率编码方向的采样数）F、相位编码步长数（相位编码方向的采样数）P。

然后，时序生成部720使用摄像参数决定倾斜磁场波形，根据预先保存在存储装置72等中的脉冲时序来生成摄像时序（步骤S1102）。测量控制部730根据摄像时序向时序产生器进行指示，进行测量，将数据填充到k空间中（步骤S1103）。然后，图像重构部740一边进行体动修正一边重构图像（步骤S1104）。

在本实施方式中，作为回波信号测量时的扫描区域，代替在现有的混合径向法中使用的长方形的叶片，而使用扇形的叶片、即具有与k空间的低空间频率区域共通的区域（重叠的区域）的扇形叶片。另外，在多个该扇形叶片中扫描k空间。各扇形叶片内的扫描轨迹被决定为将希望的TE的回波信号配置在k空间的低空间频率区域中，并且其形状是钟摆状。

时序生成部720决定实现这样的测量的倾斜磁场波形，生成摄像时序。测量控制部730根据所生成的摄像时序，控制各部件而执行上述测量。另外，图像重构部740使用多个扇形叶片的共通的区域（重叠区域）的数据，计算每个叶片的被检测体10的体动量，一边进行体动修正一边重构图像。

接着，说明上述步骤S1102的时序生成部720进行的摄像时序生成处理。

首先，说明现有的将非正交系测量的混合径向法和FSE法组合所得的摄像时序（称为现有法）和由此测量的k空间上的叶片。图4（A）是现有法的摄像时序210。另外，图5（A）是按照图4（A）所示的摄像时序210测量的k空间上的长方形的区域（长方形叶片）310。在此，RF、Gx、Gy的各个轴分别表示RF脉冲、2轴方向的倾斜磁场脉冲的施加定时。

在FSE法中，在施加1次激励RF脉冲211后到下一个激励RF脉冲211的施加为止的时间TR的期间，施加多个再收敛RF脉冲212，在每次施加再收敛RF脉冲时取得回波信号。这时，向各回波信号赋予不同的相位编码。

将其与混合径向法组合，在一次的TR期间（拍摄），将图5（A）所示的包含k空间的原点的长方形的区域（长方形叶片）310作为单位区域，在每个拍摄中使长方形叶片310与k空间的kx轴所成的角度（旋转角θ）变化地重复对单位区域内进行测量的单位测量，来测量k空间全体。在现有法中，决定倾斜磁场脉冲213、214的波形使得其得以实现。如上述那样，由操作者将1个长方形叶片310的测量（1次拍摄）的重复次数B（B是自然数）设定为摄像参数。

在此，将在第b（b是满足1≤b≤B的自然数）次重复中测量的长方形叶片（第b个长方形叶片）设为310（b），将该长方形叶片310（b）的频率编码方向设为kx（b）轴，将相位编码方向设为ky（b）轴。将k空间的x轴（kx轴）和kx（b）轴所成的角度设为长方形叶片310（b）的旋转角度θ（b）。如上述那样，各长方形的叶片310（b）的kx（b）轴方向的采样数F以及ky（b）轴方向的采样数P分别由操作者设定为摄像参数（频率编码数和相位编码数）。

与此相对，在本实施方式中，在FSE的一个拍摄期间中，作为单位测量，对与该长方形叶片310（b）相同面积（相同采样数）的扇形叶片、即具有与k空间的低空间频率区域共通的区域的扇形叶片（单位区域）内进行测量。另外，在每个拍摄中使该扇形叶片与kx轴所成的角度（旋转角度θ）变化而重复测量，来测量k空间全体。这时，决定旋转角度θ使得各扇形叶片在圆周方向上不重叠。图4（B）表示实现它的本实施方式的摄像时序220。图6（A）表示本实施方式的扇形叶片320的整体的形状，图5（B）表示其配置。

本实施方式的扇形叶片320，本实施方式的扇形叶片320如图6（A）所示，具备收集用于重构图像的数据的扇形区域322、收集用于修正扇形叶片320之间的体动用的数据的体动检测区域323。

扇形区域322是以k空间的原点为中心的半径R的圆的被2个半径R和2个半径R之间的圆弧围住的区域，一部分与体动检测区域323重叠。如图5（B）所示，扇形区域322的中心角为φ（rad）。第b个扇形叶片320（b）的by（b）轴设为将扇形叶片320（b）的扇形区域322（b）的中心角φ进行二分割的线（中心线）上。kx（b）轴为与ky（b）轴垂直的方向。另外，扇形叶片320（b）的旋转角度θ（b）为kx（b）轴与kx轴所成的角度。

如图6（B）所示，体动检测区域323在每个拍摄的全部扇形叶片320所共通的重复区域中，被设定为以k空间的原点为中心的预先确定的低空间频率区域。在本实施方式中，作为一个例子，为半径C的圆形区域。

本实施方式的摄像时序220也如图4（B）所示那样，基本上与现有的FSE同样地，在施加一次激励RF脉冲211后到施加下一个激励RF脉冲211为止的时间TR的期间，施加多个再收敛RF脉冲212，在每次施加再收敛RF脉冲时取得回波信号。这时，如上述那样，在一次的TR期间中，对具有图6所示的形状的扇形叶片320内进行测量，并且在每个TR中使扇形叶片320的旋转角度θ变化而重复测量来测量k空间全体，决定Gx轴和Gy轴方向的倾斜磁场脉冲223、224的波形。进而，决定Gx轴和Gy轴方向的倾斜磁场脉冲223、224的波形使得在各扇形叶片320内，将希望的对比度的回波信号配置在k空间的低空间频率区域中。

以下，使用图7说明本实施方式的时序生成部720决定的上述那样的倾斜磁场形状，和生成摄像时序的摄像时序生成处理的流程。

首先，作为叶片参数，时序生成部720使用摄像参数计算并决定扇形区域322的半径R、中心角φ以及总叶片数B_fan（叶片参数决定：步骤S1201）。总叶片数B_fan是为了对k空间全体进行一次测量所需要的扇形叶片320的个数。

首先，为了使最终的全部采样面积与基于长方形叶片310的全部采样面积相同，将半径R设为F/2（R=F/2）。

接着，如下这样决定中心角φ。

在本实施方式中，代替现有的长方形叶片310，在每一个拍摄中对扇形叶片320内进行测量，为此与长方形叶片310的采样数同样，根据由操作者设定为摄像参数的频率编码步长数F和相位编码步长数P来决定各扇形叶片320内采样数，即FP。

首先，考虑以下的情况，即在一个拍摄中，只测量扇形叶片320的扇形区域322的面积S。如上述那样，扇形区域与长方形叶片310的面积FP相同。

因此，使用采样数F、采样数P，用以下的公式（1）表示扇形叶片320的扇形区域322的面积S。

[公式1]

$S=\mathrm{\&pi;}{\left(\frac{F}{2}\right)}^2\frac{\mathrm{\&phi;}}{2\mathrm{\&pi;}}=\mathrm{FP}---\left(1\right)$

如果对φ求解公式（1），则用以下的公式（2）表示中心角φ。

[公式2]

$\mathrm{\&phi;}=\frac{8P}{F}---\left(2\right)$

此外，在本实施方式中，将扇形叶片320的扇形区域322决定为使得在圆周方向上没有重叠地测量k空间的半径R的区域。因此，如以下的公式（3）所示，根据扇形区域322测量的k空间内的采样面积N_fan与和k空间内接的圆（半径R=F/2）的面积相等。

[公式3]

$N_{\mathrm{fan}}={\mathrm{\&pi;R}}^2=\frac{\mathrm{\&pi;}}{4}{F}^2---\left(3\right)$

由此，用以下的公式（4）表示扇形叶片的总叶片数B_fan。

[公式4]

$B_{\mathrm{fan}}=\frac{N_{\mathrm{fan}}}{S}=\frac{N_{\mathrm{fan}}}{\mathrm{FP}}=\frac{{\mathrm{\&pi;R}}^2}{\mathrm{FP}}=\frac{\mathrm{\&pi;}}{4}\frac{F}{P}---\left(4\right)$

因此，是1≤b≤B_fan。这是将2π除以中心角φ所得的。

另外，根据0≤θ（b）<2π，用以下的公式（5）表示第b个扇形叶片320（b）的旋转角度θ（b）。

[公式5]

$\mathrm{\&theta;}\left(b\right)=(b-1)\&times;\frac{2\mathrm{\&pi;}}{B_{\mathrm{fan}}}+\mathrm{\&Phi;}---\left(5\right)$

在此，Φ是规定b=1的叶片旋转角度的常数。

接着，时序生成部720根据在一个拍摄内取得的回波信号数，在半径方向上对扇形叶片320的扇形区域322进行分割，决定各回波信号扇形叶片320内的配置区域（步骤S1202）。

在此，将扇形叶片320内的回波数设为E（E是自然数；E=P）。如图8所示，将E个回波信号分别配置在对扇形叶片320的扇形区域322在半径方向上进行E分割而得到的面积相互相等的区域（分割区域）321（n）（n是满足1≤n≤E的自然数）中。以下，将第b个的各扇形叶片320（b）的扇形区域322（b）的分割区域表示为321（b，n）。

在本实施方式中，进行控制，使得希望的对比度的回波信号配置在k空间的低空间频率区域中，除此以外的回波信号配置在高空间频率区域中。因此，进行控制，使得将希望的定时（实效TE）的回波信号配置在分割区域321中的k空间的原点附近的低空间频率区域中。

扇形叶片320（b）内的每个分割区域321（b，n）的回波信号配置顺序根据希望的对比度而改变。在此，举例说明例如从k空间的中心侧按顺序取得数据的中心（centric）配置的情况。即，将配置第一个回波信号的分割区域321（b，1）设为扇形叶片320（b）的扇形区域322（b）的半径r（b，1），将配置第n个回波信号的分割区域321（b，n）设为扇形叶片320（b）的在半径r（b，n）和半径r（b，n-1）之间的区域。此外，在中心配置中，r（b，E）是扇形叶片320（b）的扇形区域322（b）的半径R（b）。

根据分割区域321（b，n）的面积相互相等，所以用以下的公式（6）表示半径r（b，n）。

[公式6]

$r(b,n)=R\left(b\right)\sqrt{\frac{n}{E}}---\left(6\right)$

接着，时序生成部720决定各分割区域321（b，n）内的扫描轨迹（步骤S1203）。

在图9（A）所示的长方形叶片310（b）中，如图9（B）所示，一个回波信号成为直线状的轨迹（直线状轨迹510）。在本实施方式中，如图9（C）所示，在扇形叶片320（b）的扇形区域322（b）内的各分割区域321（b，n）中设定扫描轨迹。各分割区域321（b，n）内的扫描轨迹如图9（D）所示，是一连串的轨迹，为钟摆状的轨迹（钟摆状轨迹520），即交替地具有与扇形叶片320（b）的扇形区域322（b）的圆弧同心圆的圆弧状部分、以及与连接相邻的2个圆弧状部分的半径方向平行的直线状部分。

将扇形叶片320（b）的扇形区域322（b）的在分割区域321（b，n）中的每个回波信号的频率编码方向的切换数设为M（b，n）。图10表示切换数M（b，n）=3的情况下设定在分割区域321（b，n）中的钟摆状轨迹520的例子。如本图所示那样，钟摆状轨迹520由交替地重复M次（在此为3次）的、圆弧部分（圆弧状轨迹）521和直线部分（直线状轨迹522）构成。

使用扇形叶片320（b）的扇形区域322（b）的分割区域321（b，n）的半径r（b，n）和切换数M（b，n），用以下的公式（7）表示与直线状轨迹522的长度相当的圆弧状轨迹521的间隔A（b，n）。

[公式7]

$A(b,n)=\frac{r(b,n)-r(b,n-1)}{M(b,n)}---\left(7\right)$

此外，各分割区域321（b，n）中的钟摆状轨迹520的采样长度L（b，n）是切换数M（b，n）个圆弧状轨迹521和直线状轨迹522的总和，所以用以下的公式（8）来表示。

[公式8]

$L(b,n)=\underset{m=1}{\overset{M(b,n)}{\mathrm{\&Sigma;}}}(r(b,n-1)+(m-1)A(b,n))\mathrm{\&phi;}\left(b\right)+\underset{m=1}{\overset{M(b,n)}{\mathrm{\&Sigma;}}}A(b,n)---\left(8\right)$

在此，m是满足1≤m≤M的整数。扇形叶片320（b）内的钟摆状轨迹520的总采样长度与长方形叶片310（b）的频率编码方向的采样数（采样长度）F相等，因此用公式（9）表示扇形叶片320（b）的扇形区域322（b）的分割区域321（b，n）的切换数M（b，n）。

[公式9]

$M(b,n)=\frac{(\mathrm{\&phi;}-2)(r(b,n)-r(b,n-1))+2F}{\mathrm{\&phi;}(r(b,n)+r(b,n-1))}---\left(9\right)$

根据以上，用以下的公式（10）表示扇形叶片320（b）的扇形区域322（b）的分割区域321（b，n）的第m（满足1≤m≤M的自然数）次切换时的圆弧状轨迹521。其中，-φ/2≤t≤φ/2。

[公式10]

$\left\{\begin{array}{c}\mathrm{kx}\left(b\right)=\mathrm{mA}(b,n)\sin \left({(-1)}^{m}t\right)\\ \mathrm{ky}\left(b\right)=\mathrm{mA}(b,n)\cos \left({(-1)}^{m}t\right)\end{array}\right.---\left(10\right)$

另外，用以下的公式（11）表示扇形叶片320（b）的扇形区域322（b）的分割区域321（b，n）中的第m次切换时的直线状轨迹522。其中，(m-1)A(n)≤t≤mA(n)。

[公式11]

$\left\{\begin{array}{c}\mathrm{kx}\left(b\right)=t\sin \left({(-1)}^m\frac{\mathrm{\&phi;}\left(b\right)}{2}\right)\\ \mathrm{ky}\left(b\right)=-\mathrm{tos}\left({(-1)}^m\frac{\mathrm{\&phi;}\left(b\right)}{2}\right)\end{array}\right.---\left(11\right)$

在本实施方式中，如上述那样，对扇形叶片320（b）之间的被检测体10的运动进行修正，因此在各扇形叶片320（b）的k空间中心部分的体动检测区域323（b）中取得体动检测用的数据。

本实施方式的体动检测区域323是以k空间原点为中心的圆。如果将其半径设为C，则用以下的公式（12）表示体动检测区域323。

[公式12]

kx(b)²+ky(b)²≤C²      (12)

本实施方式的扇形叶片320（b）除了扇形区域322（b）以外，还具备体动检测区域323（b）。因此，各分割区域321（b）的半径追加该体动检测区域323（b），而从公式（6）所示的r（b，n）如公式（13）所示的r_mc（b，n）那样变化。此外，公式（13）中的r（b，n）表示在公式（6）中规定的各分割区域321（b，n）的半径。

[公式13]

$r_{\mathrm{mc}}(b,n)=\sqrt{{r(b,n)}^2-(\frac{2\mathrm{\&pi;}}{\mathrm{\&phi;}}-1)C^2}---\left(13\right)$

另外，第一个回波信号的总采样长度L（b，l）是圆状的体动检测区域323（b）内的轨迹的M（b，n）次的圆弧状轨迹521和直线状轨迹522的总和，用以下的公式（14）表示。在此，ω表示体动检测区域323（b）的采样点之间的间隔（像素/采样）。

[公式14]

$L(b,1)=\mathrm{\&omega;T}+\underset{m=1}{\overset{M}{\mathrm{\&Sigma;}}}(m-1)A(b,n)\mathrm{\&phi;}+\underset{m=1}{\overset{M}{\mathrm{\&Sigma;}}}A(b,n)---\left(14\right)$

包含体动检测区域323（b）的分割区域321（b，l）中的轨迹的总采样长度与长方形叶片的采样长度F相等。因此，用以下的公式（15）表示M（b，l）。

[公式15]

$M(b,1)=\frac{(\mathrm{\&phi;}-2)r_{\mathrm{mc}}\left(1\right)+2(F-\mathrm{\&omega;T})}{\mathrm{\&phi;}{r}_{\mathrm{mc}}\left(1\right)}---\left(15\right)$

如果决定了各回波信号的配置区域和扫描轨迹，则时序生成部720追加每个扇形叶片320（b）的旋转角度θ（b），而决定总叶片数B_fan部分的各拍摄的倾斜磁场形状（步骤S1204）。

在此，时序生成部720首先使用各扇形叶片320（b）的旋转角度θ（b），生成以下的公式（16）所示的旋转矩阵Rot（b）。

[公式16]

$\mathrm{Rot}\left(b\right)=\left[\begin{array}{cc}\cos \mathrm{\&theta;}\left(b\right)& -\sin \mathrm{\&theta;}\left(b\right)\\ \sin \mathrm{\&theta;}\left(b\right)& \cos \mathrm{\&theta;}\left(b\right)\end{array}\right]---\left(16\right)$

另外，将其适用于各回波信号的扫描轨迹，计算与各扇形叶片320（b）一致地旋转后的扫描轨迹，由此决定2轴（Gx，Gy）方向的倾斜磁场脉冲223和224的波形。

此外，通过在一个回波信号的读出中一边切换频率编码方向一边在相位编码方向上施加尖峰信号来实现上述钟摆状轨迹520。

通过以上的步骤，本实施方式的时序生成部720使用由操作者设定的摄像参数，决定测量的倾斜磁场波形，使得在各扇形叶片320内的k空间的低空间频率区域中配置希望的对比度的回波信号。

本实施方式的测量控制部730根据按照具有上述倾斜磁场波形的摄像时序而从信号处理系统70输出的指示，使各部件动作。即，使各扇形叶片320（b）分别从kx轴旋转旋转角度θ（b）而进行测量，这样重复总叶片数B_fan次。这时，在各扇形叶片320（b）内，进行控制使得在k空间的低空间频率区域中配置希望的对比度的回波信号。

接着，说明上述步骤S1104的本实施方式的图像重构部740的伴随着体动修正的图像重构处理。图11（A）是本实施方式的图像重构处理的处理流程。

此外，本实施方式的图像重构部740如图2所示，具备体动修正部750，其使用各扇形叶片320（b）的体动检测区域323（b）的数据（体动检测用数据），进行叶片间体动修正。叶片间体动修正的处理为，对各扇形叶片320（b）的体动检测用数据进行比较，检测作为被检测体10的体动量的并进量和旋转量的至少一方，为了抵消这些，使各扇形叶片320（b）的叶片数据并进、旋转而进行修正。

图像重构部740检测各扇形叶片320（b）中的被检测体10的体动量，因此只使用各扇形叶片320（b）的体动检测区域323内的k空间数据，分别重构图像（步骤S1301）。这时，图像重构部740通过将填充在各扇形叶片320（b）的体动检测区域323中的k空间数据在k空间的正交坐标系的格子点上进行网格化，并进行二维傅立叶变换而重构图像。此外，重构的图像能够计算体动量即可，因此可以是低空间分辨率的图像。将根据各扇形叶片320（b）的体动检测区域323（b）的k空间数据重构出的图像称为体动检测图像（b）。

接着，图像重构部740进行体动量计算处理，即在体动修正部750计算作为各体动检测图像（b）的相对于基准图像的体动量的并进移动量和旋转移动量（步骤S1302）。此外，基准图像例如为根据第一个扇形叶片320（1）的体动检测区域323（1）重构出的体动检测图像。计算出的每个体动检测图像（b）的并进移动量和旋转移动量分别为得到体动检测图像的扇形叶片320（b）的并进移动量Δd（b）和旋转移动量Δq（b）。

接着，图像重构部740使用各扇形叶片320（b）的扇形区域322（b）的k空间数据，用希望的分辨率来重构图像（步骤S1303）。这时，图像重构部740通过将填充在各扇形叶片320（b）中的k空间数据在k空间的正交坐标系的格子点上进行网格化，并进行二维傅立叶变换来重构图像。将根据各扇形叶片320（b）重构出的图像称为图像（b）。

接着，图像重构部740进行体动修正处理，使每个图像（b）并进和旋转各个扇形区域322（b）的并进移动量Δd（b）和旋转移动量Δq（b）（步骤S1304）。

图像重构部740对修正后的各图像（b）进行合成后生成重构图像（步骤S1305）。进行合成后生成重构图像的处理首先对修正后的各图像（b）进行二维傅立叶变换，设为叶片k空间数据。对各扇形叶片重复进行该操作，最后对全部扇形叶片进行网格化，并进行二维傅立叶变换，生成重构图像。

本实施方式的图像重构部740通过以上的步骤一边进行体动修正一边重构图像。

如以上说明的那样，根据本实施方式，由于是非正交系测量，所以能够减少伪影。进而，只将希望的对比度的信号配置在k空间的决定对比度的低空间频率区域中。另外，将实效TE以外的回波信号配置在低空间频率区域以外。因此，不混入其他对比度的信号，就能够高精度地得到希望的对比度的图像。

另外，根据本实施方式，各扇形叶片在k空间的低空间频率区域中具有共通的区域。另外，使用该共通的区域的数据，计算每个扇形叶片的相对于基准叶片的体动量，进行修正使得抵消它。因此，根据本实施方式，也进行叶片间的体动的修正。因此，也能够降低因叶片间的体动造成的伪影。

因此，根据本实施方式，能够伪影少地高速地得到具有希望的对比度的高质量的图像。

另外，用上述公式（4）表示填充本实施方式的半径为R（=F/2）的圆状k空间所需要的总叶片数B_fan。另一方面，在现有的混合径向法的长方形叶片310的情况下，通过公式（17）求出填充半径为R（=F/2）的圆状k空间所需要的总叶片数B_rec。

[公式17]

$B_{\mathrm{rec}}=\frac{\mathrm{\&pi;}}{2}\frac{F}{P}---\left(17\right)$

此外，用公式（18）表示k空间全体的采样面积N_rec。

[公式18]

$N_{\mathrm{rec}}=\mathrm{FP}\&times;B_{\mathrm{rec}}=\frac{\mathrm{\&pi;}}{2}{F}^2---\left(18\right)$

这样，如果对公式（4）和公式（17）进行比较，则使用扇形叶片320的情况与使用长方形叶片310的情况相比，扫描相同k空间的面积（半径R的圆状k空间）所需要的叶片数只需1/2，能够进一步缩短测量时间。因此，根据本实施方式，能够提高k空间的填充效率，进一步缩短测量时间。

进而，在基于扇形叶片的测量的情况下，通过填充椭圆状的k空间，能够进行在长方形叶片中无法实现的矩形视野测量。用以下的公式（19）表示该情况下的采样面积N_fan。

[公式19]

$N_{\mathrm{fan}}=\mathrm{\&pi;R}{R}^{\&prime;}\&le;\frac{\mathrm{\&pi;}}{4}{F}^2---\left(19\right)$

用以下的公式（20）表示该情况下的扇形叶片320（b）的总叶片数B_fan。

[公式20]

$B_{\mathrm{fan}}=\frac{N_{\mathrm{fan}}}{S}=\frac{N_{\mathrm{fan}}}{\mathrm{FP}}=\frac{\mathrm{\&pi;R}{R}^{\&prime;}}{\mathrm{FP}}\&le;\frac{\mathrm{\&pi;}}{4}\frac{F}{P}---\left(20\right)$

如果对公式（17）和公式（19）进行比较，则在使用了扇形叶片的测量中，即使组合了矩形视野测量，也能够减少叶片数，能够缩短测量时间。

此外，在本实施方式，举例说明与FSE时序组合的情况。

但是，本实施方式通过将每个拍摄的回波信号的扫描区域设为扇形叶片，如上述那样决定扫描顺序、扫描轨迹，从而得到上述效果。因此，在本实施方式中组合的脉冲时序只要是在施加一个激励脉冲后在TR期间取得多个回波信号的脉冲时序即可，能够与时序种类、对比度无关地进行应用。

另外，在上述实施方式中，在由测量控制部730测量全部k空间后，由图像重构部740一边进行体动修正一边重构图像，但并不限于此。例如，构成也可以为，如果将根据第一个扇形叶片320（1）的体动检测区域323（1）重构出的图像设为基准图像，则针对第2个以后的扇形叶片320（b），在每次测量各扇形叶片320（b）时，一边进行体动修正一边重构图像。

另外，在上述实施方式中，图像重构部740具备体动修正部750，在图像重构时使用重构图像进行体动修正，但体动修正并不限于此。例如也可以构成为使用图像重构前的体动检测区域323的k空间数据计算体动量，对扇形区域322的k空间数据进行体动修正。图11（b）表示该情况下进行体动修正的图像重构处理的处理流程。

体动修正部750使用体动检测区域323（b）的k空间数据，计算旋转移动量和并进移动量作为该扇形叶片320（b）的k空间数据的体动量（步骤S1401）。这时，作为成为基准的扇形叶片320（b）的体动检测区域323（b）的k空间数据的旋转角度来检测旋转移动量。另外，作为与成为基准的k空间数据的相位差来检测并进移动量。另外，体动修正部750对各k空间数据进行修正使得抵消所得到的体动量（步骤S1402）。

然后，图像重构部740使用修正后的k空间数据进行网格化处理而重构图像（步骤S1403）。

另外，规定体动检测区域323的半径C也可以根据测量部位的运动的特性等进行变化。进而，也可以构成为与体动修正部750计算出的体动量相应地进行变化。即，在体动量大时，为了提高精度而增大体动检测区域323，在体动量小时减小体动检测区域323。能够任意地确定称为使半径C变化的基准的体动量、半径C的变化量。

进而，也可以构成为在每个TR中、即在每次测量第二个扇形叶片以后的一个扇形叶片320（b）时，计算体动量，如上述那样决定半径C。

此外，各扇形叶片320（b）的体动检测区域323（b）内的扫描轨迹也可以如图12所示那样为螺旋状（螺旋状扫描）。用以下的公式（21）表示这时的体动检测区域323内的扫描轨迹（螺旋状轨迹523）。

[公式21]

$\left\{\begin{array}{c}\mathrm{kx}\left(b\right)=R\left(t\right)\cos t\\ \mathrm{ky}\left(b\right)=R\left(t\right)\sin t\end{array}\right.---\left(21\right)$

在此，t是采样数（0≤t≤T），R（t）是相对于k空间原点的距离（0≤R（t）≤C）。

另外，体动检测区域323的形状并不限于圆形。例如，也可以是由与kx轴平行的2边和与ky轴平行的2边组成的长方形。图13（A）表示体动检测区域323是长方形的情况下的体动检测区域323的形状。在此，将体动检测区域323（b）设为1边为X、另一边为Y的像素的长方形。另外，图13（B）和图13（C）表示体动检测区域323是长方形的情况下的扫描轨迹。如本图所示，在体动检测区域323的形状是长方形的情况下，扫描轨迹为蛇形的轨迹（蛇形轨迹524）。

此外，各扇形叶片320（b）中的体动检测区域323（b）根据叶片旋转角度θ（b）如图13（B）和图13（C）所示那样进行旋转。

如图13（A）所示，在体动检测区域323是长方形的情况下，根据蛇形轨迹524在各扇形叶片320（b）中对同一格子点上进行采样。

因此，在体动检测时容易进行扇形叶片320（b）之间的比较。

也可以构成为与体动检测区域323的形状是圆形时同样，体动检测区域323是长方形时的各边的像素数X、Y根据测量部位的运动的大小而进行变化。

另外，在上述实施方式中，作为回波信号的配置顺序，举例说明了中心配置，但回波信号的配置顺序并不限于此。决定将希望的TE（实效TE）时的回波信号配置在最接近k空间的原点的分割区域321（1）中。通过这样的结构，能够得到希望的对比度的图像。

例如，也可以是按照从外周侧开始按顺序取得数据的非中心的数据配置。如果考虑在TR期间取得的回波信号的取得顺序e（1≤e≤E），则分割区域编号n（1≤n≤N）在中心配置中为n=e，但在非中心配置中为n=E-e+1。

另外，也可以适用回波偏移。在该情况下，时序生成部720通过上述同样的方法，决定分割区域和扫描轨迹，在步骤S1204决定各拍摄的倾斜磁场形状时，考虑回波配置顺序来决定。

进而，本实施方式也可以适用多对比度测量。多对比度测量是在一次测量中同时取得2个以上不同的对比度的图像的测量方法。例如，通过在一次脉冲时序中取得2个以上的回波信号，将它们填充到2个以上的k空间中来实现。

在TR长且同时取得质子密度强调图像（PDW）和T2强调图像（T2W）时经常使用多对比度测量。PDW使用短的实效TE，因此回波信号配置使用中心配置。另外，T2W使用比较长的TE，因此使用非中心或回波偏移。此外，根据T2W的实效TE自动计算使用回波偏移的情况下的回波偏移数。

这样，在将本测量适用于多对比度测量的情况下，使用与各图像的TE对应的最适合的回波信号配置。

另外，也可以将本实施方式与扫描不同形状的叶片内的测量进行组合。

例如，在需要提高时间方向的分辨率的动态测量的情况下，需要提高k空间的低频率部分的数据取得频度。

因此，首先使用本实施方式的扇形叶片320，进行成为基准的基础测量，取得k空间的半径R的圆的内部区域全体的数据。然后，在每个预定的时间重复进行基于在以k空间的原点为中心的同心圆状的圆周进行分割的同心圆状区域即同心圆叶片的测量。这时，对于一次取得的多个同心圆叶片，每次只对k空间的低空间频率区域进行测量，对于高空间频率区域，每次测量不同的区域。另外，在各次不足的区域使用在基础测量中得到的k空间的数据。

通过这样的结构，能够在每个预定的时间得到只更新了k空间的低频率的数据的k空间数据。另外，在初次以外，一次的测量区域比全部区域小，因此能够缩短测量时间。因此，在相同的时间内，能够增加重复次数，还能够提高动态测量的时间方向的分辨率。

<第二实施方式>

接着，说明适用本发明的第二实施方式。在第一实施方式中，针对全部扇形叶片，在k空间的低空间频率区域中设置共通的体动检测区域。在本实施方式中，使各扇形叶片在圆周方向上重叠而进行测量，将该重叠区域设为体动检测区域。

本实施方式的MRI装置基本上与第一实施方式相同。另外，信号处理系统70所实现的各功能结构也与第一实施方式相同，基于这些功能的摄影处理也相同。

但是，如上述那样，在每个TR中测量的叶片的形状不同。因此，通过决定这样的叶片的形状的时序生成部进行的叶片参数决定处理不同。以下，针对本实施方式，着眼于与第一实施方式不同的结构的时序生成部720进行的叶片参数决定处理进行说明。

首先，说明本实施方式的一次TR中的测量区域即叶片的形状。

在本实施方式中，在一次TR中测量扇形形状的叶片（扇形叶片）330。

本实施方式的扇形叶片330如图14（A）所示，具有与第一实施方式的扇形叶片320的扇形区域322相同的形状。

如图14（B）所示，本实施方式的扇形叶片330是以k空间的原点为中心的半径R的圆的由2个半径R和两个半径R之间的圆弧围住的区域。扇形叶片330的中心角为φ（rad）。第b个扇形叶片330（b）的ky（b）轴在对扇形叶片330（b）的中心角φ进行二分割的线（中心线）上。kx（b）轴在与ky（b）轴垂直的方向。另外，扇形叶片330（b）的旋转角度θ（b）为kx（b）轴与kx轴所成的角。

在第一实施方式中，如图15（A）所示，在摄像时序中，各扇形叶片320的扇形区域322与相邻的扇形叶片320的扇形区域322不重叠，决定扇形叶片320（b）的中心角φ（rad）和旋转角度θ（b）从而对k空间的半径R的区域全体进行测量。

另一方面，在本实施方式中，如图15（B）所示，各扇形叶片330（b）具有在圆周方向上与相邻的扇形叶片330重叠的区域（重叠区域331），并且决定其中心角φ（rad）和旋转角度θ（b）从而对k空间的半径R的区域全体进行测量。

按照相对于叶片中心角φ的比例来规定重叠区域331的中心角α。

例如预先将比例设定为10%等。

图16表示本实施方式的时序生成部720进行的叶片参数决定处理的流程。

本实施方式的时序生成部720在决定叶片参数之前，先受理该重叠区域331的中心角α的设定（步骤S2101）。

然后，时序生成部720使用摄像参数作为叶片参数，计算并决定扇形区域322的半径R、中心角φ和总叶片数B_fan（叶片参数决定：步骤S2102），从而实现具有设定的中心角α的重叠区域331。

扇形叶片330的半径R与第一实施方式相同，为F/2。F是长方形叶片310的频率编码方向的采样数。

首先，决定没有重叠区域情况下的扇形叶片330的中心角φ和总叶片数B_fan。这些与第一实施方式的总叶片数B_fan和扇形区域332的中心角α的计算方法同样，通过上述公式（4）和公式（2）来计算。

通过加大上述没有重叠区域地配置扇形叶片330的情况下的总叶片数B_fan和叶片中心角φ的任意一方，来实现中心角α的重叠区域331。

在增加总叶片数B_fan的情况下，使用通过上述公式（2）求出的中心角φ、重叠区域331的中心角α，通过以下的公式（22）得到增加后的总叶片数B_fan’。

[公式22]

${B_{\mathrm{fan}}}^{,}=\frac{2\mathrm{\&pi;}}{\mathrm{\&phi;}-\mathrm{\&alpha;}}---\left(22\right)$

这时，通过以下的公式（23）求出各扇形叶片330（b）的旋转角度θ（b）。

[公式23]

θ(b)＝(φ-α)(b-1)      (23)

另一方面，在加大叶片中心角φ的情况下，用以下的公式（24）表示加大后的叶片中心角φ`。

[公式24]

φ′＝φ+α      (24)

这时，通过以下的公式（25）求出各扇形叶片330（b）的旋转角度θ（b）。

[公式25]

θ(b)＝φ′(b-1)      (25)

本实施方式的时序生成部720如以上那样决定总叶片数、中心角、旋转角（（B_fan’、φ、θ（b））或（B_fan、φ’、θ（b）））。

时序生成部720与第一实施方式的扇形区域322同样，在半径方向上对本实施方式的扇形叶片330进行分割，决定各回波信号在扇形叶片330内的配置区域（步骤S2103）。决定方法与第一实施方式相同。

另外，时序生成部720决定各分割区域321内的扫描轨迹（步骤S2104）。在本实施方式中，不具备体动检测区域323，因此本实施方式的扫描轨迹与第一实施方式的考虑体动检测区域323之前的扇形区域322内的钟摆状轨迹520相同。

另外，如果决定了各回波信号的配置区域和扫描轨迹，时序生成部720则追加每个扇形叶片330（b）的旋转角θ（b），而通过与第一实施方式相同的方法，决定B_fan或B_fan’次的各拍摄的倾斜磁场形状（步骤S2105）。

接着，说明本实施方式的图像重构部740进行体动修正的图像重构处理。本实施方式的图像重构处理基本上与第一实施方式的图像重构处理相同。

即，在本实施方式中，图像重构部740为了检测扇形叶片330（b）中的被检测体的体动量，只使用各重叠区域331内的k空间数据，分别重构图像。这时，图像重构部740通过将填充在各重叠区域331中的k空间数据在k空间的正交坐标系的格子点上进行网格化，并进行二维傅立叶变换，由此重构图像。此外，重构的图像能够计算体动量即可，因此可以是低空间分辨率的图像。将通过各扇形叶片330（b）的重叠区域331的k空间数据重构出的图像称为体动检测图像b。

接着，图像重构部740进行体动量计算处理，即在体动修正部750计算并进移动量和旋转移动量作为各体动检测图像（b）的相对于基准图像的体动量。此外，基准图像例如为由第一个扇形叶片330（1）和相邻的扇形叶片330（b）的重叠区域331重构出的体动检测图像（1）。

在此，在本实施方式中，重叠区域331是只在相邻的扇形叶片330之间共通的区域，因此从与得到基准图像的扇形叶片330相邻的扇形叶片330开始按顺序地计算相对于相邻的扇形叶片的体动量，与相对于该相邻的扇形叶片的基准图像的体动量相加，由此得到相对于基准图像的体动量。

即，如果设第b个扇形叶片330（b）和第b-1个扇形叶片330（b-1）之间的体动量为ΔD_b（ΔDb=（Δx、Δy、Δz、Δα、Δβ、Δγ）；x、y、z是并进量，α、β、γ是旋转量），则第（b+1）个扇形叶片330（b+1）中的外表的体动量ΔD_b+1表示为D_b+1=D_b+ΔD_b。因此，根据以下的公式（26），通过此前测量出的扇形叶片330的体动量来求出全部扇形叶片330的外表的体动量。

[公式26]

$\begin{array}{c}{D}_{b+1}=D_b+{\mathrm{\&Delta;D}}_b\\ =(D_{b-1}+{\mathrm{\&Delta;D}}_{b-1})+{\mathrm{\&Delta;D}}_b\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ =\underset{m=1}{\overset{b}{\mathrm{\&Sigma;}}}{\mathrm{\&Delta;D}}_m\end{array}---\left(26\right)$

在此，D_i=0。

此后的图像重构部740进行的使用各扇形叶片330（b）的体动量的体动修正与第一实施方式相同。即，使用各扇形叶片330（b）的k空间数据重构图像，使用各扇形叶片的体动量进行体动修正，得到重构图像。或者，在各扇形叶片330（b）的k空间数据上进行体动修正，得到重构图像。

如以上说明的那样，根据本实施方式，与第一实施方式同样，能够伪影少地高速地得到具有希望的对比度的高质量的图像。

此外，在本实施方式中，与第一实施方式同样也可以构成为，在每次测量各扇形叶片330时进行伴随着体动修正的图像重构。

另外，在本实施方式中，与第一实施方式同样也可以构成为，体动修正部750使用重叠区域331的k空间数据计算体动量，进行体动修正。在该情况下，与第一实施方式同样，图像重构部740根据体动修正后的k空间数据重构图像。

另外，在本实施方式中，与第一实施方式的体动检测区域323同样也可以构成为，根据计算出的体动量的大小来调整重叠区域331的中心角α。即，在体动量大时，为了提高精度而变更中心角α使得重叠区域331变大，在体动量小时，进行变更使得中心角α减小。预先确定体动量的大小的判别基准值、中心角α的变化量。

另外，在也可以构成为在每个TR中计算体动量，根据其结果使下一个TR的重叠区域331的中心角α发生变化。

另外，本实施方式的扇形叶片330也可以还具备第一实施方式的体动检测区域323。

在该情况下，图像重构部740的体动修正部750使用体动检测区域323和重叠区域331的两个区域的数据来检测体动量。

此外，在上述各实施方式中，构成为MRI装置100的信号处理系统70根据摄像条件计算出实现基于扇形叶片的k空间扫描的倾斜磁场波形，但并不限于此。也可以构成为在与MRI装置100可进行数据的收发、与MRI装置独立的信息处理装置上计算倾斜磁场波形。

另外，在上述各实施方式中，在摄像时，时序生成部720在每次设定摄像参数时计算扇形叶片的参数，计算倾斜磁场形状，但并不限于此。例如，也可以构成为预先对每个有可能使用的摄像参数计算倾斜磁场形状，与摄像参数对应地作为数据库保存在存储装置72等中。在该情况下，如果在摄影时受理摄影参数，则在上述步骤S1102中，时序生成部720参照数据库，抽出与受理的摄影参数对应存储的倾斜磁场形状，生成摄影时序。

另外，在上述各实施方式中，将得到重构图像用的数据的区域设为以k空间的原点为中心的半径R的圆的由2个半径R和两个半径R之间的圆弧围住的扇形的区域，但并不限于此。只要是由以k空间的原点为起始点的2个线段和将该2个线段的另一个端点连接起来的线围住的区域即可。

符号说明

10：被检测体；20：静磁场产生系统；30：倾斜磁场产生系统；31：倾斜磁场线圈；32：倾斜磁场电源；40：时序产生器；50：发送系统；51：发送线圈；52：合成器；53：调制器；54：高频放大器；60：接收系统；61：接收线圈；62：信号放大器；63：正交相位检波器；64：A/D变频器；70：信号处理系统；71：CPU；72：存储装置；73：外部存储装置；74：显示装置；75：输入装置；100：MRI装置；210：摄像时序；211：激励RF脉冲；212：再收敛RF脉冲；213：倾斜磁场脉冲；220：摄像时序；223：倾斜磁场脉冲；310：长方形叶片；320：扇形叶片；321：分割区域；322：扇形区域；323：体动检测区域；330：扇形叶片；331：重叠区域；332：扇形区域；510：直线状轨迹；520：钟摆状轨迹；521：圆弧状轨迹；522：直线状轨迹；523：螺旋状轨迹；524：蛇形轨迹；710：受理部；720：时序生成部；730：测量控制部；740：图像重构部；750：体动修正部。

Claims (15)

1.一种磁共振成像装置，具备：静磁场产生单元，其产生静磁场；磁场施加单元，其向配置在上述静磁场中的被检测体的希望的摄像区域施加倾斜磁场和高频磁场；检测单元，其检测来自上述希望的摄像区域的回波信号；测量控制单元，其控制上述磁场施加单元和上述检测单元，测量上述回波信号从而取得k空间内的预定的区域的数据；图像重构单元，其使用上述k空间的数据来重构上述摄像区域的图像，该磁共振成像装置的特征在于，

上述测量控制单元进行控制，使得在每次上述单位测量时使上述单位区域以上述k空间的原点为中心旋转预先决定的旋转角度，从而重复在1个TR内测量单位区域的多个数据的单位测量，

在上述单位测量中，将希望的对比度的回波信号配置在上述k空间的低空间频率区域中，

上述单位区域具备通过以上述k空间的原点为起始点的2个线段和将该2个线段的另一个端点连接的线而进行包围的第一区域、和在相邻的上述单位区域之间重叠的区域，

上述图像重构单元具备：体动修正单元，其在使用上述重叠的区域的数据重构上述图像时，进行体动修正。

2.根据权利要求1所述的磁共振成像装置，其特征在于，

上述2个线段分别是以上述k空间的原点为中心的圆的半径，

将上述端点连接起来的线是上述圆的圆弧。

3.根据权利要求2所述的磁共振成像装置，其特征在于，

上述重叠的区域是设定在k空间的低空间频率区域中的与全部单位区域形状相同的区域。

4.根据权利要求3所述的磁共振成像装置，其特征在于，

决定上述各单位区域的旋转角，使得在以上述k空间的原点为中心的圆的圆周方向上不重复地配置上述第一区域。

5.根据权利要求1所述的磁共振成像装置，其特征在于，

决定上述各单位区域的旋转角，使得将上述第一区域配置为在以上述k空间的原点为中心的圆的圆周方向上具有预先确定的中心角的重叠区域。

6.根据权利要求2所述的磁共振成像装置，其特征在于，

具备：

摄像条件受理单元，其从操作者受理摄像条件；

摄像时序生成单元，其根据所受理的上述摄像条件来生成摄像时序，

上述测量控制单元根据上述摄像时序进行上述控制，

上述摄像时序生成单元具备：

参数决定单元，其决定以上述k空间的原点为中心的圆的半径、上述2个线段所成的角度的中心角、以及上述旋转角分别不同的单位区域的个数的总单位区域数；

区域决定单元，其决定每个上述单位测量的各回波信号的配置区域；

扫描轨迹决定单元，其决定所决定的上述各配置区域内的扫描轨迹；

波形决定单元，其根据各单位区域与上述k空间的预先确定的坐标轴所成的角度的旋转角度，决定每个上述单位测量的扫描轨迹，决定上述摄像时序的上述倾斜磁场波形。

7.根据权利要求6所述的磁共振成像装置，其特征在于，

以上述k空间的原点为中心的圆的半径由频率编码数来决定，

上述中心角由上述圆的直径和该单位区域内的总数据数来决定，

上述总单位区域数在将2π除以上述中心角所得的值以上。

8.根据权利要求6所述的磁共振成像装置，其特征在于，

上述各配置区域内的扫描轨迹是钟摆状。

9.根据权利要求1所述的磁共振成像装置，其特征在于，

上述体动修正单元在重构上述重叠的区域的数据所得的图像上计算各单位区域的图像的体动量，修正根据各单位区域的数据重构出的图像使得抵消计算出的体动量。

10.根据权利要求1所述的磁共振成像装置，其特征在于，

上述体动修正单元使用上述重叠的区域的数据，在k空间上计算各单位区域的体动量，对各单位区域的数据进行修正使得抵消计算出的体动量，

上述图像重构单元根据上述修正后的数据来重构图像。

11.根据权利要求1所述的磁共振成像装置，其特征在于，

根据上述体动修正单元计算出的体动量来决定上述重叠的区域。

12.根据权利要求3所述的磁共振成像装置，其特征在于，

上述重叠的区域是以k空间的原点为中心的圆内的区域。

13.根据权利要求3所述的磁共振成像装置，其特征在于，

上述重叠的区域是以k空间的原点为中心，而各边与kx轴和ky轴的任意一个平行的长方形的区域。

14.一种重构图像取得方法，是磁共振成像装置的重构图像取得方法，其特征在于，包括：

取得单位区域的数据的单位测量步骤；

使上述单位区域以上述原点为中心旋转预先确定的旋转角而重复进行上述单位测量步骤的重复测量步骤；

根据在上述重复测量步骤中得到的数据重构图像的图像重构步骤，

上述单位区域具备：通过以上述k空间的原点为起始点的2个线段和将该2个线段的另一个端点连接的线而进行包围的第一区域、和在相邻的上述单位区域之间重叠的区域，

上述图像重构步骤包括：在根据每个上述单位区域的重叠的区域的数据重构出的图像上计算该每个单位区域的体动量的体动量计算步骤；

使用该单位区域的上述体动量来修正根据每个上述单位区域的第一区域的数据重构出的图像的体动修正步骤；

对每个上述单位区域的修正后的图像进行合成的图像合成步骤。

15.一种重构图像取得方法，是磁共振成像装置的重构图像取得方法，其特征在于，包括：

取得单位区域的数据的单位测量步骤；

使上述单位区域以上述原点为中心旋转预先确定的旋转角而重复进行上述单位测量步骤的重复测量步骤；

根据在上述重复测量步骤中得到的数据重构图像的图像重构步骤，

上述单位区域具备：通过以上述k空间的原点为起始点的2个线段和将该2个线段的另一个端点连接的线而进行包围的第一区域、和在相邻的上述单位区域之间重叠的区域，

上述图像重构步骤包括：

在每个上述单位区域的重叠的区域的数据上计算该每个单位区域的体动量的体动量计算步骤；

使用该单位区域的上述体动量来修正每个上述单位区域的第一区域的数据的体动修正步骤；

根据每个上述单位区域的修正后的数据重构图像的图像重构步骤。

Images