一种基于面模型的断面网格生成方法技术

本发明专利技术公开了一种基于面模型的断面网格生成方法，包括：针对每个三角形，按逆时针顺序将每条边与裁切面进行求交；根据顶点及顶点间的邻接关系，将顶点分成一个或多个首尾相连的闭合多边形点集；计算p沿着指向闭合多边形内部偏移得到的新顶点q，由q组成的新顶点集合记为Q；对集合Q，每次取其中内角最小的顶点q，将其前后相邻的两个顶点X和Y组成三角形，并根据分裂系数μ决定是否取q对边的中点W，并将该三角形分裂成XqW和qYW两个三角形；获取内层网格中的所有三角形顶点，并以原始集合Q中的顶点作为约束，调整裁切面的三角形分布和大小；得到器官模型的切面网格。本发明专利技术具有处理速度快、准确性好、处理效率高等特点。

A surface mesh generation method based on surface model

The invention discloses a section mesh generation method based on a surface model, which comprises: for each triangle, each edge is intersected with a cutting surface in counterclockwise order; vertices are divided into one or more closed polygonal point sets connected by one or more ends according to the adjacency relationship between vertices and vertices; and P is calculated to close more points along the direction; The new vertex Q obtained by the internal migration of the edge is marked as Q by the new vertex set composed of Q. For the set Q, the vertex Q with the smallest internal angle is taken each time, and the two adjacent vertices X and Y are formed into a triangle. According to the splitting coefficient, the midpoint W of the q-pair edge is determined, and the triangle is divided into two triangles X Q W and qYW. All the triangle vertices in the inner grid are obtained, and the triangle distribution and size of the cutting surface are adjusted with the constraint of the vertices in the original set Q. The invention has the advantages of fast processing speed, good accuracy and high processing efficiency.

【技术实现步骤摘要】
一种基于面模型的断面网格生成方法
本专利技术涉及计算机医学辅助诊疗系统中的医学图像可视化处理方法，尤其涉及一种基于面模型的断面网格生成方法。
技术介绍
随着计算机技术的不断发展，计算机辅助诊断(CAD)在临床中扮演着越来越重要的角色。医生采集患者CT图像，通过计算机图像处理技术将感兴趣的器官和组织以二值图像的形式分割出来，然后采用图形技术进行三维重建，从而得到器官和组织的三维网格模型。相对于传统的二维图像诊断，借助病变器官和组织的三维模型，医生可以多角度、全方位地观察和了解患者的病变情况，从而更加准确地对患者的病症和病情进行诊断和评估，进而提出更加有针对性和行之有效的治疗方案。在计算机辅助诊断的整个过程中，由于某些组织器官反应在CT和MRI图像上的特性相似，在三维重建中经常会出现不同的组织器官粘连和某些器官非正常“增生”的现象。因此，在三维重建之后，需要对模型进行后续处理，其中，将粘连器官进行切割，对异常“增生”网格进行裁切切除是非常重要的后处理手段。裁切之后的网格模型会在裁切部位出现“孔洞”，需要对“孔洞”进行三角化生成网格进行添补。目前已经存在的网格裁切面生成方法主要有以下几种：第一种，3DMax、Maya等网格处理工具。该类工具软件提供孔洞添补工具，可以自动识别网格存在的孔洞并进行自动添补。由于三维重建之后的三维网格模型比较复杂，裁切之后生成的孔洞边缘信息可能存在冗余和错误，直接采用该类工具进行处理可能会出现错误结果甚至导致程序崩溃，导致该方法无法使用。第二种，基于耳裁切算法的裁切面生成。由于该方法无法再裁切部分(孔洞)内部生成新的顶点，因此采用该方法生成的裁切面网格多为较狭长的三角形。在原网格三角形数量密度较大时，裁切面的三角形密度会过少，可能导致面模型后续处理上的问题。在裁切面与原网格的夹角较大时，连接部分的法向量计算会出现突变，导致渲染效果较差。请参照图5，传统裁切面生成算法生成的网格简单地将裁切生成的顶点相连，裁切面与被裁切网格的交界处由于法线计算的畸变导致渲染效果极差。
技术实现思路
本专利技术要解决的技术问题在于，针对现有技术的不足，提供一种处理速度快、准确性好、处理效率高的基于面模型的断面网格生成方法。为解决上述技术问题，本专利技术采用如下技术方案。一种基于面模型的断面网格生成方法，其包括有：步骤S1，导入3D医学图像，对3D医学图像中需要离断的器官模型进行裁切，得到裁切面；步骤S2，针对每个三角形，按逆时针顺序将每条边与裁切面进行求交，获取裁切面的顶点并生成顶点之间的邻接关系；步骤S3，根据顶点及顶点间的邻接关系，将顶点分成一个或多个首尾相连的闭合多边形点集；步骤S4，将步骤S3获取的每一个闭合多边形点集记为P，取闭合多边形点集中的每一个顶点记为p，取内角平分线的方向向量记为d，计算p沿着指向闭合多边形内部偏移得到的新顶点q，则有q＝p+u*d，其中，u为数值小于0.01的标量，之后由q组成的新顶点集合记为Q，将集合P和Q中的顶点根据位置和拓扑关系进行连接，得到外层三角形带；步骤S5，对步骤S4获得的集合Q，每次取其中内角最小的顶点q，将其前后相邻的两个顶点X和Y组成三角形，并根据分裂系数μ决定是否取q对边的中点W，并将该三角形分裂成XqW和qYW两个三角形，其中，μ＝f*Lavg/Lxy，f为细分参数，Lavg为边的平均长度，Lxy为边XY的长度；步骤S6，获取步骤S5得到的内层网格中的所有三角形顶点，并以原始集合Q中的顶点作为约束，调整裁切面的三角形分布和大小；步骤S7，得到器官模型的切面网格。优选地，所述步骤S5中，若分裂(μ>1)，则用W替代q参加下一轮计算，否则将q从集合Q中删除。优选地，所述步骤S5中，更新该步骤所修改的顶点内角及其左右相邻的两个顶点的邻接关系，之后进行下一轮计算，直至集合Q中的顶点数目为3，将剩余的3个顶点组合成三角形。优选地，所述步骤S6中，利用最小二乘法调整裁切面的三角形分布和大小。本专利技术公开的基于面模型的断面网格生成方法中，其主要针对裁切之后的网格进行裁切面的三角形网格生成。首先获取裁切得到的顶点集合及其相互之间的拓扑顺序关系。然后对顶点集合进行预处理，根据顶点及其邻接关系，将顶点分成一个或多个的闭合多边形顶点集合。对于集合中的每一个闭合多边形点集P，将其中的每一个顶点向孔洞的内部进行长度为d的偏移并生成新的顶点集合Q。将集合P和Q中的顶点根据位置关系生成三角形带。对顶点集Q，每次选取其中内角最小的顶点p，并将p和其前后相邻的两个顶点Y和Z组成三角形，并根据分裂系数μ决定是否取三角形的最长边的中点W并将该三角形分裂成两个三角形YpW和pZW。若分裂，则最长边的中点W代替顶点p参加下一轮计算，否则将该顶点移除集合Q。每一轮计算完成后需要更新本轮修改的顶点的内角及其左右相邻的两个顶点的邻接关系。当顶点集合Q中顶点数量为3时，将剩余的3个顶点组成三角形，计算结束。最后，将上述过程中所有生成的三角形进行最小二乘法优化，并以固定原始集合Q中的所有顶点作为约束条件，求所有裁切面网格上所有顶点的最优位置。相比现有技术而言，本专利技术对于根据CT图像重建后的网格，将“黏连”的器官和“增生”的组织切除是一种常见的操作，切割面的网格生成速度和效果对于后续的网格处理和观察有着重大的意义，可以有效地提高计算机辅助诊断的效率。附图说明图1为本专利技术断面网格生成方法的流程图。图2为优选实施例中裁切过后的胆囊模型图。图3为生成了两层三角形带的结果图。图4为裁切面网格生成结果图及优化结果图。图5为传统裁切面生成算法的运行结果图。具体实施方式下面结合附图和实施例对本专利技术作更加详细的描述。本专利技术公开了一种基于面模型的断面网格生成方法，请参照图1，其包括有：步骤S1，导入3D医学图像，对3D医学图像中需要离断的器官模型进行裁切，得到裁切面；步骤S2，针对每个三角形，按逆时针顺序将每条边与裁切面进行求交，获取裁切面的顶点并生成顶点之间的邻接关系；步骤S3，根据顶点及顶点间的邻接关系，将顶点分成一个或多个首尾相连的闭合多边形点集；步骤S4，将步骤S3获取的每一个闭合多边形点集记为P，取闭合多边形点集中的每一个顶点记为p，取内角平分线的方向向量记为d，计算p沿着指向闭合多边形内部偏移得到的新顶点q，则有q＝p+u*d，其中，u为数值小于0.01的标量，之后由q组成的新顶点集合记为Q，将集合P和Q中的顶点根据位置和拓扑关系进行连接，得到外层三角形带；步骤S5，对步骤S4获得的集合Q，每次取其中内角最小的顶点q，将其前后相邻的两个顶点X和Y组成三角形，并根据分裂系数μ决定是否取q对边的中点W，并将该三角形分裂成XqW和qYW两个三角形，其中，μ＝f*Lavg/Lxy，f为细分参数，Lavg为边的平均长度，Lxy为边XY的长度；步骤S6，获取步骤S5得到的内层网格中的所有三角形顶点，并以原始集合Q中的顶点作为约束，调整裁切面的三角形分布和大小；步骤S7，得到器官模型的切面网格。上述方法中，其主要针对裁切之后的网格进行裁切面的三角形网格生成。首先获取裁切得到的顶点集合及其相互之间的拓扑顺序关系。然后对顶点集合进行预处理，根据顶点及其邻接关系，将顶点分成一个或多个的闭合多边形顶点集合。对于集本文档来自技高网...

【技术保护点】
1.一种基于面模型的断面网格生成方法，其特征在于，包括有：步骤S1，导入3D医学图像，对3D医学图像中需要离断的器官模型进行裁切，得到裁切面；步骤S2，针对每个三角形，按逆时针顺序将每条边与裁切面进行求交，获取裁切面的顶点并生成顶点之间的邻接关系；步骤S3，根据顶点及顶点间的邻接关系，将顶点分成一个或多个首尾相连的闭合多边形点集；步骤S4，将步骤S3获取的每一个闭合多边形点集记为P，取闭合多边形点集中的每一个顶点记为p，取内角平分线的方向向量记为d，计算p沿着指向闭合多边形内部偏移得到的新顶点q，则有q＝p+u*d，其中，u为数值小于0.01的标量，之后由q组成的新顶点集合记为Q，将集合P和Q中的顶点根据位置和拓扑关系进行连接，得到外层三角形带；步骤S5，对步骤S4获得的集合Q，每次取其中内角最小的顶点q，将其前后相邻的两个顶点X和Y组成三角形，并根据分裂系数μ决定是否取q对边的中点W，并将该三角形分裂成XqW和qYW两个三角形，其中，μ＝f*Lavg/Lxy，f为细分参数，Lavg为边的平均长度，Lxy为边XY的长度；步骤S6，获取步骤S5得到的内层网格中的所有三角形顶点，并以原始集合Q中的顶点作为约束，调整裁切面的三角形分布和大小；步骤S7，得到器官模型的切面网格。...

【技术特征摘要】
1.一种基于面模型的断面网格生成方法，其特征在于，包括有：步骤S1，导入3D医学图像，对3D医学图像中需要离断的器官模型进行裁切，得到裁切面；步骤S2，针对每个三角形，按逆时针顺序将每条边与裁切面进行求交，获取裁切面的顶点并生成顶点之间的邻接关系；步骤S3，根据顶点及顶点间的邻接关系，将顶点分成一个或多个首尾相连的闭合多边形点集；步骤S4，将步骤S3获取的每一个闭合多边形点集记为P，取闭合多边形点集中的每一个顶点记为p，取内角平分线的方向向量记为d，计算p沿着指向闭合多边形内部偏移得到的新顶点q，则有q＝p+u*d，其中，u为数值小于0.01的标量，之后由q组成的新顶点集合记为Q，将集合P和Q中的顶点根据位置和拓扑关系进行连接，得到外层三角形带；步骤S5，对步骤S4获得的集合Q，每次取其中内角最小的顶点q，将其前后相邻的两个顶点X和Y组成三角形，并根据分裂系数μ决定是否取q对边的中点W...

【专利技术属性】
技术研发人员：叶建平，范应方，
申请(专利权)人：深圳市一图智能科技有限公司，
类型：发明
国别省市：广东,44