基于时间信号曲线的灌注动态影像的分割方法及系统与流程

本发明涉及一种灌注成像数据处理技术，具体涉及一种基于时间信号曲线的灌注动态影像的分割方法。本发明还涉及一种基于时间信号曲线的灌注动态影像的分割系统。

背景技术：

随着医学影像技术及计算机技术的发展，在医疗领域内以观察脑组织血流灌注状况为目的的动态脑灌注成像和定量分析技术开始应用于临床医学辅助诊断。通过对动态脑灌注检查所获得的一个时间影像序列的量化分析，可以测量出感兴趣区域脑组织的血流灌注量，从而绘制出对应的时间密度曲线。根据中心容积原理对这些曲线进行进一步分析、计算，可得到患者的相关脑血流动力学参数，如脑血流量（cbf）、脑血容量（cbv）、平均通过时间（mtt）、达峰时间（ttp）、最大剩余功能时间（tmax）等。按照一定的色彩分布对所得到的参数进行映像，就可以分别绘制出具有临床辅助诊断意义的功能影像。

目前，脑灌注成像后处理方法包括以下步骤：

第一步，寻找到动脉和静脉的位置，得到动脉流入曲线（arterialinputfunction，aif）和静脉流出曲线（veinoutputfunction，vof）；

第二步，求解反卷积以得到组织特性曲线；

第三步，由曲线的最大值、曲线下面积等进一步得到脑血流量（cbf）、脑血容量（cbv）、平均通过时间（mtt）、达峰时间（ttp）、最大剩余功能时间（tmax）等参数图像。

中国发明专利文献cn105701815a公开了一种磁共振灌注成像后处理方法及系统，通过对动脉输入函数和求解矩阵进行加权优化，将非线性问题转化为线性求解，加快了后处理速度。

这种后处理方法依赖于通过简化的生理模型将灌注数据降维约化为若干参数图，再从参数图像中通过目测、勾选区域测值、阈值处理等方式推测病变区域。但是，参数图像作为灌注数据的约化结果，实质丢失了大量有用信息，导致仅根据参数图像确定病变区域时往往不能得到稳定、有效、准确的结果。另外，参数图计算过程中，需要对大小血管的血红细胞含量等参数进行假设，还需要对动脉静脉进行识别以测定输入输出曲线，这些假设和识别均容易引入误差。此外，在对参数图像进行目测、勾选区域测值、阈值处理等二次处理过程中，还需要受过训练的医疗技术人员参与，过程耗时昂贵且容易产生人为误差。

技术实现要素：

本发明所要解决的技术问题是提供一种基于时间信号曲线的灌注动态影像的分割方法，它可以在完整地保留灌注数据全部信息的基础上，通过对灌注数据进行后处理，将灌注数据上的所有体素按照特征参数进行分类并标注，从而得到具有分类结果的三维影像数据。

为解决上述技术问题，本发明基于时间信号曲线的灌注动态影像的分割方法的技术解决方案为，包括以下步骤：

（1）、分类器训练步骤；

（1.1）、训练数据读入步骤；

读入灌注影像训练数据s；所述灌注影像训练数据s包含有脑组织所有体素集合v上的空间坐标的信息，以及采集时长t中的时间点t，标注为：

其中，s为灌注影像训练数据，

为空间坐标，

v为体素集合，

t为时间点，

t为采集时长；

所述灌注影像训练数据s为磁共振或计算机断层扫描的灌注成像所得到的训练数据，且训练数据经过预处理。

（1.2）、训练数据特征提取步骤；遍历灌注影像训练数据s的体素集合v，对所有空间坐标的时间信号曲线分别进行特征提取，使时间信号曲线以特征参数表达，得到一个特征提取函数；

所述特征提取的方法为高斯gaussian函数、gamma函数、梅尔倒谱系数mfcc、感知线性预测系数plp中的一种；

训练集生成步骤；将灌注影像训练数据s中的体素分为多个类别，对每一类别进行标注，从而生成训练集；

所述训练集生成步骤中的灌注影像训练数据s中的体素分为六类：0骨骼、1正常脑组织、2缺血脑组织、3梗死脑组织、4动静脉血管、5脑脊液；

（1.3）、分类器生成步骤；使用训练集训练分类器f，得到任意空间坐标体素的特征参数与其分类c的映射关系；

所述分类器f的训练采用随机森林或支持向量机。

（2）、分类器使用步骤；

（2.1）、待分析数据读入步骤；

读入灌注影像待分析数据，标注为：

其中，为灌注影像待分析数据；

所述灌注影像待分析数据为磁共振或计算机断层扫描的灌注成像所得到的待分析数据，且待分析数据经过预处理。

（2.2）、待分析数据特征提取步骤；遍历灌注影像待分析数据的体素集合v，对所有空间坐标的时间信号曲线使用特征提取函数g分别进行特征提取，得到；

（2.3）、分类标注步骤；使用分类器训练步骤所得到的分类器f，遍历灌注影像待分析数据体素集合v，对所有空间坐标的体素，按照其特征参数进行分类并标注，得到具有分类结果的三维影像数据。

所述步骤2.3中所得到的具有分类结果的三维影像数据中任一体素值为：

步骤2.4，图像后处理步骤，对具有分类结果的三维影像数据进行图像后处理，去除散落的奇异点。

步骤2.5，量化具有分类结果的三维影像数据中各类组织的体积。

本发明还提供一种基于时间信号曲线的灌注动态影像的分割系统，其技术解决方案为，包括分类器训练模块和分类器使用模块；

分类器训练模块包括训练数据读入单元、训练数据特征提取单元、训练集生成单元、分类器生成单元；

训练数据读入单元，被配置为读入灌注影像训练数据s；训练数据读入单元将灌注影像训练数据s标注为：

其中，s为灌注影像训练数据，

为空间坐标，

v为体素集合，

t为时间点，

t为采集时长；

训练数据特征提取单元，被配置为遍历灌注影像训练数据s的体素集合v，对所有空间坐标的时间信号曲线分别进行特征提取，使时间信号曲线以特征参数表达，得到一个特征提取函数；

训练集生成单元，被配置为将灌注影像训练数据s中的体素分为多个类别，对每一类别进行标注，从而生成训练集；

分类器生成单元，被配置为使用训练集训练分类器f，得到任意空间坐标体素的特征参数与其分类c的映射关系；

分类器使用模块包括待分析数据读入单元、待分析数据特征提取单元、分类标注单元；

待分析数据读入单元，被配置为读入灌注影像待分析数据，标注为：

其中，为灌注影像待分析数据；

待分析数据特征提取单元，被配置为遍历灌注影像待分析数据的体素集合v，对所有空间坐标的时间信号曲线使用特征提取函数g分别进行特征提取，得到；

分类标注单元，被配置为使用分类器训练步骤所得到的分类器f，遍历灌注影像待分析数据的体素集合v，对所有空间坐标的体素，按照其特征参数进行分类并标注，得到具有分类结果的三维影像数据。

在另一实施例中，所述分割系统还包括图像后处理单元，被配置为对具有分类结果的三维影像数据进行图像后处理，去除散落的奇异点。

在另一实施例中，所述分割系统还包括分割结果量化单元，被配置为量化具有分类结果的三维影像数据中各类组织的体积。

在另一实施例中，所述训练数据读入单元还包括数据预处理单元，数据预处理单元对磁共振或计算机断层扫描的灌注成像所得到的训练数据进行预处理，得到所述灌注影像训练数据s。

在另一实施例中，所述待分析数据读入单元还包括数据预处理单元，数据预处理单元对磁共振或计算机断层扫描的灌注成像所得到的待分析数据进行预处理，得到所述灌注影像待分析数据。

本发明可以达到的技术效果是：

本发明直接利用灌注数据中各体素的时间信号曲线的特征参数，对灌注影像数据中的所有体素进行分类标注，以实现对不同组织的分类，与现有技术基于时间信号曲线先计算出cbf、cbv、mtt、tmax等参数图，再进行分割的方法相比，能够完整地保留灌注数据全部信息。

本发明所得到的具有分类结果的三维影像数据，可以作为中间结果，供专业医疗技术人员参考。

本发明直接从所有体素的时间信号曲线上提取特征参数，并通过对各体素的特征参数进行分类标注，能够更好利用原始灌注数据的丰富信息对任一体素进行类别判断，避免了现有技术基于不完备生理模型的计算参数图过程造成的信息流失及误差引入。

本发明将计算量较大的训练任务提前完成，训练完毕的分类器投入使用后具有使用计算资源少，计算速度快等优点。本发明相较于传统方案参数图时的反卷积计算更加快捷。

本发明相较传统方案，还能够节省专业医疗技术人员对参数图像的勾画测量、阈值处理等步骤，医疗技术人员可以直接从本发明的分类结果与原始影像结合进行阅读。

附图说明

本领域的技术人员应理解，以下说明仅是示意性地说明本发明的原理，所述原理可按多种方式应用，以实现许多不同的可替代实施方式。这些说明仅用于示出本发明的教导内容的一般原理，不意味着限制在此所公开的发明构思。

结合在本说明书中并构成本说明书的一部分的附图示出了本发明的实施方式，并且与上文的总体说明和下列附图的详细说明一起用于解释本发明的原理。

下面结合附图和具体实施方式对本发明作进一步详细的说明：

图1是本发明基于时间信号曲线的灌注动态影像的分割方法的框架示意图；

如图2为采用本发明对灌注影像数据中的缺血梗死区域进行标注及体积计算的结果示意图。

如图3为采用本发明对灌注影像数据中的动脉进行标注的结果示意图；

图4为灌注影像数据中各组织的时间信号曲线图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。除非另外定义，此处使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本文中使用的“包括”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同，而不排除其他元件或者物件。

如图1所示，本发明基于时间信号曲线的灌注动态影像的分割方法，包括以下步骤：

1、分类器训练步骤；

1.1，训练数据读入步骤；

读入磁共振(mri)或计算机断层扫描（ct）的灌注成像所得到的训练数据，对训练数据进行预处理，得到经处理的灌注影像训练数据s；

预处理方法包括：逐时间点进行运动矫正、背景去除、图像去噪；

灌注影像训练数据s包含有脑组织所有体素集合v上的空间坐标的信息，以及采集时长t中的时间点t，标注为：

其中，s为灌注影像训练数据，

为空间坐标，

v为体素集合，

t为时间点，

t为采集时长；

1.2、训练数据特征提取步骤；遍历灌注影像训练数据s的体素集合v，对所有空间坐标的时间信号曲线分别进行特征提取，使时间信号曲线以特征参数表达，得到一个特征提取函数；

特征提取的方法可以采用现有技术中任意一种方法，如高斯gaussian函数、gamma函数、使用梅尔倒谱系数mfcc、感知线性预测系数plp等；实际训练时可迭代尝试多种特征提取方法以获得最佳表现；

训练集生成步骤；将灌注影像训练数据s中的体素分为多个类别，对每一类别进行标注，从而生成训练集；

例如：将缺血性卒中的灌注影像训练数据s中的体素根据其不同的特征参数分为六类：0骨骼、1正常脑组织、2缺血脑组织、3梗死脑组织、4动静脉血管、5脑脊液，此时；

当然，根据不同的应用场景，也可以将灌注影像训练数据分为其它类别；

1.3、分类器生成步骤；使用训练集训练分类器f，得到任意空间坐标体素的特征参数与其分类c的映射关系；

分类器训练时可以采用随机森林（randomforest），支持向量机（supportvectormachine）等方案，在实际训练时可迭代尝试不同方案以获得最佳表现。

2、分类器使用步骤；

2.1，待分析数据读入步骤；

读入待分析的磁共振(mri)或计算机断层扫描（ct）的灌注影像数据，对待分析数据进行预处理，得到经处理的灌注影像待分析数据，标注为：

其中，为灌注影像待分析数据，

为空间坐标，

v为体素集合，

t为时间点，

t为采集时长；

2.2、待分析数据特征提取步骤；遍历灌注影像待分析数据的体素集合v，对所有空间坐标的时间信号曲线使用特征提取函数g分别进行特征提取，得到灌注影像待分析数据中所有空间坐标的时间信号曲线的特征参数；

2.3、分类标注步骤；使用分类器训练步骤所得到的分类器f，遍历灌注影像待分析数据体素集合v，对所有空间坐标的体素，按照其特征参数进行分类并标注，得到一个具有分类结果的三维影像数据（如图2、图3所示），则该三维影像数据中任一体素值为：

2.4，图像后处理步骤；因同类组织均相链接，因此可利用组织的形态信息（如缺血组织应集中在一块区域）进行开闭操作及连接域操作，去除散落的奇异点；

2.5，得到分割结果；之后还可进一步地量化各类组织的体积，如计算缺血体积为：

缺血体积=分类为缺血脑组织2的体素个数×体素的体积。

由于磁共振(mri)和计算机断层扫描（ct）的灌注成像所得到的数据可以理解为一个三维的电影，灌注影像数据中的每一个体素（voxel）都可以描述为一条随时间变化的曲线（即如图4所示的时间信号曲线），造影剂在各组织中的代谢及通过速率的不同，各组织上体素的时间信号曲线也具有不同形态，因此不同组织上体素的时间信号曲线具有不同的特征参数。

本发明直接利用上述特性，根据不同的特征参数对不同组织上的体素进行分类标注，从而分辨出不同组织。本发明大大减少了数据处理量，加快了数据分析速度，从根本上解决了磁共振扫描速度慢的问题。

另外，本发明直接对灌注影像数据中的各体素进行分类标注，准确度更高。

本发明基于时间信号曲线的灌注动态影像的分割系统，包括分类器训练模块和分类器使用模块；

分类器训练模块包括训练数据读入单元、训练数据特征提取单元、训练集生成单元、分类器生成单元；

训练数据读入单元，被配置为读入磁共振(mri)或计算机断层扫描（ct）的灌注成像所得到的训练数据；进一步地，训练数据读入单元还包括数据预处理单元，数据预处理单元对所输入的训练数据进行预处理，得到经处理的灌注影像训练数据s；训练数据读入单元将灌注影像训练数据s标注为：

其中，s为灌注影像训练数据，

为空间坐标，

v为体素集合，

t为时间点，

t为采集时长；

训练数据特征提取单元，被配置为遍历灌注影像训练数据s的体素集合v，对所有空间坐标的时间信号曲线分别进行特征提取，使时间信号曲线以特征参数表达，得到一个特征提取函数；

训练集生成单元，被配置为将灌注影像训练数据s中的体素分为多个类别，对每一类别进行标注，从而生成训练集；

分类器生成单元，被配置为使用训练集训练分类器f，得到任意空间坐标体素的特征参数与其分类c的映射关系；

分类器使用模块包括待分析数据读入单元、待分析数据特征提取单元、分类标注单元、图像后处理单元、分割结果量化单元；

待分析数据读入单元，被配置为读入待分析的磁共振(mri)或计算机断层扫描（ct）的灌注影像数据，对待分析数据进行预处理，得到经处理的灌注影像待分析数据，标注为：

其中，为灌注影像待分析数据；

为空间坐标，

v为体素集合，

t为时间点，

t为采集时长；

待分析数据特征提取单元，被配置为遍历灌注影像待分析数据的体素集合v，对所有空间坐标的时间信号曲线使用特征提取函数g分别进行特征提取，得到；

分类标注单元，被配置为使用分类器训练步骤所得到的分类器f，遍历灌注影像待分析数据的体素集合v，对所有空间坐标的体素，按照其特征参数进行分类并标注，得到具有分类结果的三维影像数据，则该三维影像数据中任一体素值为：

该具有分类结果的三维影像数据能够直观地显示出各组织的分布，从而分辨出缺血组织、动静脉血管等；

图像后处理单元，被配置为利用组织的形态信息进行开闭操作及连接域操作，去除散落的奇异点；

分割结果量化单元，被配置为量化各类组织的体积。

显然，本领域的技术人员可以对本发明进行各种改动和变形，而不脱离本发明的精神和范围。这样，倘若本发明的这些修改属于本发明权利要求及其同等技术的范围之内，则本发明也意图包含这些改动和变形在内。