一种基于X射线冠脉造影图像的多角度血管重建方法与流程

本发明属于数字图像处理与医学图像的交叉领域，可应用于临床研究中x射冠脉造影图像分析，尤其是涉及一种基于x射线冠脉造影图像的多角度血管重建方法。

背景技术：

目前，国内大多数医院都用的是x射线单臂造影系统对病人做x射线造影，通过旋转造影臂得到一个对应于不同造影角度的造影图像序列。单臂造影可以很方便对病人进行不同角度的造影，但是对重建来说个缺点就是我们无法得到同一时刻的不同视角的造影图像，这给重建带来很大的困难。要重建血管的真实三维空间结构，需要得到血管至少两个不同角度的投影信息。传统方法首先提取出血管的骨架，然后通过不同视角空间约束关系，对不同视角投影图像的血管像素点进行正确匹配并重建，当血管几何形变不大，几何关系明显的时候，才能大体上恢复出血管的三维空间结构。

血管树从整体上可以看作在空间中弯曲延展的管状系统，其骨架是具有树状结构的连续空间曲线，反映了血管树的整体形态特征。对于心血管双视角造影图像的匹配问题，传统的方法是通过外极线约束来找两幅造影图像中匹配的血管点对，但是外极线约束极易发生误匹配。这种方法误差较大，在血管缺失或者重叠的情况下，并不能保证双侧造影图像中过的血管段和取样点是一一对应的关系。常见的匹配方法对图像质量要求高，耗时长，在面对不同的血管的个体差异，鲁棒性差。

鉴于传统血管匹配方法中存在的问题，很多研究者都通过在序列图像中选择处于心脏运动周期中同一时刻的造影图像作为重建的参考图像，比如舒张末期。但这种方法依然解决不了这个问题：在两个视角中，同一根血管不同段的伸缩和扭曲不同。

技术实现要素：

本发明目的是：提供一种基于x射线冠脉造影图像的多角度血管重建方法，通过分析两个角度血管的特征点，并依据特征点将整段血管分割成小段，然后依据每个小段之间的对应关系分别做三维合成，从而降低了整段血管直接合成造成的误差。

本发明的技术方案是：一种基于x射线冠脉造影图像的多角度血管重建方法，包括如下步骤：

步骤s1：基于两个不同角度的x射线冠脉造影图像，分别获得两造影图像上血管的2d结构数据及特征点；

步骤s2：根据血管的特征点将其对应的造影图像上的血管进行分段，并得到每段血管的2d结构数据；

步骤s3：将两个不同角度的x射线冠脉造影图像上、且呈映射关系的两分段血管的2d结构数据进行三维重建，并得到该分段血管的3d结构数据；

步骤s4：重复步骤s3直到所有分段血管三维重建完成，再将重建后的分段血管合并，即得到完整的三维血管。

作为优选的技术方案，步骤s1中血管的2d结构数据包括中心线、半径和角度；

作为优选的技术方案，步骤s1血管的特征点包括起始点、结束点、多个狭窄点和多个分叉点。

作为优选的技术方案，步骤s3中呈映射关系的两分段血管进行三维重建的具体方法如下：

步骤a：通过矩阵变换将两分段血管的2d结构数据分别变换到三维空间的一个正方形面上；

步骤b：将两分段血管的2d结构数据所在的正方形面与空间原点分别构成两个四棱锥体，两个四棱锥体相交的空间即为三维血管的实际限定空间；

步骤c：将两分段血管的中心线与空间原点分别构成两个3d空间曲面，两个3d空间曲面相交形成一条相交线，位于步骤b中实际限定空间内的相交线即为三维血管的中心线。

本发明的优点是：

1.本发明通过分析两个角度血管的特征点，并依据特征点将整段血管分割成小段，然后依据每个小段之间的对应关系分别做三维合成，然后再将分段血管合并，得到完整的三维血管，从而降低了整段血管直接合成造成的误差；

2.本发明解决多角度x射线冠脉造影图像的多角度血管匹配问题，为临床医学心血管疾病辅助检测提供更为有效的手段，提高了血管匹配的可靠性和精确度，从而可以有效增大血管造影三维重建的精度。

附图说明

下面结合附图及实施例对本发明作进一步描述：

图1为本发明的流程图；

图2为x射线冠脉造影图像一示意图；

图3为x射线冠脉造影图像二示意图；

图4为对应造影图像一角度的三维血管合成示意图；

图5为对应造影图像二角度的三维血管合成示意图；

图6为x角度的三维血管合成示意图；

图7为y角度的三维血管合成示意图；

图8为z角度的三维血管合成示意图；

图9为血管分叉点对应示意图；

图10为血管狭窄点对应示意图。

具体实施方式

实施例：参照图1所示：本发明基于x射线冠脉造影图像的多角度血管重建方法，包括如下步骤：

步骤s1：基于两个不同角度的x射线冠脉造影图像(参照图2和图3所示)，分别获得两造影图像上血管的2d结构数据及特征点，2d结构数据包括中心线、半径和角度；特征点包括起始点、结束点、多个狭窄点和多个分叉点；

步骤s2：根据血管的特征点将其对应的造影图像上的血管进行分段，并得到每段血管的2d结构数据；

步骤s3：将两个不同角度的x射线冠脉造影图像上、且呈映射关系的两分段血管的2d结构数据进行三维重建，并得到该分段血管的3d结构数据；呈映射关系的两分段血管进行三维重建的具体方法如下：

步骤a：通过矩阵变换将两分段血管的2d结构数据分别变换到三维空间的一个正方形面上；

步骤b：将两分段血管的2d结构数据所在的正方形面与空间原点分别构成两个四棱锥体，两个四棱锥体相交的空间即为三维血管的实际限定空间；

步骤c：将两分段血管的中心线与空间原点分别构成两个3d空间曲面，两个3d空间曲面相交形成一条相交线，位于步骤b中实际限定空间内的相交线即为三维血管的中心线；

步骤s4：重复步骤s3直到所有分段血管三维重建完成，再将重建后的分段血管合并，即得到完整的三维血管，参照图4至图8为多角度血管合成示意图。

参照图9，该图为血管分叉点对应示意图，a1b1是平面a上的一条血管中心线，其中c1是该血管上的一个分叉点；a2b2是平面b上的一条血管中心线，其中c2是该血管上的一个分叉点；经过人工确认后，在这两条中心线上，点与点之间是相互对应的，a1对应a2，b1对应b2，c1对应c2，此时，分叉点将血管分成两段，a1c1对应a2c2，c1b1对应c2b2；同上所述在每一小段之间进行线性插值得到n个点，并将平面a和平面b上的分段血管一一对应。

参照图10，该图为血管狭窄点对应示意图，a1b1是平面a上的一条血管中心线，其中d1是该血管上的一个狭窄点；a2b2是平面b上的一条血管中心线，其中d2是该血管上的一个狭窄点；经过人工确认后，在这两条中心线上，点与点之间是相互对应的，a1对应a2，b1对应b2，d1对应d2，此时，狭窄点将血管分成两段，a1d1对应a2d2，d1b1对应d2b2；同上所述在每一小段之间进行线性插值得到n个点，并将平面a和平面b上的分段血管一一对应。

通常在一条血管上会同时出现n个(n>＝1)分叉点和狭窄点，所以会将该血管分成n+1段，然后将两x射线冠脉造影图像上的每一段血管一一对应。

上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。