基于kriging插值的大型料堆體积计算方法.docVIP

PAGE PAGE 3 基于kriging插值的大型料堆体积计算方法 陈春，杨为民，张朋新 （东北大学 信息科学与工程学院 ，辽宁 沈阳 110819） 摘要：本文提出一种定量测量物料堆体积的方法，将料堆剖分成若干个小矩形柱体，通过计算每个矩形柱体的体积，从而得到料堆的体积。该方法包括： 原始数据点的整理、 边界点的提取与确定、投影区域划分、插值计算Z值和计算体积5个步骤。并在此基础上开发一套测量系统，通过实验验证了算法的可行性和正确性。 关键词：克里金插值；料堆体积测量；计算机辅助测量 0．引言 如何准确测量大型物料堆的体积一直是工程领域的重点。采用手工测量方法工作繁重，时间长，结果误差大。 随着计算机视觉技术的迅速发展，人们开始研究将此技术运用于测量一个物体的体积，在进行三维重建时采用的方法主要建立在三角网法基础上[1][2]。而常用的一般三角网法未考虑地性线的骨架作用，精度不高。故在此，本文提出一种新的定量测量物料堆体积测量方法，通过黑白网格在料堆上人工设置标识点，并在计算机中对每张照片上料堆的标志点进行标识后，综合各标志点在不同照片上的位置，进行立体恢复，从而计算出各标志点的真实空间坐标置X、Y 、Z 值； 利用所得标志点的坐标，计算出料堆的体积。 1.算法设计1.1算法的基本思想 根据计算得出各标志点的真实空间坐标值，求出其投影XY平面，并将该平面划分为n 个小矩形，利用克里金插值法分别对其进行矩形顶点（X，Y）插值求取Z值，计算每一个小矩形柱的体积，则整个料堆的体积等于所有区域体积之和： 1.2算法实现的步骤 (1)整理原始数据点。①删除距离过近的点 ； ② 归一化。将所有数据点的XY Z坐标变换到区域，，内。③设置零层点，当原始数据点的Z 坐标小于“ 零层点最小误差” 时，认为该点是零层点。 (2)对料堆在XY平面的投影进行边界提取。因为料堆本身的特点，Z=0的点不一定就是边界。故应该先寻找边界。算法步骤如下： ①求取中心坐标（,）。，。②以中心坐标（，）为原点，以中心坐标和X取最大值点连线为初始搜索径向量（向量夹角可变，变化范围（0~360）），旋转向量，计算点到远点的距离，在向量上距离最短的点即为边界点，记录。③判断搜索的数据是否是边界点。以搜索到原始数据点为中心，用②所示方法搜索得到其附近三维点，计算该点和附近三维点的空间距离，选择距离较小的三个点。计算其偏导数，不为0 ，则该点为边界点。 (3)划分插值点，将投影划分为N个小矩形。①分别对获取的离散边界点采用拉格然日插值法绘制上下边界曲线，以便于获得要插值的二维平面点的变化范围。②将料堆的投影平面划分为若干个小区间，如图所示。将X轴分成M等分，在y轴上以d为步长，确定每一个要插值的二维点。由于料堆边界的不规则性，如图可见，投影并不是被划分为很规则的矩形，在靠近边界曲线处，会出现梯形或者三角形，在计算体积是还会再具体讨论。 图一 投影划分示意图 (4)区域曲面拟合 为保证每一个插值点z值都足够准确，我们选用克里金插值法来获取z值。克里金插值法在插值过程中根据某种优化准则函数来动态的决定变量的数值，内插函数时刻处于最佳状态。克里金方法的插值公式为 (1) 式中：处的估计值 ； 处的观测值；克里金权重系数； n观测点个数。 在该算法中， 将研究对象称为区域化变量，区域化变量在整个研究区域内满足二阶平稳假设： (1)的数学期望存在且等于常数：。 (2)的协方差存在且至与两点之间的相对位置有关。 或满足本征假设 ： (3) 。 (4) 增量的方差存在且平稳 ：。 依据无偏性要求 ： 经推导可得 ： (2) 在无偏条件下使估计方差达到最小，即 ： (3) 其中：为拉格朗日乘子 。 可得求解权系数的方程组 ： (4) 求出诸权系数后，如此，可求出未采样点处的属性值。 （5）计算体积 根据插值的划分区间投影形状，矩形柱体积的计算分为以下几种情况： 投影为矩形，如示意图二: (5) 投影为三角形，如示意图三: (6) 投影为梯形，如示意图四: (7) 图二 示意图 图三 示意图 图Ⅳ 示意图 如此，计算出每一个小区间的体积相加之后，即可计算出料堆的体积。可得到料堆体积公式 (8) 2 试验 为了验