用于患者特异性血流建模的方法和系统与流程

分案说明

本申请是申请日为2011年7月29日、国家申请号为201180048878.7、发明名称为“用于患者特异性血流建模的方法和系统”的发明专利申请的分案申请。

优先权

本申请要求2010年8月12日提交的美国临时申请号61/401,462、2010年8月20日提交的美国临时申请号61/401,915、2010年8月26日提交的美国临时申请号61/402,308、2010年8月27日提交的美国临时申请号61/402,345，以及2010年10月1日提交的美国临时申请号61/404,429的优先权益，所述美国临时申请以引用的方式整体并入本文中。

实施方案包括用于流体流量建模的方法和系统，以及更具体来说用于患者特异性血流建模的方法和系统。

背景技术：

冠状动脉病可在向心脏提供血液的血管中产生冠状动脉病变，如狭窄(血管异常变窄)。因此，流向心脏的血流可能受到限制。罹患冠状动脉病的患者可在身体用力期间经历胸痛(称为慢性稳定绞痛)或当患者在休息时经历不稳定绞痛。较严重的疾病表现可引起心肌梗塞或心脏病发作。

需要提供与冠状动脉病变有关的较精确数据，例如尺寸、形状、位置、功能意义(例如病变是否影响血流)等。可对罹患胸痛和/或显示冠状动脉病症状的患者进行一种或多种测试，所述测试可提供一些有关冠状动脉病变的间接证据。举例来说，非侵入性测试可包括心电图、由血液测试进行生物标记评价、平板测试(treadmilltest)、心回波描记术(echocardiography)、单正电子发射计算机断层摄影术(spect)，和正电子发射断层摄影术(pet)。然而，这些非侵入性测试通常不提供冠状动脉病变的直接评估或评估血液流动速率。非侵入性测试可通过寻找以下的变化来提供冠状动脉病变的间接证据：心脏电活动(例如使用心电描记术(ecg))、心肌运动(例如使用应力心回波描记术)、心肌灌注(例如使用pet或spect)，或代谢变化(例如使用生物标记)。

举例来说，解剖学数据可使用冠状动脉计算断层摄影血管造影术(ccta)以非侵入方式获得。ccta可用于对患有胸痛的患者进行成像，并且涉及使用计算机断层摄影术(ct)技术以在静脉内输注对比剂之后对心脏和冠状动脉进行成像。然而，ccta也不能提供关于冠状动脉病变的功能意义(例如病变是否影响血流)的直接信息。另外，因为ccta纯粹是一种诊断性测试，所以它不能用以预测以下的变化：在其它生理学状态(例如运动)下的冠状动脉血流、压力或心肌灌注，它也不能用以预测介入的结果。

因此，患者还可需要侵入性测试，如诊断性心导管插入术以使冠状动脉病变可视化。诊断性心导管插入术可包括进行常规冠状动脉血管造影术(cca)以收集关于冠状动脉病变的解剖学数据，所述解剖学数据是由医生用动脉的尺寸和形状图像来提供。然而，cca不提供用于评估冠状动脉病变的功能意义的数据。举例来说，医生在没有确定病变是否具有功能意义的情况下不能诊断冠状动脉病变是否有害。因此，cca已引起已称为“眼狭窄反射(oculostenoticreflex)”的现象，即一些介入的心脏病专家对用cca发现的每一病变而言不论病变是否具有功能意义都插入血管内支架。因此，cca可引起对患者不必要的操作，所述不必要的操作可对患者造成增加的风险并且可对患者引起不必要的健康护理成本。

在诊断性心导管插入术期间，冠状动脉病变的功能意义可通过测量所观察的病变的血流储备分数(fractionalflowreserve；ffr)以侵入方式评估。ffr定义为在冠状动脉血流量增大(例如由静脉内施用腺苷诱导)的条件下，病变下游的平均血压除以病变上游的平均血压(例如主动脉压)的比率。血压可通过将压力线插入患者体内来测量。因此，可在诊断性心导管插入术的初始成本和风险已产生之后依据所确定的ffr作出治疗病变的决定。

因此，需要一种用于以非侵入方式评估冠状动脉解剖结构、心肌灌注和冠状动脉流的方法。所述方法和系统可能有利于诊断疑患有冠状动脉疾病的患者并且为所述患者订治疗计划的心脏病专家。另外，需要一种在不能直接测量的条件(例如运动)下预测冠状动脉流量和心肌灌注，并且预测对冠状动脉血流和心肌灌注的医药治疗、介入治疗和手术治疗的结果的方法。

应了解，前述一般描述和以下详细描述都仅具有示例性和说明性并且不限制公开内容。

技术实现要素：

根据一个实施方案，用于确定患者心血管信息的系统包括至少一个计算机系统，所述计算机系统被配置成接收关于患者心脏几何形状的患者特异性数据并且依据患者特异性数据创建呈现患者心脏的至少一部分的三维模型。至少一个计算机系统进一步被配置成创建有关患者心脏血流特征的基于物理学的模型，并且依据三维模型和基于物理学的模型来确定患者心脏内的血流储备分数。

根据另一实施方案，使用至少一个计算机系统的用于确定患者心血管信息的方法包括将关于患者心脏的几何形状的患者特异性数据输入至少一个计算机系统中，以及使用至少一个计算机系统依据患者特异性数据来创建呈现患者心脏的至少一部分的三维模型。方法还包括使用至少一个计算机系统创建有关患者心脏血流特征的基于物理学的模型，以及使用至少一个计算机系统依据三维模型和基于物理学的模型来确定患者心脏内的血流储备分数。

根据另一实施方案，提供一种用于至少一个计算机系统的非临时性计算机可读介质，其含有用于进行一种用于确定患者特异性心血管信息的方法的计算机可执行编程指令。方法包括接收关于患者心脏几何形状的患者特异性数据，以及依据患者特异性数据创建呈现患者心脏的至少一部分的三维模型。方法还包括创建有关患者心脏中血流特征的基于物理学的模型，以及依据三维模型和基于物理学的模型来确定患者心脏内血流储备分数。

根据另一实施方案，用于为患者订治疗计划的系统包括至少一个计算机系统，所述计算机系统被配置成接收关于患者解剖结构的几何形状的患者特异性数据，并且依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型。至少一个计算机系统进一步被配置成依据有关患者解剖结构的三维模型和基于物理学的模型来确定关于患者解剖结构内血流特征的第一信息、修改三维模型，以及依据经过修改的三维模型来确定关于患者解剖结构内血流特征的第二信息。

根据另一实施方案，提供一种用于计算机系统的非临时性计算机可读介质，其含有用于进行一种为患者订治疗计划方法的计算机可执行编程指令。方法包括接收关于患者解剖结构的几何形状的患者特异性数据，以及依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型。方法还包括依据有关患者解剖结构的三维模型和基于物理学的模型来确定关于患者解剖结构内血流特征的第一信息，以及依据所需的患者解剖结构几何形状变化来确定关于患者解剖结构内血流特征的第二信息。

根据另一实施方案，使用计算机系统的用于为患者订治疗计划的方法包括将关于患者解剖结构的几何形状的患者特异性数据输入至少一个计算机系统中，以及使用至少一个计算机系统、依据患者特异性数据来创建呈现患者解剖结构的至少一部分的三维模型。方法还包括使用至少一个计算机系统、依据有关患者解剖结构的三维模型和基于物理学的模型来确定关于患者解剖结构内血流特征的第一信息。方法还包括使用至少一个计算机系统来修改三维模型，以及使用至少一个计算机系统、依据经过修改的三维模型来确定关于患者解剖结构内血流特征的第二信息。

根据另一实施方案，用于为患者订治疗计划的系统包括至少一个计算机系统，所述计算机系统被配置成接收关于患者解剖结构的几何形状的患者特异性数据，并且依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型。至少一个计算机系统还被配置成依据三维模型和关于患者生理条件的信息来确定关于患者解剖结构内血流特征的第一信息、修改患者的生理条件，以及依据经过修改的患者生理条件来确定关于患者解剖结构内血流特征的第二信息。

根据另一实施方案，提供一种用于计算机系统的非临时性计算机可读介质，其含有用于为患者订治疗计划的方法的计算机可执行编程指令。方法包括接收关于患者解剖结构的几何形状的患者特异性数据，以及依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型。方法还包括依据三维模型和关于患者生理条件的信息来确定关于患者解剖结构内血流特征的第一信息，以及依据所需的患者生理条件变化来确定关于患者解剖结构内血流特征的第二信息。

根据另一实施方案，使用至少一个计算机系统的用于为患者订治疗计划的方法包括将关于患者解剖结构的几何形状的患者特异性数据输入至少一个计算机系统中，以及使用至少一个计算机系统、依据患者特异性数据来创建呈现患者解剖结构的至少一部分的三维模型。方法还包括使用至少一个计算机系统、依据三维模型和有关患者生理条件的信息来确定关于患者解剖结构内血流特征的第一信息。方法还包括使用至少一个计算机系统来修改患者生理条件，以及使用至少一个计算机系统、依据经过修改的患者生理条件来确定关于患者解剖结构内血流特征的第二信息。

根据另一实施方案，用于确定患者特异性心血管信息的系统包括至少一个计算机系统，所述计算机系统被配置成接收关于患者解剖结构的几何形状的患者特异性数据，并且依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型。至少一个计算机系统还被配置成确定与经过患者解剖结构的部分的总流量相关的总阻力，并且依据有关患者解剖结构的三维模型、基于物理学的模型以及所确定的总阻力来确定关于患者解剖结构内血流特征的信息。

根据另一实施方案，用于使用至少一个计算机系统来确定患者特异性心血管信息的方法包括将关于患者解剖结构的几何形状的患者特异性数据输入至少一个计算机系统中，以及使用至少一个计算机系统、依据患者特异性数据来创建呈现患者解剖结构的至少一部分的三维模型。方法还包括使用至少一个计算机来确定与经过患者解剖结构的部分的总流量相关的总阻力，以及使用至少一个计算机、依据有关患者解剖结构的三维模型、基于物理学的模型以及所确定的总阻力来确定关于患者解剖结构内血流特征的信息。

根据另一实施方案，提供一种用于计算机系统的非临时性计算机可读介质，其含有用于进行一种用于确定患者特异性心血管信息的方法的计算机可执行编程指令。方法包括接收关于患者解剖结构的几何形状的患者特异性数据，以及依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型。方法还包括确定与经过患者解剖结构的部分的总流量相关的总阻力，以及依据有关患者解剖结构的三维模型、基于物理学的模型以及所确定的总阻力来确定关于患者解剖结构内血流特征的信息。

根据另一实施方案，一种用于使用网址来提供患者特异性心血管信息的系统包括至少一个计算机系统，所述计算机系统被配置成允许远程用户访问网址、接收关于患者解剖结构几何形状的至少一部分的患者特异性数据、依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型，并且依据患者的三维模型和生理条件来确定关于患者解剖结构内血流特征的信息。至少一个计算机系统还被配置成使用网址向远程用户传达关于患者解剖结构的至少所述部分的第一个三维模拟的显示信息。三维模拟包括所确定的关于血流特征的信息。

根据另一实施方案，一种用于使用网址来提供患者特异性心血管信息的方法包括使用至少一个计算机系统允许远程用户访问网址，以及使用至少一个计算机系统接收关于患者解剖结构的几何形状的患者特异性数据。方法还包括使用至少一个计算机系统、依据患者特异性数据来创建呈现患者解剖结构的至少一部分的三维模型，以及使用至少一个计算机系统、依据三维模型和患者的生理条件来确定关于患者解剖结构内血流特征的信息。方法还包括使用至少一个计算机系统向使用网址的远程用户传达关于患者解剖结构的至少所述部分的第一个三维模拟的显示信息。三维模拟包括所确定的关于血流特征的信息。

根据另一实施方案，提供一种用于计算机系统的非临时性计算机可读介质，其含有用于进行一种用于使用网址来提供患者特异性心血管信息的方法的计算机可执行编程指令。方法包括允许远程用户访问网址、接收关于患者解剖结构的几何形状的患者特异性数据，以及依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型。方法还包括依据有关患者解剖结构的三维模型和基于物理学的模型来确定关于患者解剖结构内血流特征的信息，以及向使用网址的远程用户传达关于患者解剖结构的至少所述部分的第一个三维模拟的显示信息。三维模拟包括所确定的关于血流特征的信息。

根据另一实施方案，用于确定患者特异性的时变心血管信息的系统包括至少一个计算机系统，所述计算机系统被配置成在不同时刻接收关于患者解剖结构的至少一部分的几何形状的时变患者特异性数据，并且依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型。至少一个计算机系统还被配置成依据有关患者解剖结构的三维模型和基于物理学的模型来确定关于患者解剖结构内血流特征随时间变化的信息。

根据另一实施方案，用于使用至少一个计算机系统来确定患者特异性的时变心血管信息的方法包括使用至少一个计算机系统在不同时刻接收关于患者解剖结构几何形状的时变患者特异性数据。方法还包括使用至少一个计算机系统、依据患者特异性数据来创建呈现患者解剖结构的至少一部分的三维模型。方法还包括使用至少一个计算机系统、依据有关患者解剖结构的三维模型和关于基于物理学的模型的信息来确定关于患者解剖结构内血流特征随时间变化的信息。

根据另一实施方案，提供一种用于计算机系统的非临时性计算机可读介质，其含有用于进行一种用于确定患者特异性的时变心血管信息的方法的计算机可执行编程指令。方法包括在不同时刻接收随时间变化的关于患者解剖结构几何形状的患者特异性数据、依据患者特异性数据创建呈现患者解剖结构的至少一部分的三维模型，以及依据有关患者解剖结构的三维模型和关于基于物理学的模型的信息来确定关于患者解剖结构内血流特征随时间变化的信息。

根据另一实施方案，用于确定患者心血管信息的系统包括至少一个计算机系统，所述至少一个计算机系统被配置成接收关于患者解剖结构的至少一部分的几何形状和至少一种物质性质的患者特异性数据。解剖结构包括血管的至少一部分。至少一个计算机系统还被配置成依据患者特异性数据来创建呈现患者解剖结构的三维模型，并且依据三维模型和患者的生理条件来确定关于患者解剖结构内血流特征的信息。至少一个计算机系统还被配置成鉴别血管内斑块的位置。

根据另一实施方案，用于使用至少一个计算机系统来确定患者心血管信息的方法包括使用至少一个计算机系统来接收关于患者解剖结构的至少一部分的几何形状和至少一种物质性质的患者特异性数据。解剖结构包括血管的至少一部分。方法还包括使用至少一个计算机系统、依据患者特异性数据来创建呈现患者解剖结构的三维模型，以及使用至少一个计算机系统、依据三维模型和患者的生理条件来确定关于患者解剖结构内血流特征的信息。方法还包括使用至少一个计算机系统鉴别血管内的斑块。

根据另一实施方案，提供一种用于计算机系统的非临时性计算机可读介质，其含有用于进行一种用于确定患者心血管信息的方法的计算机可执行编程指令。方法包括接收关于患者解剖结构的至少一部分的几何形状和至少一种物质性质的患者特异性数据。解剖结构包括血管的至少一部分。方法还包括依据患者特异性数据来创建呈现患者解剖结构的三维模型、依据三维模型和患者的生理条件来确定关于患者解剖结构内血流特征的信息，以及鉴别血管内斑块的位置。

根据另一实施方案，用于确定患者心血管信息的系统包括至少一个计算机系统，所述至少一个计算机系统被配置成接收关于患者解剖结构的至少一部分的几何形状的患者特异性数据。解剖结构包括多条动脉的至少一部分和与多条动脉的至少一部分连接的组织。至少一个计算机系统还被配置成依据患者特异性数据创建呈现患者解剖结构的三维模型、将呈现组织的三维模型的至少一部分分成片段，并且依据三维模型和患者的生理条件来确定关于与至少一个片段相关的血流特征的信息。

根据另一实施方案，用于使用至少一个计算机系统来确定患者心血管信息的方法包括使用至少一个计算机系统来接收关于患者解剖结构的至少一部分的几何形状的患者特异性数据。解剖结构包括多条动脉的至少一部分和与多条动脉的至少一部分连接的组织。方法还包括使用至少一个计算机系统、依据患者特异性数据来创建呈现患者解剖结构的三维模型，以及使用至少一个计算机系统来扩展三维模型以形成扩充模型。方法还包括使用至少一个计算机系统将呈现组织的扩充模型的至少一部分分成片段，以及使用至少一个计算机系统、依据扩充模型和患者的生理条件来确定关于与至少一个片段相关的血流特征的信息。

根据另一实施方案，提供一种用于计算机系统的非临时性计算机可读介质，其含有用于进行一种用于确定患者心血管信息的方法的计算机可执行编程指令。方法包括接收关于患者解剖结构的至少一部分的几何形状的患者特异性数据。解剖结构包括多条动脉的至少一部分和与多条动脉的至少一部分连接的组织。方法还包括依据患者特异性数据创建呈现患者解剖结构的三维模型、将呈现组织的三维模型的至少一部分分成片段，以及依据有关解剖结构的三维模型和基于物理学的模型来确定关于与至少一个片段相关的血流特征的信息。

根据另一实施方案，用于确定患者心血管信息的系统包括至少一个计算机系统，所述至少一个计算机系统被配置成接收关于患者脑部几何形状的患者特异性数据。至少一个计算机系统还被配置成依据患者特异性数据来创建呈现患者脑部的至少一部分的三维模型，并且依据有关患者脑部的三维模型和基于物理学的模型来确定关于患者脑内血流特征的信息。

根据另一实施方案，用于使用至少一个计算机系统来确定患者特异性心血管信息的方法包括将关于患者多根脑动脉的至少一部分的几何形状的患者特异性数据输入至少一个计算机系统中。方法还包括使用至少一个计算机系统、依据患者特异性数据来创建呈现患者脑动脉的至少所述部分的三维模型，以及使用至少一个计算机系统、依据有关患者脑动脉的三维模型和基于物理学的模型来确定关于患者脑动脉内血流特征的信息。

根据另一实施方案，提供一种用于至少一个计算机系统的非临时性计算机可读介质，其含有用于进行一种用于确定患者特异性心血管信息的方法的计算机可执行编程指令。方法包括接收关于患者脑部几何形状的患者特异性数据、依据患者特异性数据来创建呈现患者脑部的至少一部分的三维模型，以及依据有关患者脑部的三维模型和基于物理学的模型来确定关于患者脑内血流特征的信息。

其它实施方案和优势将在以下说明书的部分中加以阐述，并且部分地将根据说明书而显而易见，或可通过所述公开内容的实践得到学习。实施方案和优势将借助于下文详细指出的要素和组合来实现和达到。

附图说明

并入本说明书中并且构成本说明书的一部分的随附图式说明几个实施方案，并且连同说明书一起用以解释公开内容的原理。

图1是根据一个示例性实施方案的一种系统的示意图，所述系统用于提供有关特定患者的冠状动脉血流的各种信息；

图2是根据一个示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的血流的各种信息；

图3是展示图2的方法的子步骤的流程图；

图4展示根据一个示例性实施方案以非侵入方式从患者获得的成像数据；

图5展示使用图4的成像数据产生的示例性三维模型；

图6展示一层图4的成像数据的一部分，包括用于形成第一初始模型的种子；

图7展示通过扩展图6的种子形成的第一初始模型的一部分；

图8展示根据一个示例性实施方案的经过修剪的实体模型；

图9展示当患者休息时的示例性计算ffr(cffr)模型；

图10展示当患者处于最大充血时的示例性cffr模型；

图11展示当患者处于最大运动时的示例性cffr模型；

图12展示根据一个示例性实施方案，提供用于形成集中参数模型的经过修剪的实体模型的一部分；

图13展示提供用于形成集中参数模型的图12的经过修剪的实体模型的中线的一部分；

图14展示提供用于形成集中参数模型的依据图12的经过修剪的实体模型形成的片段；

图15展示提供用于形成集中参数模型的以电阻器替代的图14的片段；

图16展示根据一个示例性实施方案的示例性集中参数模型，所述集中参数模型呈现在实体模型边界的流入和流出处的上游结构和下游结构；

图17展示依据图8的实体模型准备的三维网；

图18和图19展示图17的三维网的部分；

图20展示患者解剖结构的包括血流信息的模型，其中模型上的某些点由个别参考标签鉴别；

图21是在主动脉中以及在图20中鉴别的一些点处的模拟血压随时间变化的图；

图22是在图20中鉴别的各点处的模拟血流随时间变化的图；

图23是根据一个示例性实施方案的定稿报道；

图24是根据一个示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的冠状动脉血流的各种信息；

图25展示根据一个示例性实施方案的经过修改的cffr模型，其依据通过加宽左前降支(lad)动脉的一部分和lcx动脉的一部分创建的实体模型来确定；

图26展示根据一个示例性实施方案，在加宽lad动脉的一部分和左旋支(lcx)动脉的一部分之后，经过修改的模拟血流模型的一个实施例；

图27是根据一个示例性实施方案的一种方法的流程图，所述方法用于使用降阶模型来模拟各种治疗选择；

图28是根据另一示例性实施方案的一种方法的流程图，所述方法用于使用降阶模型来模拟各种治疗选择；

图29是根据一个示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的心肌灌注的各种信息；

图30是根据另一示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的心肌灌注的各种信息；

图31展示根据一个示例性实施方案的一种患者特异性模型，所述模型提供有关心肌灌注的各种信息；

图32是根据另一示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的心肌灌注的各种信息；

图33是沿血管壁积累的斑块的截面图；

图34展示根据一个示例性实施方案的一种患者特异性模型，所述模型提供有关斑块脆弱性的各种信息；

图35是根据一个示例性实施方案的一种方法的流程图，所述方法用于提供有关评估特定患者的斑块脆弱性、心肌容量风险和心肌灌注风险的各种信息；

图36是展示根据一个示例性实施方案的信息的示意图，所述信息获自图35的方法；

图37是脑动脉图；

图38是根据一个示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的颅内血流和颅外血流的各种信息；

图39是根据一个示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的大脑灌注的各种信息；

图40是根据另一示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的大脑灌注的各种信息；

图41是根据另一示例性实施方案的一种方法的流程图，所述方法用于提供有关特定患者的大脑灌注的各种信息；以及

图42是根据一个示例性实施方案的一种方法的流程图，所述方法用于提供有关评估特定患者的斑块脆弱性、大脑容量风险和大脑灌注风险的各种信息。

具体实施方式

现将详细参考示例性实施方案，所述示例性实施方案的实施例在随附图式中得到说明。只要可能，相同参考数字将遍及图式被使用以指相同或类似部分。这个描述是根据以下大纲来组织：

i.概述

ii.获得并预处理患者特异性解剖学数据

iii.依据所得解剖学数据来创建三维模型

iv.使模型为分析做准备以及确定边界条件

a.使模型为分析做准备

b.确定边界条件

i.确定降阶模型

ii.示例性集中参数模型

c.创建三维网

v.进行计算分析以及输出结果

a.进行计算分析

b.显示血压、流量和cffr的结果

c.验证结果

d.用于提供冠状动脉血流信息的系统和方法的另一实施方案

vi.提供患者特异性治疗计划

a.使用降阶模型来比较不同治疗选择

vii.其它结果

a.评估心肌灌注

b.评估斑块脆弱性

viii.其它应用

a,为颅内和颅外血流建模

i.评估大脑灌注

ii.评估斑块脆弱性

i.概述

在一个示例性实施方案中，方法和系统使用以非侵入方式从患者取得的信息来确定有关特定患者的血流的各种信息。所确定的信息可涉及患者冠状动脉血管结构中的血流。或者，如下文更详细地描述，所确定的信息可涉及患者血管结构(如颈动脉血管结构、周边血管结构、腹部血管结构、肾血管结构和大脑血管结构)的其它区域中的血流。冠状动脉血管结构包括在大动脉到小动脉、毛细管、小静脉、静脉等的范围内的复杂血管网状结构。冠状动脉血管结构使血液循环到心脏并在心脏内循环，并且包括主动脉2(图5)，其向多条主要冠状动脉4(图5)(例如左前降支(lad)动脉、左旋支(lcx)动脉、右冠状(rca)动脉等)提供血液，其还可分成动脉的分支，或主动脉2和主要冠状动脉4下游的其它类型血管的分支。因此，示例性方法和系统可确定有关主动脉、主要冠状动脉和/或主要冠状动脉下游的其它冠状动脉或血管内血流的各种信息。尽管下文讨论了主动脉和冠状动脉(和由此延伸的分支)，但所公开的方法和系统也可适用于其它类型血管。

在一个示例性实施方案中，通过公开的方法和系统确定的信息可包括但不限于主动脉、主要冠状动脉和/或主要冠状动脉下游的其它冠状动脉或血管中各种位置的各种血流特征或参数，如血液流动速度、压力(或其比率)、流动速率和ffr。这信息可用以确定病变是否具有功能意义和/或是否要治疗病变。这信息可使用以非侵入方式从患者获得的信息来确定。因此，是否要治疗病变的决定可在没有与侵入性程序相关的成本和风险下作出。

图1展示根据一个示例性实施方案的一种系统的方面，所述系统用于提供有关特定患者的冠状动脉血流的各种信息。患者解剖结构的三维模型10可使用以非侵入方式从患者获得的数据创建，如下文较详细描述。其它患者特异性信息也可以非侵入方式获得。在一个示例性实施方案中，由三维模型10表示的患者解剖结构的部分可包括主动脉的至少一部分，和与主动脉连接的主要冠状动脉(和由此延伸或发源的分支)的邻近部分。

有关冠状动脉血流的各种生理规律或关系20可例如从实验数据推导，如下文较详细描述。使用三维解剖模型10和所推导的生理规律20，有关冠状动脉血流的多个等式30可如下文较详细描述来确定。举例来说，等式30可使用任何数值方法(例如有限差分方法、有限体积方法、光谱方法、晶格波耳兹曼(latticeboltzmann)方法、基于颗粒的方法、水平集(levelset)方法、有限元方法等)来确定和求解。等式30可能是可求解的以便确定关于在由模型10表示的解剖结构中各种点的患者解剖结构中冠状动脉血流的信息(例如压力、速度、ffr等)。

等式30可使用计算机40来求解。依据已求解的等式，计算机40可输出一个或多个图像或模拟结果，其指示有关由模型10表示的患者解剖结构中血流的信息。举例来说，图像可包括模拟血压模型50、模拟血液流量或速度模型52、计算ffr(cffr)模型54等，如下文将更详细地描述。模拟血压模型50、模拟血流模型52和cffr模型54提供关于以下的信息：在由模型10表示的患者解剖结构中沿三维的各种位置的各别压力、速度和cffr。在冠状动脉血流增大(例如常规地由静脉内施用腺苷所诱导)的条件下，cffr可计算为模型10中具体位置的血压除以主动脉中(例如在模型10的流入边界处)的血压的比率。

在一个示例性实施方案中，计算机40可包括一个或多个存储指令的非临时性计算机可读存储设备，所述指令当由处理器、计算机系统等执行时可进行任一本文所述的用于提供有关患者血流的各种信息的行为。计算机40可包括台式计算机或便携式计算机、工作站、服务器、个人数字助理或任何其它计算机系统。计算机40可包括处理器，只读存储器(rom),随机存取存储器(ram),用于连接外围设备(例如输入设备，输出设备，存储设备等)的输入/输出(i/o)适配器，用于连接输入设备，如键盘，鼠标，触摸屏，语音输入，和/或其他设备的用户界面适配器，用于连接计算机40至网络的通讯适配器，用于连接计算机40至显示器的显示器适配器等。举例来说，显示器可用以显示三维模型10和/或通过求解等式30产生的任何图像，如模拟血压模型50、模拟血流模型52和/或cffr模型54。

图2展示根据另一示例性实施方案的一种方法的方面，所述方法用于提供有关特定患者的血流的各种信息。方法可包括获得患者特异性解剖学数据，如关于以下的信息：患者解剖结构(例如主动脉的至少一部分和与主动脉连接的主要冠状动脉的邻近部分(和由此延伸的分支))，以及预处理所述数据(步骤100)。患者特异性解剖学数据可以非侵入方式，例如通过ccta获得，如下文将描述。

患者解剖结构的三维模型可依据所得解剖学数据创建(步骤200)。举例来说，三维模型可为上文关于图1描述的患者解剖结构的三维模型10。

可使三维模型为分析做准备，并且可确定边界条件(步骤300)。举例来说，上文关于图1描述的患者解剖结构的三维模型10可被修剪并且离散化成体积网(volumetricmesh)，例如有限元网或有限体积网。体积网可用以产生上文关于图1描述的等式30。

边界条件也可被赋值并且并入上文关于图1描述的等式30中。边界条件提供关于三维模型10在它边界处的信息，例如流入边界322(图8)、流出边界324(图8)、血管壁边界326(图8)，等。流入边界322可包括流被引入三维模型的解剖结构中所经过的边界，如在接近主动脉根的主动脉末端(例如图16中所示的a端)。各流入边界322可例如用速度、流动速率、压力或其它特征的指定值或值域、通过将心脏模型和/或集中参数模型与边界偶合等来赋值。流出边界324可包括流被向外引出三维模型的解剖结构所经过的边界，如在接近主动脉弓的主动脉末端(例如图16中所示的b端)，以及主要冠状动脉和由此延伸的分支的下游末端(例如图16中所示的a-m端)。各流出边界可例如通过偶合集中参数或分布(例如一维波传播)模型来赋值，所述模型将在下文详细描述。流入和/或流出边界条件的指定值可通过以非侵入方式测量患者生理学特征来确定，患者生理学特征如但不限于心输出量(来自心脏的血流体积)、血压、心肌质量，等。血管壁边界可包括三维模型10的主动脉、主要冠状动脉和/或其它冠状动脉或血管的实体边界。

计算分析可使用所准备的三维模型和所确定的边界条件来进行(步骤400)以确定患者的血流信息。举例来说，计算分析可用等式30并且使用上文关于图1描述的计算机40来进行，以产生上文关于图1描述的图像，如模拟血压模型50、模拟血流模型52和/或cffr模型54。

方法还可包括使用结果来提供患者特异性治疗选择(步骤500)。举例来说，在步骤200中创建的三维模型10和/或在步骤300中赋值的边界条件可得到调整，以便为一种或多种治疗(例如将冠状动脉血管内支架放在三维模型10中呈现的冠状动脉之一中)或其它治疗选择建模。接着，计算分析可如上在步骤400中所述来进行，以产生新图像，如更新版本的血压模型50、血流模型52和/或cffr模型54。如果采用所述治疗选择，那么这些新图像可用以确定血液流动速度和压力的变化。

本文公开的系统和方法可并入医师可得的软件工具中，以提供定量冠状动脉中的血流并且评估冠状动脉疾病的功能意义的非侵入性手段。另外，医师可使用软件工具来预测医药治疗、介入治疗和/或手术治疗对冠状动脉血流的影响。软件工具可预防、诊断、管理和/或治疗心血管系统的其它部分的疾病，所述其它部分包括颈部动脉(例如颈动脉)、头部动脉(例如脑动脉)、胸部动脉、腹部动脉(例如腹主动脉和它的分支)、臂动脉或腿动脉(例如大腿动脉和腿弯部动脉)。软件工具可具有互动性以使医师能够为患者产生最优的个人化疗法。

举例来说，软件工具可至少部分并入医师或其它用户使用的计算机系统，例如图1中所示的计算机40中。计算机系统可接收以非侵入方式从患者获得的数据(例如用以创建三维模型10的数据、用以应用边界条件或进行计算分析的数据，等)。举例来说，数据可由医师输入或可从能够存取和提供所述数据的另一来源(如放射实验室或其它医药实验室)接收。数据可经由网络或用于传递数据的其它系统传输，或直接传入计算机系统。软件工具可使用所述数据以产生并显示三维模型10或其它模型/网和/或任何模拟或其它结果(通过求解上文关于图1描述的等式30来确定)，如模拟血压模型50、模拟血流模型52和/或cffr模型54。因此，软件工具可进行步骤100-500。在步骤500中，医师可向计算机系统提供其它输入以选择可能的治疗选择，并且计算机系统可向医师显示依据所选可能的治疗选择的新模拟。此外，图2中所示的步骤100-500中每一个都可使用单独软件包或软件模块来进行。

或者，软件工具可作为基于网络的服务或其它服务(例如与医师分开的实体提供的服务)的部分来提供。服务提供者可例如操作基于网络的服务，并且可提供网络入口或其它基于网络的应用程序(例如在服务器或由服务提供者操作的其它计算机系统上运行)，所述网络入口或其它基于网络的应用程序为医师或其它用户经由网络或在计算机系统之间传递数据的其它方法可得。举例来说，以非侵入方式从患者获得的数据可提供给服务提供者，并且服务提供者可使用所述数据来产生三维模型10或其它模型/网和/或任何模拟或其它结果(通过求解上文关于图1描述的等式30来确定)，如模拟血压模型50、模拟血流模型52和/或cffr模型54。接着，基于网络的服务可传输有关三维模型10或其它模型/网和/或模拟的信息，以便三维模型10和/或模拟可在医师的计算机系统上向医师显示。因此，基于网络的服务可进行步骤100-500和用于提供患者特异性信息的下述任何其它步骤。在步骤500中，医师可提供其它输入，例如以选择可能的治疗选择或对计算分析进行其它调整，并且输入可被传输到由服务提供者操作的计算机系统(例如经由网络入口)。基于网络的服务可依据所选的可能治疗选择来产生新的模拟或其它结果，并且可将有关新模拟的信息传回医师，以便新模拟可向医师显示。

应了解，本文所述的步骤中一个或多个可由以下进行：一个或多个人操作员(例如心脏病专家或其它医师、患者、提供基于网络的服务或其它由第三方提供的服务的服务提供者的员工、其它用户等)，或一个或多个由所述人操作员使用的计算机系统，如台式计算机或手提式计算机、工作站、服务器、个人数字助理，等。可经由网络或传递数据的其它方法来连接计算机系统。

图3展示示例性方法的其它方面，所述示例性方法用于提供有关特定患者的血流的各种信息。图3中所示的方面可并入软件工具中，所述软件工具可至少部分并入计算机系统中和/或作为基于网络的服务的部分。

ii.获得并预处理患者特异性解剖学数据

如上关于图2中所示的步骤100所述，示例性方法可包括获得患者特异性解剖学数据，如关于患者心脏的信息，以及预处理所述数据。在一个示例性实施方案中，步骤100可包括以下步骤。

最初，可选择患者。举例来说，当医师确定需要关于患者冠状动脉血流的信息时，例如，如果患者经历与冠状动脉疾病相关的症状，如胸痛、心脏病发作等，那么患者可由医师选择。

可获得患者特异性解剖学数据，如关于患者心脏的几何形状的数据，所述几何形状例如为患者主动脉的至少一部分、与主动脉连接的主要冠状动脉(和由此延伸的分支)的邻近部分，以及心肌。患者特异性解剖学数据可以非侵入方式，例如使用非侵入性成像方法获得。举例来说，ccta是一种成像方法，其中用户可操作计算机断层摄影(ct)扫描仪来观察和创建结构图像，所述结构例如是心肌、主动脉、主要冠状动脉和与其连接的其它血管。ccta数据可随时间变化，例如来展示血管形状在心动周期内的变化。ccta可用以产生患者心脏的图像。举例来说，可获得64层ccta数据，例如有关患者心脏的64层的数据，并且将所述数据收集到三维图像中。图4展示由64层ccta数据产生的三维图像120的一个实施例。

或者，其它非侵入性成像方法，如磁共振成像(mri)或超声波(us)，或侵入性成像方法，如数字减影血管造影(dsa)可用以产生患者解剖学结构的图像。成像方法可涉及向患者静脉内注射对比剂以实现解剖学结构的鉴别。所得成像数据(例如由ccta、mri等提供)可由第三方供应者提供，如由放射学实验室或心脏病专家，由患者医师等提供。

其它患者特异性解剖学数据还可以非侵入方式从患者确定。举例来说，可测量生理学数据，如患者的血压、基线心率、身高、体重、血细胞比容、心搏量等。血压可为患者肱动脉中的血压(例如使用压力袖套)，如最大(收缩)压力和最小(舒张)压力。

如上所述获得的患者特异性解剖学数据可经安全通信线(例如经由网络)转移。举例来说，数据可被转移到服务器或用于进行计算分析(例如上文在步骤400中描述的计算分析)的其它计算机系统。在一个示例性实施方案中，数据可被转移到服务器或由提供基于网络的服务的服务提供者操作的其它计算机系统。或者，数据可被转移到由患者的医师或其它用户操作的计算机系统。

回头参考图3，可检查所转移的数据以确定数据是否为可接受的(步骤102)。所述确定可由用户和/或由计算机系统进行。举例来说，所转移的数据(例如ccta数据和其它数据)可由用户和/或由计算机系统验证，例如以确定ccta数据是否完整(例如包括主动脉和主要冠状动脉的足够部分)并且对应于正确患者。

所转移的数据(例如ccta数据和其它数据)还可得到预处理和评估。预处理和/或评估可由用户和/或计算机系统进行，并且可包括例如检查ccta数据中的重合失调(misregistration)、不一致或模糊、检查ccta数据中所示的血管内支架、检查可能妨碍血管内腔可视性的其它人工制品、检查患者的结构(例如主动脉、主要冠状动脉和其它血管)与其它部分之间的足够对比，等。

可评价所转移的数据以依据上述验证、预处理和/或评估来确定数据是否可接受。在上述验证、预处理和/或评估期间，用户和/或计算机系统也许能校正数据的某些错误或问题。然而，如果有太多错误或问题，那么可确定数据为不可接受，并且用户和/或计算机系统可产生拒绝报告，所述拒绝报告解释错误或问题迫使拒绝所转移的数据。任选地，可进行新的ccta扫描和/或可再次从患者测量上述生理学数据。如果确定所转移的数据为可接受的，那么所述方法可往下进行到下述步骤202。

因此，图3中所示并且上文描述的步骤102可被视为图2的步骤100的子步骤。

iii.依据所得解剖学数据来创建三维模型

如上文关于图2中所示的步骤200所述，示例性方法可包括依据所得解剖学数据来创建三维模型。在一个示例性实施方案中，步骤200可包括以下步骤。

可使用ccta数据来产生冠状血管的三维模型。图5展示使用ccta数据产生的三维模型220的表面的一个实施例。举例来说，模型220可包括例如主动脉的至少一部分、一根或多根与主动脉的所述部分连接的主要冠状动脉的至少邻近部分，一根或多根与主要冠状动脉连接的分支的至少邻近部分，等。主动脉、主要冠状动脉和/或分支的建模部分可互连并且为枝状以使没有哪部分从模型220的其余部分断开。形成模型220的方法称为分割。

回头参考图3，计算机系统可自动分割主动脉的至少一部分(步骤202)并且分割心肌(或其它心脏组织，或与待建模的动脉连接的其它组织)的至少一部分(步骤204)。计算机系统还可分割与主动脉连接的主要冠状动脉的至少一部分。在一个示例性实施方案中，计算机系统可允许用户选择一个或多个冠状动脉根或起始点(步骤206)以分割主要冠状动脉。

分割可使用各种方法进行。分割可由计算机系统依据用户输入或在无用户输入下自动进行。举例来说，在一个示例性实施方案中，用户可向计算机系统提供输入以产生第一初始模型。举例来说，计算机系统可向用户显示从ccta数据产生的三维图像120(图4)或其层。三维图像120可包括变化强度的光的部分。举例来说，较亮区域可指示主动脉、主要冠状动脉和/或分支的内腔。较暗区域可指示患者心脏的心肌和其它组织。

图6展示可向用户显示的三维图像120的层222的一部分，并且层222可包括具有相对亮度的区域224。计算机系统可允许用户通过添加一种或多种种子226来选择相对亮度的区域224，并且种子226可充当用于分割主要冠状动脉的冠状动脉根或起始点。在用户作出指令时，计算机系统则可使用种子226作为起始点以形成第一初始模型。用户可将种子226加入主动脉和/或个别主要冠状动脉中一根或多根中。任选地，用户还可将种子226加入与主要冠状动脉连接的分支中一根或多根中。或者，计算机系统可例如使用提取的中线信息来自动放置种子。计算机系统可确定图像120的强度值，其中种子226已被放置并且可如下形成第一初始模型：沿具有相同强度值(或在集中于所选强度值的强度值的范围或临限值内)的图像120的部分扩展种子226。因此，这种分割方法可被称作“基于临限值的分割”。

图7展示通过扩展图6的种子226形成的第一初始模型的一部分230。因此，用户输入种子226作为计算机系统开始形成第一初始模型的起始点。可重复这种方法直到整个目标部分(例如主动脉和/或主要冠状动脉的部分)被分割为止。或者，第一初始模型可由计算机系统在没有用户输入下产生。

或者，分割可使用一种称为“基于边缘的分割”的方法进行。在一个示例性实施方案中，基于临限值的割方法和基于边缘的分割方法可如下文将描述来进行，以形成模型220。

第二初始模型可使用基于边缘的分割方法来形成。用这种方法可使主动脉和/或主要冠状动脉的内腔边缘定位。举例来说，在一个示例性实施方案中，用户可向计算机系统提供输入，例如，如上所述的种子226，以产生第二初始模型。计算机系统可沿图像120的部分扩展种子226直到达到边缘为止。内腔边缘可例如由用户以视觉定位，和/或由计算机系统定位(例如在强度值变化在设定临限值以上的位置)。基于边缘的分割方法可由计算机系统和/或用户进行。

心肌或其它组织也可在步骤204中依据ccta数据来分割。举例来说，可分析ccta数据来确定心肌(例如左心室和/或右心室)内表面和外表面的位置。表面位置可依据在ccta数据中，心肌相比心脏其它结构的对比度(例如相对明暗度)来确定。因此，可确定心肌的几何形状。

必要时可检查和/或校正主动脉、心肌和/或主要冠状动脉的分割(步骤208)。检查和/或校正可由计算机系统和/或用户进行。举例来说，在一个示例性实施方案中，如果有任何错误，例如，如果模型220中主动脉、心肌和/或主要冠状动脉的任何部分被遗漏或不精确，那么计算机系统可自动检查分割，并且用户可手动校正分割。

举例来说，可比较上述第一初始模型和第二初始模型以确保主动脉和/或主要冠状动脉的分割是精确的。第一初始模型和第二初始模型之间的差异的任何区域都可加以比较以校正分割并且形成模型220。举例来说，模型220可为第一初始模型和第二初始模型之间的平均。或者，可进行上述分割方法中的仅一种，并且通过所述方法形成的初始模型可用作模型220。

可计算心肌质量(步骤240)。计算可由计算机系统进行。举例来说，心肌体积可依据如上所述确定的心肌表面位置来计算，并且所计算的心肌体积可乘以心肌密度来计算心肌质量。心肌密度可预设。

模型220(图5)的各种血管(例如主动脉、主要冠状动脉等)的中线可得到确定(步骤242)。在一个示例性实施方案中，所述确定可由计算机系统自动进行。

必要时可检查和/或校正在步骤242中确定的中线(步骤244)。检查和/或校正可由计算机系统和/或用户进行。举例来说，在一个示例性实施方案中，如果有任何错误，例如，如果任何中线被遗漏或不精确，那么计算机系统可自动检查中线，并且用户可手动校正中线。

可检测到钙或斑块(导致血管变窄)(步骤246)。在一个示例性实施方案中，计算机系统可自动检测斑块。举例来说，斑块可在三维图像120中检测到并且从模型220去除。斑块可三维图像120中的鉴别出，因为斑块显现为比主动脉、冠状动脉和/或分支的内腔甚至更亮的区域。因此，斑块可由计算机系统检测为具有低于设定值的强度值，或可由用户视觉检测到。在检测到斑块之后，计算机系统可从模型220去除斑块，以便斑块不会被认为是血管中内腔或开放空间的部分。或者，计算机系统可使用不同于主动脉、主要冠状动脉和/或分支的颜色、色调(shading)或其它视觉指示物来指示模型220上的斑块。

计算机系统也可自动分割检测到的斑块(步骤248)。举例来说，斑块可依据ccta数据来分割。可分析ccta数据以依据ccta数据中斑块相比心脏其它结构的对比(例如相对明暗度)来定位斑块(或其表面)。因此，也可确定斑块的几何形状。

必要时可检查和/或校正斑块的分割(步骤250)。检查和/或校正可由计算机系统和/或用户进行。举例来说，在一个示例性实施方案中，如果有任何错误，例如，如果任何斑块被遗漏或展示不精确，那么计算机系统可自动检查分割，并且用户可手动校正分割。

计算机系统可自动分割与主要冠状动脉连接的分支(步骤252)。举例来说，分支可使用类似于用于分割主要冠状动脉(例如图6和图7中所示)并且上文关于步骤206描述的方法来分割。计算机系统也可使用与上文关于步骤248和步骤250所述类似的方法自动分割经过分割的分支中的斑块。或者，分支(和其中所含的任何斑块)可与主要冠状动脉同时被分割(例如在步骤206中)。

必要时可检查和/或校正分支的分割(步骤254)。检查和/或校正可由计算机系统和/或用户进行。举例来说，在一个示例性实施方案中，如果有任何错误，例如，如果模型220中分支的任何部分被遗漏或不精确，那么计算机系统可自动检查分割，并且用户可手动校正分割。

如果定位任何重合失调、血管内支架或其它人工制品(例如在步骤102中检查ccta数据期间)，那么都可校正模型220(步骤256)。所述校正可由用户和/或由计算机系统进行。举例来说，如果定位重合失调或其它人工制品(例如影响内腔可视性的不一致、模糊、人工制品，等)，那么可检查和/或校正模型220以避免血管截面积的人为变化或错误变化(例如人为变窄)。如果定位血管内支架，那么可检查和/或校正模型220以指示血管内支架的位置和/或例如依据血管内支架的尺寸来校正定位血管内支架的血管的截面积。

也可独立检查模型220的分割(步骤258)。所述检查可由用户和/或由计算机系统进行。举例来说，用户和/或计算机系统也许能鉴别模型220的某些错误，如可校正错误和/或可能需要模型220至少部分重做或重分割的错误。如果鉴别出所述错误，那么分割可确定为不可接受，并且某些步骤，例如步骤202-208、步骤240-256中一个或多个会取决于错误而被重复。

如果模型220的分割被独立地验证为可接受，那么任选地，模型220可被输出并且光滑化(步骤260)。所述光滑化可由用户和/或由计算机系统进行。举例来说，隆起、点或其它不连续部分可得到光滑化。可将模型220输出到准备用于计算分析等的单独软件模块。

因此，可将图3中展示并且上文描述的步骤202-208和步骤240-260视为图2的步骤200的子步骤。

iv.使模型为分析做准备以及确定边界条件

如上文关于图2中所示的步骤300所述，示例性方法可包括使模型为分析做准备以及确定边界条件的模型。在一个示例性实施方案中，步骤300可包括以下步骤。

a.使模型为分析做准备

回头参考图3，也可确定模型220(图5)的各种血管(例如主动脉、主要冠状动脉和/或分支)的截面积(步骤304)。在一个示例性实施方案中，所述确定可由计算机系统进行。

可修剪(步骤306)模型220(图5)并且可产生实体模型。图8展示依据类似于图5中所示的模型220的模型所准备的经过修剪的实体模型320的一个实施例。实体模型320是三维患者特异性几何模型。在一个示例性实施方案中，修剪可由计算机系统在有或无用户输入下进行。可修剪流入边界322和流出边界324中每一个，以便表面形成垂直于在步骤242中确定的中线的各别边界。流入边界322可包括流被引入模型320的解剖结构中所经过的边界，如在主动脉的上游端，如图8中所示。流出边界324可包括流被向外引出模型320的解剖结构中所经过的边界，如在主动脉的下游端和主要冠状动脉和/或分支的下游端。

b.确定边界条件

可提供边界条件以描述在模型(例如图8的三维实体模型320)的边界处在发生什么。举例来说，边界条件可涉及至少一个与患者建模解剖结构相关、例如在建模解剖结构的边界处的血流特征，并且血流特征可包括血液流动速度、压力、流动速率、ffr等。通过适当确定边界条件，可进行计算分析以确定在模型内各种位置的信息。现将描述边界条件和用于确定所述边界条件的方法的实施例。

在一个示例性实施方案中，所确定的边界条件可将由实体模型320表示的血管的部分的上游结构和下游结构简化成一维或二维降阶模型。用于确定边界条件的一组示例性等式和其它细节公开于例如美国专利申请公布号2010/0241404和美国临时申请号61/210,401中，这两个专利的标题都是“心血管系统的患者特异性血液动力学(patient-specifichemodynamicsofthecardiovascularsystem)”并且以引用的方式整体并入本文中。

边界条件可取决于患者生理条件而变化，因为经过心脏的血流可取决于患者生理条件而不同。举例来说，ffr通常在充血生理条件下测量，充血一般在患者例如因应激等而经历心脏血流增大时发生。ffr是在最大应激条件下冠状动脉压与主动脉压的比率。充血也可以药理学方式，例如用腺苷诱导。图9-11展示计算ffr(cffr)模型的实施例，其指示模型320中冠状动脉压与主动脉压的比率变化，所述变化取决于患者生理条件(在休息时、在最大充血下，或在最大运动下)。图9展示在整个模型320中，当患者休息时，冠状动脉压与主动脉压的比率的最小变化。图10展示在整个模型320中，当患者经历最大充血时，冠状动脉压与主动脉压的比率的较大变化。图11展示在整个模型320中，当患者经历最大运动时，冠状动脉压与主动脉压的比率的甚至更大变化。

回头参考图3，可确定充血条件的边界条件(步骤310)。在一个示例性实施方案中，腺苷的作用可使冠状动脉阻力降为1/5-1、主动脉血压降低约0-20％并且使心跳速率增大约0-20％来建模。举例来说，腺苷的作用可使冠状动脉阻力降为1/4、主动脉血压降低约10％并且使心跳速率增大约10％来建模。尽管在示例性实施方案中确定了充血条件的边界条件，但应了解，也可确定其它生理状态的边界条件，其它生理状态如为休息、不同程度充血、不同程度运动、用力、应激或其它条件。

边界条件提供关于三维实体模型320在它边界处的信息，例如流入边界322、流出边界324、血管壁边界326等，如图8中所示。血管壁边界326可包括模型320的主动脉、主要冠状动脉和/或其它冠状动脉或血管的实体边界。

各流入边界322或流出边界324可被赋予速度、流动速率、压力或其它血流特征的指定值或值域。或者，各流入边界322或流出边界324可通过使心脏模型与边界偶合、集中参数或分布(例如一维波传播)模型、另一类型的一维或二维模型或其它类型模型来赋值。特定边界条件可依据以下来确定：例如从所得患者特异性信息确定的流入边界322或流出边界324的几何形状，或其它测量参数，如心输出量、血压、在步骤240中计算的心肌质量，等。

i.确定降阶模型

与实体模型320连接的上游结构和下游结构可呈现为降阶模型，所述降阶模型呈现上游结构和下游结构。举例来说，图12-15展示根据一个示例性实施方案的一种方法的方面，所述方法用于从在流出边界324之一的三维患者特异性解剖学数据来准备集中参数模型。方法可与图2和图3中所示的方法分开进行以及在图2和图3中所示的方法之前进行。

图12展示主要冠状动脉之一或由此延伸的分支的实体模型320的一部分330，并且图13展示图12中所示的部分330的在步骤242中确定的中线的部分。

部分330可分成片段332。图14展示片段332的一个实施例，片段332可由部分330形成。片段332的长度选择可由用户和/或计算机系统进行。取决于例如片段332的几何形状，片段332的长度可变化。使用各种技术可来分割部分330。举例来说，可在一个或多个单独片段332中提供患病部分(例如截面相对较窄的部分)、病变和/或狭窄(血管异常变窄)。患病部分和狭窄可例如通过沿中线长度测量截面积以及计算局部最小截面积来鉴别。

片段332可由电路图来近似，所述电路图包括一个或多个(线性或非线性)电阻器334和/或其它电路元件(例如电容器、感应器等)。图15展示以一系列线性和非线性电阻器334替代的片段332的一个实施例。电阻器334的个别电阻可例如依据在相应片段332上的估计流量和/或压力来确定。

阻力可为恒定、线性或非线性的，例如取决于经过相应片段332的估计流动速率。对复杂几何形状(如狭窄)来说，阻力可随流动速率而变化。各种几何形状的阻力可依据计算分析(例如有限差分方法、有限体积方法、光谱方法、晶格波耳兹曼方法、基于颗粒的方法、水平集方法、等几何方法或有限元方法，或其它计算流体动力学(cfd)分析技术)来确定，并且在不同流量和压力条件下进行的计算分析的多重解可用以得出患者特异性阻力、血管特异性阻力和/或病变特异性阻力。结果可用以确定可被建模的任何片段的各种类型特征和几何形状的阻力。因此，如上所述得出患者特异性阻力、血管特异性阻力和/或病变特异性阻力可允许计算机系统识别并评价较复杂几何形状，如不对称狭窄、多种病变、在分叉和分支和扭曲血管等处的病变。

还可包括电容器，并且电容可例如依据相应片段的血管壁的弹性来确定。可包括感应器，并且感应系数可例如依据惯性作用来确定，所述惯性作用与流经相应片段的血体积的加速度或减速度有关。

电阻、电容、感应系数和与集中参数模型中所用的其它电组件相关的其它变量的个别值可依据来自许多患者的数据得出，并且类似血管几何形状可具有类似值。因此，经验模型可从大量患者特异性数据产生，从而创建对应于特定几何特征的值的文库，所述值文库可在未来分析中应用于类似患者。几何形状可在两个不同血管片段之间匹配以自动选择来自先前模拟的患者的片段332的值。

ii.示例性集中参数模型

或者，代替进行上文关于图12-15描述的步骤，可预设集中参数模型。举例来说，图16展示集中参数模型340、350、360的实施例，所述集中参数模型呈现在实体模型320的流入边界322和流出边界324处的上游结构和下游结构。a端位于流入边界322，并且a-m端和b端位于流出边界。

集中参数心脏模型340可用以确定在实体模型320的流入边界322的a端的边界条件。集中参数心脏模型340可用以呈现在充血条件下来自心脏的血流。集中参数心脏模型340包括各种参数(例如pla、rav、lav、rv-art、lv-art和e(t))，所述各种参数可依据关于患者的已知信息来确定，所述已知信息例如是关于主动脉压、患者的收缩血压和舒张血压(例如在步骤100中测定)、患者的心输出量(来自心脏的血流的体积，例如依据在步骤100中测定的患者的心搏量和心率来计算)，和/或用实验方式确定的常数。

集中参数冠状动脉模型350可用以确定在实体模型320的流出边界324的a-m端的边界条件，所述实体模型320位于主要冠状动脉和/或由此延伸的分支的下游末端。集中参数冠状动脉模型350可用以呈现在充血条件下经过a-m端离开建模血管的血流。集中参数冠状动脉模型350包括各种参数(例如ra、ca、ra-micro、cim和rv)，所述各种参数可依据关于患者的已知信息来确定，所述已知信息例如是计算心肌质量(例如在步骤240中确定)和a-m端终端阻抗(例如依据在步骤304中确定的血管在a-m端的截面积来确定)。

举例来说，计算心肌质量可用以估计经过多个流出边界324的基线(休息)平均冠状动脉流。这种关系可依据用实验方式得出的生理规律(例如图1的生理规律20)，所述生理规律将平均冠状动脉流q与心肌质量m(例如在步骤240中确定)关联为q∝qom^α，其中α是预设标度指数并且qo是预设常数。在基线(休息)条件下流出边界324的总冠状动脉流q以及患者的血压(例如在步骤100中测定)则可用以依据预设的用实验方式得出的等式来确定在流出边界324的总阻力r。

总阻力r可依据a-m端的各别截面积(例如在步骤304中测定)在a-m端分布。这种关系可依据用实验方式得出的生理规律(例如图1的生理规律20)，所述生理规律将a-m端的各别阻力关联为ri∝ri,odi^β，其中ri是在i-th出口的流动阻力，并且ri,o是预设常数，di是出口直径，并且β是预设幂律指数，例如在-3与-2之间，冠状动脉流为-2.7，大脑流为-2.9，等。经过个别a-m端的冠状动脉流以及在个别a-m端的平均压力(例如依据血管a-m端的个别截面积来确定，如在步骤304中测定)可用以确定集中参数冠状动脉模型350在相应a-m端的阻力的和(例如ra+ra-micro+rv)。其它参数(例如ra/ra-micro、ca、cim)可为用实验方式确定的常数。

弹性贮器(windkessel)模型360可用以确定在实体模型320的流出边界324的b端的边界条件，所述实体模型320位于朝着主动脉弓的主动脉的下游端。弹性贮器模型360可用以呈现在充血条件下经过b端离开建模主动脉的血流。弹性贮器模型360包括各种参数(例如rp、rd和c)，所述各种参数可依据关于患者的已知信息来确定，所述已知信息例如是上文关于集中参数心脏模型340描述的患者的心输出量、上文关于集中参数冠状动脉模型350描述的基线平均冠状动脉流、主动脉压(例如依据主动脉b端的截面积确定，如在步骤304中确定)，和/或用实验方式确定的常数。

边界条件，例如集中参数模型340、350、360(或其中包括的任一常数)或其它降阶模型可依据其它因素做调整。举例来说，如果患者的流量与血管尺寸比因在生理应激下血管扩张能力相对减小而较低，那么可调整(例如增大)阻力值。如果患者患有糖尿病、处于医药治疗、以往已经历心脏事件等，那么也可调整阻力值。

替代性集中参数模型或分布模型、一维网络模型可用以呈现实体模型320下游的冠状血管。使用mri、ct、pet或spect进行的心肌灌注成像可用以对所述模型的参数赋值。而且，替代性成像来源，例如磁共振血管造影术(mra)、回顾性电影门控或前瞻性电影门控计算机断层摄影血管造影术(cta)等可用以对所述模型的参数赋值。回顾性电影门控可与图像处理方法组合以在心动周期内获得心室腔容量变化，以便为集中参数心脏模型的参数赋值。

使用集中参数模型340、350、360或其它降阶一维或二维模型来简化患者解剖结构的一部分使得计算分析(例如下述图3的步骤402)较快速地进行，特别是当计算分析进行多次时，如当除未经处理的状态(例如图2和图3的步骤400)以外，还评价可能的治疗选择(例如图2的步骤500)，同时使最终结果维持高精度时。

在一个示例性实施方案中，边界条件的确定可由计算机系统依据用户输入(如在步骤100中获得的患者特异性生理学数据)来进行。

c.创建三维网

回头参考图3，可依据在步骤306中产生的实体模型320产生三维网(步骤312)。图17-19展示三维网380的一个实施例，所述三维网380是依据在步骤306中产生的实体模型320来准备。网380沿实体模型320的表面和在实体模型320的整个内部包括多个节点382(网点或格点)。如图18和图19中所示，网380可创建为具有四面体元(具有形成节点382的点)。或者，可使用具有其它形状的元，例如六面体或其它多面体、曲线元，等。在一个示例性实施方案中，节点382的数目可为百万，例如五百万至五千万。节点382的数目随着网380变细而增大。节点382数目越多，则可在模型320内更多点提供信息，但计算分析可能要运行更久，因为节点382的较大数目会增大待求解等式(例如图1中所示的等式30)的数目。在一个示例性实施方案中，网380的产生可由计算机系统在有或无用户输入(例如规定节点382的数目、元形状，等)下进行。

回头参考图3，可验证网380和所确定的边界条件(步骤314)。所述验证可由用户和/或由计算机系统进行。举例来说，用户和/或计算机系统也许能鉴别网380和/或边界条件的某些错误，例如，如果网380扭曲或不具有足够空间分辨率、如果边界条件不足以进行计算分析、如果在步骤310中确定的阻力不正确等，那么所述错误需要网380和/或边界条件被重做。如果是这样的话，那么可确定网380和/或边界条件为不可接受，并且可重复步骤304-314中一个或多个。如果网380和/或边界条件确定为可接受，那么所述方法可往下进行到下述步骤402。

另外，用户可检查所得患者特异性信息或其它测量参数，如心输出量、血压、身高、体重、在步骤240中计算的心肌质量得到正确输入和/或正确计算。

因此，图3中所示并且上文描述的步骤304-314可被视为图2的步骤300的子步骤。

v.进行计算分析以及输出结果

如上文关于图2中所示的步骤400所述，示例性方法可包括进行计算分析以及输出结果。在一个示例性实施方案中，步骤400可包括以下步骤。

a.进行计算分析

参考图3，计算分析可由计算机系统进行(步骤402)。在一个示例性实施方案中，步骤402可持续数分钟至数小时，取决于例如网380中节点382(图17-19)的数目，等。

分析涉及产生一系列等式，所述等式描述产生网380的模型320中的血流。如上所述，在示例性实施方案中，所需信息涉及在充血条件下经过模型320的血流的模拟。

分析还涉及使用计算机系统来使用数值方法来求解血流的三维等式。举例来说，数值方法可为已知方法，如有限差分方法、有限体积方法、光谱方法、晶格波耳兹曼方法、基于颗粒的方法、水平集方法、等几何方法或有限元方法，或其它计算流体动力学(cfd)数值技术。

使用这些数值方法，可将血液建模为牛顿流体、非牛顿流体或多相流体。在步骤100中测量的患者的血细胞比容或其它因素可用以确定用于并入分析中的血液粘度。可假设血管壁是刚性或顺从的。在后一情况下，壁动力学的等式，例如弹性动力学等式可连同血流等式一起求解。在步骤100中获得的时变三维成像数据可用作输入，以便为血管形状在心动周期内的变化建模。用于进行计算分析的一组示例性等式和步骤更详细地公开于以下中：例如标题为“methodforpredictivemodelingforplanningmedicalinterventionsandsimulatingphysiologicalconditions”的美国专利号6,236,878，和标题都为“patient-specifichemodynamicsofthecardiovascularsystem”的美国专利申请公布号2010/0241404和美国临时申请号61/210,401，所有这些专利都以引用的方式整体并入本文中。

使用所准备的模型和边界条件的计算分析可确定呈现三维实体模型320的网380的节点382中每一处的血液流量和压力。举例来说，计算分析的结果可包括节点382中每一处的各种参数的值，所述各种参数如但不限于各种血流特征或参数，如血液流动速度、压力、流动速率，或计算参数，如cffr，如下所述。参数也可在三维实体模型320上内插。因此，计算分析的结果可向用户提供通常可能以侵入方式确定的信息。

回头参考图3，可验证计算分析的结果(步骤404)。所述验证可由用户和/或由计算机系统进行。举例来说，用户和/或计算机系统也许能鉴别出某些错误，例如，如果节点382的数目不足导致信息不足、如果节点382的数目过量导致分析费时太久，那么结果是需要网380和/或边界条件重做或修正。

如果计算分析的结果在步骤404中确定为不可接受，那么用户和/或计算机系统可确定例如是否以及如何修正或改进在步骤306中产生的实体模型320和/或在步骤312中产生的网380、是否以及如何修正在步骤310中确定的边界条件，或是否对用于计算分析的任一输入进行其它修正。接着，可依据所确定的修正或改进来重复一个或多个上述步骤，例如步骤306-314、402和404。

b.显示血压、流量和cffr的结果

回头参考图3，如果计算分析的结果在步骤404中确定为可接受，那么计算机系统可输出计算分析的某些结果。举例来说，计算机系统可显示依据计算分析结果产生的图像，如上文关于图1描述的图像，如模拟血压模型50、模拟血流模型52和/或cffr模型54。如上所述，例如因为在步骤310中确定的边界条件是关于充血条件所确定，所以这些图像指示在模拟充血条件下的模拟血压、血流和cffr。

模拟血压模型50(图1)展示在整个患者解剖结构中的局部血压(例如以毫米汞柱或mmhg表示)，如由图17-19的在模拟充血条件下的网380表示。计算分析可确定网380的各节点382处的局部血压，并且模拟血压模型50可将相应颜色、色调或其它视觉指示物赋予各别压力，以使模拟血压模型50可以视觉方式指示在整个模型50中的压力变化，而不用必须规定各节点382的个别值。举例来说，图1中所示的模拟血压模型50展示，对这个具体患者来说，在模拟充血条件下，主动脉中的压力一般可为均一并且较高的(如由较暗色调指示)，并且压力随着血液向下游流到主要冠状动脉中并且流到分支中而逐渐并不断地降低(如由朝着分支的下游端，色调逐渐并连续变亮所示)。模拟血压模型50可伴有指示特定血压数值的标度，如图1中所示。

在一个示例性实施方案中，模拟血压模型50可以颜色形式提供，并且可使用色谱来指示整个模型50中的压力变化。依序从最高压力到最低压力，色谱可包括红色、橙色、黄色、绿色、蓝色、靛蓝和紫色。举例来说，上限(红色)可指示约110mmhg或更大(或80mmhg、90mmhg、100mmhg等)，并且下限(紫色)可指示约50mmhg或更小(或20mmhg、30mmhg、40mmhg等)，其中绿色指示约80mmhg(或约在上下限之间一半的其它值)。因此，对一些患者来说，模拟血压模型50可将大部分或所有主动脉展示为红色或朝着光谱高端的其它颜色，并且朝着冠状动脉和由此延伸的分支的远端，颜色可在整个光谱上逐渐变化(例如朝着光谱下端(下至紫色))。对一个具体患者来说，冠状动脉的远端可具有不同颜色，例如从红色到紫色之间的任何颜色，这取决于对各别远端测定的局部血压。

模拟血流模型52(图1)展示在整个患者解剖结构中的局部血液速度(例如以厘米/秒或cm/s表示)，这由图17-19的在模拟充血条件下的网380表示。计算分析可确定网380的各节点382处的局部血液速度，并且模拟血流模型52可将相应颜色、色调或其它视觉指示物赋予各别速度，以使模拟血流模型52可以视觉方式指示在整个模型52中的速度变化，而不用必须规定各节点382的个别值。举例来说，图1中所示的模拟血流模型52展示，对这个具体患者来说，在模拟充血条件下，在主要冠状动脉和分支的某些区域中速度一般较高(如由图1中区域53中的较暗色调指示)。模拟血流模型52可伴有指示特定血液速度数值的标度，如图1中所示。

在一个示例性实施方案中，模拟血流模型52可以颜色形式提供，并且可使用色谱来指示整个模型52中的速度变化。依序从最高速度到最低速度，色谱可包括红色、橙色、黄色、绿色、蓝色、靛蓝和紫色。举例来说，上限(红色)可指示约100(或150)cm/s或更大，并且下限(紫色)可指示约0cm/s，其中绿色指示约50cm/s(或约在上下限之间一半的其它值)。因此，对一些患者来说，模拟血流模型52可将大部分或所有主动脉展示为朝着光谱下端的颜色(例如绿色至紫色)的混合色，并且在测定血液速度增大的某些位置，颜色可在整个光谱上逐渐变化(例如朝着光谱高端(上至红色))。

cffr模型54(图1)展示在整个患者解剖结构中的局部cffr，这由图17-19的在模拟充血条件下的网380表示。如上所述，cffr可计算为通过计算分析(例如模拟血压模型50中所示)确定的在具体节点382的局部血压除以在主动脉中(例如在流入边界322(图8))的血压的比率。计算分析可确定网380的各节点382处的cffr，并且cffr模型54可将相应颜色、色调或其它视觉指示物赋予各别cffr值，以使cffr模型54可以视觉方式指示在整个模型54中的cffr变化，而不用必须规定各节点382的个别值。举例来说，图1中所示的cffr模型54展示，对这个具体患者来说，在模拟充血条件下，主动脉中的cffr一般可为均一的并且为约1.0，并且随着血液向下游流到主要冠状动脉中并且流到分支中，cffr逐渐并不断降低。cffr模型54也可指示在整个cffr模型54中某些点处的cffr值，如图1中所示。cffr模型54可伴有指示特定cffr数值的标度，如图1中所示。

在一个示例性实施方案中，cffr模型54可以颜色形式提供，并且可使用色谱来指示整个模型54中的压力变化。依序从最低cffr(指示功能上有意义的病变)到最高cffr，色谱可包括红色、橙色、黄色、绿色、蓝色、靛蓝和紫色。举例来说，上限(紫色)可指示cffr为1.0，并且下限(红色)可指示约0.7(或0.75或0.8)或更小，其中绿色指示约0.85(或约在上下限之间一半的其它值)。举例来说，下限可依据用于确定cffr测量值是否指示功能上有意义的病变或可能需要介入的其它特征的下限(例如0.7、0.75或0.8)来确定。因此，对一些患者来说，cffr模型54可将大部分或所有主动脉展示为紫色或朝着光谱高端的其它颜色，并且朝着冠状动脉和由此延伸的分支的远端，颜色可在整个光谱上逐渐变化(例如朝着光谱高端(上至红色到紫色之间的任何颜色))。对一个具体患者来说，冠状动脉的远端可具有不同颜色，例如从红色到紫色之间的任何颜色，这取决于对各别远端测定的局部cffr值。

在确定cffr已降到下限以下的程度之后，可评估动脉或分支以定位功能上有意义的病变，所述下限是用于确定功能上有意义的病变或可能需要介入的其它特征的存在。计算机系统或用户可依据动脉或分支的几何形状(例如使用cffr模型54)来定位功能上有意义的病变。举例来说，可如下定位功能上有意义的病变：找到在具有局部最小cffr值的cffr模型54的位置附近(例如上游)定位的变窄或狭窄。计算机系统可向用户指示或显示包括功能上有意义的病变的cffr模型54(或其它模型)的部分。

也可依据计算分析结果产生其它图像。举例来说，计算机系统可提供关于具体主要冠状动脉的其它信息，例如图20-22中所示。举例来说，如果具体冠状动脉包括最低cffr，那么冠状动脉可由计算机系统选择。或者，用户可选择具体冠状动脉。

图20展示患者解剖结构的模型，其包括计算分析结果，其中模型上的某些点由个别参考标签(例如lm、lad1、lad2、lad3等)鉴别。在图21中所示的示例性实施方案中，在lad动脉中提供点，所述lad动脉是在模拟充血条件下具有对这个具体患者来说最低的cffr的主要冠状动脉。

图21和图22展示在一些或所有这些点(例如lm、lad1、lad2、lad3，等)处和/或在模型上某些其它位置(例如在主动脉中，等)处某些变量随时间变化的图。图21是主动脉中和在图20中指示的点lad1、lad2和lad3处的压力(例如以毫米汞柱或mmhg表示)随时间变化的图。所述图的顶部曲线指示主动脉中的压力，从上往下第二条曲线指示在点lad1的压力，从上往下第三条曲线指示在点lad2的压力，并且底部曲线指示在点lad3的压力。图22是在图20中指示的点lm、lad1、lad2和lad3处的流量(例如以立方厘米/秒或cc/s表示)随时间变化的图。另外，可提供其它图，如在一些或所有这些点和/或在其它点处剪应力随时间变化的图。所述图的顶部曲线指示在点lm的流量，从上往下第二条曲线指示在点lad1的流量，从上往下第三条曲线指示在点lad2的流量，并且底部曲线指示在点lad3的流量。也可提供展示这些变量沿具体主要冠状动脉和/或由此延伸的分支的长度的变化的图：例如血压、流量、速度或cffr。

任选地，上述各种图和其它结果可在报告中定稿(步骤406)。举例来说，上述图像和其它信息可插入具有设置模板的文件中。模板可预设并且对多个患者来说是通用的，并且可用于向医师和/或患者报道计算分析结果。在计算分析完成之后，文件或报告可由计算机系统自动完成。

举例来说，定稿报告可包括图23中所示的信息。图23包括图1的cffr模型54，并且还包括概述信息，如主要冠状动脉和由此延伸的分支中每一个中的最低cffr值。举例来说，图23指示lad动脉中的最低cffr值是0.66，lcx动脉中的最低cffr值是0.72，rca动脉中的最低cffr值是0.80。其它概述信息可包括患者的名字、患者的年龄、患者的血压(bp)(例如在步骤100中获得)、患者的心率(hr)(例如在步骤100中获得)，等。定稿报告还可包括如上述产生的图像和其它信息的版本，医师或其它用户可得到所述版本以确定其它信息。可使由计算机系统产生的图像有一定格式以使医师或其它用户在任一点上定位光标，以确定任一上述变量的值，所述变量例如是在所述点的血压、速度、流量、cffr等。

定稿报告可被传输到医师和/或患者。定稿报告可使用任何已知通信方法(例如无线网络或有线网络，通过邮件等)传输。或者，可通知医师和/或患者定稿报告可下载或提取得到。接着，医师和/或患者可登录到基于网络的服务中经由安全通信线来下载定稿报告。

c.验证结果

回头参考图3，可独立地验证计算分析的结果(步骤408)。举例来说，用户和/或计算机系统也许能鉴别出计算分析结果(例如在步骤406中产生的图像和其它信息)的某些错误，所述错误需要任一上述步骤重做。如果鉴别出所述错误，那么计算分析结果可确定为不可接受，并且可重复某些步骤，例如步骤100、200、300、400、子步骤102、202-208、240-260、304-314和402-408等。

因此，图3中所示并且上文描述的步骤402-408可被视为图2的步骤400的子步骤。

用于验证计算分析结果的另一方法可包括使用另一方法测量来自患者的结果中包括的任一变量，例如血压、速度、流量、cffr等。在一个示例性实施方案中，可测量(例如以侵入方式)变量，接着与通过计算分析确定的结果相比。举例来说，ffr可例如使用插入如上所述的患者中的压力线来测定，测定点在由实体模型320和网380表示的患者解剖结构内的一个或多个点。可将在一个位置测量的ffr与在同一位置的cffr相比，并且可在多个位置进行比较。任选地，计算分析和/或边界条件可依据比较做调整。

d.用于提供冠状动脉血流信息的系统和方法的另一实施方案

图24中展示一种方法600的另一实施方案，所述方法600用于提供有关特定患者的冠状动脉血流的各种信息。方法600可在上述计算机系统中实施，所述计算机系统例如是用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。方法600可使用一个或多个输入610来进行，并且可包括依据输入610来产生一个或多个模型620、依据输入610和/或模型620对一个或多个条件630赋值，以及依据模型620和条件630得出一个或多个解640。

输入610可包括患者的主动脉、冠状动脉(和由此延伸的分支)和心脏的医学成像数据611，如ccta数据(例如在图2的步骤100中获得)。输入610还可包括患者臂血压的测量值612和/或其它测量值(例如在图2的步骤100中获得)。测量值612可以非侵入方式获得。输入610可用以产生模型620和/或确定下述条件630。

如上所述，一个或多个模型620可依据输入610来产生。举例来说，方法600可包括依据成像数据611来产生患者解剖结构(例如主动脉、冠状动脉和由此延伸的分支)的一个或多个患者特异性三维几何模型(步骤621)。举例来说，几何模型可为在图3的步骤306中产生的图8的实体模型320，和/或在图3的步骤312中产生的图17-19的网380。

回头参考图24，方法600还可包括产生一个或多个基于物理学的血流模型(步骤622)。血流模型可包括涉及以下的模型：经过在步骤621中产生的患者特异性几何模型的血流、心脏和主动脉循环、远端冠状循环等。血流模型可涉及至少一个与患者经过建模的解剖结构(例如血液流动速度、压力、流动速率、ffr等)相关的血流特征。血流模型可赋值为在三维几何模型的流入边界322和流出边界324的边界条件。血流模型可包括降阶模型或其它边界条件(上文关于图3的步骤310所述)，例如集中参数心脏模型340、集中参数冠状动脉模型350、弹性贮器模型360，等。

如上所述，一个或多个条件630可依据输入610和/或模型620来确定。条件630包括针对在步骤622(和图3的步骤310)中确定的边界条件所计算的参数。举例来说，方法600可包括如下确定条件：依据成像数据611(例如在图3的步骤240中确定)来计算患者特异性心室质量或心肌质量(步骤631)。

方法600可包括如下确定条件：使用在步骤631中计算的心室质量或心肌质量、休息时冠状动脉流依据关系q＝qom^α(其中α是预设标度指数，m是心室质量或心肌质量，并且qo是预设常数)进行计算(例如上文在图3的步骤310中关于确定集中参数模型所述)(步骤632)。或者，所述关系可具有形式q∝qom^α，如上文在图3的步骤310中关于确定集中参数模型所述。

方法600还可包括如下确定条件：使用在步骤632中计算的所得冠状动脉流量和患者的测量血压612、总的休息时冠状动脉阻力进行计算(例如上文在图3的步骤310中关于确定集中参数模型所述)(步骤633)。

方法600还可包括如下确定条件：使用在步骤633和模型620中计算的总的休息时冠状动脉阻力、个别冠状动脉(和由此延伸的分支)的个别阻力进行计算(步骤634)。举例来说，如上文关于图3的步骤310所述，依据个别冠状动脉和分支的远端的尺寸(例如从步骤621中产生的几何模型确定)，并且依据关系r＝rod^β(其中r是在具体远端的流动阻力，并且ro是预设常数，d是尺寸(例如所述远端的直径)，并且β是预设幂律指数，如上文关于确定图3的步骤310中的集中参数模型所述)，在步骤633中计算的总的休息时冠状动脉阻力可分布到个别冠状动脉和分支。

回头参考图24，方法600可包括依据患者的一个或多个身体条件来调整边界条件(步骤635)。举例来说，在步骤631-634中确定的参数可依据解640是否意图模拟休息、不同程度的充血、不同程度的运动或用力、不同药物治疗等来修改。依据输入610、模型620和条件630，计算分析可例如如上文关于图3的步骤402所述来进行以确定解640，所述解640包括关于患者的冠状动脉血流在步骤635中选择的身体条件下的信息(步骤641)。现将描述可从解640提供的信息的实例。

组合的患者特异性解剖学(几何)模型和生理学(基于物理学的)模型可用以确定不同药物治疗或生活方式变化(例如停止吸烟、改变膳食或增加身体活动)对冠状动脉血流的作用，所述药物治疗或生活方式变化会改变心率、心搏量、血压或冠状动脉微循环功能。所述信息可用以最优化医学疗法或避免可能危险的药物治疗后果。组合模型还可用以确定替代形式和/或不同程度的暴露于可能的外来力(例如当踢足球时、在太空飞行期间、当斯库巴潜水(scubadiving)时、在坐飞机飞行期间，等)的身体活动或风险对冠状动脉血流的作用。所述信息可用以鉴别对特定患者来说可能安全并有效的身体活动的类型和程度。组合模型还可用以预测经皮冠状动脉介入术对冠状动脉血流的可能益处以选择最优的介入策略，和/或预测冠状动脉旁路移植术对冠状动脉血流的可能益处以选择最优的手术策略。

组合模型还可用以说明动脉疾病负担增大对冠状动脉血流的可能的有害作用，并且使用机械或现象学疾病进展模型或经验数据来预测进展中疾病何时可对流向心脏肌肉的血流引起损害。所述信息可实现确定“保证期(warrantyperiod)”，其中使用非侵入性成像观察到最初没有热力学有意义的疾病的患者可能不被预期需要医药疗法、介入疗法或手术疗法，或者，如果继续有不利因素，那么可确定出现进展的比率。

组合模型还可用以说明由冠状动脉疾病负担减小引起的对冠状动脉血流的可能的有益作用，并且使用机械或现象学疾病进展模型或经验数据来预测疾病消退何时可引起经过冠状动脉到心脏肌肉的血流增加。所述信息可用以指导医学管理程序，所述医学管理程序包括但不限于改变膳食、增加身体活动、开出士他汀(statin)或其它医药的处方，等。

vi.提供患者特异性治疗计划

如上文关于图2中所示的步骤500所述，示例性方法可包括提供患者特异性治疗计划。在一个示例性实施方案中，步骤500可包括以下步骤。尽管图3不展示以下步骤，应了解这些步骤可与图3中所示的步骤(例如在步骤406或408之后)结合进行。

如上所述，图1和图23中所示的cffr模型54指示处于未治疗状态并且在模拟充血条件下的整个患者解剖结构中的cffr值，所述解剖结构由图17-19的网380表示。医师可使用这个信息来指定对患者的治疗，如增加运动、改变膳食、开医药处方、对建模解剖结构的任何部分或心脏的其它部分进行手术(例如冠状动脉旁路移植术、插入一个或多个冠状动脉血管内支架，等)，等。

为确定指定哪种治疗，计算机系统可用以预测如何依据所述治疗来改变从计算分析确定的信息。举例来说，某些治疗(如插入血管内支架或其它手术)可引起建模解剖结构的几何形状改变。因此，在一个示例性实施方案中，可修正在步骤306中产生的实体模型320以指示一个或多个插有血管内支架的内腔的加宽。

举例来说，图1和图23中所示的cffr模型54指示lad动脉中的最低cffr值是0.66，lcx动脉中的最低cffr值是0.72，rca动脉中的最低cffr值是0.80。如果cffr值例如小于0.75，那么可建议进行治疗。因此，计算机系统可建议用户修正实体模型320以指示lad动脉和lcx动脉的加宽，以模拟在这些冠状动脉中插入血管内支架。用户可被提示选择对应于模拟血管内支架的位置和尺寸的加宽(例如长度和直径)的位置和量。或者，加宽的位置和量可由计算机系统依据各种因素(如cffr值小于0.75的节点位置、血管显著变窄的位置、常规血管内支架的尺寸，等)自动确定。

图25展示依据实体模型确定的经过修改的cffr模型510的一个实施例，上述实体模型是通过加宽lad动脉在位置512的一部分和lcx动脉在位置514的一部分来产生。在一个示例性实施方案中，任一上述步骤，例如步骤310-314和402-408都可使用经过修改的实体模型重复。在步骤406中，定稿报告可包括有关未经过治疗的患者(例如没有血管内支架)的信息，如图23中所示的信息，和有关患者的模拟治疗的信息，如图25和图26中所示的信息。

图25包括经过修改的cffr模型510并且还包括概述信息，如与建议治疗相关的经过修改的实体模型的主要冠状动脉和由此延伸的分支中的最低cffr值。举例来说，图25指示lad动脉(和它的下游分支)中的最低cffr值是0.78，lcx动脉(和它的下游分支)中的最低cffr值是0.78，rca动脉(和它的下游分支)中的最低cffr值是0.79。因此，未经过治疗的患者(没有血管内支架)的cffr模型54和建议治疗(有血管内支架插入)的cffr模型510的比较指示建议治疗可使lad动脉中的最小cffr从0.66增大到0.78，并且使lcx动脉中的最小cffr从0.72增大到0.76，同时使rca动脉中的最小cffr将有从0.80到0.79的最小降低。

图26展示经过修改的模拟血流模型520的一个实施例，所述模型是如上所述在lad动脉在位置512的加宽部分和lcx动脉在位置514的加宽部分之后确定。图26还包括概述信息，如与建议治疗相关的经过修改的实体模型的主要冠状动脉和由此延伸的分支中各种位置的血流值。举例来说，图26指示未经过治疗的患者(没有血管内支架)和经过治疗的患者(有血管内支架插入)的lad动脉中四个位置lad1、lad2、lad3和lad4以及lcx动脉中两个位置lcx1和lcx2的血流值。图26还指示未经过治疗的状态和经过治疗的状态之间的血流值变化百分比。因此，未经过治疗的患者的模拟血流模型52和建议治疗的模拟血流模型520的比较指示，在所有位置lad1-lad4、lcx1和lcx2，取决于位置，建议治疗可将经过lad动脉和lcx动脉的流增大9％至19％。

也可在未经过治疗的状态和经过治疗的状态(如冠状动脉血压)之间比较其它信息。依据这个信息，医师可与患者讨论是否继续建议治疗选择。

其它治疗选择也可涉及以不同方式修改实体模型320。举例来说，冠状动脉旁路移植术可涉及在实体模型320中创建新的内腔或通路，并且去除病变还可涉及加宽内腔或通道。其它治疗选择可能不涉及修改实体模型320。举例来说，例如由于血管收缩、扩张、心率降低等，因此增加运动或用力、改变膳食或其它生活方式变化、医药处方等可涉及改变在步骤310中确定的边界条件。举例来说，患者的心率、心输出量、心搏量、血压、冠状动脉微循环功能、集中参数模型的配置等可取决于指定医药、所采用的运动(或其它用力)的类型和频率、所采用的生活方式变化的类型(例如停止吸烟、改变膳食等)，从而以不同方式影响在步骤310中确定的边界条件。

在一个示例性实施方案中，经过修改的边界条件可用实验方式使用来自许多患者的数据来确定，并且类似治疗选择可能需要以类似方式修改边界条件。经验模型可从大量患者特异性数据产生，从而创建边界条件或用于计算边界条件的函数的文库，所述文库对应于可在未来分析中应用于类似患者的特定治疗选择。

在修改边界条件之后，上述步骤，例如步骤312、314和402-408可使用经过修改的边界条件来重复，并且在步骤406中，定稿报告可包括有关未经过治疗的患者的信息，如图23中所示的信息，以及有关患者的模拟治疗的信息，如图25和图26中所示的信息。

或者，可向医师、患者或其它用户提供允许与三维模型(例如图8的实体模型320)互动的用户界面。模型320可分成用户可选片段，所述片段可由用户编辑来反映一个或多个治疗选择。举例来说，用户可选择具有狭窄(或梗塞，例如急性阻塞)的部分，并且调整所述部分来去除狭窄，用户可将一个部分加入模型320中以用作旁路等。用户还可被提示指定其它治疗选择和/或生理学参数，这可改变上文确定的边界条件，例如改变心输出量、心率、心搏量、血压、运动或用力程度、充血程度、药物治疗，等。在一个替代性实施方案中，计算机系统可确定或建议一个治疗选择。

用户界面可允许与三维模型320互动以允许用户模拟狭窄(或梗塞，例如急性梗塞)。举例来说，用户可选择用于包括狭窄的部分，并且计算机系统可用以预测从计算分析确定的信息将如何依据狭窄的添加来变化。因此，本文所述的方法可用以预测使动脉梗塞的影响。

用户界面还可允许与三维模型320互动以模拟受损动脉或动脉去除，所述动脉去除可例如在某些手术程序中发生，如当去除癌性肿瘤时发生。模型还可得到修改以模拟防止血流经过某些动脉的作用，以便预测并行通路向患者提供足够血流的可能。

a.使用降阶模型来比较不同治疗选择

在一个示例性实施方案中，计算机系统可允许用户用降阶模型替代三维实体模型320或网380来较快速地模拟各种治疗选择。图27展示根据一个示例性实施方案，有关使用降阶模型来模拟各种治疗选择的方法700的示意图。方法700可在上述计算机系统中实施。

一个或多个呈现血流或其它参数的患者特异性模拟血流模型可从上述计算分析输出(步骤701)。举例来说，模拟血流模型可包括图1的模拟血压模型50、图1的模拟血流模型52、图1的cffr模型54等，这些模型是使用上文描述并且在图2和图3中所示的方法提供。如上所述，模拟血流模型可包括患者的解剖结构的三维几何模型。

功能信息可从模拟血流模型提取以指定降阶模型的条件(步骤702)。举例来说，功能信息可包括使用上述计算分析确定的血压、流量或速度信息。

可提供降阶(例如零维或一维)模型来替代用以产生患者特异性模拟血流模型(在步骤701中产生)的三维实体模型320，并且降阶模型可用以确定关于患者的冠状动脉血流的信息(步骤703)。举例来说，降阶模型可为如上文关于图3的步骤310所述产生的集中参数模型。因此，集中参数模型是患者的解剖结构的简化模型，其可用以确定关于患者的冠状动脉血流的信息，而不用必须对与图17-19的网380相关的较复杂等式系统求解。

从对步骤703中的降阶模型求解所确定的信息接着可被绘制或外推到患者的解剖结构的三维实体模型(例如实体模型320)中(步骤704)，并且必要时用户可对降阶模型进行变动以模拟各种治疗选择和/或改为患者的生理学参数，所述生理学参数可由用户选择(步骤705)。可选生理学参数可包括心输出量、运动或用力程度、充血程度、药物治疗类型，等。可选治疗选择可包括去除狭窄、增加旁路，等。

接着，降阶模型可依据治疗选择和/或用户选择的生理学参数加以修改，并且经过修改的降阶模型可用以确定关于患者的冠状动脉血流的信息(步骤703)，所述信息与所选治疗选择和/或生理学参数相关。从对步骤703中的降阶模型求解所确定的信息接着可被绘制或外推到患者的解剖结构的三维实体模型320中，以预测所选治疗选择和/或生理学参数对患者解剖结构的冠状动脉血流的影响(步骤704)。

对各种不同治疗选择和/或生理学参数来说，可重复步骤703-705以比较所预测的各种治疗选择彼此之间以及对关于未经过治疗的患者的冠状动脉血流的信息的影响。因此，各种治疗选择和/或生理学参数的预测结果可针对彼此并且针对关于未经过治疗的患者的信息来评价，而不用必须使用三维网380来重新进行较复杂分析。代之以，可使用降阶模型，降阶模型可允许用户更容易和快速地分析并比较不同治疗选择和/或生理学参数。

图28展示根据一个示例性实施方案，用于使用降阶模型来模拟各种治疗选择的示例性方法的其它方面。方法700可在上述计算机系统中实施。

如上文关于图3的步骤306所述，患者特异性几何模型可依据患者的成像数据来产生(步骤711)。举例来说，成像数据可包括在图2的步骤100中获得的ccta数据，并且几何模型可为在图3的步骤306中产生的图8的实体模型320，和/或在图3的步骤312中产生的图17-19的网380。

计算分析可使用患者特异性三维几何模型来进行，例如上文关于图3的步骤402所述，以确定关于患者的冠状动脉血流的信息(步骤712)。计算分析可输出一种或多种三维患者特异性模拟血流模型，所述模型呈现血流或其它参数，例如图1的模拟血压模型50、图1的模拟血流模型52、图1的cffr模型54，等。

可依据模型的解剖特征(步骤713)分割模拟血流模型(例如上文关于图14所述)。举例来说，从主要冠状动脉延伸的分支可提供于个别片段中(步骤714)，具有狭窄或患病区域的部分可提供于个别片段中(步骤716)，并且在分支和具有狭窄或患病区域的部分之间的部分可提供于个别片段中(步骤715)。可在分割模拟血流模型时提供不同分辨度，以使各血管可包括多个短的离散片段或较长片段，例如包括整根血管。而且，各种技术可提供用于分割模拟血流模型，包括产生中线以及依据所产生的中线制作切面图，或检测分支点以及依据检测到的中线制作切面图。患病部分和狭窄可例如通过沿中线长度测量截面积以及计算局部最小截面积来鉴别。步骤711-716可被视为图27的步骤701的子步骤。

片段可用集中参数模型的组件替代，所述组件如为电阻器、电容器、感应器等，如上文关于图15所述。电阻、电容、感应系数和与集中参数模型中所用的其它电组件相关的其它变量的个别值可来源于步骤712中提供的模拟血流模型。举例来说，对在分支和具有狭窄或患病区域的部分之间的分支和部分来说，来源于模拟血流模型的信息可用以将线性阻力赋予相应片段(步骤717)。对具有复杂几何形状的部分(如狭窄或患病区域)来说，阻力可随流动速率而变化。因此，多种计算分析可用以获得各种流量和压力条件的模拟血流模型以得出这些复杂几何形状的患者特异性、血管特异性和病变特异性阻力函数，如上文关于图15所述。因此，对具有狭窄或患病区域的部分来说，来源于这多种计算分析的信息或来源于先前数据的模型可用以将非线性流量依赖性阻力赋予相应片段(步骤718)。步骤717和步骤718可被视为图27的步骤702的子步骤。

使用在步骤717和步骤718中确定的阻力，可产生降阶(例如零维或一维)模型(步骤719)。举例来说，降阶模型可为如上文关于图3的步骤310所述产生的集中参数模型。因此，集中参数模型是患者的解剖结构的简化模型，其可用以确定关于患者的冠状动脉血流的信息，而不用必须对与图17-19的网380相关的较复杂等式系统求解。

可提供一种用户界面，其允许用户与在步骤719中创建的降阶模型互动(步骤720)。举例来说，用户可选择并编辑降阶模型的不同片段以模拟不同治疗选择和/或可编辑各种生理学参数。举例来说，介入，如插入血管内支架以修复患病区域可通过减小插有血管内支架的片段的阻力来建模。形成旁路可通过增加具有平行于患病片段的低阻力的片段来建模。

经过修改的降阶模型可被求解以确定关于在治疗和/或在步骤720中选择的生理学参数变化下患者的冠状动脉血流的信息(步骤721)。在步骤721中确定的各片段中流量和压力的解值接着可与在步骤712中确定的三维解相比，并且任何差异都可通过调整片段的阻力函数(例如在步骤717和步骤718中确定)以及求解降阶模型(例如步骤721)直到解匹配为止而减到最小。因此，可创建降阶模型，接着用简化等式组对所述降阶模型求解，所述简化等式组允许相对快速的计算(例如与完整三维模型相比)，并且可用以求出流动速率和压力，这结果可密切接近完整三维计算解的结果。降阶模型允许相对快速迭代以为各种不同治疗选择建模。

从对步骤721中的降阶模型求解所确定的信息接着可被绘制或外推到患者的解剖结构的三维实体模型(例如实体模型320)(步骤722)。步骤719-722可能类似于图27的步骤703-705，并且必要时可由用户重复以模拟治疗选择和/或生理学参数的不同组合。

或者，不从三维模型沿片段计算阻力(例如上文关于步骤717和步骤718所述)，沿中线相隔一定间隔的流量和压力可被指定到集中参数模型或一维模型中。有效阻力或损失系数可在边界条件和指定流量和压力的限制下求出。

而且，在个别片段上的流动速率和压力梯度可用以使用来源于降阶模型的解(例如上关于步骤721所述)来计算心外膜冠状动脉阻力。心外膜冠状动脉阻力可计算为心外膜冠状动脉(从医学成像数据重建的患者特异性模型中包括的冠状动脉和由此延伸的分支的部分)的等效阻力。这可具有临床意义，所述临床意义是在于解释为什么患有冠状动脉扩散性动脉粥样硬化的患者可显示局部缺血的症状(血液供给的限制)。而且，在模拟的药理学诱发性充血或不同运动强度的条件下每单位心肌组织体积(或质量)的流量和/或每单位心脏作功的流量可使用来自降阶模型的数据来计算。

因此，三维血流建模的精度可与一维和集中参数模型技术中固有的计算简单性和相对速度组合。三维计算方法可用以用数字方式得出患者特异性一维模型或集中参数模型，所述模型嵌有在正常片段、狭窄、接合点(junction)和其它解剖特征上压力损失的以数字方式得出的经验模型。可向患有心血管疾病的患者提供改善的诊断，并且医学治疗、介入治疗和手术治疗的计划可较快地进行。

而且，三维计算流体动力学技术的精度可与血流的集中参数和一维模型的计算简单性和进行能力组合。三维几何和生理学模型可自动分解成降阶一维模型或集中参数模型。三维模型可用以计算经过正常片段、狭窄和/或分支的血流的线性或非线性血液动力学作用，并且可用以设置经验模型的参数。一维模型或集中参数模型可更有效并快速地求出患者特异性模型中的血液流量和压力，并且显示集中参数或一维解的结果。

降阶患者特异性解剖学和生理学模型可用以确定不同药物治疗或生活方式变化(例如停止吸烟、改变膳食或增加身体活动)对冠状动脉血流的作用，所述药物治疗或生活方式变化会改变心率、心搏量、血压或冠状动脉微循环功能。所述信息可用以最优化医学疗法或避免可能危险的药物治疗后果。降阶模型还可用以确定替代形式和/或不同程度的暴露于可能的外来力(例如当踢足球时、在太空飞行期间、当斯库巴潜水时、在坐飞机飞行期间，等)的身体活动或风险对冠状动脉血流的作用。所述信息可用以鉴别对特定患者来说可能安全并有效的身体活动的类型和程度。降阶模型还可用以预测经皮冠状动脉介入术对冠状动脉血流的可能益处以选择最优的介入策略，和/或预测冠状动脉旁路移植术对冠状动脉血流的可能益处以选择最优的手术策略。

降阶模型还可用以说明动脉疾病负担增大对冠状动脉血流的可能的有害作用，并且使用机械或现象学疾病进展模型或经验数据来预测进展中疾病何时可对流向心脏肌肉的血流引起损害。所述信息可实现确定“保证期”，其中使用非侵入性成像观察到最初没有热力学有意义的疾病的患者可能不被预期需要医药疗法、介入疗法或手术疗法，或者，如果继续有不利因素，那么可确定出现进展的比率。

降阶模型还可用以说明由冠状动脉疾病负担减小引起的对冠状动脉血流的可能的有益作用，并且使用机械或现象学疾病进展模型或经验数据来预测疾病消退何时可引起经过冠状动脉到心脏肌肉的血流增加。所述信息可用以指导医学管理程序，所述医学管理程序包括但不限于改变膳食、增加身体活动、开出士他汀或其它医药的处方，等。

降阶模型还可并入血管造影系统中以允许医师在心导管插入实验室中检查患者的同时现场计算治疗选择。模型可记录到与血管造影显示器相同的取向，从而允许冠状动脉的现场血管造影图与模拟血流解具有并行或重叠的结果。当在程序期间进行观察时，医师可计划并改变治疗计划，从而允许在作出医学决定之前有相对快速的反馈。医师可以侵入方式获得压力、ffr或血流测量值，并且在进行预测性模拟之前，测量值可用以进一步改进模型。

降阶模型还可并入医学成像系统或工作站中。如果来源于先前患者特异性模拟结果的文库，那么降阶模型可与几何分割算法结合使用以在完成成像扫描之后相对快速地求出血流信息。

降阶模型还可用以对新型医学疗法对大的患者群体的有效性或治疗选择对大的患者群体的成本/益处进行建模。多种患者特异性集中参数模型(例如数百、数千或更多)的数据库可提供在相对较短时间长度内求解的模型。相对快速的迭代和最优化可提供用于药物、疗法或临床试验模拟或设计。调适所述模型来呈现治疗、患者的药物反应或手术介入可允许估计待获得的有效性，而无需进行可能高成本并且可能冒险的大规模临床试验。

vii.其它结果

a.评估心肌灌注

可计算其它结果。举例来说，计算分析可提供定量心肌灌注(经过心肌的血流)的结果。如由于局部缺血(血液供给的限制)、疤痕或其它心脏问题，因此定量心肌灌注可有助于鉴别心肌血流减少的区域。

图29展示根据一个示例性实施方案，有关方法800的示意图，所述方法800用于提供有关特定患者的心肌灌注的各种信息。方法800可在上述计算机系统中实施，所述计算机系统例如是用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。

方法800可使用一个或多个输入802来进行。输入802可包括患者的主动脉、冠状动脉(和由此延伸的分支)和心脏的医学成像数据803，如ccta数据(例如在图2的步骤100中获得)。输入802还可包括从患者测量的其它生理学数据804，如患者的臂血压、心率和/或其它测量值(例如在图2的步骤100中获得)。其它生理学数据804可以非侵入方式获得。输入802可用以进行下述步骤。

患者心肌组织的三维几何模型可依据成像数据803产生(步骤810)，并且几何模型可分成片段或体积(步骤812)。举例来说，图31展示三维几何模型846，所述三维几何模型846包括分到片段842中的患者心肌组织的三维几何模型838。个别片段842的尺寸和位置可依据以下来确定：冠状动脉(和由此延伸的分支)的流出边界324(图8)的位置、个别片段842(例如邻近血管)中的血管或与各别片段842(例如邻近血管)连接的血管的尺寸，等。将几何心肌模型838分到片段842中可使用各种已知方法进行，所述已知方法如为快速行进法、通用快速行进法(generalizedfastmarchingmethod)、水平集方法、扩散等式、支配经过多孔介质的流的等式，等。

三维几何模型还可包括患者的主动脉和冠状动脉(和由此延伸的分支)的一部分，所述部分可依据成像数据803来建模(步骤814)。举例来说，图31的三维几何模型846包括患者的主动脉和冠状动脉(和由此延伸的分支)的三维几何模型837，以及在步骤810中创建的患者的心肌组织的三维几何模型838。

回头参考图29，例如上文关于图3的步骤402所述，可进行计算分析以确定解，所述解包括关于在由用户确定的身体条件下患者的冠状动脉血流的信息(步骤816)。举例来说，身体条件可包括休息、所选程度的充血、所选程度的运动或用力或其它条件。所述解可提供在步骤814中建模的患者的解剖结构中各种位置并且在指定身体条件下的信息，如血液流量和压力。计算分析可使用来源于集中参数模型或一维模型的在流出边界324(图8)的边界条件进行。可产生一维模型来填充片段842，如下文关于图30所述。

依据在步骤816中确定的血流信息，可计算向在步骤812中创建的心肌的各别片段842中的血流灌注(步骤818)。举例来说，灌注可如下计算：将来自流出边界324(图8)各出口的流量除以出口所灌注的分割心肌的体积。

在步骤818中确定的心肌各别片段的灌注可显示在步骤810或步骤812中产生的心肌几何模型(例如图31中所示的患者心肌组织的三维几何模型838)上(步骤820)。举例来说，图31展示几何模型838的心肌的片段842可用不同色调或颜色说明以指示血流灌注到各别片段842中。

图30展示根据一个示例性实施方案，有关方法820的另一示意图，所述方法820用于提供有关特定患者的心肌灌注的各种信息。方法820可在上述计算机系统中实施，所述计算机系统例如是用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。

方法820可使用一个或多个输入832来进行，输入832可包括患者的主动脉、冠状动脉(和由此延伸的分支)和心脏的医学成像数据833，如ccta数据(例如在图2的步骤100中获得)。输入832可用以进行下述步骤。

患者心肌组织的三维几何模型可依据成像数据833创建(步骤835)。模型还可包括患者的主动脉和冠状动脉(和由此延伸的分支)的一部分，所述部分还可依据成像数据803来创建。举例来说，如上所述，图31展示三维几何模型836，所述三维几何模型836包括患者的主动脉和冠状动脉(和由此延伸的分支)的几何模型837，以及患者心肌组织的几何模型838。步骤835可包括上述图29的步骤810和步骤814。

回头参考图30，几何心肌模型838可分成体积或片段842(步骤840)。步骤840可包括上述图29的步骤812。如上所述，图31展示三维几何模型846，所述三维几何模型846包括分到片段842中的患者心肌组织的几何模型838。

回头参考图30，几何模型846可被修改来包括冠状动脉树中的下一代分支857(步骤855)。分支857(图31中虚线所示)的位置和尺寸可依据冠状动脉(和由此延伸的分支)的中线来确定。中线可例如依据成像数据833来确定(步骤845)。一种算法还可用以依据形态计量学模型(用以预测在流出边界324(图8)的已知出口下游的血管位置和尺寸的模型)和/或与血管尺寸有关的生理学分支规律来确定分支857的位置和尺寸(步骤850)。形态计量学模型可扩充到几何模型837中包括的冠状动脉(和由此延伸的分支)的下游端，并且提供于心外膜表面(心脏组织的外层)上或含在心肌壁的几何模型838内。

心肌可依据在步骤855中创建的分支857进一步分割(步骤860)。举例来说，图31展示可分成子体积或子段862的片段842。

可在子段862中创建其它分支857，并且子段862可进一步分割成较小片段867(步骤865)。创建分支以及将体积分成子段的步骤可重复到获得体积尺寸和/或分支尺寸的所需分解为止。已在步骤855和步骤865中扩充为包括新分支857的模型846接着可用以计算向子段(如在步骤865中产生的子段867)中的冠状动脉血流和心肌灌注。

因此，扩充模型可用以进行上述计算分析。计算分析结果可将有关血流的信息从患者特异性冠状动脉模型(例如图31的模型837)提供到所产生的形态计量学模型(包括在步骤855和步骤865产生的分支)中，所述信息可延伸到在步骤865中产生的灌注子段867的每一个中。计算分析可使用静态心肌灌注体积或动态模型进行，所述模型并有来自偶合心脏机械学模型的数据。

图32展示根据一个示例性实施方案，有关方法870的另一示意图，所述方法870用于提供有关特定患者的心肌灌注的各种信息。方法870可在上述计算机系统中实施，所述计算机系统例如是用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。

方法870可使用一个或多个输入872来进行。输入872可包括患者的主动脉、冠状动脉(和由此延伸的分支)和心脏的医学成像数据873，如ccta数据(例如在图2的步骤100中获得)。输入872还可包括从患者测量的其它生理学数据874，如患者的臂血压、心率和/或其它测量值(例如在图2的步骤100中获得)。其它生理学数据874可以非侵入方式获得。输入872还可包括从患者测量(例如使用ct、pet、spect，等)的心脏灌注数据875。输入872可用以进行下述步骤。

患者的主动脉和冠状动脉(和由此延伸的分支)的三维几何模型可依据成像数据873来创建(步骤880)。举例来说，图31展示患者的主动脉和冠状动脉(和由此延伸的分支)的三维几何模型837。步骤880可类似于上述图29的步骤814。

计算分析例如可如上文关于图3的步骤402所述来进行以确定解，所述解包括关于在由用户确定的身体条件下患者的冠状动脉血流的信息(步骤882)。举例来说，身体条件可包括休息、所选程度的充血、所选程度的运动或用力或其它条件。所述解可提供在步骤880中建模的患者的解剖结构中各种位置并且在指定身体条件下的信息，如血液流量和压力。步骤882可类似于上述图29的步骤816。

而且，患者心肌组织的三维几何模型可依据成像数据873创建(步骤884)。举例来说，如上所述，图31展示三维几何模型836，其包括患者心肌组织的三维几何模型838(例如在步骤884中创建)和患者的主动脉和冠状动脉(和由此延伸的分支)的三维几何模型837(例如在步骤880中创建)。步骤884可类似于上述图29的步骤810。

几何模型可分成片段或子体积(步骤886)。举例来说，图31展示几何模型846，所述几何模型846包括患者的分成片段842的心肌组织的模型838。步骤886可类似于上述图29的步骤812。

依据在步骤882中确定的血流信息，可计算向在步骤886中创建的心肌的各别片段842中的血流灌注(步骤888)。步骤888可类似于上述图29的步骤818。

心肌各别片段的计算灌注可显示在步骤884或步骤886中产生的心肌几何模型(例如图31中所示的患者心肌组织的三维几何模型838)上(步骤890)。举例来说，图31展示几何模型838的心肌的片段842可用不同色调或颜色说明以指示血流灌注到各别片段842中。步骤890可类似于上述图29的步骤820。

在步骤890中绘制于心肌三维几何模型上的模拟灌注数据可与测量的心脏灌注数据875相比(步骤892)。可例如在基于体元的心肌表示法或不同的心肌离散表示法(例如有限元网)上进行比较。比较可使用心肌三维表示上各种颜色和/或色调来指示模拟灌注数据和测量灌注数据的差异。

在步骤880中创建的三维几何模型的出口处的边界条件可调整，以减小模拟灌注数据和测量灌注数据之间的误差(步骤894)。举例来说，为减小误差，可调整边界条件，以使向一个区域(例如片段842、862或867)(其中模拟灌注低于测量灌注)供料的血管的指定流动阻力可减小。可调整边界条件的其它参数。或者，可修改模型的分支结构。举例来说，如上文关于图30和图31所述，可扩充在步骤880中创建的几何模型来创建形态计量学模型。边界条件和/或形态计量学模型的参数可使用参数估计或数据同化方法来凭经验或系统地调整，所述方法如为标题为“methodfortuningpatient-specificcardiovascularsimulations”的美国专利申请公布号2010/0017171中所述的方法，或其它方法。

可重复图32的步骤882、步骤888、步骤890、步骤892、步骤894和/或其它步骤，例如直到模拟灌注数据和测量灌注数据之间的误差低于预定临限值为止。因此，计算分析可使用有关解剖学信息、冠状动脉血流信息和心肌灌注信息的模型来进行。所述模型可适用于诊断目的并且适用于预测医药疗法、介入疗法或手术疗法的益处。

因此，在休息和/或应激条件下的冠状动脉血流和心肌灌注可在从三维医学成像数据构建的患者特异性几何模型中模拟。测量的心肌灌注数据可与模拟心肌灌注结果组合使用，以调整边界条件直到模拟心肌灌注结果在给定容许量内匹配测量的心肌灌注数据为止(例如上文关于图32所述)。可提供较精确的患者特异性冠状动脉血流计算，并且心脏病专家可能够预测在测量数据可能不可得的情况下(如当模拟患者在最大运动或用力、模拟治疗或其它条件下时)的冠状动脉血流和心肌灌注。

左心室和/或右心室心肌的患者特异性三维模型可分成灌注片段或子体积。而且，从医学成像数据确定的冠状动脉的患者特异性三维几何模型可与其余冠状动脉树的一部分的形态计量学模型组合，所述其余冠状动脉树是在心外膜表面上或含于由灌注子体积表示的左心室和/或右心室心肌壁中以形成扩充模型。可计算扩充模型中给定(例如患病)位置下游的总心肌体积的百分比。还可计算在扩充模型中给定(例如患病)位置的总心肌血流的百分比。扩充模型可用以计算冠状动脉血流和心肌灌注。还可修改冠状动脉血流模型直到模拟灌注在指定容许量内匹配测量的灌注数据为止。

b.评估斑块脆弱性

计算分析还可提供定量作用于斑块的患者特异性生物力学力的结果，所述斑块可在患者的主动脉和冠状动脉(和由此延伸的分支)中积累，例如冠状动脉粥样硬化斑块。生物力学力可由脉动压力、流量和心脏运动引起。

图33展示沿血管壁902积累的斑块900的一个实施例，所述血管壁902如为主要冠状动脉之一或由此延伸的分支之一的壁。斑块上游端和下游端的压力和/或表面积的差异可产生力904，所述力904至少沿血流方向作用于斑块900，例如由血液流经血管引起。另一力906可至少沿朝着并且垂直于血管壁902的方向作用于斑块900的表面。力906可由流经血管的血液的血压所引起。又一力908可至少沿血流方向作用于斑块900的表面，并且可归因于在休息、运动等期间的血液动力学力。

结果还可评估斑块破裂(例如当在血管壁上累积的斑块变得不稳定并且脱落或破开时)的风险以及可能受所述破裂影响的心肌体积。可在各种模拟生理条件(如休息、运动等)下评估结果。斑块破裂风险可定义为模拟斑块应力与斑块强度的比率，所述斑块强度是使用来源于ccta或mri的物质组成数据(例如在图2的步骤100中测定)来估计。

举例来说，图34展示计算分析可输出的结果的一个实施例。结果可包括图31的三维几何模型846，所述三维几何模型846可包括患者的主动脉和冠状动脉(和由此延伸的分支)的三维几何模型837，以及患者的分成片段842的心肌组织的三维几何模型838。结果还可指示冠状动脉(由此延伸的分支)之一中的位置910，其中斑块可被确定为脆弱的，并且位置910可依据斑块破裂风险评估(将在下文更详细地描述)和/或依据用户输入而被鉴别出。而且，如图34中所示，由于在位置910鉴别出的斑块破裂，因此心肌片段912(片段842的多个的)可鉴别为具有高的低灌注概率。

图35和图36是根据一个示例性实施方案展示方法920的方面的示意图，所述方法920用于提供有关评估特定患者的斑块脆弱性、心肌体积风险和心肌灌注风险的各种信息。方法920可在上述计算机系统中实施，所述计算机系统例如是用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。方法920可使用一个或多个输入922来进行，并且可包括依据输入922来产生一个或多个模型930、依据一个或多个模型930进行一种或多种生物力学分析940，以及依据模型930和生物力学分析940提供各种结果。

输入922可包括患者的主动脉、冠状动脉(和由此延伸的分支)和心脏的医学成像数据923，如ccta数据(例如在图2的步骤100中获得)。输入922还可包括从患者测量的其它生理学数据924，如患者的臂血压、心率和/或其它测量值(例如在图2的步骤100中获得)。其它生理学数据924可以非侵入方式获得。输入922可用以产生模型930和/或进行下述生物力学分析940。

如上所述，一个或多个模型930可依据输入922来产生。举例来说，方法920可包括产生血液动力学模型932，所述血液动力学模型932包括在整个患者解剖结构的三维几何模型中各种位置的计算血液流量和压力信息。患者解剖结构的模型可使用以下创建：医学成像数据923，例如在图3的步骤306中产生的图8的实体模型320，和/或在图3的步骤312中产生的图17-19的网380，并且在一个示例性实施方案中，血液动力学模型932可为模拟血压模型50(图1)、模拟血流模型52(图1)、cffr模型54(图1)，或在进行计算分析之后产生的其它模拟，例如上文关于图3的步骤402所述。固体力学模型(包括流体结构互动模型)可以用已知数值方法的计算分析求解。斑块和血管的性质可建模成线性或非线性、各向同性或各向异性。解可提供斑块以及斑块和血管之间的界面的应力和应变。在图36中所示的示例性实施方案中，血液动力学模型932是cffr模型54。

方法920可包括使用血液动力学模型932如下来进行生物力学分析940：计算在各种生理状态(如休息、不同程度的运动或用力等)下因血液动力学力而作用于斑块腔面的压力906(图33)和剪应力908(图33)(步骤942)。压力906和剪应力908可依据来自血液动力学模型932的信息(例如血液压力和流量)来计算。

任选地，方法920还可包括产生用于从四维成像数据定量血管变形的几何分析模型934，所述四维成像数据例如是在心动周期的多个阶段(如收缩阶段和舒张阶段)获得的成像数据。成像数据可使用各种已知成像方法获得。几何分析模型934可包括关于在心动周期的不同阶段，例如由心脏运动造成的血管位置、变形、取向和尺寸的信息。举例来说，可通过几何分析模型934模拟患者的主动脉、冠状动脉(和由此延伸的分支)和斑块的各种类型的变形，如纵向变长(伸长)或缩短、扭曲(扭转)、径向膨胀或压缩和弯曲。

方法920可包括使用几何分析模型934如下进行生物力学分析940：计算患者的主动脉、冠状动脉(和由此延伸的分支)和斑块由心脏诱导性脉动压力造成的各种变形特征，如纵向变长(伸长)或缩短、扭曲(扭转)、径向膨胀或压缩，和弯曲等(步骤944)。这些变形特征可依据来自几何分析模型934的信息来计算，所述信息例如是在心动周期的多个阶段内，血管位置、取向和尺寸的改变。

计算变形特征可通过确定所建模几何形状(例如患者的主动脉、冠状动脉(和由此延伸的分支)、斑块等的几何形状)的中线或表面网来简化。为确定所建模几何形状在不同阶段之间的变化，分支口、钙化病变和软斑块可用作陆标(landmark)。在不具有陆标的区域中，沿所建模几何形状的长度的截面积概况可用以鉴别在两个图像帧之间的相应位置(以“记录(register)”两个图像帧)。依据原始图像数据的可变形记录算法可用以提取三维变形场。计算的三维变形场接着可投射到与所建模几何形状(例如血管长度)对准的曲线轴以计算变形场的切线分量和法向分量。所建模几何形状(例如血管长度)、分支分开的角度和收缩和心脏舒张之间曲率的所得差异可用以确定血管经历的应变。

方法920还可包括产生用于从医学成像数据923确定斑块组成和性质的斑块模型936。举例来说，斑块模型936可包括关于斑块的密度和其它物质性质的信息。

方法920还可包括产生血管壁模型938，所述血管壁模型938是用于计算关于斑块、血管壁和/或斑块和血管壁之间的界面的信息。举例来说，血管壁模型938可包括关于应力和应变的信息，所述信息可依据斑块模型936中包括的斑块组成和性质、在步骤942中计算的压力906和剪应力908和/或在步骤944中计算的变形特征来计算。

方法920可包括使用血管壁模型938如下进行生物力学分析940：计算血液动力学力和心脏运动诱导性应变造成的斑块上应力(例如急性应力或积累应力)(步骤946)。举例来说，可计算作用于斑块的流量诱导力904(图33)。由血液动力学力和心脏运动诱导性应变造成的斑块上应力或力可依据来自血管壁模型938的信息(例如斑块上应力和应变)来计算。

方法920可包括依据上述模型930中一个或多个和生物力学分析940中一个或多个来确定其它信息。

可计算斑块破裂脆弱性指数(步骤950)。斑块破裂脆弱性指数可例如依据总血液动力学应力、应力频率、应力方向和/或斑块强度或其它性质来计算。举例来说，包围目标斑块的区域可与斑块的三维模型930(如斑块模型936)分离。斑块强度可从斑块模型936中提供的物质性质确定。由脉动压力、流量和心脏运动造成的目标斑块上的血液动力学应力和组织应力可在模拟基线和运动(或用力)条件下通过使用以下来计算：在步骤946中预先计算的血液动力学应力和运动诱导性应变。斑块脆弱性可依据斑块应力与斑块强度的比率来评估。

还可计算心肌容量风险指数(mvri)(步骤952)。mvri可定义为在动脉树中给定位置，受斑块破裂和血管梗塞(闭塞或阻塞)影响的总心肌容量的百分比。mvri可依据给定斑块下游血管提供的心肌的部分来计算，这可相对于下游血管的尺寸来考虑斑块尺寸，以及依据三维血液动力学解，斑块可流入不同血管中的概率。

心肌可被建模并分成由血液动力学模拟中的各血管提供的片段842(例如关联图30的步骤835和步骤840所述)。几何模型可被修改为包括冠状动脉树中分支857的下一代(例如关联图30的步骤855所述)，并且心肌可被进一步分割(例如关联图30的步骤860所述)。可在子段862中产生其它分支857，并且子段862可进一步分割成较小片段867(例如关联图30的步骤865所述)。如先前描述的生理学关系可用以关联血管尺寸与成比例量的所提供的心肌。

可确定破裂斑块所行进的可能路径。血液动力学解可用以确定斑块片段或栓塞物可流入不同下游血管中的机会百分比。

可将破裂斑块的尺寸与下游血管的尺寸相比以确定斑块可在哪里最终产生流障碍。这个信息可与脆弱性指数组合以提供可能受破裂斑块影响的心肌容量的概率图。mvri在各可能的受影响的片段都可被赋值。图34展示片段912的一个实施例，在所述片段912处，在远端血管中位置910的脆弱斑块影响小面积心肌的概率较高。

还可计算心肌灌注风险指数(mpri)(步骤954)。mpri可定义为在动脉树中给定位置，受斑块破裂和血管梗塞影响的总心肌血流的百分比。举例来说，相比lad动脉的邻近部分中的斑块破裂，lad动脉的远端部分中的斑块破裂将得到较低mvri和较低mpri。然而，如果受进料血管中的脆弱斑块影响的心肌容量的一部分不可行(例如由于可在心肌梗塞之后形成的疤痕组织)，那么这些指数可不同。因此，mpri指示向心肌片段的可能灌注损失，而非如mvri指示的受影响容量。可计算向图31的各片段842、862或867的灌注率，并且灌注损失可依据脆弱性指数、血液动力学解，和斑块和血管的尺寸来计算。

因此，可计算由脉动血压、脉动血流、脉动血液剪应力和/或脉动心脏运动造成的斑块应力，并且斑块强度可依据医学成像数据来估计，并且可定量有关斑块脆弱性、心肌容量风险和心肌灌注风险的指数。

viii.其它应用

上述实施方案与评估关于患者的冠状动脉血流的信息相关。或者，实施方案还可针对身体其它区域的血流来调适，所述身体其它区域如但不限于颈动脉、周边动脉、腹部动脉、肾动脉、大腿动脉、腿弯部动脉和脑动脉。

a.为颅内和颅外血流建模

现将描述有关脑动脉的实施方案。众多疾病可影响颅外或颅内动脉中的血液流量和压力或受颅外或颅内动脉中的血液流量和压力的影响。颅外(例如颈动脉和脊椎)动脉中的动脉粥样硬化疾病可限制血液流向脑部。动脉粥样硬化疾病的严重表现可导致短暂性缺血性发作或缺血性中风。颅内或颅外动脉中的动脉瘤疾病可造成以下风险：导致缺血性中风的栓塞或导致出血性中风的动脉瘤破裂。其它病状(如头外伤、高血压、头颈癌、动静脉畸形、直立耐受不能等)也可影响脑血流量。此外，脑血流量减小可诱发如昏厥的症状，或影响慢性神经病症，如继阿兹海默氏病(alzheimer'sdisease)或帕金森氏病(parkinson'sdisease)之后的痴呆。

已知患有或怀疑患有颅外或颅内动脉疾病的患者通常可接受一种或多种以下非侵入性诊断性测试：us、mri、ct、pet。然而，这些测试不能有效提供用于大多数患者的颅外和颅内动脉的解剖学数据和生理学数据。

图37是脑动脉图，所述脑动脉包括颅内(在颅骨内)和颅外(在颅骨外)动脉。用于确定关于患者特异性颅内和颅外血流的信息的方法一般可类似于如上所述用于确定关于患者特异性冠状动脉血流的信息的方法。

图38是展示方法1000的方面的示意图，所述方法1000用于提供有关特定患者的颅内血流量和颅外血流的各种信息。方法1000可在计算机系统中实施，所述计算机系统例如类似于用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。方法1000可使用一个或多个输入1010来进行，并且可包括依据输入1010来产生一个或多个模型1020、依据输入1010和/或模型1020对一个或多个条件1030赋值，以及依据模型1020和条件1030得出一个或多个解1040。

输入1010可包括患者的颅内和颅外动脉(例如患者的主动脉、颈动脉(图37中所示)、脊椎动脉(图37中所示)和脑部)的医学成像数据1011，如ccta数据(例如以与上文关于图2的步骤100所述类似的方式获得)。输入1010还可包括测量1012患者的臂血压、颈动脉血压(例如使用张力测定法)和/或其它测量(例如以与上文关于图2的步骤100所述类似的方式获得)。测量1012可以非侵入方式获得。输入1010可用以产生模型1020和/或确定下述条件1030。

如上所述，一个或多个模型1020可依据输入1010来产生。举例来说，方法1000可包括依据成像数据1011来产生患者的颅内和颅外动脉的一个或多个患者特异性三维几何模型(步骤1021)。三维几何模型1021可使用如上关于产生图8的实体模型320和图17-19的网380所述的类似方法来产生。举例来说，与图3的步骤306和步骤312类似的步骤可用以产生呈现患者的颅内和颅外动脉的三维实体模型和网。

回头参考图38，方法1000还可包括产生一个或多个基于物理学的血流模型(步骤1022)。举例来说，血流模型可为呈现以下的模型：经过在步骤1021中产生的患者特异性几何模型的流、心脏和主动脉循环、远端颅内和颅外循环，等。在三维几何模型1021的流入边界和/或流出边界，血流模型可包括如上文关于图3的步骤310所述的降阶模型，例如集中参数模型或分布(一维波传播)模型等。或者，流入边界和/或流出边界可被赋予速度、流动速率、压力或其它特征等的各别指定值或值域。作为另一替代方案，流入边界可与心脏模型(例如包括主动脉弓)偶合。流入边界和/或流出边界的参数可调整来匹配所测量或选择的生理条件，生理条件包括但限于心输出量和血压。

如上所述，一个或多个条件1030可依据输入1010和/或模型1020来确定。条件1030包括针对在步骤1022(和图3的步骤310)中确定的边界条件所计算的参数。举例来说，方法1000可包括如下确定条件：依据成像数据1011(例如以与上文关于图3的步骤240所述类似的方式获得)来计算患者特异性脑容量或头容量(步骤1031)。

方法1000可包括如下确定条件：使用在步骤1031中计算的脑容量或头容量、依据关系q＝qom^α的休息脑血流量q(其中α是预设标度指数，m是从脑容量或头容量确定的脑质量，并且qo是预设常数)进行计算(例如类似于上文关于在图3的步骤310中确定集中参数模型描述的生理关系)(步骤1032)。或者，所述关系可具有形式q∝qom^α，如上文在图3的步骤310中关于确定集中参数模型所述。

方法1000还可包括如下确定条件：使用在步骤1032中计算的所得冠状动脉流量和患者的测量血压1012、总休息大脑阻力进行计算(例如类似于上文关于在图3的步骤310中确定集中参数模型描述的方法)(步骤1033)。举例来说，在基线(休息)条件下在步骤1032中确定的在三维几何模型1021的流出边界的总脑血流量q以及测量血压1012可用以依据预设的用实验方式得出的等式来确定在流出边界的总阻力r。阻力、电容、感应系数和与集中参数模型中所用的各种电组件相关的其它变量可并入边界条件中(例如上文在图3的步骤310中关于确定集中参数模型所述)。

方法1000还可包括如下确定条件：使用在步骤1033和模型1020中计算的总休息大脑阻力、个别颅内和颅外动脉的个别阻力进行计算(步骤1034)。举例来说，类似于上文关于图3的步骤310描述的方法，在步骤1033中计算的总休息大脑阻力r可依据个别颅内和颅外动脉的远端的尺寸(例如从在步骤1021中产生的几何模型确定)，并且依据关系r＝rod^β(其中r是在具体远端的流动阻力，并且ro是预设常数，d是尺寸(例如远端的直径)，并且β是预设幂律指数，如上文关于在图3的步骤310中确定集中参数模型所述)分布到个别颅内和颅外动脉。

回头参考图38，方法1000可包括依据患者的一个或多个身体条件来调整边界条件(步骤1035)。举例来说，在步骤1031-1034中确定的参数可依据以下来修改：解1040是否意图模拟休息、不同程度的应激、不同程度的压力感受性反应或其它自主反馈控制、不同程度的充血、不同程度的运动、用力、高血压或低血压、不同药物治疗、体位变化和/或其它条件。参数(例如有关在流出边界的边界条件的参数)还可依据以下来调整：颅内和颅外动脉的舒张血管能力(血管变宽的能力)，例如归因于微血管功能障碍或内皮健康。

依据输入1010、模型1020和条件1030，可进行计算分析，例如上文关于图3的步骤402所述，以便确定解1040，所述解1040包括关于患者的冠状动脉血流在步骤1035中选择的身体条件下的信息(步骤1041)。可从解1040提供的信息的实例可类似于上文关于图1和图21-24提供的实例，例如模拟血压模型、模拟血流模型，等。结果还可用以确定例如流动速率、总脑流量、血管壁剪应力、牵引力或剪切力，所述力作用于血管壁或动脉粥样硬化斑块或动脉瘤、颗粒/血液滞留时间、血管壁活动、血液剪切速率，等。这些结果还可用以分析离开血管系统中特定区域的栓塞物在哪里可很可能因血液循环而移动。

计算机系统可允许用户模拟几何形状的各种变化。举例来说，模型1020，例如在步骤1021中产生的患者特异性几何模型可被修改来预测动脉梗塞(例如急性梗塞)的影响。在一些手术程序中，如当去除癌性肿瘤时，一根或多根颅外动脉可被损伤或去除。因此，在步骤1021中产生患者特异性几何模型还可修改来模拟防止血液流到一根或一根以上颅外动脉中的影响，以预测并行通路用于向患者提供足够血流的可能。

计算机系统可允许用户模拟各种治疗选择的结果，所述治疗选择如为介入修复或手术修复例如急性梗塞。通过用如上文关于图27和图28所述的降阶模型替代如上所述的呈现颅内和颅外动脉的三维实体模型或网，模拟可较快速地进行。因此，降阶模型，如一维模型或集中参数模型可较有效并快速地求出患者特异性模型中的血液流量和压力并且显示解的结果。

特定患者对舒张血管刺激物的反应可依据患者休息时的血液动力学信息或依据不同疾病病况的基于群体的数据来预测。举例来说，在基线(休息)中，用依据幂律和脑质量(例如上文关于步骤1032所述)赋值的流动分布运行模拟(例如上文在步骤1041中所述)。可调整阻力值(例如在步骤1033和步骤1034中确定)以允许足够的灌注。或者，将来自有如糖尿病、药物治疗和以往心脏事件的因素的患者群体的数据用以对不同阻力赋值。单独或与血液动力学信息(例如壁剪应力，或流量和血管尺寸的关系)组合的在休息条件下的阻力的调整可用以确定远端大脑血管扩张的剩余能力。需要阻力减小来满足休息流量要求的患者，或流量与血管尺寸比较高的患者可具有减小的在生理应激下进一步扩张其血管的能力。

在脑动脉的个别片段上的流动速率和压力梯度(例如在步骤1041中确定)可用以计算大脑动脉阻力。大脑动脉阻力可计算为患者特异性几何模型中包括的颅外和颅内动脉的部分的等效阻力，所述患者特异性几何模型是从医学成像数据产生(例如在步骤1021中产生)。大脑动脉阻力的临床意义可在于解释为什么在颅外和/或颅内动脉中患有扩散性动脉粥样硬化的患者可显示昏厥(暂时性意识丧失或姿势，例如昏晕)或局部缺血(血液供给限制)的症状。

而且，例如依据在步骤1041中确定的流量信息以及在步骤1031中计算的脑组织容量或质量，可计算在基线或改变的生理学条件下每单位脑组织容量(或质量)的流量。这个计算可适用于理解血流减少对慢性神经病症的影响。这个计算还可适用于选择或改进医学疗法，例如抗高血压剂的剂量。其它结果可包括定量外伤、脑震荡、外部生理应激、过量重力、失重、太空飞行、深海减压(例如弯曲)等的影响。

组合的患者特异性解剖学(几何)模型和生理学(基于物理学的)模型可用以确定不同药物治疗或生活方式变化(例如停止吸烟、改变膳食或增加身体活动)对脑动脉血流的作用，所述药物治疗或生活方式变化会改变心率、心搏量、血压或大脑微循环功能。组合模型还可用以确定替代形式和/或不同程度的暴露于可能的外来力(例如当踢足球时、在太空飞行期间、当斯库巴潜水时、在坐飞机飞行期间，等)的身体活动或风险对脑动脉血流的作用。所述信息可用以鉴别对特定患者来说可能安全并有效的身体活动的类型和程度。组合模型也可用以预测经皮介入对脑动脉血流的可能益处以选择最优的介入策略，和/或预测颈动脉内膜切除术或外部颈动脉到内部颈动脉旁路移植对脑动脉血流的可能益处以选择最优的手术策略。

组合模型还可用以说明动脉疾病负担增大对脑动脉血流的可能的有害作用，并且使用机械或现象学疾病进展模型或经验数据来预测进展中疾病何时可对流向脑部的血流引起损害。所述信息可实现确定“保证期”，其中使用非侵入性成像观察到最初没有热力学有意义的疾病的患者可能不被预期需要医药疗法、介入疗法或手术疗法，或者，如果继续有不利因素，那么可确定出现进展的比率。

组合模型也可用以说明由疾病负担减小引起的对脑动脉血流的可能的有益作用，并且使用机械或现象学疾病进展模型或经验数据来预测疾病消退何时可引起到脑部的血流增加。所述信息可用以指导医学管理程序，所述医学管理程序包括但不限于改变膳食、增加身体活动、开出士他汀或其它医药的处方，等。

组合模型还可用以预测动脉梗塞的影响。在一些手术程序(如去除癌性肿瘤)中，一些颅外动脉可被损伤或去除。模拟防止血液流到颅外动脉之一的作用可允许预测并行通路向特定患者提供足够血流的可能。

i.评估大脑灌注

可计算其它结果。举例来说，计算分析可提供定量大脑灌注(经过大脑的血流)的结果。定量大脑灌注可有助于鉴别脑血流量减少的区域。

图39展示根据一个示例性实施方案，有关方法1050的示意图，所述方法1050用于提供有关特定患者的大脑灌注的各种信息。方法1050可在上述计算机系统中实施，所述计算机系统例如类似于用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。

方法1050可使用一个或多个输入1052来进行。输入1052可包括患者的颅内和颅外动脉(例如患者的主动脉、颈动脉(图37中所示)、脊椎动脉(图37中所示)和脑部)的医学成像数据1053，如ccta数据(例如以与上文关于图2的步骤100所述类似的方式获得)。输入1052还可包括从患者测量的其它生理学数据1054，如患者的臂血压、心率和/或其它测量值(例如以与上文关于图2的步骤100所述类似的方式获得)。其它生理学数据1054可以非侵入方式获得。输入1052可用以进行下述步骤。

患者脑组织的三维几何模型可依据成像数据1053创建(步骤1060)，并且几何模型可分成片段或体积(步骤1062)(例如用与如上文关于图29-32所述类似的方式)。个别片段的尺寸和位置可依据以下来确定：颅内和颅外动脉的流出边界的位置、在各别片段(例如邻近血管)中或与各别片段(例如邻近血管)连接的血管的尺寸，等。将几何模型分成片段可使用各种已知方法进行，所述已知方法如为快速行进法、通用快速行进法、水平集方法、扩散等式、经过多孔介质的等式支配流，等。

三维几何模型还可包括患者的颅内和颅外动脉的一部分，所述部分可依据成像数据1053来建模(步骤1064)。举例来说，在步骤1062和步骤1064中，可创建包括脑组织和颅内和颅外动脉的三维几何模型。

计算分析例如可如上文关于图3的步骤402所述来进行以确定解，所述解包括关于在由用户确定的身体条件下患者的脑血流量的信息(步骤1066)。举例来说，身体条件可包括休息、不同程度的应激、不同程度的压力感受性反应或其它自主反馈控制、不同程度的充血、不同程度的运动或用力、不同药物治疗、体位变化和/或其它条件。所述解可提供在步骤1064中建模的患者的解剖结构中各种位置并且在指定身体条件下的信息，如血液流量和压力。计算分析可使用来源于集中参数模型或一维模型的在流出边界的边界条件来进行。可产生一维模型来填充脑组织的片段，如下文关于图40所述。

依据在步骤1066中确定的血流信息，可计算向在步骤1062中创建的脑部的各别片段中的血流灌注(步骤1068)。举例来说，灌注可如下计算：将来自流出边界各出口的流量除以出口所灌注的被分割脑部的容量。

在步骤1068中确定的脑部各别片段的灌注可显示在步骤1060或步骤1062中产生的脑部几何模型上(步骤1070)。举例来说，在步骤1060中创建的几何模型中所示的脑部片段可用不同色调或颜色说明以指示血流灌注到各别片段中。

图40展示根据一个示例性实施方案，有关方法1100的另一示意图，所述方法1100用于提供有关特定患者的大脑灌注的各种信息。方法1100可在上述计算机系统中实施，所述计算机系统例如类似于用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。

方法1100可使用一个或多个输入1102来进行，输入1102可包括患者的主动脉、颈动脉(图37中所示)、脊椎动脉(图37中所示)和脑部的医学成像数据1103，如ccta数据(例如以与上文关于图2的步骤100所述类似的方式获得)。输入1102可用以进行下述步骤。

患者脑组织的三维几何模型可依据成像数据1103创建(步骤1110)。模型还可包括患者的主动脉、颈动脉(图37中所示)和脊椎动脉(图37中所示)的一部分，所述部分还可依据成像数据1103创建。举例来说，如上所述，可创建包括脑组织和颅内和颅外动脉的三维几何模型。步骤1110可包括上述图39的步骤1060和步骤1064。

在步骤1110中创建的几何脑组织模型可分成体积或片段(步骤1112)。步骤1112可包括上述图39的步骤1062。几何脑组织模型还可进一步修改来包括大脑树中分支的下一代(步骤1118)(例如用与如上文关于图29-32所述类似的方式)。分支的位置和尺寸可依据颅内和颅外动脉的中线来确定。中线可例如依据成像数1103来确定(步骤1114)。一种算法还可用以依据形态计量学模型(用以预测在流出边界的已知出口下游的血管位置和尺寸的模型)和/或与血管尺寸有关的生理学分支规律来确定分支的位置和尺寸(步骤1116)。形态计量学模型可扩充到几何模型中包括的颅内和颅外动脉的下游端，并且提供于脑组织外层上或含在脑组织几何模型内。

脑部可依据在步骤1118中创建的分支来进一步片段(步骤1120)(例如用与如上文关于图29-32所述类似的方式)。可在子段中创建其它分支，并且子段可进一步分割成较小片段(步骤1122)(例如用与如上文关于图29-32所述类似的方式)。创建分支以及将体积分成子段的步骤可重复到获得体积尺寸和/或分支尺寸的所需分解为止。已在步骤1118和步骤1122中扩充为包括新分支的几何模型接着可用以计算向子段(如在步骤1122中产生的子段)中的脑血流量和大脑灌注。

因此，扩充模型可用以进行上述计算分析。计算分析结果可将有关血流的信息从患者特异性脑动脉模型提供到所产生的形态计量学模型(包括在步骤1118和步骤1122产生的分支)中，所述信息可延伸到在步骤1122中产生的灌注子段的每一个中。

图41展示根据一个示例性实施方案，有关方法1150的另一示意图，所述方法1150用于提供有关特定患者的大脑灌注的各种信息。方法1150可在上述计算机系统中实施，所述计算机系统例如是用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。

方法1150可使用一个或多个输入1152来进行。输入1152可包括患者的主动脉、颈动脉(图37中所示)、脊椎动脉(图37中所示)和脑部的医学成像数据1153，如ccta数据(例如以与上文关于图2的步骤100所述类似的方式获得)。输入1152还可包括从患者测量的其它生理学数据1154，如患者的臂血压、心率和/或其它测量值(例如在图2的步骤100中获得)。其它生理学数据1154可以非侵入方式获得。输入1152还可包括从患者测量(例如使用ct、pet、spect、mri，等)的脑灌注数据1155。输入1152可用以进行下述步骤。

患者的颅内和颅外动脉的三维几何模型可依据成像数据1153来创建(步骤1160)。步骤1160可类似于上述图39的步骤1064。

计算分析例如可如上文关于图3的步骤402所述来进行以确定解，所述解包括关于在由用户确定的身体条件下患者的脑血流量的信息(步骤1162)。举例来说，身体条件可包括休息、不同程度的应激、不同程度的压力感受性反应或其它自主反馈控制、不同程度的充血、不同程度的运动或用力、不同药物治疗、体位变化和/或其它条件。所述解可提供在步骤1160中建模的患者的解剖结构中各种位置并且在指定身体条件下的信息，如血液流量和压力。步骤1162可类似于上述图39的步骤1066。

而且，患者脑组织的三维几何模型可依据成像数据1153创建(步骤1164)。举例来说，在步骤1160和步骤1164中，可创建包括脑组织和颅内和颅外动脉的三维几何模型。步骤1164可类似于上述图39的步骤1060。

几何模型可分成片段或子体积(步骤1166)。步骤1166可类似于上述图39的步骤1062。

依据在步骤1162中确定的血流信息，可计算向在步骤1166中创建的脑组织的各别片段中的血流灌注(步骤1168)。步骤1168可类似于上述图39的步骤1068。

脑组织各别片段的计算灌注可显示在步骤1164或步骤1166中产生的脑组织几何模型上(步骤1170)。步骤1170可类似于上述图39的步骤1070。

在步骤1170中绘制于脑组织三维几何模型上的模拟灌注数据可与测量的大脑灌注数据1155相比(步骤1172)。比较可使用脑组织三维表示上各种颜色和/或色调来指示模拟灌注数据和测量灌注数据的差异。

在步骤1160中创建的三维几何模型的出口处的边界条件可调整，以减小模拟灌注数据和测量灌注数据之间的误差(步骤1174)。举例来说，为减小误差，可调整边界条件，以使向一个区域(例如在步骤1166中创建的片段)(其中模拟灌注低于测量灌注)供料的血管的指定流动阻力可减小。可调整边界条件的其它参数。或者，可修改模型的分支结构。举例来说，如上文关于图40所述，可扩充在步骤1160中创建的几何模型来创建形态计量学模型。边界条件和/或形态计量学模型的参数可使用参数估计或数据同化方法来凭经验或系统地调整，所述方法如为标题为“methodfortuningpatient-specificcardiovascularsimulations”的美国专利申请公布号2010/0017171中所述的方法，或其它方法。

图41的步骤1162、1168、1170、1172、1174和/或其它步骤可重复，例如直到模拟灌注数据和测量灌注数据之间的误差低于预定临限值为止。因此，计算分析可使用有关解剖学信息、脑血流量信息和大脑灌注信息的模型来进行。所述模型可适用于诊断目的并且适用于预测医药疗法、介入疗法或手术疗法的益处。

因此，可计算在基线条件或改变的生理学状态下的颅外和颅内动脉血流和大脑灌注。大脑灌注数据可与模拟大脑灌注结果组合使用，以调整颅内动脉血流计算值的边界条件，直到模拟大脑灌注结果在给定容许量内匹配测量的大脑灌注数据为止。因此，可提供较精确的患者特异性颅外和颅内动脉血流计算，并且医师可预测在测量数据可能不可得的情况下(例如某些身体条件，如运动、用力、体位变化或模拟治疗)的脑动脉血流和大脑灌注。脑部的患者特异性三维模型可分成灌注片段或子体积，并且可确定患者是否接收到适当最小的灌注到脑部的各种区域。

颅内动脉的患者特异性三维几何模型可从医学成像数据产生，并且与由灌注片段或子体积表示的其余颅内动脉树的一部分的形态计量学模型(例如上文关于图40所述)组合以形成扩充模型。可计算扩充模型中给定(例如患病)位置下游的总脑容量(或质量)的百分比。而且，可计算在扩充模型中给定(例如患病)位置的总脑血流量的百分比。另外，功能成像研究(例如功能性磁共振成像(fmri))、灌注ct或mri中记录的不足接着可追踪到进料血管中的疾病、解剖学变体、受损自调节机制、低血压或其它病状，从而可适用于患有缺血性中风、昏厥、直立耐受不能、外伤或慢性神经病症的患者。

ii.评估斑块脆弱性

计算分析还可提供定量作用于斑块的患者特异性生物力学力的结果，所述斑块可在患者的颅内和颅外动脉(例如颈动脉粥样硬化斑块)中积累。生物力学力可由脉动压力、流量和颈部运动引起。

图42是根据一个示例性实施方案展示方法1200的方面的示意图，所述方法1200用于提供有关评估特定患者的斑块脆弱性、大脑容量风险和大脑灌注风险的各种信息。方法1200可在上述计算机系统中实施，所述计算机系统例如类似于用以实施上文描述并在图3中展示的步骤中一个或多个的计算机系统。方法1200可使用一个或多个输入1202来进行，并且可包括依据输入1202来产生一个或多个模型1210、依据一个或多个模型1210进行一种或多种生物力学分析1220，以及依据模型1220和生物力学分析1210提供各种结果。

输入1202可包括患者的颅内和颅外动脉(例如患者的主动脉、颈动脉(图37中所示)、脊椎动脉(图37中所示)和脑部)的医学成像数据1203，如ccta数据(例如以与上文关于图2的步骤100所述类似的方式获得)。输入1202还可包括从患者测量的其它生理学数据1204，如患者的臂血压、心率和/或其它测量值(例如以与上文关于图2的步骤100所述类似的方式获得)。其它生理学数据1204可以非侵入方式获得。输入1202可用以产生模型1210和/或进行下述生物力学分析1220。

如上所述，一个或多个模型1210可依据输入1202来产生。举例来说，方法1200可包括产生血液动力学模型1212，所述血液动力学模型1212包括在整个患者解剖结构的三维几何模型中各种位置的计算血液流量和压力信息。患者解剖结构的模型可使用医学成像数据1203创建，并且在一个示例性实施方案中，血液动力学模型1212可为模拟血压模型、模拟血流模型或在进行计算分析之后产生的其它模拟，例如上文关于图3的步骤402所述。固体力学模型(包括流体结构互动模型)可以用已知数值方法的计算分析求解。斑块和血管的性质可建模成线性或非线性、各向同性或各向异性。解可提供斑块以及斑块和血管之间的界面的应力和应变。用于产生血液动力学模型1212的步骤可类似于用于产生上述图35的血液动力学模型932的步骤。

方法1200可包括使用血液动力学模型1212如下来进行生物力学分析1220：计算在各种生理状态(如休息、不同程度的运动或用力等)下因血液动力学力而作用于斑块腔面的压力和剪应力(步骤1222)。压力和剪应力可依据来自血液动力学模型1212信息(例如血液压力和流量)来计算。步骤1222可类似于上述图35的步骤942。

任选地，方法1200还可包括以与上文关于图35的几何分析模型934所述类似的方式，产生用于从四维成像数据定量血管变形的几何分析模型，所述四维成像数据例如是在心动周期的多个阶段(如收缩阶段和舒张阶段)获得的成像数据。方法1200还可包括以与上文关于图35的步骤944所述类似的方式，使用几何分析模型如下进行生物力学分析1220：计算患者的颅内和颅外动脉和斑块由心脏诱导性脉动压力造成的各种变形特征，如纵向变长(伸长)或缩短、扭曲(扭转)、径向膨胀或压缩，和弯曲等。

方法1200还可包括产生用于从医学成像数据1203确定斑块组成和性质的斑块模型1214。举例来说，斑块模型1214可包括关于斑块的密度和其它物质性质的信息。

方法1200还可包括产生血管壁模型1216，所述血管壁模型1216是用于计算关于斑块、血管壁和/或斑块和血管壁之间的界面的信息。举例来说，血管壁模型1216可包括关于应力和应变的信息，所述信息可依据斑块模型1214中包括的斑块组成和性质和在步骤1220中计算的压力和剪应力来计算。任选地，应力和应变还可使用如上所述计算的变形特征来计算。用于产生斑块模型1214和/或血管壁模型1216的步骤可类似于上述图35的用于产生斑块模型936和/或血管壁模型938的步骤。

方法1200可包括使用血管壁模型1216如下进行生物力学分析1220：计算血液动力学力和颈部活动诱导性应变造成的斑块上应力(例如急性应力或积累应力)(步骤1224)。举例来说，可计算作用于斑块的流量诱导力904(图33)。由血液动力学力和颈部活动诱导性应变造成的斑块上应力或力可依据来自血管壁模型1216的信息(例如斑块上应力和应变)来计算。步骤1224可类似于上述图35的步骤946。

方法1200可包括依据上述模型1210中一个或多个和生物力学分析1220中一个或多个来确定其它信息。

可计算斑块破裂脆弱性指数(步骤1230)。斑块破裂脆弱性指数可例如依据血液动力学应力、应力频率、应力方向和/或斑块强度或其它性质来计算。举例来说，包围目标斑块的区域可与斑块的三维模型1210(如斑块模型1214)分离。斑块强度可从斑块模型1214中提供的物质性质确定。由脉动压力、流量和颈部运动造成的目标斑块上的血液动力学应力和组织应力可在模拟基线和运动(或用力)条件下通过使用以下来计算：在步骤1224中预先计算的血液动力学应力和运动诱导性应变。斑块脆弱性可依据斑块应力与斑块强度的比率来评估。步骤1230可类似于上述图35的步骤950。举例来说，可计算位于颅外动脉中的斑块的斑块破裂脆弱性指数以用于中风评估。

还可计算大脑容量风险指数(cvri)(步骤1232)。cvri可定义为在动脉树中给定位置，血管受斑块破裂或栓塞和梗塞(闭塞或阻塞)影响的总脑容量的百分比。cvri可依据给定斑块下游血管提供的脑部的部分来计算，这可相对于下游血管的尺寸来考虑斑块尺寸，以及依据三维血液动力学解，斑块可流入不同血管中的概率。cvri可在患病状态评估，或在介入之前或之后评估。步骤1232可类似于上述图35的步骤952。

脑组织可被建模并分成由血液动力学模拟中的各血管提供的片段(例如关联图40的步骤1110和步骤1112所述)。几何模型可被修改为包括大脑树中分支的下一代(例如关联图40的步骤1118所述)，并且脑组织可被进一步分割(例如关联图40的步骤1120所述)。可在子段中创建其它分支，并且子段可进一步分割成较小片段(例如关联图40的步骤1122所述)。如先前描述的生理学关系可用以关联血管尺寸与成比例量的所提供的脑组织。

可确定破裂斑块所行进的可能路径。血液动力学解可用以确定斑块片段或栓塞物可流入不同下游血管中的机会百分比。

可将破裂斑块的尺寸与下游血管的尺寸相比以确定斑块可在哪里最终产生流障碍。这个信息可与脆弱性指数组合以提供可能受破裂斑块影响的脑组织的容量的概率图。cvri在各可能的受影响的片段都可被赋值。

还可计算大脑灌注风险指数(cpri)(步骤1234)。cpri可定义为在动脉树中给定位置，受斑块破裂和血管梗塞影响的总脑血流量的百分比。cpri指示向脑组织片段的可能灌注损失，而非如cvri指示的受影响容量。举例来说，颈动脉斑块破裂或栓塞的影响可取决于患者的威利斯环(circleofwillis)(图37中所示)的几何形状而变化，并且因解剖结构的这些差异而可得到不同cvri和cpri值。可计算向脑组织各片段的灌注率，并且灌注损失可依据脆弱性指数、血液动力学解，和斑块和血管的尺寸来计算。cpri可在患病状态评估，或在介入之前或之后评估。步骤1234可类似于上述图35的步骤954。

因此，可评估作用于颈动脉粥样硬化斑块的生物力学力，所述生物力学力由脉动压力、脉动血流和/或任选地颈部运动引起。可定量斑块所经历的由脉动压力、脉动血流和/或任选地颈部运动引起的总应力。解可考虑作用于斑块或作用于斑块和血管壁之间的界面的患者特异性血液动力学应力的多个来源。而且，斑块强度可依据医学成像数据来估计，并且可定量有关斑块脆弱性、大脑容量风险和大脑灌注风险的指数。

通过确定如下文所述的颅外和颅内动脉的解剖学和生理学数据，可预测特定患者的动脉或器官层面在各种身体条件下的血流变化。此外，可提供其它信息，如短暂性缺血性发作、缺血性中风或动脉瘤破裂的风险、作用于动脉粥样硬化斑块或动脉瘤的力、医学介入或手术疗法对颅内或颅外血流、压力、壁应力或脑灌注的预测影响。可定量颅内或颅外动脉中的血流、压力和壁应力，以及总脑灌注和局部脑灌注，并且可确定疾病的功能意义。

除定量从成像数据构建的三维几何模型中的血流(例如上文在步骤1212中描述)以外，模型还可被修改来模拟疾病进展或消退或医药介入、经皮介入或手术介入的影响。在一个示例性实施方案中，动脉粥样硬化进展可如下建模：随时间对解进行迭代，例如通过求出剪应力或颗粒滞留时间，以及依据血液动力学因素和/或患者特异性生物化学测量值来调适几何模型以使动脉粥样硬化斑块发展得到进展。此外，血流、心率、血压和其它生理学变量的变化对颅外和/或颅内动脉血流或大脑灌注的影响可通过改变边界条件来建模，并且用以计算这些变量随时间变化的积累作用。

任何实施方案中阐述的任何方面都可与本文所述的任何其它实施方案一起使用。本文所述的每一设备和装置都可用于任何适合的医学程序，可进展经过任何适合的身体内腔和体腔，并且可用于对任何适合的身体部分进行成像。

将对所属领域的技术人员显而易知的是，各种修改和变化可在公开系统和方法中进行，而不背离公开内容的范畴。考虑本文公开的公开内容的说明书和实践后，所属领域的技术人员将显而易知其它实施方案。旨在将说明书和实施例视为仅具有示例性，公开内容的真实范畴和精神由以下权利要求指示。