基于脑血管CT影像的图像处理方法及图像处理装置与流程

本发明涉及一种灌注成像数据处理技术，具体涉及一种基于脑血管ct影像的图像处理方法。本发明还涉及一种基于脑血管ct影像的图像处理装置。

背景技术：

脑动脉的急性或慢性闭塞使得部分脑组织得不到充足的血液供应而坏死，从而形成缺血型卒中。在缺血型卒中的急性期影像学检查中，常使用计算机断层扫描动态灌注成像技术(computertomographyperfusionimaging，ctp)来计算cbf(cerebralbloodflow)、cbv(cerebralbloodvolume)、tmax(timetomaxenhancement)、ttp(timetopeak)等参数图，再进一步对这些参数图进行阈值处理，得到所谓的“影像学组织窗”，即核心梗死区(一般定义为cbf小于全脑中位数cbf的30％的区域)和缺血区(一般定义为tmax大于6秒的区域)，为医生提供诊疗决策的信息。

然而该方法依赖于覆盖全脑组织的ctp扫描，该扫描需要128排(256层)并配有摇篮床的ct扫描仪。由于此扫描仪价格昂贵，实践中并没有得到很好的普及。大部分医疗机构在急诊时，ct扫描还局限于传统的平扫像及增强像两种，无法开展ctp灌注成像，因而得不到组织窗参数。

目前也有一些基于人工智能技术的研究，希望仅使用ct平扫或ct增强像两种常规扫描即可得到梗死核心和/或缺血区。如brainomix公司的e-aspects软件，其使用了ct平扫影像进行梗死核心区的分割(参见nagels,jolyo,pfaffj,etal.e-aspectsderivedacuteischemicvolumesonnon-contrast-enhancedcomputedtomographyimages:[j].internationaljournalofstroke,2019)。omano等人发表了一篇使用cta影像使用3dcnn技术进行梗死核心区的分割(参见o.,t.,salli,e.,savolainen,s.,&kangasniemi,m.(2019).3dconvolutionalneuralnetworksappliedtoctangiographyinthedetectionofacuteischemicstroke.europeanradiologyexperimental,3(1).https://doi.org/10.1186/s41747-019-0085-6)。

但上述研究还处于理论阶段，并不能应用于实践。

技术实现要素：

本发明所要解决的技术问题是提供一种基于脑血管ct影像的图像处理方法，它可以利用常规ct扫描仪的扫描结果得到梗死核心区域和缺血区域。

为解决上述技术问题，本发明基于脑血管ct影像的图像处理方法的技术解决方案为，包括以下步骤：

步骤1，对深度神经网络进行训练；采用ctp灌注数据，根据ctp灌注数据计算出核心梗死和缺血区的标注参数图；

步骤2，采用常规的ct扫描仪进行三种方式的扫描，得到ncct、cta和cta+影像；

步骤2.1，进行ct平扫得到ncct影像数据，并重建为所述ctp灌注数据的层厚；

步骤2.2，注射造影剂，进行ct增强扫描得到cta影像数据，并重建为所述ctp灌注数据的层厚；

步骤2.3，等待3～8秒后，再进行一次ct扫描得到cta+影像数据；

步骤3，对ncct、cta和cta+影像进行自动fov的裁剪和配准，以使ncct、cta和cta+影像与训练深度神经网络时所输入的ctp影像的fov大致相当；并进行层方向上的重采样，以调整至步骤1神经网络训练中使用训练数据的层厚；

步骤4，对cta及cta+图像进行图像增强变换；

步骤5，将三种影像数据ncct、cta和cta+输入已训练好的深度神经网络中的三个通道中，神经网络的输出即为缺血区和梗死区的分割标注。

在另一实施例中，所述步骤1中对深度神经网络进行训练的方法是：

抽出ctp的三期时相；以第一期作为增强前的ct平扫像，以ctp动脉期的ctp期相作为cta增强像，以ctp动脉期后3～8秒的ctp期相作为cta+第二增强像；

将ctp的三期时相输入深度神经网络；

使深度神经网络根据输入的ctp灌注数据进行多次计算，输出训练结果。

在另一实施例中，所述训练结果为根据该ctp数据计算得到参数图标注的核心梗死区和缺血区。

在另一实施例中，所述输入的ctp灌注数据包括不同梗死大小和不同梗死位置的ctp灌注数据。

在另一实施例中，所述输入的ctp灌注数据的层厚与常规cta扫描的层厚保持一致。

在另一实施例中，所述输入的ctp灌注数据的层厚为1mm层厚。

在另一实施例中，在所述步骤1.2的cta扫描时，使扫描对象进入ct机架并同时进行螺旋扫描；在所述步骤1.3中的cta+扫描时，使扫描对象退出ct机架并同时进行螺旋扫描。

本发明还提供一种基于脑血管ct影像的图像处理装置，其技术解决方案为，包括：

深度神经网络训练模块，被配置为采用ctp灌注数据，根据ctp灌注数据计算出核心梗死和缺血区的标注参数图；

扫描影像获取模块，被配置为采用常规的ct扫描仪进行三种方式的扫描，得到ncct、cta和cta+影像；所述ncct影像为通过ct平扫所得到的ncct影像数据，并重建为所述ctp灌注数据的层厚；所述cta影像为注射造影剂后通过ct增强扫描所得到的cta影像数据，并重建为所述ctp灌注数据的层厚；所述cta+影像为等待3～8秒后再进行一次ct扫描所得到的cta+影像数据；

视场角裁剪和配准模块，被配置为对ncct、cta和cta+影像进行自动fov的裁剪和配准，使ncct、cta和cta+影像与训练深度神经网络时所输入的ctp影像的fov一致；

图像增强变换模块，被配置为对cta及cta+图像进行图像增强变换；

影像数据输入模块，被配置为将三种影像数据ncct、cta和cta+输入深度神经网络训练模块所训练好的深度神经网络中的三个通道中，深度神经网络的输出即为缺血区和梗死区的分割标注。

在另一实施例中，还包括重采样模块，被配置为对ncct、cta和cta+影像进行层方向上的重采样，以调整至与深度神经网络训练模块中使用训练数据的层厚。

在另一实施例中，所述深度神经网络训练模块包括：

ctp三期时相抽出模块，被配置为以第一期作为增强前的ct平扫像，以ctp动脉期的ctp期相作为cta增强像，以ctp动脉期后3～8秒的ctp期相作为cta+第二增强像；

输入模块，被配置为将ctp的三期时相输入深度神经网络；

训练模块，被配置为使深度神经网络根据所述输入模块输入的ctp灌注数据进行多次计算，输出训练结果。

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

本发明使用常规ct平扫、ct增强，以及第二期ct增强共三次扫描，能够计算得到所谓的影像学组织窗，即梗死核心区域和缺血区域。

本发明在cta扫描之后，间隔3～8秒后再进行一次cta+扫描，即多使用了一期cta+延迟扫描。通过增加该延迟扫描，能够得到准确的影像学组织窗参数，即梗死核心区域和缺血区域。

附图说明

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

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

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

图1是深度神经网络训练的输入及输出示意图；

图2是采用常规的ct扫描仪所得到的ncct、cta和cta+影像的示意图；

图3是缺血脑组织与正常脑组织随时间变化的图像信号。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。除非另外定义，此处使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本文中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。“包括”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同，而不排除其他元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电性的连接，不管是直接的还是间接的。“上”、“下”、“左”、“右”等仅用于表示相对位置关系，当被描述对象的绝对位置改变后，则该相对位置关系也可能相应地改变。

如图1所示，本发明基于脑血管ct影像的图像处理方法，包括以下步骤：

步骤1，对深度神经网络进行训练；

训练方法可以是采用ctp灌注数据，根据ctp灌注数据计算出核心梗死和缺血区标注；

具体地，首先，抽出ctp的三期时相；

以第一期作为增强前的ct平扫像(即ncct)，选取ctp动脉期即一段动脉中信号达到最大值时的ctp期相作为cta增强像(血管增强像)，以ctp动脉期后3～8秒的ctp期相作为cta+第二增强像；

然后，将ctp的三期时相输入深度神经网络；

最后，使深度神经网络根据输入的ctp灌注数据进行多次计算，输出训练结果；该训练结果为根据该ctp数据计算得到的核心梗死区和缺血区标注图，即ctp定义下的核心梗死区(一般定义为cbf小于全脑中位数cbf的30％的区域)和缺血区(一般定义为tmax大于6秒的区域)，如图1所示；

深度神经网络进行计算的过程中，所采用的ctp灌注数据作为训练数据，可以包括不同梗死大小和不同梗死位置的ctp灌注数据；所采用的ctp灌注数据的层厚与常规cta扫描的层厚保持一致，优选为1mm左右层厚的ctp灌注数据。

深度神经网络可以是三通道3dunet深度神经网络；

步骤2，采用常规的ct扫描仪进行三种方式的扫描，得到ncct、cta和cta+影像；

步骤2.1，进行ct平扫(即ncct)并重建为1mm左右层厚；

步骤2.2，注射造影剂，进行ct增强扫描(即cta扫描)并重建为1mm左右层厚；

步骤2.3，等待3～8秒后，再进行一次ct扫描(即cta+扫描)，扫描及重建参数与步骤2.2中所述cta扫描一致；

优选地，在步骤1.2的cta扫描时，使病人床进入ct机架并同时进行螺旋扫描；在步骤1.3中的cta+扫描时，使病人床退出ct机架并同时进行螺旋扫描；工作流程上，多增加的这次扫描带来的工作量比较少。

步骤3，对ncct、cta和cta+影像进行自动fov(thefieldofview，视场角)的裁剪和配准，以使ncct、cta和cta+影像与训练深度神经网络时所输入的ctp影像的fov大致相当；并进行层方向上的重采样，以调整至步骤1神经网络训练中使用训练数据的层厚；

步骤4，对cta及cta+图像进行图像增强变换；

由于ctp扫描与cta扫描所使用的造影剂浓度与注射速率不相同，因此需要对cta进行图像增强变换，以贴近训练深度神经网络所用ctp动脉期的效果；

步骤5，将三种扫描数据ncct、cta和cta+输入已训练好的深度神经网络中的三个通道中，神经网络的输出即为缺血区和梗死区的分割标注。

经大量研究发现，当造影剂进入大动脉后，造影剂进入脑组织的缺血部分和健康部分的延迟时间不同，二者之间的间隔时间为3～8秒，如图3所示；由于延迟时间不同，扫描结果也会有明显的灰度改变，因此在间隔3～8秒之后，脑组织中的缺血部分和健康部分才会因造影剂进入延迟不同而发生明显的灰度改变。

传统的ct扫描仪的扫描方式是仅使用ncct，或仅使用ncct和cta影像数据，无法体现出造影剂进入脑组织的缺血部分和健康部分的延迟时间的不同，在cta扫描时，造影剂很可能还没有进入脑组织或刚进入脑组织，此时缺血部分和健康部分的脑组织都还没有得到信号增强，因此无法得到影像学组织窗参数。

本发明利用了造影剂进入脑组织的缺血部分和健康部分的延迟时间不同的特性，在常规的cta扫描之后，间隔3～8秒后再进行一次cta+扫描，通过cta扫描与cta+延迟扫描的比对，能够直观地发现梗死核心区域和缺血区域。

值得注意的是，目前国内外对脑组织区域影像的前期研究中也采用多期cta扫描，即在cta动脉期扫描后，追加平衡期、静脉期等2～3个期相的扫描，以判断动脉的延迟充盈及静脉的引流情况。这种多期cta扫描一般用于评估各cta期相中血管的显影多少，再通过各种目测量表以对侧枝循环进行评估。本发明直接使用cta及cta+的图像的脑组织区域影像进行脑组织梗死或缺血的评估，而并不利用血管影像。因此本发明与现有的多期cta扫描之间具有显著的区别。

本发明无需依赖于覆盖全脑组织的ctp扫描，仅使用常规的ct扫描仪就能够得到影像学组织窗参数。

本发明适用于无全脑灌注扫描能力的ct扫描仪。

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