CN106918794B - 磁共振系统及成像的方法 - Google Patents

磁共振系统及成像的方法 Download PDF

Info

Publication number
CN106918794B
CN106918794B CN201510992357.7A CN201510992357A CN106918794B CN 106918794 B CN106918794 B CN 106918794B CN 201510992357 A CN201510992357 A CN 201510992357A CN 106918794 B CN106918794 B CN 106918794B
Authority
CN
China
Prior art keywords
field
radio frequency
echo signal
target
magnetic resonance
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.)
Active
Application number
CN201510992357.7A
Other languages
English (en)
Other versions
CN106918794A (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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201510992357.7A priority Critical patent/CN106918794B/zh
Publication of CN106918794A publication Critical patent/CN106918794A/zh
Application granted granted Critical
Publication of CN106918794B publication Critical patent/CN106918794B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/387Compensation of inhomogeneities

Abstract

本发明公开了一种磁共振系统成像的方法,具体包括:在多通道发射模式下,根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化以产生目标B1场,所述目标B1场同时满足B1场均匀性和射频发射功率优化;将经过优化的幅度和相位参数加载到射频脉冲控制器,驱动多个射频发射通道产生目标B1场;在所述目标B1场下激发受检者待扫描部位,获取待扫描部位的磁共振信号;对所述磁共振信号进行傅里叶变换,获取待扫描部位的磁共振图像。本发明通过采集自旋回波信号和刺激回波信号强度反映B1场的分布,可快速校正得到均匀B1场并可对射频发射功率校准。此外,本发明还提出一种磁共振系统。

Description

磁共振系统及成像的方法
【技术领域】
本发明涉及磁共振成像技术领域,尤其涉及射频场优化相关的技术。
【背景技术】
磁共振成像(Magnetic Resonance Imaging,MRI)技术是一种随着计算机技术、电子电路技术、超导体技术的发展而产生的一种非介入式探测技术,在当今的临床诊断和医学研究中有着极其重要的意义。磁共振成像的基本过程是:发射单元产生垂直于扫描设备主磁场的射频场,射频脉冲能量使人体内共振的原子核自旋,此时的人体相当于接收机,被准确地调谐在共振频率上,射频脉冲停止后,人体被照射部分的质子经过弛豫过程逐渐回到原来的平衡位置,同时发射出与刺激脉冲波频率相同的射频信号,这一过程称为自由感应衰减振荡。病灶内的氢原子具有不同于正常组织内的氢原子的“弛豫”时间常数,这样可以获得人体内各处质子密度和质子周围的化学环境的图像,医生可通过图像准确诊断组织的早期病变。自磁共振成像技术应用以来,为满足人们对成像高分辨率、高信噪比的需求,扫描设备的主磁场场强不断提升,对应的射频脉冲载波频率也随之提高。然而,上述操作的结果是短波长的射频场与负载的耦合效应增强,导致射频能量在空间上的不均匀传递以及对空间内质子的不均匀激发,成像结果不够准确反映真实信息,从而无法提供用于临床诊断和研究的可靠参照。因此,射频(B1)场的均匀性是决定磁共振性能的重要指标。
在扫描区内均匀的射频线圈发射磁场,才能使线圈在共振频率处获得高增益Q值,从而获得清晰的图像。B1场的不均匀可以造成不同空间位置处的翻转角差异,而无论是自旋回波或者梯度回波序列,实际的翻转角常常低于理论数值,通常会造成特定范围内的信号减小;或者在损毁梯度回波序列中,由于序列的T1对比主要取决于翻转角,可能造成图像对比特征的改变。此外,180°或90°饱和脉冲的翻转角差异,也会造成所抑制组织的不完全。因此,高场MRI中B1场的均匀性是一个关键性技术问题,直接影响到图像的均匀性以及对比度。现有技术中,B1场均匀性的获取方法主要包括:扩大线圈的几何尺寸和多通道并行发射技术。其中,扩大线圈的几何尺寸在提高B1场均匀性的同时,会加大线圈损耗、减小填充因数,而且由于线圈内所包含的组织量增多,还会形成额外的噪声来源,虽然B1场均匀性增强却伴随信噪比(Signal Noise Ratio,S/N)降低;多通道并行发射技术主要根据所有单个发射通道的B1场调节各个通道的幅值以及相位优化B1场的均匀性,然而其调节过程比较耗时,难以满足临床序列扫描前进行预扫描系统校正的要求。此外,MRI系统中为了使发射的射频脉冲达到需要的翻转角,还需要再在扫描前对体线圈的射频发射功率进行校准,高场射频激发场的不均匀使得传统的校准方法得出的电压值只能代表激发区域的平均水平,无法对高场下的B1场的均匀性问题提供帮助。鉴于此,有必要对现有B1场校准方法进行改进。
【发明内容】
本发明所要解决的技术问题是提供一种可快速、准确实现B1场均匀性优化的方法。
本发明解决上述技术问题所采用的技术方案为一种磁共振系统成像的方法,所述磁共振系统包括若干个射频发射通道,具体包括以下步骤:
在多通道发射模式下,根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化以产生目标B1场,所述目标B1场同时满足B1场均匀性和射频发射功率优化;
将经过优化的幅度和相位参数加载到射频脉冲控制器,驱动多个射频发射通道产生目标B1场;
在所述目标B1场下激发受检者待扫描部位,获取待扫描部位的磁共振信号;
对所述磁共振信号进行傅里叶变换,获取待扫描部位的磁共振图像。
进一步地,所述根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化的具体过程为:
(a)在磁共振系统空载状态下,将初始校正得到的若干个射频发射通道的幅度、相位参数和射频发射功率初始值加载到射频脉冲控制器;
(b)驱动所述多个射频发射通道并采用刺激回波序列获取成像区域内自旋回波信号强度SSE和刺激回波信号强度SSTE
(c)根据自旋回波信号强度SSE和刺激回波信号强度SSTE判断当前B1场是否为目标B1场,如果否,则调整所述若干个射频发射通道的幅度、相位参数,并返回步骤(b);否则,结束优化过程。
进一步地,所述根据自旋回波号强度SSE和刺激回波信号强度SSTE判断当前B1场是否为目标B1场的具体过程为:
根据自旋回波信号强度SSE和刺激回波信号强度SSTE获取当前翻转角θ;
获取目标B1场所对应的目标翻转角θm,计算所述当前翻转角θ与所述目标翻转角θm的差值是否在设定范围内,如果是,则判定当前B1场为目标B1场;如果否,则判定当前B1场不是目标B1场。
进一步地,所述目标翻转角θm=90°。
进一步地,所述刺激回波含有三个射频脉冲,对应翻转角分别为θ1、θ2、θ3,且θ1=θ2=θ3=θ,0≤θ≤180°,则
Figure BDA0000889996690000041
其中θ表示当前翻转角。
进一步地,所述刺激回波含有三个射频脉冲,对应翻转角分别为θ1、θ2、θ3,且θ1=θ3=θ,θ2=2θ,0≤θ≤90°,则
Figure BDA0000889996690000042
其中θ表示当前翻转角。
进一步地,所述根据自旋回波信号强度SSE和刺激回波信号强度SSTE判断当前B1场是否为目标B1场的具体过程为:
获取当前B1场下采集的刺激回波信号强度SSTE与自旋回波信号强度SSE的比值i,
Figure BDA0000889996690000043
获取目标B1场对应的刺激回波信号强度SSTE与自旋回波信号强度SSE的参考比值j,0≤j≤2;
如果i与j的差值在设定范围内,则判定当前B1场为目标B1场;否则,则判定当前B1场不是目标B1场。
进一步地,还包括对目标B1场对应的射频发射功率PTar进行归一化,获取参考射频发射功率PRef,Pref=K×PTar,其中比例因子
Figure BDA0000889996690000044
Figure BDA0000889996690000045
其中B1Tar表示目标B1场的幅度值,B1Ref表示参考B1场的幅度值,θm表示目标B1场对应的目标翻转角,θref表示参考B1场对应的参考翻转角。
本发明还提供一种磁共振系统,包括:
体发射线圈模块,包含多个射频线圈发射通道;
优化模块,用于根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化;
驱动模块,用于根据多个射频发射通道优化后的幅度和相位参数驱动所述多个射频发射通道产生目标B1场,所述目标B1场同时满足B1场均匀性和射频发射功率优化;
信号采集模块,用于在所述目标B1场下采集待扫描部位的磁共振信号;
重建模块,用于对所述磁共振信号进行傅里叶变换,获取待扫描区域的磁共振图像。
进一步地,所述优化模块根据待成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE获取当前B1场对应的当前翻转角θ,根据当前翻转角θ与目标B1场对应的目标翻转角θm的差值确定是否对多个射频发射通道的幅度和相位参数进行优化;
或所述优化模块获取当前B1场下采集的刺激回波信号强度SSTE与自旋回波信号强度SSE的比值i,根据i和目标B1场对应的刺激回波信号强度与自旋回波信号强度的参考比值j,确定是否对射频发射通道的幅度和相位参数进行优化,其中,0≤i≤2,0≤j≤2。
与现有技术相比,本发明的优点在于:通过采用刺激回波获取成像区域内的自旋回波和刺激回波信号,进而获得翻转角,而翻转角与B1场存在线性对应关系,因此通过采集自旋回波信号和刺激回波信号的强度可反映B1场的情况;调节多个射频通道的幅度和相位参数可获得与之对应的自旋回波信号和刺激回波信号,通过检测上述两种信号可快速校正得到均匀B1场,提高图像的均匀性和对比度;采用迭代优化方法实现B1场均匀性校准的同时也可对射频发射功率校准,保证了校准值的一致性和准确性,在一定程度上避免局部高场射频能量吸收率过高。
【附图说明】
图1为本发明磁共振系统成像的方法流程图;
图2为射频发射通道的幅度和相位参数优化流程图;
图3为在成像区域内采集自旋回波和刺激回波信号示意图;
图4a为射频发射通道幅相参数优化前得到的B1场分布图;
图4b为采用本发明方法优化的射频发射通道幅相参数得到的B1场分布图;
图5为射频线圈发射通道结构示意图。
【具体实施方式】
为使本发明的上述目的、特征和优点能够更为明显易懂,下面结合附图和实施例对本发明的具体实施方式做详细的说明。
磁共振系统主要包括扫描控制系统、梯度系统、射频系统和超导磁体,其中,超导磁体、梯度系统、射频系统设置在扫描腔的周围,超导磁体在超导开关控制下产生均匀稳定的静磁场B0的主磁场,主磁场的强度均匀性、扫描腔的大小等都会影响磁共振信号的信噪比和图像分辨率;梯度系统由梯度线圈、梯度控制器、数模转换器、梯度放大器和梯度冷却系统等组成,提供系统线性度满足要求的、可快速开关的梯度场,根据需要动态地在主磁场附加一个X、Y、Z正交的三维空间线性变化的梯度磁场,一个方向的梯度用于射频脉冲选择性的激发一个层面内质子的自旋,第二个梯度对沿层面内一个方向的MR信号进行频率编码,第三个梯度对沿层面内另一个方向的MR信号进行相位空间编码,使受检者在不同位置具有不同的共振频率,实现成像体素的选层和空间编码的功能,上述三个梯度场分别由三个正交的梯度交流线圈产生,每一组线圈要求有一个独立的电流驱动放大器供电,每组线圈由两个电流方向相反的同轴线圈组成,以产生其轴向的最大梯度磁场;射频系统由射频线圈、发射接收系统、射频功率系统等组成,是磁共振系统中实施射频激励并接收和处理射频信号的功能单元,其不仅要根据扫描序列的要求发射各种翻转角的射频波,还要接收成像区域内氢质子的共振信号,射频系统中的射频线圈主要在射频控制器的作用下,提供扫描序列所需要的各种射频脉冲,为了保证图像空间的均匀性,产生的B1场在整个成像区域内要尽可能均匀。
通常情况下,射频激发后,线圈中感应出的信号正比于磁化矢量的横向分量,横向分量与翻转角直接相关,而翻转角又与B1场成比例,因此,通过计算成像区域内翻转角的大小可获得当前B1场的分布情况。为实现B1场均匀性,本发明一种磁共振系统成像的方法,其中磁共振系统包括若干个射频发射通道,具体包括以下步骤:
S10、在多通道发射模式下,根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化以产生目标B1场,该目标B1场既可满足B1场均匀性也可对射频发射功率优化;
S20、将经过优化的幅度和相位参数加载到射频脉冲控制器中,驱动多个射频发射通道,其中射频通道与射频发射线圈一一对应,从而产生目标B1场;
S30、在目标B1场下激发受检者待扫描部位,获取待扫描部位的磁共振信号;
S40、对获取的磁共振信号进行傅里叶变换,获取待扫描部位的磁共振图像。
在本发明方法中,根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化具体过程如下:
(a)在磁共振系统空载状态下,将初始校正得到的若干个射频发射通道的幅度、相位参数和射频发射功率初始值PIni加载到射频脉冲控制器。磁共振系统包括若干个射频发射通道(通道数大于等于2),每个射频发射通道都可单独发射射频场,多个射频发射通道形成合成矢量B1场,合成场与各通道发射场之间满足矢量合成关系,即各通道独立B1场通过矢量叠加形成了合成的系统B1场。在此具体实施例中,以双通道发射为例,记射频通道A产生的独立B1场为
Figure BDA0000889996690000085
射频通道B产生的独立B1场为
Figure BDA0000889996690000086
则B1场与单通道
Figure BDA0000889996690000087
之间的合成关系满足:
Figure BDA0000889996690000088
其中αA、αB为复数,且
Figure BDA0000889996690000081
Figure BDA0000889996690000082
a0、a1分别为通道A和通道B产生的射频脉冲的幅值比例因子,βA,βB分别表示通道A和通道B产生的射频脉冲的相位因子。
(b)驱动多个射频发射通道并采用刺激回波序列获取激发区域内自旋回波信号强度SSE和刺激回波信号强度SSTE。本具体实施例中,根据步骤(a)中每个通道的幅度和相位参数,如图3所示,Gss为选层梯度,在选层梯度时间范围内,采用包含三个射频脉冲的刺激回波(STEAM)序列,三个射频脉冲对应的翻转角分别为θ1、θ2、θ3,且满足θ1=θ2=θ3=θ,θ表示当前翻转角,0≤θ≤180°,,τ1和τ2是射频脉冲到采集窗的时间间隔,通过数据采集窗口ADC获取的成像区域内自旋回波和刺激回波信号的强度是成像区域内每一点信号值的总和,且成像区域内任一点的坐标用(x,y)表示,待成像区域内自旋回波以及刺激回波信号强度可表示如下:
Figure BDA0000889996690000083
Figure BDA0000889996690000084
其中,ρ0表示纵向磁化,T1为氢质子纵向弛豫时间,T2为氢质子横向弛豫时间。
(c)根据自旋回波信号强度SSE和刺激回波信号强度SSTE判断当前B1场是否为目标 B1场,如果否,则调整若干个射频发射通道的幅度、相位参数,并返回步骤(b);否则,结束优 化过程。翻转角与B1场的分布满足线性关系式,以双射频通道为例翻转角其中,Tp表示脉冲持续的时间,A(t)为t时刻施加的脉冲波 形,γ为旋磁比,由上述公式可知在脉冲波形恒定的情况下,翻转角与B1成线性关系,因此, 翻转角与B1场存在对应关系,通过翻转角可间接反映B1场的分布情况。本实施例中获取自旋 回波和刺激回波信号强度的比值,并忽略弛豫项,可得
Figure BDA0000889996690000092
因此,可 计算得到当前翻转角
利用B1场与翻转角的线性对应关系,获取目标B1场所对应的目标翻转角为θm,0≤θm≤180°,计算当前翻转角θ与目标翻转角θm的差值是否在设定范围内,具体为:根据优化需求选定收敛域值ε11为设定的阈值),计算当前翻转角θ与目标翻转角θm的差值是否小于ε1,如果满足该条件,则判定此时的翻转角θ达到目标翻转角θm;否则,判定当前B1场不是目标B1场,此时需调整射频发射通道A、B的幅度、相位参数(调整射频通道的幅度和相位参数可改变翻转角),并返回步骤(b),直至判定当前翻转角θ与目标翻转角θm的差值位于设定范围内(最理想的情况为调整为当前翻转角θ与目标翻转角相等),并结束优化过程。
如图4a所示为射频发射通道幅相参数优化前得到的B1场分布图,如图4b所示为采用本发明方法对射频发射通道幅相参数优化后得到的B1场分布图,两图都以翻转角表示。在未优化之前,由于受双通道相位相互抵消的影响,B1场分布不均匀,具体表现为部分边缘区域信号强,另一部分边缘区域信号较弱,整体分布不均匀;而采用优化后的参数得到的B1场均匀性明显改善,可避免出现边缘区域相位相互抵消的情况。
需要说明的是,刺激回波选择序列的不同,通过计算得到的当前翻转角公式也会不同。在实施例二中,刺激回波含有三个射频脉冲,对应翻转角分别为θ1、θ2、θ3,且θ1=θ3=θ,θ2=2θ,0≤θ≤90°,θ表示当前射频发射功率翻转角。如实施例一中所述,自旋回波以及刺激回波信号强度可如下表示:
Figure BDA0000889996690000101
Figure BDA0000889996690000102
在忽略弛豫项的情况下,自旋回波和刺激回波信号强度的比值为:
Figure BDA0000889996690000103
根据自旋回波信号强度SSE和刺激回波信号强度SSTE获得的当前翻转角θ表达式为:
Figure BDA0000889996690000104
目标B1场所对应的目标翻转角为θm,0≤θm≤90°。计算当前翻转角θ与目标翻转角θm的差值是否在设定范围内,具体为:根据优化需求选定收敛域值ε2(设定的阈值),计算当前翻转角θ与目标翻转角θm的差值是否小于ε2,如果满足该条件,则判定此时的翻转角θ达到目标翻转角θm;否则,判定当前B1场不是目标B1场,此时需调整射频发射通道A、B的幅度、相位参数(调整射频通道的幅度和相位参数可改变翻转角),并返回步骤(b),直至判定当前翻转角θ与目标翻转角θm的差值位于设定范围内(最理想的情况为调整为当前翻转角θ与目标翻转角θm相等),并结束优化迭代过程。
在上述两个实施例中,调整各个射频发射通道的幅度、相位参数可改变B1场的分布,B1场分布的变化通过计算当前发射功率对应的翻转角θ体现,而计算检测得到的自旋回波信号强度SSE和刺激回波信号强度SSTE可得到当前发射功率达到的翻转角θ,因此,也可直接通过计算自旋回波信号强度SSE和刺激回波信号强度SSTE的比值判断当前射频场是否达到目标场。在实施例三中,根据初始校正分配得到的每个通道的幅度和相位参数,采用包含三个射频脉冲的刺激回波(STEAM)序列,三个射频脉冲对应的翻转角分别为θ1、θ2、θ3,且满足θ1=θ2=θ3=θ,θ表示当前翻转角,0≤θ≤180°,如实施例一所述,自旋回波信号强度SSE和刺激回波信号强度SSTE与当前射频发射功率翻转角θ的关系满足:
Figure BDA0000889996690000111
若当前射频B1场调整为目标场,即当前翻转角θ=θm,0≤θm≤180°,则此时应满足关系
Figure BDA0000889996690000112
因此,根据自旋回波信号强度SSE和刺激回波信号强度SSTE判断当前B1场是否为目标B1场的也可按照如下过程实现过程:
获取当前B1场下采集的刺激回波信号强度SSTE与自旋回波信号强度SSE的比值i,
Figure BDA0000889996690000113
获取目标B1场对应的刺激回波信号强度与自旋回波信号强度的参考比值j,j根据具体情况设定,本实施例中
Figure BDA0000889996690000114
如果i与j的差值在设定范围ε3内,则判定当前B1场为目标B1场;否则,则判定当前B1场不是目标B1场,本实施例中设定目标函数Cost=|i-j|,判定目标函数是否为ε3(达到设定范围),如果目标函数为ε3(处于设定范围)结束优化过程;如果目标函数不为ε3(处于设定范围),则调整各个射频通道的幅度和相位参数,根据调整后的参数重新获取一组自旋回波信号强度SSE和刺激回波信号强度SSTE,并重新判断目标函数的值,直至判定目标函数为ε3(处于设定范围),获取此时射频发射功率的翻转角,并根据翻转角与射频发射功率的线性关系,获取目标发射场中各个射频通道经过调整后的最优幅值和相位参数。
在实施例四中,根据初始校正分配得到的每个通道的幅度和相位参数,含有三个射频脉冲的刺激回波序列,三个射频脉冲对应翻转角分别为θ1、θ2、θ3,且θ1=θ3=θ,θ2=2θ,0≤θ≤90°,θ表示当前射频发射功率翻转角,如实施例二所述,自旋回波信号强度SSE和刺激回波信号强度SSTE与当前射频发射功率翻转角θ的关系满足:
Figure BDA0000889996690000121
若当前射频B1场调整为目标场,当前射频发射功率翻转角θ=θm,0≤θm≤90°,则此时应满足关系
Figure BDA0000889996690000122
直接判断刺激回波信号强度SSTE与自旋回波信号强度SSE两者的比值即可判定当前射频场是否调整到目标场,具体为:获取当前B1场下采集的刺激回波信号强度SSTE与自旋回波信号强度SSE的比值i,
Figure BDA0000889996690000123
获取目标B1场对应的刺激回波信号强度与自旋回波信号强度的参考比值j,j根据具体情况设定,本实施例中j=cosθm,0≤j≤1;
如果i与j的差值在设定范围ε4内,则判定当前B1场为目标B1场;否则,则判定当前B1场不是目标B1场,本实施例汇总设定目标函数Cost=|i-j|,判定目标函数是否为ε4(或达到设定范围),如果目标函数为ε4(或处于设定范围)结束优化过程;如果目标函数不为ε4(或处于设定范围),则调整各个射频通道的幅度和相位参数,根据调整后的参数重新获取一组自旋回波信号强度SSE和刺激回波信号强度SSTE,并重新判断目标函数的值,直至判定目标函数为ε4(或处于设定范围),获取此时射频发射功率的翻转角,并根据翻转角与射频发射功率的线性关系,获取目标发射场中各个射频通道经过调整后的最优幅值和相位参数。
本发明的磁共振成像方法,在通过调整各通道相位、幅度参数获得目标B1场后,还包括,对目标B1场对应的射频发射功率PTar进行归一化获取参考射频发射功率PRef,Pref=K×PTar,其中比例因子
Figure BDA0000889996690000131
Figure BDA0000889996690000132
B1Tar表示目标B1场的幅度值,B1Ref表示参考B1场的幅度值,θm目标B1场对应的目标翻转角,θref是参考B1场对应的参考翻转角。在上述成像方法基础上,本发明还提供一种磁共振系统,包括磁体、梯度线圈、射频线圈和控制系统,其中,磁体在其工作孔径内产生匀强磁场,即主磁场或静磁场B0(static magneticfield);磁体的空间内还设置有梯度线圈、匀场线圈、射频线圈从而形成容纳受检者的有效孔径,梯度线圈为磁共振系统提供线性度要求的、可快速开关的梯度场;射频线圈,实施射频激励并接收和处理射频信号。
本发明的磁共振系统中的体发射线圈模块,包含多个射频线圈发射通道,如图5所示,包含有射频振荡器、频率合成器、滤波放大器、波形调制器、脉冲功率放大器、终端匹配网络及RF发射线圈等,作用是在射频脉冲控制器的作用下,控制序列发生器、脉冲发生器通过射频线圈发射通道提供扫描序列所需的各种射频脉冲,其产生的射频B1场垂直于主磁场B0,使得RF脉冲能够将其能量耦合给共振的原子核而引起质子进动。需要指出的是,不同的射频发射功率以及不同的射频发射通道加权矢量合成会形成不同的射频B1场,从而导致空间翻转角分布不同,序列的实现效果也不一样。
在扫描成像前,为保证获取B1场的均匀性,本发明的磁共振系统设置有优化模块,用于根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化,具体为:根据待成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE获取当前B1场对应的当前翻转角,根据当前翻转角与目标B1场对应的目标翻转角的差值确定是否对多个射频发射通道的幅度和相位参数进行优化;或优化模块获取当前B1场下采集的刺激回波信号强度SSTE与自旋回波信号强度SSE的比值i,0≤i≤2,根据i和目标B1场对应的刺激回波信号强度与自旋回波信号强度的参考比值j,0≤j≤2,确定是否对射频发射通道的幅度和相位参数进行优化或调整;驱动模块,即射频脉冲控制器,用于根据多个射频发射通道优化后的幅度和相位参数驱动射频发射通道产生目标B1场,目标B1场同时满足B1场均匀性和射频发射功率均匀性。
射频接收单元接收人体产生的磁共振信号,并经适当放大后供信号采集模块使用。从接收线圈中感应出的信号只有微瓦(μW)数量级的功率,这就要求它既要有很高的放大倍数,又要有很小的噪声,因此射频接收单元中还包括前置放大器,前置放大器可对1μV以下的信号发生反应,同时在工作频率附近要求有较为平坦的频率响应,并在很大范围内有足够的线性放大特性。此外,从安全性能考虑,该前置放大器设置为可接受1V左右的过载,且过载后可在小于1μs的时间内迅速恢复。信号经前置放大器放大后分别经过混频器、低通滤波器和A/D转换器后,可滤除其中混杂的交流成分,并可使模拟信号成为离散数字信号,这些数字信号被信号采集模块接收,获得待扫描部位的磁共振信号。重建模块,用于对上述磁共振信号进行傅里叶变换,获取待扫描区域的磁共振图像。利用本磁共振成像系统,可在短时间内完成对射频场的校正,且获得的图像可有效减小由于B1场不均匀带来的图像伪影,提高获取临床图像的均匀性和对比度。此外由于对B1场分布进行了优化,则在一定程度上避免局部高场射频能量吸收率SAR过高的情况存在。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (9)

1.一种磁共振系统成像的方法,所述磁共振系统包括若干个射频发射通道,具体包括以下步骤:
在多通道发射模式下,采用包含三个射频脉冲的刺激回波序列获取自旋回波信号强度SSE和刺激回波信号强度SSTE
根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化以产生目标B1场,所述目标B1场同时满足B1场均匀性和射频发射功率优化;
将经过优化的幅度和相位参数加载到射频脉冲控制器,驱动多个射频发射通道产生目标B1场;
在所述目标B1场下激发受检者待扫描部位,获取待扫描部位的磁共振信号;
对所述磁共振信号进行傅里叶变换,获取待扫描部位的磁共振图像;
在对多个射频发射通道的幅度和相位参数进行优化以产生目标B1场之后,还包括:对目标B1场对应的射频发射功率PTar进行归一化,获取参考射频发射功率PRef,Pref=K×PTar,其中比例因子
Figure FDA0002770609200000011
其中θm表示目标B1场对应的目标翻转角,θref表示参考B1场对应的参考翻转角。
2.根据权利要求1所述的磁共振系统成像的方法,其特征在于,所述根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化的具体过程为:
(a)在磁共振系统空载状态下,将初始校正得到的若干个射频发射通道的幅度、相位参数和射频发射功率初始值加载到射频脉冲控制器;
(b)驱动所述多个射频发射通道并采用刺激回波序列获取成像区域内自旋回波信号强度SSE和刺激回波信号强度SSTE
(c)根据自旋回波信号强度SSE和刺激回波信号强度SSTE判断当前B1场是否为目标B1场,如果否,则调整所述若干个射频发射通道的幅度、相位参数,并返回步骤(b);否则,结束优化过程。
3.根据权利要求2所述的磁共振系统成像的方法,其特征在于,所述根据自旋回波信号强度SSE和刺激回波信号强度SSTE判断当前B1场是否为目标B1场的具体过程为:
根据自旋回波信号强度SSE和刺激回波信号强度SSTE获取当前翻转角θ;
获取目标B1场所对应的目标翻转角θm,计算所述当前翻转角θ与所述目标翻转角θm的差值是否在设定范围内,如果是,则判定当前B1场为目标B1场;如果否,则判定当前B1场不是目标B1场。
4.根据权利要求3所述的磁共振系统成像的方法,其特征在于,所述目标翻转角θm=90°。
5.根据权利要求3所述的磁共振系统成像的方法,其特征在于,所述刺激回波的三个射频脉冲对应翻转角分别为θ1、θ2、θ3,且θ1=θ2=θ3=θ,0≤θ≤180°,则
Figure FDA0002770609200000021
其中θ表示当前翻转角。
6.根据权利要求3所述的磁共振系统成像的方法,其特征在于,所述刺激回波的三个射频脉冲对应翻转角分别为θ1、θ2、θ3,且θ1=θ3=θ,θ2=2θ,0≤θ≤90°,则
Figure FDA0002770609200000022
其中θ表示当前翻转角。
7.根据权利要求2所述的磁共振系统成像的方法,其特征在于,所述根据自旋回波信号强度SSE和刺激回波信号强度SSTE判断当前B1场是否为目标B1场的具体过程为:
获取当前B1场下采集的刺激回波信号强度SSTE与自旋回波信号强度SSE的比值i,
Figure FDA0002770609200000031
获取目标B1场对应的刺激回波信号强度SSTE与自旋回波信号强度SSE的参考比值j,0≤j≤2;
如果i与j的差值在设定范围内,则判定当前B1场为目标B1场;否则,则判定当前B1场不是目标B1场。
8.一种磁共振系统,包括:
体发射线圈模块,包含多个射频线圈发射通道;
优化模块,用于根据成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE,对多个射频发射通道的幅度和相位参数进行优化,所述自旋回波信号强度SSE和刺激回波信号强度SSTE通过采用包含三个射频脉冲的刺激回波序列获得;
驱动模块,用于根据多个射频发射通道优化后的幅度和相位参数驱动所述多个射频发射通道产生目标B1场,所述目标B1场同时满足B1场均匀性和射频发射功率优化;
信号采集模块,用于在所述目标B1场下采集待扫描部位的磁共振信号;
重建模块,用于对所述磁共振信号进行傅里叶变换,获取待扫描区域的磁共振图像;
在对多个射频发射通道的幅度和相位参数进行优化以产生目标B1场之后,所述驱动模块还用于对目标B1场对应的射频发射功率PTar进行归一化,获取参考射频发射功率PRef,Pref=K×PTar,其中比例因子
Figure FDA0002770609200000041
其中θm表示目标B1场对应的目标翻转角,θref表示参考B1场对应的参考翻转角。
9.根据权利要求8所述的磁共振系统,其特征在于,所述优化模块根据待成像区域内采集的自旋回波信号强度SSE和刺激回波信号强度SSTE获取当前B1场对应的当前翻转角θ,根据当前翻转角θ与目标B1场对应的目标翻转角θm的差值确定是否对多个射频发射通道的幅度和相位参数进行优化;
或所述优化模块获取当前B1场下采集的刺激回波信号强度SSTE与自旋回波信号强度SSE的比值i,根据i和目标B1场对应的刺激回波信号强度与自旋回波信号强度的参考比值j,确定是否对射频发射通道的幅度和相位参数进行优化,其中,0≤i≤2,0≤j≤2。
CN201510992357.7A 2015-12-25 2015-12-25 磁共振系统及成像的方法 Active CN106918794B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510992357.7A CN106918794B (zh) 2015-12-25 2015-12-25 磁共振系统及成像的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510992357.7A CN106918794B (zh) 2015-12-25 2015-12-25 磁共振系统及成像的方法

Publications (2)

Publication Number Publication Date
CN106918794A CN106918794A (zh) 2017-07-04
CN106918794B true CN106918794B (zh) 2021-01-08

Family

ID=59459556

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510992357.7A Active CN106918794B (zh) 2015-12-25 2015-12-25 磁共振系统及成像的方法

Country Status (1)

Country Link
CN (1) CN106918794B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107390149B (zh) * 2017-07-21 2019-12-13 上海联影医疗科技有限公司 梯度线圈极性的检测方法、存储介质及磁共振扫描系统
CN107656222B (zh) * 2017-09-30 2020-09-18 上海联影医疗科技有限公司 磁共振扫描策略确定方法、磁共振扫描系统
CN107942397A (zh) * 2017-12-29 2018-04-20 吉林大学 用预极化场增强信号幅度的磁共振多通道探测方法及装置
DE102018102965A1 (de) * 2018-02-09 2019-08-14 Otto-Von-Guericke-Universität Magdeburg Hybridsystem für die Durchführung einer Magnetresonanztomographie und einer Radiofrequenzablation sowie Verfahren zu dessen Betrieb
EP3591418A1 (en) * 2018-07-03 2020-01-08 Koninklijke Philips N.V. Mri method for b0-mapping
US11175367B2 (en) * 2018-08-10 2021-11-16 General Electric Company Methods and systems for estimating transmit attenuation for a magnetic resonance imaging scan
CN109696646B (zh) * 2019-01-17 2020-10-30 华东师范大学 一种快速自旋回波脉冲序列中射频脉冲的优化方法
CN110068780B (zh) * 2019-05-07 2021-03-23 上海东软医疗科技有限公司 一种磁共振成像的自适应方法、装置和磁共振成像系统
CN110261801A (zh) * 2019-06-19 2019-09-20 上海东软医疗科技有限公司 一种射频激励方法及装置
CN110361682A (zh) * 2019-07-31 2019-10-22 张忠铭 快速确定磁共振射频激发翻转角度对应能量自动校准方法
CN113093078B (zh) * 2021-03-31 2021-11-30 无锡鸣石峻致医疗科技有限公司 一种射频脉冲强度与翻转角线性关系的确定方法、装置、计算机设备及存储介质
CN115561687B (zh) * 2022-10-13 2023-04-07 兰州大学 一种优化植入医疗器械在mri中射频发热的方法和系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN85103819A (zh) * 1985-05-11 1986-11-05 菲利浦光灯制造公司 确定人体内某一部分核磁共振分布的方法及其设备
DE69431894D1 (de) * 1993-05-12 2003-01-30 Marconi Medical Sys Inc Gerät und Verfahren zum Anregen von Kernspinresonanz in einer begrenzten Unterregion einer Untersuchungsregion
CN101327123A (zh) * 2007-06-20 2008-12-24 通用电器横河医疗系统株式会社 磁共振成像装置和磁共振图像产生方法
CN102159966A (zh) * 2008-09-17 2011-08-17 皇家飞利浦电子股份有限公司 利用用于mri的多信道rf发射系统的具有rf功率正则化的rf匀场
CN103597370A (zh) * 2011-04-21 2014-02-19 马普协会 空间编码的相位对比磁共振成像
CN103814303A (zh) * 2011-07-08 2014-05-21 明尼苏达大学董事会 对径向采样图像数据集中的角向一致伪影的多步校正
CN104062611A (zh) * 2013-03-22 2014-09-24 西门子(深圳)磁共振有限公司 一种磁共振系统的射频激发方法和装置
CN104520728A (zh) * 2012-08-08 2015-04-15 皇家飞利浦有限公司 具有重影稳定性校正的多次激发磁共振成像

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4049649B2 (ja) * 2002-10-01 2008-02-20 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴撮影装置
AU2010273298B2 (en) * 2009-07-15 2014-10-23 Viewray Technologies, Inc. Method and apparatus for shielding a linear accelerator and a magnetic resonance imaging device from each other
DE102011005649B4 (de) * 2011-03-16 2013-07-04 Siemens Aktiengesellschaft Bestimmung des tatsächlichen Kippwinkels und Justierung der Transmitterspannung bei der MR-Bildgebung eines kontinuierlich verfahrenen Untersuchungsobjekts
CN104799859B (zh) * 2014-01-27 2019-01-08 西门子(深圳)磁共振有限公司 一种磁共振弥散加权成像图像的校正方法和装置
CN104515963B (zh) * 2014-12-19 2017-02-08 中国科学院深圳先进技术研究院 磁共振射频匀场系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN85103819A (zh) * 1985-05-11 1986-11-05 菲利浦光灯制造公司 确定人体内某一部分核磁共振分布的方法及其设备
DE69431894D1 (de) * 1993-05-12 2003-01-30 Marconi Medical Sys Inc Gerät und Verfahren zum Anregen von Kernspinresonanz in einer begrenzten Unterregion einer Untersuchungsregion
CN101327123A (zh) * 2007-06-20 2008-12-24 通用电器横河医疗系统株式会社 磁共振成像装置和磁共振图像产生方法
CN102159966A (zh) * 2008-09-17 2011-08-17 皇家飞利浦电子股份有限公司 利用用于mri的多信道rf发射系统的具有rf功率正则化的rf匀场
CN103597370A (zh) * 2011-04-21 2014-02-19 马普协会 空间编码的相位对比磁共振成像
CN103814303A (zh) * 2011-07-08 2014-05-21 明尼苏达大学董事会 对径向采样图像数据集中的角向一致伪影的多步校正
CN104520728A (zh) * 2012-08-08 2015-04-15 皇家飞利浦有限公司 具有重影稳定性校正的多次激发磁共振成像
CN104062611A (zh) * 2013-03-22 2014-09-24 西门子(深圳)磁共振有限公司 一种磁共振系统的射频激发方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
磁共振快速自旋回波技术及临床应用;戴捷 等;《实用医技杂志》;20040831;第11卷(第8期);第1492-1493页 *

Also Published As

Publication number Publication date
CN106918794A (zh) 2017-07-04

Similar Documents

Publication Publication Date Title
CN106918794B (zh) 磁共振系统及成像的方法
US7649353B2 (en) Magnetic resonance imaging method and apparatus using a circularly polarized RF resonator detecting right and left circularly polarized components of MR signals
US8319496B2 (en) Magnetic resonance method and apparatus for reducing RF heating in the patient
US9547064B2 (en) Method for the control of a magnetic resonance system
US9095302B2 (en) Method to generate magnetic resonance exposures
US10234528B2 (en) Method and apparatus to correct noise effects in quantitative techniques in magnetic resonance imaging
US7769425B2 (en) Magnetic resonance device and method
US7265545B2 (en) Method and apparatus for accelerated spiral-coded imaging in magnetic resonance tomography
US7282914B2 (en) Specific energy absorption rate model
US9335395B2 (en) Magnetic resonance system and method for frequency calibration of the magnetic resonance system
US10018691B2 (en) Correction of imaging methods in a magnetic resonance device
WO2014050818A1 (ja) 磁気共鳴イメージング装置およびその高周波出力測定装置
US9494662B2 (en) Magnetic resonance method and apparatus for automatic calculation of a maximum pulse-length of an excitation pulse
US6734673B2 (en) Apparatus and method for magnetic resonance imaging employing partial parallel acquisition with formation of signal combinations and solving of an equation system based thereon
CN107440718B (zh) 磁共振成像的sar值确定方法及磁共振成像装置
US11698427B2 (en) Methods and systems for estimating transmit attenuation for a magnetic resonance imaging scan
US9678188B2 (en) Magnetic resonance system and method to acquire at least two measurement data sets from an examination subject
CN110811619B (zh) 用于估计磁共振成像扫描的发射衰减的方法和系统
US20050237059A1 (en) Method and apparatus for suppression of artifacts in MRT imaging
Runderkamp et al. Whole‐liver flip‐angle shimming at 7 T using parallel‐transmit k T‐point pulses and Fourier phase‐encoded DREAM B 1+ mapping
US20220206097A1 (en) Magnetic resonance image processing method and device, and computer readable storage medium
US11698432B2 (en) Magnetic resonance imaging system, and main magnetic field correction method therefor and storage medium
CN115963439A (zh) 用于估计磁共振成像扫描的发射衰减的方法和系统
Runderkamp et al. Whole-liver flip-angle shimming at 7 T using parallel-transmit kT-point pulses and Fourier phase-encoded DREAM B
CN114966504A (zh) Mri系统、确定mri系统sar值的方法、计算机可读存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258

Applicant after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant