[发明专利]一种共面网格单元自动分类方法、装置和可读存储介质

说明书

本发明提出了一种共面网格单元自动分类方法、装置和可读存储介质，该方法包括如下步骤：提取三角形网格单元集；定义共面网格单元集合和临时网格单元集合，并初始化；从三角形网格单元集中随机选取种子网络单元，并存入共面网格单元集合；从共面网格单元集合中，选取种子网格单元；从三角形网格单元集中将与种子网格单元共节点的单元取出，存入临时网格单元集合并从中随机选取网格单元；分别计算种子网格单元和网格单元的法向量，并计算两者法向量的夹角；依次判断临时网格单元集合、共面网格单元集合和三角形网格单元集，完成对共面网格单元的分类。本发明实现了在平面区域边界处，对边界处网格选取精确控制，避免造成跨区域网格误选，效率高。

技术领域

本发明涉及医学有限元分析领域，具体而言，涉及一种共面网格单元自动分类方法、装置和可读存储介质。

背景技术

随着计算力学和医学图像的快速发展，通过血液动力学仿真模拟血液在血管中的流动情况，能够充分了解心脑血管疾病患者的血管健康状况。为了能够精确模拟血液的流动状态，需要对血管网格单元模型中血流出口端、入口端及模型表面的网格进行精确选择，施加相应的边界条件。目前，网格选取的方式，包括种子网格区域生长法、网格邻近区域扩展方法。但在平面区域边界处，这两种方法对边界处网格的选取难以精确控制，易造成跨区域的网格误选，最终需要人工进行干预，逐个移除误选网格，效率较低。

发明内容

鉴于上述问题，本发明提供了一种共面网格单元自动分类方法、装置和可读存储介质。

为解决上述技术问题，本发明采用的技术方案是：一种共面网格单元自动分类方法，包括如下步骤：S101，获取网格数据，从中提取出三角形网格单元集；

S102，定义共面网格单元集合和临时网格单元集合，并将其初始化，所述定义共面网格单元集合和临时网格单元集合均包含于三角形网格单元集；S103，从所述三角形网格单元集中随机选取1个网格单元作为种子网络单元，并将其存入共面网格单元集合；S104，从所述共面网格单元集合未遍历过的网格单元中，选取1个网格单元作为种子网格单元；S105，从所述三角形网格单元集中将与种子网格单元共节点的单元取出，并存入临时网格单元集合；S106，从临时网格单元集合中随机选取1个网格单元；S107，分别计算所述种子网格单元和网格单元的法向量，并计算两者法向量的夹角，若，则将所述网格单元存入共面网格单元集合,否则，将其放回三角形网格单元集，其中为指定阈值，用于判定2个网格单元是否共面；S108，依次判断所述临时网格单元集合是否为空，共面网格单元集合是否遍历完成，三角形网格单元集是否为空，根据判断结果重复上述步骤，完成对共面网格单元的分类。

作为优选方案，将三角形网格单元集记为E，，三角形网格节点集记为P，，三角形网格节点对应的坐标为，每个单元对应的节点关系为；则每个单元对应的空间法向量的计算公式如下：

。

作为优选方案，在步骤S102中，所述初始化为将共面网格单元集合和临时网格单元集合置为空集。

作为优选方案，在步骤S107中，所述两者法向量的夹角的计算公式如下所示：

作为优选方案，所述步骤S108具体包括：S1081，判断所述临时网格单元集合是否为空，若为空，则重复步骤S106-S108，反之，执行步骤S1082；S1082，判断所述共面网格单元集合中的网格单元是否全部遍历完成，若未完成，则重复步骤S104-S1081，反之，执行步骤S1083；S1083，判断三角形网格单元集是否为空，若不为空，则重复步骤S103- S1082，反之，则完成对共面网格单元的分类。

作为优选方案，所述判断所述共面网格单元集合中的网格单元是否全部遍历完成，具体为：判断共面网格单元集合中是否有网格单元未在步骤S104中作为种子单元被选取过。