CN113030813A - 一种磁共振t2定量成像方法及系统 - Google Patents

一种磁共振t2定量成像方法及系统 Download PDF

Info

Publication number
CN113030813A
CN113030813A CN202110215919.2A CN202110215919A CN113030813A CN 113030813 A CN113030813 A CN 113030813A CN 202110215919 A CN202110215919 A CN 202110215919A CN 113030813 A CN113030813 A CN 113030813A
Authority
CN
China
Prior art keywords
echo
sample
scanning
image
neural network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110215919.2A
Other languages
English (en)
Other versions
CN113030813B (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.)
Xiamen University
Original Assignee
Xiamen University
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 Xiamen University filed Critical Xiamen University
Priority to CN202110215919.2A priority Critical patent/CN113030813B/zh
Publication of CN113030813A publication Critical patent/CN113030813A/zh
Application granted granted Critical
Publication of CN113030813B publication Critical patent/CN113030813B/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/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4816NMR imaging of samples with ultrashort relaxation times such as solid samples, e.g. MRI using ultrashort TE [UTE], single point imaging, constant time imaging

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明提供了一种磁共振T2定量成像方法及系统,涉及磁共振成像领域。方法包括:设计多扫描多重叠回波T2定量成像脉冲序列并确定其采样参数;利用脉冲序列对待测样品进行多次扫描采样并对采样得到的信号进行处理得到待重建多扫描多重叠回波图像;根据待测样品特征生成设定量的训练样本;确定用于T2定量图像重建的深度神经网络并采用设定量的训练样本对该深度神经网络进行训练;将待测样品的待重建多扫描多重叠回波图像输入训练好的深度神经网络进行重建,得到待测样品的高分辨的T2定量图像。

Description

一种磁共振T2定量成像方法及系统
技术领域
本发明涉及磁共振成像领域,特别是涉及一种磁共振T2定量成像方法及系统。
背景技术
磁共振成像(Magnetic Resonance Imaging,MRI)是利用原子核在强磁场内发生共振产生的信号经图像重建的一种成像技术。MRI可分为定性的参数加权成像和定量的参数成像。虽然参数加权成像可以提供良好的组织对比度,但是由于图像的对比度与获得图像所用的脉冲序列以及扫描参数有很大的关系,因此在不同设备上采集到的图像组织对比度往往存在差异,难以进行横向和纵向对比。定量参数成像可以直接测量组织的特性参数值,减少仪器设备条件对成像结果的影响,因此具有重要的应用价值。
现有的T2定量成像往往需要较长的扫描时间,成像结果易受生理运动等的影响,导致图像分辨率降低。为了缩短定量成像的扫描时间,已有多种方法被提出。例如,磁共振指纹(MRF)方法可实现快速定量T2成像,但其图像重建时间长。已经提出的单扫描多重叠回波方法可实现超快速定量T2成像,然而,其图像分辨率有限,一些结构信息可能会丢失。因此,本发明用多扫描的方式改善单扫描多重叠回波采样,以提高图像分辨率。在使用多扫描的方式时,本发明设计了一个预相位编码梯度,使每次扫描采集的信号合并后可以得到样品完整的多重叠回波信号。同时,由于多扫描的方式会带来更严重的涡流伪影,影响定量成像效果,本发明提出一种适用于多扫描多重叠回波成像的涡流伪影去除方法。
发明内容
基于以上背景,本发明提供了一种磁共振T2定量成像方法及系统。通过设计多扫描多重叠回波T2定量成像脉冲序列,缩短单次扫描的采样回波链长度,在保证信噪比的同时提高了图像的分辨率;设计了一种适用于多扫描多重叠回波图像的涡流伪影去除方法,有效去除多扫描多重叠回波图像的涡流伪影;将去涡流伪影后的多扫描多重叠回波图像输入到训练好的深度神经网络中进行重建,即可获得高分辨T2定量图像。
为实现上述目的,本发明提供了如下方案:
一种磁共振T2定量成像方法,包括以下步骤:
设计磁共振T2定量成像脉冲序列,称为多扫描多重叠回波T2定量成像脉冲序列;
确定多扫描多重叠回波T2定量成像脉冲序列采样参数;
利用所述脉冲序列在设定的采样参数下对待测样品进行多次扫描采样,得到待测样品的多扫描多重叠回波信号;
对所述待测样品的多扫描多重叠回波信号进行处理,得到待测样品的待重建多扫描多重叠回波图像;
根据待测样品特征模拟生成设定量的深度神经网络训练样本,构成训练样本集;
确定用于T2定量图像重建的深度神经网络;
采用所述训练样本集对所述深度神经网络进行训练,获得训练好的深度神经网络;
将所述待测样品的待重建多扫描多重叠回波图像输入所述训练好的深度神经网络进行重建,获得待测样品的高分辨T2定量图像。
本发明还提供了一种磁共振T2定量成像系统,包括:
脉冲序列设计模块,用于设计磁共振T2定量成像脉冲序列,所述磁共振T2定量成像脉冲序列称为多扫描多重叠回波T2定量成像脉冲序列;确定多扫描多重叠回波T2定量成像脉冲序列的采样参数;
信号采集模块,用于利用所述脉冲序列在设定的采样参数下对待测样品进行多次扫描采样,得到待测样品的多扫描多重叠回波信号;
信号处理模块,用于对所述待测样品的多扫描多重叠回波信号进行处理,得到待测样品的待重建多扫描多重叠回波图像;
训练样本集生成模块,用于将所述脉冲序列输入磁共振成像仿真软件中,并且依据真实实验的非理想性添加相应的非理想项,以尽量模拟真实情况;根据所述待测样品特征生成设定量的随机模板;利用仿真软件对各随机模板进行模拟采样,得到各随机模板的多扫描多重叠回波信号;将各随机模板的多扫描多重叠回波信号重排成二维k空间信号,然后进行二维傅里叶变换,得到各随机模板的多扫描多重叠回波图像;将各随机模板的多扫描多重叠回波图像和相应的随机模板构成一个训练样本,得到设定量的训练样本,构成训练样本集。
深度神经网络确定模块,用于确定用于T2定量图像重建的深度神经网络;采用所述训练样本集对所述深度神经网络进行训练,获得训练好的深度神经网络;
T2定量图像重建模块,用于将所述待测样品的待重建多扫描多重叠回波图像输入所述训练好的深度神经网络进行重建,获得待测样品的T2定量图像。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
1)与单扫描多重叠回波成像相比,本发明设计的多扫描多重叠回波成像脉冲序列缩短了单次扫描的采样回波链长度,在保证信噪比的同时提高了图像的分辨率;
2)通过设计适用于多扫描多重叠回波图像的涡流伪影去除方法,有效去除多扫描多重叠回波图像的涡流伪影;
3)在模拟生成深度神经网络的训练样本集阶段,考虑到实际实验中重聚脉冲的非理想性,对重聚脉冲引入了不相干的随机平滑的等效重聚角度图以尽量模拟真实情况,减小了重建的T2定量图像的误差;
4)利用深度神经网络对去除涡流伪影后的多扫描多重叠回波图像进行T2定量图像重建,比传统重建方法效率更高,使用更方便。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例1中磁共振T2定量成像方法的流程图;
图2为本发明实施例1中脉冲序列的结构示意图;
图3为本发明实施例2中深度神经网络的结构示意图;
图4为本发明实施例3中磁共振T2定量成像系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提出了一种磁共振T2定量成像方法及系统。通过设计多扫描多重叠回波成像脉冲序列,有效缩短每次扫描的采样回波链长度,在保证信噪比的同时提高了图像的分辨率;通过设计适用于多扫描多重叠回波图像的涡流伪影去除方法,有效解决涡流导致的伪影,进一步提升图像质量;利用深度神经网络对去涡流伪影后的图像进行处理,可直接获得高分辨的T2定量图像。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1:
如图1所示,为本发明实施例1中磁共振T2定量成像方法的流程图。S1-S8表示方法的各步骤:
S1:设计磁共振T2定量成像脉冲序列,称为多扫描多重叠回波T2定量成像脉冲序列;
S2:确定多扫描多重叠回波T2定量成像脉冲序列的采样参数;
S3:利用所述脉冲序列在设定的采样参数下对待测样品进行多次扫描采样,得到待测样品的多扫描多重叠回波信号;
S4:对所述待测样品的多扫描多重叠回波信号进行处理,得到待测样品的待重建多扫描多重叠回波图像;
S5:根据待测样品模拟生成设定量的深度神经网络训练样本,构成训练样本集;
S6:确定用于T2定量图像重建的深度神经网络;
S7:采用所述训练样本集对所述深度神经网络进行训练,获得训练好的深度神经网络;
S8:将所述待测样品的待重建多扫描多重叠回波图像输入所述训练好的深度神经网络进行重建,获得待测样品的高分辨T2定量图像。
接下来将对上述各步骤做详细描述:
在设计成像所需的脉冲序列时,考虑到单扫描多重叠回波成像要获得高的图像分辨率需要较长的采样回波链长度,导致图像信噪比降低,因此本发明提出多扫描方法,通过多次扫描缩短脉冲序列中采样回波链长度,从而在保证图像信噪比的同时提高图像的分辨率。
如图2所示,为本发明实施例1中磁共振T2定量成像脉冲序列结构示意图。
其中,αi表示第i个射频激发脉冲的翻转角,Gi表示第i个射频激发脉冲对应的回波移位梯度;β表示重聚脉冲翻转角;Gcr表示破坏梯度;Gpre表示预相位编码梯度;Echoi表示第i个回波重聚的中心位置,TEi表示第i个回波的回波时间,其中,i=1,2,…,n,n为射频激发脉冲个数,决定了回波的个数;j表示第j次扫描,j+1表示第j+1次扫描,TR表示第j次扫描和第j+1次扫描的时间间隔,j=1,2,…,J,J为扫描次数;回波移位梯度包括频率编码维的移位梯度和相位编码维的移位梯度;
破坏梯度包括频率编码维、相位编码维和层面选择维的破坏梯度;
采样回波链包括频率编码维的梯度链和相位编码维的梯度链;
其中,m1+m2+…+mn+mn+1+2n=N,N为采样回波链中相位编码步数,ml、m2、…、mn和mn+1分别表示图2中各虚线括号内包含的相位编码步数;预相位编码梯度施加在相位编码维方向,用以改变每次扫描的k空间填充起点,每次扫描取值不同。
确定多扫描多重叠回波成像脉冲序列的采样参数,具体包括:
确定射频激发脉冲的个数n和各射频激发脉冲的翻转角αi
确定重聚脉冲的翻转角β;
确定各回波的重叠方式,即各回波在k空间中的位置,从而确定各射频激发脉冲之间的时间间隔比例以及各射频激发脉冲之后的回波移位梯度的大小比例;
确定脉冲序列扫描次数、成像视野、成像矩阵以及回波间隔,从而确定采样回波链长度与采样回波链中频率编码维和相位编码维的梯度值、预相位编码梯度值、各射频脉冲之间的时间间隔与各回波移位梯度值;
确定脉冲序列重复时间TR。
在利用所述脉冲序列得到多扫描多重叠回波图像后,本发明考虑到:由于所述脉冲序列的采样回波链中频率编码维上梯度的快速切换易产生涡流,使采集的信号产生相位累积误差,从而导致多扫描多重叠回波图像出现涡流伪影,降低图像质量,为此,本发明提出了一个相位校正方法,具体方案如下:
设定一个待定系数的线性相位模型,并设置其中各个参数的取值范围和取值步长;
将所述待测样品的多扫描多重叠回波信号重排成二维k空间信号S(kx,ky),其中kx为k空间频率编码维坐标,ky为k空间相位编码维坐标;
对S(kx,ky)沿频率编码维进行一维傅里叶变换,得到I(x,ky),其中x为图像域频率编码维坐标;
依照设定的线性相位模型向I(x,ky)添加相位,得到I′(x,ky),其中沿频率编码维方向添加的相位与x成线性关系,沿相位编码维方向添加的相位大小相等但同一次扫描采集的相邻行之间添加的相位符号相反,从而构成奇偶相位差;
对I′(x,ky)沿相位编码维进行一维傅里叶变换,得到添加相位后待测样品的多扫描多重叠回波图像;
将添加相位后的多扫描多重叠回波图像与将添加相位后的多扫描多重叠回波图像按扫描次数进行循环位移和叠加后得到的图像进行比较,计算添加相位后的多扫描多重叠回波图像的涡流伪影程度;
在设定的次数内重复以上过程,找到在线性相位模型参数取值范围内,添加相位后涡流伪影程度最小的多扫描多重叠回波图像,该图像即为涡流伪影被最大程度地消除了的待测样品的待重建多扫描多重叠回波图像。
在确定深度神经网络前,本发明考虑到深度神经网络需要大量的训练样本,而真实的训练样本难以获得,因此,本发明利用模拟方法生成深度神经网络所需的训练样本,具体步骤如下:
将所述脉冲序列输入磁共振成像仿真软件中,依据真实实验的非理想性添加相应的非理想项,以尽量模拟真实情况;
根据待测样品特征生成设定量的随机模板;
利用仿真软件对各随机模板进行模拟采样,得到各随机模板的多扫描多重叠回波信号;将各随机模板的多扫描多重叠回波信号重排成二维k空间信号,然后进行二维傅里叶变换,得到各随机模板的多扫描多重叠回波图像;
将各随机模板的多扫描多重叠回波图像和相应的随机模板构成一个训练样本,得到设定量的训练样本,构成训练样本集。
此外,本发明考虑到实际实验中脉冲序列的重聚脉冲往往角度不准,因此在利用磁共振成像仿真软件对随机模板进行模拟采样前,对重聚脉冲引入不相干的随机平滑的等效重聚角度图,以保证获得的训练样本尽可能与真实情况相符,降低最终重建结果的误差。
接下来对本发明中确定用于T2定量图像重建的深度神经网络的方法做具体描述。方法包括:
确定深度神经网络的网络结构;
确定深度神经网络的输入和输出通道数;
确定训练深度神经网络的损失函数。
训练深度神经网络时,将训练样本集分批次输入深度神经网络进行迭代训练,每次训练网络计算出损失函数的值,根据该值自动调整神经网络的参数值使损失函数的值减小,重复上述训练直至损失函数的值不再减小,保存深度神经网络参数。
深度神经网络训练完成后,将去涡流伪影后的多扫描多重叠回波图像作为深度神经网络的输入,直接得到高分辨的T2定量图像。
基于上述,本发明在实施例1中设计好了成像所需的脉冲序列,训练好了图像重建所需的深度神经网络,并提出了一些为提升T2定量图像分辨率的方法。
在对真实样品进行T2定量成像时,可直接基于本发明提出的磁共振T2定量成像方法,利用设计好的脉冲序列对样品进行采样得到多扫描多重叠回波图像,对多扫描多重叠回波图像做去涡流伪影处理,利用训练好的深度神经网络对去涡流伪影后的多扫描多重叠回波图像进行处理,得到高分辨T2定量图像。
实施例2:
如图3所示,为本发明实施例2中深度神经网络的结构示意图。其中的1表示去除涡流伪影后的多扫描多重叠回波图像,2表示重建得到的T2定量图像。
该深度神经网络是一个五层的U-Net网络,可以看作是对称的编-解码器结构,包括4个编码处理单元(图3的左半部分)和4个解码处理单元(图3的右半部分);
每个编码处理单元包括依次串联的第一卷积模块和最大池化模块;每个最大池化模块做下采样处理;每个解码处理单元包括依次串联的反卷积模块和第二卷积模块;每个反卷积模块做上采样处理。
需要说明的是,本发明中深度神经网络不限于U-Net,可以是其它网络;编码处理单元和解码处理单元的数量也不限定为4个,可为任意正整数个。只要能够实现本发明如上述的深度神经网络的功能,都在本发明的保护范围内。
更具体的,本发明实施例2中的第一卷积模块包括依次串联的设定量的第一卷积层和一个Relu激活函数;最大池化模块包括一个最大池化层;
第二卷积模块包括依次串联的一个反卷积层和设定量的第二卷积层;
第三卷积模块包括一个第三卷积层;
第一卷积层和第二卷积层中包含设定量的卷积核,第三卷积层中包含一个卷积核;卷积核的大小均为3×3、步长均为1、padding形式均为same;
最大池化层中包含设定量的卷积核,卷积核的大小为2×2、步长为2、padding形式为same;
反卷积层中包含设定量的卷积核,卷积核的大小为2×2、步长为2。
另外,本发明实施例2中还对编码处理单元的第一卷积层中卷积核的数量以及解码处理单元的第二卷积层中卷积核的数量进行了限定。
其中,后一个编码处理单元的第一卷积层中卷积核的数量是当前编码处理单元的第一卷积层中卷积核的数量的二倍;
后一个解码处理单元的第二卷积层中卷积核的数量是当前解码处理单元的第二卷积层中卷积核的数量的二分之一。
图3中的64、128、256、512和1024表示的就是卷积层中卷积核的数量。
而且,由图3可以看出,随着下采样的进行,编码处理单元的侧面积逐渐减小为上一个编码处理单元的二分之一,表示编码处理单元中特征图的面积依次减小为上一个特征图面积的二分之一;随着上采样的进行,解码处理单元的侧面积逐渐增大为上一个解码处理单元的二倍,表示解码处理单元中特征图的面积依次增大为上一个特征图面积的二倍。
需要说明的是,本发明中各卷积模块中卷积层的数量可为任意正整数个,各卷积层的卷积核数量也不仅限于实施例2中所示的64、128、256、512和1024,只要编码处理单元与解码处理单元的数量相等,且后一个编码处理单元的第一卷积层中卷积核的数量是当前编码处理单元的第一卷积层中卷积核的数量的二倍;后一个解码处理单元的第二卷积层中卷积核的数量是当前解码处理单元的第二卷积层中卷积核的数量的二分之一,卷积层中卷积核的各种数量设置均在本发明的保护范围内。
深度神经网络的损失函数为:
Figure BDA0002953781620000091
其中,L表示损失函数的值,M表示深度神经网络训练样本集中训练样本的数量,||||表示范数,f()表示深度神经网络的映射关系,xk表示第k个训练样本中的图像,W和b表示深度神经网络的参数,yk表示第k个训练样本中的随机模板,ychange表示将yk中小于设定阈值的值置为设定阈值后的矩阵,F表示范数的类型,λ表示约束项系数,▽表示梯度算子,ymask表示对yk使用Canny算子求得的图像边缘信息;
采用指数衰减的方式调整学习率以控制所述深度神经网络的训练次数:设定初始学习率为0.0001,每达到设定训练次数后,所述学习率降低一次。
训练深度神经网络时,将训练样本集分批次输入深度神经网络进行迭代训练,每次训练网络计算出损失函数的值,根据该值自动调整W和b的值使损失函数的值减小,重复上述训练直至损失函数的值不再减小,保存深度神经网络参数W和b。
实施例3:
如图4所示,本发明还提供了一种磁共振T2定量成像系统,该系统基于如实施例1中的磁共振T2定量成像方法实现。
具体的,系统包括:
脉冲序列设计模块,用于设计磁共振T2定量成像脉冲序列,该脉冲序列称为多扫描多重叠回波成像脉冲序列;确定多扫描多重叠回波成像脉冲序列的采样参数;
信号采集模块,用于利用脉冲序列在设定的采样参数下对待测样品进行多次扫描采样,得到待测样品的多扫描多重叠回波信号;
信号处理模块,用于对待测样品的多扫描多重叠回波信号进行处理,得到待测样品的待重建多扫描多重叠回波图像;
训练样本集生成模块,用于将脉冲序列输入磁共振成像仿真软件中,并且依据真实实验的非理想性添加相应的非理想项,以尽量模拟真实情况;根据待测样品特征生成设定量的随机模板;利用仿真软件对各随机模板进行模拟采样,得到各随机模板的多扫描多重叠回波信号;将各随机模板的多扫描多重叠回波信号重排成二维k空间信号,然后进行二维傅里叶变换,得到各随机模板的多扫描多重叠回波图像;将各随机模板的多扫描多重叠回波图像和相应的随机模板构成一个训练样本,得到设定量的训练样本,构成训练样本集。
深度神经网络确定模块,用于确定用于T2定量图像重建的深度神经网络;采用训练样本集对所述深度神经网络进行训练,获得训练好的深度神经网络;
T2定量图像重建模块,用于将待测样品的待重建多扫描多重叠回波图像输入训练好的深度神经网络进行重建,获得待测样品的T2定量图像。
综上所述,本发明中提供的磁共振T2定量成像方法及系统,具备以下技术效果:
1)与单扫描多重叠回波成像相比,本发明设计的多扫描多重叠回波成像脉冲序列,缩短了单次扫描的采样回波链长度,在保证信噪比的同时提高了图像的分辨率;
2)通过设计适用于多扫描多重叠回波图像的涡流伪影去除方法,有效去除多扫描多重叠回波图像的涡流伪影;
3)在模拟生成深度神经网络的训练样本集阶段,考虑到实际实验中重聚脉冲的非理想性,对重聚脉冲引入了不相干的随机平滑的等效重聚角度图以尽量模拟真实情况,减小了重建的T2定量图像的误差;
4)利用深度神经网络对去除涡流伪影后的多扫描多重叠回波图像进行T2定量图像重建,比传统重建方法效率更高,使用更方便。
因此,利用本发明提供的方法及系统,能够获得高分辨率的磁共振T2定量图像。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (7)

1.一种磁共振T2定量成像方法,其特征在于,包括以下步骤:
设计磁共振T2定量成像脉冲序列,称为多扫描多重叠回波T2定量成像脉冲序列;
确定多扫描多重叠回波T2定量成像脉冲序列的采样参数;
利用所述脉冲序列在设定的采样参数下对待测样品进行多次扫描采样,得到待测样品的多扫描多重叠回波信号;
对所述待测样品的多扫描多重叠回波信号进行处理,得到待测样品的待重建多扫描多重叠回波图像;
根据待测样品特征模拟生成设定量的深度神经网络训练样本,构成训练样本集;
确定用于T2定量图像重建的深度神经网络;
采用所述训练样本集对所述深度神经网络进行训练,获得训练好的深度神经网络;
将所述待测样品的待重建多扫描多重叠回波图像输入所述训练好的深度神经网络进行重建,获得待测样品的T2定量图像。
2.根据权利要求1所述的磁共振T2定量成像方法,其特征在于,所述多扫描多重叠回波T2定量成像脉冲序列包括:
翻转角为αi的射频激发脉冲和对应的回波移位梯度Gi、翻转角为β的重聚脉冲、破坏梯度Gcr、预相位编码梯度Gpre、采样回波链;
其中,i=1,2,...,n;n为射频激发脉冲个数,决定了回波的个数;
每个射频激发脉冲都与层面选择维的层选梯度相结合进行层选;
回波移位梯度施加在频率编码维和相位编码维方向;
破坏梯度施加在频率编码维、相位编码维和层面选择维方向;
采样回波链由分别作用在频率编码维和相位编码维的梯度链组成,频率编码维的梯度链由一系列正负梯度构成,相位编码维的梯度链由一系列面积相等的梯度构成;
预相位编码梯度施加在相位编码维方向,用以改变每次扫描的k空间填充起点,每次扫描取值不同。
3.根据权利要求2所述的磁共振T2定量成像方法,其特征在于,确定多扫描多重叠回波T2定量成像脉冲序列的采样参数,具体包括:
确定射频激发脉冲的个数n和各射频激发脉冲的翻转角αi
确定重聚脉冲的翻转角β;
确定各回波的重叠方式,即各回波在k空间中的位置,从而确定各射频激发脉冲之间的时间间隔比例以及各射频激发脉冲之后的回波移位梯度的大小比例;
确定脉冲序列扫描次数、成像视野、成像矩阵以及回波间隔,从而确定采样回波链长度与采样回波链中频率编码维和相位编码维的梯度值、预相位编码梯度值、各射频脉冲之间的时间间隔与各回波移位梯度值;
确定脉冲序列重复时间TR。
4.根据权利要求1所述的磁共振T2定量成像方法,其特征在于,对所述待测样品的多扫描多重叠回波信号进行处理,得到待测样品的待重建多扫描多重叠回波图像,具体包括:
设定一个待定系数的线性相位模型,并设置其中各个参数的取值范围和取值步长;
将所述待测样品的多扫描多重叠回波信号重排成二维k空间信号S(kx,ky),其中kx为k空间频率编码维坐标,ky为k空间相位编码维坐标;
对S(kx,ky)沿频率编码维进行一维傅里叶变换,得到I(x,ky),其中x为图像域频率编码维坐标;
依照设定的线性相位模型向I(x,ky)添加相位,得到I′(x,ky),其中沿频率编码维方向添加的相位与x成线性关系,沿相位编码维方向添加的相位大小相等但同一次扫描采集的相邻行之间添加的相位符号相反,从而构成奇偶相位差;
对I′(x,ky)沿相位编码维进行一维傅里叶变换,得到添加相位后待测样品的多扫描多重叠回波图像;
将添加相位后的多扫描多重叠回波图像与将添加相位后的多扫描多重叠回波图像按扫描次数进行循环位移和叠加后得到的图像进行比较,计算添加相位后的多扫描多重叠回波图像的涡流伪影程度;
在设定的次数内重复以上过程,找到在线性相位模型参数取值范围内,添加相位后涡流伪影程度最小的多扫描多重叠回波图像,该图像即为涡流伪影被最大程度地消除了的待测样品的待重建多扫描多重叠回波图像。
5.根据权利要求1所述的磁共振T2定量成像方法,其特征在于,根据待测样品特征模拟生成设定量的深度神经网络训练样本,构成训练样本集,具体包括:
将所述脉冲序列输入磁共振成像仿真软件中,并且依据真实实验的非理想性添加相应的非理想项,以尽量模拟真实情况;
根据所述待测样品特征生成设定量的随机模板;
利用仿真软件对各随机模板进行模拟采样,得到各随机模板的多扫描多重叠回波信号;
将各随机模板的多扫描多重叠回波信号重排成二维k空间信号,然后进行二维傅里叶变换,得到各随机模板的多扫描多重叠回波图像;
将各随机模板的多扫描多重叠回波图像和相应的随机模板构成一个训练样本,得到设定量的训练样本,构成训练样本集。
6.根据权利要求1所述的磁共振T2定量成像方法,其特征在于,所述确定用于T2定量图像重建的深度神经网络,具体包括:
确定深度神经网络的网络结构;
确定深度神经网络的输入和输出通道数;
确定训练深度神经网络的损失函数。
7.一种磁共振T2定量成像系统,其特征在于,包括:
脉冲序列设计模块,用于设计磁共振T2定量成像脉冲序列,所述磁共振T2定量成像脉冲序列称为多扫描多重叠回波T2定量成像脉冲序列;确定多扫描多重叠回波成像脉冲序列的采样参数;
信号采集模块,用于利用所述脉冲序列在设定的采样参数下对待测样品进行多次扫描采样,得到待测样品的多扫描多重叠回波信号;
信号处理模块,用于对所述待测样品的多扫描多重叠回波信号进行处理,得到待测样品的待重建多扫描多重叠回波图像;
训练样本集生成模块,用于将所述脉冲序列输入磁共振成像仿真软件中,并且依据真实实验的非理想性添加相应的非理想项,以尽量模拟真实情况;根据所述待测样品特征生成设定量的随机模板;利用仿真软件对各随机模板进行模拟采样,得到各随机模板的多扫描多重叠回波信号;将各随机模板的多扫描多重叠回波信号重排成二维k空间信号,然后进行二维傅里叶变换,得到各随机模板的多扫描多重叠回波图像;将各随机模板的多扫描多重叠回波图像和相应的随机模板构成一个训练样本,得到设定量的训练样本,构成训练样本集;
深度神经网络确定模块,用于确定用于T2定量图像重建的深度神经网络;采用所述训练样本集对所述深度神经网络进行训练,获得训练好的深度神经网络;
T2定量图像重建模块,用于将所述待测样品的待重建多扫描多重叠回波图像输入所述训练好的深度神经网络进行重建,获得待测样品的T2定量图像。
CN202110215919.2A 2021-02-26 2021-02-26 一种磁共振t2定量成像方法及系统 Active CN113030813B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110215919.2A CN113030813B (zh) 2021-02-26 2021-02-26 一种磁共振t2定量成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110215919.2A CN113030813B (zh) 2021-02-26 2021-02-26 一种磁共振t2定量成像方法及系统

Publications (2)

Publication Number Publication Date
CN113030813A true CN113030813A (zh) 2021-06-25
CN113030813B CN113030813B (zh) 2022-02-22

Family

ID=76461751

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110215919.2A Active CN113030813B (zh) 2021-02-26 2021-02-26 一种磁共振t2定量成像方法及系统

Country Status (1)

Country Link
CN (1) CN113030813B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113598716A (zh) * 2021-08-24 2021-11-05 南方医科大学 一种光声与磁共振联合成像方法
CN115166612A (zh) * 2022-07-21 2022-10-11 厦门大学 基于深度学习的超快速磁共振成像高精度仿真方法及系统

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1933772A (zh) * 2004-03-08 2007-03-21 内尔科尔普里坦贝内特公司 对组织中静脉与动脉混合血脉动光学检测的方法和装置
CN101470178A (zh) * 2007-12-29 2009-07-01 西门子(中国)有限公司 一种抑制残余运动伪影的方法及装置
CN101657732A (zh) * 2007-03-23 2010-02-24 居里研究所 在nmr弛豫测量成像中确定相位修正振幅的方法
CN104865545A (zh) * 2014-02-21 2015-08-26 上海联影医疗科技有限公司 平面回波成像方法及装置
CN105232045A (zh) * 2015-11-26 2016-01-13 厦门大学 基于双回波的单扫描定量磁共振扩散成像方法
CN109581319A (zh) * 2019-01-24 2019-04-05 西安电子科技大学 基于多扫描递归的海杂波多普勒偏移和带宽估计方法
CN109741409A (zh) * 2018-11-30 2019-05-10 厦门大学 回波平面成像涡流伪影的无参考扫描校正方法
WO2019224800A1 (en) * 2018-05-25 2019-11-28 Mahajan Vidur Method and system for simulating and constructing original medical images from one modality to other modality
CN110782452A (zh) * 2019-11-05 2020-02-11 厦门大学 一种t2定量图像成像方法及系统
CN113075598A (zh) * 2020-01-03 2021-07-06 上海联影医疗科技股份有限公司 磁共振系统的磁场分布检测方法、磁共振系统及成像方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1933772A (zh) * 2004-03-08 2007-03-21 内尔科尔普里坦贝内特公司 对组织中静脉与动脉混合血脉动光学检测的方法和装置
CN101657732A (zh) * 2007-03-23 2010-02-24 居里研究所 在nmr弛豫测量成像中确定相位修正振幅的方法
CN101470178A (zh) * 2007-12-29 2009-07-01 西门子(中国)有限公司 一种抑制残余运动伪影的方法及装置
CN104865545A (zh) * 2014-02-21 2015-08-26 上海联影医疗科技有限公司 平面回波成像方法及装置
CN105232045A (zh) * 2015-11-26 2016-01-13 厦门大学 基于双回波的单扫描定量磁共振扩散成像方法
WO2019224800A1 (en) * 2018-05-25 2019-11-28 Mahajan Vidur Method and system for simulating and constructing original medical images from one modality to other modality
CN109741409A (zh) * 2018-11-30 2019-05-10 厦门大学 回波平面成像涡流伪影的无参考扫描校正方法
CN109581319A (zh) * 2019-01-24 2019-04-05 西安电子科技大学 基于多扫描递归的海杂波多普勒偏移和带宽估计方法
CN110782452A (zh) * 2019-11-05 2020-02-11 厦门大学 一种t2定量图像成像方法及系统
CN113075598A (zh) * 2020-01-03 2021-07-06 上海联影医疗科技股份有限公司 磁共振系统的磁场分布检测方法、磁共振系统及成像方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JUN ZHANG,等: "Robust Single-Shot T2 Mapping via Multiple Overlapping-Echo Acquisition and Deep Neural Network", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 *
蔡淑惠,等: "时空编码核磁共振方法及其应用", 《厦门大学学报(自然科学版)》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113598716A (zh) * 2021-08-24 2021-11-05 南方医科大学 一种光声与磁共振联合成像方法
CN113598716B (zh) * 2021-08-24 2022-04-15 南方医科大学 一种光声与磁共振联合成像方法
CN115166612A (zh) * 2022-07-21 2022-10-11 厦门大学 基于深度学习的超快速磁共振成像高精度仿真方法及系统

Also Published As

Publication number Publication date
CN113030813B (zh) 2022-02-22

Similar Documents

Publication Publication Date Title
CN111656392B (zh) 用于合成磁共振图像的系统和方法
EP3534170A1 (en) System and method for sparse image reconstruction
CN108375746B (zh) 一种相位反卷折方法与设备
KR101627394B1 (ko) 자화율 가중 이미징 및 자화율 매핑(swim)을 이용한 핵 자기 공명 영상 생성 방법
CN110807492B (zh) 一种磁共振多参数同时定量成像方法及系统
CN113096208B (zh) 基于双域交替卷积的神经网络磁共振图像的重建方法
CN113030813B (zh) 一种磁共振t2定量成像方法及系统
CN113994225A (zh) 用于磁共振成像中的图像重建的系统和方法
JPH1170095A (ja) イメージングのための装置及び方法
Chen et al. Self‐calibrating wave‐encoded variable‐density single‐shot fast spin echo imaging
CN113298902B (zh) 卷褶视野磁共振图像的重建方法、计算机设备及存储介质
Pawar et al. A deep learning framework for transforming image reconstruction into pixel classification
CN110782452B (zh) 一种t2定量图像成像方法及系统
EP3953725A1 (en) Parameter map determination for time domain magnetic resonance
JP2002000581A (ja) 磁気共鳴画像化を使用して得られた複素画像の合成
KR102163337B1 (ko) 언더샘플링 방향 변화를 통한 다중영상획득 mri의 고속화 방법 및 이를 위한 장치
CN116725515A (zh) 一种磁共振快速成像方法
US11965945B2 (en) Magnetic resonance system and shimming method and imaging method thereof
CN113971706B (zh) 一种快速磁共振智能成像方法
Gutierrez et al. Reducing the complexity of model-based MRI reconstructions via sparsification
JP2023069890A (ja) 磁気共鳴イメージング装置、画像処理装置、及び、画像処理方法
CN115137347A (zh) 三维超短回波时间磁共振指纹成像的髓鞘定量成像方法
Shimron et al. CORE‐PI: Non‐iterative convolution‐based reconstruction for parallel MRI in the wavelet domain
US20200333414A1 (en) Method and system for generating synthethic images with switchable image contrasts
EP1176555A2 (en) Image processing method and apparatus, recording medium, and imaging apparatus

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
GR01 Patent grant
GR01 Patent grant