CN112419277B - 三维血管中心线合成方法、***及存储介质 - Google Patents

Abstract

本申请提供了三维血管中心线合成方法、***及存储介质，包括：获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息；从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线；根据每幅所述冠状动脉二维造影图像的图像信息，将每条所述二维血管中心线投影到三维空间内，合成所述三维血管中心线。本申请的三维血管中心线合成方法摆脱了对多角度冠状动脉二维造影图像成像时间相同、空间相同的要求，在保证精度的前提下，适用范围更广。

Description

三维血管中心线合成方法、***及存储介质

技术领域

本发明涉及冠状动脉医学技术领域，特别是涉及三维血管中心线合成方法、***及存储介质。

背景技术

人体血液中的脂类及糖类物质在血管壁上的沉积将在血管壁上形成斑块，继而导致血管狭窄；特别是发生在心脏冠脉附近的血管狭窄将导致心肌供血不足，诱发冠心病、心绞痛等病症，对人类的健康造成严重威胁。据统计，我国现有冠心病患者约1100万人，心血管介入手术治疗患者数量每年增长大于10％。

冠脉造影CAG、计算机断层扫描CT等常规医用检测手段虽然可以显示心脏冠脉血管狭窄的严重程度，但是并不能准确评价冠脉的缺血情况。为提高冠脉血管功能评价的准确性，1993年Pijls提出了通过压力测定推算冠脉血管功能的新指标——血流储备分数(Fractional Flow Reserve，FFR)，经过长期的基础与临床研究，FFR已成为冠脉狭窄功能性评价的金标准。

血流储备分数(FFR)通常是指心肌血流储备分数，定义为病变冠脉能为心肌提供的最大血流与该冠脉完全正常时最大供血流量之比，研究表明，在冠脉最大充血状态下，血流量的比值可以用压力值来代替。即FFR值的测量可在冠脉最大充血状态下，通过压力传感器对冠脉远端狭窄处的压力和冠脉狭窄近端压力进行测定继而计算得出。

现有存在的问题在于多角度冠状动脉二维造影图像成像的过程中，大部分的DSA成像设备，在同一时间只能拍同一角度的冠状动脉二维造影图像成像，也就是说，多个不同角度的冠状动脉二维造影图像成像必须由不同时间得到。不符合直接三维血管中心线合成方法对图像的要求，因此合成结果大多失真。

发明内容

本发明提供了一种三维血管中心线合成方法、***及存储介质，以解决三维血管中心线合成要求多角度冠状动脉二维造影图像成像时间和空间必须都相同的问题。

为实现上述目的，第一方面，本申请提供了一种三维血管中心线的合成方法，包括：

获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息；

从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线；

根据每幅所述冠状动脉二维造影图像的图像信息，将每条所述二维血管中心线投影到三维空间内，合成所述三维血管中心线。

可选地，上述的三维血管中心线的合成方法，所述获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息的方法包括：

获取至少两组拍摄角度不同的冠状动脉二维造影图像组；

读取每组所述冠状动脉二维造影图像组的图像信息，包括拍摄角度和探测距离；

根据所述探测距离，分别从每组所述冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。

可选地，上述的三维血管中心线的合成方法，所述根据每幅所述冠状动脉二维造影图像的图像信息，将每条所述二维血管中心线投影到三维空间内，合成所述三维血管中心线的方法包括：

以心脏为坐标原点，建立三维坐标系；

获取每幅所述感兴趣的二维造影图像的左右角度α、前后角度β，以及人体与平板探测器之间的距离S，所述二维血管中心线上各点的坐标为(x，y)；

将各点(x，y)投影到三维空间内，获得一系列三维坐标点P，坐标为(x’、y’、z’)；

将放射源投影到所述三维空间内形成放射点R；

将每个所述三维坐标点P均与所述放射点R进行连线，从PR连线上获得所述三维血管中心线上的点，将所述三维血管中心线上的点依次连接获得所述三维血管中心线。

可选地，上述的三维血管中心线的合成方法，所述将各点(x，y)投影到三维空间内，获得一系列三维坐标点P，坐标为(x’、y’、z’)的方法，包括：

将各点(x，y)绕y轴旋转获得(x’、y’、z’)系列点，具体公式为：

将(x’、y’、z’)系列点绕x轴旋转获得一系列三维坐标点P，坐标为(x’、y’、z’)；

可选地，上述的三维血管中心线的合成方法，所述将将放射源投影到所述三维空间内形成放射点R的方法，包括：

获得每幅所述感兴趣的二维造影图像中人体与放射源之间的距离S’；

根据公式

其中，R在三维空间内的坐标为(a，b，c)。

可选地，上述的三维血管中心线的合成方法，所述将每个所述三维坐标点P均与所述放射点R进行连线，从PR连线上获得所述三维血管中心线上的点，将所述三维血管中心线上的点依次连接获得所述三维血管中心线的方法，包括：

将从同一幅所述感兴趣的二维造影图像中获得的一系列三维坐标点P和放射点R对应连接，获得多条PR直线；

获取所述血管的同一位置的两条PR直线之间的最小距离的点，分别为A点和B点；

将所述A点与所述B点连线，获取所述AB线段的中点作为所述三维血管的中心线上的点；

将获得的一系列的所述三维血管的中心线上的点依次连接，获得所述三维血管中心线。

可选地，上述的三维血管中心线的合成方法，所述分别从每幅所述冠状动脉二维造影图像中提取一条二维血管中心线的方法包括：

读取冠状动脉二维造影图像；

获取感兴趣的血管段；

拾取所述感兴趣的血管段的起始点、种子点和结束点；

分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割，得到至少两个局部血管区域图；

从每个所述局部血管区域图中提取至少一条血管局部路径线；

将每个所述局部血管区域图上相对应的血管局部路径线连接，获得至少一条所述血管路径线；

选取一条所述血管路径线作为所述二维血管中心线。

可选地，上述的三维血管中心线的合成方法，所述从每个所述局部血管区域图上提取至少一条血管局部路径线的方法包括：

对所述局部血管区域图做图像增强处理，得到对比强烈的粗略血管图；

对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线。

可选地，上述的三维血管中心线的合成方法，所述对所述局部血管区域图做图像增强处理，得到对比强烈的粗略血管图的方法，包括：

在每幅所述局部血管区域图中，以所述感兴趣的血管段作为前景，其他区域作为背景，强化所述前景，弱化所述背景，得到对比强烈的所述粗略血管图。

可选地，上述的三维血管中心线的合成方法，所述对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线的方法包括：

对所述粗略血管图进行网格划分；

沿着所述起始点至所述结束点的血管延伸方向，搜索所述起始点与周边n个网格上的交叉点的最短时间路径作为第二个点，搜索所述第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点，所述第三个点重复上述步骤，直至最短时间路径到达结束点，其中，n为大于等于1的正整数；

按照搜索顺序，从所述起始点至所述结束点的血管延伸方向连线，获得至少一条血管局部路径线。

可选地，上述的三维血管中心线的合成方法，所述选取一条所述血管路径线作为所述二维血管中心线的方法包括：

如果血管路径线为两条或两条以上，则对每条血管路径线从所述起始点至所述结束点所用的时间求和；

取用时最少的所述血管路径线作为所述二维血管中心线。

第二方面，本申请提供了一种三维血管中心线合成***，用于上述的三维血管中心线的合成方法，包括：图像读取装置，与所述图像读取装置连接的二维血管中心线提取装置和所述三维血管中心线获取装置，所述二维血管中心线提取装置与所述三维血管中心线获取装置连接；

所述图像读取装置，用于获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息；

所述二维血管中心线提取装置，用于接收所述图像读取装置发送的冠状动脉二维造影图像，从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线；

所述三维血管中心线获取装置，用于接收所述二维血管中心线提取装置发送的二维血管中心线，以及接收所述图像读取装置发送的冠状动脉二维造影图像和图像信息，根据每幅所述冠状动脉二维造影图像的图像信息，将每条所述二维血管中心线投影到三维空间内，合成所述三维血管中心线。

可选地，上述的三维血管合成***，所述图像读取装置包括：依次连接的二维图像获取结构、图像信息获取结构、筛选感兴趣图像结构，所述筛选感兴趣图像结构与所述二维图像获取结构连接；

所述二维图像获取结构，用于获取至少两组拍摄角度不同的冠状动脉二维造影图像组；

所述图像信息获取结构，用于读取所述二维图像获取结构传递的每组所述冠状动脉二维造影图像组的图像信息，包括拍摄角度和探测距离；

所述筛选感兴趣图像结构，用于根据所述探测距离，分别从每组所述冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。

第三方面，本申请提供了一种冠状动脉分析***，包括：上述的三维血管合成***。

第四方面，本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的三维血管的合成方法。

本申请实施例提供的方案带来的有益效果至少包括：

本申请提供了一种三维血管的合成方法，摆脱了对多角度冠状动脉二维造影图像成像时间相同、空间相同的要求，在保证精度的前提下，适用范围更广。该方法适应于现有部署在实际场景中的成像设备。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1为本申请的三维血管中心线的合成方法的流程图；

图2为本申请的S100的流程图；

图3为本申请的S200的流程图；

图4为本申请的S260的流程图；

图5为本申请的S262的流程图；

图6为本申请的S270的流程图；

图7为本申请的S300的流程图；

图8为本申请的S350的流程图；

图9为本申请的三维血管中心线合成***的结构框图；

图10为本申请的图像读取装置100的结构框图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

以下将以图式揭露本发明的多个实施方式，为明确说明起见，许多实务上的细节将在以下叙述中一并说明。然而，应了解到，这些实务上的细节不应用以限制本发明。也就是说，在本发明的部分实施方式中，这些实务上的细节是非必要的。此外，为简化图式起见，一些习知惯用的结构与组件在图式中将以简单的示意的方式绘示之。

实施例1：

如图1所示，本申请为了解决上述问题，本申请提供了一种三维血管中心线的合成方法，包括：

S100，获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息；

S200，从每幅感兴趣的二维造影图像中提取一条二维血管中心线；

S300，根据每幅冠状动脉二维造影图像的图像信息，将每条二维血管中心线投影到三维空间内，合成三维血管中心线。

实施例2：

本申请为了解决上述问题，本申请提供了一种三维血管中心线的合成方法，包括：

S100，如图2所示，获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息，包括：

S110，获取至少两组拍摄角度不同的冠状动脉二维造影图像组；

S120，读取每组冠状动脉二维造影图像组的图像信息，包括拍摄角度和探测距离；

S130，根据探测距离，分别从每组冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。

S200，如图3所示，从每幅感兴趣的二维造影图像中提取一条二维血管中心线，包括：

S210，读取冠状动脉二维造影图像；

S220，获取感兴趣的血管段；

S230，拾取感兴趣的血管段的起始点、种子点和结束点；

S240，分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割，得到至少两个局部血管区域图；

S250，从每个局部血管区域图中提取至少一条血管局部路径线；

S260，如图4所示，将每个局部血管区域图上相对应的血管局部路径线连接，获得至少一条血管路径线，包括：

S261，对局部血管区域图做图像增强处理，得到对比强烈的粗略血管图，包括：在每幅局部血管区域图中，以感兴趣的血管段作为前景，其他区域作为背景，强化前景，弱化背景，得到对比强烈的粗略血管图。

S262，如图5所示，对粗略血管图做网格划分，沿着起始点至结束点方向，提取至少一条血管局部路径线，包括：

S2621，对粗略血管图进行网格划分；

S2622，沿着起始点至结束点的血管延伸方向，搜索起始点与周边n个网格上的交叉点的最短时间路径作为第二个点，搜索第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点，第三个点重复上述步骤，直至最短时间路径到达结束点，其中，n为大于等于1的正整数；

S2623，按照搜索顺序，从起始点至结束点的血管延伸方向连线，获得至少一条血管局部路径线。

S270，如图6所示，选取一条血管路径线作为二维血管中心线，包括：

S271，如果血管路径线为两条或两条以上，则对每条血管路径线从起始点至结束点所用的时间求和；

S272，取用时最少的血管路径线作为二维血管中心线。

S300，如图7所示，根据每幅冠状动脉二维造影图像的图像信息，将每条二维血管中心线投影到三维空间内，合成三维血管中心线，包括：

S310，以心脏为坐标原点，建立三维坐标系；

S320，获取每幅感兴趣的二维造影图像的左右角度α、前后角度β，以及人体与平板探测器之间的距离S，二维血管中心线上各点的坐标为(x，y)；

S330，将各点(x，y)投影到三维空间内，获得一系列三维坐标点P，坐标为(x’、y’、z’)，包括：

将各点(x，y)绕y轴旋转获得(x’、y’、z’)系列点，具体公式为：

将(x’、y’、z’)系列点绕x轴旋转获得一系列三维坐标点P，坐标为(x’、y’、z’)；

S340，将放射源投影到三维空间内形成放射点R，包括：

获得每幅感兴趣的二维造影图像中人体与放射源之间的距离S’；

根据公式

其中，R在三维空间内的坐标为(a，b，c)。

S350，如图8所示，将每个三维坐标点P均与放射点R进行连线，从PR连线上获得三维血管中心线上的点，将三维血管中心线上的点依次连接获得三维血管中心线，包括：

S351，将从同一幅感兴趣的二维造影图像中获得的一系列三维坐标点P和放射点R对应连接，获得多条PR直线；

S352，获取血管的同一位置的两条PR直线之间的最小距离的点，分别为A点和B点；

S353，将A点与B点连线，获取AB线段的中点作为三维血管的中心线上的点；

S354，将获得的一系列的三维血管的中心线上的点依次连接，获得三维血管中心线。

如图9所示，本申请提供了一种三维血管中心线合成***，用于上述的三维血管中心线的合成方法，包括：图像读取装置100，与图像读取装置100连接的二维血管中心线提取装置200和三维血管中心线获取装置300，二维血管中心线提取装置200与三维血管中心线获取装置300连接；图像读取装置100，用于获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息；二维血管中心线提取装置200用于接收图像读取装置发送的冠状动脉二维造影图像，从每幅感兴趣的二维造影图像中提取一条二维血管中心线；三维血管中心线获取装置300用于接收二维血管中心线提取装置200发送的二维血管中心线，以及接收图像读取装置100发送的冠状动脉二维造影图像和图像信息，根据每幅冠状动脉二维造影图像的图像信息，将每条二维血管中心线投影到三维空间内，合成三维血管中心线。

如图10所示，本申请的一个实施例中，图像读取装置100包括：依次连接的二维图像获取结构110、图像信息获取结构120、筛选感兴趣图像结构130，筛选感兴趣图像结构130与二维图像获取结构110连接；二维图像获取结构110用于获取至少两组拍摄角度不同的冠状动脉二维造影图像组；图像信息获取结构120用于读取二维图像获取结构传递的每组冠状动脉二维造影图像组的图像信息，包括拍摄角度和探测距离；筛选感兴趣图像结构130用于根据探测距离，分别从每组冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。

本申请提供了一种冠状动脉分析***，包括：上述的三维血管合成***。

本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的三维血管的合成方法。

所属技术领域的技术人员知道，本发明的各个方面可以实现为***、方法或计算机程序产品。因此，本发明的各个方面可以具体实现为以下形式，即：完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等)，或硬件和软件方面结合的实施方式，这里可以统称为“电路”、“模块”或“***”。此外，在一些实施例中，本发明的各个方面还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式，该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或***的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。

例如，可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件，可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作***执行的多个软件指令。在本发明的示例性实施例中，由数据处理器来执行如本文的根据方法和/或***的示例性实施例的一个或多个任务，诸如用于执行多个指令的计算平台。可选地，该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器，例如，磁硬盘和/或可移动介质。可选地，也提供了一种网络连接。可选地也提供显示器和/或用户输入设备，诸如键盘或鼠标。

可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的***、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项：

具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行***、装置或者器件使用或者与其结合使用。

计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行***、装置或者器件使用或者与其结合使用的程序。

计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括(但不限于)无线、有线、光缆、RF等等，或者上述的任意合适的组合。

例如，可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码，包括诸如Java、Smalltalk、C++等面向对象编程语言和常规过程编程语言，诸如"C"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络--包括局域网(LAN)或广域网(WAN)-连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。

应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。

也可以把这些计算机程序指令存储在计算机可读介质中，这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作，从而，存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(article of manufacture)。

还可将计算机程序指令加载到计算机(例如，冠状动脉分析***)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程，使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。

本发明的以上的具体实例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (12)

1.一种三维血管中心线的合成方法，其特征在于，包括：

获取至少两组拍摄角度不同的冠状动脉二维造影图像组；

读取每组所述冠状动脉二维造影图像组的图像信息，包括拍摄角度和探测距离；

根据所述探测距离，分别从每组所述冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像；

从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线；

根据每幅所述冠状动脉二维造影图像的图像信息，将每条所述二维血管中心线投影到三维空间内，合成所述三维血管中心线的方法包括：

以心脏为坐标原点，建立三维坐标系；

获取每幅所述感兴趣的二维造影图像的左右角度α、前后角度β，以及人体与平板探测器之间的距离S，所述二维血管中心线上各点的坐标为(x，y)；

将各点(x，y)投影到三维空间内，获得一系列三维坐标点P，坐标为(x’、y’、z’)；

将放射源投影到所述三维空间内形成放射点R；

将每个所述三维坐标点P均与所述放射点R进行连线，从PR连线上获得所述三维血管中心线上的点，将所述三维血管中心线上的点依次连接获得所述三维血管中心线；

所述将各点(x，y)投影到三维空间内，获得一系列三维坐标点P，坐标为(x’、y’、z’)的方法，包括：

将各点(x，y)绕y轴旋转获得(x’、y’、z’)系列点，具体公式为：

将(x’、y’、z’)系列点绕x轴旋转获得一系列三维坐标点P，坐标为(x’、y’、z’)，具体公式为：

2.根据权利要求1所述的一种三维血管中心线的合成方法，其特征在于，所述将放射源投影到所述三维空间内形成放射点R的方法，包括：

获得每幅所述感兴趣的二维造影图像中人体与放射源之间的距离S’；

根据公式

其中，R在三维空间内的坐标为(a，b，c)。

3.根据权利要求2所述的一种三维血管中心线的合成方法，其特征在于，所述将每个所述三维坐标点P均与所述放射点R进行连线，从PR连线上获得所述三维血管中心线上的点，将所述三维血管中心线上的点依次连接获得所述三维血管中心线的方法，包括：

将从同一幅所述感兴趣的二维造影图像中获得的一系列三维坐标点P和放射点R对应连接，获得多条PR直线；

获取所述血管的同一位置的两条PR直线之间的最小距离的点，分别为A点和B点；

将所述A点与所述B点连线，获取AB线段的中点作为所述三维血管的中心线上的点；

将获得的一系列的所述三维血管的中心线上的点依次连接，获得所述三维血管中心线。

4.根据权利要求1所述的一种三维血管中心线的合成方法，其特征在于，所述分别从每幅所述冠状动脉二维造影图像中提取一条二维血管中心线的方法包括：

读取冠状动脉二维造影图像；

获取感兴趣的血管段；

拾取所述感兴趣的血管段的起始点、种子点和结束点；

分别对起始点、种子点、结束点的相邻两点间的二维造影图像进行分割，得到至少两个局部血管区域图；

从每个所述局部血管区域图中提取至少一条血管局部路径线；

将每个所述局部血管区域图上相对应的血管局部路径线连接，获得至少一条所述血管路径线；

选取一条所述血管路径线作为所述二维血管中心线。

5.根据权利要求4所述的一种三维血管中心线的合成方法，其特征在于，所述从每个所述局部血管区域图上提取至少一条血管局部路径线的方法包括：

对所述局部血管区域图做图像增强处理，得到对比强烈的粗略血管图；

对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线。

6.根据权利要求5所述的一种三维血管中心线的合成方法，其特征在于，所述对所述局部血管区域图做图像增强处理，得到对比强烈的粗略血管图的方法，包括：

在每幅所述局部血管区域图中，以所述感兴趣的血管段作为前景，其他区域作为背景，强化所述前景，弱化所述背景，得到对比强烈的所述粗略血管图。

7.根据权利要求6所述的一种三维血管中心线的合成方法，其特征在于，所述对所述粗略血管图做网格划分，沿着所述起始点至所述结束点方向，提取至少一条血管局部路径线的方法包括：

对所述粗略血管图进行网格划分；

沿着所述起始点至所述结束点的血管延伸方向，搜索所述起始点与周边n个网格上的交叉点的最短时间路径作为第二个点，搜索所述第二个点与周边n个网格上的交叉点的最短时间路径作为第三个点，所述第三个点重复上述步骤，直至最短时间路径到达结束点，其中，n为大于等于1的正整数；

按照搜索顺序，从所述起始点至所述结束点的血管延伸方向连线，获得至少一条血管局部路径线。

8.根据权利要求7所述的一种三维血管中心线的合成方法，其特征在于，所述选取一条所述血管路径线作为所述二维血管中心线的方法包括：

如果血管路径线为两条或两条以上，则对每条血管路径线从所述起始点至所述结束点所用的时间求和；

取用时最少的所述血管路径线作为所述二维血管中心线。

9.一种三维血管中心线合成***，用于权利要求1～8任一项所述的三维血管中心线的合成方法，其特征在于，包括：图像读取装置，与所述图像读取装置连接的二维血管中心线提取装置和所述三维血管中心线获取装置，所述二维血管中心线提取装置与所述三维血管中心线获取装置连接；

所述图像读取装置，用于获取至少两幅拍摄角度不同的冠状动脉二维造影图像的图像信息；

所述二维血管中心线提取装置，用于接收所述图像读取装置发送的冠状动脉二维造影图像，从每幅所述感兴趣的二维造影图像中提取一条二维血管中心线；

所述三维血管中心线获取装置，用于接收所述二维血管中心线提取装置发送的二维血管中心线，以及接收所述图像读取装置发送的冠状动脉二维造影图像和图像信息，根据每幅所述冠状动脉二维造影图像的图像信息，将每条所述二维血管中心线投影到三维空间内，合成所述三维血管中心线。

10.根据权利要求9所述的一种三维血管中心线合成***，其特征在于，所述图像读取装置包括：依次连接的二维图像获取结构、图像信息获取结构、筛选感兴趣图像结构，所述筛选感兴趣图像结构与所述二维图像获取结构连接；

所述二维图像获取结构，用于获取至少两组拍摄角度不同的冠状动脉二维造影图像组；

所述图像信息获取结构，用于读取所述二维图像获取结构传递的每组所述冠状动脉二维造影图像组的图像信息，包括拍摄角度和探测距离；

所述筛选感兴趣图像结构，用于根据所述探测距离，分别从每组所述冠状动脉二维造影图像中选取一幅感兴趣的二维造影图像。

11.一种冠状动脉分析***，其特征在于，包括：权利要求9～10任一项所述的一种三维血管中心线合成***。

12.一种计算机存储介质，其特征在于，计算机程序被处理器执行时实现权利要求1～8任一项所述的一种三维血管中心线的合成方法。