CN111537930B - 基于梯度波形调整的磁共振参数成像方法与设备 - Google Patents

Abstract

本申请适用于磁共振技术领域，提供了一种基于梯度波形调整的磁共振参数成像方法与设备，其中，一种基于梯度波形调整的磁共振参数成像方法，通过获取预先设定的扫描序列参数与控制参数，根据扫描序列参数与控制参数，可以确定待调整信号梯度波形的目标面积值，其中，待调整信号梯度波形包括平台波段与渐变波段，按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，对渐变波段进行平滑调整得到新的渐变波段，降低了原有渐变波段所表征的信号拐点的落差，无需对梯度场切换过程中产生的噪音进行采集，再匹配相应的降噪信号进行降噪操作，简化了降噪方案且节约了降噪成本。

Description

基于梯度波形调整的磁共振参数成像方法与设备

技术领域

本申请属于磁共振技术领域，尤其涉及一种基于梯度波形调整的磁共振参数成像方法、装置、磁共振成像设备及计算机可读存储介质。

背景技术

磁共振成像技术(Magnetic Resonance Imaging，MRI)是通过对静磁场中的目标体施加特定频率的射频脉冲，使该目标体中的氢质子受到激励而发生磁共振现象。在停止施加射频脉冲后，目标体中的质子在弛豫过程中产生磁共振MR信号，通过对该MR信号的接收、空间编码和图像重建等处理过程，即可获得可用的MR信号。由于磁共振过程中，每一个信号都含有全层的信息，因此需要对磁共振信号进行空间定位编码，即频率编码和相位编码。具体地，因为接收线圈采集到的MR信号实际是带有空间编码信息的无线电波，属于模拟信号，所以需要经过模数转换变成数字信息，再将数字信息填充到K空间，最后得到相应的数字点阵。其中，K空间与磁共振信号的空间定位息息相关，该K空间也叫傅里叶空间，是带有空间定位编码信息的MR信号原始数字信息的填充空间。每一幅MR图像都有其相应的K空间数据点阵，通过对K空间的数据进行傅里叶转换，就能对原始数字数据中的空间定位编码信息进行解码，分解出不同频率、相位和幅度的MR信号，不同的频率和相位代表不同的空间位置，而幅度则代表MR信号强度。把不同频率、相位及信号强度的MR数字信号分配到相应的像素中，就得到MR图像数据，也即重建出了MR图像。傅里叶变换就是把K空间的原始数据点阵转变成MR图像点阵的过程。

在使用磁共振成像技术进行扫描成像的过程中，需要在静磁场环境中施加一个梯度场，该梯度场用于磁共振成像过程中配合射频脉冲的激励实现成像区域的选择，以及对成像的目标体上产生的MR信号进行空间位置的编码。现有的磁共振成像设备中，用于制造梯度场环境的梯度场组件中，通过三组梯度线圈不断地在接通、断开状态下切换，并在空间中，如坐标X、Y、Z三个方向上产生梯度磁场，从而构建出梯度场环境。由于磁共振成像设备在工作过程中，梯度场组件必须快速地到达最大功率，而且相位编码梯度与层面选择梯度必须在读出梯度接通前快速地断开。有时梯度场的极性还会快速地切换，而在快速切换过程中，梯度线圈中的金属丝之间会产生剧烈地震动，进而产生较大的噪声。

虽然现有技术中在解决磁共振成像过程中的噪音问题时，可以通过使用带阻滤波器来抑制梯度波形声压级高的频带成分，但是该方案需要采集梯度线圈中的声压，然后计算频率响应函数进而滤除梯度线圈在切换梯度场时说产生的特定频段声音。由此可见，现有的磁共振成像技术中，在降低梯度场切换过程中的噪音时，存在降噪方案较复杂且成本较高的问题。

发明内容

有鉴于此，本申请实施例提供了一种基于梯度波形调整的磁共振参数成像方法、装置、磁共振成像设备及计算机可读存储介质，以解决现有的磁共振成像技术中，在降低梯度场切换过程中的噪音时，存在降噪方案较复杂且成本较高的问题。

本申请实施例的第一方面提供了一种基于梯度波形调整的磁共振参数成像方法，包括：

获取预先设定的扫描序列参数与控制参数；其中，所述控制参数用于描述待调整信号梯度波形；

根据所述扫描序列参数与所述控制参数，确定所述待调整信号梯度波形的目标面积值；其中，所述待调整信号梯度波形包括平台波段与渐变波段；

按照预设幅度调整参数调整所述平台波段的波形幅度，得到新的平台波段；

基于所述目标面积值与所述新的平台波段的第一面积值，对所述渐变波段进行平滑调整，得到新的渐变波段；其中，所述新的渐变波段的第二面积值与所述第一面积值之和等于所述目标面积值；

基于所述新的渐变波段与所述新的平台波段组成的目标信号梯度波形，控制所述梯度场。

进一步的，所述控制参数包括用于描述所述待调整信号梯度波形的梯度函数；

所述扫描序列参数包括：K空间尺寸、K空间单元尺寸、原子核的旋磁比、扫描视野、带宽、采样时间，以及与所述采样时间关联的采样点个数。

进一步的，所述根据所述扫描序列参数与所述控制参数，确定所述待调整信号梯度波形的目标面积值，包括：

通过以下公式测算出所述待调整信号梯度波形的目标面积值；

k＝N·Δk；

其中，k(t)为所述采样时间为t时刻的K空间位置；γ为所述原子核的旋磁比；G(t′)为所述梯度函数；k为所述K空间尺寸；N为所述采样点个数；Δk为所述K空间单元尺寸；FOV为所述扫描视野；A为所述目标面积值；BW为所述带宽。

进一步的，所述按照预设幅度调整参数调整所述平台波段的波形幅度，得到新的平台波段，包括：

根据所述梯度函数确定所述平台波段的平台幅度值；

按照预设幅度调整参数调整所述平台幅度值，得到新的平台幅度值；其中，所述新的平台波段的幅度值等于所述平台幅度值与所述调整参数之和；

根据所述新的平台幅度值得到所述新的平台波段。

进一步的，所述采样时间包括平台波段持续时间与渐变波段持续时间；所述采样点个数包括与所述平台波段持续时间对应的第一采样点个数，以及与所述渐变波段持续时间对应的第二采样点个数；

所述基于所述目标面积值与所述新的平台波段的第一面积值，对所述渐变波段进行平滑调整，得到新的渐变波段，包括：

获取所述第一采样点个数；

将所述第一采样点个数与所述新的平台幅度值的乘积，识别为所述第一面积值；

测算所述目标面积值与所述第一面积值之差，得到调整面积值；

基于所述调整面积值对所述渐变波段进行平滑调整，得到新的渐变波段。

进一步的，所述渐变波段包括在所述渐变波段持续时间内连续的多个渐变点；

所述基于所述调整面积值对所述渐变波段进行平滑调整，得到新的渐变波段，包括：

获取所述渐变波段持续时间；

基于所述渐变波段持续时间确定多个渐变点；

通过以下公式调整每个所述渐变点的幅度值，得到多个新的渐变点；

X(t)＝1-exp(-w*t)；

G(t)＝G0*Xn(t)；

其中，X(t)为所述新的渐变点；t为所述渐变波段持续时间中的时刻；w为调节因子，且0<|w|<1；Xn(t)表征X(t)的归一化结果；G0为所述渐变点的幅度值；G(t)为所述新的渐变点的幅度值；多个所述新的渐变点组成新的渐变波段，且多个所述新的渐变点的幅度值之和等于所述调整面积值。

进一步的，所述方法还包括：

确定与所述目标信号梯度波形对应的目标梯度函数；

将所述目标梯度函数与所述控制参数关联存储至预设数据库中。

本申请实施例的第二方面提供了一种基于梯度波形调整的磁共振参数成像装置，包括：

第一获取单元，用于获取预先设定的扫描序列参数与控制参数；其中，所述控制参数用于描述待调整信号梯度波形；

第一确定单元，用于根据所述扫描序列参数与所述控制参数，确定所述待调整信号梯度波形的目标面积值；其中，所述待调整信号梯度波形包括平台波段与渐变波段；

第一调整单元，用于按照预设幅度调整参数调整所述平台波段的波形幅度，得到新的平台波段；

第二调整单元，用于基于所述目标面积值与所述新的平台波段的第一面积值，对所述渐变波段进行平滑调整，得到新的渐变波段；其中，所述新的渐变波段的第二面积值与所述第一面积值之和等于所述目标面积值；

执行单元，用于基于所述新的渐变波段与所述新的平台波段组成的目标信号梯度波形，控制所述梯度场。

本申请实施例的第三方面提供了一种磁共振成像设备，包括存储器、处理器以及存储在所述存储器中并可在所述磁共振设备上运行的计算机程序，所述处理器执行所述计算机程序时实现第一方案提供的基于梯度波形调整的磁共振参数成像方法的各步骤。

本申请实施例的第四方面提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现第一方案提供的基于梯度波形调整的磁共振参数成像方法的各步骤。

本申请实施例的第五方面提供了一种计算机程序产品，当计算机程序产品在磁共振设备上运行时，使得磁共振设备执行上述第一方面中任一项所述的基于梯度波形调整的磁共振参数成像方法的各步骤。

实施本申请实施例提供的一种基于梯度波形调整的磁共振参数成像方法、装置、磁共振成像设备及计算机可读存储介质具有以下有益效果：

本申请实施例提供的一种基于梯度波形调整的磁共振参数成像方法，通过获取预先设定的扫描序列参数与控制参数，由于控制参数能够用于描述待调整信号梯度波形，因此根据扫描序列参数与控制参数，可以确定待调整信号梯度波形的目标面积值，且待调整信号梯度波形包括平台波段与渐变波段，由于该渐变波段所表征的信号拐点落差越大，梯度场切换过程中的噪音就越大，因此按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，对渐变波段进行平滑调整得到新的渐变波段，进而降低了原有渐变波段所表征的信号拐点的落差，最后基于新的渐变波段与新的平台波段组成的目标信号梯度波形控制梯度场，即可实现在磁共振成像过程中对切换梯度场所导致的噪音进行降噪处理，无需对梯度场切换过程中产生的噪音进行采集，再匹配相应的降噪信号进行降噪操作，简化了降噪方案且节约了降噪成本。

此外，按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，调整渐变波段得到新的渐变波段，使得新的渐变波段与新的平台波段组成的目标信号梯度波形，保持有待调整信号梯度波形的目标面积值不变的情况下，尽可能保留平台波段，也即实现最大化地保留控制参数同时，降低梯度切换剧烈带来的噪声问题。

附图说明

为了更清楚地说明本申请实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本申请实施例提供的一种基于梯度波形调整的磁共振参数成像方法的实现流程图；

图2是本申请整体方案实现原理示意图；

图3是本申请另一实施例提供的一种基于梯度波形调整的磁共振参数成像方法的实现流程图；

图4是本申请实施例提供的一种基于梯度波形调整的磁共振参数成像装置结构框图；

图5是本申请另一实施例提供的一种磁共振成像设备的结构框图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。

请参阅图1，图1是本申请实施例提供的一种基于梯度波形调整的磁共振参数成像方法的实现流程图。本实施例中，基于梯度波形调整的磁共振参数成像方法用于在磁共振成像过程中的梯度场切换控制，其执行主体为磁共振成像设备。

如图1所示的基于梯度波形调整的磁共振参数成像方法包括以下步骤：

S11：获取预先设定的扫描序列参数与控制参数；其中，所述控制参数用于描述待调整信号梯度波形。

在步骤S11中，扫描序列参数用于描述采集窗口，也即进行数据采集的时间段。由于采集窗口与控制参数在时序上存在关联关系，也即控制参数所描述的待调整信号梯度波形与采集窗口之间在时序上存在重叠区域，因此该重叠区域为数据采集的时间段。控制参数为控制梯度场的控制参数，也即控制磁共振成像过程中梯度场切换的控制参数，因此该控制参数所描述的梯度波形为待调整信号梯度波形。

需要说明的是，在磁共振成像过程中，通过构建一个静磁场环境，再对静磁场中的目标体如人体，施加某种特定频率的射频脉冲，使目标体中的氢质子受到激励而发生磁共振现象。在静磁场环境下，还需要构建一个梯度场，该梯度场用于在磁共振成像过程中做选择激励以选择成像区域，以及对成像的目标体上产生的MR信号进行空间位置的编码。现有的磁共振成像设备中，用于制造梯度场环境的梯度场组件中，通过三组梯度线圈不断地在接通、断开状态下切换，并在空间中产生梯度磁场。在本申请的所有实施例中，控制参数所描述的待调整信号梯度波形，是用于控制梯度磁场切换的信号波形。

图2示出了整体方案实现原理示意图。如图2所示，在实际应用中，由于控制参数所呈现出来的信号梯度波形存在较为剧烈的拐点，如图2中的点P1与点P2，而该剧烈拐点所对应的参数信号就是导致噪音的原因，因此为了降低梯度场切换过程中的噪音，需要对该待调整信号梯度波形中的剧烈拐点进行调整。

在本实施例中，获取扫描序列参数与控制参数，由于控制参数能够用于描述待调整信号梯度波形，因此将该控制参数进行图像化，能够更多好地对控制参数进行分析和调整，也是为对整体方案进行改进提供途径和基础。应当理解的是，在本申请的所有实施例中，对待调整信号梯度波形进行调整，实际上是对控制参数进行调整。

至于何时获取预先设定的扫描序列参数与控制参数，可以包括但不仅限于以下两个场景。

场景1：若在进行磁共振成像的过程中，检测到用于选择成像区域的预设指令时，获取预先设定的扫描序列参数与控制参数。

例如，在对目标体进行磁共振成像的过程中，当通过选择不同的信号激励进而对目标体的成像区域进行选择时，获取预先设定的扫描序列参数与控制参数。

场景2：若检测到用于对控制参数进行调整的预设指令，则获取预先设定的扫描序列参数与控制参数。

例如，在向磁共振成像设备导入控制方案后，通过触发用于对控制参数进行调整的预设指令，进而获取预先设定的扫描序列参数与控制参数。

应当理解的是，在实际应用中，由于磁共振成像过程的自动化程度较高，因此对控制参数进行调整，可以是在使用磁共振成像设备的过程中进行调整，或者是在对磁共振成像设备进行调试时进行调整。预先设定的扫描序列参数与控制参数，可以预先存储在磁共振成像设备的数据库中，获取扫描序列参数与控制参数时，可以根据控制指令所对应的磁共振成像策略，从数据库中获取相应的扫描序列参数与控制参数。

S12：根据所述扫描序列参数与所述控制参数，确定所述待调整信号梯度波形的目标面积值；其中，所述待调整信号梯度波形包括平台波段与渐变波段。

在步骤S12中，待调整信号梯度波形能够用于表征控制参数与时间之间的对应关系，也即不同的时间点对应不同的幅度值。连续的时间段内，多个离散的采样点的幅度值在坐标上相连即形成了待调整信号梯度波形，且各幅度值在连续的时间段内进行累加，即可得到该待调整信号梯度波形的目标面积值。

需要说明的是，扫描序列参数用于描述采集窗口，该采集窗口与待调整信号梯度波形在时序上存在重叠，而待调整信号梯度波形与采集窗口在时序上重叠的区域，用于表征有效的数据采集时段和数据内容，也即用于表征采集到的数据在K空间上的位置。

在本实施例中，根据扫描序列参数与控制参数确定待调整信号梯度波形的目标面积值，该目标面积值并非待调整信号梯度波形的完整面积值，而是待调整信号梯度波形与采集窗口在时序重叠部分的面积值，也即待调整信号梯度波形的中有效的数据采集波段。

如图2所示，待调整信号梯度波形的平台波段a的时序与扫描序列参数所描述的采集窗口a’的时序完全重叠，也即采集窗口所表征的有效采集数据和有效采集时段，与待调整信号梯度波形的平台波段a对应。

如图2所示，在本实施例中，待调整信号梯度波形201中的平台波段a与渐变波段(a1、a2)是由采样时间进行区分，也即渐变波段(a1、a2)对应的采样时间为时刻1至时刻2，以及时刻3至时刻4，平台波段a对应的采样时段为时刻2至时刻3。

如图2所示，由于时刻0至时刻1属于准备时刻，时刻4至时刻5为结束时刻，因此不纳入采样窗口范围内。由于平台波段a的前后两个波段a1与a2，都包含有两个临界值点，也即变化剧烈的拐点P1与拐点P2，因此在对待调整信号梯度波形进行调整时，波段a1与波段a2被识别为渐变波段，平台波段则为a。

可以理解的是，无论是中平台波段还是渐变波段，都可以是由点组成的线。在待调整信号梯度波形中，渐变波段则由平台波段的前后两段包含有临界值端点的波段所组成，也即可以将渐变波段视为平台波段的前后两段波段，且每段渐变波段中都包含有临界值端点。

在实际应用中，待调整信号梯度波形的具体形状取决于控制参数，也即控制参数是一系列梯度幅度值与时间之间的关系数据，在坐标系中，可以将横坐标代表时间，纵坐标代表梯度幅度值，即可绘制出待调整信号梯度波形。由于控制参数可以在坐标系中绘制出相应的待调整信号梯度波形，因此控制参数中的部分数据可以用描述待调整信号梯度波形的梯度函数所代替。

作为本实施例一种可能实现的方式，控制参数包括用于描述所述待调整信号梯度波形的梯度函数；扫描序列参数包括：K空间尺寸、K空间单元尺寸、原子核的旋磁比、扫描视野、带宽、采样时间，以及与所述采样时间关联的采样点个数。步骤S12具体可以包括：

通过以下公式测算出所述待调整信号梯度波形的目标面积值；

k＝N·Δk；

其中，k(t)为所述采样时间为t时刻的K空间位置；γ为所述原子核的旋磁比；G(t′)为所述梯度函数；k为所述K空间尺寸；N为所述采样点个数；Δk为所述K空间单元尺寸；FOV为所述扫描视野；A为所述目标面积值；BW为所述带宽。

在本实施例中，K空间尺寸为采样点个数K空间单元尺寸的乘积，K空间单元尺寸Δk则与扫描视野FOV之间为倒数关系；N为采样点个数，且N为大于0的整数。

需要说明的是，K空间也叫傅里叶空间，是带有空间定位编码信息的MR信号原始数字数据的填充空间，每一幅MR图像都有其相应的K空间数据点阵，对K空间的数据进行傅里叶转换，就能对原始数字数据中的空间定位编码信息进行解码，分解出不同频率、相位和幅度的MR信号，不同的频率和相位代表不同的空间位置，而幅度则代表MR信号强度，也即傅里叶变换就是把K空间的原始数据点阵转变成磁共振图像点阵的过程。由于待调整信号梯度波形与采集窗口在时序上重叠的区域，用于表征有效的数据采集时段和数据内容，也即用于表征采集到的数据在K空间上的位置，通过上述公式，结合扫描序列参数中所包括的K空间尺寸、K空间单元尺寸、原子核的旋磁比、扫描视野、带宽、采样时间，以及与采样时间关联的采样点个数之间换算关系，即可确定出待调整信号梯度波形的目标面积值。

在本申请的所有实施例中，由于待调整信号梯度波形的目标面积值用于表征采集到的数据在K空间上的位置，因此为了确保数据在K空间上的位置保持不变，在对待调整信号梯度波形进行调整，所得到的新的信号梯度波形的面积值与待调整信号梯度波形的目标面积值相等。

S13：按照预设幅度调整参数调整所述平台波段的波形幅度，得到新的平台波段。

在步骤S13中，预设幅度调整参数可以为平台波段的波形幅度的增量或者倍数，也即调整平台波段的波形幅度，是在平台波段原有的波形幅度的基础上增加其波形幅度值的大小。

需要说明的是，当预设幅度调整参数为平台波段的波形幅度的增量时，该增量为大于0的增量数值，当预设幅度调整参数为平台波段的波形幅度的倍数时，该倍数为大于1的倍数。

在本申请的所有实施例中，对待调整信号梯度波形进行调整，具体是对待调整信号梯度波形中具有剧烈拐点的波段进行平滑处理，而为了保证调整后的新的梯形波段所围成的面积值与待调整信号梯度波形的目标面积值相等，因此按照预设幅度调整参数调整平台波段的波形幅度，是为了通过增加平台波段的波形幅度，进而补偿因平滑处理后损失的部分面积。

如图2所示，新的平台波段b的波形幅度大于原有的平台波段a的波形幅度。且由于对待调整信号梯度波形201进行调整，将具有剧烈拐点的波段进行平滑处理后，损失部分面积，通过配置得到新的平台波段b的波形，因为其幅度大于原有的平台波段a的波形幅度，能够补偿因平滑处理后损失的部分面积。

作为本实施例一种可能实现的方式，步骤S13可以包括：

根据所述梯度函数确定所述平台波段的平台幅度值；按照预设幅度调整参数调整所述平台幅度值，得到新的平台幅度值；其中，所述新的平台波段的幅度值等于所述平台幅度值与所述调整参数之和；根据所述新的平台幅度值得到所述新的平台波段。

在本实施例中，新的平台幅度值为原有的平台幅度值与调整参数之和。由于梯度场的采样过程属于离散型采样，也即待调整信号梯度波形中的平台波段与渐变波段分别对应有采样时间，也即采样点，且梯度函数能够用于描述整个待调整信号梯度波形，因此采样时间或采样点个数，即可确定出平台波段和渐变波段，结合梯度函数则可确定出的平台波段的平台幅度值。

需要说明的是，在本申请的所有实施例中，对待调整信号梯度波形进行调整的过程，是先对平台波段进行幅度值调整，再基于调整后得到的新的平台波段，对渐变波段进行调整。由于新的平台波段中，每个采样时间或采样点对应的幅度值都高于原始幅度值，因此新的平台波段与未调整的渐变波段之间存在落差，为了确保调整后梯度波形能够较大程度地保留待调整信号梯度波形的原有特性，因此对渐变波段进行平滑调整，也即使其能够与新的平台波段之间平滑且连贯行程新的信号梯度波形。

S14：基于所述目标面积值与所述新的平台波段的第一面积值，对所述渐变波段进行平滑调整，得到新的渐变波段；其中，所述新的渐变波段的第二面积值与所述第一面积值之和等于所述目标面积值。

在步骤S14中，目标面积值为待调整信号梯度波形的面积。第一面积值为新的平台波段的在时间与幅度坐标轴上的面积，也是新的平台波段的所有采样点的幅度值累加。

需要说明的是，第一面积值为新的平台波段的面积值，也是新的平台波段上每个采样点对应幅度值的累加之和；第二面积值为新的渐变波段的面积值，也是新的渐变波段上每个采样点对应幅度值的累加之和。为了保持待调整信号梯度波形所表征的K空间的数据位置不变，新的平台波段与新的渐变波段组成的新的信号梯度波形，其面积必须与待调整信号梯度波形的目标面积，因此在对渐变波段进行平滑调整时，以新的渐变波段的第二面积值与第一面积值之和等于目标面积，作为限制条件对渐变波段进行平滑调整。

作为本实施例一种可能实现的方式，所述采样时间包括平台波段持续时间与渐变波段持续时间；所述采样点个数包括与所述平台波段持续时间对应的第一采样点个数，以及与所述渐变波段持续时间对应的第二采样点个数；步骤S14可以包括：

获取所述第一采样点个数；将所述第一采样点个数与所述新的平台幅度值的乘积，识别为所述第一面积值；测算所述目标面积值与所述第一面积值之差，得到调整面积值；基于所述调整面积值对所述渐变波段进行平滑调整，得到新的渐变波段。

在本实施例中，由于第一面积值为新的平台波段的在时间与幅度坐标轴上的面积，而新的平台波段上各个点的幅度值均大于原有平台波段上各个点的幅度值，因此为了确保调整后梯度波形所围成很的面积与待调整信号梯度波形的目标面积一致，在对渐变波段进行平滑调整时，需要考虑目标面积与新的平台波段的第一面积值之间的差值。也即对渐变波段进行平滑调整后，新的渐变波段上所有采样点的幅度值的累加之和，等于目标面积与第一面积值之间的差值。

在实际应用中，由于已经获悉了对渐变波段进行平滑调整的限制条件，也即新的渐变波段的第二面积值与第一面积值之和等于目标面积，因此在其他实施例中，可以通过配置相应的条件函数或者方程组，描述该限制条件，进而实现对渐变波段进行平滑调整。

作为本实施例一种可能实现的方式，所述渐变波段包括在所述渐变波段持续时间内连续的多个渐变点；所述基于所述调整面积值对所述渐变波段进行平滑调整，得到新的渐变波段，包括：

获取所述渐变波段持续时间；

基于所述渐变波段持续时间确定多个渐变点；

通过以下公式调整每个所述渐变点的幅度值，得到多个新的渐变点；

X(t)＝1-exp(-w*t)；

G(t)＝G0*Xn(t)；

其中，X(t)为所述新的渐变点；t为所述渐变波段持续时间中的时刻；w为调节因子，且0<|w|<1；Xn(t)表征X(t)的归一化结果；G0为所述渐变点的幅度值；G(t)为所述新的渐变点的幅度值；多个所述新的渐变点组成新的渐变波段，且多个所述新的渐变点的幅度值之和等于所述调整面积值。

需要说明的是，X(t)为新的渐变点，w的取值范围与渐变点位于平台波段相对位置有关，也即多个新的渐变点随时间呈上升趋势，或者多个新的渐变点随时间呈下降趋势，与渐变点位于平台波段相对位置有关。

在本实施例中，当X(t)的t时刻处于平台波段到来之前，调节因子w取值范围在0与1之间，也即0<w<1；当X(t)的t时刻处于平台波段到来之后，调节因子w取值范围在-1与0之间，也即-1<w<0。

如图2所示，在本申请的所有实施例中，待调整信号梯度波形201中的渐变波段(a1、a2)可以为包含临界拐点(P1、P2)对应采样点的波段，且该渐变波段(a1、a2)可以与平台波段相同，也即波段a1与波段a2都与平台波段a的幅度值相同。与待调整信号梯度波形201中的渐变波段不同的是，新的渐变波段(b1、b2)上每个新的渐变点的幅度值都不相同，如图2所示，波段b1和波段b2上所有采样时间点对应的幅度值均不相同，通过对待调整信号梯度波形201进行调整，得到由新的平台波段b与新的渐变波段(b1、b2)组成的目标信号梯度波形202，其面积与待调整信号梯度波形201的目标面积相同。

如图2所示，由于在a1波段中的渐变点均在平台波段到来之前，也即对a1波段的调整应当是随时序逐渐上升的趋势调整，进而得到b1波段，因此当X(t)的t时刻处于在a1波段中，调节因子w取值范围在0与1之间；由于在a2波段中的渐变点均在平台波段到来之后，也即对a2波段的调整应当是随时序逐渐下降的趋势调整，进而得到b2波段，因此当X(t)的t时刻处于在a2波段中，调节因子w取值范围在-1与0之间。

需要说明的是，本实施例在对待调整信号梯度波形中的渐变波段进行调整时，只是将部分平台波段划定为渐变波段，连同待调整信号梯度波形中的拐点进行平滑调整，能够在保证大部分平台波段不变的条件下，实现对渐变波段的平滑处理，消除信号梯度波形中的剧烈拐点，使得新的平台波段与新的渐变波段之间能够平滑过渡，从消除梯度场切换过程中的噪音。

S15：基于所述新的渐变波段与所述新的平台波段组成的目标信号梯度波形，控制所述梯度场。

在步骤S15中，目标信号梯度波形的面积与待调整信号梯度波形的目标面积相等，因此目标信号梯度波形的面积所表征的K空间的位置，与目标面积所表征的K空间的位置为同一位置。

需要说明的是，基于新的渐变波段与新的平台波段组成的目标信号梯度波形，控制梯度场，是利用目标信号梯度波形对应的新的控制参数，控制梯度场进行切换。

在本申请的所有实施例中，由于控制参数是控制梯度场切换的控制参数，对控制参数所描述的待调整信号梯度波形进行调整，实际上是对该控制参数进行调整，因此基于调整后得到新的渐变波段与新的平台波段组成的目标信号梯度波形，对应调整后的控制参数，故基于新的渐变波段与新的平台波段组成的目标信号梯度波形控制梯度场，就是基于目标信号梯度波形对应的调整后的控制参数控制梯度场进行切换操作。

以上可以看出，本实施例提供的一种基于梯度波形调整的磁共振参数成像方法，通过获取预先设定的扫描序列参数与控制参数，由于控制参数能够用于描述待调整信号梯度波形，因此根据扫描序列参数与控制参数，可以确定待调整信号梯度波形的目标面积值，且待调整信号梯度波形包括平台波段与渐变波段，由于该渐变波段所表征的信号拐点落差越大，梯度场切换过程中的噪音就越大，因此按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，对渐变波段进行平滑调整得到新的渐变波段，进而降低了原有渐变波段所表征的信号拐点的落差，最后基于新的渐变波段与新的平台波段组成的目标信号梯度波形控制梯度场，即可实现在磁共振成像过程中对切换梯度场所导致的噪音进行降噪处理，无需对梯度场切换过程中产生的噪音进行采集，再匹配相应的降噪信号进行降噪操作，简化了降噪方案且节约了降噪成本。

此外，按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，调整渐变波段得到新的渐变波段，使得新的渐变波段与新的平台波段组成的目标信号梯度波形，保持有待调整信号梯度波形的目标面积值不变的情况下，尽可能保留平台波段，也即实现最大化地保留控制参数同时，降低梯度切换剧烈带来的噪声问题。

请参阅图3，图3是本申请另一实施例提供的一种基于梯度波形调整的磁共振参数成像方法的实现流程图。相对于图1对应的实施例，本实施例提供的基于梯度波形调整的磁共振参数成像方法在步骤S15之后还包括S21～S22。详述如下：

S21：确定与所述目标信号梯度波形对应的目标梯度函数。

S22：将所述目标梯度函数与所述控制参数关联存储至预设数据库中。

在本实施例中，目标梯度函数是基于目标信号梯度波形所得出的函数。也即，将新的渐变波段与新的平台波段进行重组后，得到目标信号梯度波，基于该目标信号梯度波即可得出目标梯度函数。

需要说明的是，由于待调整信号梯度波形对应的梯度函数G(t^′)中，采样时间点对应有变化剧烈的拐点，因此容易导致梯度场切换时产生较大噪音。而通过对待调整信号梯度波形进行调整后，将平台波形的幅度值增加，将渐变波形平滑调整，最后将新的渐变波段与新的平台波段进行重组后，得到目标信号梯度波，基于该目标信号梯度波形即可得出目标梯度函数，由于该目标梯度函数所表征的梯度波形为目标信号梯度波形，且该波形中没有剧烈的幅度值拐点，因此不会在控制梯度场切换时产生较大的噪声。

应当理解的是，将目标梯度函数与控制参数关联存储至预设数据库中，当再次需要对控制参数对应的待调整信号梯度波形进行调整时，可以直接根据控制参数从预设数据库中获取相应的目标梯度函数，也即调整之后梯度波形所对应的目标梯度函数，即可直接进行基于梯度波形调整的磁共振参数成像。

以上可以看出，本实施例提供的一种基于梯度波形调整的磁共振参数成像方法，通过获取预先设定的扫描序列参数与控制参数，由于控制参数能够用于描述待调整信号梯度波形，因此根据扫描序列参数与控制参数，可以确定待调整信号梯度波形的目标面积值，且待调整信号梯度波形包括平台波段与渐变波段，由于该渐变波段所表征的信号拐点落差越大，梯度场切换过程中的噪音就越大，因此按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，对渐变波段进行平滑调整得到新的渐变波段，进而降低了原有渐变波段所表征的信号拐点的落差，最后基于新的渐变波段与新的平台波段组成的目标信号梯度波形控制梯度场，即可实现在磁共振成像过程中对切换梯度场所导致的噪音进行降噪处理，无需对梯度场切换过程中产生的噪音进行采集，再匹配相应的降噪信号进行降噪操作，简化了降噪方案且节约了降噪成本。

此外，按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，调整渐变波段得到新的渐变波段，使得新的渐变波段与新的平台波段组成的目标信号梯度波形，保持有待调整信号梯度波形的目标面积值不变的情况下，尽可能保留平台波段，也即实现最大化地保留控制参数同时，降低梯度切换剧烈带来的噪声问题。

此外，通过确定与目标信号梯度波形对应的目标梯度函数，再将改目标梯度函数与控制参数关联存储至预设数据库中，使得在下次采用与控制参数一直的信号对梯度场进行控制时，无需重复梯度波形的调整步骤，即可从预设数据库中查找目标信号梯度波形或者与其对应的目标梯度函数进行梯度场切换控制。

请参阅图4，图4是本申请实施例提供的一种基于梯度波形调整的磁共振参数成像装置的结构框图。本实施例中该基于梯度波形调整的磁共振参数成像装置包括的各单元用于执行图1与图3对应的实施例中的各步骤。具体请参阅图1与图3以及图1与图3所对应的实施例中的相关描述。为了便于说明，仅示出了与本实施例相关的部分。参见图4，基于梯度波形调整的磁共振参数成像装置400包括：第一获取单元41、第一确定单元42、第一调整单元43、第二调整单元44以及执行单元45。其中：

第一获取单元41，用于获取预先设定的扫描序列参数与控制参数；其中，所述控制参数用于描述待调整信号梯度波形。

第一确定单元42，用于根据所述扫描序列参数与所述控制参数，确定所述待调整信号梯度波形的目标面积值；其中，所述待调整信号梯度波形包括平台波段与渐变波段。

第一调整单元43，用于按照预设幅度调整参数调整所述平台波段的波形幅度，得到新的平台波段。

第二调整单元44，用于基于所述目标面积值与所述新的平台波段的第一面积值，对所述渐变波段进行平滑调整，得到新的渐变波段；其中，所述新的渐变波段的第二面积值与所述第一面积值之和等于所述目标面积值。

执行单元45，用于基于所述新的渐变波段与所述新的平台波段组成的目标信号梯度波形，控制所述梯度场。

作为本申请一实施例，所述控制参数包括用于描述所述待调整信号梯度波形的梯度函数；所述扫描序列参数包括：K空间尺寸、K空间单元尺寸、原子核的旋磁比、扫描视野、带宽、采样时间，以及与所述采样时间关联的采样点个数。

作为本申请一实施例，第一确定单元42具体用于，通过以下公式测算出所述待调整信号梯度波形的目标面积值；

k＝N·Δk；

其中，k(t)为所述采样时间为t时刻的K空间位置；γ为所述原子核的旋磁比；G(t′)为所述梯度函数；k为所述K空间尺寸；N为所述采样点个数；Δk为所述K空间单元尺寸；FOV为所述扫描视野；A为所述目标面积值；BW为所述带宽。

作为本申请一实施例，第一调整单元43具体用于，根据所述梯度函数确定所述平台波段的平台幅度值；按照预设幅度调整参数调整所述平台幅度值，得到新的平台幅度值；其中，所述新的平台波段的幅度值等于所述平台幅度值与所述调整参数之和；根据所述新的平台幅度值得到所述新的平台波段。

作为本申请一实施例，所述采样时间包括平台波段持续时间与渐变波段持续时间；所述采样点个数包括与所述平台波段持续时间对应的第一采样点个数，以及与所述渐变波段持续时间对应的第二采样点个数。

第二调整单元44具体用于，获取所述第一采样点个数；将所述第一采样点个数与所述新的平台幅度值的乘积，识别为所述第一面积值；测算所述目标面积值与所述第一面积值之差，得到调整面积值；基于所述调整面积值对所述渐变波段进行平滑调整，得到新的渐变波段。

作为本申请一实施例，所述渐变波段包括在所述渐变波段持续时间内连续的多个渐变点；第二调整单元44具体用于，

获取所述渐变波段持续时间；

基于所述渐变波段持续时间确定多个渐变点；

通过以下公式调整每个所述渐变点的幅度值，得到多个新的渐变点；

X(t)＝1-exp(-w*t)；

G(t)＝G0*Xn(t)；

其中，X(t)为所述新的渐变点；t为所述渐变波段持续时间中的时刻；w为调节因子，且0<|w|<1；Xn(t)表征X(t)的归一化结果；G0为所述渐变点的幅度值；G(t)为所述新的渐变点的幅度值；多个所述新的渐变点组成新的渐变波段，且多个所述新的渐变点的幅度值之和等于所述调整面积值。

作为本申请一实施例，基于梯度波形调整的磁共振参数成像装置400，还包括：第二确定单元46与存储单元47。具体地：

第二确定单元46，用于确定与所述目标信号梯度波形对应的目标梯度函数。

存储单元47，用于将所述目标梯度函数与所述控制参数关联存储至预设数据库中。

以上可以看出，本实施例提供的方案，通过获取预先设定的扫描序列参数与控制参数，由于控制参数能够用于描述待调整信号梯度波形，因此根据扫描序列参数与控制参数，可以确定待调整信号梯度波形的目标面积值，且待调整信号梯度波形包括平台波段与渐变波段，由于该渐变波段所表征的信号拐点落差越大，梯度场切换过程中的噪音就越大，因此按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，对渐变波段进行平滑调整得到新的渐变波段，进而降低了原有渐变波段所表征的信号拐点的落差，最后基于新的渐变波段与新的平台波段组成的目标信号梯度波形控制梯度场，即可实现在磁共振成像过程中对切换梯度场所导致的噪音进行降噪处理，无需对梯度场切换过程中产生的噪音进行采集，再匹配相应的降噪信号进行降噪操作，简化了降噪方案且节约了降噪成本。

此外，按照预设幅度调整参数调整待调整信号梯度波形的平台波段的波形幅度，再基于目标面积值与新的平台波段的第一面积值，调整渐变波段得到新的渐变波段，使得新的渐变波段与新的平台波段组成的目标信号梯度波形，保持有待调整信号梯度波形的目标面积值不变的情况下，尽可能保留平台波段，也即实现最大化地保留控制参数同时，降低梯度切换剧烈带来的噪声问题。

此外，通过确定与目标信号梯度波形对应的目标梯度函数，再将改目标梯度函数与控制参数关联存储至预设数据库中，使得在下次采用与控制参数一直的信号对梯度场进行控制时，无需重复梯度波形的调整步骤，即可从预设数据库中查找目标信号梯度波形或者与其对应的目标梯度函数进行梯度场切换控制。

图5是本申请另一实施例提供的一种磁共振成像设备的结构框图。如图5所示，该实施例的磁共振成像设备5包括：处理器50、存储器51以及存储在所述存储器51中并可在所述处理器50上运行的计算机程序52，例如基于梯度波形调整的磁共振参数成像方法的程序。处理器50执行所述计算机程序52时实现上述各个基于梯度波形调整的磁共振参数成像方法各实施例中的步骤，例如图1所示的S11至S15。或者，所述处理器50执行所述计算机程序52时实现上述图4对应的实施例中各单元的功能，例如，图4所示的单元41至45的功能，具体请参阅图5对应的实施例中的相关描述，此处不赘述。

示例性的，所述计算机程序52可以被分割成一个或多个单元，所述一个或者多个单元被存储在所述存储器51中，并由所述处理器50执行，以完成本申请。所述一个或多个单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序52在所述磁共振成像设备5中的执行过程。例如，所述计算机程序52可以被分割成第一获取单元、第一确定单元、第一调整单元、第二调整单元以及执行单元，各单元具体功能如上所述。

所述磁共振设备可包括，但不仅限于，处理器50、存储器51。本领域技术人员可以理解，图5仅仅是磁共振成像设备5的示例，并不构成对磁共振设备5的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述磁共振设备还可以包括输入输出设备、网络接入设备、总线等。

所称处理器50可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

所述存储器51可以是所述磁共振成像设备5的内部存储单元，例如磁共振设备5的硬盘或内存。所述存储器51也可以是所述磁共振成像设备5的外部存储设备，例如所述磁共振成像设备5上配备的插接式硬盘，智能存储卡(Smart Media Card，SMC)，安全数字(Secure Digital，SD)卡，闪存卡(Flash Card)等。进一步地，所述存储器51还可以既包括所述磁共振成像设备5的内部存储单元也包括外部存储设备。所述存储器51用于存储所述计算机程序以及所述磁共振设备所需的其他程序和数据。所述存储器51还可以用于暂时地存储已经输出或者将要输出的数据。

以上所述实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。

Claims (10)

1.一种基于梯度波形调整的磁共振参数成像方法，其特征在于，包括：

获取预先设定的扫描序列参数与控制参数；其中，所述控制参数用于描述待调整信号梯度波形；

根据所述扫描序列参数与所述控制参数，确定所述待调整信号梯度波形的目标面积值；其中，所述待调整信号梯度波形包括平台波段与渐变波段；

按照预设幅度调整参数调整所述平台波段的波形幅度，得到新的平台波段，具体为：在平台波段原有的波形幅度的基础上增加其波形幅度值的大小；

基于所述目标面积值与所述新的平台波段的第一面积值，对所述渐变波段进行平滑调整，得到新的渐变波段；其中，所述新的渐变波段的第二面积值与所述第一面积值之和等于所述目标面积值；

基于所述新的渐变波段与所述新的平台波段组成的目标信号梯度波形，控制所述梯度场，所述梯度波形调整前后的时间周期不变。

2.根据权利要求1所述的基于梯度波形调整的磁共振参数成像方法，其特征在于，所述控制参数包括用于描述所述待调整信号梯度波形的梯度函数；

所述扫描序列参数包括：K空间尺寸、K空间单元尺寸、原子核的旋磁比、扫描视野、带宽、采样时间，以及与所述采样时间关联的采样点个数。

3.根据权利要求2所述的基于梯度波形调整的磁共振参数成像方法，其特征在于，所述根据所述扫描序列参数与所述控制参数，确定所述待调整信号梯度波形的目标面积值，包括：

通过以下公式测算出所述待调整信号梯度波形的目标面积值；

k＝N·Δk；

其中，k(t)为所述采样时间为t时刻的K空间位置；γ为所述原子核的旋磁比；G(t′)为所述梯度函数；k为所述K空间尺寸；N为所述采样点个数；Δk为所述K空间单元尺寸；FOV为所述扫描视野；A为所述目标面积值；BW为所述带宽。

4.根据权利要求2所述的基于梯度波形调整的磁共振参数成像方法，其特征在于，所述按照预设幅度调整参数调整所述平台波段的波形幅度，得到新的平台波段，包括：

根据所述梯度函数确定所述平台波段的平台幅度值；

按照预设幅度调整参数调整所述平台幅度值，得到新的平台幅度值；其中，所述新的平台波段的幅度值等于所述平台幅度值与所述调整参数之和；

根据所述新的平台幅度值得到所述新的平台波段。

5.根据权利要求4所述的基于梯度波形调整的磁共振参数成像方法，其特征在于，所述采样时间包括平台波段持续时间与渐变波段持续时间；所述采样点个数包括与所述平台波段持续时间对应的第一采样点个数，以及与所述渐变波段持续时间对应的第二采样点个数；

所述基于所述目标面积值与所述新的平台波段的第一面积值，对所述渐变波段进行平滑调整，得到新的渐变波段，包括：

获取所述第一采样点个数；

将所述第一采样点个数与所述新的平台幅度值的乘积，识别为所述第一面积值；

测算所述目标面积值与所述第一面积值之差，得到调整面积值；

基于所述调整面积值对所述渐变波段进行平滑调整，得到新的渐变波段。

6.根据权利要求5所述的基于梯度波形调整的磁共振参数成像方法，其特征在于，所述渐变波段包括在所述渐变波段持续时间内连续的多个渐变点；

所述基于所述调整面积值对所述渐变波段进行平滑调整，得到新的渐变波段，包括：

获取所述渐变波段持续时间；

基于所述渐变波段持续时间确定多个渐变点；

通过以下公式调整每个所述渐变点的幅度值，得到多个新的渐变点；

X(t)＝1-exp(-w*t)；

G(t)＝G0*Xn(t)；

其中，X(t)为所述新的渐变点；t为所述渐变波段持续时间中的时刻；w为调节因子，且0<|w|<1；Xn(t)表征X(t)的归一化结果；G0为所述渐变点的幅度值；G(t)为所述新的渐变点的幅度值；多个所述新的渐变点组成新的渐变波段，且多个所述新的渐变点的幅度值之和等于所述调整面积值。

7.根据权利要求1至6任一项所述的基于梯度波形调整的磁共振参数成像方法，其特征在于，所述方法还包括：

确定与所述目标信号梯度波形对应的目标梯度函数；

将所述目标梯度函数与所述控制参数关联存储至预设数据库中。

8.一种基于梯度波形调整的磁共振参数成像装置，其特征在于，包括：

第一获取单元，用于获取预先设定的扫描序列参数与控制参数；其中，所述控制参数用于描述待调整信号梯度波形；

第一确定单元，用于根据所述扫描序列参数与所述控制参数，确定所述待调整信号梯度波形的目标面积值；其中，所述待调整信号梯度波形包括平台波段与渐变波段；

第一调整单元，用于按照预设幅度调整参数调整所述平台波段的波形幅度，得到新的平台波段，具体为：在平台波段原有的波形幅度的基础上增加其波形幅度值的大小；

第二调整单元，用于基于所述目标面积值与所述新的平台波段的第一面积值，对所述渐变波段进行平滑调整，得到新的渐变波段；其中，所述新的渐变波段的第二面积值与所述第一面积值之和等于所述目标面积值；

执行单元，用于基于所述新的渐变波段与所述新的平台波段组成的目标信号梯度波形，控制所述梯度场，所述梯度波形调整前后的时间周期不变。

9.一种磁共振成像设备，其特征在于，所述磁共振设备包括存储器、处理器以及存储在所述存储器中并可在所述磁共振设备上运行的计算机程序，所述处理器执行所述计算机程序时实现如权利要求1至7任一项所述基于梯度波形调整的磁共振参数成像方法的步骤。

10.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述基于梯度波形调整的磁共振参数成像方法的步骤。

Images