CN103142229B - 扩散峭度张量成像的高阶张量特征参数提取方法 - Google Patents
扩散峭度张量成像的高阶张量特征参数提取方法 Download PDFInfo
- Publication number
- CN103142229B CN103142229B CN201310056881.4A CN201310056881A CN103142229B CN 103142229 B CN103142229 B CN 103142229B CN 201310056881 A CN201310056881 A CN 201310056881A CN 103142229 B CN103142229 B CN 103142229B
- Authority
- CN
- China
- Prior art keywords
- tensor
- diffusion
- kurtosis
- eigenvalue
- sigma
- 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
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及医疗器械技术领域,为提出一种新的实用的人脑白质结构异常的磁共振成像参数提取方法,通过一定的特征提取方法对组织的各向异性进行分析,达到客观评价病变程度的目的,本发明采取的技术方案是,扩散峭度张量成像的高阶张量特征参数提取方法,包括如下步骤:通过受试者在磁共振扫描仪上采集组织沿多个方向的扩散加权信号,拟合得到反映组织内水分子扩散分布概率密度函数特征的二阶及四阶张量,通过一定的张量分析得到由于组织结构异常引起的水分子扩散分布异常的病变特征,主要有扩散系数和峭度系数两个不同角度的参数。本发明主要应用于医疗器械设计制造。
Description
技术领域
本发明涉及医疗器械技术领域,具体讲,涉及扩散峭度张量成像的高阶张量特征参数提取方法。
背景技术
阿尔兹海默症(Alzheimer’sdisease,AD)是老年痴呆症中最普遍最重要的一种,也是严重影响老年人健康的三大疾病之一。随着目前我国乃至直接全球社会老龄化程度的加剧,此类疾病给社会带来的经济负担和社会负担也随之大大的加重。此类疾病依然还没有能治愈的药物或者其它方法,而且该疾病也是有着发病周期长具有隐匿性,早期的诊断对于该疾病的早处理有着十分重要的意义。目前针对阿尔兹海默症的诊断研究开展较多,也有一定的研究进展,但是依然缺乏统一标准,国内外学者也积极采用了多种方法进行研究。
发明内容
本发明旨在克服现有技术的不足,提出一种新的实用的人脑白质结构异常的磁共振成像参数提取方法,通过一定的特征提取方法对组织的各向异性进行分析,达到客观评价病变程度的目的,为此,本发明采取的技术方案是,扩散峭度张量成像的高阶张量特征参数提取方法,包括如下步骤:通过受试者在磁共振扫描仪上采集组织沿多个方向的扩散加权信号,拟合得到反映组织内水分子扩散分布概率密度函数特征的二阶及四阶张量,通过一定的张量分析得到由于组织结构异常引起的水分子扩散分布异常的病变特征,主要有扩散系数和峭度系数两个不同角度的参数。
采集组织沿多个方向的扩散加权信号进一步细化为:
1.1数据的采集
在3.0T的磁共振扫面仪器上,采用单次激发平面回波(SE-EPI)序列进行30个敏感梯度方向3个扩散敏感因子b值(b=0,b1000和b2000s/mm2)的扩散加权DWI信号采集。
1.2数据的预处理
首先采用FSL(theFMRIBSoftwareLibrary)软件的FDT对数据进行涡流矫正和头动矫正操作,以及SUSAN对数据进行初步的去噪处理,其次在张量拟合再对数据进行非局部均值滤波操作。
采用非局部的均值滤波(NonLocalMeansFilter)方法,对每个像素点xi处理如下式(5-1):
这里,I代表待处理图像;v(xj)代表像素xj的灰度值;NL(v)(xi)代表像素v(xj)经过非局部均值滤波的结果;w(xi,xj)代表恢复v(xi)对应于v(xj)的权值;其中i和j均代表像素的索引编号。
这里,Z(i)是标准化的常数Z(i)=∑jw(i,j),θ是使用残差方法估计的噪声标准偏差,h是一个滤波参数,距离d表示为如下: 这里N=card(Ni)=card(Nj),yk和zk分别代表xi和yj领域Ni和Nj的第k个点,对于灰度级图像,Δ(v(yk),v(zk))=|(v(yk),v(zk))||。
张量的拟合方案:
S(b)是添加了梯度磁场的信号强度,即水分子在某个方向上的扩散情况,b值是扩散磁敏感因子;S(0)是未加梯度磁场的信号强度,及水分子的背景扩散情况;g是梯度敏感方向的单位向量,D是对应组织的扩散张量,是一个3×3的二维矩阵;基于量化水分子真实扩散分布偏离高斯扩散分布程度的峭度的扩散峭度成像(DKI,diffusionkurtosisimaging)有如下模型,如式(5-3)
这里,
其中,Dapp和Kapp分别为组织在梯度敏感方向g上的表观扩散系数和表观峭度系数;对应的张量即为D和W,分别对应一个二维三阶的实对称矩阵D和一个四维三阶的实对称矩阵W;MD是对应组织部位的平均扩散系数,由扩散张量矩阵D计算得来的;g依然是扩散敏感梯度方向,gi、gj、gk、gl是g的第i、j、k、l个分量元素,同上标T表示转置;
上述式(5-3)通过两个不同b值的信号S(b1)和S(b2),b1=1000s/mm2、b2=2000s/mm2以及一个b=0的信号;联立得到对应b1和b2的关于未知变量Dapp和Kapp的二元二次方程组如式(5-5),可解得对应30个方向上的Dapp和Kapp;
D的9个元素中仅有6个未知量,而W的81个元素中只有15个位置变量,分别对(5-4-1)和(5-4-2)展开,可得到如下式,
D11表示矩阵D的元素,其下标表示该元素在矩阵中的位置,与之类似W3321。通过(5-4)、(5-6)、(5-7)式采用非线性最小二乘算法即可得到最优的张量矩阵D和W,即对应的张量D和W。
特征提取方案:
张量的分解算法
首先是扩散张量矩阵D的分解,依然采用对角化分解方法,如下式,
这里,λi和vi是对应的特征值和特征向量(λ1≥λ2≥λ3),且3个vi两两正交,分别对应三个特征值的方向;其中作为特征值的λi代表组织内三分方向上的扩散系数大小,且最大特征值对应方向极为组织内神经纤维的主要走向,另两个特征值反映神经纤维髓鞘的特性;
其次对于高阶的峭度张量W,转化为解决此优化问题,
即得到如下方程组,
带入x=(x1,x2,x3)展开如下,
其解的情况如下,
如果W2111=W3111=0,即x2=x3=0,则λ=W1111,且对应特征向量为x=(±1,0,0);若x2≠0,x3=0,则另t=x1/x2,消元得到如下方程,
取其实根t,则对应的特征向量和特征值为,
若x3≠0,则另u=x1/x3,v=x2/x3,消元得到如下方程,
取其实根对(u,v),则特征向量和特征值为,
由此得到W的特征向量xi和特征值λi,1≤i≤n,且3≤n≤13代表特征值的个数。
特征参数计算
根据5.4.1的峭度张量特征分解结果给出峭度张量的提取特征,
平均峭度系数MK,meankurtosiscoefficient,
MK=mean(λi)
峭度各向异性KA,kurtosisanisotropy,
对每个i,将xi按着三个正交的vi′方向进行分组,即寻找其对应的最优i‘,
得到三组X,分别对每组Xi′向量对应的λXi′求平均,则
Ki′=mean(λXi′)(i′=1,2,3)
径向的峭度系数RK,radialkurtosiscoefficient,
RK=K1
轴向的峭度系数AK,axialkurtosiscoefficient,
AK=(K2+K3)/2。
本发明的技术特点及效果:
本发明提出了基于峭度信息的表征脑白质异变特征的参数提取方法。在采集MR信号并经过与处理之后,采用本发明所述的方法进行计算,可以得出包括扩散张量参数在内的以及峭度张量峭度特征参数的多种参数图像。利用这些参数,通过ROI统计分析、分类算法等可以获取反映白质退化程度的明显表征,对阿尔兹海默症早期的病变特征具有明显的更敏感的反映。如附图,及为各种参数反映阿尔兹海默症病人病变部位的图像。
附图说明
图1图像降噪示意图。
图2方法步骤流程图。
图3提取参数图像。
具体实施方式
本发明旨在提出一种新的实用的人脑白质结构异常的磁共振成像参数提取方法,通过多方向检测的DWI信号,拟合得到扩散位移分布的方差张量和峭度张量,通过一定的特征提取方法对组织的各向异性进行分析,达到客观评价病变程度的目的。
本发明提出了基于磁共振扩散峭度张量成像(MR-DKI,MR-diffusionkurtosistensorimaging)理论的反映组织结构特征的参数图像提取方法。其技术流程是:通过受试者在磁共振扫描仪上采集组织沿多个方向的扩散加权信号,拟合得到反映组织内水分子扩散分布概率密度函数特征的二阶及四阶张量,通过一定的张量分析得到由于组织结构异常引起的水分子扩散分布异常的病变特征,主要有扩散系数和峭度系数两个不同角度的参数。通过上述特征参数的分析,可以用于临床AD病人的诊断依据。
磁共振扫描仪可以对生物组织呈现多种模态的检测和信息获取,水分子的扩散特性是生物组织内的十分普遍的现象也能充分反映组织的围观结构特性。磁共振扩散加权成像(MR-DWI,diffusionweightedimaging)是检测活体组织内水分子扩散情况的技术,DKI技术就是通过多个方向上的DWI信号采集,获取组织内水分子的扩散分布情况,进而研究组织的微观结构特性。同时,本专利利用了DKI理论模型中的高阶峭度张量(KT,kurtosistensor),通过敏感的四阶张量,能更加敏感和精细的得到组织微观结构的异常,为阿尔兹海默症的进一步诊断提供有效的依据。
本发明方法步骤流程图如附图2。
1.1数据的采集方案
在3.0T的磁共振扫面仪器上,采用单次激发平面回波(SE-EPI)序列进行30个敏感梯度方向3个扩散敏感因子b值(b=0,b=1000和b=2000s/mm2)的扩散加权DWI信号采集。另外其它参数设置为:重复时间TR(repetitiontime)=10500ms,回波时间TE(ehcotime)=103ms,重复次数为Average(orNEXT)=1,数据获取时间TA(AcquisitionTime)=11’14’’,噪声水平为30,获取图像矩阵大小为128×128,视野FOV=230×230mm2,层厚为1.8mm,全脑采集73层,层间无间隔。因为对每个被试,将采集得到一个大小为128×128×73×61的矩阵数据集,采集时间为约11分钟。
1.2数据的预处理方案
首先采用英国牛津大学的FSL(theFMRIBSoftwareLibrary)软件的FDT对数据进行涡流矫正和头动矫正操作,以及SUSAN对数据进行初步的去噪处理。其次在张量拟合再对数据进行非局部均值滤波操作,本发明采用非局部的均值滤波方法,对每个像素点xi处理如下式(5-1):(见图1)
这里,I代表待处理图像;xi,xj均代表I中的像素点,xi为待处理像素点,0<i≤card(I),0<j≤card(I),且i≠j,i、j∈Z(Z代表整数,card(N)代表域N的像素点个数);v(xj)代表像素xj的灰度值;w(xi,xj)代表恢复v(xi)对应于v(xj)的权值。
这里,Z(i)是标准化的常数Z(i)=∑jw(i,j),是使用为残差方法估计的噪声标准偏差,h是一个滤波参数。距离d表示为如下: 这里Ni和Nj分别代表像素点xi和xj对应的领域;v(Ni),v(Nj)分别代表领域Ni和Nj的所有像素点的灰度值矩阵,N=card(Ni)=card(Nj);yk和zk分别代表xi和xj领域Ni和Nj的第k个点,则v(yk),v(zk)分别代表其对应灰度值;对于灰度级图像,Δ(v(yk),v(zk))=|(v(yk),v(zk))|。
1.3张量的拟合方案
扩散峭度成像(DKI,diffusionkurtosisimaging)技术是基于传统的扩散张量成像(DTI,diffusiontensorimaging)推广发展而来的,经典的扩散张量成像模型如式S(b)=S(0)e-b(xTDx),其中S(b)是添加了梯度磁场的信号强度(b=1000,2000s/mm2),即水分子在某个方向上的扩散情况,b值是扩散磁敏感因子;S(0)是未加梯度磁场的信号强度,及水分子的背景扩散情况;g是梯度敏感方向的单位向量,D是对应组织的扩散张量,是一个3×3的二维矩阵。但是一个张量D获取的组织结构信息对于微观病变特性很难识别,所以新提出的基于量化水分子真实扩散分布偏离高斯扩散分布程度的峭度的扩散峭度成像(DKI)技术有如下模型,如式(5-3)
这里,
其中,Dapp和Kapp分别为组织在某方向(这里指数据采集的梯度敏感方向g)上的表观扩散系数和表观峭度系数;Dij代表张量矩阵D的第i行第j列的元素;Wijkl代表张量矩阵W的第一维第i,第二维第j,第三维k和第四维第l的元素,对应的张量即为D和W,分别对应一个二维三阶的实对称矩阵D和一个四维三阶的实对称矩阵W;MD是对应组织部位的平均扩散系数,由扩散张量矩阵D计算得来的;g依然是扩散敏感梯度方向,gi、gj、gk、gl是g的第i、j、k、l个分量元素,同上标T表示转置;
上述式(5-3)通过两个不同b值的信号S(b1)和S(b2)(b1=1000s/mm2、b2=2000s/mm2)以及一个b0=0的信号,这里b0,b1,b2为b值的三个不同的取值;联立得到对应b1和b2的关于未知变量Dapp和Kapp的二元二次方程组如式(5-5),可解得对应30个方向上的Dapp和Kapp。
因为张量D和W,对应的矩阵是实对称的,D的9个元素中仅有6个未知量,而W的81个元素中只有15个位置变量。因此后分别对(5-4-1)和(5-4-2)展开,可得到如下式,
通过(5-4)、(5-6)、(5-7)式采用非线性最小二乘算法即可得到最优的张量矩阵D和W,即对应的张量D和W。
1.4特征提取方案
1.4.1张量的分解算法
首先是扩散张量矩阵D的分解,依然采用对角化分解方法,如下式,
这里,λi和vi是对应D的特征值和特征向量(λ1≥λ2≥λ3),vi是特征分解λi所对应的1×3的特征向量,这里i=1,2,3,且3个vi两两正交,分别对应三个特征值的方向;其中其中作为特征值的λi代表组织内三分方向上的扩散系数大小,且最大特征值对应方向极为组织内神经纤维的主要走向,另两个特征值反映神经纤维髓鞘的特性。
其次对于高阶的峭度张量W,我们转化为解决此优化问题,
即得到如下方程组,
带入x=(x1,x2,x3)展开如下,
这里诸如W1222均为张量矩阵的某一个独立元素值,即Wijkl各小标在(1,2,3)取值得到;上式其解的情况如下,
如果W2111=W3111=0,即x2=x3=0,则λ=W1111,且对应特征向量为x=(±1,0,0)。
若x2≠0,x3=0,则另t=x1/x2,消元得到如下方程,
取其实根t,则对应的特征向量和特征值为,
λ=Wx4
若x3≠0,则另u=x1/x3,v=x2/x3,消元得到如下方程,
取其实根对(u,v),则特征向量和特征值为,
由此得到W的特征向量xi和特征值λi,1≤i≤n,且3≤n≤13代表特征值的个数。
1.4.2特征参数计算
经典的扩散张量成像(DTI)算法的常用参数有各向异性FA(fractionalanisotropy),平均扩散系数MD(meandiffusivity),轴向扩散系数AD(axialdiffusivity)和径向扩散系数RD(radialdiffusivity)。本发明主要根据5.4.1的峭度张量特征分解结果给出峭度张量的提取特征,
平均峭度系数MK,meankurtosiscoefficient,
MK=mean(λi)
峭度各向异性KA,kurtosisanisotropy,
对每个i,将xi对应三个正交的vi′方向进行分组,即求夹角寻找其对应的最优i‘,这里i‘=1,2,3,代表扩散张量D的三个特征值或特征向量,参见1.4.1,
得到特征向量的三个分组X,分别对每组Xi’向量对应的λXi′求平均,则
Ki′=mean(λXi′)(i’=1,2,3)
径向的峭度系数RK,radialkurtosiscoeffocient,
RK=K1
轴向的峭度系数AK,axialkurtosiscoefficient,
AK=(K2+K3)/2
1.5特征的验证方案
分别对更具现有阿尔兹海默症最新诊断标准,获取阿尔兹海默症(AD)、轻度认知障碍(MCI)和控制组(ND)三组病人MR数据,采用本发明所述的方法进行特征提取,得到扩散信息参数和峭度信息参数若干。然后对阿尔兹海默症发病的几个显著部位(如颞叶、额叶等)的ROI区域或者体积进行统计病理分析等。
本发明针对阿尔兹海默症发病周期长,具有隐匿性,且诊断标准主观化的特点,提出了利用磁共振扩散结构成像技术进行客观诊断的新方法。基于峭度张量得到峭度信息相关参数,一次作为脑白质病变程度的输入参数,实现对患阿尔兹海默症高风险的老年人进行客观有效的诊断识别,具有重大的社会意义。
Claims (4)
1.一种扩散峭度张量成像的高阶张量特征参数提取方法,其特征是,包括如下步骤:通过受试者在磁共振扫描仪上采集组织沿多个方向的扩散加权信号,拟合得到反映组织内水分子扩散分布概率密度函数特征的二阶及四阶张量,通过一定的张量分析得到扩散系数和峭度系数两个不同角度的参数;对应张量矩阵的特征值和特征向量,采集组织沿多个方向的扩散加权信号进一步细化为:
1.1数据的采集
在3.0T的磁共振扫面仪器上,采用单次激发平面回波(SE-EPI)序列进行30个敏感梯度方向3个扩散敏感因子b值的扩散加权DWI信号采集;
1.2数据的预处理
首先采用FSL(theFMRIBSoftwareLibrary)软件的FDT对数据进行涡流矫正和头动矫正操作,以及SUSAN对数据进行初步的去噪处理,其次在张量拟合再对数据进行非局部均值滤波操作;
采用非局部的均值滤波(NonLocalMeansFilter)方法,对每个像素点Xp处理如下式(5-1):
这里,I代表待处理图像;v(Xp)代表像素Xp的灰度值;NL(v)(Xp)代表像素Xp的灰度值v(XP)经过非局部均值滤波的结果;w(Xp,Yq)代表恢复v(Xp)对应于v(Xq)的权值;其中p和q均代表像素的索引编号:
这里,Z(p)是标准化的常数Z(p)=∑qw(p,q),是使用残差方法估计的噪声标准偏差,h是一个滤波参数,距离d表示为如下:这里N=card(Np)=card(Nq),和分别代表像素Xp和像素Yq领域Np和Nq中的第个像素点,对于灰度级图像,
2.如权利要求1所述的扩散峭度张量成像的高阶张量特征参数提取方法,其特征是对应张量矩阵的特征值和特征向量,张量的拟合方案:
S(b)是添加了梯度磁场的信号强度,即水分子在某个方向上的扩散情况,b值是扩散磁敏感因子;S(0)是未加梯度磁场的信号强度,及水分子的背景扩散情况;g是梯度敏感方向的单位向量,D是对应组织的扩散张量,是一个3×3的二维矩阵;基于量化水分子真实扩散分布偏离高斯扩散分布程度的峭度的扩散峭度成像(DKI,diffusionkurtosisimaging)有如下模型,如式(5-3)
这里,
其中,Dapp和Kapp分别为组织在梯度敏感方向g上的表观扩散系数和表观峭度系数;对应的张量即为D和W,分别对应一个二维三阶的实对称矩阵D和一个四维三阶的实对称矩阵W;MD是对应组织部位的平均扩散系数,由扩散张量矩阵D计算得来的;g依然是扩散敏感梯度方向,gi、gj、gk、gl是向量g的第i、j、k、1个分量元素,i/j/k/l=1,2,3,上标T表示转置;
上述式(5-3)通过两个不同b值的信号S(b1)和S(b2),b1=1000s/mm2、b2=2000s/mm2以及一个b=0的信号;联立得到对应b1和b2的关于未知变量Dapp和Kapp的二元二次方程组如式(5-5),可解得对应30个方向上的Dapp和Kapp;
D的9个元素中仅有6个未知量,而W的81个元素中只有15个位置变量,分别对(5-4-1)和(5-4-2)展开,得到下式,
D11、D12、D13、D22、D23、D33表示矩阵D的元素,其中D12=D21、D23=D32、D13=D31,W3321与之类似;通过(5-4-1)、(5-4-2)、(5-6)、(5-7)式采用非线性最小二乘算法即可得到最优的张量矩阵D和W,即对应的张量D和W。
3.如权利要求2所述的扩散峭度张量成像的高阶张量特征参数提取方法,其特征是对应张量矩阵的特征值和特征向量,特征提取方案:
张量的分解算法
首先是扩散张量矩阵D的分解,依然采用对角化分解方法,如下式,
这里,和是对应D的特征值和特征向量,L1≥L2≥L3,且3个两两正交,分别对应三个特征值的方向;其中作为特征值的代表组织内三分方向上的扩散系数大小,且最大特征值对应方向极为组织内神经纤维的主要走向,另两个特征值反映神经纤维髓鞘的特性;
其次对于高阶的峭度张量W,转化为解决此优化问题,
即得到如下方程组,
带入x=(x1,x2,x3)展开如下,
其解的情况如下,
如果W2111=W3111=0,即x2=x3=0,则λ=W1111,且对应特征向量为x=(±1,0,0);若x2≠0,x3=0,则另t=x1/x2,消元得到如下方程,
取其实根t,则对应的特征向量和特征值为,
λ=Wx4
若x3≠0,则另u=x1/x3,v=x2/x3,消元得到如下方程,
取其实根对(u,v),则特征向量和特征值为,
λ=Wx4
由此得到W的特征向量xw和特征值λw,1≤w≤n,且3≤n≤13代表特征值的个数。
4.如权利要求3所述的扩散峭度张量成像的高阶张量特征参数提取方法,其特征是对应张量矩阵的特征值和特征向量,特征参数计算:
根据5-4-2的峭度张量特征分解结果给出峭度张量的提取特征,
平均峭度系数MK(meankurtosiscoefficient):
MK=mean(λw),
峭度各向异性KA(kurtosisanisotropy):
对每个w,将W的特征向量xw按着三个正交方向vw’进行分组,即寻找w其对应的最优w’,
得到三组向量V1、V2、V3,分别对每组Vw’向量对应的λVw’求平均,则
Kw‘=mean(λVw’)
径向的峭度系数RK(radialkurtosiscoefficient):RK=K1
轴向的峭度系数AK(axialkurtosiscoefficient):AK=(K2+K3)/2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310056881.4A CN103142229B (zh) | 2013-02-22 | 2013-02-22 | 扩散峭度张量成像的高阶张量特征参数提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310056881.4A CN103142229B (zh) | 2013-02-22 | 2013-02-22 | 扩散峭度张量成像的高阶张量特征参数提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103142229A CN103142229A (zh) | 2013-06-12 |
CN103142229B true CN103142229B (zh) | 2015-12-02 |
Family
ID=48540780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310056881.4A Active CN103142229B (zh) | 2013-02-22 | 2013-02-22 | 扩散峭度张量成像的高阶张量特征参数提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103142229B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11835611B2 (en) * | 2017-04-06 | 2023-12-05 | The United States Of America, As Represented By The Secretary, Department Of Health And Human Services | Isotropic generalized diffusion tensor MRI |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5946800B2 (ja) * | 2013-07-22 | 2016-07-06 | 株式会社日立製作所 | 磁気共鳴イメージング装置、画像処理装置、画像処理方法、及び画像処理プログラム |
CN104207776A (zh) * | 2014-08-22 | 2014-12-17 | 南昌大学 | 一种综合性磁共振成像装置及方法 |
CN105574849A (zh) * | 2015-11-25 | 2016-05-11 | 天津大学 | 一种基于扩散峭度张量的白质微结构特征可视化方法 |
CN105842642B (zh) * | 2016-03-17 | 2019-04-05 | 天津大学 | 基于峰度张量分数各向异性微结构特征提取方法与装置 |
CN107219483B (zh) * | 2017-04-22 | 2019-11-26 | 天津大学 | 一种基于扩散峰度成像的径向峰度各项异性定量方法 |
CN109256023B (zh) * | 2018-11-28 | 2020-11-24 | 中国科学院武汉物理与数学研究所 | 一种肺部气道微结构模型的测量方法 |
CN111351813B (zh) * | 2020-03-17 | 2021-09-24 | 无锡鸣石峻致医疗科技有限公司 | 一种基于非均匀场磁共振系统的表观扩散系数测量方法 |
CN114155225B (zh) * | 2021-12-07 | 2023-06-06 | 浙江大学 | 一种定量测量脑白质髓鞘内外水分子交换速率的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1885058A (zh) * | 2005-06-20 | 2006-12-27 | 西门子公司 | 借助磁共振确定扩散张量的系数的方法和装置 |
CN102202572A (zh) * | 2008-08-07 | 2011-09-28 | 纽约大学 | 用于提供实时扩散峭度成像的系统、方法和计算机可存取介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7995825B2 (en) * | 2001-04-05 | 2011-08-09 | Mayo Foundation For Medical Education | Histogram segmentation of FLAIR images |
-
2013
- 2013-02-22 CN CN201310056881.4A patent/CN103142229B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1885058A (zh) * | 2005-06-20 | 2006-12-27 | 西门子公司 | 借助磁共振确定扩散张量的系数的方法和装置 |
CN102202572A (zh) * | 2008-08-07 | 2011-09-28 | 纽约大学 | 用于提供实时扩散峭度成像的系统、方法和计算机可存取介质 |
Non-Patent Citations (3)
Title |
---|
Diffusion Kurtosis Imaging:The Quantification of Non-Gaussian Water Diffusion by Means of Magnetic Resonance Imaging;Jens H.Jensen,et al.,;《Magnetic Resonance in Medicine》;20051231;第53卷;第1432-1440页 * |
基于DKI的脑神经结构特性参数的优化提取研究;陈元园等;《2011年学术年会论文摘要》;20111231;第33页 * |
扩散峭度成像(DKI)在中枢神经系统的应用;曾丁巳等;《临床放射学杂志》;20111231;第30卷(第9期);第1400-1402页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11835611B2 (en) * | 2017-04-06 | 2023-12-05 | The United States Of America, As Represented By The Secretary, Department Of Health And Human Services | Isotropic generalized diffusion tensor MRI |
Also Published As
Publication number | Publication date |
---|---|
CN103142229A (zh) | 2013-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103142229B (zh) | 扩散峭度张量成像的高阶张量特征参数提取方法 | |
Schilling et al. | Challenges in diffusion MRI tractography–Lessons learned from international benchmark competitions | |
CN107658018B (zh) | 一种基于结构连接和功能连接的融合脑网络构建方法 | |
CN112529894B (zh) | 一种基于深度学习网络的甲状腺结节的诊断方法 | |
CN108288070B (zh) | 一种神经指纹提取分类方法及系统 | |
WO2024083057A1 (zh) | 基于多模态磁共振成像的图卷积神经网络疾病预测系统 | |
Assaf et al. | The CONNECT project: combining macro-and micro-structure | |
CN103886328B (zh) | 基于脑网络模块结构特征的功能磁共振影像数据分类方法 | |
Ahmed et al. | Single volume image generator and deep learning-based ASD classification | |
Nath et al. | Tractography reproducibility challenge with empirical data (TraCED): the 2017 ISMRM diffusion study group challenge | |
Kumar et al. | Fiberprint: A subject fingerprint based on sparse code pooling for white matter fiber analysis | |
CN105842642B (zh) | 基于峰度张量分数各向异性微结构特征提取方法与装置 | |
CN104267361A (zh) | 基于结构特征的自适应定量磁化率分布图复合重建的方法 | |
CN110232332A (zh) | 动态功能连接局部线性嵌入特征的提取及脑状态分类方法和系统 | |
CN104574298A (zh) | 一种基于互信息的多b值扩散权重图像的降噪方法 | |
CN109920550A (zh) | 一种基于dMRI研究青少年肌阵挛性癫痫的方法 | |
CN117172294B (zh) | 一种稀疏脑网络的构建方法、系统、设备和存储介质 | |
WO2024083058A1 (zh) | 一种大脑纤维束异常区域精准定位系统 | |
CN112305006A (zh) | 一种基于核磁共振信号的肿瘤组织标本的识别方法及系统 | |
CN107731283A (zh) | 一种基于多子空间建模的图像预选系统 | |
Karimi et al. | A machine learning-based method for estimating the number and orientations of major fascicles in diffusion-weighted magnetic resonance imaging | |
Yang et al. | Diagnosis of Parkinson’s disease based on 3D ResNet: The frontal lobe is crucial | |
CN114842969A (zh) | 一种基于关键纤维束的轻度认知障碍症评估方法 | |
CN104899884A (zh) | 一种用于早期帕金森症预测的综合分析方法 | |
CN112990266B (zh) | 多模态脑影像数据处理的方法、装置、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |