CN115564820A - 基于贪婪投影三角化的体积确定方法、系统、设备及介质 - Google Patents

基于贪婪投影三角化的体积确定方法、系统、设备及介质 Download PDF

Info

Publication number
CN115564820A
CN115564820A CN202211404837.3A CN202211404837A CN115564820A CN 115564820 A CN115564820 A CN 115564820A CN 202211404837 A CN202211404837 A CN 202211404837A CN 115564820 A CN115564820 A CN 115564820A
Authority
CN
China
Prior art keywords
point cloud
cloud data
target
point
curved surface
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
CN202211404837.3A
Other languages
English (en)
Other versions
CN115564820B (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.)
Wuhan Institute of Technology
Original Assignee
Wuhan Institute of Technology
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 Wuhan Institute of Technology filed Critical Wuhan Institute of Technology
Priority to CN202211404837.3A priority Critical patent/CN115564820B/zh
Publication of CN115564820A publication Critical patent/CN115564820A/zh
Application granted granted Critical
Publication of CN115564820B publication Critical patent/CN115564820B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • 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/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于贪婪投影三角化的体积确定方法、系统、设备及介质,涉及体积计算技术领域,方法包括:获取多个针对于目标对象的原始点云数据,对每个所述原始点云数据进行预处理,得到多个目标点云数据;根据多个所述目标点云数据进行点云补全,得到完整点云数据,所述完整点云数据包括多个点云补全后的点云数据;利用贪婪投影三角化算法,对所述完整点云数据进行曲面重建,得到针对于所述目标对象的目标曲面模型,所述目标曲面模型包含多个目标面,每个所述目标面包括多个目标点,每个所述目标点表征一个点云数据;根据所述目标曲面模型包含的多个所述目标面,确定所述目标对象的体积。

Description

基于贪婪投影三角化的体积确定方法、系统、设备及介质
技术领域
本发明涉及体积计算技术领域,具体涉及一种基于贪婪投影三角化的体积确定方法、系统、设备及介质。
背景技术
近年来,随着新能源汽车的兴起,人们越来越关注其中电池的重要组成部分——锂元素。锂元素主要生产于盐湖,但除了锂元素,人类已经从盐湖中大量开采盐,以及钾、镁、硼、溴等工业原料。盐作为一个国家的重要资源,不仅是人们日常生活的必需品,同时也是工业上的重要原材料,因此,合理规划并利用盐湖产盐区中的盐也就显得十分必要。
合理规划并利用产盐区中的盐的前提是获取其中盐堆的体积,将盐堆作为需要进行体积计算的目标对象,传统的获取盐堆体积的方法主要是通过人工对盐堆外形进行规整,然后拉尺测量盐堆的尺寸数据,最后根据测得的尺寸数据确定盐堆的体积。然而,在对盐堆进行规整过程中,往往会存在盐堆表面不水平等问题,在进行拉尺测量过程中,也容易存在测量误差,最终导致盐堆体积的计算精度不高。并且,一般情况下,产盐区中盐堆体积大、数量多,利用传统方法获取盐堆的体积不仅耗费人力,而且也难以快速、精确地测量出盐堆的体积,导致不能很好地满足工业需求。
因此,目前亟需一种确定目标对象体积的效率及准确率高的方法,使之能够快速、精确地计算出目标对象的体积。
发明内容
本发明所要解决的技术问题为:现有方法确定目标对象体积的人力成本高,效率及准确率低。为解决该技术问题,本发明提供了一种基于贪婪投影三角化的体积确定方法、系统、设备及介质。
本发明解决上述技术问题的技术方案如下:
一种基于贪婪投影三角化的体积确定方法,包括:
步骤S1,获取多个针对于目标对象的原始点云数据,对每个所述原始点云数据进行预处理,得到多个目标点云数据;
步骤S2,根据多个所述目标点云数据进行点云补全,得到完整点云数据,所述完整点云数据包括多个点云补全后的点云数据;
步骤S3,利用贪婪投影三角化算法,对所述完整点云数据进行曲面重建,得到针对于所述目标对象的目标曲面模型,所述目标曲面模型包含多个目标面,每个所述目标面包括多个目标点,每个所述目标点表征一个点云数据;
步骤S4,根据所述目标曲面模型包含的多个所述目标面,确定所述目标对象的体积。
本发明的有益效果是:通过对原始点云数据进行预处理,可降低点云密度,缩短后续曲面重建时间,提高目标对象体积的计算效率;通过对目标点云数据进行点云补全,可消除激光雷达只能采集目标对象表面点云的缺陷,从而得到完整点云数据,便于后续构建针对于目标对象的目标曲面模型;通过对完整点云数据使用贪婪投影三角化的曲面重建方法,可生成更精确的曲面模型,保证了目标对象体积计算的准确率;基于目标曲面模型中包含的目标点,利用几何体计算公式即可确定目标对象的体积,计算速度快,且计算精度高。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述步骤S1中,对每个所述原始点云数据进行预处理,得到多个目标点云数据,包括:
对于每个所述原始点云数据,通过VoxelGrid滤波器对所述原始点云数据进行下采样处理,得到所述原始点云数据对应的目标点云数据。
采用上述进一步方案的有益效果是:通过利用VoxelGrid滤波器对原始点云数据进行下采样处理,既能筛除不属于目标对象的点云噪声,又能降低点云密度、缩短后续曲面重建时间,提高了目标对象体积的计算效率。
进一步,每个所述目标点云数据包含各自在预先建立的多维空间坐标系中对应的点云位置信息;
所述步骤S2包括:
从多个所述目标点云数据中提取多个边界点云数据,每个所述边界点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
根据多个所述目标点云数据和多个所述边界点云数据,确定表面点云数据,每个所述表面点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
对于每个所述表面点云数据,以所述目标对象的底部作为投影方向,对所述表面点云数据进行垂直投影,得到所述表面点云数据对应的底部点云数据,所述底部点云数据表征对应的所述表面点云数据投影到所述目标对象底面的点的点云数据,每个所述底部点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
对于每个所述边界点云数据,根据所述边界点云数据对应的点云位置信息和预设的步长,得到所述边界点云数据对应的多个四周点云数据;
根据多个所述目标点云数据、多个所述底部点云数据和多个所述四周点云数据,得到完整点云数据;
所述目标对象底面的各个点在所述多维空间坐标系中对应有第一点位置信息,所述第一点位置信息包括横坐标值、纵坐标值和高度值,每个所述边界点云数据各自对应的点云位置信息包括点云横坐标、点云纵坐标和点云竖坐标;
所述根据所述边界点云数据对应的点云位置信息和预设的步长,得到所述边界点云数据对应的多个四周点云数据,包括:
将所述边界点云数据对应的点云横坐标作为第一坐标值,将所述边界点云数据对应的点云纵坐标作为第二坐标值;
获取所述目标对象底面中横坐标值等于所述第一坐标值、纵坐标值等于所述第二坐标值的点的高度值,将所述高度值作为目标高度值;
根据所述边界点云数据对应的点云竖坐标、所述目标高度值和所述步长,确定所述边界点云数据对应的四周点云数据的个数;
根据所述边界点云数据对应的点云位置信息、所述步长和所述个数,得到多个四周点云数据,每个所述四周点云数据包含各自在所述多维空间坐标系中对应的点云位置信息。
采用上述进一步方案的有益效果是:通过对目标点云数据进行点云补全,可消除激光雷达只能采集目标对象表面点云的缺陷,从而得到完整点云数据,为后续构建针对于目标对象的目标曲面模型奠定了基础。
进一步,所述步骤S3包括:
对于所述完整点云数据中的每个所述点云数据,通过预设的最邻近参数计算所述点云数据对应的协方差矩阵,根据所述协方差矩阵,计算所述点云数据对应的法向量;
对于所述完整点云数据中的每个所述点云数据,根据所述点云数据对应的所述协方差矩阵和所述法向量,确定所述点云数据对应的法向量的方向;
利用贪婪投影三角化算法,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型。
采用上述进一步方案的有益效果是:利用贪婪投影三角化算法对完整点云数据进行曲面重建,保证了构建出的目标曲面模型的完整性,且构建出的目标曲面模型精度高。
进一步,所述利用贪婪投影三角化算法,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型,包括:
步骤A1,确定迭代次数
Figure 274777DEST_PATH_IMAGE001
步骤A2,在超参数空间中随机选择参数集合,所述贪婪投影三角化算法使用所述参数集合,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到所述参数集合对应的初始曲面模型;
步骤A3,重复
Figure 758628DEST_PATH_IMAGE002
次所述步骤A2,得到
Figure 600683DEST_PATH_IMAGE001
个初始曲面模型;
步骤A4,对于每个所述参数集合,确定所述参数集合对应的所述初始曲面模型与预先构建好的标准曲面模型的相似度;
步骤A5,根据所述相似度,确定目标参数集合;
步骤A6,所述贪婪投影三角化算法使用所述目标参数集合,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型。
采用上述进一步方案的有益效果是:通过利用随机搜索算法进行超参数调优,能更详细地搜索曲面重建过程中的参数,便于得到更好的曲面模型。
进一步,每个所述目标点包含各自在预先建立的多维空间坐标系中对应的点云位置信息,每个所述目标点各自对应的点云位置信息包括目标点横坐标、目标点纵坐标和目标点竖坐标;
所述步骤S4包括:
根据多个所述目标点,构建多个目标面,每个所述目标面由多个所述目标点的连线组成;
对于每个所述目标面,确定组成所述目标面的每个所述目标点各自对应的点云位置信息;
对于每个所述目标面,以所述目标对象的底部作为投影方向,对所述目标面进行垂直投影,得到所述目标面对应的多个底面点,所述底面点表征对应的所述目标点投影到所述目标对象底面的点,每个所述底面点包含各自在所述多维空间坐标系中对应的第二点位置信息,所述第二点位置信息包括横坐标值、纵坐标值和高度值;
对于每个目标面,根据组成所述目标面的多个所述目标点和所述目标面对应的多个所述底面点,构建几何体,根据所述目标点对应的点云位置信息和所述底面点对应的第二点位置信息,计算所述几何体的体积,根据所述几何体的体积确定所述目标面对应的体积;
将多个所述目标面对应的体积相加,得到所述目标对象的体积。
采用上述进一步方案的有益效果是:通过根据目标点构建目标面,对目标面进行垂直投影,再基于投影结果计算目标面的体积,从而确定目标对象的体积,计算速度快且计算精度高,避免了因目标对象表面不水平而导致计算精度低的问题。
为解决上述技术问题,本发明还提供了一种基于贪婪投影三角化的体积确定系统,包括:
第一数据处理模块,用于获取多个针对于目标对象的原始点云数据,对每个所述原始点云数据进行预处理,得到多个目标点云数据;
第二数据处理模块,用于根据多个所述目标点云数据进行点云补全,得到完整点云数据,所述完整点云数据包括多个点云补全后的点云数据;
模型构建模块,用于利用贪婪投影三角化算法,对所述完整点云数据进行曲面重建,得到针对于所述目标对象的目标曲面模型,所述目标曲面模型包含多个目标面,每个所述目标面包括多个目标点,每个所述目标点表征一个点云数据;
体积确定模块,根据所述目标曲面模型包含的多个所述目标面,确定所述目标对象的体积。
为解决上述技术问题,本发明还提供了一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如前所述的基于贪婪投影三角化的体积确定方法。
为解决上述技术问题,本发明还提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如前所述的基于贪婪投影三角化的体积确定方法。
附图说明
图1为本发明中基于贪婪投影三角化的体积确定方法的流程示意图;
图2为本发明实施例中不规则几何体的示意图;
图3为本发明中基于贪婪投影三角化的体积确定系统的结构示意图;
图4为本发明中电子设备的结构示意图。
具体实施方式
以下对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。为了便于对本申请方案的理解,对本申请中涉及的名词和原理进行解释,在本发明中,目标对象可以是盐堆,也可以是矿石、物料等可获取其点云数据、且自身的底面具有一定面积的实体,实体由多个面闭合包裹形成,每个面由多个点组成。组成所述目标对象的各个点在多维空间坐标系中对应有位置信息,所述完整点云数据包含的每个点云数据在多维空间坐标系中对应有位置信息(包括点云横坐标、点云纵坐标和点云竖坐标),所述目标对象的底面表征所述完整点云数据中,最小的点云竖坐标对应的所述点云数据在所述目标对象上所在的面,所述目标对象的底面即为所述目标对象的底部。
实施例一
为解决现有技术中存在的技术问题,本实施例提供了一种基于贪婪投影三角化的体积确定方法,如图1所示,包括:
步骤S1,获取多个针对于目标对象的原始点云数据,对每个所述原始点云数据进行预处理,得到多个目标点云数据;
步骤S2,根据多个所述目标点云数据进行点云补全,得到完整点云数据,所述完整点云数据包括多个点云补全后的点云数据;
步骤S3,利用贪婪投影三角化算法,对所述完整点云数据进行曲面重建,得到针对于所述目标对象的目标曲面模型,所述目标曲面模型包含多个目标面,每个所述目标面包括多个目标点,每个所述目标点表征一个点云数据;
步骤S4,根据所述目标曲面模型包含的多个所述目标面,确定所述目标对象的体积。
本实施例中,以盐堆作为目标对象,通过激光雷达在盐堆上方采集多个原始点云数据,每个所述原始点云数据包含X、Y、Z三个维度的位置信息;在通过激光雷达获取原始点云数据后,人工限定点云区域,将不属于点云区域的原始点云数据(即不属于盐堆部分的原始点云数据)删除,得到最终用于进行预处理的原始点云数据。
其中,所述步骤S1中,对每个所述原始点云数据进行预处理,得到多个目标点云数据,包括:
对于每个所述原始点云数据,通过VoxelGrid滤波器对所述原始点云数据进行下采样处理,得到所述原始点云数据对应的目标点云数据。本实施例中,所述VoxelGrid滤波器进行下采样处理的体素大小为1cm×1cm×1cm的立方体。
其中,每个所述目标点云数据包含各自在预先建立的多维空间坐标系中对应的点云位置信息;本实施例中,所述多维空间坐标系为三维坐标系;
所述步骤S2包括:
从多个所述目标点云数据中提取多个边界点云数据,每个所述边界点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
根据多个所述目标点云数据和多个所述边界点云数据,确定表面点云数据,每个所述表面点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
对于每个所述表面点云数据,以盐堆的底部作为投影方向,对所述表面点云数据进行垂直投影,得到所述表面点云数据对应的底部点云数据,所述底部点云数据表征对应的所述表面点云数据投影到盐堆底面的点的点云数据,每个所述底部点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
对于每个所述边界点云数据,根据所述边界点云数据对应的点云位置信息和预设的步长,得到所述边界点云数据对应的多个四周点云数据;
根据多个所述目标点云数据、多个所述底部点云数据和多个所述四周点云数据,得到完整点云数据。
其中,对于所述边界点云数据的提取,本实施例采用边界提取算法BoundaryEstimation提取点云边界,得到多个边界点云数据,以补全盐堆不能被激光雷达获取到的部分点云。其中,BoundaryEstimation中表示某个点的最近邻点的数量的参数的值,可以为20;
对于所述表面点云数据的确定,本实施例通过遍历所述目标点云数据,将不属于所述边界点云数据的目标点云数据确定为表面点云数据;
盐堆底面的各个点在所述多维空间坐标系中对应有第一点位置信息,所述第一点位置信息包括横坐标值、纵坐标值和高度值,每个所述边界点云数据各自对应的点云位置信息包括点云横坐标、点云纵坐标和点云竖坐标;
对于所述完整点云数据的确定,本实施例通过将多个所述目标点云数据、多个所述底部点云数据和多个所述四周点云数据合并的方式,确定完整点云数据,目的是得到一个盐堆的三维模型,而不是只有盐堆表面的曲面模型;
所述根据所述边界点云数据对应的点云位置信息和预设的步长,得到所述边界点云数据对应的多个四周点云数据,包括:
将所述边界点云数据对应的点云横坐标作为第一坐标值,将所述边界点云数据对应的点云纵坐标作为第二坐标值;
获取盐堆底面中横坐标值等于所述第一坐标值、纵坐标值等于所述第二坐标值的点的高度值,将所述高度值作为目标高度值;
根据所述边界点云数据对应的点云竖坐标、所述目标高度值和所述步长,确定所述边界点云数据对应的四周点云数据的个数;
根据所述边界点云数据对应的点云位置信息、所述步长和所述个数,得到多个四周点云数据,每个所述四周点云数据包含各自在所述多维空间坐标系中对应的点云位置信息,每个所述四周点云数据各自对应的点云位置信息包括点云横坐标、点云纵坐标和点云竖坐标。
例如,将所述步长设置为1mm,表示沿盐堆底面的方向,高度每间隔1mm,便生成1个边界点云数据对应的四周点云数据,直至生成的四周点云数据对应的点云竖坐标小于或等于基于该边界点云数据对应的点云位置信息所确定的盐堆底面的点的高度值。
其中,所述步骤S3包括:
对于所述完整点云数据中的每个所述点云数据,通过预设的最邻近参数计算所述点云数据对应的协方差矩阵,根据所述协方差矩阵,计算所述点云数据对应的法向量;
对于所述完整点云数据中的每个所述点云数据,根据所述点云数据对应的所述协方差矩阵和所述法向量,确定所述点云数据对应的法向量的方向;
利用贪婪投影三角化算法,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于盐堆的目标曲面模型。
其中,所述利用贪婪投影三角化算法,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于盐堆的目标曲面模型,包括:
步骤A1,确定迭代次数
Figure 995892DEST_PATH_IMAGE001
步骤A2,在超参数空间中随机选择参数集合,所述贪婪投影三角化算法使用所述参数集合,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到所述参数集合对应的初始曲面模型;
步骤A3,重复
Figure 872581DEST_PATH_IMAGE003
次所述步骤A2,得到
Figure 757360DEST_PATH_IMAGE001
个初始曲面模型;
步骤A4,对于每个所述参数集合,确定所述参数集合对应的所述初始曲面模型与预先构建好的标准曲面模型的相似度;
步骤A5,根据所述相似度,确定目标参数集合;
步骤A6,所述贪婪投影三角化算法使用所述目标参数集合,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于盐堆的目标曲面模型。
本实施例中,利用随机搜索算法RandomizedSearch在超参数空间中随机选择参数集合。具体地,RandomizedSearch通过在超参数空间中搜索几十至几百个点,通过比对曲面重建后的初始曲面模型和标准曲面模型,将相似度最高的参数集合作为最优的参数集合(即目标参数集合),根据最优的参数集合,从而确定目标曲面模型。
其中,每个所述目标点包含各自在所述多维空间坐标系中对应的点云位置信息,每个所述目标点各自对应的点云位置信息包括目标点横坐标、目标点纵坐标和目标点竖坐标;所述目标面可以为三角面;
以三角面为例,所述步骤S4包括:
根据多个所述目标点,通过所述贪婪投影三角化算法构建多个三角面,每个所述三角面由3个所述目标点的连线组成;
对于每个所述三角面,确定组成所述三角面的3个所述目标点各自对应的点云位置信息;
对于每个所述三角面,以盐堆的底部作为投影方向,对所述三角面进行垂直投影,得到所述三角面对应的3个底面点,所述底面点表征对应的所述目标点投影到盐堆底面的点,每个所述底面点包含各自在所述多维空间坐标系中对应的第二点位置信息,所述第二点位置信息包括横坐标值、纵坐标值和高度值;
对于每个三角面,通过将组成所述三角面的3个所述目标点连接成三角形、将所述三角面对应的3个所述底面点连接成三角形、将每个所述目标点与各自对应的所述底面点连接,构建出所述三角面对应的几何体,根据所述目标点对应的点云位置信息和所述底面点对应的第二点位置信息,计算所述几何体的体积,所述三角面对应的体积等于所述几何体的体积;
将多个所述三角面对应的体积相加,得到盐堆的体积。
其中,所述几何体可以是规则的三棱柱,也可能是不规则几何体。例如,组成某一 三角面的3个目标点(每个目标点即对应完整点云数据中包含的一个经点云补全后的点云 数据)分别为H、G、F,将目标点H、G、F分别向盐堆底面垂直投影,得到3个底面点,分别为A、B、 C,如图2所示。将目标点H、G、F连接成三角形,将底面点A、B、C连接成三角形,将目标点H与对 应的底面点A连接、将目标点G与对应的底面点B连接、将目标点F与对应的底面点C连接,从 而构建出该三角面对应的不规则几何体H-G-F-C-B-A。在多维空间坐标系中,目标点H、G、F, 以及底面点A、B、C各自对应有位置信息,根据各个目标点、各个底面点的坐标位置,通过添 加辅助线的方式,新增两个辅助点,分别为D、E,不规则几何体H-G-F-C-B-A即由三棱锥H- DFG、三棱锥G-DEF和三棱柱ABC-DEF组成。对于三棱锥H-DFG,通过第一公式计算三棱锥H- DFG的体积
Figure 24394DEST_PATH_IMAGE004
,第一公式为:
Figure 787076DEST_PATH_IMAGE005
上式中,
Figure 518271DEST_PATH_IMAGE006
表示由三棱锥H-DFG的四个顶点坐标,通过行列式计算三棱锥H-DFG的 体积;
Figure 511635DEST_PATH_IMAGE007
表示三棱锥H-DFG各顶点在所述多维空间坐标系 中对应的位置信息;
三棱锥G-DEF体积
Figure 328281DEST_PATH_IMAGE008
的计算方法与三棱锥H-DFG体积
Figure 65293DEST_PATH_IMAGE004
的计算方法相同,相同之 处不再赘述;
通过第二公式计算三棱柱ABC-DEF的体积
Figure 650995DEST_PATH_IMAGE009
,第二公式为:
Figure 641692DEST_PATH_IMAGE010
上式中,
Figure 883317DEST_PATH_IMAGE011
表示三棱柱ABC-DEF的底面积,
Figure 220757DEST_PATH_IMAGE012
表示三棱柱ABC-DEF的高;
通过第三公式计算不规则几何体H-G-F-C-B-A的体积
Figure 926545DEST_PATH_IMAGE013
,第三公式为:
Figure 261712DEST_PATH_IMAGE014
不规则几何体H-G-F-C-B-A的体积
Figure 52950DEST_PATH_IMAGE013
即为该三角面对应的体积。
实施例二
基于与实施例一中基于贪婪投影三角化的体积确定方法相同的原理,本实施例提供了一种基于贪婪投影三角化的体积确定系统,如图3所示,包括:
第一数据处理模块,用于获取多个针对于目标对象的原始点云数据,对每个所述原始点云数据进行预处理,得到多个目标点云数据;
第二数据处理模块,用于根据多个所述目标点云数据进行点云补全,得到完整点云数据,所述完整点云数据包括多个点云补全后的点云数据;
模型构建模块,用于利用贪婪投影三角化算法,对所述完整点云数据进行曲面重建,得到针对于所述目标对象的目标曲面模型,所述目标曲面模型包含多个目标面,每个所述目标面包括多个目标点,每个所述目标点表征一个点云数据;
体积确定模块,根据所述目标曲面模型包含的多个所述目标面,确定所述目标对象的体积。
其中,所述第一数据处理模块包括数据预处理单元,所述数据预处理单元用于通过VoxelGrid滤波器对每个所述原始点云数据进行下采样处理,得到每个所述原始点云数据各自对应的目标点云数据。
其中,所述第二数据处理模块包括:
第一处理单元,用于根据多个所述目标点云数据和多个所述边界点云数据,确定表面点云数据,每个所述表面点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
第二处理单元,用于确定每个所述表面点云数据各自对应的底部点云数据,即:
对于每个所述表面点云数据,以所述目标对象的底部作为投影方向,对所述表面点云数据进行垂直投影,得到所述表面点云数据对应的底部点云数据,所述底部点云数据表征对应的所述表面点云数据投影到所述目标对象底面的点的点云数据,每个所述底部点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
第三处理单元,用于确定每个所述边界点云数据各自对应的四周点云数据,即:
对于每个所述边界点云数据,根据所述边界点云数据对应的点云位置信息和预设的步长,得到所述边界点云数据对应的多个四周点云数据;
第四处理单元,用于根据多个所述目标点云数据、多个所述底部点云数据和多个所述四周点云数据,得到完整点云数据;
对于每个所述边界点云数据,所述第三处理单元具体用于:
将所述边界点云数据对应的点云横坐标作为第一坐标值,将所述边界点云数据对应的点云纵坐标作为第二坐标值;
获取所述目标对象底面中横坐标值等于所述第一坐标值、纵坐标值等于所述第二坐标值的点的高度值,将所述高度值作为目标高度值;
根据所述边界点云数据对应的点云竖坐标、所述目标高度值和所述步长,确定所述边界点云数据对应的四周点云数据的个数;
根据所述边界点云数据对应的点云位置信息、所述步长和所述个数,得到多个四周点云数据,每个所述四周点云数据包含各自在所述多维空间坐标系中对应的点云位置信息。
其中,所述模型构建模块具体用于:
对于所述完整点云数据中的每个所述点云数据,通过预设的最邻近参数计算所述点云数据对应的协方差矩阵,根据所述协方差矩阵,计算所述点云数据对应的法向量;
对于所述完整点云数据中的每个所述点云数据,根据所述点云数据对应的所述协方差矩阵和所述法向量,确定所述点云数据对应的法向量的方向;
利用贪婪投影三角化算法,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型。
其中,所述目标模型构建模块在利用贪婪投影三角化算法,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型时,包括:
第五处理单元,用于确定迭代次数
Figure 131765DEST_PATH_IMAGE015
第六处理单元,用于在超参数空间中随机选择参数集合,所述贪婪投影三角化算法使用所述参数集合,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到所述参数集合对应的初始曲面模型;
第七处理单元,用于重复执行
Figure 193524DEST_PATH_IMAGE016
次所述第六处理单元的处理过程,得到
Figure 27488DEST_PATH_IMAGE015
个初始曲面模型;
第八处理单元,用于确定每个所述参数集合对应的相似度;具体地,对于每个所述参数集合,确定所述参数集合对应的所述初始曲面模型与预先构建好的标准曲面模型的相似度;
第九处理单元,用于根据所述相似度,确定目标参数集合;
第十处理单元,用于所述贪婪投影三角化算法使用所述目标参数集合,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型。
其中,所述体积确定模块具体用于:
根据多个所述目标点,构建多个目标面,每个所述目标面由多个所述目标点的连线组成;
对于每个所述目标面,确定组成所述目标面的每个所述目标点各自对应的点云位置信息;
对于每个所述目标面,以所述目标对象的底部作为投影方向,对所述目标面进行垂直投影,得到所述目标面对应的多个底面点,所述底面点表征对应的所述目标点投影到所述目标对象底面的点,每个所述底面点包含各自在所述多维空间坐标系中对应的第二点位置信息,所述第二点位置信息包括横坐标值、纵坐标值和高度值;
对于每个目标面,根据组成所述目标面的多个所述目标点和所述目标面对应的多个所述底面点,构建几何体,根据所述目标点对应的点云位置信息和所述底面点对应的第二点位置信息,计算所述几何体的体积,根据所述几何体的体积确定所述目标面对应的体积;
将多个所述目标面对应的体积相加,得到所述目标对象的体积。
实施例三
基于与实施例一中基于贪婪投影三角化的体积确定方法相同的原理,本实施例提供了一种电子设备,如图4所示,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如实施例一所述的基于贪婪投影三角化的体积确定方法。
实施例四
基于与实施例一中基于贪婪投影三角化的体积确定方法相同的原理,本实施例提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如实施例一所述的基于贪婪投影三角化的体积确定方法。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (10)

1.一种基于贪婪投影三角化的体积确定方法,其特征在于,包括:
步骤S1,获取多个针对于目标对象的原始点云数据,对每个所述原始点云数据进行预处理,得到多个目标点云数据;
步骤S2,根据多个所述目标点云数据进行点云补全,得到完整点云数据,所述完整点云数据包括多个点云补全后的点云数据;
步骤S3,利用贪婪投影三角化算法,对所述完整点云数据进行曲面重建,得到针对于所述目标对象的目标曲面模型,所述目标曲面模型包含多个目标面,每个所述目标面包括多个目标点,每个所述目标点表征一个点云数据;
步骤S4,根据所述目标曲面模型包含的多个所述目标面,确定所述目标对象的体积。
2.根据权利要求1所述的方法,其特征在于,所述步骤S1中,对每个所述原始点云数据进行预处理,得到多个目标点云数据,包括:
对于每个所述原始点云数据,通过VoxelGrid滤波器对所述原始点云数据进行下采样处理,得到所述原始点云数据对应的目标点云数据。
3.根据权利要求1所述的方法,其特征在于,每个所述目标点云数据包含各自在预先建立的多维空间坐标系中对应的点云位置信息;
所述步骤S2包括:
从多个所述目标点云数据中提取多个边界点云数据,每个所述边界点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
根据多个所述目标点云数据和多个所述边界点云数据,确定表面点云数据,每个所述表面点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
对于每个所述表面点云数据,以所述目标对象的底部作为投影方向,对所述表面点云数据进行垂直投影,得到所述表面点云数据对应的底部点云数据,所述底部点云数据表征对应的所述表面点云数据投影到所述目标对象底面的点的点云数据,每个所述底部点云数据包含各自在所述多维空间坐标系中对应的点云位置信息;
对于每个所述边界点云数据,根据所述边界点云数据对应的点云位置信息和预设的步长,得到所述边界点云数据对应的多个四周点云数据;
根据多个所述目标点云数据、多个所述底部点云数据和多个所述四周点云数据,得到完整点云数据。
4.根据权利要求3所述的方法,其特征在于,所述目标对象底面的各个点在所述多维空间坐标系中对应有第一点位置信息,所述第一点位置信息包括横坐标值、纵坐标值和高度值,每个所述边界点云数据各自对应的点云位置信息包括点云横坐标、点云纵坐标和点云竖坐标;
所述根据所述边界点云数据对应的点云位置信息和预设的步长,得到所述边界点云数据对应的多个四周点云数据,包括:
将所述边界点云数据对应的点云横坐标作为第一坐标值,将所述边界点云数据对应的点云纵坐标作为第二坐标值;
获取所述目标对象底面中横坐标值等于所述第一坐标值、纵坐标值等于所述第二坐标值的点的高度值,将所述高度值作为目标高度值;
根据所述边界点云数据对应的点云竖坐标、所述目标高度值和所述步长,确定所述边界点云数据对应的四周点云数据的个数;
根据所述边界点云数据对应的点云位置信息、所述步长和所述个数,得到多个四周点云数据,每个所述四周点云数据包含各自在所述多维空间坐标系中对应的点云位置信息。
5.根据权利要求1所述的方法,其特征在于,所述步骤S3包括:
对于所述完整点云数据中的每个所述点云数据,通过预设的最邻近参数计算所述点云数据对应的协方差矩阵,根据所述协方差矩阵,计算所述点云数据对应的法向量;
对于所述完整点云数据中的每个所述点云数据,根据所述点云数据对应的所述协方差矩阵和所述法向量,确定所述点云数据对应的法向量的方向;
利用贪婪投影三角化算法,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型。
6.根据权利要求5所述的方法,其特征在于,所述利用贪婪投影三角化算法,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型,包括:
步骤A1,确定迭代次数
Figure 898079DEST_PATH_IMAGE001
步骤A2,在超参数空间中随机选择参数集合,所述贪婪投影三角化算法使用所述参数集合,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到所述参数集合对应的初始曲面模型;
步骤A3,重复
Figure 607015DEST_PATH_IMAGE002
次所述步骤A2,得到
Figure 959499DEST_PATH_IMAGE001
个初始曲面模型;
步骤A4,对于每个所述参数集合,确定所述参数集合对应的所述初始曲面模型与预先构建好的标准曲面模型的相似度;
步骤A5,根据所述相似度,确定目标参数集合;
步骤A6,所述贪婪投影三角化算法使用所述目标参数集合,根据所述完整点云数据中每个所述点云数据对应的法向量和所述法向量的方向进行曲面重建,得到针对于所述目标对象的目标曲面模型。
7.根据权利要求1所述的方法,其特征在于,每个所述目标点包含各自在预先建立的多维空间坐标系中对应的点云位置信息,每个所述目标点各自对应的点云位置信息包括目标点横坐标、目标点纵坐标和目标点竖坐标;
所述步骤S4包括:
根据多个所述目标点,构建多个目标面,每个所述目标面由多个所述目标点的连线组成;
对于每个所述目标面,确定组成所述目标面的每个所述目标点各自对应的点云位置信息;
对于每个所述目标面,以所述目标对象的底部作为投影方向,对所述目标面进行垂直投影,得到所述目标面对应的多个底面点,所述底面点表征对应的所述目标点投影到所述目标对象底面的点,每个所述底面点包含各自在所述多维空间坐标系中对应的第二点位置信息,所述第二点位置信息包括横坐标值、纵坐标值和高度值;
对于每个目标面,根据组成所述目标面的多个所述目标点和所述目标面对应的多个所述底面点,构建几何体,根据所述目标点对应的点云位置信息和所述底面点对应的第二点位置信息,计算所述几何体的体积,根据所述几何体的体积确定所述目标面对应的体积;
将多个所述目标面对应的体积相加,得到所述目标对象的体积。
8.一种基于贪婪投影三角化的体积确定系统,其特征在于,包括:
第一数据处理模块,用于获取多个针对于目标对象的原始点云数据,对每个所述原始点云数据进行预处理,得到多个目标点云数据;
第二数据处理模块,用于根据多个所述目标点云数据进行点云补全,得到完整点云数据,所述完整点云数据包括多个点云补全后的点云数据;
模型构建模块,用于利用贪婪投影三角化算法,对所述完整点云数据进行曲面重建,得到针对于所述目标对象的目标曲面模型,所述目标曲面模型包含多个目标面,每个所述目标面包括多个目标点,每个所述目标点表征一个点云数据;
体积确定模块,根据所述目标曲面模型包含的多个所述目标面,确定所述目标对象的体积。
9.一种电子设备,其特征在于,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如权利要求1至7中任一项所述的基于贪婪投影三角化的体积确定方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至7中任一项所述的基于贪婪投影三角化的体积确定方法。
CN202211404837.3A 2022-11-10 2022-11-10 基于贪婪投影三角化的体积确定方法、系统、设备及介质 Active CN115564820B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211404837.3A CN115564820B (zh) 2022-11-10 2022-11-10 基于贪婪投影三角化的体积确定方法、系统、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211404837.3A CN115564820B (zh) 2022-11-10 2022-11-10 基于贪婪投影三角化的体积确定方法、系统、设备及介质

Publications (2)

Publication Number Publication Date
CN115564820A true CN115564820A (zh) 2023-01-03
CN115564820B CN115564820B (zh) 2023-03-21

Family

ID=84769905

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211404837.3A Active CN115564820B (zh) 2022-11-10 2022-11-10 基于贪婪投影三角化的体积确定方法、系统、设备及介质

Country Status (1)

Country Link
CN (1) CN115564820B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116188803A (zh) * 2023-04-23 2023-05-30 武汉工程大学 一种基于极坐标的边界点云提取方法、系统、设备及介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6996505B1 (en) * 2000-06-21 2006-02-07 Raindrop Geomagic, Inc. Methods, apparatus and computer program products for automatically generating nurbs models of triangulated surfaces using homeomorphisms
CN104794748A (zh) * 2015-03-17 2015-07-22 上海海洋大学 基于Kinect视觉技术的三维空间地图构建方法
CN106504332A (zh) * 2016-10-19 2017-03-15 未来科技(襄阳)有限公司 三维点云的曲面重建方法和装置
CN110176061A (zh) * 2019-04-30 2019-08-27 中科恒运股份有限公司 一种三维重建中人体表面重构方法
CN111265343A (zh) * 2020-01-10 2020-06-12 大连理工大学 假肢接受腔自适应匹配与自动生成方法
CN112950771A (zh) * 2021-04-02 2021-06-11 中国铁建重工集团股份有限公司 一种围岩结构的检测方法、系统及相关组件
CN113776429A (zh) * 2021-07-21 2021-12-10 武汉市异方体科技有限公司 一种车载物料体积测量方法
CN114419130A (zh) * 2021-12-22 2022-04-29 中国水利水电第七工程局有限公司 一种基于图像特征和三维点云技术的散料体积测量方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6996505B1 (en) * 2000-06-21 2006-02-07 Raindrop Geomagic, Inc. Methods, apparatus and computer program products for automatically generating nurbs models of triangulated surfaces using homeomorphisms
CN104794748A (zh) * 2015-03-17 2015-07-22 上海海洋大学 基于Kinect视觉技术的三维空间地图构建方法
CN106504332A (zh) * 2016-10-19 2017-03-15 未来科技(襄阳)有限公司 三维点云的曲面重建方法和装置
CN110176061A (zh) * 2019-04-30 2019-08-27 中科恒运股份有限公司 一种三维重建中人体表面重构方法
CN111265343A (zh) * 2020-01-10 2020-06-12 大连理工大学 假肢接受腔自适应匹配与自动生成方法
CN112950771A (zh) * 2021-04-02 2021-06-11 中国铁建重工集团股份有限公司 一种围岩结构的检测方法、系统及相关组件
CN113776429A (zh) * 2021-07-21 2021-12-10 武汉市异方体科技有限公司 一种车载物料体积测量方法
CN114419130A (zh) * 2021-12-22 2022-04-29 中国水利水电第七工程局有限公司 一种基于图像特征和三维点云技术的散料体积测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
杨玉泽;林文树;: "基于激光点云数据的树木枝叶分割和三维重建" *
王婷玉: "《采空区三角形投影体积算法及其应用》", 《中国地质灾害与防治学报》 *
胡燕威: "《基于激光雷达的空间物体三维建模与体积计算》", 《中国硕士学位优秀论文全文数据库》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116188803A (zh) * 2023-04-23 2023-05-30 武汉工程大学 一种基于极坐标的边界点云提取方法、系统、设备及介质
CN116188803B (zh) * 2023-04-23 2023-09-15 武汉工程大学 一种基于极坐标的边界点云提取方法、系统、设备及介质

Also Published As

Publication number Publication date
CN115564820B (zh) 2023-03-21

Similar Documents

Publication Publication Date Title
US8300920B2 (en) System and method for measuring form and position tolerances of an object
JP4900204B2 (ja) 物体認識方法
CN114998338B (zh) 一种基于激光雷达点云的矿山开采量计算方法
CN103278115B (zh) 一种基于dem计算淤地坝淤积量的方法及系统
CN115372989A (zh) 基于激光雷达的越野自动小车长距离实时定位系统及方法
CN112215958B (zh) 一种基于分布式计算的激光雷达点云数据投影方法
CN104616349A (zh) 基于局部曲面变化因子的散乱点云数据精简处理方法
CN115564820B (zh) 基于贪婪投影三角化的体积确定方法、系统、设备及介质
CN112819962A (zh) 数字图像相关中非均匀网格划分及局部网格疏密方法
TW201344633A (zh) 曲面量測系統及方法
CN112017219B (zh) 一种激光点云配准方法
CN107492120A (zh) 点云配准方法
CN103315770A (zh) 基于三维超声图像的目标参数处理方法及系统
CN112232248A (zh) 一种多线LiDAR点云数据平面特征的提取方法及装置
CN107993242B (zh) 基于机载LiDAR点云数据缺失区域边界提取方法
Maurya et al. Performance of greedy triangulation algorithm on reconstruction of coastal dune surface
CN107748834A (zh) 一种计算起伏观测面磁场的快速、高精度数值模拟方法
Liu et al. Accurate volume calculation driven by Delaunay triangulation for coal measurement
CN115131571A (zh) 一种基于点云预处理六领域的建筑物局部特征点识别方法
CN115375875A (zh) 基于三维形貌的粮堆凹凸形状特征识别方法、设备及介质
Sun et al. Datum feature extraction and deformation analysis method based on normal vector of point cloud
CN116740303B (zh) 水下目标的离散网格生成方法
Liang et al. Hole-filling algorithm for airborne lidar point cloud data
Zhao et al. Design and development of system for mapping granary on basis of WPF 3D application
CN115797573B (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