血流储备分数模拟参数的定制、校准和/或训练的制作方法

以下总体涉及成像，并且更具体而言涉及血流储备分数模拟(ffr)的定制、校准和/或训练，并且针对ffr-计算机扫描断层(ffr-ct)的特定应用进行了描述，并且也适用于x射线ffr。

背景技术：

血流储备分数(ffr)是导管实验室(cathlab)中的一项侵入性测量，通过ffr指数来量化由于钙化或软斑块引起的冠状动脉病变的血流动力学显著性。该指标指示冠状支架的功能严重性，其根据在冠状血管造影术期间进行的压力测量被计算并且被定义为在充血条件下远端血压(狭窄后面)相比于近端血压(靠近心门)。也就是，ffr指标表达在存在支架时沿血管向下的最大流量相比于假设没有狭窄时的最大流量。ffr值是介于0与1之间的绝对值，其中，值0.50指示给定的狭窄会导致血压下降50％。

侵入性ffr流程要求将导管插入股动脉或桡动脉中，并且将导管推进到狭窄处，其中在由影响血管几何结构、顺应性和阻力和/或其他特性的各种药剂促进的情况期间，在导管尖端的传感器感测在狭窄两端的压力。非侵入性方法通过计算流体动力学(cfd)模拟来根据心脏的ct图像数据(例如，根据造影增强的计算机断层摄影冠状动脉血管造影(ccta))来估计ffr指标，在该模拟中，模拟了通过冠状动脉的血流和血压。这包括使用ccta图像数据以导出冠状动脉树的几何模型、从其中提取特征并且根据特征来确定针对该模拟的边界条件。

ffr-ct模拟软件已经在临床站点部署，具有预先设定的、工厂经调整的算法参数(例如，用于根据ccta图像来分割冠状动脉)，或针对控制通过冠状动脉的模拟血流的外周阻力的比例因子。然而，使用模式、个人习惯和/或患者人群可能因站点而异。例如，当分割血管时，一个临床医师可以在血管壁的内表面内进行分割，而另一临床医师则在血管壁上进行分割。在另一示例中，相对于较低的图像质量的图像，在较高的图像质量的图像中血管壁的实际边缘可以更清晰。ffr-ct模拟结果的准确性可以因临床医师而异、因站点而异等等。

技术实现要素：

本文所描述的方面解决了上面提到的问题和其他问题。

在一个方面中，一种计算系统包括具有计算机可执行指令的计算机可读存储介质，所述计算机可执行指令包括具有分割器和边界条件确定器的生物物理模拟器。所述计算系统还包括处理器，所述处理器被配置为执行所述生物物理模拟器以利用经适配的冠状动脉树分割和经适配的边界条件中的至少一个以及心脏成像数据来计算血流储备分数指标。

在另一个方面中，一种计算机可读存储介质编码有计算机可读指令，所述计算机可读指令在由计算系统的处理器执行时，使所述处理器进行以下操作：接收心脏成像数据，并且执行生物物理模拟器以利用经适配的冠状动脉树分割和经适配的边界条件中的至少一个以及心脏成像数据来计算血流储备分数指标。

在另一个方面中，一种方法包括：接收心脏成像数据，并且利用经适配的冠状动脉树分割和经适配的边界条件中的至少一个以及心脏成像数据来计算血流储备分数指标。

本领域技术人员通过阅读并理解所附的描述，将认识到本申请的其他方面。

附图说明

本发明可以采取各种部件和部件的布置以及各种步骤和步骤布置的形式。附图仅出于图示优选实施例的目的，并且不应被解释为限制本发明。

图1示意性地图示了系统，其包括计算系统和成像系统。

图2图示了示例生物物理模拟器。

图3图示了冠状动脉的示例分割。

图4图示了冠状动脉的示例解剖模型。

图5图示了根据本文的实施例的示例方法。

图6图示了根据本文的实施例的另一示例方法。

图7图示了根据本文的实施例的另一示例方法。

具体实施方式

图1示意性地示出了包括成像系统102(例如，ct扫描器)的系统100。在变型中，系统100包括x射线成像器。成像系统102包括体固定机架104和旋转机架106，旋转机架104由固定机架104可旋转地支撑并且围绕检查区域108关于z-轴旋转。对象支撑件110，例如卧榻，在检查区域108中支撑目标或对象。

辐射源112，例如x-射线管，由旋转机架106可旋转地支撑，随旋转机架106旋转，并且发出穿过检查区域108的辐射。辐射敏感探测器阵列114对向一角度弧，跨检查区域1088与辐射源112相对。辐射敏感探测器阵列114探测穿过检查区域108的辐射并生成指示其的(一个或多个)电信号(投影数据)。重建器116重建投影数据，生成指示定位于检查区域108中的对象或目标的被扫描部分的体积图像数据。

在该示例中，计算系统118用作操作者控制台。控制台118包括处理器120(例如，微处理器、中央处理单元等)和计算机可读存储介质122，其排除瞬态介质，并且包括诸如物理存储设备之类的瞬态介质，等等。控制台118还包括(一个或多个)人类可读输出设备(例如显示监视器)和(一个或多个)输入设备(例如鼠标、键盘等)。

计算机可读存储介质122包括用于生物物理模拟器126的指令124。处理器120被配置为执行指令124和/或允许操作者经由图形用户界面(gui)或其他方式与扫描器102交互和/或操作扫描器102的软件。处理器120可以额外地或备选地执行由载波、信号和/或其他暂态介质承载的(一条或多条)计算机可读指令。

在变型中，生物物理模拟器126是另一计算系统的一部分，其与控制台118和系统100分离。在这种情况下，另一计算系统与控制台118的相似之处在于它包括处理器、计算机可读存储介质、输入设备和输出设备，但是它不包括允许操作者与扫描器102交互和/或操作扫描器102的软件。

生物物理模拟器126被配置为至少处理体积图像数据以执行生物物理模拟。关于ffr，生物物理模拟器确定基于ffr指标的ccta图像数据(或x射线血管造影)。ffr指标可以经由显示监视器显示、存储、传送到另一设备等。如下面更详细描述的，生物物理模拟器126包括一个或多个反馈回路，例如，分割反馈回路和/或边界条件反馈回路。反馈回路可以用于针对个体临床医师、站点等的训练、模拟定制和/或校准和/或预测的目的。

因此，生物物理模拟器126可以被部署在临床站点并且针对个体用户和/或站点进行调谐，例如，用于根据ccta图像数据(或x射线血管造影)来分割冠状动脉树，或针对控制通过冠状动脉的模拟血流的外周阻力的比例因子。在一个实例中，这减轻了因站点而异的使用模式、个人习惯和/或患者人群的不同结果。因此，相对于其中没有利用本文所描述的方法和/或省略了本文所描述的生物物理模拟器126的配置，可以改进由不同的临床医师、不同的站点等导致的ffr-ct(或x射线ffr)模拟的准确性。

图2示意性地示出了结合ffr-ct的生物物理模拟器126的示例。在该示例中，生物物理模拟器126包括具有分割反馈回路204的分割器202、具有边界条件反馈回路208的边界条件确定器206和流动模拟器210。在变型中，省略了分割反馈回路204。在另一变型中，省略了边界条件反馈回路208。生物物理模拟器126从成像系统100、数据存储库(例如，放射信息系统(ris)、图片存档及通信系统(pacs)等)和/或其他装置接收ccta成像数据作为输入。

分割器202采用分割算法来根据ccta成像数据分割冠状动脉树。可以自动地(例如，机器学习等)或半自动地(例如，利用用户协助)执行分割。在一个实例中，分割包括从其识别和/或提取冠状动脉中心线和/或管腔几何结构(例如，直径、周长、横截面积等)。分割可以基于体素强度、物体形状和/或其他特性。图3示出了个体血管的一部分的分割，示出了血管腔的相对血管壁302，并且图4示出了经分割的冠状动脉树400。

在以下文章中讨论了用于根据ccta成像数据来提取冠状动脉树的合适方法的示例：zheng等人的“robustandaccuratecoronaryarterycenterlineextractioninctabycombiningmodel-drivenanddata-drivenapproaches”medimagecomputassistinterv，2013；16(pt3)：74-81；ecabert等人的“segmentationoftheheartandgreatvesselsinctimagesusingamodel-basedadaptationframework”，medimageanal，2011dec；15(6)：863-76以及freiman等人的“improvingccta-basedlesions’hemodynamicsignificanceassessmentbyaccountingforpartialvolumemodelinginautomaticcoronarylumensegmentation”，medphys，2017mar；44(3)：1040-1049。本文也预期其他方法。

在一个实例中，ccta成像数据是训练集，并且专家们已经执行了准确的冠状动脉树分割，所述分割导致与已知的“真实情况”侵入性ffr测量相匹配的对ffr指标的模拟，从而为训练集提供了参考分割。分割器202根据训练集来分割冠状动脉树，并且用户利用分割器202的工具自由地操纵冠状动脉树分割。分割器202将经用户调整的冠状动脉树分割与参考分割进行比较。分割器202经由反馈回路204提供反馈，所述反馈指示经用户调整的冠状动脉树分割与参考分割之间的任何差异。反馈可以指示随后的经调整的冠状动脉树分割是否更接近参考分割、用于改进分割的推荐和/或关于分割的其他信息。

在另一实例中，ccta成像数据再次是训练集。在该实例中，分割器202学习经用户调整的冠状动脉树分割与参考分割之间的差异，并且将该偏差存储为用户特异的校准数据。分割器202利用用户特异的校准数据，例如以自动地根据处于评估或分析中的患者ccta成像数据(非训练数据集)来适配随后的经用户调整的冠状动脉树分割。为此，分割器202经由反馈回路204基于用户特异的校准数据视觉地呈现经适配的经调整的冠状动脉树分割。为了适配，分割器202可以移动分割的一个或多个点。例如，分割器202可以增加或减小所分割的冠状动脉血管的直径、移动所分割的冠状动脉血管壁等。

用户可以接受、拒绝和/或修改视觉呈现的经适配经调整的冠状动脉树分割。用户还可以重做分割。用户可以定期地使用训练集(和/或另一训练集)来更新其用户特异的校准数据。例如，随着时间的推移，用户相对于参考创建更为准确的经调整的冠状动脉树分割的能力可以随着经验的增加而增加。在这种实例中，当前的用户特异校准数据可能会过度校正。为此，用户可以处理训练集(和/或其他训练集)，并且分割器202可以基于由此的偏差来更新用户特异的校准数据。额外地或者替代地，分割器202可以在启动时或以其他方式向用户呈现训练集以更新其用户特异的校准数据。

边界条件确定器206根据经用户调整的冠状动脉树分割和/或分割器202经适配的经用户调整的冠状动脉树分割来确定用于血管中血流的计算流体动力学模拟的边界条件。其中一种方法采用参数集总模型。该模型包括使用非线性阻力的中心线表示，其中元素指示流入边界条件和流出边界条件，并且元素表示树分割传递函数，其中包括反映血管几何结构(例如直径、周长、横截面面积等)和/或水力效应(例如，表1中的1-6)的一系列线性和非线性阻力元素。

表1.水力效应，其中a＝血管的局部横截面面积，p＝血管的周长，l＝血管的长度，r＝血管的半径，以及κ＝血管的曲率，ρ＝血液的密度，并且μ＝血液的粘度。

在nickisch等人的“learningpatient-specificlumpedmodelsforinteractivecoronarybloodflowsimulations”(medicalimagecomputingandcomputer-assistedintervention——miccai2015：第18届国际会议，lncs，vol.9350，2015，vol.9350，第433-441页)中讨论了集总模型的示例。导出边界条件的示例在于2014年6月30日提交的、题为“enhancedpatient’sspecificmodellingforffr-ct”的ep14174891.3中进行了描述，其整体通过引用的方式并入本文。本文也预期其他方法。

在所示的示例中，在ffr-ct模拟之后收集的患者特异的数据被提供给边界条件确定器206。这种数据的非限制性示例包括但不限于血球容积、存在/没有糖尿病、急性冠状动脉综合征和血压。数据可以由用户提供、从电子患者记录中提取和/或以其他方式取回。经由反馈回路208使用患者特异的数据来更新边界条件。例如，针对先前已经接受过ffr-ct模拟的患者进行的侵入性ffr测量结果可以用于改进未来的模拟。此外，已记录的结果、附加的生物物理和/或功能测量可以用于重新调整边界条件。

以下提供了如何利用患者特异的数据的非限制性示例。血球容积是红细胞体积百分比的量度。通常，血球容积升高意味着血液粘度更高。可以使用经验曲线来调整血液的粘度，例如，在40％-45％之间的线性伸缩。糖尿病的存在通常意味着更坚硬的血管壁和更高的心肌阻力。可以例如通过将默认值增加10％来调整对应的阻力边界条件。在处置的临界ffr阈值为0.80的情况下，对于患有急性冠状动脉综合征(acs)的患者，可以将该阈值提高至例如0.85，以确保对边缘患者进行治疗。

可以使用针对特定患者的额外的侵入性参考ffr测量来适配任何模型参数(例如，心肌阻力)，以便(在给定当前的分割的情况下)获得参考与模型预测之间更好的匹配。这可以通过梯度下降、参数搜索等完成。考虑到先前获得的模型预测和药物处置，可以以类似的方式使用已记录的患者结果(例如，心脏事件或生存数据)来适配参数。例如，如果经由ct-ffr评估认为病变是不显著的，并且该病变引起了大的心脏事件，则可以重新考虑该评估，并且可以更新诸如ffr阈值的参数以更好地匹配结果。

额外地或者替代地，对于一个或多个输入参数，边界条件确定器206还针对输入参数中的一个或多个预定变化来确定边界条件。例如，对于所测量的血压x，边界条件确定器206还可以针对至少为(x)(±0.z)的血压值确定边界条件，其中z表示作为百分比的公差，例如，1％、2％、5％、10％、25％等。这允许计算系统118计算输入参数的ffr值以及输入值周围的(一个或多个)值。计算出的ffr值将指示参数值如何影响ffr-ct模拟的结果。

在参数值(例如，所测量的血压)的置信度低，并且值的微小改变对ffr值有显著影响(例如，改变了从观察患者到处置患者的建议)的情况下，用户可以选择获取并使用更准确的值，而不是使用当前值。替代地，用户可以继续了解该参数值如何影响ffr值。替代地，用户可以在不使用该参数值的情况下继续和/或调整受此参数影响的边界条件。在参数值的置信度高的情况下，用户可以继续使用该值。

血流模拟器利用边界条件执行血流模拟，并且生成并输出ffr值。可以例如使用计算流体动力学(cfd)方法和/或其他方法来完成血流模拟。计算ffr值的示例在于2013年5月10日提交的题为“determinationofafractionalflowreserve(ffr)valueforastenosisofavessel”的美国专利us2015/0092999a1，于2013年10月24日提交的题为“fractionalflowreserve(ffr)index”的美国专利us2015/0282765a1中进行了描述，通过引用将其整体并入本文。

图5示出了根据本文所描述的实施例的示例方法。

在502处，根据ccta成像数据的训练集来分割冠状动脉血管。

在504处，接收到指示对分割的调整的用户输入，并将其施加到经分割的血管，从而产生经用户调整的冠状动脉血管分割。

在506处，将所述经用户调整的冠状动脉血管分割与针对所述训练集的参考血管分割进行比较。

在508处，视觉地呈现分割之间的差异。所述差异可以是图形的和/或数字的。

该方法非常适合用于冠状动脉树分割训练。

图6示出了根据本文所描述的实施例的另一示例方法。

在602处，根据ccta成像数据的训练集来分割冠状动脉树。

在604处，接收到指示对分割的调整的用户输入，并将其施加到经分割的树，从而产生经用户调整的冠状动脉树分割。

在606处，将所述经用户调整的冠状动脉树分割与针对心脏成像数据的训练集的参考树分割进行比较。

在608处，将其之间的偏差存储为用户特异的校准。

在610处，根据受评估的患者的ccta成像数据分割冠状动脉树。

在612处，接收到指示对分割的调整的用户输入，并将其施加到经分割的树，从而产生经用户调整的冠状动脉树分割。

在614处，将用户特异的校准施加到经用户调整的冠状动脉血管分割，从而为患者产生经适配的经用户调整的冠状动脉树分割。

在616处，使用所述经适配的经用户调整的冠状动脉树分割来生成ffr指数。

图7示出了根据本文所描述的实施例的另一示例方法。

在702处，根据患者的ccta成像数据来分割冠状动脉树。在一个实例中，如结合图6所描述地适配冠状动脉树。

在704处，基于冠状动脉树来确定边界条件。

在706处，利用边界条件来确定ffr值。

在708处，基于患者特异的数据和/或其他数据来更新边界条件。

以上可以通过编码或嵌入在计算机可读存储介质上的计算机可读指令的方式来实现，所述计算机可读指令当由(多个)计算机处理器执行时，使得(多个)处理器执行所描述的动作。额外地或者替代地，计算机可读指令中的至少一个由不是计算机可读存储介质的信号、载波或其他瞬态介质来承载。此外，应该认识到，动作的排序不是限制性的。因此，本文考虑了其他排序。另外地，可以省略一个或多个动作和/或可以包括一个或多个附加动作。

尽管已经在附图和前述描述中详细说明和描述了本发明，但是这样的说明和描述应被认为是说明性或示例性的而非限制性的；本发明不限于所公开的实施例。本领域技术人员通过研究附图、公开内容和所附权利要求书，在实践所要求保护的发明时可以理解并实现所公开的实施例的其他变型。

在权利要求书中，词语“包括”不排除其他元件或步骤，并且词语“一”或“一个”不排除多个。单个处理器或其他单元可以实现权利要求书中引述的若干项的功能。在互不相同的从属权利要求中引述某些措施的仅有事实并不指示着不能有利地使用这些措施的组合。

计算机程序可以存储在/分发在合适的介质上，例如，与其他硬件一起提供或作为其一部分提供的光学存储介质或固态介质，但是也可以以其他形式发布，例如，经由互联网或其他有线或无线电信系统。权利要求书中的任何附图标记都不应被解释为限制范围。