用于基于医学图像数据对心血管循环建模的方法和设备.pdf

上传人:t**** 文档编号:8076705 上传时间:2019-12-25 格式:PDF 页数:14 大小:447.98KB
返回 下载 相关 举报
摘要
申请专利号:

CN201110340679.5

申请日:

20110916

公开号:

CN102525443B

公开日:

20160511

当前法律状态:

有效性:

有效

法律详情:

IPC分类号:

A61B5/026,G06T7/20

主分类号:

A61B5/026,G06T7/20

申请人:

西门子公司

发明人:

P·沙马,B·乔治斯库,V·米哈列夫,T·陈,D·科马尼丘

地址:

德国慕尼黑

优先权:

61/383478,61/384382,13/226779

专利代理机构:

中国专利代理(香港)有限公司

代理人:

臧永杰;卢江

PDF下载: PDF下载
内容摘要

本发明公开了一种非侵入性患者特定的冠状动脉病评估的方法和系统。从医学图像数据产生冠状动脉的解剖学模型。基于医学图像数据中造影剂传播的时间空间表示来估计冠状动脉中血液的速度。在冠状动脉的解剖学模型中采用计算流体动力学(CFD)模拟来模拟血液流动,所述计算流体动力学采用冠状动脉中血液的所估计的速度作为边界条件。

权利要求书

1.一种用于基于医学图像数据对心血管循环建模的方法,包括:从医学图像数据中产生冠状动脉的解剖学模型;基于医学图像数据中造影剂传播的时间空间表示估计冠状动脉中血液的速度;和采用计算流体动力学(CFD)模拟来模拟冠状动脉的解剖学模型中的血液流动,该计算流体动力学模拟采用冠状动脉中血液的所估计的速度作为边界条件。 2.根据权利要求1所述的方法,其中医学图像数据包括四维医学图像数据,并且从医学图像数据产生冠状动脉的解剖学模型的步骤包括:从四维医学图像数据产生冠状动脉的四维解剖学模型。 3.根据权利要求2所述的方法,其中从四维医学图像数据产生冠状动脉的四维解剖学模型的步骤包括:在四维医学图像数据的多个帧中的每一个中分割感兴趣的冠状动脉;和在四维医学图像数据的多个帧中的每一个中产生所分割的感兴趣的冠状动脉的三维模型。 4.根据权利要求1所述的方法,其中基于医学图像数据中造影剂传播的时间空间表示估计冠状动脉中血液的速度的步骤包括:基于从医学图像数据提取的时间强度曲线来产生造影剂传播的时间空间表示;和估计造影剂传播的时间空间表示的斜率。 5.根据权利要求4所述的方法,其中基于从医学图像数据提取的时间强度曲线来产生造影剂传播的时间空间表示的步骤包括:在冠状动脉长度上随时间映射医学图像数据中的点的反向强度大小。 6.根据权利要求1所述的方法,其中采用计算流体动力学(CFD)模拟来模拟冠状动脉的解剖学模型中的血液流动的步骤,该计算流体动力学模拟采用冠状动脉中血液的所估计的速度作为边界条件,包括:解纳维叶-斯托克斯方程,其由具有所估计的血液速度作为入口速度边界条件的冠状动脉的解剖学模型所约束。 7.根据权利要求1所述的方法,其中采用计算流体动力学(CFD)模拟来模拟冠状动脉的解剖学模型中的血液流动的步骤,该计算流体动力学模拟采用冠状动脉中血液的所估计的速度作为边界条件,包括:在冠状动脉上随时间模拟压力和速度值。 8.根据权利要求1所述的方法,还包括:基于计算流体动力学(CFD)模拟计算一个或多个血液动力学参数。 9.根据权利要求8所述的方法,其中一个或多个血液动力学参数包括流动速率、压力降以及壁剪切应力中的至少一个。 10.根据权利要求8所述的方法,其中一个或多个血液动力学参数包括冠状动脉血流储备(CFR)以及血流储备分数(FFR)中的至少一个。 11.一种用于基于医学图像数据对心血管循环建模的设备,包括:用于从医学图像数据中产生冠状动脉的解剖学模型的装置;用于基于医学图像数据中造影剂传播的时间空间表示估计冠状动脉中血液的速度的装置;和用于采用计算流体动力学(CFD)模拟来模拟冠状动脉的解剖学模型中的血液流动的装置,该计算流体动力学模拟采用冠状动脉中血液的所估计的速度作为边界条件。 12.根据权利要求11所述的设备,其中用于基于医学图像数据中造影剂传播的时间空间表示估计冠状动脉中血液的速度的装置包括:用于基于从医学图像数据提取的时间强度曲线来产生造影剂传播的时间空间表示的装置;和用于估计造影剂传播的时间空间表示的斜率的装置。 13.根据权利要求12所述的设备,其中用于基于从医学图像数据提取的时间强度曲线来产生造影剂传播的时间空间表示的装置包括:用于在冠状动脉长度上随时间映射医学图像数据中的点的反向强度大小的装置。 14.根据权利要求11所述的设备,其中用于采用使用冠状动脉中血液的所估计的速度作为边界条件的计算流体动力学(CFD)模拟来模拟冠状动脉的解剖学模型中的血液流动的装置包括:用于解纳维叶-斯托克斯方程的装置,该方程由具有血液的所估计的速度作为入口速度边界条件的冠状动脉的解剖学模型所约束。 15.根据权利要求11所述的设备,还包括:用于基于计算流体动力学(CFD)模拟计算一个或多个血液动力学参数的装置。

说明书

本申请要求2010年9月16日提交的美国临时申请No.61/383,478和2010年9月20日提交的美国临时申请No.61/384,382的权益,其公开通过引用结合于此。

技术领域

本发明涉及采用医学图像来对心血管循环建模,并且更具体而言,涉及基于四维医学图像数据和数字模拟的非侵入性患者特定的冠状动脉病评估。

背景技术

心脏病是引起美国男人和女人死亡的主要原因,并且占不少于全球死亡的30%。尽管近年来医学进展已经在诊断和治疗复杂心脏病方面提供了重大的改进,但是早产儿发病率和死亡率(morality)的发生仍然巨大。对此的一个原因是缺乏对患者特定的参数的在体内和体外精确的估计,所述参数准确地表征了解剖学、生理学、以及血液动力学,所有这些在心血管疾病的发展中起重要作用。

基于技术(例如,计算机断层成像(CT)、血管造影术等)的医学成像通常在临床实践中用于表征冠状动脉中狭窄的严重程度。然而,这些技术只提供了解剖学评估,其对于临床决定做出经常是不充分的。特别是,冠状动脉狭窄的严重程度的解剖学评估经常导致高估或低估,高估或低估两者是不希望的。狭窄严重程度的高估可能导致不必要的介入和再狭窄的继发风险,而低估将有可能导致不治疗。精确功能性评估可能需要压力和/或流量的测量,所述测量以侵入性的方式来确定。

发明内容

本发明提供一种基于四维医学图像数据的非侵入性患者特定的冠状动脉病评估的方法和系统。具体而言,本发明的实施例采用非侵入性计算流体动力学(CFD(ComputationalFluidDynamics))提供流量和压力测量,其基于采用从四维医学图像数据得到的患者特定的边界条件的方法。本发明的实施例还提供基于四维医学图像数据测量患者的冠状动脉血流储备(CFR(coronaryflowreserve))的非侵入性方法。

在本发明的一个实施例中,冠状动脉的解剖学模型从医学图像数据产生。基于医学图像数据中造影剂传播的时间空间表示来估计冠状动脉中血液的速度。之后在冠状动脉的解剖学模型中采用计算流体动力学(CFD)模拟来模拟血液流动,该CFD模拟采用冠状动脉中血液的所估计的速度作为边界条件。

在本发明的另一个实施例中,接收在静止(rest)期间获取的第一医学图像数据序列和在充血期间获取的第二医学图像数据序列。从第一医学图像数据序列产生冠状动脉的第一解剖学模型,以及从第二医学图像数据序列产生冠状动脉的第二解剖学模型。基于第一医学图像数据序列中造影剂传播的时间空间表示来估计静止状态期间冠状动脉中血液的最大速度,和基于第一医学图像数据序列中造影剂传播的时间空间表示来估计冠状动脉中血液的最大充血速度。为确定血液流动速率,静止状态期间的最大速度映射为静止状态期间的平均速度,而充血时的最大速度映射为充血时的平均速度。基于静止流动速率和充血期间流动速率来计算冠状动脉的冠状动脉血流储备。

参照下面详细描述和附图,本发明的这些和其他优点对于本领域普通技术人员来说将变得显而易见。

附图说明

图1图示了根据本发明的实施例用于患者特定的冠状动脉病评估的方法;

图2图示了图1的方法的步骤的示例性结果;

图3图示了根据本发明的实施例基于造影剂传播的时间空间表示的速度估计;

图4图示了根据本发明的实施例基于医学图像数据非侵入性评估冠状动脉血流储备的方法;

图5图示了从血管中血液速度分布图确定血管中的流动速率;

图6是能够实施本发明的计算机的高级框图。

具体实施方式

本发明涉及非侵入性冠状动脉病评估,其采用从诸如计算机断层成像(CT)、磁共振成像(MRI)、以及超声波心动描记术数据的容积数据序列对心脏的患者特定的建模。这样的容积数据序列,这里也称为四维图像数据或四维图像,是在一段时期内获得以覆盖一个或多个心动周期的序列,其中每一帧均为三维图像(容积)。在此对本发明的实施例进行描述以便给出冠状动脉病评估方法的可视理解。数字图像经常由一个或多个对象(或形状)的数字表示来组成。对象的数字表示在这里经常在识别和操纵对象方面来描述。这样的操纵是在计算机系统的存储器或其他电路/硬件中完成的虚拟操纵。因此,要理解的是,本发明的实施例可在计算机系统内采用存储于计算机系统内的数据来执行。

基于用于冠心病功能性评估的分析的CFD技术通常基于冠状动脉的简化几何结构,具有得自人口范围数据的一般边界条件。这使得这样的技术不适于全面的患者特定的冠状动脉病评估,诸如在冠状动脉狭窄情况下的狭窄严重程度的评估。然而,在本发明的实施例中,基于CFD的非侵入性方法使用用于流量和几何结构(geometry)两者的患者特定的边界条件,其得自医学图像数据,例如高分辨CT数据。

本发明的实施例提供了用于非侵入性功能性冠状动脉病评估的方法和系统,其基于四维医学图像数据,诸如高分辨CT数据,结合采用计算流体动力学(CFD)建模和模拟的根本的患者特定的血液动力学分析。为了根本的血液动力学分析以产生要用于功能性评估的患者特定的参数,用于感兴趣的冠状动脉的四维(三维+时间)患者特定的几何模型从医学图像数据确定。基于通过造影剂传播的时间空间表示的造影剂传播分析的图像被执行来鲁棒地重新获得在感兴趣的冠状动脉上随时间的速度分布图(profile)。在通过从造影剂传播得出的速度分布图确定的入口边界条件下,在感兴趣的冠状动脉中执行患者特定的CFD模拟,以及从CFD模拟中得出血液动力学参数,以表征狭窄的程度。

图1图示了根据本发明的实施例用于患者特定的冠状动脉病评估的方法。图1的方法将表示患者冠状区域的图像数据转换为冠状动脉的患者特定的解剖学模型,并使用患者特定的冠状动脉模型来模拟冠状动脉心脏中血液流动。

参照图1,在步骤102,接收医学图像数据。具体而言,接收至少一个图像数据序列。图像数据序列可为在一定时期内获取的三维图像(容积)序列。例如,这样的四维图像数据(三维+时间)可在一个完整心搏周期上获取。一个或多个图像数据序列可采用各种医学成像模态来接收。例如,根据本发明的各种实施例,可以接收三维CT数据、二维动态血管造影术数据、和/或旋转血管造影术数据的序列,然而本发明并不限制于此。图像数据可直接从一个或多个图像获取装置接收,例如CT扫描仪或X射线装置。也可以例如从计算机系统的存储器或存储装置、或一些其他的计算机可读存储介质加载以前存储的图像数据。医学图像序列可以表示造影剂经过冠状动脉的传播。图2示例性地图示了图1的方法的步骤的结果。如图2中所示,图像202表示了高分辨CT容积的帧。

返回图1,在步骤104,从所接收的医学图像数据产生冠状动脉的患者特定的解剖学模型。患者特定的解剖学模型为使用四维医学图像数据产生的冠状动脉的四维(三维+时间)几何模型。为了产生冠状动脉的患者特定的解剖学模型,感兴趣的冠状动脉在四维图像数据的每一帧中被分割。感兴趣的冠状动脉采用任何冠状动脉分割方法在每一帧内被分割。例如,感兴趣的冠状动脉可以采用美国公布专利申请No.2010/0067760中描述的方法在CT容积内被分割,其中所述美国公布专利申请通过引用结合于此。之后产生几何表面模型用于每一帧中感兴趣的分割的冠状动脉。例如,用于冠状动脉的解剖学建模的方法在美国专利No.7,860,290和美国专利No.7/953,266中得以描述,两者均通过引用结合于此。这导致感兴趣的冠状动脉的解剖学模型,其表示冠状动脉随时间变化的解剖学。如图2中所示,图像204表示冠状动脉分割,图像206表示冠状动脉的解剖学建模,导致三维解剖学模型。

返回图1,在步骤106,基于所接收的医学图像数据中造影剂传播估计冠状动脉中血液流动的速度。具体而言,一旦产生了冠状树分割,从所接收的医学图像数据的时间强度曲线得出造影剂传播的时间空间表示。造影剂传播的时间空间表示被用于重新获得感兴趣的冠状动脉上随时间的速度分布图。

造影传播的时间空间表示通过结合强度分布图的可用测量能够鲁棒地估计速度分布图。图3图示了根据本发明的实施例基于造影剂传播的时间空间表示的速度估计。如图3中所示,图像300表示输入二维透视图像的时间强度曲线。图像300的时间强度曲线表示在血管长度上随时间二维透视图像内的点的反向强度大小(invertedintensitymagnitude)。该信息提供了造影剂传播的时间空间表示,如图310中所示。斜率312对于造影剂传播的时间空间表示被估计,以及该斜率312为血管中血液速度的一阶近似。之后该速度的近似可以被用于设定用于随后CFD建模的边界条件。血管中速度的这样的近似可对于所有的感兴趣的冠状动脉执行,导致感兴趣的冠状动脉的估计速度分布图。如图2中所示,图像208和210表示造影传播的时间空间分析,其用于产生时间下标的(time-indexed)速度分布图,导致用于CFD建模的速度边界条件。

返回图1,在步骤108,在冠状动脉中采用具有患者特定的边界条件的CFD来模拟血液流动。具体而言,CFD建模使用从造影剂传播所估计的冠状动脉流动速度分布图,来设立用于CFD建模的边界条件。血液被建模为牛顿流体,以及在刚性壁假定下通过用数字解离散化的纳维叶-斯托克斯(Navier-Stokes)方程(连续性和动量方程)来获得速度场。离散化的纳维叶-斯托克斯方程用来增量地模拟随时间在冠状动脉内的血液流动的速度和压力。冠状动脉的患者特定的解剖学也被输入给CFD建模,以便约束基于患者特定的解剖学的血液流动模拟。

几个血液动力学参数已经被建议用于冠状动脉病的功能性评估,诸如流动速率和压力降用于评估严重程度和狭窄,以及壁剪切应力用于斑块形成。然而,这样的参数是以前在得自人口范围数据的边界条件下基于冠状动脉的简化几何结构被计算的。根据本发明的有益实施例,医学图像数据,例如高分辨CT数据,不仅用于提供解剖学模型,而且用来估计患者特定的边界条件,用于经由CFD模拟提取这些血液动力学参数,以及采用这些参数用于冠状动脉病的功能性评估。如图2中所示,图像210表示基于CFD的血液动力学分析。

返回图1,在步骤110,输出患者特定的血液动力学参数。基于CFD模拟计算患者特定的血液动力学参数。具体而言,CFD模拟产生随时间经过冠状动脉的血液流动的速度和压力的模拟值。这些模拟值可用来计算各种血液动力学参数,它们可用来评估冠状动脉病。例如,流动速率和压力降可用于评估严重程度和狭窄,以及壁剪切应力可用于评估斑块形成。

用于冠状动脉狭窄功能性评估的另一类型的参数为血流储备参数,诸如冠状动脉血流储备(CFR)以及血流储备分数(FFR(fractionalflowreserve))。CFR定义为冠状动脉中最大充血流量与静止时同一冠状动脉中流量的比。FFR定义为狭窄血管中最大血流量与正常血管中最大血流量的比,并且用来表征狭窄的严重程度。在临床实践中,基于压力/流量的测量用来确定这些血流储备。因此,根据本发明的实施例,权利要求1的方法可用于估计这些血流储备。

根据本发明的实施例,医学图像数据,例如高分辨CT数据,可用于冠状动脉的解剖学建模,并用于基于造影剂传播确定最大速度。这随后被用于执行患者特定的CFD分析,其结果用于拟合程序中,该拟合程序将最大速度映射为平均速度,以便确定计算CFD值所必要的流动速率。

图4图示了根据本发明实施例的基于医学图像数据非侵入性评估冠状动脉血流储备的方法。如图4中所示,在步骤402,对于静止和充血接收医学图像数据。具体而言,接收在静止期间获取的第一图像数据序列和接收在最大充血流量条件期间获取的第二图像数据序列。这些图像数据序列可以为在一定时期内获取的三维图像(容积)序列。图像数据序列可采用各种医学成像模态来获取。例如,根据本发明的各种实施例,可以接收三维CT数据、二维动态血管造影术数据、和/或旋转血管造影术数据的序列,然而本发明并不限制于此。可以从一个或多个图像获取装置直接接收图像数据,例如CT扫描仪或X射线装置。也能够例如从计算机系统的存储器或存储装置、或一些其他的计算机可读存储介质加载以前存储的图像数据。医学图像的序列可以表示造影剂经过冠状动脉的传播。

在步骤404,对于所接收的图像序列的每一个接收感兴趣的冠状动脉的解剖学模型。应理解的是,步骤404可类似如上文所述的图1的步骤104被执行。在步骤406,基于造影剂传播对于静止图像数据和充血图像数据两者估计最大速度。具体而言,每个图像数据序列的最大速度可采用造影剂传播的时间空间表示来估计,如上文结合图1的步骤106所述。

在步骤406,采用患者特定的CFD模拟,将静止时的最大速度映射为平均静止速度,并且将充血时的最大速度映射为平均充血速度。图5图示了从血管中血液的速度分布图确定血管中的流动速率。如图5中所示,对于给定血管500中的血液流动速率的精确评估,在血管500的整个横截面S上需要速度信息u。通过在表面区域S上对流量分布图ρu(其中ρ为血液的密度)积分可以确定血液流动速率Q。因此,如图5中所示,质量流动速率Q也可被确定为Q=ρAU,其中U为平均血液速度以及A为横截面S的面积。非侵入性血液速度测量记录给定横截面上的最大速度Vmax,因此需要假定速度分布图来确定血液流动速率。在一些情况中,流动速率通过假定抛物线速度分布图从最大速度值来确定。根据该准则,平均速度被确定为0.5Vmax。然而,该假定并未考虑患者特定的几何结构和血液动力学。

代替于做关于血管速度分布图的假定,CFD模拟可用来获得将最大速度Vmax映射到平均速度Vavg的更实际的关系。对于小血管,包括冠状动脉,如下关系被用来将Vmax映射到Vavg:

V avg = f ( V max , W ) = V max ( 1 + pW q ) 2 , - - - ( 1 ) ]]>

其中W为沃姆斯莱(Womersley)数,定义为其中f是脉动流的频率,R为特征长度尺度(冠状动脉的横截面面积),ρ为密度,以及μ为血液的动态粘性。

参数p和q采用拟合程序来确定,所述拟合程序使用来自在有缺陷的边界条件下执行的一系列CFD模拟的数据,其中代替于规定速度分布图边界条件,只指定流值。具体而言,一系列CFD模拟被执行,其中不同流动速率被指定,以及诸如W、f、R、ρ和/或μ的参数的值是变化的。每个模拟都产生用于Vmax的值以及用于Vavg的值。基于用于每个模拟的值以及产生的用于Vmax和Vavg的值,映射参数p和q通过将这些参数拟合为由模拟产生的数据来确定。例如,这些参数可拟合为使用非线性最小二乘法拟合算法的数据。

一旦基于CFD模拟确定了映射参数p和q,步骤404中确定的最大速度Vmax便可以使用沃姆斯莱数和映射参数映射为平均速度Vavg,如方程(1)中所示。这对于静止图像数据和充血图像数据单独地被执行,导致平均静止血液速度和平局充血血液速度。

在步骤408,基于平均静止血液速度和平均充血血液速度计算CFR。具体而言,CFR可以被计算为:

CFR = Q hyp Q rest = ρA hyp V hyp ρA rest V rest - - - ( 2 ) ]]>

其中Vhyp和Vrest为得自的平均速度值,以及Vhyp和Vrest为冠状动脉分别在充血时和静止时的横截面面积。

在缺乏第二医学图像序列(即充血期间获取的图像)的情况下,类似上文图4中所呈现的方法,仍可以被应用。在这种情况下,充血血液流动分析可使用一组不同的出口边界条件(即对冠状动脉床(coronarybed)的流动特征建模的边界条件)来执行。由于冠状动脉中流动的多数阻力存在于微脉管循环中,充血条件通过适当地改变出口边界条件中的阻力/阻抗值来模拟。适于模拟充血条件的阻力/阻抗值的其他细节在S.Mantero等人的“TheCoronaryBedanditsRoleintheCardiovascularSystem:AReviewandanIntroductorySingle-BracnhModel”(JournalofBiomedicalEngineering,Volume14,Issue2,March1992,Pages109-116)中描述。

上文所述的用于非侵入性冠状动脉病评估的方法可在计算机上使用众所周知的计算机处理器、存储单元、存储装置、计算机软件以及其他的组件来实施。这种计算机的高级框图示于图6中。计算机602包含处理器604,其通过执行限定全部操作的计算机程序指令控制计算机602的这样的操作。计算机程序指令可以存储于存储装置612(例如,磁盘)中,并且在需要执行计算机程序指令时加载到存储器610中。因此,图1和4的方法的步骤可由存储于存储器6l0和/或存储装置612中的计算机程序指令定义,并由执行计算机程序指令的处理器604来控制。图像获取装置620,例如CT扫描装置、X射线获取装置等,可连接至计算机602来输入图像数据至计算机602。将图像获取装置620和计算机602作为一个装置来实施是可能的。图像获取装置620和计算机602通过网络无线通信也是可能的。计算机602还包括一个或多个网络接口606,用于与其他装置经由网络通信。计算机602还包括其他的输入/输出装置608,其能够使得用户与计算机602交互(例如,显示器、键盘、鼠标、扬声器、按钮等)。这种输入/输出装置608可与一组计算机程序结合作为注解工具(annotationtool)来使用用以注解从图像获取装置620接收的容积。本领域技术人员将意识到,实际的计算机的实施还能够包含其他组件,以及图6是用于说明性目的的这种计算机的一些组件的高级表示。

应将前述的详细描述理解为在每个方面说明性的和示例性的,但不是限制性的,并且这里所公开的本发明的范围不由详细描述所确定,而由根据专利法许可的完整宽度所解释的权利要求书所确定。应理解的是,这里示出和描述的实施例仅说明本发明的原理,并且本领域技术人员可在不偏离本发明的范围和精神的情况下实施各种修改。本领域技术人员能够在不偏离本发明的范围和精神的情况下执行各种其他特征组合。

用于基于医学图像数据对心血管循环建模的方法和设备.pdf_第1页
第1页 / 共14页
用于基于医学图像数据对心血管循环建模的方法和设备.pdf_第2页
第2页 / 共14页
用于基于医学图像数据对心血管循环建模的方法和设备.pdf_第3页
第3页 / 共14页
点击查看更多>>
资源描述

《用于基于医学图像数据对心血管循环建模的方法和设备.pdf》由会员分享,可在线阅读,更多相关《用于基于医学图像数据对心血管循环建模的方法和设备.pdf(14页珍藏版)》请在专利查询网上搜索。

1、(10)授权公告号 (45)授权公告日 (21)申请号 201110340679.5 (22)申请日 2011.09.16 61/383478 2010.09.16 US 61/384382 2010.09.20 US 13/226779 2011.09.07 US A61B 5/026(2006.01) G06T 7/20(2006.01) (73)专利权人 西门子公司 地址 德国慕尼黑 (72)发明人 P沙马 B乔治斯库 V米哈列夫 T陈 D科马尼丘 (74)专利代理机构 中国专利代理(香港)有限公 司 72001 代理人 臧永杰 卢江 US 2005/0020903 A1,2005.01。

2、.27, US 2006/0083415 A1,2006.04.20, US 2007/0248250 A1,2007.10.25, US 7693315 B2,2010.04.06, Taylor.CA.Patient-Specific Modeling of Cardiovascular Mechanics.Annu Rev Biomed Eng .2009, 第 11 卷 109-134. (54) 发明名称 用于基于医学图像数据对心血管循环建模的 方法和设备 (57) 摘要 本发明公开了一种非侵入性患者特定的冠状 动脉病评估的方法和系统。从医学图像数据产生 冠状动脉的解剖学模型。基于医。

3、学图像数据中造 影剂传播的时间空间表示来估计冠状动脉中血液 的速度。在冠状动脉的解剖学模型中采用计算流 体动力学 (CFD) 模拟来模拟血液流动, 所述计算 流体动力学采用冠状动脉中血液的所估计的速度 作为边界条件。 (30)优先权数据 (51)Int.Cl. (56)对比文件 审查员 王舒妍 (19)中华人民共和国国家知识产权局 (12)发明专利 权利要求书2页 说明书6页 附图5页 CN 102525443 B 2016.05.11 CN 102525443 B 1/2 页 2 1.一种用于基于医学图像数据对心血管循环建模的方法, 包括 : 从医学图像数据中产生冠状动脉的解剖学模型 ; 基。

4、于医学图像数据中造影剂传播的时间空间表示估计冠状动脉中血液的速度 ; 和 采用计算流体动力学 (CFD) 模拟来模拟冠状动脉的解剖学模型中的血液流动, 该计算 流体动力学模拟采用冠状动脉中血液的所估计的速度作为边界条件。 2.根据权利要求 1 所述的方法, 其中医学图像数据包括四维医学图像数据, 并且从医 学图像数据产生冠状动脉的解剖学模型的步骤包括 : 从四维医学图像数据产生冠状动脉的四维解剖学模型。 3.根据权利要求 2 所述的方法, 其中从四维医学图像数据产生冠状动脉的四维解剖学 模型的步骤包括 : 在四维医学图像数据的多个帧中的每一个中分割感兴趣的冠状动脉 ; 和 在四维医学图像数据的。

5、多个帧中的每一个中产生所分割的感兴趣的冠状动脉的三维 模型。 4.根据权利要求 1 所述的方法, 其中基于医学图像数据中造影剂传播的时间空间表示 估计冠状动脉中血液的速度的步骤包括 : 基于从医学图像数据提取的时间强度曲线来产生造影剂传播的时间空间表示 ; 和 估计造影剂传播的时间空间表示的斜率。 5.根据权利要求 4 所述的方法, 其中基于从医学图像数据提取的时间强度曲线来产生 造影剂传播的时间空间表示的步骤包括 : 在冠状动脉长度上随时间映射医学图像数据中的点的反向强度大小。 6.根据权利要求 1 所述的方法, 其中采用计算流体动力学 (CFD) 模拟来模拟冠状动脉 的解剖学模型中的血液流。

6、动的步骤, 该计算流体动力学模拟采用冠状动脉中血液的所估计 的速度作为边界条件, 包括 : 解纳维叶 - 斯托克斯方程, 其由具有所估计的血液速度作为入口速度边界条件的冠状 动脉的解剖学模型所约束。 7.根据权利要求 1 所述的方法, 其中采用计算流体动力学 (CFD) 模拟来模拟冠状动脉 的解剖学模型中的血液流动的步骤, 该计算流体动力学模拟采用冠状动脉中血液的所估计 的速度作为边界条件, 包括 : 在冠状动脉上随时间模拟压力和速度值。 8.根据权利要求 1 所述的方法, 还包括 : 基于计算流体动力学 (CFD) 模拟计算一个或多个血液动力学参数。 9.根据权利要求 8 所述的方法, 其中。

7、一个或多个血液动力学参数包括流动速率、 压力 降以及壁剪切应力中的至少一个。 10.根据权利要求 8 所述的方法, 其中一个或多个血液动力学参数包括冠状动脉血流 储备 (CFR) 以及血流储备分数 (FFR) 中的至少一个。 11.一种用于基于医学图像数据对心血管循环建模的设备, 包括 : 用于从医学图像数据中产生冠状动脉的解剖学模型的装置 ; 用于基于医学图像数据中造影剂传播的时间空间表示估计冠状动脉中血液的速度的 装置 ; 和 权 利 要 求 书 CN 102525443 B 2 2/2 页 3 用于采用计算流体动力学 (CFD) 模拟来模拟冠状动脉的解剖学模型中的血液流动的 装置, 该计。

8、算流体动力学模拟采用冠状动脉中血液的所估计的速度作为边界条件。 12.根据权利要求 11 所述的设备, 其中用于基于医学图像数据中造影剂传播的时间空 间表示估计冠状动脉中血液的速度的装置包括 : 用于基于从医学图像数据提取的时间强度曲线来产生造影剂传播的时间空间表示的 装置 ; 和 用于估计造影剂传播的时间空间表示的斜率的装置。 13.根据权利要求 12 所述的设备, 其中用于基于从医学图像数据提取的时间强度曲线 来产生造影剂传播的时间空间表示的装置包括 : 用于在冠状动脉长度上随时间映射医学图像数据中的点的反向强度大小的装置。 14.根据权利要求 11 所述的设备, 其中用于采用使用冠状动脉。

9、中血液的所估计的速度 作为边界条件的计算流体动力学 (CFD) 模拟来模拟冠状动脉的解剖学模型中的血液流动 的装置包括 : 用于解纳维叶 - 斯托克斯方程的装置, 该方程由具有血液的所估计的速度作为入口速 度边界条件的冠状动脉的解剖学模型所约束。 15.根据权利要求 11 所述的设备, 还包括 : 用于基于计算流体动力学 (CFD) 模拟计算一个或多个血液动力学参数的装置。 权 利 要 求 书 CN 102525443 B 3 1/6 页 4 用于基于医学图像数据对心血管循环建模的方法和设备 0001 本申请要求 2010 年 9 月 16 日提交的美国临时申请 No.61/383,478 和。

10、 2010 年 9 月 20 日提交的美国临时申请 No.61/384,382 的权益, 其公开通过引用结合于此。 技术领域 0002 本发明涉及采用医学图像来对心血管循环建模, 并且更具体而言, 涉及基于四维 医学图像数据和数字模拟的非侵入性患者特定的冠状动脉病评估。 背景技术 0003 心脏病是引起美国男人和女人死亡的主要原因, 并且占不少于全球死亡的 30。 尽管近年来医学进展已经在诊断和治疗复杂心脏病方面提供了重大的改进, 但是早产儿发 病率和死亡率 (morality) 的发生仍然巨大。对此的一个原因是缺乏对患者特定的参数的 在体内和体外精确的估计, 所述参数准确地表征了解剖学、 生。

11、理学、 以及血液动力学, 所有 这些在心血管疾病的发展中起重要作用。 0004 基于技术(例如, 计算机断层成像(CT)、 血管造影术等)的医学成像通常在临床实 践中用于表征冠状动脉中狭窄的严重程度。 然而, 这些技术只提供了解剖学评估, 其对于临 床决定做出经常是不充分的。特别是, 冠状动脉狭窄的严重程度的解剖学评估经常导致高 估或低估, 高估或低估两者是不希望的。狭窄严重程度的高估可能导致不必要的介入和再 狭窄的继发风险, 而低估将有可能导致不治疗。精确功能性评估可能需要压力和 / 或流量 的测量, 所述测量以侵入性的方式来确定。 发明内容 0005 本发明提供一种基于四维医学图像数据的非。

12、侵入性患者特定的冠状动脉 病评估的方法和系统。具体而言, 本发明的实施例采用非侵入性计算流体动力学 (CFD(Computational Fluid Dynamics) 提供流量和压力测量, 其基于采用从四维医学图 像数据得到的患者特定的边界条件的方法。本发明的实施例 还提供基于四维医学图像数 据测量患者的冠状动脉血流储备 (CFR(coronary flow reserve) 的非侵入性方法。 0006 在本发明的一个实施例中, 冠状动脉的解剖学模型从医学图像数据产生。基于医 学图像数据中造影剂传播的时间空间表示来估计冠状动脉中血液的速度。 之后在冠状动脉 的解剖学模型中采用计算流体动力学 。

13、(CFD) 模拟来模拟血液流动, 该 CFD 模拟采用冠状动 脉中血液的所估计的速度作为边界条件。 0007 在本发明的另一个实施例中, 接收在静止 (rest) 期间获取的第一医学图像数据 序列和在充血期间获取的第二医学图像数据序列。 从第一医学图像数据序列产生冠状动脉 的第一解剖学模型, 以及从第二医学图像数据序列产生冠状动脉的第二解剖学模型。基于 第一医学图像数据序列中造影剂传播的时间空间表示来估计静止状态期间冠状动脉中血 液的最大速度, 和基于第一医学图像数据序列中造影剂传播的时间空间表示来估计冠状动 脉中血液的最大充血速度。为确定血液流动速率, 静止状态期间的最大速度映射为静止状 说。

14、 明 书 CN 102525443 B 4 2/6 页 5 态期间的平均速度, 而充血时的最大速度映射为充血时的平均速度。基于静止流动速率和 充血期间流动速率来计算冠状动脉的冠状动脉血流储备。 0008 参照下面详细描述和附图, 本发明的这些和其他优点对于本领域普通技术人员来 说将变得显而易见。 附图说明 0009 图 1 图示了根据本发明的实施例用于患者特定的冠状动脉病评估的方法 ; 0010 图 2 图示了图 1 的方法的步骤的示例性结果 ; 0011 图 3 图示了根据本发明的实施例基于造影剂传播的时间空间表示的速度估计 ; 0012 图 4 图示了根据本发明的实施例基于医学图像数据非侵。

15、入性评估冠状动脉血流 储备的方法 ; 0013 图 5 图示了从血管中血液速度分布图确定血管中的流动速率 ; 0014 图 6 是能够实施本发明的计算机的高级框图。 具体实施方式 0015 本发明涉及非侵入性冠状动脉病评估, 其采用从诸如计算机断层成像 (CT)、 磁共 振成像 (MRI)、 以及超声波心动描记术数据的容积数据序列对 心脏的患者特定的建模。这 样的容积数据序列, 这里也称为四维图像数据或四维图像, 是在一段时期内获得以覆盖一 个或多个心动周期的序列, 其中每一帧均为三维图像 ( 容积 )。在此对本发明的实施例进 行描述以便给出冠状动脉病评估方法的可视理解。数字图像经常由一个或多。

16、个对象 ( 或形 状 ) 的数字表示来组成。对象的数字表示在这里经常在识别和操纵对象方面来描述。这样 的操纵是在计算机系统的存储器或其他电路 / 硬件中完成的虚拟操纵。因此, 要理解的是, 本发明的实施例可在计算机系统内采用存储于计算机系统内的数据来执行。 0016 基于用于冠心病功能性评估的分析的 CFD 技术通常基于冠状动脉的简化几何结 构, 具有得自人口范围数据的一般边界条件。这使得这样的技术不适于全面的患者特定的 冠状动脉病评估, 诸如在冠状动脉狭窄情况下的狭窄严重程度的评估。 然而, 在本发明的实 施例中, 基于CFD的非侵入性方法使用用于流量和几何结构(geometry)两者的患者。

17、特定的 边界条件, 其得自医学图像数据, 例如高分辨 CT 数据。 0017 本发明的实施例提供了用于非侵入性功能性冠状动脉病评估的方法和系统, 其基 于四维医学图像数据, 诸如高分辨 CT 数据, 结合采用计算流体动力学 (CFD) 建模和模拟的 根本的患者特定的血液动力学分析。 为了根本的血液动力学分析以产生要用于功能性评估 的患者特定的参数, 用于感兴趣的冠状动脉的四维 ( 三维 + 时间 ) 患者特定的几何模型从 医学图像数据确定。 基于通过造影剂传播的时间空间表示的造影剂传播分析的图像被执行 来鲁棒地重新获得在感兴趣的冠状动脉上随时间的速度分布图 (profile)。在通过从造影 剂。

18、传播得出的速度分布图确定的入口边界条件下, 在感兴趣的冠状动脉中执行患者特定的 CFD 模拟, 以及从 CFD 模拟中得出血液动力学参数, 以表征狭窄的程度。 0018 图 1 图示了根据本发明的实施例用于患者特定的冠状动脉病评估的方法。图 1 的 方法将表示患者冠状区域的图像数据转换为冠状动脉的患者特定的解剖学模型, 并使用患 者特定的冠状动脉模型来模拟冠状动脉心脏中血液流动。 说 明 书 CN 102525443 B 5 3/6 页 6 0019 参照图 1, 在步骤 102, 接收医学图像数据。具体而言, 接收至少一个图像数据序 列。图像数据序列可为在一定时期内获取的三维图像 ( 容积 。

19、) 序列。例如, 这样的四维图 像数据 ( 三维 + 时间 ) 可在一个完整心搏周期上获取。一个或多个图像数据序列可采用各 种医学成像模态来接收。例如, 根据本发明 的各种实施例, 可以接收三维 CT 数据、 二维动 态血管造影术数据、 和 / 或旋转血管造影术数据的序列, 然而本发明并不限制于此。图像数 据可直接从一个或多个图像获取装置接收, 例如 CT 扫描仪或 X 射线装置。也可以例如从计 算机系统的存储器或存储装置、 或一些其他的计算机可读存储介质加载以前存储的图像数 据。医学图像序列可以表示造影剂经过冠状动脉的传播。图 2 示例性地图示了图 1 的方法 的步骤的结果。如图 2 中所示。

20、, 图像 202 表示了高分辨 CT 容积的帧。 0020 返回图 1, 在步骤 104, 从所接收的医学图像数据产生冠状动脉的患者特定的解剖 学模型。患者特定的解剖学模型为使用四维医学图像数据产生的冠状动脉的四维 ( 三维 + 时间 ) 几何模型。为了产生冠状动脉的患者特定的解剖学模型, 感兴趣的冠状动脉在四维 图像数据的每一帧中被分割。 感兴趣的冠状动脉采用任何冠状动脉分割方法在每一帧内被 分割。例如, 感兴趣的冠状动脉可以采用美国公布专利申请 No.2010/0067760 中描述的方 法在 CT 容积内被分割, 其中所述美国公布专利申请通过引用结合于此。之后产生几何表面 模型用于每一帧。

21、中感兴趣的分割的冠状动脉。例如, 用于冠状动脉的解剖学建模的方法在 美国专利No.7,860,290和美国专利No.7/953,266中得以描述, 两者均通过引用结合于此。 这导致感兴趣的冠状动脉的解剖学模型, 其表示冠状动脉随时间变化的解剖学。如图 2 中 所示, 图像 204 表示冠状动脉分割, 图像 206 表示冠状动脉的解剖学建模, 导致三维解剖学 模型。 0021 返回图 1, 在步骤 106, 基于所接收的医学图像数据中造影剂传播估计冠状动脉中 血液流动的速度。 具体而言, 一旦产生了冠状树分割, 从所接收的医学图像数据的时间强度 曲线得出造影剂传播的时间空间表示。 造影剂传播的时。

22、间空间表示被用于重新获得感兴趣 的冠状动脉上随时间的速度分布图。 0022 造影传播的时间空间表示通过结合强度分布图的可用测量能够鲁棒地估计速度 分布图。图 3 图示了根据本发明的实施例基于造影剂传播的时间空间表示的速度估计。如 图 3 中所示, 图像 300 表示输入二维透视图像的时间强度曲线。图像 300 的时间强度曲 线表示在血管长度上随时间二维透视图像内的点的反向强度大小 (inverted intensity magnitude)。该信息提供了造影剂传播的时间空间表示, 如图 310 中所示。斜率 312 对于 造影剂传播的时间空间表示被估计, 以及该斜率 312 为血管中血液速度的。

23、一阶近似。之 后该速度的近似可以被 用于设定用于随后 CFD 建模的边界条件。血管中速度的这样的 近似可对于所有的感兴趣的冠状动脉执行, 导致感兴趣的冠状动脉的估计速度分布图。 如图 2 中所示, 图像 208 和 210 表示造影传播的时间空间分析, 其用于产生时间下标的 (time-indexed) 速度分布图, 导致用于 CFD 建模的速度边界条件。 0023 返回图 1, 在步骤 108, 在冠状动脉中采用具有患者特定的边界条件的 CFD 来模拟 血液流动。具体而言, CFD 建模使用从造影剂传播所估计的冠状动脉流动速度分布图, 来设 立用于 CFD 建模的边界条件。血液被建模为牛顿流。

24、体, 以及在刚性壁假定下通过用数字解 离散化的纳维叶 - 斯托克斯 (Navier-Stokes) 方程 ( 连续性和动量方程 ) 来获得速度场。 离散化的纳维叶 - 斯托克斯方程用来增量地模拟随时间在冠状动脉内的血液流动的速度 说 明 书 CN 102525443 B 6 4/6 页 7 和压力。冠状动脉的患者特定的解剖学也被输入给 CFD 建模, 以便约束基于患者特定的解 剖学的血液流动模拟。 0024 几个血液动力学参数已经被建议用于冠状动脉病的功能性评估, 诸如流动速率和 压力降用于评估严重程度和狭窄, 以及壁剪切应力用于斑块形成。 然而, 这样的参数是以前 在得自人口范围数据的边界条。

25、件下基于冠状动脉的简化几何结构被计算的。 根据本发明的 有益实施例, 医学图像数据, 例如高分辨 CT 数据, 不仅用于提供解剖学模型, 而且用来估计 患者特定的边界条件, 用于经由 CFD 模拟提取这些血液动力学参数, 以及采用这些参数用 于冠状动脉病的功能性评估。如图 2 中所示, 图像 210 表示基于 CFD 的血液动力学分析。 0025 返回图 1, 在步骤 110, 输出患者特定的血液动力学参数。基于 CFD 模拟计算患者 特定的血液动力学参数。具体而言, CFD 模拟产生随时间经过冠状动脉的血液流动的速度 和压力的模拟值。这些模拟值可用来计算各种血液动力学参数, 它们可用来评估冠。

26、状动脉 病。 例如, 流动速率和压力降可用于评估严重程度和狭窄, 以及壁剪切应力可用于评估斑块 形成。 0026 用于冠状动脉狭窄功能性评估的另一类型的参数为血流储备参数, 诸如冠状动脉 血流储备 (CFR) 以及血流储备分数 (FFR(fractional flow reserve)。CFR 定义为冠状动 脉中最大充血流量与静止时同一冠状动脉中流量的比。FFR 定义为狭窄血管中最大血流量 与正常血管中最大血流量的比, 并且用来表征狭窄的严重程度。在临床实践中, 基于压力 / 流量的测量用来确定这些血流储备。因此, 根据本发明的实施例, 权利要求 1 的方法可用于 估计这些血流储备。 0027。

27、 根据本发明的实施例, 医学图像数据, 例如高分辨 CT 数据, 可用于冠状动脉的解 剖学建模, 并用于基于造影剂传播确定最大速度。这随后被用于执行患者特定的 CFD 分析, 其结果用于拟合程序中, 该拟合程序将最大速度映射为平均速度, 以便确定计算 CFD 值所 必要的流动速率。 0028 图 4 图示了根据本发明实施例的基于医学图像数据非侵入性评估冠状动脉血流 储备的方法。如图 4 中所示, 在步骤 402, 对于静止和充血接收医学图像数据。具体而言, 接 收在静止期间获取的第一图像数据序列和接收在最大充血流量条件期间获取的第二图像 数据序列。这些图像数据序列可以为在一定时期内获取的三维图。

28、像 ( 容积 ) 序列。图像数 据序列可采用各种医学成像模态来获取。 例如, 根据本发明的各种实施例, 可以接收三维CT 数据、 二维动态血管造影术数据、 和 / 或旋转血管造影术数据的序列, 然而本发明并不限制 于此。可以从一个或多个图像获取装置直接接收图像数据, 例如 CT 扫描仪或 X 射线装置。 也能够例如从计算机系统的存储器或存储装置、 或一些其他的计算机可读存储介质加载以 前存储的图像数据。医学图像的序列可以表示造影剂经过冠状动脉的传播。 0029 在步骤 404, 对于所接收的图像序列的每一个接收感兴趣的冠状动脉的解剖学模 型。应理解的是, 步骤 404 可类似如上文所述的图 1。

29、 的步骤 104 被执行。在步骤 406, 基于 造影剂传播对于静止图像数据和充血图像数据两者估计最大速度。具体而言, 每个图像数 据序列的最大速度可采用造影剂传播的时间空间表示来估计, 如上文结合图 1 的步骤 106 所述。 0030 在步骤 406, 采用患者特定的 CFD 模拟, 将静止时的最大速度映射为平均静止速 度, 并且将充血时的最大速度映射为平均充血速度。图 5 图示了从血管中血液的速度分布 说 明 书 CN 102525443 B 7 5/6 页 8 图确定血管中的流动速率。如图 5 中所示, 对于给定血管 500 中的血液流动速率的精确 评估, 在血管 500 的整个横截面。

30、 S 上需要速度信息 u。通过在表面区域 S 上对流量分布图 u(其中为血液的密度)积分可以确定血液流动速率Q。 因此, 如图5中所示, 质量流动 速率 Q 也可被确定为 Q AU, 其中 U 为平均血液速度以及 A 为横截面 S 的面积。非侵入 性血液速度测量记录给定横截面上的最大速度 Vmax, 因此需要假定速度分布图来确定血液 流动速率。在一些情况中, 流动速率通过假定抛物线速度分布图从最大速度值来确定。根 据该准则, 平均速度被确定为 0.5Vmax。然而, 该假定并未考虑患 者特定的几何结构和血液 动力学。 0031 代替于做关于血管速度分布图的假定, CFD 模拟可用来获得将最大速。

31、度 Vmax映射 到平均速度 Vavg的更实际的关系。对于小血管, 包括冠状动脉, 如下关系被用来将 V max映射 到 Vavg: 0032 0033 其中 W 为沃姆斯莱 (Womersley) 数, 定义为 其中 f 是脉动流的频率, R 为特征长度尺度 ( 冠状动脉的横截面面积 ), 为密度, 以及 为血液的动态粘性。 0034 参数 p 和 q 采用拟合程序来确定, 所述拟合程序使用来自在有缺陷的边界条件下 执行的一系列 CFD 模拟的数据, 其中代替于规定速度分布图边界条件, 只指定流值。具体而 言, 一系列 CFD 模拟被执行, 其中不同流动速率被指定, 以及诸如 W、 f、 R。

32、、 和 / 或 的参 数的值是变化的。每个模拟都产生用于 Vmax的值以及用于 V avg的值。基于用于每个模拟的 值以及产生的用于Vmax和Vavg的值, 映射参数p和q通过将这些参数拟合为由模拟产生的数 据来确定。例如, 这些参数可拟合为使用非线性最小二乘法拟合算法的数据。 0035 一旦基于 CFD 模拟确定了映射参数 p 和 q, 步骤 404 中确定的最大速度 Vmax便可以 使用沃姆斯莱数和映射参数映射为平均速度Vavg, 如方程(1)中所示。 这对于静止图像数据 和充血图像数据单独地被执行, 导致平均静止血液速度和平局充血血液速度。 0036 在步骤 408, 基于平均静止血液速。

33、度和平均充血血液速度计算 CFR。具体而言, CFR 可以被计算为 : 0037 0038 其中 Vhyp和 V rest为得自 的平均速度值, 以及 Vhyp和 V rest为冠 状动脉分别在充血时和静止时的横截面面积。 0039 在缺乏第二医学图像序列 ( 即充血期间获取的图像 ) 的情况下, 类似上文图 4 中 所呈现的方法, 仍可以被应用。 在这种情况下, 充血血液流动分析可使用一组不同的出口边 界条件 ( 即对冠状动脉床 (coronary bed) 的流动特 征建模的边界条件 ) 来执行。由于冠 状动脉中流动的多数阻力存在于微脉管循环中, 充血条件通过适当地改变出口边界条件中 的阻。

34、力 / 阻抗值来模拟。适于模拟充血条件的阻力 / 阻抗值的其他细节在 S.Mantero 等人 的 “The Coronary Bed and its Role in the Cardiovascular System : A Review and an Introductory Single-Bracnh Model” (Journal of Biomedical Engineering, Volume 14, 说 明 书 CN 102525443 B 8 6/6 页 9 Issue 2, March 1992, Pages109-116) 中描述。 0040 上文所述的用于非侵入性冠状动脉。

35、病评估的方法可在计算机上使用众所周知的 计算机处理器、 存储单元、 存储装置、 计算机软件以及其他的组件来实施。这种计算机的高 级框图示于图6中。 计算机602包含处理器604, 其通过执行限定全部操作的计算机程序指 令控制计算机 602 的这样的操作。计算机程序指令可以存储于存储装置 612( 例如, 磁盘 ) 中, 并且在需要执行计算机程序指令时加载到存储器 610 中。因此, 图 1 和 4 的方法的步骤 可由存储于存储器 6l0 和 / 或存储装置 612 中的计算机程序指令定义, 并由执行计算机程 序指令的处理器 604 来控制。图像获取装置 620, 例如 CT 扫描装置、 X 射。

36、线获取装置等, 可 连接至计算机 602 来输入图像数据至计算机 602。将图像获取装置 620 和计算机 602 作为 一个装置来实施是可能的。 图像获取装置620和计算机602通过网络无线通信也是可能的。 计算机602还包括一个或多个网络接口606, 用于与其他装置经由网络通信。 计算机602还 包括其他的输入 / 输出装置 608, 其能够使得用户与计算机 602 交互 ( 例如, 显示器、 键盘、 鼠标、 扬声器、 按钮等 )。这种输入 / 输出装置 608 可与一组计算机程序结合作为注解工具 (annotation tool) 来使用用以注解从图像获取装置 620 接收的容积。本领域。

37、技术人员将 意识到, 实际的计算机的实施还能够包含其他组件, 以及图 6 是用于说明性目的的这种计 算机的一些组件的高级表示。 0041 应将前述的详细描述理解为在每个方面说明性的和示例性的, 但不是限制性的, 并且这里所公开的本发明的范围不由详细描述所确定, 而由根据专利法许可的完整宽度所 解释的权利要求书所确定。 应理解的是, 这里示出和描述的实施例仅说明本发明的原理, 并 且本领域技术人员可在不偏离本发明的范围和精神的情况下实施各种修改。 本领域技术人 员能够在不偏离本发明的范围和精神的情况下执行各种其他特征组合。 说 明 书 CN 102525443 B 9 1/5 页 10 图 1 说 明 书 附 图 CN 102525443 B 10 2/5 页 11 图 2 说 明 书 附 图 CN 102525443 B 11 3/5 页 12 图 3A 图 3B 说 明 书 附 图 CN 102525443 B 12 4/5 页 13 图 4 图 5 说 明 书 附 图 CN 102525443 B 13 5/5 页 14 图 6 说 明 书 附 图 CN 102525443 B 14 。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 >


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1