首页 / 专利库 / 诊断设备和程序 / 磁共振成像 / 弥散张量成像 / 融合概率型和确定型纤维束追踪的纤维束筛选方法

融合概率型和确定型纤维束追踪的纤维束筛选方法

阅读:401发布:2020-06-09

专利汇可以提供融合概率型和确定型纤维束追踪的纤维束筛选方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种融合概率型和确定型 纤维 束追踪的纤维束筛选方法,主要解决了 现有技术 所提取白质微结构特征真实性、准确性低的问题。其实现步骤是:1、读取DTI图像和FA图像,2、构建张量模板,3、提取张量模板的感兴趣纤维束,4、构建纤维束参数化模型,5、提取个体DTI图像的感兴趣纤维束,6、构建纤维束 频率 图谱,7、筛选感兴趣纤维束,本发明基于概率型纤维束追踪结果构建纤维束频率图谱,利用该图谱对纤维束参数化后的感兴趣纤维束进行筛选,使得纤维束参数化所提取的白质微结构特征精细化,有效克服了现有技术所提取白质微结构特征真实性,准确性低问题。,下面是融合概率型和确定型纤维束追踪的纤维束筛选方法专利的具体信息内容。

1.一种融合概率型和确定型纤维束追踪的纤维束筛选方法,其特征在于,构建纤维束频率图谱,筛选感兴趣纤维束,该方法包括如下步骤:
(1)读取弥散张量成像DTI图像和各向异性FA图像:
读入两组图像,一组为标准图像作为对照图像组,另一组为待测图像组,所述的对照图像组中包含格式均为.nii的44人的脑部弥散张量成像DTI图像和各向异性FA图像;所述的待测图像组中包含格式均为.nii的52人的脑部弥散张量成像DTI图像和各向异性FA图像,同一人的脑部弥散张量成像DTI图像和各向异性FA图像所处空间一致,组成该个体的个体空间;
(2)构建张量模板:
利用张量模板构建方法,构建张量模板,该张量模板所处的图像空间称为张量模板空间;
(3)提取张量模板的感兴趣纤维束:
利用追踪软件,对所构建张量模板进行确定型纤维束追踪,提取感兴趣纤维束;
(4)构建纤维束参数化模型:
利用纤维束参数化方法,构建纤维束参数化模型;
(5)提取个体弥散张量成像DTI图像的感兴趣纤维束:
(5a)利用追踪软件,对两组图像中每张弥散张量成像DTI图像上的感兴趣纤维束进行概率型纤维束追踪,得到纤维连接分布图;
(5b)将纤维连接分布图中所有不为零的纤维连接分布值,按照从大到小排列,选取排序中前千分之一的体素,将所选体素标记为1,纤维连接分布图中的其余体素标记为0,组成感兴趣纤维束图像;
(6)构建纤维束频率图谱:
(6a)将两组图像中每张各向异性FA图像配准到所构建的张量模板,得到其个体空间到张量模板空间的的配准参数;
(6b)利用配准参数,将每幅感兴趣纤维束图像映射到张量模板空间;
(6c)在张量模板空间下将感兴趣纤维束图像上的非零体素标记为1,得到标准化的感兴趣纤维束图像;
(6d)求所有标准化的感兴趣纤维束图像中对应位置上体素灰度值的平均值,将所有体素灰度值的平均值组成一张平均图像;
(6e)将平均图像上灰度值小于0.7的体素标记为0,其余体素上的灰度值不变得到所构建的纤维束频率图谱;
(7)筛选感兴趣纤维束:
(7a)寻找纤维束频率图谱中灰度值为零的点,去掉纤维束参数化后的感兴趣纤维束上对应坐标位置的点,得到筛选后的纤维束;
(7b)利用输入的两组图像中每张各向异性FA图像的个体空间到张量模板空间的配准参数,计算对应的逆配准参数;
(7c)利用对应的逆配准参数,将筛选后的纤维束上所有点的坐标,映射到每个个体空间;
(7d)利用三线性插值公式,插值计算出输入的两组图像中每张各向异性FA图像上对应个体空间下筛选后的纤维束上所有点的各向异性FA值;
(7e)在输入的两组图像中每张各向异性FA图像的个体空间下,求筛选后的纤维束中所有纤维中每个对应位置上点的各向异性FA值的平均值,将所有对应位置上点的各向异性FA值的平均值组成FA值平均值向量,输入的两组图像对应得到两组FA值平均值向量。
2.根据权利要求1所述的基于白质纤维束的白质微结构特征筛选方法,其特征在于,步骤(2)中所述张量模板构建方法的步骤如下:
第一步,利用图像配准软件,将对照图像组中格式为.nii的每张各向异性FA图像配准到格式为.nii的标准模板,得到配准参数;
第二步,在标准空间下的所有图像中选取一幅未选图像,将所选图像分别非线性配准到其余未选的每幅图像上,得到所选图像配准到其余的每幅图像上的配准参数,将所有配准参数进行平均,得到所选图像到其余图像上的平均配准参数;
第三步,将配准参数与平均配准参数融合,得到所选各向异性FA图像到标准空间下其余各向异性FA图像的变换参数;
第四步,将所选各向异性FA图像的变换参数作用到对应的弥散张量成像DTI图像上,得到标准化后的弥散张量成像DTI图像;
第五步,判断是否已选取完标准空间下的所有图像,若是,则执行第六步,否则,执行第二步;
第六步,求所有标准化后的弥散张量成像DTI图像中对应位置上体素灰度值的平均值;
第七步,将所有体素灰度的平均值组成一张平均图像,将该平均图像作为构建好的张量模板。
3.根据权利要求1所述的基于白质纤维束的白质微结构特征筛选方法,其特征在于,步骤(4)中所述纤维束参数化方法的步骤如下:
第一步,按照下式,计算感兴趣纤维束中每条纤维的纤维密度
其中,Dl表示感兴趣纤维束中第l条纤维的纤维密度,nl表示感兴趣纤维束中第l条纤维上点的总数,∑表示求和操作,i表示感兴趣纤维束中第l条纤维上点的序号,Mli表示感兴趣纤维束中第l条纤维上第i个点的纤维轨迹数目;
第二步,选取感兴趣纤维束中所有纤维的纤维密度中的最大值,作为原型纤维;
第三步,计算原型纤维上每个点处的切向量,对切向量做垂线,垂线与其他纤维相交得到原型纤维上每个点的匹配点的坐标。
4.根据权利要求1所述的基于白质纤维束的白质微结构特征筛选方法,其特征在于,步骤(7d)中所述三线性插值公式如下:
其中,f(xi,yi,zi)表示筛选后的纤维束上(xi,yi,zi)点的各向异性FA值,(xi,yi,zi)表示筛选后的纤维束上第i个点处的坐标, 分别表示沿x,y,z轴正方向上与(xi,yi,zi)点最近的整数坐标值, 分别表示沿x,y,z轴负方向上与(xi,yi,zi)点最近的整数坐标值;
分别 表示各向 异性FA图 像上
点的各向异性FA值。

说明书全文

融合概率型和确定型纤维束追踪的纤维束筛选方法

技术领域

[0001] 本发明属于物理技术领域,更进一步涉及影像处理技术领域中的一种融合概率型和确定型纤维束追踪的纤维束筛选方法。本发明利用所构建纤维束频率图谱,对经过纤维束参数化分析后的纤维束进行筛选,可以降低纤维束参数化分析所提取白质微结构特征的误差。

背景技术

[0002] 目前沿白质纤维束分析的主要思路之一是在张量模板图像上提取感兴趣纤维束,在对张量模板图像上的感兴趣纤维束进行定量化后,将该纤维束映射到个体脑部弥散张量成像DTI(Diffusion Tensor Imaging)图像上,实现不同个体间纤维束的对应。
[0003] 西安电子科技大学在其拥有的专利技术“一种基于健康人群的白质纤维参数化模型构建方法”(申请号:2014108358310,授权公告号:CN104537711 B)中公开了一种基于白质纤维束的纤维束参数化分析方法。该方法首先在张量模板上用确定型纤维束追踪方法提取出感兴趣纤维束,然后在感兴趣纤维束中选取一条代表性纤维,通过点匹配算法,在该代表性纤维上的点与其他纤维上的点之间建立对应关系。最后,将点匹配后的白质纤维束映射到个体图像中,提取纤维束所有点对应位置的白质扩散特性值,得到个体的白质微结构特征。该方法存在的不足之处在于,在张量模板上提取出感兴趣纤维束,所提取纤维束中的某些纤维,在个体图像上不存在对应纤维,导致所提取白质微结构特征真实性较低。并且,将点匹配后的纤维束映射到个体图像的过程中,可能由于映射误差导致白质纤维束映射到灰质区域,影响所提取白质微结构特征的准确性。

发明内容

[0004] 本发明的目的是针对上述现有方法的不足,提供了一种融合概率型和确定型纤维束追踪的纤维束筛选方法。本发明解决现有技术所提取白质微结构特征真实性、准确性较低的问题。
[0005] 实现本发明目的的思路是,对每个个体的DTI图像进行概率型纤维束追踪,基于概率性纤维束追踪结果构建纤维束频率图谱,利用该图谱对纤维束参数化后的感兴趣纤维束进行筛选,融合确定型和概率型纤维束追踪,使得纤维束参数化分析所提取的白质微结构特征真实可靠,并且提高了所提取白质微结构的准确性。
[0006] 本发明的方法包括如下步骤:
[0007] (1)读取DTI图像和FA图像:
[0008] 读入两组图像,一组为标准图像作为对照图像组,另一组为待测图像组,所述的对照图像组中包含格式均为 .nii的44人的脑部DTI图像和各向异性FA(FractionalAnisotropy)图像;所述的待测图像组中包含格式均为.nii的52人的脑部DTI图像和FA图像,同一人的脑部DTI图像和FA图像所处空间一致,组成该个体的个体空间;
[0009] (2)构建张量模板:
[0010] 利用张量模板构建方法,构建张量模板,该张量模板所处的图像空间称为张量模板空间;
[0011] (3)提取张量模板的感兴趣纤维束:
[0012] 利用追踪软件,对所构建张量模板进行确定型纤维束追踪,提取感兴趣纤维束;
[0013] (4)构建纤维束参数化模型:
[0014] 利用纤维束参数化方法,构建纤维束参数化模型;
[0015] (5)提取个体DTI图像的感兴趣纤维束:
[0016] (5a)利用追踪软件,对两组图像中每张DTI图像上的感兴趣纤维束进行概率型纤维束追踪,得到纤维连接分布图;
[0017] (5b)将纤维连接分布图中所有不为零的纤维连接分布值,按照从大到小排列,选取排序中前千分之一的体素,将所选体素标记为1,纤维连接分布图中的其余体素标记为0,组成感兴趣纤维束图像;
[0018] (6)构建纤维束频率图谱:
[0019] (6a)将两组图像中每张FA图像配准到所构建的张量模板,得到其个体空间到张量模板空间的的配准参数;
[0020] (6b)利用配准参数,将每幅感兴趣纤维束图像映射到张量模板空间;
[0021] (6c)在张量模板空间下将感兴趣纤维束图像上的非零体素标记为1,得到标准化的感兴趣纤维束图像;
[0022] (6d)求所有标准化的感兴趣纤维束图像中对应位置上体素灰度值的平均值,将所有体素灰度值的平均值组成一张平均图像;
[0023] (6e)将平均图像上灰度值小于0.7的体素标记为0,其余体素上的灰度值不变得到所构建的纤维束频率图谱;
[0024] (7)筛选感兴趣纤维束:
[0025] (7a)寻找纤维束频率图谱中灰度值为零的点,去掉纤维束参数化后的感兴趣纤维束上对应坐标位置的点,得到筛选后的纤维束;
[0026] (7b)利用输入的两组图像中每张FA图像的个体空间到张量模板空间的配准参数,计算对应的逆配准参数;
[0027] (7c)利用对应的逆配准参数,将筛选后的纤维束上所有点的坐标,映射到每个个体空间;
[0028] (7d)利用三线性插值公式,插值计算出输入的两组图像中每张FA图像上对应个体空间下筛选后的纤维束上所有点的FA值;
[0029] (7e)在输入的两组图像中每张FA图像的个体空间下,求筛选后的纤维束中所有纤维中每个对应位置上点的FA值的平均值,将所有对应位置上点的FA值的平均值组成FA值平均值向量,输入的两组图像对应得到两组FA值平均值向量;
[0030] 本发明与现有的技术相比具有以下优点:
[0031] 第一,由于本发明的方法采用构建纤维束频率图谱的方法,克服了现有技术所提取白质微结构特征真实性较低的不足,使得本发明所筛选后纤维束的白质微结构特征具有真实性。
[0032] 第二,由于本发明的方法采用构建纤维束参数化模型,筛选感兴趣纤维束的方法,克服了现有技术所提取白质微结构特征准确性低的不足,使得本发明提高了纤维束参数化分析所提取白质微结构特征的准确性。附图说明
[0033] 图1是本发明方法的流程图
[0034] 图2是本发明所构建的张量图谱;
[0035] 图3是本发明构建的纤维束频率图谱。

具体实施方式

[0036] 下面结合附图,对本发明作进一步的详细描述。
[0037] 参照附图1,对本发明的方法做进一步的详细描述。
[0038] 步骤1,读取DTI图像和FA图像。
[0039] 读入两组图像,一组为标准图像作为对照图像组,另一组为待测图像组,所述的对照图像组中包含格式均为.nii的44人的脑部DTI图像和FA图像;所述的待测图像组中包含格式均为.nii的52人的脑部DTI图像和FA图像,同一人的脑部DTI图像和FA图像所处空间一致,组成该个体的个体空间。
[0040] 步骤2,构建张量模板。
[0041] 利用FA图像和ICBM标准模板,构建得到张量图谱,得到的结果如图2所示。
[0042] 利用张量模板构建方法,构建张量模板,该张量模板所处的图像空间称为张量模板空间。
[0043] 所述张量模板构建方法的步骤如下:
[0044] 第1步,利用FSL软件,将对照图像组中格式为.nii的每张FA图像配准到格式为.nii的国际大脑图谱协会ICBM(International Consortium for Brain Mapping)标准模板,得到配准参数。
[0045] 第2步,在标准空间下的所有图像中选取一幅未选图像,将所选图像分别非线性配准到其余未选的每幅图像上,得到所选图像配准到其余的每幅图像上的配准参数,将所有配准参数进行平均,得到所选图像到其余图像上的平均配准参数。
[0046] 第3步,利用FSL软件,将配准参数与平均配准参数融合,得到所选FA图像到标准空间下其余FA图像的变换参数。
[0047] 第4步,利用FSL软件,将所选FA图像的变换参数作用到对应的DTI图像上,得到标准化后的DTI图像。
[0048] 第5步,判断是否已选取完标准空间下的所有图像,若是,则执行第6步,否则,执行第2步。
[0049] 第6步,求所有标准化后的DTI图像中对应位置上体素灰度值的平均值。
[0050] 第7步,将所有体素灰度的平均值组成一张平均图像,将该平均图像作为构建好的张量模板。
[0051] 步骤3,提取张量模板的感兴趣纤维束。
[0052] 利用exploreDTI软件,对所构建张量模板进行确定型纤维束追踪,提取感兴趣纤维束。
[0053] 步骤4,构建纤维束参数化模型。
[0054] 利用纤维束参数化方法,构建纤维束参数化模型。
[0055] 所述纤维束参数化方法的步骤如下:
[0056] 按照下式,计算感兴趣纤维束中每条纤维的纤维密度
[0057]
[0058] 其中,Dl表示感兴趣纤维束中第l条纤维的纤维密度,nl表示感兴趣纤维束中第l条纤维上点的总数,∑表示求和操作,i表示感兴趣纤维束中第l条纤维上点的序号,Mli表示感兴趣纤维束中第l条纤维上第i个点的纤维轨迹数目;
[0059] 选取感兴趣纤维束中所有纤维的纤维密度中的最大值,作为原型纤维;
[0060] 计算原型纤维上每个点处的切向量,对切向量做垂线,垂线与其他纤维相交得到原型纤维上每个点的匹配点的坐标。
[0061] 步骤5,提取个体DTI图像的感兴趣纤维束。
[0062] 利用FSL软件,对两组图像中每张DTI图像上的感兴趣纤维束进行概率型纤维束追踪,得到纤维连接分布图。
[0063] 将纤维连接分布图中所有不为零的纤维连接分布值,按照从大到小排列,选取排序中前千分之一的体素,将所选体素标记为1,纤维连接分布图中的其余体素标记为0,组成感兴趣纤维束图像。
[0064] 步骤6,构建纤维束频率图谱。
[0065] 利用所有标准化的感兴趣纤维束图像,构建纤维束频率图谱,得到的结果如图3所示。
[0066] 利用FSL软件,将两组图像中每张FA图像配准到所构建的张量模板,得到其个体空间到张量模板空间的的配准参数。
[0067] 利用配准参数,将每幅感兴趣纤维束图像映射到张量模板空间。
[0068] 在张量模板空间下将感兴趣纤维束图像上的非零体素标记为1,得到标准化的感兴趣纤维束图像。
[0069] 求所有标准化的感兴趣纤维束图像中对应位置上体素灰度值的平均值,将所有体素灰度值的平均值组成一张平均图像。
[0070] 将平均图像上灰度值小于0.7的体素标记为0,其余体素上的灰度值不变得到所构建的纤维束频率图谱。
[0071] 步骤7,筛选感兴趣纤维束。
[0072] 寻找纤维束频率图谱中灰度值为零的点,去掉纤维束参数化后的感兴趣纤维束上对应坐标位置的点,得到筛选后的纤维束。
[0073] 利用FSL软件,对输入的两组图像中每张FA图像的个体空间到张量模板空间的配准参数,计算对应的逆配准参数
[0074] 利用对应的逆配准参数,将筛选后的纤维束上所有点的坐标,映射到每个个体空间。
[0075] 利用三线性插值公式,插值计算出输入的两组图像中每张FA图像上对应个体空间下筛选后的纤维束上所有点的FA值。
[0076] 所述的三线性插值公式如下:
[0077]
[0078] 其中,f(xi,yi,zi)表示筛选后的纤维束上(xi,yi,zi)点的FA值,(xi,yi,zi)表示筛选后的纤维束上第i个点处的坐标, 分别表示沿x,y,z轴正方向上与(xi,yi,zi)点最近的整数坐标值, 分别表示沿x,y,z轴负方向上与(xi,yi,zi)点最近的整数坐标值;分别表示FA图像

点的FA值。
[0079] 在输入的两组图像中每张FA图像的个体空间下,求筛选后的纤维束中所有纤维中每个对应位置上点的FA值的平均值,将所有对应位置上点的FA值的平均值组成FA值平均值向量,输入的两组图像对应得到两组FA值平均值向量。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈