CN113610916B - 基于点云数据的不规则物体体积测定方法及*** - Google Patents
基于点云数据的不规则物体体积测定方法及*** Download PDFInfo
- Publication number
- CN113610916B CN113610916B CN202110674138.XA CN202110674138A CN113610916B CN 113610916 B CN113610916 B CN 113610916B CN 202110674138 A CN202110674138 A CN 202110674138A CN 113610916 B CN113610916 B CN 113610916B
- Authority
- CN
- China
- Prior art keywords
- point cloud
- axis
- slice
- curvature
- volume
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000001788 irregular Effects 0.000 title claims abstract description 39
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 27
- 238000005315 distribution function Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 37
- 238000010276 construction Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 238000001914 filtration Methods 0.000 description 8
- 239000002699 waste material Substances 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 6
- 239000002910 solid waste Substances 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011449 brick Substances 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000010440 gypsum Substances 0.000 description 1
- 229910052602 gypsum Inorganic materials 0.000 description 1
- 239000010813 municipal solid waste Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
本发明提供了基于点云数据的不规则物体体积测定方法和***，将被测不规则物体的点云数据的信息投影到空间系Z轴，确定z轴上的概率密度函数，进而确定z轴投影上的曲率函数，利用曲率函数计算z轴投影上各点的曲率；根据各点的曲率、设定的基准间距、基准曲率和相对间距系数确定z轴上相邻切片间距较优长度列表；按照列表对点云进行切片，依次得到z轴方向点云台体；对每一切片横截面进行边界检测并计算各个点云切片的横截面面积；根据各个点云切片的横截面面积以及各个台体的高度确定各台体的体积，将全部台体加和得到被测不规则物体的最终体积。本发明测定过程简单，结果可靠、高效、精度可控。
Description
技术领域
本发明涉及体积计算技术领域，特别涉及一种基于点云数据的针对固体废物等不规则物体的体积计算算法及***。
背景技术
固体废物是指人类在生产生活活动中丢弃的固体或泥状的物质，是环境的一大污染源。在建筑领域，由于近年来城镇化速度的加快，越来越多的建筑物被拆除，并随之产生大量的建筑垃圾，如塑料、木头、石膏板、红砖、混凝土块、瓷砖石材、编织物等，因为其建筑垃圾体量大，类别多，仅仅靠人工分拣难度巨大且经费高昂，为了避免资源浪费与环境污染，国内外已经开始研究智能机器人分选***，而机器人对固体废物的抓握角度，决定了其是否能够稳固的将建筑垃圾投递到对应类别的箱子中，物体的体积、质心、重量等参数，是决定机器人抓握角度的重要因素。
固体建筑垃圾最大的特点就是材料不同且形状各异，无法根据材料或者形状的某一特性去计算其体积，因此有必要设计出一种针对物体不规则特性的体积测定方法。而传统的基于点云的体积计算方式，例如切片法，是按照某一步长连续均匀切片的，这种方法因为没有考虑到点云在z轴的概率密度分布，从而导致在计算过程中损失了一定精度，同样，它的体积计算模块中，按照柱体体积计算方式，没有考虑到上下两个横截面的面积差异，进一步导致点云体积计算精度降低。
发明内容
本发明的目的就是为了弥补现有的点云体积测定方法的缺陷，提出了基于点云数据的不规则物体体积测定方法及***，计算过程简单，计算精度可控，计算结果可靠，提高了对固废的体积计算的准确度。
本发明提供一种基于点云数据的不规则物体体积测定方法，包括：对被测不规则物体进行扫描获取其点云数据，可选地，对获取的点云数据进行滤波去噪预处理；点云滤波方法包含高斯滤波、均匀采样滤波、体素滤波、统计滤波等方法，这里采用体素滤波方法对原始点云进行滤波，因为体素格滤波器可以达到向下采样同时不破坏点云本身几何结构的功能。滤波后得到的点云数据是较为均匀平滑的，用来为计算点云体积做准备。
根据点云数据中的三维坐标的数据，将点云数据的信息投影到空间系Z轴，确定z轴上的概率密度函数，根据得到点云在的z轴投影上的概率密度函数进而确定z轴投影上的曲率函数，利用曲率函数计算z轴投影上各点的曲率；根据各点的曲率、设定的基准间距、基准曲率和相对间距系数确定z轴上相邻切片间距较优长度列表；按照相邻切片间距较优长度列表对点云进行切片，依次得到z轴方向点云台体；
对每一切片横截面进行边界检测并计算各个点云切片的横截面面积；根据各个点云切片的横截面面积以及各个台体的高度确定各台体的体积，将全部台体的体积加和得到被测不规则物体的最终体积。
进一步地，所述z轴上的概率密度函数表示如下：
使用投影方法得到z轴上的分布函数d(z)，z∈[zmin，zmax]，
其中zmin为点云在的z轴投影上的最小坐标值，zmax为点云在的z轴投影上的最大坐标值；
计算概率密度函数表示如下：
其中D为点的总数，满足：
进一步地，所述曲率函数表示如下：
苴中：
Δpk(z)＝p(z+k·Δz)-p(z+(k-1)Δz)，
c(z)为曲率函数，Δz为概率密度函数p(z)的变量z的单位步长符号“×”表示叉积，“·”表示点乘，“||||”表示模长；其中k＝1，2，3，Δpk(z)表示概率密度函数值p(z+(k-1)Δz)与p(z+k·Δz)的差。
进一步地，z轴上相邻切片间距较优长度列表确定方法如下：
根据下式得到z轴方向上的间距函数：
h(z)＝hb+(c(z)-cb)s
其中c(z)为曲率函数，hb为基准间距、cb为基准曲率、s为相对间距系数；
相邻切片间距较优长度列表确定方法如下：
步骤s21：令z＝zmin；
步骤s22：在相邻切片间距较优长度列表l尾部添加z轴坐标值z；
步骤s23：令z＝z+h(z)；
步骤s24：若z＜zmax，返回步骤s22，否则继续；
步骤s25：在列表l尾部添加z轴坐标值zmax，结束；
zmin为点云在的z轴投影上的最小坐标值，
从下到上数第k个切片的z轴坐标表示如下：
其中li为相邻切片间距较优长度列表l的第i个元素；
其中K为切片数量，且满足
HK＝zmax
进一步地，在z轴坐标为HK点截取的各个点云切片的横截面面积的计算方法如下：
将点云切片垂直于Z轴映射到XOY坐标系，得到二维平面图与轮廓点集：
{(xi，yi)，1≤i≤n}，
按照下式计算所述多边形轮廓的面积，作为对应点云切片横截面的面积；
其中，i为点云切片横截面顶点的符号系数，n为多边形的顶点数量，且
xn+1＝x1，yn+1＝y1。
进一步地，按照下式计算从下至上数第i个台体的体积
其中hi为从下至上数第i个台体的高度，Si为从下至上数第i个切片的面积；且：
hi＝Hk-Hk-1
照下式计算所述点云的体积：
其中n为横截面的数量，台体的数量为n一1。
本发明还提供了基于点云数据的不规则物体体积测定***，包括：数据获取模块、切片模块、边界检测模块和体积计算模块；
数据获取模块，用于对被测不规则物体进行扫描获取其点云数据；
切片模块，用于根据点云数据中的三维坐标的数据，将点云数据的信息投影到空间系Z轴，确定z轴上的概率密度函数，根据得到点云在的z轴投影上的概率密度函数进而确定z轴投影上的曲率函数，利用曲率函数计算z轴投影上各点的曲率；根据各点的曲率、设定的基准间距、基准曲率和相对间距系数确定z轴上相邻切片间距较优长度列表；按照相邻切片间距较优长度列表对点云进行切片，依次得到z轴方向点云台体；
边界检测模块，用于对每一切片横截面进行边界检测并计算各个点云切片的横截面面积；
体积计算模块，用于根据各个点云切片的横截面面积以及各个台体的高度确定各台体的体积，将全部台体加和得到被测不规则物体的最终体积。
本发明所取得的有益技术效果：本发明中提出了一种改进的切片法，该方法首先获得不规则物体的点云信息，将点云信息在z轴方向上投影，根据概率密度函数的曲率进行不等间距的连续切片，从而获得对应于点云主体的一系列离散的点云切片。克服了等间距的切片方式会造成形状体积变化大的部分切片偏少、变化不大的部分反倒切了较多的片的问题。另一方面这种根据切片顺序和多边形的区域逐个搜索点云切片的轮廓边界，从而计算出单个切片点云截面的面积，利用台体计算公式获得相邻切片间部分的体积信息，通过这种方法，将计算出来的切片间部分的体积相加，最后得到整个点云的体积信息。整个算法过程简洁，结果可靠，计算方法高效，精度可控，有效地克服了点云等距分割方法的缺点，提高了原始切片方法的计算精度，为不规则物体的体积计算提供了一个很好的解决方案。
附图说明
图1是家庭中常用的沙发的点云数据；
图2是本发明实施例提供的基于点云数据的不规则物体体积测定方法步骤示意图；
图3是本发明实施例中不规则物的z轴方向点云数据的概率密度函数；
图4是本发明实施例中将不规则物从下至上划分为多个台体的示意图。
具体实施方式
本发明提供的基于点云数据的不规则物体体积测定方法，能够应用在建筑垃圾智能分选***场景下，且物体体积计算是智能分拣***中非常重要的步骤，下面将结合所附图例来解释说明本发明所具有的效果，旨在帮助读者能够更清晰的理解本发明的本质。
实施例：基于点云数据的不规则物体体积测定方法，流程示意图如图2所示，具体包括：
步骤s1，将采集的包含大量散列点、孤立点的原始点云数据进行预处理操作。具体实施例中可选地，步骤s1具体包括：采用体素滤波方法去除原始点云中的毛刺和明显离群点，得到均匀平滑的点云数据。本实施例采用的是家庭中常用的沙发的点云数据，如图1所示。
步骤s2，根据滤波后的点云中的三维坐标的数据，将滤波后的点云数据的信息投影到空间系Z轴，得到点在z轴上的分布函数，进一步计算z轴上的概率密度函数。
本实施例中，使用投影方法得到z轴上的分布函数：d(z)，z∈[zmin，zmax]，令点的总数为D，概率密度函数表示如下：
满足：
根据得到点云在的z轴投影上的概率密度函数(如图3所示)，计算出该概率密度函数上的点在z轴方向的曲率变化函数。
按照下式曲率计算公式计算z轴投影上各点的曲率：
其中：
Δpk(z)＝p(z+k·Δz)-p(z+(k-1)Δz)，
c(z)为得到的曲率函数，Δz为概率密度函数p(z)的变量z的单位步长，符号“×”表示叉积，“·”表示点乘，“||||”表示模长。。
设定基准间距、基准曲率与相对间距系数，根据设计出的间距计算公式计算出z轴方向上的间距函数。
设定基准间距hb、基准曲率cb和相对间距系数s，根据下式得到z轴方向上的间距函数：
h(z)＝hb+(c(z)-cb)s
该函数的意义为：若以z轴上z0为某一片切片的高度，则该切片与相邻上方切片的较优间距为h(z0)。
从z轴最低点开始，运行设计出的间距获知算法，得到相邻切片间距较优长度列表。按照如下间距获知算法，得到相邻切片间距较优长度列表l：
步骤s21：令z＝zmin。
步骤s22：在列表l尾部添加z轴坐标值z。
步骤s23：令z＝z+h(z)。
步骤s24：若z＜zmax，返回步骤s22，否则继续。
步骤s25：在列表l尾部添加z轴坐标值zmax，结束。
进一步地，根据得到的相邻切片间距较优长度列表l所述的分割步长进行不等间距分割，得到多个垂直于z轴的横截面。从下到上数第k个切片的z轴坐标表示如下：
其中K为切片数量，且满足
HK＝zmax。
根据得到的相邻切片间距较优长度列表，从z轴最低点开始，以垂直于z轴方向为横截面，对被测物体进行切片，得到多个垂直于z轴的点集。根据步骤s21-s25所述算法可以获得整个点云的切分位置，即z轴坐标集合。
按照步骤s2获得的相邻切片间距较优长度列表对点云进行切片，点云整体被切割为多个台体，如图4所示。
步骤s3，使用基于α-shape边界检测算法，获取点云切片横截面的轮廓边界。
步骤s4，得到切分z轴坐标集合后，还需要体积计算公式中另一个关键的参数，即点云切片横截面面积S。本步骤确定点云切片横截面轮廓的边界点；根据所述的横截面轮廓边界，计算横截面的面积。
具体为：按顺序依次连接点云切片的边界点，取点云切片的轮廓边界，得到多边形，且不同的点云切片轮廓不相同，轮廓边界为任意多边形。
求解任意多边形的面积，步骤包括：
将点云切片垂直于Z轴映射到XOY坐标系，得到二维平面图与轮廓点集：
{(xi，yi)，1≤i≤n}
一般情况下，由于物体获得的点云图像是不规则的，所以不同点云切片的轮廓边界得到的多边形为任意多边形。求解任意多边形的面积，首先需要知道多边形顶点的坐标信息，即点云中三维坐标的X、Y信息，此处将点云的截面垂直于z轴投影到XOY平面。
按照下式计算所述多边形轮廓的面积，作为对应点云切片横截面的面积。
其中，i为点云切片横截面顶点的符号系数，n为多边形的顶点数量。
步骤s5，根据所述横截面，将被测物体从下至上分解为多个台体，对每个台体使用台体体积计算公式获取体积，再全部加和得到被测物体的最终体积。
用所述横截面，根据所述的相邻切片间距较优长度列表l所述的分割步长，将点云整体分割为多个台体，根据所述横截面的面积，计算所述台体的体积，步骤包括：
按照下式计算从下至上数第i个台体的体积
其中hi为从下至上数第i个台体的高度，Si为从下至上数第i个切片的面积。
按照下式计算所述点云的体积：
其中n为横截面的数量，即台体的数量为n-1。
计算出台体体积ΔVi，其中i为从下至上数第i个台体的体积，例如，图中S1到S10为第一个台体，S10到S11为第二个台体。。。，以此类推。
则点云体积为所以台体体积相加结果，计算公式为：
其中n为横截面的数量，台体的数量为n-1。
本发明提供的方法有效地克服了原始点云分割方法的缺点，提高了原始切片方法的计算精度，为建筑垃圾不规则固废物体的体积计算提供了一个很好的解决方案。
与以上实施例提供的基于点云数据的不规则物体体积测定方相对应地，本实施例提供了基于点云数据的不规则物体体积测定***，包括：数据获取模块、切片模块、边界检测模块和体积计算模块；
所述数据获取模块，用于对被测不规则物体进行扫描获取其点云数据；
所述切片模块，用于根据点云数据中的三维坐标的数据，将点云数据的信息投影到空间系Z轴，确定z轴上的概率密度函数，根据得到点云在的z轴投影上的概率密度函数进而确定z轴投影上的曲率函数，利用曲率函数计算z轴投影上各点的曲率；根据各点的曲率、设定的基准间距、基准曲率和相对间距系数确定z轴上相邻切片间距较优长度列表；按照相邻切片间距较优长度列表对点云进行切片，依次得到z轴方向点云台体；
所述边界检测模块，用于对每一切片横截面进行边界检测并计算各个点云切片的横截面面积；
所述体积计算模块，用于根据各个点云切片的横截面面积以及各个台体的高度确定各台体的体积，将全部台体加和得到被测不规则物体的最终体积。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的***，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
本领域内的技术人员应明白，本申请的实施例可提供为方法、***、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。
Claims (6)
1.基于点云数据的不规则物体体积测定方法，其特征在于，包括：
对被测不规则物体进行扫描获取其点云数据；
根据点云数据中的三维坐标数据，将点云数据的信息投影到空间系Z轴，确定z轴上的概率密度函数，进而确定z轴投影上的曲率函数，利用曲率函数计算z轴投影上各点的曲率；根据各点的曲率、设定的基准间距、基准曲率和相对间距系数确定z轴上相邻切片间距较优长度列表；按照相邻切片间距较优长度列表对点云进行切片，依次得到z轴方向点云台体；
对每一切片横截面进行边界检测并计算各个点云切片的横截面面积；根据各个点云切片的横截面面积以及各个台体的高度确定各台体的体积，将全部台体的体积加和得到被测不规则物体的最终体积；
z轴上相邻切片间距较优长度列表确定方法如下：
根据下式得到z轴方向上的间距函数：
其中/>为曲率函数，/>基准间距、/>基准曲率、/>相对间距系数；
相邻切片间距较优长度列表确定方法如下：
步骤s21：令；
步骤s22：在相邻切片间距较优长度列表尾部添加z轴坐标值/>；
步骤s23：令；
步骤s24：若，返回步骤s22，否则继续；
步骤s25：在列表尾部添加z轴坐标值/>，结束；
为点云在的z轴投影上的最小坐标值；
从下到上数第个切片的z轴坐标表示如下：
其中/>相邻切片间距较优长度列表/>的第i个元素；为切片数量，且满足/>。
2.根据权利要求1所述的基于点云数据的不规则物体体积测定方法，其特征在于，所述z轴上的概率密度函数表示如下：
使用投影方法得到z轴上的分布函数，其中/>为点云在的z轴投影上的最小坐标值，/>为点云在的z轴投影上的最大坐标值；
计算概率密度函数表示如下：
其中D为点的总数，满足：/>。
3.根据权利要求1所述的基于点云数据的不规则物体体积测定方法，其特征在于，所述曲率函数表示如下：
其中：，/>为曲率函数，/>为概率密度函数的变量/>的单位步长符号“/>”表示叉积，“/>”表示点乘，“/>”表示模长；k=1，2，3，/>表示概率密度函数值/>与/>的差。
4.根据权利要求1所述的基于点云数据的不规则物体体积测定方法，其特征在于，在z轴坐标为点截取的各个点云切片的横截面面积的计算方法如下：
将点云切片垂直于Z轴映射到XOY坐标系，得到二维平面图与轮廓点集：
，
轮廓为多边形轮廓；
按照下式计算所述多边形轮廓的面积，作为对应点云切片横截面的面积；
其中，/>为点云切片横截面顶点的符号系数，/>为多边形的顶点数量，且/>。
5.根据权利要求1所述的基于点云数据的不规则物体体积测定方法，其特征在于，按照下式计算从下至上数第i个台体的体积
其中/>为从下至上数第i个台体的高度，/>为从下至上数第i个切片的面积；且：/>，/>为第k个切片的z轴坐标；
照下式计算所述点云的体积：
其中/>为横截面的数量，台体的数量为/>。
6.基于点云数据的不规则物体体积测定***，其特征在于，包括：数据获取模块、切片模块、边界检测模块和体积计算模块；
所述数据获取模块，用于对被测不规则物体进行扫描获取其点云数据；
所述切片模块，用于根据点云数据中的三维坐标的数据，将点云数据的信息投影到空间系Z轴，确定z轴上的概率密度函数，根据得到点云在的z轴投影上的概率密度函数进而确定z轴投影上的曲率函数，利用曲率函数计算z轴投影上各点的曲率；根据各点的曲率、设定的基准间距、基准曲率和相对间距系数确定z轴上相邻切片间距较优长度列表；按照相邻切片间距较优长度列表对点云进行切片，依次得到z轴方向点云台体；
z轴上相邻切片间距较优长度列表确定方法如下：
根据下式得到z轴方向上的间距函数：
其中/>为曲率函数，/>基准间距、/>基准曲率、/>相对间距系数；
相邻切片间距较优长度列表确定方法如下：
步骤s21：令；
步骤s22：在相邻切片间距较优长度列表尾部添加z轴坐标值/>；
步骤s23：令；
步骤s24：若，返回步骤s22，否则继续；
步骤s25：在列表尾部添加z轴坐标值/>，结束；
为点云在的z轴投影上的最小坐标值；
从下到上数第个切片的z轴坐标表示如下：
其中/>相邻切片间距较优长度列表/>的第i个元素；为切片数量，且满足/>；
所述边界检测模块，用于对每一切片横截面进行边界检测并计算各个点云切片的横截面面积；
所述体积计算模块，用于根据各个点云切片的横截面面积以及各个台体的高度确定各台体的体积，将全部台体的体积加和得到被测不规则物体的最终体积。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110674138.XA CN113610916B (zh) | 2021-06-17 | 2021-06-17 | 基于点云数据的不规则物体体积测定方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110674138.XA CN113610916B (zh) | 2021-06-17 | 2021-06-17 | 基于点云数据的不规则物体体积测定方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113610916A CN113610916A (zh) | 2021-11-05 |
| CN113610916B true CN113610916B (zh) | 2024-04-12 |
Family
ID=78303565
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110674138.XA Active CN113610916B (zh) | 2021-06-17 | 2021-06-17 | 基于点云数据的不规则物体体积测定方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113610916B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115201778B (zh) * | 2022-09-09 | 2023-02-03 | 广州小鹏自动驾驶科技有限公司 | 不规则障碍物检测方法、车辆及计算机可读存储介质 |
| CN117291845B (zh) * | 2023-11-27 | 2024-03-19 | 成都理工大学 | 一种点云地面滤波方法、***、电子设备及存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011153624A2 (en) * | 2010-06-11 | 2011-12-15 | Ambercore Software Inc. | System and method for manipulating data having spatial coordinates |
| CN109671155A (zh) * | 2018-12-21 | 2019-04-23 | 北京林业大学 | 基于点云数据的表面网格重建方法、***及相关设备 |
| CN109697731A (zh) * | 2018-12-21 | 2019-04-30 | 北京林业大学 | 基于点云数据的植物冠层体积计算方法、***及相关设备 |
| WO2019179168A1 (zh) * | 2018-03-22 | 2019-09-26 | 深圳光峰科技股份有限公司 | 投影畸变校正方法、装置、***及存储介质 |
| WO2020118970A1 (zh) * | 2018-12-10 | 2020-06-18 | 华东理工大学 | 一种复杂表面加工的在机点云检测及补偿方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8396284B2 (en) * | 2007-10-23 | 2013-03-12 | Leica Geosystems Ag | Smart picking in 3D point clouds |
-
2021
- 2021-06-17 CN CN202110674138.XA patent/CN113610916B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011153624A2 (en) * | 2010-06-11 | 2011-12-15 | Ambercore Software Inc. | System and method for manipulating data having spatial coordinates |
| WO2019179168A1 (zh) * | 2018-03-22 | 2019-09-26 | 深圳光峰科技股份有限公司 | 投影畸变校正方法、装置、***及存储介质 |
| WO2020118970A1 (zh) * | 2018-12-10 | 2020-06-18 | 华东理工大学 | 一种复杂表面加工的在机点云检测及补偿方法 |
| CN109671155A (zh) * | 2018-12-21 | 2019-04-23 | 北京林业大学 | 基于点云数据的表面网格重建方法、***及相关设备 |
| CN109697731A (zh) * | 2018-12-21 | 2019-04-30 | 北京林业大学 | 基于点云数据的植物冠层体积计算方法、***及相关设备 |
Non-Patent Citations (2)
| Title |
|---|
| 基于四面体的不规则表面文物体积计算;张小青;朱光;侯妙乐;黄鹤;;测绘通报(10);全文 * |
| 基于机载LiDAR数据的林木冠层投影面积与体积测量;张卫正;董寿银;王国飞;裘正军;吴翔;何勇;;农业机械学报(01);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113610916A (zh) | 2021-11-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113610916B (zh) | 基于点云数据的不规则物体体积测定方法及*** | |
| WO2021000719A1 (zh) | 一种基于三维点云的小曲率薄壁零件机器人加工边界提取方法 | |
| CN108986048B (zh) | 基于线激光扫描三维点云快速复合滤波处理方法 | |
| CN110599506B (zh) | 一种复杂异形曲面机器人三维测量的点云分割方法 | |
| CN107464223B (zh) | 一种基于切片的点云孔洞修补方法 | |
| CN111598780B (zh) | 一种适用于机载LiDAR点云的地形自适应插值滤波方法 | |
| CN103310481A (zh) | 一种基于模糊熵迭代的点云精简方法 | |
| KR101549155B1 (ko) | 라이다 자료를 활용한 구조물의 직선경계 추출방법 | |
| CN112132840A (zh) | 一种车载行道树点云分类与特征信息提取方法 | |
| Fateeva et al. | System analysis of the quality of meshes in HBIM | |
| CN116258857A (zh) | 一种面向室外树木激光点云分割与提取方法 | |
| CN115795632A (zh) | 一种基于标记点云的自动几何孪生方法和*** | |
| CN103934569A (zh) | 一种基于选择性激光烧结的分层切片方法 | |
| CN113763529A (zh) | 一种基于三维扫描的变电站建模方法 | |
| CN107292100B (zh) | 一种岩石结构面分形维数的计算方法 | |
| CN115880354B (zh) | 一种基于点云自适应切片的计算树冠体积的方法 | |
| CN115854883B (zh) | 一种适用于长钢筋笼的自动化检测方法及装置 | |
| CN107545596B (zh) | 一种点云模型最优切割平面的提取方法 | |
| Kampel et al. | Classification of archaeological fragments using profile primitives | |
| CN116452604A (zh) | 一种复杂变电站场景分割方法、设备及存储介质 | |
| CN116805413A (zh) | 一种基于三维点云的油茶苗期表型特征自动计算方法 | |
| CN113283495B (zh) | 一种集料颗粒分档方法及其装置 | |
| CN112634458B (zh) | 一种基于特征显著性评价的口腔模型点云精简方法 | |
| CN112365534B (zh) | 一种基于单目相机三维重建的大型煤堆体积测量方法 | |
| CN106056124A (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 |