本发明专利技术属于激光雷达点云处理和矿山监测领域，具体技术方案为：一种基于激光雷达点云的矿山开采量计算方法，首先，基于激光雷达系统对矿山进行定期三维地形测绘，获取矿山原始点云数据；然后对激光雷达点云数据进行去噪、配准等预处理；其次检测点云局部空洞面并利用邻域信息内插三维数据点进行修补，前后两期点云数据形成某时段不规则疏密多变矿体的三维模型；最后，设置离散化步长对矿体三维模型进行块体离散化并提取块体体积，各单元块体体积累加得到整个矿山变化体积，即该时段的矿山开采量，本发明专利技术能够对矿山开挖方量进行快速估算，该方法具备测量高效安全、灵活便捷、即时性高等特点，可应用于三维变化监测、工程测管、智慧矿山建设等领域。慧矿山建设等领域。慧矿山建设等领域。
【技术实现步骤摘要】
一种基于激光雷达点云的矿山开采量计算方法
[0001]本专利技术属于激光雷达点云处理和矿山监测领域，具体涉及一种基于激光雷达点云的矿山开采量计算方法。
技术介绍
[0002]矿产资源是人类赖以生存和经济发展的重要物资基础，开采量动态监测可为其科学高效的生产管理提供依据，对矿产资源的可持续发展与智慧矿山建设意义重大。矿山开采尤其是露天开采覆盖范围大且地形复杂多变，导致开采量测算难度大、效率低、成本高。如何快速准确地获取矿山地形变化并计算矿山开采量成为了一个亟待解决的问题。
[0003]已有矿山开采量计算方法主要包括两种，一种是基于人工测绘的开采量计算，此类方法通过地测人员手动测量获取的开采区域工作面三维坐标生成地形图，勘测效率低且易受人为因素影响，导致开采量计算速度与精度不稳定；另一种是基于无人机可见光影像的开采量计算，此类方法通过获取时序无人机影像后重建三维模型，采集效率高于人工测绘方式，但由于露天矿地形起伏较大且边界持续移动，三维模型高程的真实性难以得到保障，依然无法在精度上满足矿山开采量计算的工程需求。
[0004]激光雷达系统是一种集激光，全球定位系统(GPS)和惯性导航系统(INS)三种技术于一身的系统，可高效、精确、无损地获得地表三维信息。将激光雷达系统引入矿山开采量动态监测中，可通过定期获取矿山的精细三维点云数据，构建某时期不规则疏密多变矿体的三维模型，快速准确获取矿山开挖过程中的真实地形地貌变化，为开采量计算提供可靠的数据支撑。
技术实现思路
[0005]本专利技术为了满足矿山开采量计算快速准确的需求，提供一种基于激光雷达点云的矿山开采量计算方法，通过对矿山进行定期三维地形测绘更新并基于精细三维点云数据计算矿山变化体积，从而实现矿山开采过程中的动态开采量计算。
[0006]为实现上述目的，本专利技术所采用的技术方案为：一种基于激光雷达点云的矿山开采量计算方法，具体步骤如下：一、激光雷达点云获取，对矿山进行矿区勘探后确定监测范围，根据实际监测区域大小与高差选取激光雷达系统并制定数据采集方案，完成外业数据采集，经计算机软件处理后解算得到激光雷达点云数据。
[0007]进一步是，步骤一包括以下步骤：S11、激光雷达点云数据采集：对矿山的地形起伏、开采边界、典型地物进行勘探后确定矿山监测区域，根据监测区域的大小和高差选取续航时间与探测距离符合要求的激光雷达系统，结合系统具体的性能参数确定测量距离、扫描速度与扫描角度，对矿山实施外业数据采集。
[0008]S12、数据解算：基于原始惯导数据、原始激光测距数据（由激光扫描仪扫描采集得
到）、GPS数据通过计算机联合解算得到激光雷达点云数据。
[0009]二、点云预处理：对原始激光雷达点云数据进行去噪、配准等预处理。
[0010]进一步是，步骤二包括以下步骤：S21、点云去噪：对原始激光雷达点云采用稀疏离群点移除（Statistical Outlier Removal，SOR）方法进行去噪，该方法计算每个点到其所有临近点的平均距离，假设得到的结果是一个高斯分布，平均距离在标准范围（由全局距离平均值和方差定义）之外的点，被定义为离群点并从数据集中去除掉。
[0011]S22、点云配准：对开采前后两期点云数据采用迭代最近点（Iterative Closest Point, ICP）自动配准算法，ICP算法通过找到两组点云集合中距离最近的点对，根据估计的变换关系（旋转与平移）来计算距离最近点对经过变换之后的误差，不断的迭代直至使得设定的目标函数达到最小值，得到最优平移矩阵和旋转矩阵，以确定最终的变换关系，目标函数如下：式中，R表示旋转矩阵，T表示平移矩阵，k为待配准点云的点数, ɑ
i
与b
i
分别对应参考点云与待配准点云。
[0012]最终将待配准的点云与参考点云的数据对进行一一对应，即完成两期点云数据配准。
[0013]三、点云空洞面修补：在三维空间中根据点云数据构建初始三角网格，设定最小空洞面积阈值P，自动检测面积大于阈值P的点云空洞网格。利用邻域信息对空洞网格进行三维数据点内插，点云空洞面修补完成后得到矿山精细三维点云数据，前后两期点云数据形成不规则疏密多变矿体的三维模型。
[0014]S31、空洞面检测：矿山由于其本身反射率有限且存在部分水域，数据采集中往往存在三维点云局部空洞的问题，尤其是在露天矿的矿底部与侧面，开采量计算有赖于完整精细的点云数据，需要对空洞面进行检测与修补。基于原始三维点云数据构建初始三角网络，遍历计算所有三角网格的面积，设定最小空洞面积阈值P，自动检测面积大于阈值P的网格为需要修补的空洞面。
[0015]其中，初始三角网络构建的基本步骤是：1）、以点云数据左下角任一点为起始点；2）、找出与起始点最近的数据点相互连接形成三角形的一条边作为基线；3）、在基线右侧进行扩展，搜索与基线距离最近的点构成起始三角形的第三点；4）、基线的两个端点与第三点相连，成为新的基线；5）、继续搜索与新的基线距离最近的第三点构成新的三角形，直至所有的点都被包含在三角网格中，停止搜索，初始三角网络构建完毕。
[0016]S32、空洞面修补：对需要修补的三角网格的每条边取中点，分别为点，连接三点后将空洞网格分割为四个新的三角形，分别取四个三角形的重心作为新的离散点内插至空洞网格中，重心坐标为。
[0017]再次计算四个新的三角形的面积，仍大于阈值P的三角网格被划分为空洞网格，重复进行三角形分割与离散点内插，直到整体三角网格中的每一个网格面积均小于P，停止内插，得到精细矿山三维点云数据。
[0018]空洞面修补后，前后两期点云数据形成不规则疏密多变矿体的三维模型。
[0019]四、开采量计算：根据矿体三维模型的面积大小与点密度确定离散化步长，对其进行单元块体离散化，获取各块体的网格面积与块体相对高度，计算各单元块体体积，将整个矿体所有单元块体的体积叠加获取整个矿山变化体积，即该时段的矿山开采量。
[0020]进一步是，步骤四包括以下步骤：S41、单元块体离散化：结合点云密度与测区大小确定离散化步长，将不规则疏密多变矿体三维模型划分成个单元块体。
[0021]S42、矿山开采量计算：计算各单元块体的顶面积S
t
和底面积S
d
，取顶面积与底面积的平均值作为网格面积S
m
，计算公式如下：其中，i为各单元块体的序号，；三维模型各单元块体对应的顶面与底面点云的高程值的平均高度作为单元块体高程值的最大值Z
max
与最小值Z
min
，对每个单元块体对应的顶面高程值与底面高程值做差得到块体相对高度H，计算公式如下：每个单元块体利用网格面积与块体相对高度相乘得到离散化块体的体积，经计算并累加求和得到矿体变化体积，即为该时段矿山的开采量，计算公式如下：。
[0022]与现有技术相比，本专利技术的有益效果是：一、本专利技术将激光雷达探测手段引入到矿山开本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种基于激光雷达点云的矿山开采量计算方法，其特征在于，具体步骤如下：一、激光雷达点云获取：对矿山进行矿区勘探后确定监测范围，根据实际监测区域大小与高差选取激光雷达系统并制定数据采集方案，完成外业数据采集，经计算机软件处理后解算得到激光雷达点云数据；二、点云预处理：对原始激光雷达点云数据进行去噪、配准预处理；三、点云空洞面修补：在三维空间中根据点云数据构建初始三角网格，设定最小空洞面积阈值P，自动检测面积大于阈值P的点云空洞网格，利用邻域信息对空洞网格进行三维数据点内插，点云空洞面修补完成后得到矿山精细三维点云数据，前后两期点云数据形成不规则疏密多变矿体的三维模型；四、开采量计算：根据矿体三维模型的面积大小与点密度确定离散化步长，对其进行单元块体离散化，获取各块体的网格面积与块体相对高度，计算各单元块体的体积，将整个矿体所有单元块体的体积叠加获取整个矿山变化体积，即该时段的矿山开采量。2.根据权利要求1所述的一种基于激光雷达点云的矿山开采量计算方法，其特征在于，在步骤一中，包括以下两个步骤：S11、激光雷达点云数据采集：对矿山的地形起伏、开采边界、典型地物进行勘探后确定矿山监测区域，根据监测区域的大小和高差选取续航时间与探测距离符合要求的激光雷达系统，结合系统具体的性能参数确定测量距离、扫描速度与扫描角度，对矿山实施外业数据采集；S12、数据解算：基于原始惯导数据、原始激光测距数据、GPS数据通过计算机联合解算得到激光雷达点云数据。3.根据权利要求2所述的一种基于激光雷达点云的矿山开采量计算方法，其特征在于，所述原始激光测距数据由激光扫描仪扫描采集得到。4.根据权利要求1所述的一种基于激光雷达点云的矿山开采量计算方法，其特征在于，在步骤二中，包括以下两个步骤：S21、点云去噪：对原始激光雷达点云采用稀疏离群点移除方法进行去噪，该方法计算每个点到其所有临近点的平均距离，假设得到的结果是一个高斯分布，平均距离在标准范围之外的点，被定义为离群点并从数据集中去除掉；S22、点云配准：对开采前后两期点云数据采用迭代最近点自动配准算法，ICP算法通过找到两组点云集合中距离最近的点对，根据估计的变换关系来计算距离最近点对经过变换之后的误差，不断的迭代直至使得设定的目标函数达到最小值，得到最优平移矩阵和旋转矩阵，以确定最终的变换关系，目标函数如下：式中，R表示旋转矩阵，T表示平移矩阵，k为待配准点云的点数，
ɑ
i
为参考点云；b
i
【专利技术属性】
技术研发人员：李晓方，杨风暴，王强，高敏，高佳锋，张付岗，马泽亮，赵永鑫，胡向荣，安开敏，刘一杰，曹艳君，
申请(专利权)人：山西阳光三极科技股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。