CN116862816A - 一种基于三维形态学的模型膨胀方法、装置及应用 - Google Patents

一种基于三维形态学的模型膨胀方法、装置及应用 Download PDF

Info

Publication number
CN116862816A
CN116862816A CN202310945642.8A CN202310945642A CN116862816A CN 116862816 A CN116862816 A CN 116862816A CN 202310945642 A CN202310945642 A CN 202310945642A CN 116862816 A CN116862816 A CN 116862816A
Authority
CN
China
Prior art keywords
expansion
voxel
dimensional array
dimensional
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310945642.8A
Other languages
English (en)
Other versions
CN116862816B (zh
Inventor
蔡林凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Yazhi Medical Technology Co ltd
Original Assignee
Hangzhou Yazhi Medical Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Yazhi Medical Technology Co ltd filed Critical Hangzhou Yazhi Medical Technology Co ltd
Priority to CN202310945642.8A priority Critical patent/CN116862816B/zh
Publication of CN116862816A publication Critical patent/CN116862816A/zh
Application granted granted Critical
Publication of CN116862816B publication Critical patent/CN116862816B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • G06T5/30Erosion or dilatation, e.g. thinning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20036Morphological image processing

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Image Generation (AREA)

Abstract

本申请提出了一种基于三维形态学的模型膨胀方法、装置及应用,包括以下步骤:构建包围盒围绕待膨胀模型,对其进行采样得到固定大小的体素集合;将体素集合保存到原始三维数组,并计算每个体素与其三角面片邻近体素中心点的距离,若最短邻近距离小于设定距离,则认为邻近体素包含三角面片将其标记为1,否则为0;使用立体结构元遍历膨胀三维数组,计算膨胀系数,若大于膨胀阈值则停止膨胀;基于膨胀三维数组进行模型重建,得到待膨胀模型的膨胀结果。本方案可以将模型体素化后进行膨胀来改变模型的形状解决膨胀过程中的自相交问题。

Description

一种基于三维形态学的模型膨胀方法、装置及应用
技术领域
本申请涉及三维模型领域,特别是涉及一种基于三维形态学的模型膨胀方法、装置及应用,能有效避免模型在膨胀过程中产生的自相交情况。
背景技术
膨胀是形态学中的一种基本操作,用于扩大图像中物体的边界或填充空洞,它通过对图像中的每个像素应用立体结构元素,根据立体结构元素的形状和大小将像素集合向外部蔓延,常规的二维膨胀操作通过扫描原图,找到像素值为0的点并将选择好的结构原点移动到该点,判断该立体结构元素所覆盖范围内的像素值是否至少存在一个像素值为1的点,如果存在,该点的像素值膨胀为1,如果立体结构元素所覆盖范围内像素值全为0,则该点不进行膨胀处理保持不变,遍历原图中其他的像素值为0的点,重复以上操作完成膨胀。
在二维场景中膨胀有着广泛的应用,是形态学图像处理中不可缺少的基本操作,可以通过膨胀处理来填充二维图像中的空洞部分,将空洞部分中的像素点与周围的像素点重新连接起来以完善图像的形状,膨胀操作可以使物体的边界向外蔓延,从而扩大物体的区域,这在图像分割、目标检测和形态学重建等任务中经常被用来扩大和增强物体的边界,还可以通过膨胀操作来去除图像中的小尺寸噪声,即通过立体结构元素对图像中的噪声进行膨胀,使其连接成较大的连通域然后再进行腐蚀操作以去除噪声。
但是,在膨胀过程中会存在模型的不同部分相互重叠或交叉的自相交情况,这可能导致模型的几何形状不准确造成模型的失真或无效的结果,自相交问题通常在模型形状复杂或存在尖端、凸起的情况下更加容易发生。
发明内容
本申请实施例提供了一种基于三维形态学的模型膨胀方法、装置及应用,可以通过先将模型体素化再膨胀的方式来避免模型膨胀过程中出现的自相交问题。
第一方面,本申请实施例提供了一种基于三维形态学的模型膨胀方法,所述方法包括:
获取待膨胀模型,构建包围盒将所述待膨胀模型包围,并以固定大小的体素对所述包围盒进行采样得到体素集合;
将体素集合中的每一体素根据三维坐标保存到原始三维数组中,识别所述待膨胀模型内的每一三角面片,并将所述原始三维数组中所有包含三角面片的体素值记为1,不包含三角面片的体素值记为0得到模型体素三维数组;
设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中,定义由多个体素构建的立体结构元,使用所述立体结构元的中心体素依次遍历所述膨胀三维数组中的每一体素,遍历过程中若立体结构元所覆盖范围内的体素全部为0,则被遍历体素保持不变,若立体结构元所覆盖范围内至少存在一体素值为1的体素,则将被遍历体素的体素值设定为1进行膨胀,所述被遍历体素为与立体结构元中心体素重合的体素;
实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀,并基于膨胀三维数组进行模型重建得到待膨胀模型的膨胀结果。
第二方面,本申请实施例提供了一种基于三维形态学的模型膨胀装置,包括:
获取模块:获取待膨胀模型,构建包围盒将所述待膨胀模型包围,并以固定大小的体素对所述包围盒进行采样得到体素集合;
定义模块:将体素集合中的每一体素根据三维坐标保存到原始三维数组中,识别所述待膨胀模型内的每一三角面片,并将所述原始三维数组中所有包含三角面片的体素值记为1,不包含三角面片的体素值记为0得到模型体素三维数组;
膨胀模块:设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中,定义由多个体素构建的立体结构元,使用所述立体结构元的中心体素依次遍历所述膨胀三维数组中的每一体素,遍历过程中若立体结构元所覆盖范围内的体素全部为0,则被遍历体素保持不变,若立体结构元所覆盖范围内至少存在一体素值为1的体素,则将被遍历体素的体素值设定为1进行膨胀,所述被遍历体素为与立体结构元中心体素重合的体素;
重建模块:实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀,并基于膨胀三维数组进行模型重建得到待膨胀模型的膨胀结果。
第三方面,本申请实施例提供了一种电子装置,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器被设置为运行所述计算机程序以执行一种基于三维形态学的模型膨胀方法。
第四方面,本申请实施例提供了一种可读存储介质,所述可读存储介质中存储有计算机程序,所述计算机程序包括用于控制过程以执行过程的程序代码,所述过程包括根据一种基于三维形态学的模型膨胀方法。
本发明的主要贡献和创新点如下:
本申请实施例通过构建包围盒来将待膨胀模型包围,并定义固定大小的体素来进行采样从而使得待膨胀模型体素化,将待膨胀模型体素化后可以准确地表示模型的形状,并且不存在曲面或边缘相交的情况;
本方案通过对体素进行膨胀操作,并通过增加邻近体素来改变模型的形状从而解决自相交问题,膨胀操作可以增加体素之间的间隙,使得模型的曲面或边缘在离散空间中保持分离,并且不再相交;
本方案采用轴包围盒进行体素化,轴包围盒体素化可以利用轴对齐的包围盒来表示整个场景或模型的范围。这样,只需要存储包围盒的最小和最大坐标,不需要存储每个具体的体素状态,因此内存占用较小;轴包围盒体素化仅需要对整个包围盒进行均匀划分,而不需要进行复杂的递归判断和分割,因此构建时间较短;由于轴包围盒体素化是根据整个包围盒进行简单均匀划分的,查询特定体素所处的位置只需要根据体素的空间坐标进行简单计算,因此查询时间较短。
本申请的一个或多个实施例的细节在以下附图和描述中提出,以使本申请的其他特征、目的和优点更加简明易懂。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是根据本申请实施例的一种基于三维形态学的模型膨胀方法的流程图;
图2是根据本申请实施例的一种包围盒包围待膨胀模型的示意图;
图3是根据本申请实施例的原始三维数组的示意图;
图4是根据本申请实施例的膨胀三维数组的示意图;
图5是根据本申请实施例的立体结构元示意图;
图6是根据本申请实施例的一种基于三维形态学的模型膨胀装置的结构框图;
图7是根据本申请实施例的电子装置的硬件结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本说明书一个或多个实施例相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本说明书一个或多个实施例的一些方面相一致的装置和方法的例子。
需要说明的是:在其他实施例中并不一定按照本说明书示出和描述的顺序来执行相应方法的步骤。在一些其他实施例中,其方法所包括的步骤可以比本说明书所描述的更多或更少。此外,本说明书中所描述的单个步骤,在其他实施例中可能被分解为多个步骤进行描述;而本说明书中所描述的多个步骤,在其他实施例中也可能被合并为单个步骤进行描述。
实施例一
本申请实施例提供了一种基于三维形态学的模型膨胀方法,具体地,参考图1,所述方法包括:
获取待膨胀模型,构建包围盒将所述待膨胀模型包围,并以固定大小的体素对所述包围盒进行采样得到体素集合;
将体素集合中的每一体素根据三维坐标保存到原始三维数组中,识别所述待膨胀模型内的每一三角面片,并将所述原始三维数组中所有包含三角面片的体素值记为1,不包含三角面片的体素值记为0得到模型体素三维数组;
设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中,定义由多个体素构建的立体结构元,使用所述立体结构元的中心体素依次遍历所述膨胀三维数组中的每一体素,遍历过程中若立体结构元所覆盖范围内的体素全部为0,则被遍历体素保持不变,若立体结构元所覆盖范围内至少存在一体素值为1的体素,则将被遍历体素的体素值设定为1进行膨胀,所述被遍历体素为与立体结构元中心体素重合的体素;
实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀,并基于膨胀三维数组进行模型重建得到待膨胀模型的膨胀结果。
在本方案中,在“构建包围盒将所述待膨胀模型包围”步骤中,定义一个三维坐标系,使用所述待膨胀模型在三维坐标系中每一维度坐标的最大值和最小值来构建包围盒。
具体的,如图2所示,本方案使用AABB包围盒(Axis-aligned bounding box)将所述待膨胀模型包围起来,通过对所述待膨胀模型中的所有顶点进行扫描得到三维坐标系中每一维度坐标的最大值与最小值分别为xmin、xmax、ymin、ymax、zmin、zmax,通过这些最大值与最小值构建包围盒可以精准的将所述带膨胀包围盒进行包围。
具体的,所述包围盒为轴对齐包围盒,在本方案中的包围盒为长方体包围盒,使用轴对齐包围盒来包围待膨胀模型的好处为:
1.内存占用较小:轴包围盒体素化可以利用轴对齐的包围盒来表示整个场景或模型的范围。这样,只需要存储包围盒的最小和最大坐标,不需要存储每个具体的体素状态,因此内存占用较小;
2.构建时间较短:轴包围盒体素化仅需要对整个包围盒进行均匀划分,而不需要进行复杂的递归判断和分割,因此构建时间较短。
3.查询时间较短:由于轴包围盒体素化是根据整个包围盒进行简单均匀划分的,查询特定体素所处的位置只需要根据体素的空间坐标进行简单计算,因此查询时间较短。
具体的,假设体素的固定大小为N*N*N,则使用固定大小的体素在包围盒内进行均匀的采样得到体素集合,所述体素集合是用离散的网格来表示包围盒内的每一部分,这样做的好处是方便后续的处理。
在“将体素集合中的每一体素根据三维坐标保存到原始三维数组中”步骤中,以所述待膨胀模型在三维坐标系中其中一维的最大值与最小值的差值和体素大小来决定原始三维数组的大小。
具体的,所述原始三维数组如图3所示,记所述待膨胀模型在三维坐标系中x轴的最大值为xmax,最小值为xmin,体素大小为N*N*N,则原始三维数组的大小为:
((xmax-xmin)/N)×((xmax-xmin)/N)×((xmax-xmmin)/N)定义原始三维数组来保存包围盒中的所有体素可以方便随机读取与后续的计算。
在“识别所述待膨胀模型内的每一三角面片”步骤中,获取所述待膨胀模型中每一三角面片在原始三维数组中的位置,获取每一三角面片的至少一邻近体素,计算每一三角面片与对应的每一邻近体素中心点的邻近距离,若最短的邻近距离小于设定距离,则认为最短的邻近距离对应的邻近体素包含三角面片。
在“计算每一三角面片与对应的每一邻近体素中心点的邻近距离”步骤中,获取三角面片的三个顶点坐标,使用三角面片的三个顶点坐标计算三角面片的法向量,获取邻近体素的中心点坐标,计算三角面片的法向量与对应邻近体素中心点坐标的内积得到中心点坐标到三角面片的有向距离,使用中心点坐标减去有向距离与三角面片法向量的乘积得到垂足点坐标,若所述垂足点在三角面片内则邻近体素与三角面片的邻近距离为邻近体素中心点到垂足点的距离,若所述垂足点不在三角面片内,则获取垂足点与三角面片的最短直线,使用最短直线与有向距离构成一个直角三角形,直角三角形的斜边为邻近体素与三角面片的邻近距离。
在一些具体实施例中,要想计算每一三角面片与对应的每一邻近体素中心点的邻近距离,要先通过以下步骤计算邻近体素中心点在三角面片所在平面的垂足点p1:
1、确定三角面片的三个顶点坐标,假设分别为A(x1,y1,z1),B(x2,y2,z2),C(x3,y3,z3)。
2、计算三角面片的法向量,该法向量可以通过以下公式计算:N=(AB×AC),其中N为法向量,×表示向量的叉积运算,AB表示从点A指向点B的向量,AC表示从点A指向点C的向量。
3、将邻近体素中心点P(x,y,z)与三角面片法向量N进行内积运算得到点P到三角面片的有向距离d=N·P。
4、根据Q=P-d*N计算垂足点的坐标,P为邻近体素中心点,d为点P到三角面片的有向距离,N为三角面片的法向量。
根据垂足点的坐标来计算得到每一邻近体素中心点与对应三角面片的邻近距离。
具体的,通过垂足点来判断哪个体素是包含三角面片的体素可以在原始三维数组集合中精确的找到包含三角面片的体素,方面后续进行膨胀。
在“设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中”步骤中,膨胀三维数组的边界相较于原始三维数组的边界向外膨胀,把模型体素三维数组中的体素复制到膨胀三维数组后将膨胀三维数组中剩余体素位置的体素值设置为0。
在本方案中,膨胀三维数组的边界相较于原始三维数组的边界向外膨胀的膨胀程度正相关于立体结构元的立体尺寸。
具体的,膨胀三维数组的如图4所示,由于本方案的目的是对待膨胀模型进行膨胀,所以在膨胀之前要对原始三维数组进行扩大。
具体的,在本方案中所述膨胀三维数组的长宽高相交与原始三维数组均扩大了两个单位大小。
本方案所使用的立体结构元如图5所示,立体结构元为3×3×3的立方体立体结构元,所述立体结构元的8个顶点所在体素为0,其他19个体素都是1,也就是图5中的黑色部分为0,白色部分为1,这样设置的好处在于所述立体结构元适用于不同方向和形状的待膨胀模型,而且在进行膨胀操作时,对角线方向上的膨胀会与水平和垂直方向上的膨胀保持一致,保持待膨胀模型的整体对称性。
具体的,使用立体结构元的中心体素依次遍历所述膨胀三维数组进行膨胀,其膨胀过程用以下公式进行表示:
其中,B为立体结构元,A为待膨胀模型,该公式表示将立体结构元B的中心体素平移到待膨胀模型的(x,y,z)位置。
示例性的,立体结构元是一个立方体,立体结构元的中心体素遍历每一个体素,立体结构元会覆盖以该体素为中心的27个体素,如果立体结构元为1的19个网格所对应的体素中体素值至少有一个为1,则被遍历的那个体素的体素值膨胀为1,否则保持不变。
在“实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀”步骤中,将所述膨胀后膨胀三维数组与膨胀前膨胀三维数组的比值乘以最大膨胀倍数的倒数作为膨胀系数。
示例性的,本方案中的最大膨胀倍数设置为3,在膨胀系数的计算公式为:
其中,最大膨胀倍数的倒数为1/3,Q为膨胀后膨胀三维数组,P为膨胀前三维数组。
具体的,本方案通过定义膨胀阈值的方式来控制膨胀的大小。
在一些具体实施例中,本方案使用Meshing Cubes算法基于膨胀三维数组进行模型重建的步骤如下:
1、设置等值面数值,将膨胀三维数组中每一体素的顶点数值与等值面数值进行比较判断得到顶点状态表,所述顶点状态表表示每一体素的顶点在等值面内或等值面外;
2、通过顶点状态表找出与等值面相交的体素边,并通过线性插值发来计算相交点的坐标;
3、通过中心差分法计算体素八个顶点的法向量,再通过线性插值得到每个三角面片各个顶点的法向量;
4、使用每一体素的顶点坐标和法向量连接三角面片得到等值面的三角形网格。
具体的,生成的三角形网格可能包含不合理的拓扑结构、表面缺陷和过度细分等问题,所以需要对其进行优化处理,优化处理的过程包括拓扑优化、表面光滑和去除不需要的三角形等操作,以提高网格的质量和可视化效果。
具体的,将优化后的网格进行渲染和可视化来完成模型的重建,在本方案中可以使用3D渲染引擎或者其他可视化工具来展示重建的模型。
实施例二
基于相同的构思,参考图6,本申请还提出了一种基于三维形态学的模型膨胀装置,包括:
获取模块:获取待膨胀模型,构建包围盒将所述待膨胀模型包围,并以固定大小的体素对所述包围盒进行采样得到体素集合;
定义模块:将体素集合中的每一体素根据三维坐标保存到原始三维数组中,识别所述待膨胀模型内的每一三角面片,并将所述原始三维数组中所有包含三角面片的体素值记为1,不包含三角面片的体素值记为0得到模型体素三维数组;
膨胀模块:设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中,定义由多个体素构建的立体结构元,使用所述立体结构元的中心体素依次遍历所述膨胀三维数组中的每一体素,遍历过程中若立体结构元所覆盖范围内的体素全部为0,则被遍历体素保持不变,若立体结构元所覆盖范围内至少存在一体素值为1的体素,则将被遍历体素的体素值设定为1进行膨胀,所述被遍历体素为与立体结构元中心体素重合的体素;
重建模块:实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀,并基于膨胀三维数组进行模型重建得到待膨胀模型的膨胀结果。
实施例三
本实施例还提供了一种电子装置,参考图7,包括存储器404和处理器402,该存储器404中存储有计算机程序,该处理器402被设置为运行计算机程序以执行上述任一项方法实施例中的步骤。
具体地,上述处理器402可以包括中央处理器(CPU),或者特定集成电路(ApplicationSpecificIntegratedCircuit,简称为ASIC),或者可以被配置成实施本申请实施例的一个或多个集成电路。
其中,存储器404可以包括用于数据或指令的大容量存储器404。举例来说而非限制,存储器404可包括硬盘驱动器(HardDiskDrive,简称为HDD)、软盘驱动器、固态驱动器(SolidStateDrive,简称为SSD)、闪存、光盘、磁光盘、磁带或通用串行总线(UniversalSerialBus,简称为USB)驱动器或者两个或更多个以上这些的组合。在合适的情况下,存储器404可包括可移除或不可移除(或固定)的介质。在合适的情况下,存储器404可在数据处理装置的内部或外部。在特定实施例中,存储器404是非易失性(Non-Volatile)存储器。在特定实施例中,存储器404包括只读存储器(Read-OnlyMemory,简称为ROM)和随机存取存储器(RandomAccessMemory,简称为RAM)。在合适的情况下,该ROM可以是掩模编程的ROM、可编程ROM(ProgrammableRead-OnlyMemory,简称为PROM)、可擦除PROM(ErasableProgrammableRead-OnlyMemory,简称为EPROM)、电可擦除PROM(ElectricallyErasableProgrammableRead-OnlyMemory,简称为EEPROM)、电可改写ROM(ElectricallyAlterableRead-OnlyMemory,简称为EAROM)或闪存(FLASH)或者两个或更多个以上这些的组合。在合适的情况下,该RAM可以是静态随机存取存储器(StaticRandom-AccessMemory,简称为SRAM)或动态随机存取存储器(DynamicRandomAccessMemory,简称为DRAM),其中,DRAM可以是快速页模式动态随机存取存储器404(FastPageModeDynamicRandomAccessMemory,简称为FPMDRAM)、扩展数据输出动态随机存取存储器(ExtendedDateOutDynamicRandomAccessMemory,简称为EDODRAM)、同步动态随机存取内存(SynchronousDynamicRandom-AccessMemory,简称SDRAM)等。
存储器404可以用来存储或者缓存需要处理和/或通信使用的各种数据文件,以及处理器402所执行的可能的计算机程序指令。
处理器402通过读取并执行存储器404中存储的计算机程序指令,以实现上述实施例中的任意一种基于三维形态学的模型膨胀方法。
可选地,上述电子装置还可以包括传输设备406以及输入输出设备408,其中,该传输设备406和上述处理器402连接,该输入输出设备408和上述处理器402连接。
传输设备406可以用来经由一个网络接收或者发送数据。上述的网络具体实例可包括电子装置的通信供应商提供的有线或无线网络。在一个实例中,传输设备包括一个网络适配器(Network Interface Controller,简称为NIC),其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中,传输设备406可以为射频(Radio Frequency,简称为RF)模块,其用于通过无线方式与互联网进行通讯。
输入输出设备408用于输入或输出信息。在本实施例中,输入的信息可以是待膨胀模型、立体结构元等,输出的信息可以是膨胀结果等。
可选地,在本实施例中,上述处理器402可以被设置为通过计算机程序执行以下步骤:
S101、获取待膨胀模型,构建包围盒将所述待膨胀模型包围,并以固定大小的体素对所述包围盒进行采样得到体素集合;
S102、将体素集合中的每一体素根据三维坐标保存到原始三维数组中,识别所述待膨胀模型内的每一三角面片,并将所述原始三维数组中所有包含三角面片的体素值记为1,不包含三角面片的体素值记为0得到模型体素三维数组;
S103、设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中,定义由多个体素构建的立体结构元,使用所述立体结构元的中心体素依次遍历所述膨胀三维数组中的每一体素,遍历过程中若立体结构元所覆盖范围内的体素全部为0,则被遍历体素保持不变,若立体结构元所覆盖范围内至少存在一体素值为1的体素,则将被遍历体素的体素值设定为1进行膨胀,所述被遍历体素为与立体结构元中心体素重合的体素;
S104、实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀,并基于膨胀三维数组进行模型重建得到待膨胀模型的膨胀结果。
需要说明的是,本实施例中的具体示例可以参考上述实施例及可选实施方式中所描述的示例,本实施例在此不再赘述。
通常,各种实施例可以以硬件或专用电路、软件、逻辑或其任何组合来实现。本发明的一些方面可以以硬件来实现,而其他方面可以以可以由控制器、微处理器或其他计算设备执行的固件或软件来实现,但是本发明不限于此。尽管本发明的各个方面可以被示出和描述为框图、流程图或使用一些其他图形表示,但是应当理解,作为非限制性示例,本文中描述的这些框、装置、系统、技术或方法可以以硬件、软件、固件、专用电路或逻辑、通用硬件或控制器或其他计算设备或其某种组合来实现。
本发明的实施例可以由计算机软件来实现,该计算机软件由移动设备的数据处理器诸如在处理器实体中可执行,或者由硬件来实现,或者由软件和硬件的组合来实现。包括软件例程、小程序和/或宏的计算机软件或程序(也称为程序产品)可以存储在任何装置可读数据存储介质中,并且它们包括用于执行特定任务的程序指令。计算机程序产品可以包括当程序运行时被配置为执行实施例的一个或多个计算机可执行组件。一个或多个计算机可执行组件可以是至少一个软件代码或其一部分。另外,在这一点上,应当注意,如图7中的逻辑流程的任何框可以表示程序步骤、或者互连的逻辑电路、框和功能、或者程序步骤和逻辑电路、框和功能的组合。软件可以存储在诸如存储器芯片或在处理器内实现的存储块等物理介质、诸如硬盘或软盘等磁性介质、以及诸如例如DVD及其数据变体、CD等光学介质上。物理介质是非瞬态介质。
本领域的技术人员应该明白,以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本申请范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请的保护范围应以所附权利要求为准。

Claims (10)

1.一种基于三维形态学的模型膨胀方法,其特征在于,包括以下步骤:
获取待膨胀模型,构建包围盒将所述待膨胀模型包围,并以固定大小的体素对所述包围盒进行采样得到体素集合;
将体素集合中的每一体素根据三维坐标保存到原始三维数组中,识别所述待膨胀模型内的每一三角面片,并将所述原始三维数组中所有包含三角面片的体素值记为1,不包含三角面片的体素值记为0得到模型体素三维数组;
设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中,定义由多个体素构建的立体结构元,使用所述立体结构元的中心体素依次遍历所述膨胀三维数组中的每一体素,遍历过程中若立体结构元所覆盖范围内的体素全部为0,则被遍历体素保持不变,若立体结构元所覆盖范围内至少存在一体素值为1的体素,则将被遍历体素的体素值设定为1进行膨胀,所述被遍历体素为与立体结构元中心体素重合的体素;
实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀,并基于膨胀三维数组进行模型重建得到待膨胀模型的膨胀结果。
2.根据权利要求1所述的一种基于三维形态学的模型膨胀方法,其特征在于,在“构建包围盒将所述待膨胀模型包围”步骤中,定义一个三维坐标系,使用所述待膨胀模型在三维坐标系中每一维度坐标的最大值和最小值来构建包围盒。
3.根据权利要求1所述的一种基于三维形态学的模型膨胀方法,其特征在于,在“识别所述待膨胀模型内的每一三角面片”步骤中,获取所述待膨胀模型中每一三角面片在原始三维数组中的位置,获取每一三角面片的至少一邻近体素,计算每一三角面片与对应的每一邻近体素中心点的邻近距离,若最短的邻近距离小于设定距离,则认为最短的邻近距离对应的邻近体素包含三角面片。
4.根据权利要求3所述的一种基于三维形态学的模型膨胀方法,其特征在于,在“计算每一三角面片与对应的每一邻近体素中心点的邻近距离”步骤中,获取三角面片的三个顶点坐标,使用三角面片的三个顶点坐标计算三角面片的法向量,获取邻近体素的中心点坐标,计算三角面片的法向量与对应邻近体素中心点坐标的内积得到中心点坐标到三角面片的有向距离,使用中心点坐标减去有向距离与三角面片法向量的乘积得到垂足点坐标,若所述垂足点在三角面片内则邻近体素与三角面片的邻近距离为邻近体素中心点到垂足点的距离,若所述垂足点不在三角面片内,则获取垂足点与三角面片的最短直线,使用最短直线与有向距离构成一个直角三角形,直角三角形的斜边为邻近体素与三角面片的邻近距离。
5.根据权利要求1所述的一种基于三维形态学的模型膨胀方法,其特征在于,在“设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中”步骤中,膨胀三维数组的边界相较于原始三维数组的边界向外膨胀,在将模型体素三维数组中的体素复制到膨胀三维数组后将膨胀三维数组中剩余体素位置的体素值设置为0。
6.根据权利要求1所述的一种基于三维形态学的模型膨胀方法,其特征在于,膨胀三维数组的边界相较于原始三维数组的边界向外膨胀的膨胀程度正相关于立体结构元的立体尺寸。
7.根据权利要求1所述的一种基于三维形态学的模型膨胀方法,其特征在于,在“实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀”步骤中,将所述膨胀后膨胀三维数组与膨胀前膨胀三维数组的比值乘以最大膨胀倍数的倒数作为膨胀系数。
8.一种基于三维形态学的模型膨胀装置,其特征在于,包括:
获取模块:获取待膨胀模型,构建包围盒将所述待膨胀模型包围,并以固定大小的体素对所述包围盒进行采样得到体素集合;
定义模块:将体素集合中的每一体素根据三维坐标保存到原始三维数组中,识别所述待膨胀模型内的每一三角面片,并将所述原始三维数组中所有包含三角面片的体素值记为1,不包含三角面片的体素值记为0得到模型体素三维数组;
膨胀模块:设置相较于原始三维数组膨胀的膨胀三维数组,将模型体素三维数组中所有体素根据三维坐标复制到膨胀三维数组中,定义由多个体素构建的立体结构元,使用所述立体结构元的中心体素依次遍历所述膨胀三维数组中的每一体素,遍历过程中若立体结构元所覆盖范围内的体素全部为0,则被遍历体素保持不变,若立体结构元所覆盖范围内至少存在一体素值为1的体素,则将被遍历体素的体素值设定为1进行膨胀,所述被遍历体素为与立体结构元中心体素重合的体素;
重建模块:实时计算所述膨胀三维数组的膨胀系数,若所述膨胀三维数组的膨胀系数大于等于膨胀阈值则停止膨胀,并基于膨胀三维数组进行模型重建得到待膨胀模型的膨胀结果。
9.一种电子装置,包括存储器和处理器,其特征在于,所述存储器中存储有计算机程序,所述处理器被设置为运行所述计算机程序以执行权利要求1-7任一所述的一种基于三维形态学的模型膨胀方法。
10.一种可读存储介质,其特征在于,所述可读存储介质中存储有计算机程序,所述计算机程序包括用于控制过程以执行过程的程序代码,所述过程包括根据权利要求1-7任一所述的一种基于三维形态学的模型膨胀方法。
CN202310945642.8A 2023-07-28 2023-07-28 一种基于三维形态学的模型膨胀方法、装置及应用 Active CN116862816B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310945642.8A CN116862816B (zh) 2023-07-28 2023-07-28 一种基于三维形态学的模型膨胀方法、装置及应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310945642.8A CN116862816B (zh) 2023-07-28 2023-07-28 一种基于三维形态学的模型膨胀方法、装置及应用

Publications (2)

Publication Number Publication Date
CN116862816A true CN116862816A (zh) 2023-10-10
CN116862816B CN116862816B (zh) 2024-01-26

Family

ID=88221598

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310945642.8A Active CN116862816B (zh) 2023-07-28 2023-07-28 一种基于三维形态学的模型膨胀方法、装置及应用

Country Status (1)

Country Link
CN (1) CN116862816B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101676576B1 (ko) * 2015-08-13 2016-11-15 삼성에스디에스 주식회사 3차원 모델의 복셀화 및 각 복셀의 속성 할당 장치 및 방법
US20170177746A1 (en) * 2015-12-17 2017-06-22 Fanuc Corporation Model generating device, position and orientation calculating device, and handling robot device
CN108022247A (zh) * 2016-11-03 2018-05-11 北京大学口腔医学院 基于牙周膜影像学解剖特征提取活体牙三维牙根形态方法
CN110826122A (zh) * 2019-10-12 2020-02-21 中广核工程有限公司 一种核电三维布置设计模型体素化方法及系统
CN112927334A (zh) * 2021-02-22 2021-06-08 中铁二院工程集团有限责任公司 一种基于gpu的三维模型快速体素化方法
CN113470180A (zh) * 2021-05-25 2021-10-01 杭州思看科技有限公司 三维网格重建方法、装置、电子装置和存储介质
CN113610784A (zh) * 2021-07-23 2021-11-05 湖北英库科技有限公司 一种肝段划分方法、系统、设备及存储介质
CN113793418A (zh) * 2021-09-24 2021-12-14 中国船舶工业系统工程研究院 基于线性几何体拟合的建筑物三维模型简化算法
CN114612318A (zh) * 2022-02-16 2022-06-10 西北大学 基于文物ct图像轮廓线的三维建模方法、系统及设备
CN114663637A (zh) * 2022-04-24 2022-06-24 杭州雅智医疗技术有限公司 三维牙齿模型倒凹区域填充方法、装置及应用
CN115830230A (zh) * 2022-11-24 2023-03-21 武汉大学 一种多视图几何构筑物三维建模系统及方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101676576B1 (ko) * 2015-08-13 2016-11-15 삼성에스디에스 주식회사 3차원 모델의 복셀화 및 각 복셀의 속성 할당 장치 및 방법
US20170177746A1 (en) * 2015-12-17 2017-06-22 Fanuc Corporation Model generating device, position and orientation calculating device, and handling robot device
CN108022247A (zh) * 2016-11-03 2018-05-11 北京大学口腔医学院 基于牙周膜影像学解剖特征提取活体牙三维牙根形态方法
CN110826122A (zh) * 2019-10-12 2020-02-21 中广核工程有限公司 一种核电三维布置设计模型体素化方法及系统
CN112927334A (zh) * 2021-02-22 2021-06-08 中铁二院工程集团有限责任公司 一种基于gpu的三维模型快速体素化方法
CN113470180A (zh) * 2021-05-25 2021-10-01 杭州思看科技有限公司 三维网格重建方法、装置、电子装置和存储介质
CN113610784A (zh) * 2021-07-23 2021-11-05 湖北英库科技有限公司 一种肝段划分方法、系统、设备及存储介质
CN113793418A (zh) * 2021-09-24 2021-12-14 中国船舶工业系统工程研究院 基于线性几何体拟合的建筑物三维模型简化算法
CN114612318A (zh) * 2022-02-16 2022-06-10 西北大学 基于文物ct图像轮廓线的三维建模方法、系统及设备
CN114663637A (zh) * 2022-04-24 2022-06-24 杭州雅智医疗技术有限公司 三维牙齿模型倒凹区域填充方法、装置及应用
CN115830230A (zh) * 2022-11-24 2023-03-21 武汉大学 一种多视图几何构筑物三维建模系统及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
XUAN SUN等: "Compositional Structure Recognition of 3D Building Models Through Volumetric Analysis", IEEE ACCESS, vol. 6, pages 33953 - 33968, XP011686514, DOI: 10.1109/ACCESS.2018.2842721 *
段煜: "基于激光点云的室外场景三维重建", 中国优秀硕士学位论文全文数据库 基础科学辑, no. 2, pages 005 - 232 *
温佩芝等: "三维重建网格模型的缺陷孔洞识别与修复方法", 计算机应用研究, no. 4, pages 280 - 284 *

Also Published As

Publication number Publication date
CN116862816B (zh) 2024-01-26

Similar Documents

Publication Publication Date Title
CN111581776B (zh) 一种基于几何重建模型的等几何分析方法
CN110033519B (zh) 基于隐式函数的三维建模方法、装置、系统及存储介质
WO2021142996A1 (zh) 基于图像分割的点云去噪方法、系统、装置和存储介质
CN101231761B (zh) 一种保持外观特征的网格模型简化方法
CN110084894B (zh) 三维模型的局部放大展示方法、装置与电子设备
Wang et al. Surface reconstruction from unoriented point clouds by a new triangle selection strategy
US20130265298A1 (en) Graphic processing method and apparatus
JP3265879B2 (ja) 3次元直交格子データの生成装置
CN111802978A (zh) 清扫控制方法、存储介质及扫地机器人
CN112561788A (zh) 一种bim模型的二维展开方法及纹理贴图方法、装置
US6901310B2 (en) Method and system for approximately reproducing the surface of a workpiece
US7324105B1 (en) Neighbor and edge indexing
US10937236B1 (en) Mesh smoothing for visual quality and analysis improvement
Branch et al. A hole-filling algorithm for triangular meshes using local radial basis function
CN116862816B (zh) 一种基于三维形态学的模型膨胀方法、装置及应用
CN110555910A (zh) 一种网格切割方法、装置、终端及存储介质
Chong et al. Automatic mesh-healing technique for model repair and finite element model generation
Amiri et al. Connectivity maps for subdivision surfaces
CN113592976B (zh) 地图数据的处理方法、装置、家用电器和可读存储介质
CN114155327A (zh) 一种智能多分辨率特征优化的三维重建方法及系统
CN112967396B (zh) 基于镜向反射的3d模型球面保面积参数化方法及系统
US7196703B1 (en) Primitive extension
Fang et al. 3D shape recovery of complex objects from multiple silhouette images
CN111932566B (zh) 一种模型轮廓图生成方法、装置和系统
CN114119407A (zh) 四角网格数据的去噪方法、装置、设备及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant