磁共振成像装置以及翻转角决定方法
磁共振成像装置以及翻转角决定方法
【专利摘要】为了在多回波序列中使考虑了信号修正的图像的SNR最大，决定多回波摄像序列的多个再收敛高频磁场脉冲的翻转角，使用反映所取得的多个回波信号的信号修正后的图像的SNR的指标，依照预先确定的最优化法，使确定各再收敛RF脉冲的翻转角的信息变化，重复计算指标，决定图像的SNR最大的翻转角。
【专利说明】磁共振成像装置以及翻转角决定方法
【技术领域】
[0001]本发明涉及ー种测定来自被检测体的氢、磷等的核磁共振(以下称为“NMR”)信号来对核密度分布、迟豫时间分布等进行图像化的核磁共振成像(以下称为“MRI ”)技术，特别涉及ー种使用多回波摄像的VRFA (可变再调焦翻转角)的成像技木。
【背景技术】
[0002]MRI装置是测定构成被检测体特别是人体的组织的原子核自旋所产生的NMR信号(回波信号)，二维或三维地对其头部、腹部、四肢等的形态或功能进行图像化(摄像)的装置。在摄像中，对NMR信号赋予根据倾斜磁场而不同的相位编码，并且进行频率编码，作为时间序列数据来测量。对测量到的回波信号进行二维或三维傅立叶变换，来重新构成图像。
[0003]依照预先确定的摄像序列測量回波信号。该摄像序列有各种各样，但其中，有以下的序列,在施加一次激励高频磁场(RF)脉冲后，一边改变相位编码量一边多次施加再收敛高频磁场(RF)脉冲来收集回波信号，谋求高速化。将这样的摄像称为多回波摄像，将在各再收敛RF脉冲后测量的回波信号群称为回波链。另外，将在施加一次的RF脉冲后得到的回波信号的个数称为回波链数(ETL:Echo Train Length)。
[0004]在多回波摄像中，有时改变再收敛RF脉冲的FA(翻转角:Flip Angle)。将这样的摄像序列称为VRFA多回波序列。在VRFA多回波序列中，为了达到预定的目的而决定各再收敛RF脉冲的FA。
[0005]例如，有决定FA来提高`图像的分辨率的技术(例如參照专利文献I)。在专利文献I所公开的技术中，决定FA以使回波链内的各回波信号的信号強度恒定。这是因为如果各回波信号的信号強度恒定，则能够抑制信号衰减，即使使ETL变长，图像也不会模糊(分辨率提高)。这时的FA典型的是在一旦成为小的值(FAmin)后会向最后的FA平滑地増大。在此，将k空间中心的FA定义为FAcenter,将最后的FA定义为FAmax。
[0006]此外，在本说明书中，假设回波信号的信号强度忽略由于相位编码、频率编码造成的强度变化。另外，对于信号強度的值，将激励的瞬间的假想信号强度标准化为I。
[0007]但是，画质不只由分辨率决定，还由对运动的耐性、SNR (信噪比)等各种要素決定。特别是SNR，在决定画质的好坏上是主要的。在专利文献I所公开的技术中，以使信号強度恒定为目的，因此所得到的图像的SNR并不一定最优。另外，恒定的值的决定方法是任意的，并不知道在怎样的值时成为怎样的画质。
[0008]作为提高该SNR的技术，有决定FA以便提高配置在k空间中心的回波信号的信号強度的技术(例如參照非专利文献I)。这是因为一般配置在k空间的中心的回波信号决定了 SNR。在非专利文献I中，决定FA以使在指定了 FAmin、FAcenter> FAmax之后信号平滑地变化。这是因为通过増大FAmin能够提高图像对运动的耐性，通过提高FAcenter能够提高不考虑信号修正时的SNR。
[0009]另外，有以下的技木，即在将2种数据相加(合成)计算2种对比度的图像时，决定FA使得各数据的信号强度统一(例如參照非专利文献2)。在非专利文献2中，作为2种数据，取得CSF (脑脊液)和脑实质的信号強度的平衡不同的数据。将这些数据相加，得到T2W和FLAIR的图像。这时，决定FA使得2种数据的CSF的信号強度在k空间中心的区域中统一。通过在取得2种数据的时刻使信号强度统一，在相加时不需要进行用于统一信号強度的加权，提闻了 SNR。
[0010]此外，作为根据所決定的FA计算信号強度的变化的方法，已知EPG (扩展相位图)(例如參照非专利文献3)。另外，相反，作为根据信号強度的变化计算FA的方法，已知ProspectiveEPG (例如參照非专利文献4)。
[0011]现有技术文献
[0012]非专利文献
[0013]非专利文献1:Reed F.Busse, Anja C.S.Brau, Anthony Vu, CharlesR.Michelich, Ersin Bayram, Richard Kijowski, Scott B.Reeder, and HowardA.Rowley."Effects of Refocusing Flip Angle Modulation and View Ordering in3DFast Spin Ech0."Magnetic Resonance in Medicine60:640649(2008)
[0014]非专利文献2:H.Lee, J-S.Seo, and J.Park.^SNR-Optimized AcceleratedPhase-Sensitive Dual-Acquisition Single_Slab3D Turbo Spin Echo Imaging〃Proc.1ntl.Soc.Mag.Reson.Med.19:386 (2011)
[0015]非专利文献 3:Juergen Hennig."Multieco Imaging Sequences with LowRefocusing Flip Angles^Journal of Magnetic Resonance78:397407(1988)
[0016]非专利文献4: J uergen Hennig, Matthias Weigel, and KlausSheffler."Calculation of Flip Anbgles for Echo Trains With PredefinedAmplitudes With the Extended Phase Graph (EPG)-Algorithm:Principles andApplications to Hyperecho and TRAPS Sequences."Magnetic Resonance inMedicine51:6880(2004)
【发明内容】
[0017]发明要解决的问题
[0018]一般在MRI中，以提高画质等为目的，进行各种信号修正。例如在多回波摄像中，有时修正信号的衰减(T2修正)来抑制图像的模糊。如果进行了这样的信号修正，则与回波信号強度对应地噪声提高，因此全部的回波信号的信号强度对SNR有影响。因此，在进行信号修正的情况下，在非专利文献I所公开的方法中无法进行充分的SNR的控制。
[0019]另外，在非专利文献2所公开的技术中，也只是在k空间中心的区域中统ー 2种图像的CSF的信号強度，一般没有考虑到回波链整体的信号強度来使SNR最大。因此，在进行信号修正的情况下，在非专利文献2所公开的方法中也无法进行SNR的控制。
[0020]本发明就是鉴于上述情况而提出的，其目的在于，提供ー种在VRFA多回波序列中使考虑到信号修正的SNR最大的技术。
[0021]解决问题的方案
[0022]本发明用于决定磁共振成像装置中的多回波摄像序列的、多个再收敛高频磁场脉冲的翻转角，使用反映所取得的多个回波信号的信号修正后的图像的SNR的指标，依照预先确定的最优化法，改变用于确定各再收敛RF脉冲的翻转角的信息，重复计算指标，决定图像的SNR最大的翻转角。
[0023]例如，使用还反映配置在k空间的中心以外的回波信号的信号強度的指标，决定VRFA多回波序列的再收敛RF脉冲的FA，使得SNR大致成为最大。根据不只是反映配置在k空间的中心的回波信号的信号強度，还反映配置在k空间的中心以外的回波信号的信号强度的指标，判别改变再收敛RF脉冲的FA时的SNR的大小。
[0024]发明效果
[0025]根据本发明，在VRFA多回波序列中，能够使考虑到信号修正的SNR最大。 【专利附图】
【附图说明】
[0026]图1是第一实施方式的MRI装置的框图。
[0027]图2是用于说明FSE序列的脉冲序列的说明图。
[0028]图3 Ca)是第一实施方式的序列作成部的功能框图，(b)是第一实施方式的变形例的序列作成部的功能框图。
[0029]图4是第一实施方式的摄像处理的流程图。
[0030]图5 (a)是用于说明第一实施方式的UI画面的说明图，(b)是用于说明第一实施方式的n画面的变形例子的说明图。
[0031]图6是第一实施方式的FA决定处理的流程图。
[0032]图7 (a)、(b)分别是通过第一实施方式的FA决定处理得到的信号强度变化形状和FA变化形状的图表。
[0033]图8 (a)、(b)是用于说明第一实施方式的FA修正部的修正处理的说明图。
[0034]图9 (a)、(b)是用于说明第一实施方式的FA修正部的修正处理的说明图。
[0035]图10是第一实施方式的摄像处理的变形例的流程图。
[0036]图11是第一实施方式的摄像处理的变形例的流程图。
[0037]图12是第一实施方式的摄像处理的变形例的流程图。
[0038]图13是第二实施方式的序列作成部的功能框图。
[0039]图14是第二实施方式的FA决定处理的流程图。
[0040]图15是第三实施方式的序列作成部的功能框图。
[0041]图16是第三实施方式的FA决定处理的流程图。
[0042]图17 (a)、(b)分别是通过第三实施方式的FA决定处理得到的信号强度变化形状和FA变化形状的图表。
[0043]图18是第四实施方式的序列作成部的功能框图。
[0044]图19是第四实施方式的FA决定处理的流程图。
[0045]图20是第四实施方式的信号强度变化形状决定处理的流程图。
[0046]图21是第四实施方式的最小噪声放大率计算处理的流程图。
[0047]图22是第四实施方式的最终FA最小化处理的流程图。
[0048]图23是用于说明第四实施方式的UI画面的ー个例子的说明图。
[0049]图24 (a)、(b)分别是通过第四实施方式的FA决定处理得到的信号强度变化形状和FA变化形状的图表。
[0050]图25是第四实施方式的修正后的FA变化形状的图表。【具体实施方式】
[0051]本发明的用于决定多回波摄像序列的多个再收敛高频磁场脉冲的翻转角的结构的概要如下所述。
[0052]ー种磁共振成像装置，具备产生静磁场的静磁场产生系统、向配置在上述静磁场中的被检测体施加倾斜磁场的倾斜磁场产生系统、发送以预定的翻转角激励上述被检测体磁化的高频磁场脉冲的高频磁场发送系统、接收上述被检测体产生的回波信号的信号接收系统、根据上述信号接收系统接收到的回波信号重新构成图像并且依照摄像序列控制上述倾斜磁场产生系统、上述高频磁场发送系统、上述信号接收系统的动作的控制系统，
[0053]上述控制系统具备:
[0054]接受部，其从用户接受决定上述摄像序列的摄像条件；序列生成部，其依照上述接受部接受的摄像条件和脉冲序列生成上述摄像序列，上述脉冲序列是在施加I个激励高频磁场脉冲后施加多个再收敛高频磁场脉冲的序列，
[0055]上述序列生成部根据上述信号接收系统接收的多个回波信号的信号強度，决定上述各再收敛闻频磁场脉冲的翻转角以使上述图像的SNR最大。
[0056]在上述磁共振成像装置中，上述序列生成部可以使用反映上述多个回波信号的信号修正后的图像的SNR的指标，决定使该图像的SNR最大的上述各再收敛闻频磁场脉冲的翻转角。
[0057]在上述磁共振成像装置中，例如上述序列生成部具备:翻转角决定部，其根据上述摄像条件决定上述各再收敛高频磁场脉冲的翻转角，上述翻转角决定部可以通过依照预先确定的最优化法使确定上述各再收敛高频磁场脉冲的翻转角的信息变化，重复计算上述指标，来决定上述翻转角。
[0058]在上述磁共振成像装置中，例如上述翻转角决定部具备:信号強度计算部，其根据确定上述各再收敛高频磁场脉冲的翻转角的信息，分别计算通过施加与各信息对应的再收敛高频磁场脉冲而得到的回波信号的信号強度，作为上述指标，可以使用根据计算出的上述各回波信号的信号強度的倒数的平方的平均的平方根得到的噪声放大率。可以决定上述各再收敛高频磁场脉冲的翻转角以使其噪声放大率最小。
[0059]在上述磁共振成像装置中，例如确定上述各再收敛高频磁场脉冲的翻转角的信息是再收敛高频磁场脉冲的翻转角。
[0060]在上述磁共振成像装置中，例如确定上述各再收敛高频磁场脉冲的翻转角的信息由比再收敛高频磁场脉冲少的数构成。
[0061]在上述磁共振成像装置中，例如确定上述各再收敛高频磁场脉冲的翻转角的信息可以是依照预定的规则间隔抽取各再收敛高频磁场脉冲的翻转角所得到的中间值列，上述翻转角决定部具备:中间值列决定部，其使上述中间值列变化来决定使上述SNR最大的中间值列；翻转角计算部，其对上述决定的中间值进行内插，计算上述各再收敛高频磁场脉冲的翻转角。
[0062]在上述磁共振成 像装置中，例如确定上述各再收敛高频磁场脉冲的翻转角的信息可以是表示通过施加各再收敛高频磁场脉冲而得到的回波信号的信号強度的变化形状的函数的參数，上述翻转角决定部具备:信号强度变化形状决定部，其使上述參数变化来决定使上述SNR最大的上述參数，决定上述信号強度的变化的形状；翻转角计算部，其根据用上述參数确定的信号強度的变化的形状，计算上述各再收敛高频磁场脉冲的翻转角。
[0063]在上述磁共振成像装置中，例如上述函数可以是渐近地逼近指数函数的函数，上述參数可以是上述指数函数的系数以及时间常数中的至少一方。
[0064]在上述磁共振成像装置中，例如上述函数可以是二重指数函数，上述參数是构成上述二重指数函数的2个指数函数各自的系数和时间常数中的ー个以上。
[0065]在上述磁共振成像装置中，上述二重指数函数也可以是从通过上述激励高频磁场脉冲激励的瞬间的虚拟的信号強度衰减的函数。
[0066]在上述磁共振成像装置中，例如可以预先确定使上述翻转角变化的范围，上述信号强度变化形状决定部使上述參数变化，以使上述翻转角收敛在上述范围内并且不饱和。
[0067]在上述磁共振成像装置中，例如作为上述指标，可以使用根据计算出的上述各回波信号的信号強度的倒数的平方的平均的平方根得到的噪声放大率，上述二重指数函数由上述信号強度的变化形状渐近地逼近的第一指数函数、表示渐近的变化的第二指数函数构成，上述參数是上述第一指数函数的时间常数，上述信号强度变化形状决定部在固定上述噪声放大率来使上述第一指数函数的时间常数变化时把最后的再收敛高频磁场脉冲的上述翻转角最小的上述第一指数函数的时间常数决定为使上述SNR为最大的上述參数。
[0068]在上述磁共振成 像装置中，例如可以预先确定上述翻转角的范围，上述翻转角决定部在计算出的翻转角不存在于上述范围内的情况下，依照预先确定的规则决定翻转角。
[0069]在上述磁共振成像装置中，例如上述最优化法也可以是缩小解的范围的探索法。
[0070]在上述磁共振成像装置中，例如上述序列生成部还可以具备:翻转角修正部，其依照预先确定的规则修正由上述翻转角决定部決定的翻转角。
[0071]在上述磁共振成像装置中，例如上述翻转角修正部也可以进行修正使得由上述翻转角决定部決定的翻转角的一部分或全部收敛在预先确定的可设定范围内。
[0072]在上述磁共振成像装置中，例如上述翻转角修正部也可以在将预先确定的上述再收敛高频磁场脉冲的翻转角修正为预先确定的值后，对该翻转角以外的最小的翻转角以后的翻转角进行修正使其平滑地变化。
[0073]在上述磁共振成像装置中，例如修正为上述预先确定的值的翻转角可以是得到配置在k空间的中心的回波信号的翻转角。
[0074]在上述磁共振成像装置中，例如还可以具备:翻转角存储部，其对每个摄像条件存储作为上述再收敛高频磁场脉冲的翻转角的使上述图像的SNR最大而决定的翻转角，上述序列生成部使用与上述接受部接受的摄像条件对应地存储在上述翻转角存储部中的翻转角来作为上述翻转角。
[0075]在上述磁共振成像装置中，例如上述序列生成部也可以在上述翻转角存储部中没有存储与上述摄像条件对应的翻转角的情况下，使用与最接近的摄像条件对应地存储的翻转角作为上述翻转角。
[0076]在上述磁共振成像装置中，例如上述序列生成部也可以在上述翻转角存储部中没有存储与上述摄像条件对应的翻转角的情况下，使用对与最接近的摄像条件对应地存储的翻转角进行内插而计算出的翻转角作为上述翻转角。
[0077]在上述磁共振成像装置中，例如在上述摄像条件中包含用于确定成为摄像对象的核素的纵向迟豫时间和横向迟豫时间的信息。
[0078]另外，一种翻转角决定方法，其决定在磁共振成像装置中使用的多回波摄像序列的多个再收敛高频磁场脉冲的翻转角，该翻转角决定方法具有:使用反映在上述磁共振成像装置中取得的多个回波信号的信号修正后的图像的SNR的指标，依照预先确定的最优化法，使确定上述各再收敛高频磁场脉冲的翻转角的信息变化，重复计算上述指标，决定上述图像的SNR成为最大的翻转角的翻转角决定步骤。
[0079]以下，详细说明各实施方式。
[0080]〈く第一实施方式》
[0081]说明应用本发明的第一实施方式。以下，在用于说明本发明的实施方式的全部图中，具有相同功能的部分赋予相同符号，省略其重复的说明。
[0082]最初，根据图1说明本实施方式的MRI装置100的概要。图1是表示本实施方式的MRI装置100的整体结构的框图。本实施方式的MRI装置100利用NMR现象得到被检测体的断层图像，具备静磁场产生系统20、倾斜磁场产生系统30、序列产生器40、高频磁场发送系统(发送系统)50、信号接收系统(接收系统)60、控制系统(序列产生器40和信号处理系统70)。
[0083]静磁场产生系统20如果是垂直磁场方式，则在被检测体10的周围的空间中在与其体轴垂直的方向上产生均匀的静磁场，如果是水平磁场方式，则在体轴方向上产生均匀的静磁场，在被检测体10的周围配置永磁体方式、常导方式或超导方式的静磁场产生源(例如磁铁)。
[0084]倾斜磁场产生系统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信号(回波信号)编码各个方向的位置信息。
[0085]序列产生器40控制倾斜磁场产生系统30、发送系统50、接收系统60，使得依照预定的摄像序列重复施加高频磁场脉冲(以下称为“RF”脉冲)和倾斜磁场脉沖。序列产生器40依照来自后述的信号处理系统70具备的CPU71的控制信号而动作，向它们发送收集被检测体10的断层图像的数据所需要的各种指令。
[0086]发送系统50为了使构成被检测体10的生物体组织的原子的原子核自旋产生核磁共振，向被检测体10照射RF脉冲，具备高频振动器(合成器)52、调制器53、高频放大器54、发送侧的高频线圈(发送线圈)51。在基于来自序列产生器40的指令的定时下通过调制器53对从合成器52输出的高频脉冲进行振幅调制，通过高频放大器54对该振幅调制后的高频脉冲进行放大，提供给接近被检测体10配置的发送线圈51，由此向被检测体10照射RF脉冲。
[0087]接收系统60检测通过构成被检测体10的生物体组织的原子核自旋的核磁共振释放的回波信号(NMR信号)，具备接收侧的高频线圈(接收线圈)61、信号放大器62、正交相位检波器63、A/D变换器64。通过从发送线圈51照射的电磁波而感应出的被检测体10的响应的回波信号被接近被检测体10配置的接收线圈61检测，在通过信号放大器62放大后，在基于来自序列产生器40的指令的定时下通过正交相位检波器63分割为正交的两个系统的信号，分别通过A/D变换器64变换为数字量，发送到信号处理系统70。
[0088]信号处理系统70进行各种数据处理、处理结果的显示和保存等，具备CPU71、存储装置72、外部存储装置73、显示装置74、输入装置75。例如使用来自接收系统60的数据，重新构成被检测体10的断层图像。另外，依照序列，向序列产生器40发送控制信号。重新构成出的断层图像被显示在显示装置74上，并且记录到存储装置72或外部存储装置73中。输入装置75用于由操作者输入MRI装置100的各种控制信息、由信号处理系统70进行的处理的控制信息，具备跟踪球或鼠标、以及键盘。该输入装置75接近显示装置74地配置，操作者ー边看着显示装置74，ー边通过输入装置75交互地控制MRI装置100的各种处理。
[0089]此外，在图1中，发送线圈51和倾斜磁场线圈31被设置得在插入被检测体10的静磁场产生系统20的静磁场空间内，如果是垂直磁场方式则与被检测体10相対，如果是水平磁场方式则围住被检测体10。另外，接收线圈61被设置得与被检测体10相对或围住被检测体10。
[0090]现在，在临床普及的MRI装置的摄像对象原子核种类是作为被检测体的主要构成物质的氢原子核(质子)。通过对与质子密度的空间分布、激励状态的迟豫时间的空间分布相关的信息进行图像化，二维或三维地拍摄人体头部、腹部、四肢等的形态或功能。
[0091]根据确定RF脉冲和倾斜磁场脉冲的施加定时的脉冲序列、指定RF脉冲和倾斜磁场脉冲的施加強度、施加定时等的參数，来决定成为信号处理系统70的CPU71向序列产生器40赋予控制信号的基础的摄像序列。预先设定脉冲序列并保存在存储装置72中。另外，根据由操作者经由输入装置75设定的摄像条件，在信号处理系统70中计算參数。
[0092]本实施方式的信号处理系统70具备:序列作成部，其根据预先保存的脉冲序列、从操作者输入的參数，作成在摄像中使用的摄像序列。通过由CPU71将预先保存在存储装置72或外部存储装置73中的程序加载到存储中并执行，来实现该序列作成部。
[0093]在详细说明本实施方式的序列作成部之前，先说明在本实施方式中使用的上述脉冲序列。在本实施方式中，作为脉冲序列，使用FSE (快速自旋回波)等多回波摄像用的脉冲序列。以下，在本实施方式中，列举使用FSE的情况为例子进行说明。此外，在本实施方式中使用的脉冲序列不限于FSE，只要是在激励脉冲(90°脉冲)后在重复时间TR的时间内施加多个再收敛(重聚焦)RF脉冲的脉冲序列即可。
[0094]在图2中表示FSE脉冲序列200的RF脉冲的施加定时、回波信号的取得定时。如本图所示那样，在施加激励RF脉冲(90°脉冲)201后，在TR期间，施加多个(在此，作为ー个例子示例6个的情況)再收敛RF脉冲202 (202^202^202^202^202^202^,,然后，在施加各再收敛RF脉冲202后，測量回波信号203 (203^203^203^203^203^2036^将得到的回波信号203的信号強度设为SS。此外，在此使用的信号強度SS如上述那样忽略因各种编码产生的强度变化。
[0095]在本实施方式中，使用使该再收敛RF脉冲202的各翻转角(FA =FA1, FA2, FA3、FA4,FA5, FA6)可变的基于VRFA的脉冲序列。
[0096]以下，在本实施方式中，设在ー个激励RF脉冲201后施加N (N是自然数)个再收敛RF脉冲，将各再收敛RF脉冲表示为再收敛RF脉冲202n (n是满足I < n < N的自然数)。按照施加顺序附加下标的n。另外，将第n个施加的再收敛RF脉冲202n的翻转角表示为FAn。另外,将第n个施加的再收敛RF脉冲202n的翻转角FAn称为第n个FA。另外，作为第n个施加的再收敛RF脉冲202n的下ー个测量的回波信号203n，将其回波编号设为n。此外，在不需要特別区别每个回波编号的情况下，称为再收敛RF脉冲202、翻转角FA、回波信号203、信号強度SS。
[0097]进而，将按顺序排列第I个~第n个的各翻转角FA的值后的列(FA值列)称为FA变化形状FAP，将按顺序排列第I个~第n个回波信号203的信号強度SS所得的列称为信号强度变化形状(SSP)。即，FA变化形状FAP由每个再收敛RF脉冲的FA构成，信号强度变化形状SSP由针对每个再收敛RF脉冲得到的每个回波信号(回波编号)的信号強度SS构成。
[0098]在本实施方式中，决定该FA变化形状(各再收敛RF脉冲202的FA值列)，以便即使在进行信号修正的情况下也使所得到的图像的SNR大致为最大。
[0099]为了实现它，本实施方式的序列作成部300如图3 (a)所示具备决定FA变化形状的FA决定部310。FA决定部310进行FA决定处理，即决定再收敛RF脉冲202的FA (FA变化形状)以便在给定的条件下图像的SNR大致为最大。
[0100]FA决定部310例如通过使噪声放大率NA为最小来实现SNR的最大化，噪声放大率NA由构成信号强度变化形状SSP的各回波信号的信号強度SS的倒数的平方的平均的平方根(RMS:均方根)来定义。该噪声放大率NA相当于依照信号強度SS进行信号修正的情况下的忽略了频率特性的噪声放大率。如上述那样，在进行信号修正的情况下，配置在k空间中心以外的回波信号的信号強度也对SNR产生影响。通过使用这样定义的噪声放大率NA，能够还考虑k空间中心以外的信号強度来进行SNR的评价。另外，也可以考虑后处理的滤波器等而考虑进行了其他加权的噪声放大率。
[0101]在计算使噪声放大率`NA为最小的FA变化形状FAP吋，使用普通的最优化法。例如预先设定使用最优化法时的FA变化形状FAP的初始值、即各再收敛RF脉冲的FA值的初始值。也可以从操作者输入初始值作为摄像条件。
[0102]作为所使用的最优化法，例如有最急下降法。例如将全部的FA值的初始值设为0度。在使用最急下降法的情况下，FA决定部310在每次使FA变化形状FAP变化吋，计算噪声放大率NA对FA变化形状FAP的偏微分(在计算中是差商)，从FA变化形状FAP减去将该偏微分乘以预先确定的正的整数所得的值，更新FA变化形状FAP。例如重复预定次数地进行FA变化形状FAP的更新。然后，将最终得到的FA变化形状FAP决定为目的的FA变化形状FAP来输出。此外，也可以对每个再收敛RF脉冲设定与偏微分相乘的正的整数。
[0103]接着，说明本实施方式的信号处理系统70的摄像时的处理流程。图4是本实施方式的摄像时的处理流程。
[0104]如果从操作者输入了摄像条件，则本实施方式的信号处理系统70接受它(步骤SllODo然后，使用所接受的摄像条件，序列作成部300使决定部310决定各再收敛RF脉冲的FA (FA变化形状FAP) (FA决定处理:步骤S1102)，使用所決定的各再收敛RF脉冲的FA (FA变化形状FAP)和摄像条件，作成摄像序列(步骤S1103)。
[0105]信号处理系统70使用所作成的摄像序列，向序列产生器40给出指令，执行摄像(步骤 S1104)。[0106]在此，在图5(a)中表示在步骤SllOl中由操作者设定Tl和T2来作为摄像条件的UI (用户接ロ)画面400。n画面400具备接受摄像条件的输入的摄像条件设定区域410。摄像条件设定区域410具备接受Tl的输入的Tl输入栏401、接受T2的输入的T2输入栏402。该UI画面400被显示在显示装置74上，操作者经由该UI画面400，使用输入装置75设定由在MRI装置100的摄像中所需要的各种參数构成的摄像条件。
[0107]接着，说明步骤S1102的FA决定部310的FA决定处理的流程。
[0108]图6是本实施方式的FA决定处理的处理流程。在此，列举根据最急下降法求出FA变化形状FAP的情况为例子进行说明。另外，根据是否重复进行了 M次(M是自然数)来进行收敛判定。
[0109]FA决定部310向针对重复次数进行计数的计数器i设定I (步骤S1201)。然后，对FA变化形状FAP设定再收敛RF脉冲的FA的初始值(步骤S1202)。例如将初始值全部设为0度。
[0110]接着，FA决定部310计算噪声放大率NA对FA变化形状FAP的偏微分(在计算中是差商)(步骤S1203)。在此，首先使用回波间隔IET (Inter EchoTime)、成为对象的原子核种类的纵向迟豫时间(Tl)和横向迟豫时间(T2)，根据FA变化形状FAP，通过非专利文献2所公开的EPG的方法，计算信号强度变化形状SSP。另外，使用构成信号强度变化形状SSP的各信号強度SS，计算噪声放大率NA，计算其偏微分。此外，在根据FA变化形状FAP计算信号强度变化形状SSP吋，也可以直接使用布洛赫方程式。
[0111]EPG的详细如非专利文献2所公开的那样，但概要如下。在EPG中，分为各散相水平下的纵向磁化成分和横向磁化成分来考虑回波信号203的状态。在此，散相水平表示由于再收敛RF脉冲间的倾斜磁场等，磁化扩散何种程度。如果再收敛RF脉冲是完全的180°脉冲，磁化完全重聚焦，则由于相 同的倾斜磁场，散相水平毎次反转而倒回，但如果不是完全的180°，则出现没有倒回的成分，留有散相水平发展了 I阶段的状态。当施加多个再收敛RF脉冲吋，散相水平积累或抵消，前进到n阶段或倒回到0阶段。
[0112]产生回波信号是在散相水平倒回0阶段吋。因此，通过保存过去的全部散相水平，能够求出散相水平恢复为0阶段的定时。可以只用整数简单地考虑FSE，使得各再收敛RF脉冲之间的倾斜磁场取得平衡，散相水平重叠而成为n阶段。即，根据回波信号203n-l的状态、FAn、Tl、T2、IET决定回波信号203n的状态。能够利用它，逐次求出回波信号203的状态(信号強度)。
[0113]如果得到偏微分，则FA决定部310使用得到的偏微分，更新FA变化形状FAP (步骤S1204)。在此，如上述那样，从更新前的FA变化形状FAP减去将偏微分乘以预定的正的整数所得的值，得到更新后的FA变化形状FAP。
[0114]然后，判别是否重复进行了 M次处理(步骤S1205)，如果重复了，则将该时刻的FA变化形状FAP决定为处理结果(步骤S1208)，结束处理。
[0115]另ー方面，如果在步骤S1205中不满M次，则将i加I (步骤S1206)，返回到步骤S1203。
[0116]此外，为了简单根据重复次数来进行收敛判定，但收敛判定并不限于此。例如，也可以在噪声放大率NA的偏微分足够小的情况下判定为已收敛。在该情况下，预先确定判定所使用的阈值，在偏微分的值成为阈值以下的情况下，判定为已收敛。[0117]另外，在FA决定处理中使用的最优化法并不限于最急下降法。例如也可以整体地使FA变化形状FAP变化，对每个FA变化形状FAP计算噪声放大率NA，将得到计算出的噪声放大率NA中的最小的NA的FA变化形状FAP作为处理結果。
[0118]另外，0度到180度以外的FA具有与从0度到180度相同的效果，因此在最优化法中可以把使FA变化的范围限制为从0度到180度。另外，也可以限制得更窄。如果提高使FA变化的范围的下限值则对运动耐性变强。如果降低使FA变化的范围的上限值，则能够抑制再收敛RF脉冲的功率。这些上限值和下限值既可以预先由MRI装置100保存，也可以由操作者输入。
[0119]图7 (a)和图7 (b)分别表示作为通过上述方法进行探索的结果所得到的使噪声放大率NA成为最小的信号强度变化形状SSP511和FA变化形状FAP512的ー个例子。在此，设为回波链数ETL为80，回波间隔IET为0.0073秒，对象的Tl为I秒，对象的T2为0.1秒。另外，构成信号强度变化形状SSP的各回波信号的信号強度SS如上述那样是将激励瞬间的假想信号强度标准化为I的值。另外，在2个图表中，纵向轴是信号強度和FA (deg),横向轴是回波编号。
[0120]如以上说明的那样，根据本实施方式，导入能够评价考虑到信号修正的SNR的噪声放大率NA，使用最优化法决定使其最小的FA变化形状。因此，根据本实施方式，即使是进行信号修正等，配置在k空间的中心以外的回波信号对SNR产生影响的多回波摄像，也能够决定得到最大的SNR的FA变化形状。
[0121 ] 另外，根据本实施方式，根据由这样决定的FA变化形状所确定的FA施加各再收敛RF脉冲，进行多回波摄像。因此，根据本实施方式，操作者只要输入摄像条件，就能够在多回波摄像中得到考虑到信号修正的SNR最大的图像。
[0122]此外，在本实施方式中，序列作成部300直接使用由FA决定部310決定的FA变化形状FAP来生成摄像序列，但并不限于此。
[0123]也可以构成为与FA的范围的限制等对应地修正由FA决定部310決定的FA变化形状FAP。
[0124]图3 (b)表示该情况下的序列作成部300的功能框图。如本图所示那样，序列作成部300除了 FA决定部310以外，还具备FA修正部311。
[0125]FA修正部311依照预先确定的规则，修正由FA决定部310決定的FA变化形状。以下，使用【专利附图】
【附图说明】具体的修正方法。
[0126]首先，列举限制FA的范围，由FA修正部311进行修正以使根据FA变化形状FAP确定的各再收敛RF脉冲的FA收敛在该范围内的情况为例子进行说明。
[0127]此外，在上述实施方式中，作为变形例子说明应用FA的范围的限制作为FA决定部310中的限制条件的例子，但也可以组合该变形例子和FA变化形状决定后的修正。进而，也可以通过FA决定部310只限制上限值，来修正只得到下限值的FA变化形状FAP等。
[0128]以下，作为ー个例子，说明以下情況，即已经由FA决定部310对上限值进行了限制，由FA修正部311对得到的FA变化形状FAP进行修正使其收敛在下限值FAL。图8 (a)是用于说明该修正方法的ー个例子的说明图。在此，将FA的上限值设为FAHJf FA的下限值设为FAL。
[0129]在该情况下， FA修正部311如本图所示那样，以上限值FAH为中心对修正前的FA变化形状FAP601进行比例缩放，收敛为到下限值FAL为止的之间，成为修正后的FA变化形状FAP'602。在此，以上限值FAL为中心的比例缩放是指将构成修正前的FA变化形状FAP601的各FA值变换为将与上限值的差(FA-FAH)乘以(FAL-FAH) / (Famin-FAH)的结果与上限值FAH相加所得的值。
[0130]序列作成部300使用FA修正部311进行了修正后的FA变化形状FAP'，作成摄像序列。
[0131]此外，使FA收敛在预定的范围内的修正方法也可以是以下这样的修正，即如图8(b)所示，从第一 FA到最小的FA (FAmin:设为第m)为止，通过上述同样的比例缩放使其收敛到下限值FAL内，对于以后的FA，将下限值以下的值的FA全部设为下限值FAL。在此，用图表603表示修正后的FA变化形状FAP'。
[0132]接着，列举FA修正部311修正预定的FA使其成为预定的值的情况为例子进行说明。在该情况下，对其他FA进行修正使其从最小的FA到最大的FA平滑地变化。图9 (a)是用于说明该修正方法的ー个例子的说明图。在此，将得到配置在k空间的中心附近的回波信号的FA设为修正为上述预定的值的FA。将得到配置在k空间的中心附近的回波信号的FA设为第n的FA，将修正后的预定的值设为FAl。其中，修正前的FA变化形状FAP内的最小FA即FAmin为第m，n>m。用图表601表示修正前的FA变化形状FAP，用图表604表示修正后的FA变化形状FAP'。
[0133]在该情况下，例如在第m~第n的FA提前(在时间方向上移动)以使成为FAl的第P个脉冲大致成为第n个后，以第n个FA为中心在FA值方向上进行比例缩放，使得与第m个FA平滑地连接。第m个以后的 FA与上述同样地，以上限值FAH为中心在FA值方向上进行比例缩放，与第n个FA平滑地连接。
[0134]此外，如图9 (b)所示，FA修正部311也可以在FA变化形状FAP的修正中，针对稳定化脉冲部分的FA，排除上限值FAH以下这样的限制。
[0135]此外，稳定化脉冲是为了使回波信号強度稳定而施加的脉冲，由从回波链的开头开始的预定个数的回波信号构成。在该情况下，从最初到FAmin (第m个)的FA以180度为中心通过上述方法在FA值方向上进行比例缩放，对于FAmin以后的FA，将下限以下的FA修正为下限值FAL。用图表605表示修正后的FA变化形状FAP'。
[0136]另外，在上述实施方式中，在毎次设定摄像条件时决定FA，但并不限于此。也可以构成为预先对Tl、T2、ETL、IET这样的在FA决定处理中需要的每个參数计算FA变化形状FAP，将计算结果作为FA变化形状数据库存储在存储装置72或外部存储装置73中。
[0137]在该情况下，序列作成部300也可以不具备FA决定部310。即，序列作成部300当设定了摄像条件吋，从FA变化形状数据库中抽出与之对应的FA变化形状FAP，使用它作成摄像序列。
[0138]在图10中表示该情况下的摄像处理的流程。在此，与图4相同的处理附加相同的符号。如本图所示那样，当从操作者输入了摄像条件时，信号处理系统70接受它(步骤SI 101)。然后，序列作成部300从FA变化形状数据库中抽出与接受的摄像条件对应的FA变化形状FAP (步骤S1112)。然后，序列作成部300使用抽出的FA变化形状FAP、摄像条件，作成摄像序列(步骤S1113)。信号处理系统70使用作成的摄像序列，向序列产生器40给出指令，执行摄像(步骤S1104)。[0139] 此外，也可以构成为在步骤S1112中在FA变化形状数据库中没有与接受的摄像条件对应的FA变化形状FAP的情况下，选择与接受的摄像条件最接近的FA变化形状FAP，或者还可以进行内插。
[0140]此外，在该情况下，也可以构成为在步骤S1112后，修正所抽出的FA变化形状。
[0141]另外，在上述实施方式中，构成为由操作者直接输入Tl、T2作为摄像条件，但操作者指定的信息并不限于此。只要是能够确定Tl、T2的信息即可。例如也可以构成为代替Tl、T2而由操作者设定对象组织。在该情况下，将每个组织的Tl、T2值保存在存储装置72或外部存储装置73中。如果指定了对象组织，则序列作成部300从存储装置72等中抽出与之对应的Tl、T2值，使用它决定FA变化形状FAP，作成摄像序列。
[0142]在该情况下，如图5 (b)所示，UI画面400在摄像条件设定区域410中具备对象组织输入栏403，来取代Tl输入栏401和T2输入栏402。
[0143]此外，在有多个对象组织的情况下，也可以采用摄像对象的全部组织的Tl和T2各自的平均值。另外，也可以使用每个组织的Tl、T2，对每个组织计算上述FA变化形状FAP，对每个组织的FA变化形状FAP进行平均，设为最終的FA变化形状FAP。在计算平均吋，也可以根据摄像对象内部的组织的分布、比例进行加权。
[0144]图11表示使用各个Tl和T2计算每个组织的FA变化形状FAP、然后取得平均时的处理流程。在此，列举组织有K种(I(是自然数)的第一实施方式的情况为例子进行说明。
[0145]如果从操作者输入了摄像条件，则本实施方式的信号处理系统70接受它(步骤S1101)。
[0146]然后，使用所接受的摄像条件，序列作成部300使FA决定部310决定各再收敛RF脉冲的FA。这时，在本变形例子中，FA决定部310对输入的每个Tl和T2进行FA决定处理，计算各自的FA变化形状(步骤S1121、S1102、S1122、S1123)。
[0147]然后，序列作成部300计算每个T1、T2的FA变化形状的平均，将它设为FA变化形状(步骤S1124)。
[0148]然后，序列作成部300使用计算出的各再收敛RF脉冲的FA (FA变化形状)、摄像条件，作成摄像序列(步骤S1103)。信号处理系统70使用作成的摄像序列，向序列产生器40给出指令，执行摄像(步骤S1104)。
[0149]另外，也可以通过预扫描来求出T1、T2。在图12中表示该情况下的摄像处理的流程。本处理基本上与第一实施方式的摄像处理相同，但如果接受了摄像条件，则信号处理系统70进行预扫描(步骤SI 131)，决定Tl、T2 (步骤SI 132)，转移到FA决定处理SI 102。
[0150]此外，在图11和图12所示的变形例子中，还可以从预先保存在存储装置72等中的数据库中，抽出与Tl和T2对应的FA变化形状，来代替FA决定处理。
[0151]〈く第二实施方式〉〉
[0152]接着，说明应用本发明的第二实施方式。在第一实施方式中，使全部再收敛RF脉冲的FA在可能的范围内变化，探索使SNR最大(使噪声放大率NA最小)的FA变化形状。即，探索的參数(探索參数)个数是在一次TR期间施加的全部再收敛RF脉冲数(回波数)N。例如，也可以把与FA之间能够相互变换的信号強度SS作为探索參数，在该情况下，探索參数个数为N。在本实施方式中，使该探索參数个数比N少。
[0153]本实施方式的MRI装置100、序列作成部300的结构基本上与第一实施方式相同。另外，本实施方式的信号处理系统70的摄像处理的流程也与第一实施方式相同。因此，省略对它们的说明。以下，主要着眼于与第一实施方式不同的结构来说明本实施方式。
[0154]在本实施方式中，如上述那样准备回波数N以下的L个(L〈N的自然数)探索參数(总称为CTRL)，使用它决定使SNR最大化的FA变化形状FAP。此外，使用EPG、ProspectiveEPG等，能够与信号强度和FA相互进行变换。因此，探索參数既可以表示FA，也可以表示信号強度。此外，ProspectiveEPG与EPG相反，是根据信号强度变化形状SSP计算FA变化形状FAP的技术。
[0155]在ProspectiveEPG中，通过从信号强度变化形状SSP的开头顺序地求出实现信号強度SS的FA，计算FA变化形状FAP。在该过程中也可能不存在能够实现信号強度SS那样的FA、或不存在于设定的范围内。将这样的状态称为FA已饱和。 [0156]在FA饱和的情况下,通过依照预先确定的规则决定FA,来计算FA变化形状FAP。例如，进行将饱和以后的FA全部设为预先确定的上限值FAH等的处理。
[0157]本实施方式的FA决定部320如图13所示那样，具备决定使SNR最大化的探索參数CTRL的參数决定部321、根据探索參数CTRL计算各FA值(FA变化形状FAP)的FA计算部 322。
[0158]以下，列举具体例子说明本实施方式的FA决定处理。在此，使用隔I个抽取构成FA变化形状FAP的N个各FA值所得的FA值列作为探索參数CTRL。将探索參数CTRL的个数设为L (L=N/2:N是偶数，L= (N+D/2:N是奇数)。
[0159]FA计算部322根据探索參数CTRL计算FA变化形状FAP。例如，在探索參数CTRL是间隔抽取构成FA变化形状FAP的N个FA值而作成的情况下，使用线性内插等方法计算FA变化形状FAP。
[0160]參数决定部321使用FA计算部322将探索參数CTRL变换为FA变化形状FAP，通过与第一实施方式相同的探索方法，决定SNR最大、即噪声放大率NA最小的探索參数CTRL。
[0161]在图14中表示本实施方式的FA决定处理的流程。在本实施方式中，也与第一实施方式同样地，说明使用最急下降法重复M次进行值的更新的情況。
[0162]首先，參数决定部321向针对重复次数进行计数的计数器i设定I (步骤S2101)。然后，对探索參数CTRL设定初始值(步骤S2102)。与第一实施方式同样地，例如将初始值全部设为0度。另外，參数决定部321使用FA计算部322，根据所得到的探索參数CTRL，作成FA变化形状FAP (步骤S2103)。
[0163]接着，參数决定部321计算噪声放大率NA对探索參数CTRL的偏微分(在计算中为差商)(步骤S2104)。在此，与第一实施方式同样地，使用IET、T1、T2，根据在步骤S2103中计算出的FA变化形状FAP，通过EPG等计算信号强度变化形状SSP。然后，使用计算出的信号强度变化形状SSP，计算噪声放大率NA及其偏微分。
[0164]然后，通过与第一实施方式相同的方法更新探索參数CTRL (步骤S2105)。即，从更新前的探索參数CTRL的各值减去将在步骤S2104中得到的偏微分的值乘以预先确定的正的整数所得的值，得到更新后的探索參数CTRL。
[0165]然后，參数决定部321判别是否重复进行了 M次处理(步骤S2106)，如果不满M次，则将计数器i加I (步骤S2107)，返回到步骤S2103。
[0166]另ー方面，在步骤S2106中重复进行了 M次的情况下，參数决定部321将该时刻的探索參数CTRL决定为使SNR最大的探索參数CTRL (步骤S2108)。然后，FA决定部320使FA计算部322根据在步骤S2108中決定的探索參数CTRL计算FA变化形状(步骤S2109)，
结束处理。
[0167]此外，在本实施方式中，探索噪声放大率NA的最小值的方法也不限于此。例如也可以与第一实施方式同样地，不通过最急下降法，而是整体地将根据得到对变化的每个探索參数CTRL计算出的噪声放大率NA中的最小的NA的探索參数CTRL求出的FA变化形状FAP作为处理結果。
[0168]另外，也根据重复次数来进行收敛判定，但并不限于此。例如，也可以与第一实施方式同样地，根据噪声放大率NA的偏微分的值来判定。
[0169]如以上说明的那样，根据本实施方式，与第一实施方式同样地，在进行信号修正等，配置在k空间的中心以外的回波信号对SNR产生影响的情况下，也能够在多回波摄像中决定使SNR最大的FA变化形状，使用该FA变化形状进行摄像，因此能够得到考虑了信号修正的SNR大的图像。
[0170]另外，本实施方式与第一实施方式相比，探索所用的參数少。因此，与第一实施方式相比，能够通过少的计算时间实现最优的SNR。
[0171]此外，降低參数个数的方法并不限于上述方法。例如，作为探索參数，既可以使用表示FA变化形状FAP的曲线的近似式的系数，也可以使用定义表示信号强度变化形状SSP的曲线的公式的系数。如果使用EPG、ProspectiveEPG等，则能够对信号强度和FA相互变换，因此探索參数既可以表示FA，也可以表示信号強度。
[0172]另外，在本实施方式中，也可以与第一实施方式同样地，预先确定在最优化法中使FA变化的范围来进行限制。另外`，也可以具备修正由FA决定部320決定的FA变化形状FAP的FA修正部311。
[0173]另外，还可以将预先作成的FA变化形状FAP保存在存储装置72等中，在摄像中从其中抽出。
[0174]并且，还可以设定摄像对象组织来代替T1、T2。另外，在摄像对象组织有多种的情况下，也可以使用T1、T2值的平均值计算FA变化形状，还可以对每个T1、T2值计算FA变化形状，计算FA变化形状的平均。另外，也可以通过预扫描得到Tl、T2。
[0175]另外，在本实施方式中，在步骤S2103中，根据探索參数CTRL计算FA变化形状FAP、信号强度变化形状SSP，使用信号强度变化形状SSP计算噪声放大率NA的偏微分，但并不限于该方法。例如，在如上述实施方式那样间隔抽取FA变化形状来作成的情况下，也可以不计算没有间隔抽取的FA变化形状而更新探索參数CTRL。即，只根据作为探索參数CTRL被抽出的FA值列计算间隔抽取出的信号强度值列。另外，根据该信号强度值列计算噪声放大率NA及其偏微分。更新探索參数CTRL。
[0176]〈く第三实施方式〉〉
[0177]接着，说明应用本发明的第三实施方式。在上述各实施方式中，作为具体例子，控制直接与FA变化形状有关的參数使噪声放大率最小化。
[0178]在本实施方式中，使信号强度变化形状变化，决定使噪声放大率最小化的信号强度变化形状，将实现该信号强度变化形状的FA变化形状作为在摄像中使用的FA变化形状。即，作为第二实施方式的探索參数，使用确定信号强度变化形状的信息。[0179]本实施方式的MRI装置100、序列作成部300具有与第二实施方式基本相同的结构。另外，本实施方式的信号处理系统70的摄像处理的流程也与第二实施方式相同。因此，在此省略它们的说明。以下，主要着眼于与第二实施方式不同的结构来说明本实施方式。
[0180]在本实施方式中，如上述那样，作为探索參数，使用确定信号强度变化形状SSP的信息。使该探索參数变化，决定使SNR最大化的信号强度变化形状SSP，根据所決定的信号强度变化形状SSP决定FA变化形状FAP。为了实现它，本实施方式的FA决定部330如图15所示，具备决定用于确定使SNR最大化的信号强度变化形状SSP的信息的信号强度变化形状决定部331、根据确定作为探索參数的信号强度变化形状SSP的信息计算FA变化形状FAP的FA计算部332。
[0181]信号强度变化形状决定部331用预定的函数表示信号强度变化形状SSP，决定作为探索參数的函数的參数(系数等)使得噪声放大率NA最小。这时，与第一实施方式、第二实施方式同样地，使用最急下降法等最优化法决定该函数的參数。
[0182]在调查了若干种情况后可知，信号强度变化形状SSP的一部分或全部是渐近地逼近指数函数的形状时，噪声放大率NA大致为最小。在本实施方式中，作为信号強度SS (t)的函数，使用以下的公式(I)表示的二重指数函数。
[0183]
【权利要求】
1.一种磁共振成像装置，具备产生静磁场的静磁场产生系统、向配置在上述静磁场中的被检测体施加倾斜磁场的倾斜磁场产生系统、发送以预定的翻转角激励上述被检测体磁化的高频磁场脉冲的高频磁场发送系统、接收上述被检测体产生的回波信号的信号接收系统、根据上述信号接收系统接收到的回波信号重新构成图像并且依照摄像序列控制上述倾斜磁场产生系统、上述高频磁场发送系统、上述信号接收系统的动作的控制系统， 上述磁共振成像装置的特征在于， 上述控制系统具备: 接受部，其从用户接受决定上述摄像序列的摄像条件； 序列生成部，其依照上述接受部接受的摄像条件和脉冲序列生成上述摄像序列， 上述脉冲序列是在施加I个激励高频磁场脉冲后施加多个再收敛高频磁场脉冲的序列， 上述序列生成部根据上述信号接收系统接收的多个回波信号的信号强度，决定上述各再收敛闻频磁场脉冲的翻转角以使上述图像的SNR最大。
2.根据权利要求1所述的磁共振成像装置，其特征在于， 上述序列生成部使用反映上述多个回波信号的信号修正后的图像的SNR的指标，决定使该图像的SNR最大的上述各再收敛闻频磁场脉冲的翻转角。
3.根据权利要求2所述的磁共振成像装置，其特征在于， 上述序列生成部具备:翻转角决定部，其根据上述摄像条件决定上述各再收敛高频磁场脉冲的翻转角，` 上述翻转角决定部通过依照预先确定的最优化法使确定上述各再收敛高频磁场脉冲的翻转角的信息变化，重复计算上述指标，来决定上述翻转角。
4.根据权利要求3所述的磁共振成像装置，其特征在于， 上述翻转角决定部具备:信号强度计算部，其根据确定上述各再收敛高频磁场脉冲的翻转角的信息，分别计算通过施加与各信息对应的再收敛高频磁场脉冲而得到的回波信号的信号强度， 作为上述指标，使用根据计算出的上述各回波信号的信号强度的倒数的平方的平均的平方根得到的噪声放大率， 决定上述各再收敛高频磁场脉冲的翻转角以使上述噪声放大率最小。
5.根据权利要求3所述的磁共振成像装置，其特征在于， 确定上述各再收敛高频磁场脉冲的翻转角的信息是上述各再收敛高频磁场脉冲的翻转角。
6.根据权利要求3所述的磁共振成像装置，其特征在于， 确定上述各再收敛高频磁场脉冲的翻转角的信息由比上述各再收敛高频磁场脉冲数少的数构成。
7.根据权利要求3所述的磁共振成像装置，其特征在于， 确定上述各再收敛高频磁场脉冲的翻转角的信息是依照预定的规则间隔抽取上述各再收敛高频磁场脉冲的翻转角所得到的中间值列， 上述翻转角决定部具备: 中间值列决定部，其使上述中间值列变化来决定使上述SNR最大的中间值列；翻转角计算部，其对上述决定的中间值进行内插，计算上述各再收敛高频磁场脉冲的翻转角。
8.根据权利要求3所述的磁共振成像装置，其特征在于， 确定上述各再收敛高频磁场脉冲的翻转角的信息是表示通过施加上述各再收敛高频磁场脉冲而得到的回波信号的信号强度的变化形状的函数的参数， 上述翻转角决定部具备:信号强度变化形状决定部，其使上述参数变化来决定使上述SNR最大的上述参数，决定上述信号强度的变化的形状； 翻转角计算部，其根据用上述参数确定的信号强度的变化的形状，计算上述各再收敛高频磁场脉冲的翻转角。
9.根据权利要求8所述的磁共振成像装置，其特征在于， 上述函数是渐近地逼近指数函数的函数， 上述参数是上述指数函数的系数以及时间常数中的至少一方。
10.根据权利要求8所述的磁共振成像装置，其特征在于， 上述函数是二重指数函数， 上述参数是构成上述二重指数函数的2个指数函数各自的系数和时间常数中的一个以上。
11.根据权利要求8所述的磁共振成像装置，其特征在于， 预先确定使上述翻转角变化`的范围， 上述信号强度变化形状决定部使上述参数变化，以使上述翻转角收敛在上述范围内并且不饱和。
12.根据权利要求10所述的磁共振成像装置，其特征在于， 作为上述指标，使用根据计算出的上述各回波信号的信号强度的倒数的平方的平均的平方根得到的噪声放大率， 上述二重指数函数由上述信号强度的变化形状渐近地逼近的第一指数函数、表示渐近的变化的第二指数函数构成， 上述参数是上述第一指数函数的时间常数， 上述信号强度变化形状决定部在固定上述噪声放大率来使上述第一指数函数的时间常数变化时把最后的再收敛高频磁场脉冲的上述翻转角最小的上述第一指数函数的时间常数决定为使上述SNR为最大的上述参数。
13.根据权利要求8所述的磁共振成像装置，其特征在于， 预先确定上述翻转角的范围， 上述翻转角决定部在计算出的翻转角不存在于上述范围内的情况下，依照预先确定的规则决定翻转角。
14.根据权利要求3~13的任意一项所述的磁共振成像装置，其特征在于， 上述序列生成部还具备:翻转角修正部，其依照预先确定的规则修正由上述翻转角决定部决定的翻转角。
15.根据权利要求1或2所述的磁共振成像装置，其特征在于， 还具备:翻转角存储部，其对每个摄像条件存储作为上述各再收敛高频磁场脉冲的翻转角的使上述图像的SNR最大而决定的翻转角，上述序列生成部使用与上述接受部接受的摄像条件对应地存储在上述翻转角存储部中的翻转角来作为上述翻转角。
16.根据权利要求1~13的任意一项所述的磁共振成像装置，其特征在于， 在上述摄像条件中包含用于确定成为摄像对象的核素的纵向迟豫时间和横向迟豫时间的信息。
17.一种翻转角决定方法，其决定在磁共振成像装置中使用的多回波摄像序列的多个再收敛高频磁场脉冲的翻转角， 该翻转角决定方法的特征在于，具有: 使用反映在上述磁共振成像装置中取得的多个回波信号的信号修正后的图像的SNR的指标， 依照预先确定的最优化法，使确定上述各再收敛高频磁场脉冲的翻转角的信息变化，重复计算上述指标，决定上述图像的SNR成为最大的翻转角的翻转角决定步骤。
【文档编号】A61B5/055GK103561645SQ201280024917
【公开日】2014年2月5日 申请日期:2012年5月23日 优先权日:2011年6月9日
【发明者】黑川真次 申请人:株式会社日立医疗器械
【专利摘要】为了在多回波序列中使考虑了信号修正的图像的SNR最大，决定多回波摄像序列的多个再收敛高频磁场脉冲的翻转角，使用反映所取得的多个回波信号的信号修正后的图像的SNR的指标，依照预先确定的最优化法，使确定各再收敛RF脉冲的翻转角的信息变化，重复计算指标，决定图像的SNR最大的翻转角。
【专利说明】磁共振成像装置以及翻转角决定方法
【技术领域】
[0001]本发明涉及ー种测定来自被检测体的氢、磷等的核磁共振(以下称为“NMR”)信号来对核密度分布、迟豫时间分布等进行图像化的核磁共振成像(以下称为“MRI ”)技术，特别涉及ー种使用多回波摄像的VRFA (可变再调焦翻转角)的成像技木。
【背景技术】
[0002]MRI装置是测定构成被检测体特别是人体的组织的原子核自旋所产生的NMR信号(回波信号)，二维或三维地对其头部、腹部、四肢等的形态或功能进行图像化(摄像)的装置。在摄像中，对NMR信号赋予根据倾斜磁场而不同的相位编码，并且进行频率编码，作为时间序列数据来测量。对测量到的回波信号进行二维或三维傅立叶变换，来重新构成图像。
[0003]依照预先确定的摄像序列測量回波信号。该摄像序列有各种各样，但其中，有以下的序列,在施加一次激励高频磁场(RF)脉冲后，一边改变相位编码量一边多次施加再收敛高频磁场(RF)脉冲来收集回波信号，谋求高速化。将这样的摄像称为多回波摄像，将在各再收敛RF脉冲后测量的回波信号群称为回波链。另外，将在施加一次的RF脉冲后得到的回波信号的个数称为回波链数(ETL:Echo Train Length)。
[0004]在多回波摄像中，有时改变再收敛RF脉冲的FA(翻转角:Flip Angle)。将这样的摄像序列称为VRFA多回波序列。在VRFA多回波序列中，为了达到预定的目的而决定各再收敛RF脉冲的FA。
[0005]例如，有决定FA来提高`图像的分辨率的技术(例如參照专利文献I)。在专利文献I所公开的技术中，决定FA以使回波链内的各回波信号的信号強度恒定。这是因为如果各回波信号的信号強度恒定，则能够抑制信号衰减，即使使ETL变长，图像也不会模糊(分辨率提高)。这时的FA典型的是在一旦成为小的值(FAmin)后会向最后的FA平滑地増大。在此，将k空间中心的FA定义为FAcenter,将最后的FA定义为FAmax。
[0006]此外，在本说明书中，假设回波信号的信号强度忽略由于相位编码、频率编码造成的强度变化。另外，对于信号強度的值，将激励的瞬间的假想信号强度标准化为I。
[0007]但是，画质不只由分辨率决定，还由对运动的耐性、SNR (信噪比)等各种要素決定。特别是SNR，在决定画质的好坏上是主要的。在专利文献I所公开的技术中，以使信号強度恒定为目的，因此所得到的图像的SNR并不一定最优。另外，恒定的值的决定方法是任意的，并不知道在怎样的值时成为怎样的画质。
[0008]作为提高该SNR的技术，有决定FA以便提高配置在k空间中心的回波信号的信号強度的技术(例如參照非专利文献I)。这是因为一般配置在k空间的中心的回波信号决定了 SNR。在非专利文献I中，决定FA以使在指定了 FAmin、FAcenter> FAmax之后信号平滑地变化。这是因为通过増大FAmin能够提高图像对运动的耐性，通过提高FAcenter能够提高不考虑信号修正时的SNR。
[0009]另外，有以下的技木，即在将2种数据相加(合成)计算2种对比度的图像时，决定FA使得各数据的信号强度统一(例如參照非专利文献2)。在非专利文献2中，作为2种数据，取得CSF (脑脊液)和脑实质的信号強度的平衡不同的数据。将这些数据相加，得到T2W和FLAIR的图像。这时，决定FA使得2种数据的CSF的信号強度在k空间中心的区域中统一。通过在取得2种数据的时刻使信号强度统一，在相加时不需要进行用于统一信号強度的加权，提闻了 SNR。
[0010]此外，作为根据所決定的FA计算信号強度的变化的方法，已知EPG (扩展相位图)(例如參照非专利文献3)。另外，相反，作为根据信号強度的变化计算FA的方法，已知ProspectiveEPG (例如參照非专利文献4)。
[0011]现有技术文献
[0012]非专利文献
[0013]非专利文献1:Reed F.Busse, Anja C.S.Brau, Anthony Vu, CharlesR.Michelich, Ersin Bayram, Richard Kijowski, Scott B.Reeder, and HowardA.Rowley."Effects of Refocusing Flip Angle Modulation and View Ordering in3DFast Spin Ech0."Magnetic Resonance in Medicine60:640649(2008)
[0014]非专利文献2:H.Lee, J-S.Seo, and J.Park.^SNR-Optimized AcceleratedPhase-Sensitive Dual-Acquisition Single_Slab3D Turbo Spin Echo Imaging〃Proc.1ntl.Soc.Mag.Reson.Med.19:386 (2011)
[0015]非专利文献 3:Juergen Hennig."Multieco Imaging Sequences with LowRefocusing Flip Angles^Journal of Magnetic Resonance78:397407(1988)
[0016]非专利文献4: J uergen Hennig, Matthias Weigel, and KlausSheffler."Calculation of Flip Anbgles for Echo Trains With PredefinedAmplitudes With the Extended Phase Graph (EPG)-Algorithm:Principles andApplications to Hyperecho and TRAPS Sequences."Magnetic Resonance inMedicine51:6880(2004)
【发明内容】
[0017]发明要解决的问题
[0018]一般在MRI中，以提高画质等为目的，进行各种信号修正。例如在多回波摄像中，有时修正信号的衰减(T2修正)来抑制图像的模糊。如果进行了这样的信号修正，则与回波信号強度对应地噪声提高，因此全部的回波信号的信号强度对SNR有影响。因此，在进行信号修正的情况下，在非专利文献I所公开的方法中无法进行充分的SNR的控制。
[0019]另外，在非专利文献2所公开的技术中，也只是在k空间中心的区域中统ー 2种图像的CSF的信号強度，一般没有考虑到回波链整体的信号強度来使SNR最大。因此，在进行信号修正的情况下，在非专利文献2所公开的方法中也无法进行SNR的控制。
[0020]本发明就是鉴于上述情况而提出的，其目的在于，提供ー种在VRFA多回波序列中使考虑到信号修正的SNR最大的技术。
[0021]解决问题的方案
[0022]本发明用于决定磁共振成像装置中的多回波摄像序列的、多个再收敛高频磁场脉冲的翻转角，使用反映所取得的多个回波信号的信号修正后的图像的SNR的指标，依照预先确定的最优化法，改变用于确定各再收敛RF脉冲的翻转角的信息，重复计算指标，决定图像的SNR最大的翻转角。
[0023]例如，使用还反映配置在k空间的中心以外的回波信号的信号強度的指标，决定VRFA多回波序列的再收敛RF脉冲的FA，使得SNR大致成为最大。根据不只是反映配置在k空间的中心的回波信号的信号強度，还反映配置在k空间的中心以外的回波信号的信号强度的指标，判别改变再收敛RF脉冲的FA时的SNR的大小。
[0024]发明效果
[0025]根据本发明，在VRFA多回波序列中，能够使考虑到信号修正的SNR最大。 【专利附图】
【附图说明】
[0026]图1是第一实施方式的MRI装置的框图。
[0027]图2是用于说明FSE序列的脉冲序列的说明图。
[0028]图3 Ca)是第一实施方式的序列作成部的功能框图，(b)是第一实施方式的变形例的序列作成部的功能框图。
[0029]图4是第一实施方式的摄像处理的流程图。
[0030]图5 (a)是用于说明第一实施方式的UI画面的说明图，(b)是用于说明第一实施方式的n画面的变形例子的说明图。
[0031]图6是第一实施方式的FA决定处理的流程图。
[0032]图7 (a)、(b)分别是通过第一实施方式的FA决定处理得到的信号强度变化形状和FA变化形状的图表。
[0033]图8 (a)、(b)是用于说明第一实施方式的FA修正部的修正处理的说明图。
[0034]图9 (a)、(b)是用于说明第一实施方式的FA修正部的修正处理的说明图。
[0035]图10是第一实施方式的摄像处理的变形例的流程图。
[0036]图11是第一实施方式的摄像处理的变形例的流程图。
[0037]图12是第一实施方式的摄像处理的变形例的流程图。
[0038]图13是第二实施方式的序列作成部的功能框图。
[0039]图14是第二实施方式的FA决定处理的流程图。
[0040]图15是第三实施方式的序列作成部的功能框图。
[0041]图16是第三实施方式的FA决定处理的流程图。
[0042]图17 (a)、(b)分别是通过第三实施方式的FA决定处理得到的信号强度变化形状和FA变化形状的图表。
[0043]图18是第四实施方式的序列作成部的功能框图。
[0044]图19是第四实施方式的FA决定处理的流程图。
[0045]图20是第四实施方式的信号强度变化形状决定处理的流程图。
[0046]图21是第四实施方式的最小噪声放大率计算处理的流程图。
[0047]图22是第四实施方式的最终FA最小化处理的流程图。
[0048]图23是用于说明第四实施方式的UI画面的ー个例子的说明图。
[0049]图24 (a)、(b)分别是通过第四实施方式的FA决定处理得到的信号强度变化形状和FA变化形状的图表。
[0050]图25是第四实施方式的修正后的FA变化形状的图表。【具体实施方式】
[0051]本发明的用于决定多回波摄像序列的多个再收敛高频磁场脉冲的翻转角的结构的概要如下所述。
[0052]ー种磁共振成像装置，具备产生静磁场的静磁场产生系统、向配置在上述静磁场中的被检测体施加倾斜磁场的倾斜磁场产生系统、发送以预定的翻转角激励上述被检测体磁化的高频磁场脉冲的高频磁场发送系统、接收上述被检测体产生的回波信号的信号接收系统、根据上述信号接收系统接收到的回波信号重新构成图像并且依照摄像序列控制上述倾斜磁场产生系统、上述高频磁场发送系统、上述信号接收系统的动作的控制系统，
[0053]上述控制系统具备:
[0054]接受部，其从用户接受决定上述摄像序列的摄像条件；序列生成部，其依照上述接受部接受的摄像条件和脉冲序列生成上述摄像序列，上述脉冲序列是在施加I个激励高频磁场脉冲后施加多个再收敛高频磁场脉冲的序列，
[0055]上述序列生成部根据上述信号接收系统接收的多个回波信号的信号強度，决定上述各再收敛闻频磁场脉冲的翻转角以使上述图像的SNR最大。
[0056]在上述磁共振成像装置中，上述序列生成部可以使用反映上述多个回波信号的信号修正后的图像的SNR的指标，决定使该图像的SNR最大的上述各再收敛闻频磁场脉冲的翻转角。
[0057]在上述磁共振成像装置中，例如上述序列生成部具备:翻转角决定部，其根据上述摄像条件决定上述各再收敛高频磁场脉冲的翻转角，上述翻转角决定部可以通过依照预先确定的最优化法使确定上述各再收敛高频磁场脉冲的翻转角的信息变化，重复计算上述指标，来决定上述翻转角。
[0058]在上述磁共振成像装置中，例如上述翻转角决定部具备:信号強度计算部，其根据确定上述各再收敛高频磁场脉冲的翻转角的信息，分别计算通过施加与各信息对应的再收敛高频磁场脉冲而得到的回波信号的信号強度，作为上述指标，可以使用根据计算出的上述各回波信号的信号強度的倒数的平方的平均的平方根得到的噪声放大率。可以决定上述各再收敛高频磁场脉冲的翻转角以使其噪声放大率最小。
[0059]在上述磁共振成像装置中，例如确定上述各再收敛高频磁场脉冲的翻转角的信息是再收敛高频磁场脉冲的翻转角。
[0060]在上述磁共振成像装置中，例如确定上述各再收敛高频磁场脉冲的翻转角的信息由比再收敛高频磁场脉冲少的数构成。
[0061]在上述磁共振成像装置中，例如确定上述各再收敛高频磁场脉冲的翻转角的信息可以是依照预定的规则间隔抽取各再收敛高频磁场脉冲的翻转角所得到的中间值列，上述翻转角决定部具备:中间值列决定部，其使上述中间值列变化来决定使上述SNR最大的中间值列；翻转角计算部，其对上述决定的中间值进行内插，计算上述各再收敛高频磁场脉冲的翻转角。
[0062]在上述磁共振成 像装置中，例如确定上述各再收敛高频磁场脉冲的翻转角的信息可以是表示通过施加各再收敛高频磁场脉冲而得到的回波信号的信号強度的变化形状的函数的參数，上述翻转角决定部具备:信号强度变化形状决定部，其使上述參数变化来决定使上述SNR最大的上述參数，决定上述信号強度的变化的形状；翻转角计算部，其根据用上述參数确定的信号強度的变化的形状，计算上述各再收敛高频磁场脉冲的翻转角。
[0063]在上述磁共振成像装置中，例如上述函数可以是渐近地逼近指数函数的函数，上述參数可以是上述指数函数的系数以及时间常数中的至少一方。
[0064]在上述磁共振成像装置中，例如上述函数可以是二重指数函数，上述參数是构成上述二重指数函数的2个指数函数各自的系数和时间常数中的ー个以上。
[0065]在上述磁共振成像装置中，上述二重指数函数也可以是从通过上述激励高频磁场脉冲激励的瞬间的虚拟的信号強度衰减的函数。
[0066]在上述磁共振成像装置中，例如可以预先确定使上述翻转角变化的范围，上述信号强度变化形状决定部使上述參数变化，以使上述翻转角收敛在上述范围内并且不饱和。
[0067]在上述磁共振成像装置中，例如作为上述指标，可以使用根据计算出的上述各回波信号的信号強度的倒数的平方的平均的平方根得到的噪声放大率，上述二重指数函数由上述信号強度的变化形状渐近地逼近的第一指数函数、表示渐近的变化的第二指数函数构成，上述參数是上述第一指数函数的时间常数，上述信号强度变化形状决定部在固定上述噪声放大率来使上述第一指数函数的时间常数变化时把最后的再收敛高频磁场脉冲的上述翻转角最小的上述第一指数函数的时间常数决定为使上述SNR为最大的上述參数。
[0068]在上述磁共振成 像装置中，例如可以预先确定上述翻转角的范围，上述翻转角决定部在计算出的翻转角不存在于上述范围内的情况下，依照预先确定的规则决定翻转角。
[0069]在上述磁共振成像装置中，例如上述最优化法也可以是缩小解的范围的探索法。
[0070]在上述磁共振成像装置中，例如上述序列生成部还可以具备:翻转角修正部，其依照预先确定的规则修正由上述翻转角决定部決定的翻转角。
[0071]在上述磁共振成像装置中，例如上述翻转角修正部也可以进行修正使得由上述翻转角决定部決定的翻转角的一部分或全部收敛在预先确定的可设定范围内。
[0072]在上述磁共振成像装置中，例如上述翻转角修正部也可以在将预先确定的上述再收敛高频磁场脉冲的翻转角修正为预先确定的值后，对该翻转角以外的最小的翻转角以后的翻转角进行修正使其平滑地变化。
[0073]在上述磁共振成像装置中，例如修正为上述预先确定的值的翻转角可以是得到配置在k空间的中心的回波信号的翻转角。
[0074]在上述磁共振成像装置中，例如还可以具备:翻转角存储部，其对每个摄像条件存储作为上述再收敛高频磁场脉冲的翻转角的使上述图像的SNR最大而决定的翻转角，上述序列生成部使用与上述接受部接受的摄像条件对应地存储在上述翻转角存储部中的翻转角来作为上述翻转角。
[0075]在上述磁共振成像装置中，例如上述序列生成部也可以在上述翻转角存储部中没有存储与上述摄像条件对应的翻转角的情况下，使用与最接近的摄像条件对应地存储的翻转角作为上述翻转角。
[0076]在上述磁共振成像装置中，例如上述序列生成部也可以在上述翻转角存储部中没有存储与上述摄像条件对应的翻转角的情况下，使用对与最接近的摄像条件对应地存储的翻转角进行内插而计算出的翻转角作为上述翻转角。
[0077]在上述磁共振成像装置中，例如在上述摄像条件中包含用于确定成为摄像对象的核素的纵向迟豫时间和横向迟豫时间的信息。
[0078]另外，一种翻转角决定方法，其决定在磁共振成像装置中使用的多回波摄像序列的多个再收敛高频磁场脉冲的翻转角，该翻转角决定方法具有:使用反映在上述磁共振成像装置中取得的多个回波信号的信号修正后的图像的SNR的指标，依照预先确定的最优化法，使确定上述各再收敛高频磁场脉冲的翻转角的信息变化，重复计算上述指标，决定上述图像的SNR成为最大的翻转角的翻转角决定步骤。
[0079]以下，详细说明各实施方式。
[0080]〈く第一实施方式》
[0081]说明应用本发明的第一实施方式。以下，在用于说明本发明的实施方式的全部图中，具有相同功能的部分赋予相同符号，省略其重复的说明。
[0082]最初，根据图1说明本实施方式的MRI装置100的概要。图1是表示本实施方式的MRI装置100的整体结构的框图。本实施方式的MRI装置100利用NMR现象得到被检测体的断层图像，具备静磁场产生系统20、倾斜磁场产生系统30、序列产生器40、高频磁场发送系统(发送系统)50、信号接收系统(接收系统)60、控制系统(序列产生器40和信号处理系统70)。
[0083]静磁场产生系统20如果是垂直磁场方式，则在被检测体10的周围的空间中在与其体轴垂直的方向上产生均匀的静磁场，如果是水平磁场方式，则在体轴方向上产生均匀的静磁场，在被检测体10的周围配置永磁体方式、常导方式或超导方式的静磁场产生源(例如磁铁)。
[0084]倾斜磁场产生系统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信号(回波信号)编码各个方向的位置信息。
[0085]序列产生器40控制倾斜磁场产生系统30、发送系统50、接收系统60，使得依照预定的摄像序列重复施加高频磁场脉冲(以下称为“RF”脉冲)和倾斜磁场脉沖。序列产生器40依照来自后述的信号处理系统70具备的CPU71的控制信号而动作，向它们发送收集被检测体10的断层图像的数据所需要的各种指令。
[0086]发送系统50为了使构成被检测体10的生物体组织的原子的原子核自旋产生核磁共振，向被检测体10照射RF脉冲，具备高频振动器(合成器)52、调制器53、高频放大器54、发送侧的高频线圈(发送线圈)51。在基于来自序列产生器40的指令的定时下通过调制器53对从合成器52输出的高频脉冲进行振幅调制，通过高频放大器54对该振幅调制后的高频脉冲进行放大，提供给接近被检测体10配置的发送线圈51，由此向被检测体10照射RF脉冲。
[0087]接收系统60检测通过构成被检测体10的生物体组织的原子核自旋的核磁共振释放的回波信号(NMR信号)，具备接收侧的高频线圈(接收线圈)61、信号放大器62、正交相位检波器63、A/D变换器64。通过从发送线圈51照射的电磁波而感应出的被检测体10的响应的回波信号被接近被检测体10配置的接收线圈61检测，在通过信号放大器62放大后，在基于来自序列产生器40的指令的定时下通过正交相位检波器63分割为正交的两个系统的信号，分别通过A/D变换器64变换为数字量，发送到信号处理系统70。
[0088]信号处理系统70进行各种数据处理、处理结果的显示和保存等，具备CPU71、存储装置72、外部存储装置73、显示装置74、输入装置75。例如使用来自接收系统60的数据，重新构成被检测体10的断层图像。另外，依照序列，向序列产生器40发送控制信号。重新构成出的断层图像被显示在显示装置74上，并且记录到存储装置72或外部存储装置73中。输入装置75用于由操作者输入MRI装置100的各种控制信息、由信号处理系统70进行的处理的控制信息，具备跟踪球或鼠标、以及键盘。该输入装置75接近显示装置74地配置，操作者ー边看着显示装置74，ー边通过输入装置75交互地控制MRI装置100的各种处理。
[0089]此外，在图1中，发送线圈51和倾斜磁场线圈31被设置得在插入被检测体10的静磁场产生系统20的静磁场空间内，如果是垂直磁场方式则与被检测体10相対，如果是水平磁场方式则围住被检测体10。另外，接收线圈61被设置得与被检测体10相对或围住被检测体10。
[0090]现在，在临床普及的MRI装置的摄像对象原子核种类是作为被检测体的主要构成物质的氢原子核(质子)。通过对与质子密度的空间分布、激励状态的迟豫时间的空间分布相关的信息进行图像化，二维或三维地拍摄人体头部、腹部、四肢等的形态或功能。
[0091]根据确定RF脉冲和倾斜磁场脉冲的施加定时的脉冲序列、指定RF脉冲和倾斜磁场脉冲的施加強度、施加定时等的參数，来决定成为信号处理系统70的CPU71向序列产生器40赋予控制信号的基础的摄像序列。预先设定脉冲序列并保存在存储装置72中。另外，根据由操作者经由输入装置75设定的摄像条件，在信号处理系统70中计算參数。
[0092]本实施方式的信号处理系统70具备:序列作成部，其根据预先保存的脉冲序列、从操作者输入的參数，作成在摄像中使用的摄像序列。通过由CPU71将预先保存在存储装置72或外部存储装置73中的程序加载到存储中并执行，来实现该序列作成部。
[0093]在详细说明本实施方式的序列作成部之前，先说明在本实施方式中使用的上述脉冲序列。在本实施方式中，作为脉冲序列，使用FSE (快速自旋回波)等多回波摄像用的脉冲序列。以下，在本实施方式中，列举使用FSE的情况为例子进行说明。此外，在本实施方式中使用的脉冲序列不限于FSE，只要是在激励脉冲(90°脉冲)后在重复时间TR的时间内施加多个再收敛(重聚焦)RF脉冲的脉冲序列即可。
[0094]在图2中表示FSE脉冲序列200的RF脉冲的施加定时、回波信号的取得定时。如本图所示那样，在施加激励RF脉冲(90°脉冲)201后，在TR期间，施加多个(在此，作为ー个例子示例6个的情況)再收敛RF脉冲202 (202^202^202^202^202^202^,,然后，在施加各再收敛RF脉冲202后，測量回波信号203 (203^203^203^203^203^2036^将得到的回波信号203的信号強度设为SS。此外，在此使用的信号強度SS如上述那样忽略因各种编码产生的强度变化。
[0095]在本实施方式中，使用使该再收敛RF脉冲202的各翻转角(FA =FA1, FA2, FA3、FA4,FA5, FA6)可变的基于VRFA的脉冲序列。
[0096]以下，在本实施方式中，设在ー个激励RF脉冲201后施加N (N是自然数)个再收敛RF脉冲，将各再收敛RF脉冲表示为再收敛RF脉冲202n (n是满足I < n < N的自然数)。按照施加顺序附加下标的n。另外，将第n个施加的再收敛RF脉冲202n的翻转角表示为FAn。另外,将第n个施加的再收敛RF脉冲202n的翻转角FAn称为第n个FA。另外，作为第n个施加的再收敛RF脉冲202n的下ー个测量的回波信号203n，将其回波编号设为n。此外，在不需要特別区别每个回波编号的情况下，称为再收敛RF脉冲202、翻转角FA、回波信号203、信号強度SS。
[0097]进而，将按顺序排列第I个~第n个的各翻转角FA的值后的列(FA值列)称为FA变化形状FAP，将按顺序排列第I个~第n个回波信号203的信号強度SS所得的列称为信号强度变化形状(SSP)。即，FA变化形状FAP由每个再收敛RF脉冲的FA构成，信号强度变化形状SSP由针对每个再收敛RF脉冲得到的每个回波信号(回波编号)的信号強度SS构成。
[0098]在本实施方式中，决定该FA变化形状(各再收敛RF脉冲202的FA值列)，以便即使在进行信号修正的情况下也使所得到的图像的SNR大致为最大。
[0099]为了实现它，本实施方式的序列作成部300如图3 (a)所示具备决定FA变化形状的FA决定部310。FA决定部310进行FA决定处理，即决定再收敛RF脉冲202的FA (FA变化形状)以便在给定的条件下图像的SNR大致为最大。
[0100]FA决定部310例如通过使噪声放大率NA为最小来实现SNR的最大化，噪声放大率NA由构成信号强度变化形状SSP的各回波信号的信号強度SS的倒数的平方的平均的平方根(RMS:均方根)来定义。该噪声放大率NA相当于依照信号強度SS进行信号修正的情况下的忽略了频率特性的噪声放大率。如上述那样，在进行信号修正的情况下，配置在k空间中心以外的回波信号的信号強度也对SNR产生影响。通过使用这样定义的噪声放大率NA，能够还考虑k空间中心以外的信号強度来进行SNR的评价。另外，也可以考虑后处理的滤波器等而考虑进行了其他加权的噪声放大率。
[0101]在计算使噪声放大率`NA为最小的FA变化形状FAP吋，使用普通的最优化法。例如预先设定使用最优化法时的FA变化形状FAP的初始值、即各再收敛RF脉冲的FA值的初始值。也可以从操作者输入初始值作为摄像条件。
[0102]作为所使用的最优化法，例如有最急下降法。例如将全部的FA值的初始值设为0度。在使用最急下降法的情况下，FA决定部310在每次使FA变化形状FAP变化吋，计算噪声放大率NA对FA变化形状FAP的偏微分(在计算中是差商)，从FA变化形状FAP减去将该偏微分乘以预先确定的正的整数所得的值，更新FA变化形状FAP。例如重复预定次数地进行FA变化形状FAP的更新。然后，将最终得到的FA变化形状FAP决定为目的的FA变化形状FAP来输出。此外，也可以对每个再收敛RF脉冲设定与偏微分相乘的正的整数。
[0103]接着，说明本实施方式的信号处理系统70的摄像时的处理流程。图4是本实施方式的摄像时的处理流程。
[0104]如果从操作者输入了摄像条件，则本实施方式的信号处理系统70接受它(步骤SllODo然后，使用所接受的摄像条件，序列作成部300使决定部310决定各再收敛RF脉冲的FA (FA变化形状FAP) (FA决定处理:步骤S1102)，使用所決定的各再收敛RF脉冲的FA (FA变化形状FAP)和摄像条件，作成摄像序列(步骤S1103)。
[0105]信号处理系统70使用所作成的摄像序列，向序列产生器40给出指令，执行摄像(步骤 S1104)。[0106]在此，在图5(a)中表示在步骤SllOl中由操作者设定Tl和T2来作为摄像条件的UI (用户接ロ)画面400。n画面400具备接受摄像条件的输入的摄像条件设定区域410。摄像条件设定区域410具备接受Tl的输入的Tl输入栏401、接受T2的输入的T2输入栏402。该UI画面400被显示在显示装置74上，操作者经由该UI画面400，使用输入装置75设定由在MRI装置100的摄像中所需要的各种參数构成的摄像条件。
[0107]接着，说明步骤S1102的FA决定部310的FA决定处理的流程。
[0108]图6是本实施方式的FA决定处理的处理流程。在此，列举根据最急下降法求出FA变化形状FAP的情况为例子进行说明。另外，根据是否重复进行了 M次(M是自然数)来进行收敛判定。
[0109]FA决定部310向针对重复次数进行计数的计数器i设定I (步骤S1201)。然后，对FA变化形状FAP设定再收敛RF脉冲的FA的初始值(步骤S1202)。例如将初始值全部设为0度。
[0110]接着，FA决定部310计算噪声放大率NA对FA变化形状FAP的偏微分(在计算中是差商)(步骤S1203)。在此，首先使用回波间隔IET (Inter EchoTime)、成为对象的原子核种类的纵向迟豫时间(Tl)和横向迟豫时间(T2)，根据FA变化形状FAP，通过非专利文献2所公开的EPG的方法，计算信号强度变化形状SSP。另外，使用构成信号强度变化形状SSP的各信号強度SS，计算噪声放大率NA，计算其偏微分。此外，在根据FA变化形状FAP计算信号强度变化形状SSP吋，也可以直接使用布洛赫方程式。
[0111]EPG的详细如非专利文献2所公开的那样，但概要如下。在EPG中，分为各散相水平下的纵向磁化成分和横向磁化成分来考虑回波信号203的状态。在此，散相水平表示由于再收敛RF脉冲间的倾斜磁场等，磁化扩散何种程度。如果再收敛RF脉冲是完全的180°脉冲，磁化完全重聚焦，则由于相 同的倾斜磁场，散相水平毎次反转而倒回，但如果不是完全的180°，则出现没有倒回的成分，留有散相水平发展了 I阶段的状态。当施加多个再收敛RF脉冲吋，散相水平积累或抵消，前进到n阶段或倒回到0阶段。
[0112]产生回波信号是在散相水平倒回0阶段吋。因此，通过保存过去的全部散相水平，能够求出散相水平恢复为0阶段的定时。可以只用整数简单地考虑FSE，使得各再收敛RF脉冲之间的倾斜磁场取得平衡，散相水平重叠而成为n阶段。即，根据回波信号203n-l的状态、FAn、Tl、T2、IET决定回波信号203n的状态。能够利用它，逐次求出回波信号203的状态(信号強度)。
[0113]如果得到偏微分，则FA决定部310使用得到的偏微分，更新FA变化形状FAP (步骤S1204)。在此，如上述那样，从更新前的FA变化形状FAP减去将偏微分乘以预定的正的整数所得的值，得到更新后的FA变化形状FAP。
[0114]然后，判别是否重复进行了 M次处理(步骤S1205)，如果重复了，则将该时刻的FA变化形状FAP决定为处理结果(步骤S1208)，结束处理。
[0115]另ー方面，如果在步骤S1205中不满M次，则将i加I (步骤S1206)，返回到步骤S1203。
[0116]此外，为了简单根据重复次数来进行收敛判定，但收敛判定并不限于此。例如，也可以在噪声放大率NA的偏微分足够小的情况下判定为已收敛。在该情况下，预先确定判定所使用的阈值，在偏微分的值成为阈值以下的情况下，判定为已收敛。[0117]另外，在FA决定处理中使用的最优化法并不限于最急下降法。例如也可以整体地使FA变化形状FAP变化，对每个FA变化形状FAP计算噪声放大率NA，将得到计算出的噪声放大率NA中的最小的NA的FA变化形状FAP作为处理結果。
[0118]另外，0度到180度以外的FA具有与从0度到180度相同的效果，因此在最优化法中可以把使FA变化的范围限制为从0度到180度。另外，也可以限制得更窄。如果提高使FA变化的范围的下限值则对运动耐性变强。如果降低使FA变化的范围的上限值，则能够抑制再收敛RF脉冲的功率。这些上限值和下限值既可以预先由MRI装置100保存，也可以由操作者输入。
[0119]图7 (a)和图7 (b)分别表示作为通过上述方法进行探索的结果所得到的使噪声放大率NA成为最小的信号强度变化形状SSP511和FA变化形状FAP512的ー个例子。在此，设为回波链数ETL为80，回波间隔IET为0.0073秒，对象的Tl为I秒，对象的T2为0.1秒。另外，构成信号强度变化形状SSP的各回波信号的信号強度SS如上述那样是将激励瞬间的假想信号强度标准化为I的值。另外，在2个图表中，纵向轴是信号強度和FA (deg),横向轴是回波编号。
[0120]如以上说明的那样，根据本实施方式，导入能够评价考虑到信号修正的SNR的噪声放大率NA，使用最优化法决定使其最小的FA变化形状。因此，根据本实施方式，即使是进行信号修正等，配置在k空间的中心以外的回波信号对SNR产生影响的多回波摄像，也能够决定得到最大的SNR的FA变化形状。
[0121 ] 另外，根据本实施方式，根据由这样决定的FA变化形状所确定的FA施加各再收敛RF脉冲，进行多回波摄像。因此，根据本实施方式，操作者只要输入摄像条件，就能够在多回波摄像中得到考虑到信号修正的SNR最大的图像。
[0122]此外，在本实施方式中，序列作成部300直接使用由FA决定部310決定的FA变化形状FAP来生成摄像序列，但并不限于此。
[0123]也可以构成为与FA的范围的限制等对应地修正由FA决定部310決定的FA变化形状FAP。
[0124]图3 (b)表示该情况下的序列作成部300的功能框图。如本图所示那样，序列作成部300除了 FA决定部310以外，还具备FA修正部311。
[0125]FA修正部311依照预先确定的规则，修正由FA决定部310決定的FA变化形状。以下，使用【专利附图】
【附图说明】具体的修正方法。
[0126]首先，列举限制FA的范围，由FA修正部311进行修正以使根据FA变化形状FAP确定的各再收敛RF脉冲的FA收敛在该范围内的情况为例子进行说明。
[0127]此外，在上述实施方式中，作为变形例子说明应用FA的范围的限制作为FA决定部310中的限制条件的例子，但也可以组合该变形例子和FA变化形状决定后的修正。进而，也可以通过FA决定部310只限制上限值，来修正只得到下限值的FA变化形状FAP等。
[0128]以下，作为ー个例子，说明以下情況，即已经由FA决定部310对上限值进行了限制，由FA修正部311对得到的FA变化形状FAP进行修正使其收敛在下限值FAL。图8 (a)是用于说明该修正方法的ー个例子的说明图。在此，将FA的上限值设为FAHJf FA的下限值设为FAL。
[0129]在该情况下， FA修正部311如本图所示那样，以上限值FAH为中心对修正前的FA变化形状FAP601进行比例缩放，收敛为到下限值FAL为止的之间，成为修正后的FA变化形状FAP'602。在此，以上限值FAL为中心的比例缩放是指将构成修正前的FA变化形状FAP601的各FA值变换为将与上限值的差(FA-FAH)乘以(FAL-FAH) / (Famin-FAH)的结果与上限值FAH相加所得的值。
[0130]序列作成部300使用FA修正部311进行了修正后的FA变化形状FAP'，作成摄像序列。
[0131]此外，使FA收敛在预定的范围内的修正方法也可以是以下这样的修正，即如图8(b)所示，从第一 FA到最小的FA (FAmin:设为第m)为止，通过上述同样的比例缩放使其收敛到下限值FAL内，对于以后的FA，将下限值以下的值的FA全部设为下限值FAL。在此，用图表603表示修正后的FA变化形状FAP'。
[0132]接着，列举FA修正部311修正预定的FA使其成为预定的值的情况为例子进行说明。在该情况下，对其他FA进行修正使其从最小的FA到最大的FA平滑地变化。图9 (a)是用于说明该修正方法的ー个例子的说明图。在此，将得到配置在k空间的中心附近的回波信号的FA设为修正为上述预定的值的FA。将得到配置在k空间的中心附近的回波信号的FA设为第n的FA，将修正后的预定的值设为FAl。其中，修正前的FA变化形状FAP内的最小FA即FAmin为第m，n>m。用图表601表示修正前的FA变化形状FAP，用图表604表示修正后的FA变化形状FAP'。
[0133]在该情况下，例如在第m~第n的FA提前(在时间方向上移动)以使成为FAl的第P个脉冲大致成为第n个后，以第n个FA为中心在FA值方向上进行比例缩放，使得与第m个FA平滑地连接。第m个以后的 FA与上述同样地，以上限值FAH为中心在FA值方向上进行比例缩放，与第n个FA平滑地连接。
[0134]此外，如图9 (b)所示，FA修正部311也可以在FA变化形状FAP的修正中，针对稳定化脉冲部分的FA，排除上限值FAH以下这样的限制。
[0135]此外，稳定化脉冲是为了使回波信号強度稳定而施加的脉冲，由从回波链的开头开始的预定个数的回波信号构成。在该情况下，从最初到FAmin (第m个)的FA以180度为中心通过上述方法在FA值方向上进行比例缩放，对于FAmin以后的FA，将下限以下的FA修正为下限值FAL。用图表605表示修正后的FA变化形状FAP'。
[0136]另外，在上述实施方式中，在毎次设定摄像条件时决定FA，但并不限于此。也可以构成为预先对Tl、T2、ETL、IET这样的在FA决定处理中需要的每个參数计算FA变化形状FAP，将计算结果作为FA变化形状数据库存储在存储装置72或外部存储装置73中。
[0137]在该情况下，序列作成部300也可以不具备FA决定部310。即，序列作成部300当设定了摄像条件吋，从FA变化形状数据库中抽出与之对应的FA变化形状FAP，使用它作成摄像序列。
[0138]在图10中表示该情况下的摄像处理的流程。在此，与图4相同的处理附加相同的符号。如本图所示那样，当从操作者输入了摄像条件时，信号处理系统70接受它(步骤SI 101)。然后，序列作成部300从FA变化形状数据库中抽出与接受的摄像条件对应的FA变化形状FAP (步骤S1112)。然后，序列作成部300使用抽出的FA变化形状FAP、摄像条件，作成摄像序列(步骤S1113)。信号处理系统70使用作成的摄像序列，向序列产生器40给出指令，执行摄像(步骤S1104)。[0139] 此外，也可以构成为在步骤S1112中在FA变化形状数据库中没有与接受的摄像条件对应的FA变化形状FAP的情况下，选择与接受的摄像条件最接近的FA变化形状FAP，或者还可以进行内插。
[0140]此外，在该情况下，也可以构成为在步骤S1112后，修正所抽出的FA变化形状。
[0141]另外，在上述实施方式中，构成为由操作者直接输入Tl、T2作为摄像条件，但操作者指定的信息并不限于此。只要是能够确定Tl、T2的信息即可。例如也可以构成为代替Tl、T2而由操作者设定对象组织。在该情况下，将每个组织的Tl、T2值保存在存储装置72或外部存储装置73中。如果指定了对象组织，则序列作成部300从存储装置72等中抽出与之对应的Tl、T2值，使用它决定FA变化形状FAP，作成摄像序列。
[0142]在该情况下，如图5 (b)所示，UI画面400在摄像条件设定区域410中具备对象组织输入栏403，来取代Tl输入栏401和T2输入栏402。
[0143]此外，在有多个对象组织的情况下，也可以采用摄像对象的全部组织的Tl和T2各自的平均值。另外，也可以使用每个组织的Tl、T2，对每个组织计算上述FA变化形状FAP，对每个组织的FA变化形状FAP进行平均，设为最終的FA变化形状FAP。在计算平均吋，也可以根据摄像对象内部的组织的分布、比例进行加权。
[0144]图11表示使用各个Tl和T2计算每个组织的FA变化形状FAP、然后取得平均时的处理流程。在此，列举组织有K种(I(是自然数)的第一实施方式的情况为例子进行说明。
[0145]如果从操作者输入了摄像条件，则本实施方式的信号处理系统70接受它(步骤S1101)。
[0146]然后，使用所接受的摄像条件，序列作成部300使FA决定部310决定各再收敛RF脉冲的FA。这时，在本变形例子中，FA决定部310对输入的每个Tl和T2进行FA决定处理，计算各自的FA变化形状(步骤S1121、S1102、S1122、S1123)。
[0147]然后，序列作成部300计算每个T1、T2的FA变化形状的平均，将它设为FA变化形状(步骤S1124)。
[0148]然后，序列作成部300使用计算出的各再收敛RF脉冲的FA (FA变化形状)、摄像条件，作成摄像序列(步骤S1103)。信号处理系统70使用作成的摄像序列，向序列产生器40给出指令，执行摄像(步骤S1104)。
[0149]另外，也可以通过预扫描来求出T1、T2。在图12中表示该情况下的摄像处理的流程。本处理基本上与第一实施方式的摄像处理相同，但如果接受了摄像条件，则信号处理系统70进行预扫描(步骤SI 131)，决定Tl、T2 (步骤SI 132)，转移到FA决定处理SI 102。
[0150]此外，在图11和图12所示的变形例子中，还可以从预先保存在存储装置72等中的数据库中，抽出与Tl和T2对应的FA变化形状，来代替FA决定处理。
[0151]〈く第二实施方式〉〉
[0152]接着，说明应用本发明的第二实施方式。在第一实施方式中，使全部再收敛RF脉冲的FA在可能的范围内变化，探索使SNR最大(使噪声放大率NA最小)的FA变化形状。即，探索的參数(探索參数)个数是在一次TR期间施加的全部再收敛RF脉冲数(回波数)N。例如，也可以把与FA之间能够相互变换的信号強度SS作为探索參数，在该情况下，探索參数个数为N。在本实施方式中，使该探索參数个数比N少。
[0153]本实施方式的MRI装置100、序列作成部300的结构基本上与第一实施方式相同。另外，本实施方式的信号处理系统70的摄像处理的流程也与第一实施方式相同。因此，省略对它们的说明。以下，主要着眼于与第一实施方式不同的结构来说明本实施方式。
[0154]在本实施方式中，如上述那样准备回波数N以下的L个(L〈N的自然数)探索參数(总称为CTRL)，使用它决定使SNR最大化的FA变化形状FAP。此外，使用EPG、ProspectiveEPG等，能够与信号强度和FA相互进行变换。因此，探索參数既可以表示FA，也可以表示信号強度。此外，ProspectiveEPG与EPG相反，是根据信号强度变化形状SSP计算FA变化形状FAP的技术。
[0155]在ProspectiveEPG中，通过从信号强度变化形状SSP的开头顺序地求出实现信号強度SS的FA，计算FA变化形状FAP。在该过程中也可能不存在能够实现信号強度SS那样的FA、或不存在于设定的范围内。将这样的状态称为FA已饱和。 [0156]在FA饱和的情况下,通过依照预先确定的规则决定FA,来计算FA变化形状FAP。例如，进行将饱和以后的FA全部设为预先确定的上限值FAH等的处理。
[0157]本实施方式的FA决定部320如图13所示那样，具备决定使SNR最大化的探索參数CTRL的參数决定部321、根据探索參数CTRL计算各FA值(FA变化形状FAP)的FA计算部 322。
[0158]以下，列举具体例子说明本实施方式的FA决定处理。在此，使用隔I个抽取构成FA变化形状FAP的N个各FA值所得的FA值列作为探索參数CTRL。将探索參数CTRL的个数设为L (L=N/2:N是偶数，L= (N+D/2:N是奇数)。
[0159]FA计算部322根据探索參数CTRL计算FA变化形状FAP。例如，在探索參数CTRL是间隔抽取构成FA变化形状FAP的N个FA值而作成的情况下，使用线性内插等方法计算FA变化形状FAP。
[0160]參数决定部321使用FA计算部322将探索參数CTRL变换为FA变化形状FAP，通过与第一实施方式相同的探索方法，决定SNR最大、即噪声放大率NA最小的探索參数CTRL。
[0161]在图14中表示本实施方式的FA决定处理的流程。在本实施方式中，也与第一实施方式同样地，说明使用最急下降法重复M次进行值的更新的情況。
[0162]首先，參数决定部321向针对重复次数进行计数的计数器i设定I (步骤S2101)。然后，对探索參数CTRL设定初始值(步骤S2102)。与第一实施方式同样地，例如将初始值全部设为0度。另外，參数决定部321使用FA计算部322，根据所得到的探索參数CTRL，作成FA变化形状FAP (步骤S2103)。
[0163]接着，參数决定部321计算噪声放大率NA对探索參数CTRL的偏微分(在计算中为差商)(步骤S2104)。在此，与第一实施方式同样地，使用IET、T1、T2，根据在步骤S2103中计算出的FA变化形状FAP，通过EPG等计算信号强度变化形状SSP。然后，使用计算出的信号强度变化形状SSP，计算噪声放大率NA及其偏微分。
[0164]然后，通过与第一实施方式相同的方法更新探索參数CTRL (步骤S2105)。即，从更新前的探索參数CTRL的各值减去将在步骤S2104中得到的偏微分的值乘以预先确定的正的整数所得的值，得到更新后的探索參数CTRL。
[0165]然后，參数决定部321判别是否重复进行了 M次处理(步骤S2106)，如果不满M次，则将计数器i加I (步骤S2107)，返回到步骤S2103。
[0166]另ー方面，在步骤S2106中重复进行了 M次的情况下，參数决定部321将该时刻的探索參数CTRL决定为使SNR最大的探索參数CTRL (步骤S2108)。然后，FA决定部320使FA计算部322根据在步骤S2108中決定的探索參数CTRL计算FA变化形状(步骤S2109)，
结束处理。
[0167]此外，在本实施方式中，探索噪声放大率NA的最小值的方法也不限于此。例如也可以与第一实施方式同样地，不通过最急下降法，而是整体地将根据得到对变化的每个探索參数CTRL计算出的噪声放大率NA中的最小的NA的探索參数CTRL求出的FA变化形状FAP作为处理結果。
[0168]另外，也根据重复次数来进行收敛判定，但并不限于此。例如，也可以与第一实施方式同样地，根据噪声放大率NA的偏微分的值来判定。
[0169]如以上说明的那样，根据本实施方式，与第一实施方式同样地，在进行信号修正等，配置在k空间的中心以外的回波信号对SNR产生影响的情况下，也能够在多回波摄像中决定使SNR最大的FA变化形状，使用该FA变化形状进行摄像，因此能够得到考虑了信号修正的SNR大的图像。
[0170]另外，本实施方式与第一实施方式相比，探索所用的參数少。因此，与第一实施方式相比，能够通过少的计算时间实现最优的SNR。
[0171]此外，降低參数个数的方法并不限于上述方法。例如，作为探索參数，既可以使用表示FA变化形状FAP的曲线的近似式的系数，也可以使用定义表示信号强度变化形状SSP的曲线的公式的系数。如果使用EPG、ProspectiveEPG等，则能够对信号强度和FA相互变换，因此探索參数既可以表示FA，也可以表示信号強度。
[0172]另外，在本实施方式中，也可以与第一实施方式同样地，预先确定在最优化法中使FA变化的范围来进行限制。另外`，也可以具备修正由FA决定部320決定的FA变化形状FAP的FA修正部311。
[0173]另外，还可以将预先作成的FA变化形状FAP保存在存储装置72等中，在摄像中从其中抽出。
[0174]并且，还可以设定摄像对象组织来代替T1、T2。另外，在摄像对象组织有多种的情况下，也可以使用T1、T2值的平均值计算FA变化形状，还可以对每个T1、T2值计算FA变化形状，计算FA变化形状的平均。另外，也可以通过预扫描得到Tl、T2。
[0175]另外，在本实施方式中，在步骤S2103中，根据探索參数CTRL计算FA变化形状FAP、信号强度变化形状SSP，使用信号强度变化形状SSP计算噪声放大率NA的偏微分，但并不限于该方法。例如，在如上述实施方式那样间隔抽取FA变化形状来作成的情况下，也可以不计算没有间隔抽取的FA变化形状而更新探索參数CTRL。即，只根据作为探索參数CTRL被抽出的FA值列计算间隔抽取出的信号强度值列。另外，根据该信号强度值列计算噪声放大率NA及其偏微分。更新探索參数CTRL。
[0176]〈く第三实施方式〉〉
[0177]接着，说明应用本发明的第三实施方式。在上述各实施方式中，作为具体例子，控制直接与FA变化形状有关的參数使噪声放大率最小化。
[0178]在本实施方式中，使信号强度变化形状变化，决定使噪声放大率最小化的信号强度变化形状，将实现该信号强度变化形状的FA变化形状作为在摄像中使用的FA变化形状。即，作为第二实施方式的探索參数，使用确定信号强度变化形状的信息。[0179]本实施方式的MRI装置100、序列作成部300具有与第二实施方式基本相同的结构。另外，本实施方式的信号处理系统70的摄像处理的流程也与第二实施方式相同。因此，在此省略它们的说明。以下，主要着眼于与第二实施方式不同的结构来说明本实施方式。
[0180]在本实施方式中，如上述那样，作为探索參数，使用确定信号强度变化形状SSP的信息。使该探索參数变化，决定使SNR最大化的信号强度变化形状SSP，根据所決定的信号强度变化形状SSP决定FA变化形状FAP。为了实现它，本实施方式的FA决定部330如图15所示，具备决定用于确定使SNR最大化的信号强度变化形状SSP的信息的信号强度变化形状决定部331、根据确定作为探索參数的信号强度变化形状SSP的信息计算FA变化形状FAP的FA计算部332。
[0181]信号强度变化形状决定部331用预定的函数表示信号强度变化形状SSP，决定作为探索參数的函数的參数(系数等)使得噪声放大率NA最小。这时，与第一实施方式、第二实施方式同样地，使用最急下降法等最优化法决定该函数的參数。
[0182]在调查了若干种情况后可知，信号强度变化形状SSP的一部分或全部是渐近地逼近指数函数的形状时，噪声放大率NA大致为最小。在本实施方式中，作为信号強度SS (t)的函数，使用以下的公式(I)表示的二重指数函数。
[0183]
【权利要求】
1.一种磁共振成像装置，具备产生静磁场的静磁场产生系统、向配置在上述静磁场中的被检测体施加倾斜磁场的倾斜磁场产生系统、发送以预定的翻转角激励上述被检测体磁化的高频磁场脉冲的高频磁场发送系统、接收上述被检测体产生的回波信号的信号接收系统、根据上述信号接收系统接收到的回波信号重新构成图像并且依照摄像序列控制上述倾斜磁场产生系统、上述高频磁场发送系统、上述信号接收系统的动作的控制系统， 上述磁共振成像装置的特征在于， 上述控制系统具备: 接受部，其从用户接受决定上述摄像序列的摄像条件； 序列生成部，其依照上述接受部接受的摄像条件和脉冲序列生成上述摄像序列， 上述脉冲序列是在施加I个激励高频磁场脉冲后施加多个再收敛高频磁场脉冲的序列， 上述序列生成部根据上述信号接收系统接收的多个回波信号的信号强度，决定上述各再收敛闻频磁场脉冲的翻转角以使上述图像的SNR最大。
2.根据权利要求1所述的磁共振成像装置，其特征在于， 上述序列生成部使用反映上述多个回波信号的信号修正后的图像的SNR的指标，决定使该图像的SNR最大的上述各再收敛闻频磁场脉冲的翻转角。
3.根据权利要求2所述的磁共振成像装置，其特征在于， 上述序列生成部具备:翻转角决定部，其根据上述摄像条件决定上述各再收敛高频磁场脉冲的翻转角，` 上述翻转角决定部通过依照预先确定的最优化法使确定上述各再收敛高频磁场脉冲的翻转角的信息变化，重复计算上述指标，来决定上述翻转角。
4.根据权利要求3所述的磁共振成像装置，其特征在于， 上述翻转角决定部具备:信号强度计算部，其根据确定上述各再收敛高频磁场脉冲的翻转角的信息，分别计算通过施加与各信息对应的再收敛高频磁场脉冲而得到的回波信号的信号强度， 作为上述指标，使用根据计算出的上述各回波信号的信号强度的倒数的平方的平均的平方根得到的噪声放大率， 决定上述各再收敛高频磁场脉冲的翻转角以使上述噪声放大率最小。
5.根据权利要求3所述的磁共振成像装置，其特征在于， 确定上述各再收敛高频磁场脉冲的翻转角的信息是上述各再收敛高频磁场脉冲的翻转角。
6.根据权利要求3所述的磁共振成像装置，其特征在于， 确定上述各再收敛高频磁场脉冲的翻转角的信息由比上述各再收敛高频磁场脉冲数少的数构成。
7.根据权利要求3所述的磁共振成像装置，其特征在于， 确定上述各再收敛高频磁场脉冲的翻转角的信息是依照预定的规则间隔抽取上述各再收敛高频磁场脉冲的翻转角所得到的中间值列， 上述翻转角决定部具备: 中间值列决定部，其使上述中间值列变化来决定使上述SNR最大的中间值列；翻转角计算部，其对上述决定的中间值进行内插，计算上述各再收敛高频磁场脉冲的翻转角。
8.根据权利要求3所述的磁共振成像装置，其特征在于， 确定上述各再收敛高频磁场脉冲的翻转角的信息是表示通过施加上述各再收敛高频磁场脉冲而得到的回波信号的信号强度的变化形状的函数的参数， 上述翻转角决定部具备:信号强度变化形状决定部，其使上述参数变化来决定使上述SNR最大的上述参数，决定上述信号强度的变化的形状； 翻转角计算部，其根据用上述参数确定的信号强度的变化的形状，计算上述各再收敛高频磁场脉冲的翻转角。
9.根据权利要求8所述的磁共振成像装置，其特征在于， 上述函数是渐近地逼近指数函数的函数， 上述参数是上述指数函数的系数以及时间常数中的至少一方。
10.根据权利要求8所述的磁共振成像装置，其特征在于， 上述函数是二重指数函数， 上述参数是构成上述二重指数函数的2个指数函数各自的系数和时间常数中的一个以上。
11.根据权利要求8所述的磁共振成像装置，其特征在于， 预先确定使上述翻转角变化`的范围， 上述信号强度变化形状决定部使上述参数变化，以使上述翻转角收敛在上述范围内并且不饱和。
12.根据权利要求10所述的磁共振成像装置，其特征在于， 作为上述指标，使用根据计算出的上述各回波信号的信号强度的倒数的平方的平均的平方根得到的噪声放大率， 上述二重指数函数由上述信号强度的变化形状渐近地逼近的第一指数函数、表示渐近的变化的第二指数函数构成， 上述参数是上述第一指数函数的时间常数， 上述信号强度变化形状决定部在固定上述噪声放大率来使上述第一指数函数的时间常数变化时把最后的再收敛高频磁场脉冲的上述翻转角最小的上述第一指数函数的时间常数决定为使上述SNR为最大的上述参数。
13.根据权利要求8所述的磁共振成像装置，其特征在于， 预先确定上述翻转角的范围， 上述翻转角决定部在计算出的翻转角不存在于上述范围内的情况下，依照预先确定的规则决定翻转角。
14.根据权利要求3~13的任意一项所述的磁共振成像装置，其特征在于， 上述序列生成部还具备:翻转角修正部，其依照预先确定的规则修正由上述翻转角决定部决定的翻转角。
15.根据权利要求1或2所述的磁共振成像装置，其特征在于， 还具备:翻转角存储部，其对每个摄像条件存储作为上述各再收敛高频磁场脉冲的翻转角的使上述图像的SNR最大而决定的翻转角，上述序列生成部使用与上述接受部接受的摄像条件对应地存储在上述翻转角存储部中的翻转角来作为上述翻转角。
16.根据权利要求1~13的任意一项所述的磁共振成像装置，其特征在于， 在上述摄像条件中包含用于确定成为摄像对象的核素的纵向迟豫时间和横向迟豫时间的信息。
17.一种翻转角决定方法，其决定在磁共振成像装置中使用的多回波摄像序列的多个再收敛高频磁场脉冲的翻转角， 该翻转角决定方法的特征在于，具有: 使用反映在上述磁共振成像装置中取得的多个回波信号的信号修正后的图像的SNR的指标， 依照预先确定的最优化法，使确定上述各再收敛高频磁场脉冲的翻转角的信息变化，重复计算上述指标，决定上述图像的SNR成为最大的翻转角的翻转角决定步骤。
【文档编号】A61B5/055GK103561645SQ201280024917
【公开日】2014年2月5日 申请日期:2012年5月23日 优先权日:2011年6月9日
【发明者】黑川真次 申请人:株式会社日立医疗器械
相关技术
网友询问留言
已有0条留言
- 还没有人留言评论。精彩留言会获得点赞！
1