一种结合曲率特征的用递归图重建分叉血管表面的方法技术

本发明专利技术公开一种结合曲率特征的用递归图重建分叉血管表面的方法，包括获取单分支血管骨架线的代数表达，结合曲线的曲率特征对骨架线进行采样；基于曲线的采样点的法向量取横截面，拟合截面上的轮廓点得到轮廓线；采用递归图的方法插补轮廓线，得到每一单分支上完整的轮廓线；将不同分支在分叉区域进行拼接，得到整个血管的完整轮廓线；将轮廓线的对应点用三角面拼接构建血管表面；选择Loop细分方法对血管表面进一步细分，使重构表面更加光滑，该方法结合曲率特征进行采样，降低了采样密度，用递归图方法进行轮廓插补，重建的血管模型更逼真，减小了计算量，提高重建效率，且重建的血管模型具备可重构性，为引入物理特性建立功能性血管模型建立基础。

Reconstruction of bifurcated vascular surface by recursion graph combined with curvature features

The invention discloses a method for reconstructing the surface of a bifurcated vessel by recursive graph in combination with curvature characteristics, which includes obtaining the algebraic expression of the skeleton line of a single branch vessel, sampling the skeleton line in combination with the curvature characteristics of the curve, taking the cross section based on the normal vector of the sampling point of the curve, fitting the outline points on the cross section to obtain the outline line; The complete contour line on each single branch is obtained by interpolating contour line with recursive graph method; the complete contour line of the whole blood vessel is obtained by splicing different branches in the bifurcated region; the corresponding points of the contour line are spliced into triangles to construct the blood vessel surface; Loop subdivision method is selected to further subdivide the blood vessel surface to make the reconstructed table. The surface is smoother. This method combines the curvature feature to sample, reduces the sampling density, uses the recursive graph method to interpolate the contour. The reconstructed vascular model is more realistic, reduces the amount of calculation, improves the reconstructed efficiency, and the reconstructed vascular model is reconfigurable. A functional vascular model is established for introducing physical characteristics. Basics.

【技术实现步骤摘要】
一种结合曲率特征的用递归图重建分叉血管表面的方法
本专利技术涉及生物医学工程和计算机视觉
，具体是一种结合曲率特征的用递归图重建分叉血管表面的方法。
技术介绍
血管结构的三维建模对于血管疾病的临床诊断、虚拟手术、预测病情发展、提高术前规划的可靠性及手术的成功率、临床解剖教学、病理研究等方面具有重要意义和广泛的应用价值。根据三维建模过程中是否基于模型假设，可以将血管造型技术分为基于无模型的方法和基于模型的方法两种。基于无模型方法是通过血管边界的体素直接构造表面，常见的方法有移动立方体(MarchingCubes)和MPU等方法。基于模型的方法一般需要先得到血管的血管中心线和管径大小等信息，再根据这些数据来构建血管表面的多边形网格，主要方法有柱面拟合、圆台面拟合、参数曲面法、细分曲面法、基于截面轮廓构造法等。基于无模型方法简单、效率较高，然而重建产生的三角面片是散乱无规则的，网格呈锯齿状，曲面具有台阶效应，存在着构造等值面的二义性，所以构造的三维模型只能用于血管形态的观测，不能反映血管的管径大小、曲率等详细信息。基于模型方法相对于无模型方法，构建过程会显得较复杂，但是这类方法构建的模型具有可重构性，拥有较多的信息量并具有较逼真的可视化效果，能够满足血管疾病的诊断、介入手术模拟和医学解剖及训练等后续处理应用。基于模型方法也有其自身的缺点，一方面是对于血管表面的重建，并不能精准地反映血管的曲率变化；另一方面是对于分叉区域的拼接，若先构造单分支管道表面再进行拼接，处理起来比较复杂，同时计算量比较大。
技术实现思路
本专利技术的目的在于克服现有技术的不足，而提供一种结合曲率特征的用递归图重建分叉血管表面的方法，该方法不仅能够精准的反映血管的曲率变化，更详细地展示血管结构信息，而且还可以降低采样密度，减小计算量，采用递归图方法可以快速实现单分支血管轮廓的插补，同时在分叉区域，针对轮廓点和交叉轮廓的处理，仍然可以减小计算量，提高模型的重建效率。实现本专利技术目的的技术方案是：一种结合曲率特征的用递归图重建分叉血管表面的方法，具体包括如下步骤：S1、获取单分支血管骨架线的代数表达，结合曲线的曲率特征对骨架线进行采样；S2、基于曲线的采样点的法向量取横截面，拟合截面上的轮廓点得到轮廓线；S3、采用递归图的方法插补步骤S2得到的轮廓线，得到每一单分支上完整的轮廓线；S4、将不同分支在分叉区域进行拼接，得到整个血管的完整轮廓线；S5、将步骤S4得到的轮廓线的对应点用三角面拼接构建血管表面；S6、选择Loop细分方法对血管表面进一步细分，从而使重构表面更加光滑。步骤S1中，所述的获取单分支血管骨架线的代数表达，结合曲线的曲率特征对骨架线进行采样，具体是：S1-1、提取血管骨架线；S1-2、按照血管的管径大小将骨架线重组为单分支；S1-3、用B样条曲线作为基元实现对每一单分支骨架点的拟合，获取单分支血管骨架线的代数表达；S1-4、结合曲线的曲率特征对曲线进行采样，采样点更能反应曲线的特征。步骤S3中，所述的采用递归图的方法插补步骤S2得到的轮廓线，得到每一单分支上完整的轮廓线；具体是：S3-1、在单分支骨架线上，选取采样点对应的横截面上的轮廓线；S3-2、采用递归图的方法，快速对轮廓线进行插补。步骤S4中，所述的将步骤S4得到的轮廓线的对应点用三角面拼接构建血管表面，具体是：S4-1、将不同分支的轮廓树在分叉区域的分叉点处简单组合；S4-2、根据轮廓线上轮廓点与不同骨架线的欧式距离；S4-3、依据不同分支的管径大小，判断查找同时存在于不同分支的轮廓点；S4-4、将同时存在于不同分支的轮廓点删除；S4-5、对于有轮廓交叉的部分，对相应轮廓进行融合；S4-6、将各分叉区域进行拼接，得到完整血管轮廓树。步骤S5中，所述的将轮廓线的对应点用三角面拼接构建血管表面，具体是：S5-1、将轮廓线按照等角度方式取轮廓点；S5-2、将不同轮廓线上的轮廓点对应起来；S5-3、依据对应轮廓点构造三角面，拼接构建血管表面。有益效果：本专利技术提供的一种结合曲率特征的用递归图重建分叉血管表面的方法，结合曲线曲率特征进行采样，一方面降低了采样密度，减小重建轮廓线的计算，另一方面更能反应血管骨架线的特征，使构造出来的模型显现更详细的血管结构信息。用递归图方法实现血管轮廓的快速重建，选取低密度采样点的轮廓线作为控制轮廓线，基于递归方法可以快速的插补轮廓线，这样的轮廓线插补不会丢失太多血管的信息，同时又可以提高效率。构建的血管模型为引进物理特性的分析的后续工作，如血管动态形变、血流分析、血管表面压力分析等，提供了良好的几何基础。附图说明图1为本专利技术实施例中的一种结合曲率特征的用递归图重建分叉血管表面的方法的流程示意图；图2为本专利技术实施例中的基于骨架线曲率特征进行采样的示意图；图3为本专利技术实施例中的基于采样点构造轮廓线的示意图；图4为本专利技术实施例中的采用递归图方法对单分支轮廓线插补的示意图；图5为本专利技术实施例中的采用递归图方法对真实单分支轮廓线进行插补后得到的单分支完整轮廓线；图6为本专利技术实施例中的分叉区域删除轮廓点进行分支拼接的示意图；图7为本专利技术实施例中的重建的模型。具体实施方式下面结合附图和实施例对本专利技术做进一步阐述，但不是对本专利技术的限定。实施例：如图1所示，一种结合曲率特征的用递归图重建分叉血管表面的方法，具体包括如下步骤：S1、获取单分支血管骨架线的代数表达，结合曲线的曲率特征对骨架线进行采样；S2、基于曲线的采样点的法向量取横截面，拟合截面上的轮廓点得到轮廓线；S3、采用递归图的方法插补步骤S2得到的轮廓线，得到每一单分支上完整的轮廓线；S4、将不同分支在分叉区域进行拼接，得到整个血管的完整轮廓线；S5、将步骤S4得到的轮廓线的对应点用三角面拼接构建血管表面；S6、选择Loop细分方法对血管表面进一步细分，从而使重构表面更加光滑。步骤S1中，所述的获取单分支血管骨架线的代数表达，结合曲线的曲率特征对骨架线进行采样，具体是：S1-1、基于拉普拉斯算子方法提取血管点云数据的骨架线；S1-2、按照血管的管径大小将骨架线重组为单分支，并对单分支骨架线进行编号，以使在分叉区域拼接时进行区分；S1-3、用B样条曲线作为基元实现对每一单分支骨架点的拟合，选择逼近方法，从单分支血管段选取n+1个数据点作为型值点，反求该曲线的控制点，采用3次的参数样条曲线，确定节点向量，根据以上四个部分拟合三次B样条曲线，获取单分支血管骨架线的代数表达S1-4、结合曲线的曲率特征对曲线进行采样，采样点更能反应曲线的特征,求得空间曲线的曲率构建采样点M与曲率k大小成正比，为自设定的数目，取3，5，10等，根据实际情况设定，那么就可以依据曲率来进行采样，在曲率大的曲线段采样点适当增多，在曲率小的曲线段采样点要减少，这样不仅可以保持曲线段的形态特征，同时可以降低采样密度，减小后续的计算量。步骤S3中，所述的采用递归图的方法插补步骤S2得到的轮廓线，得到每一单分支上完整的轮廓线，具体是：S3-1、在单分支骨架线上，选取采样点对应的横截面上的轮廓线；S3-2、以选定的轮廓线为控制轮廓，采用递归图的方法，对轮廓线进行快速插补，得到骨架线上较完整的轮廓线。步骤S4中，所述的将不同分支在分叉区域进行拼本文档来自技高网...

【技术保护点】
1.一种结合曲率特征的用递归图重建分叉血管表面的方法，其特征在于，具体包括如下步骤：S1、获取单分支血管骨架线的代数表达，结合曲线的曲率特征对骨架线进行采样；S2、基于曲线的采样点的法向量取横截面，拟合截面上的轮廓点得到轮廓线；S3、采用递归图的方法插补步骤S2得到的轮廓线，得到每一单分支上完整的轮廓线；S4、将不同分支在分叉区域进行拼接，得到整个血管的完整轮廓线；S5、将步骤S4得到的轮廓线的对应点用三角面拼接构建血管表面；S6、选择Loop细分方法对血管表面进一步细分，从而使重构表面更加光滑。

【技术特征摘要】
1.一种结合曲率特征的用递归图重建分叉血管表面的方法，其特征在于，具体包括如下步骤：S1、获取单分支血管骨架线的代数表达，结合曲线的曲率特征对骨架线进行采样；S2、基于曲线的采样点的法向量取横截面，拟合截面上的轮廓点得到轮廓线；S3、采用递归图的方法插补步骤S2得到的轮廓线，得到每一单分支上完整的轮廓线；S4、将不同分支在分叉区域进行拼接，得到整个血管的完整轮廓线；S5、将步骤S4得到的轮廓线的对应点用三角面拼接构建血管表面；S6、选择Loop细分方法对血管表面进一步细分，从而使重构表面更加光滑。2.根据权利要求1所述的一种结合曲率特征的用递归图重建分叉血管表面的方法，其特征在于，步骤S1中，所述的获取单分支血管骨架线的代数表达，结合曲线的曲率特征对骨架线进行采样，具体是：S1-1、提取血管骨架线；S1-2、按照血管的管径大小将骨架线重组为单分支；S1-3、用B样条曲线作为基元实现对每一单分支骨架点的拟合，获取单分支血管骨架线的代数表达；S1-4、结合曲线的曲率特征对曲线进行采样，采样点更能反应曲线的特征。3.根据权利要求1所述的一种结合曲率特征的用递归图重建分叉血管表面的方法，其特征在...

【专利技术属性】
技术研发人员：古天龙，韦永成，蓝如师，常亮，姜静月，罗笑南，
申请(专利权)人：桂林电子科技大学，
类型：发明
国别省市：广西,45