血流成像的制作方法

本发明涉及流体动态成像，且更特别地涉及基于流过血管或心脏结构的造影剂的动态成像的对象中血流特征的评估。相关技术的描述动态对比增强(dce)计算机断层扫描(ct)已被用于评估血管中的血流和流动压力，例如，如2019年5月16日提交的共同所有的第pct/ca2019/050668号国际pct申请(2019年11月21日公布为wo2019/218076)中所述。由dce ct得出的血流评估代表了成像扫描数秒过程中的平均值，通常大于10秒。因此，dce ct不具有足够的时间分辨率来实现4d流动成像，也就是说，以非常精细的时间分辨率跟踪所选图像体素中的血流特征，例如以小于1秒的时间间隔计算流速的变化。目前可用于非侵入性4d流动成像评估的技术是多普勒超声心动图和mri 4d血流(mri 4d flow)。多普勒超声心动图的一个潜在缺点是，如果超声束与血流射流没有很好地对准，则血流评估可能不太准确。mri 4d血流技术的一个潜在缺点是mri扫描仪的可及性有限，以及操作成本更高，且检查时间更长。最近描述的4d血流ct技术(lantz et al.(2018)intracardiac flow at4d ct:comparison with 4d flow mri；radiology,第289卷,第51-58页)使用基于navier-stokes方程的流体动力学的计算机模拟来模拟4d血流图，以预测主动脉中血流的大小和方向。使用流体动力学的计算机模拟来预测流动方向和大小存在缺点。使用计算机模拟的一个潜在缺点是，目前三维navier-stoke方程无法完全求解，即可能无法获得精确的解析解。因此，在一些情况下，精确模拟流动特性可能很困难。使用计算机模拟的另一个缺点是密集型计算机模拟是时间密集型的，并且需要大量的计算机处理，这不利于促进以精细的时间分辨率进行批量评估。这种计算机模拟需要超级计算架构来提供可以合理接近小于1秒的时间分辨率的评估，并且当在没有超级计算设施的情况下执行时，可能需要许多分钟且可能甚至几个小时来从图像采集的起点提供评估。因此，持续需要用于基于血流成像的对象血管评估的替代方法和系统。

背景技术：

技术实现思路

1、一方面，提供了用于血流成像的计算机实现的方法，包括：

2、获得包括捕获目标心血管系统中造影剂的增加阶段和下降阶段的至少一部分的多张对应图像的图像数据；

3、选择所述图像数据中的靶标体素；

4、基于所述图像数据，生成造影剂的时间增强曲线；

5、从所述时间增强曲线中选择第一和第二时间点；

6、通过减去选择的第一和第二时间点的图像数据确定所述增强的变化速率；

7、基于所述增强的变化速率确定靶标体素的血流特征。

8、其他方面，还提供了用于执行所述方法的系统和非瞬时计算机可读介质。

9、又一方面，提供了用于血流成像的系统，其包括：

10、用于存储图像数据的存储器，所述图像数据包括捕获目标心血管系统中造影剂的增加阶段和下降阶段的至少一部分的多张对应图像；

11、处理器，其被配置为基于所述图像数据生成靶标体素中造影剂的时间增强曲线；基于减去第一和第二时间点的图像数据，确定第一和第二时间点之间的增强的时间变化速率；以及基于所述增强的时间变化速率确定通过靶标体素的血流特征。

技术特征：

1.用于血流成像的计算机实现的方法，其包括：

2.如权利要求1所述的方法，其中所述时间增强曲线具有上坡或下坡。

3.如权利要求2所述的方法，其中所述时间增强曲线具有上坡和下坡。

4.如权利要求3所述的方法，其中所述图像数据包括在所述造影剂进入之前捕获所述目标心血管系统的至少一张图像。

5.如权利要求4所述的方法，还包括基于在所述造影剂进入之前捕获所述目标心血管系统的所述至少一张图像，确定所述时间增强曲线的基线值。

6.如权利要求5所述的方法，其中确定所述血流特征包括确定所述时间增强曲线下面积，以及将通过所述靶标体素的所述造影剂的持续时间确定为所述时间增强曲线中的信号强度大于所述基线值的时间间隔。

7.如权利要求6所述的方法，其中确定所述血流特征包括基于所述时间增强曲线下面积和所述持续时间，确定通过所述靶标体素的所述造影剂的密度。

8.如权利要求7所述的方法，其中确定所述血流特征包括使用雷诺输运定理确定绝对流速。

9.如权利要求8所述的方法，其中确定所述血流特征包括确定所述绝对流速的心脏引起的脉动。

10.如权利要求8所述的方法，其中确定所述血流特征包括基于多个时间间隔的绝对流速，确定所述目标心血管系统中的层流或紊流或层流-紊流-过渡血流。

11.如权利要求7所述的方法，其中确定所述血流特征包括使用雷诺输运定理确定流动加速度。

12.如权利要求11所述的方法，其中确定所述流动加速度包括获得所述图像数据期间的心电门控，并且其中所述第一时间点在心动周期的第一阶段，并且所述第二时间点在所述心动周期的第二阶段。

13.如权利要求8所述的方法，其中确定所述血流特征包括通过应用伯努利方程确定流动压力。

14.如权利要求1-5中任一项所述的方法，其中确定所述血流特征包括确定相对流速。

15.如权利要求14所述的方法，其中确定所述血流特征包括确定所述相对流速的心脏引起的脉动。

16.如权利要求14所述的方法，其中确定所述血流特征包括基于多个时间间隔的相对流速，确定所述目标心血管系统中的层流或紊流或层流-紊流-过渡血流。

17.如权利要求14所述的方法，其中确定所述血流特征包括确定所述造影剂的质量的时间变化速率。

18.如权利要求14所述的方法，其中确定所述血流特征包括通过应用伯努利方程确定流动压力。

19.如权利要求1-18中任一项所述的方法，其中所述多张对应图像为大于5张图像。

20.如权利要求1-18中任一项所述的方法，还包括从基于x射线的扫描或mri扫描，获取所述目标心血管系统的扫描数据，以及基于所述扫描数据重建图像数据。

21.如权利要求20所述的方法，其中从ct扫描获取所述扫描数据。

22.如权利要求20所述的方法，其中从mri扫描获取所述扫描数据。

23.如权利要求20所述的方法，其中从具有大于5秒的经过时间的扫描获取所述扫描数据。

24.如权利要求20-23中任一项所述的方法，还包括：

25.如权利要求1-24中任一项所述的方法，其中所述靶标体素是单一体素。

26.如权利要求1-24中任一项所述的方法，其中所述靶标体素是邻近体素的块。

27.如权利要求1-24中任一项所述的方法，其中所述靶标体素是各向同性的。

28.如权利要求1-24中任一项所述的方法，其中所述靶标体素是非各向同性的。

29.用于血流成像的系统，其包括：

30.如权利要求29所述的系统，其中所述时间增强曲线具有上坡或下坡。

31.如权利要求30所述的系统，其中所述时间增强曲线具有上坡和下坡。

32.如权利要求31所述的系统，其中所述图像数据包括在所述造影剂进入之前捕获所述目标心血管系统的至少一张图像。

33.如权利要求32所述的系统，其中所述处理器进一步被配置为基于在所述造影剂进入之前捕获所述目标心血管系统的所述至少一张图像，确定所述时间增强曲线的基线值。

34.如权利要求33所述的系统，其中确定所述血流特征包括确定所述时间增强曲线下面积，以及将通过所述靶标体素的所述造影剂的持续时间确定为所述时间增强曲线中的信号强度大于所述基线值的时间间隔。

35.如权利要求34所述的系统，其中确定所述血流特征包括基于所述时间增强曲线下面积和所述持续时间，确定通过所述靶标体素的所述造影剂的密度。

36.如权利要求35所述的系统，其中确定所述血流特征包括使用雷诺输运定理确定绝对流速。

37.如权利要求36所述的系统，其中确定所述血流特征包括确定所述绝对流速的心脏引起的脉动。

38.如权利要求36所述的系统，其中确定所述血流特征包括基于多个时间间隔的绝对流速，确定所述目标心血管系统中的层流或紊流或层流-紊流-过渡血流。

39.如权利要求35所述的系统，其中确定所述血流特征包括使用雷诺输运定理确定流动加速度。

40.如权利要求38所述的系统，其中确定所述流动加速度包括获得所述图像数据期间的心电门控，并且其中所述第一时间点在心动周期的第一阶段，且所述第二时间点在所述心动周期的第二阶段。

41.如权利要求36所述的系统，其中确定所述血流特征包括通过应用伯努利方程确定流动压力。

42.如权利要求29-33中任一项所述的系统，其中确定所述血流特征包括确定相对流速。

43.如权利要求42所述的系统，其中确定所述血流特征包括确定所述相对流速的心脏引起的脉动。

44.如权利要求42所述的系统，其中确定所述血流特征包括基于多个时间间隔的相对流速，确定所述目标心血管系统中的层流或紊流或层流-紊流-过渡血流。

45.如权利要求42所述的系统，其中确定所述血流特征包括确定所述造影剂的质量的时间变化速率。

46.如权利要求42所述的系统，其中确定所述血流特征包括通过应用伯努利方程确定流动压力。

47.如权利要求29-46中任一项所述的系统，其中所述多张对应图像为大于5张图像。

48.如权利要求29-46中任一项所述的系统，还包括x射线扫描仪或mri扫描仪，其被配置为从基于x射线的扫描或mri扫描，获取所述目标心血管系统的扫描数据，以及基于所述扫描数据重建图像数据。

49.如权利要求48所述的系统，其中从ct扫描获取所述扫描数据。

50.如权利要求48所述的系统，其中从mri扫描获取所述扫描数据。

51.如权利要求48所述的系统，其中从具有大于5秒的经过时间的扫描获取所述扫描数据。

52.如权利要求29-51中任一项所述的系统，其中所述靶标体素是单一体素。

53.如权利要求29-51中任一项所述的系统，其中所述靶标体素是邻近体素的块。

54.如权利要求29-51中任一项所述的系统，其中所述靶标体素是各向同性的。

55.如权利要求29-51中任一项所述的系统，其中所述靶标体素是非各向同性的。

56.用于血流成像的系统，其包括：

57.具现为用于执行血流成像方法的计算机可读代码的非瞬时计算机可读介质，其包括：

58.具现为用于执行权利要求1-28中任一项所述的方法的计算机可读代码的计算机可读介质。

技术总结
本文描述了用于血流成像的计算机实现的方法，包括：获得包括捕获目标心血管系统中造影剂的增加阶段和下降阶段的至少一部分的多张对应图像的图像数据；选择图像数据中的靶标体素；基于图像数据生成造影剂的时间增强曲线；从时间增强曲线中选择第一和第二时间点；通过减去选择的第一和第二时间点的图像数据确定增强的变化速率；基于增强的变化速率确定靶标体素的血流特征。还提供了用于执行所述方法的系统和非瞬时计算机可读介质。

技术研发人员：苏兆霖,李延炎
受保护的技术使用者：伦敦健康科学中心研究公司
技术研发日：
技术公布日：2024/1/14