一种边坡点云曲面重建及变形信息提取方法技术

本发明专利技术提出一种边坡点云曲面重建及变形信息提取方法，属于变形监测领域，本发明专利技术采取了改进的RBF算法，通过对经过点云预处理之后的边坡点云数据进行了曲面建模，提高了RBF算法建模效率和精度；针对边坡点云信息提取问题，提出Cyclone结合CAD提取DEM与等高线法、三维质检软件法两种方法的可行性，并将三维质检软件法提取的变形数据与传统变形监测数据进行了对比，验证了方法的正确性。

【技术实现步骤摘要】

本专利技术属于变形监测领域，具体涉及一种边坡点云曲面重建及变形信息提取方法。
技术介绍
20世纪90年代中期出现然后极为迅速地发展起来并且被人熟知的三维激光扫描技术，据有其他传统边坡变形监测方法无法比拟的优势；三维激光扫描技术获取高密度高精度的三维点云数据高效快速，对监测结果进行研究掌握其变形规律有很大帮助。在三维激光扫描应用于边坡变形监测的相关技术之中，散乱的边坡点云曲面重构技术和边坡点云变形信息提取方法选取是关键。1971年，Hardy在研究曲面的二元函数中使用了MQ基函数来构造和拟合曲面，在他研究之后许多学者围绕着他的方法又研究出大量的成果，为后续被人广泛应用的的径向基函数的出现奠定了大量的理论研究基础；屈永华建立了一种适用于地质曲面重建的径向基函数(RBF)神经网络模型，并给出了具体算法和一般曲面插值方法相比，该模型对原始数据的分布和边界条件无特别的要求，因而适合于遭到破坏或勘探数据较少的地质曲面的重建，但这种进行的点云曲面重建的算法只适用于数据量小的散乱点云数据，曲面重建的精度比较低。2010年赵小平等对边坡点云数据进行了处理，再通过扫描数据处理软件获得DEM数据边坡形态，他的研究中验证了点云预处理之后的精度，为边坡变形监测与灾害预报提供基础数据。2015年刘昌军等提出基于八叉树结构的点云与点云的精细直接比较算法，实现变形区域多个点云数据的自动比较，其研究运用的精细自动比较算法随点云的数据量增加，效率降低，不能快速高效的提取复杂大量数据的边坡点云变形信息。
技术实现思路
针对现有技术的不足，本专利技术提出一种边坡点云曲面重建及变形信息提取方法，以达到提高提取变形信息效率和正确性目的。一种边坡点云曲面重建及变形信息提取方法，包括以下步骤：步骤1、通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标；步骤2、根据实际需求间隔一段时间，再一次通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标；步骤3、对三维激光扫描仪获取的边坡三维点云数据进行去噪和拼接预处理；步骤4、基于最小包围盒的RBF点云曲面重建，构建边坡点云曲面；步骤5、对三维激光扫描仪和全站仪测得的两次点云坐标进行比较，并根据预处理过的两期三维点云数据生成DEM模型，将两个DEM模型进行比较，获得边坡变形信息。步骤1所述的通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标，具体为：三维激光扫描仪扫描时每站需布设有4～6个反射片。步骤4所述的基于最小包围盒的RBF点云曲面重建，构建边坡点云曲面，具体步骤如下：步骤4-1、输入预处理后的点云数据；步骤4-2、用坐标系旋转确定最小包围盒，具体如下：设三个坐标轴的初始转角均为0，按设定的角度步长，用三重循环分别使三个坐标由0°至360°进行旋转；在每一个旋转位置，确定包围盒的边界数据和体积，求出体积最小包围盒及相对应的转角；步骤4-3、沿坐标系x轴方向分解包围盒，若点云数据集的点数大于设定的包围盒内最多点数，并且包围盒的边长度大于设定的最小盒子边长，则获取点云数据集所有点的重心；步骤4-4、采用重心x坐标所在平面对点云数据集和包围盒进行切分，若切分后包围盒内的数据集小于设定的包围盒内最多点数，或者在x方向的包围盒边长小于等于设定的最小盒子边长，则执行步骤4-5，否则返回执行步骤4-4；步骤4-5、对切分获得的相邻包围盒及相邻的数据集求并集，若某个包围盒内的数据大于设定的包围盒内最多点数，则在y轴和z轴方向上分割包围盒，直到包围盒内的数据均小于等于设定的包围盒内最多点数，且所有包围盒的并集包含该包围盒；步骤4-6、采用径向基函数插值方法对分割后的包围盒的数据进行插值加密，根据插值后的数据构建边坡点云曲面。步骤5所述的对三维激光扫描仪和全站仪测得的两次点云坐标进行比较，并根据预处理过的两期三维点云数据生成DEM模型，将两个DEM模型进行比较，获得边坡变形信息；具体步骤如下：步骤5-1、判断三维激光扫描仪得到的反射片的坐标值与全站仪得到的反射片的坐标值的差值是否在设定范围内，若是，则三维激光扫描仪获得的数据能满足边坡变形监测的要求，否则三维激光扫描仪重新采集数据；步骤5-2、将预处理后的三维激光扫描仪两次采集的点云数据导入到Cyclone软件中；步骤5-3、利用Cyclone软件的mesh模块生成两个点云的DEM模型；步骤5-4、将两个点云的DEM模型进行插值运算，找到高程相同的点，连线生成高密度等高线图，并将其另存为dwf格式的文件；步骤5-5、用CAD打开dwf格式的文件，并对变形信息进行提取和标注，得到用于分析变形的高密度等高线文件；步骤5-6、以第一次扫描的DEM模型当做参考对象，第二次扫描的DEM模型当成结果对象，使用GeomagicQualify软件对两个DEM模型进行比较，获得边坡变形信息。本专利技术优点：本专利针对边坡点云数据模型重建问题，采取了改进的RBF算法，通过对经过点云预处理之后的边坡点云数据进行了曲面建模，提高了RBF算法建模效率和精度；针对边坡点云信息提取问题，提出Cyclone结合CAD提取DEM与等高线法、三维质检软件法两种方法的可行性，并将三维质检软件法提取的变形数据与传统变形监测数据进行了对比，验证了方法的正确性。附图说明图1为本专利技术一种实施例的边坡点云曲面重建及变形信息提取方法流程图；图2为本专利技术一种实施例的站位和控制点布局简略示意图；图3为本专利技术一种实施例的最小包围盒的3种形式示意图，其中，图(a)为AABB包围盒示意图，图(b)为0BB包围盒示意图，图(c)为凸包围盒示意图；图4为本专利技术一种实施例的确定最小包围盒的流程图；图5为本专利技术一种实施例的点云边坡曲面与点云对比图，其中，图(a)是未被改进算法构建的点云边坡曲面与点云对比图，图(b)是改进算法构建的点云边坡曲面与点云对比图；图6为本专利技术一种实施例的等高线文件示意图，其中，图(a)是第一次扫描的等高线，图(b)是第二次扫描的等高线；图7为本专利技术一种实施例的边坡变形示意图。具体实施方式下面结合附图对本专利技术一种实施例做进一步说明。本专利技术实施例中，边坡点云曲面重建及变形信息提取方法，如图1所示，包括以下步骤：步骤1、通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标；本专利技术实施例中，根据地形条件初步定出设站位置；本专利技术实施例中扫描对象为门头沟一个尾矿的边坡，环境较为复杂；尾矿底部为相对稳固的地面，可以将三维激光扫描仪架设在上面；为保证数据的完整性，同时为了避免数据重叠，总共设站6站，每一测站保证周边能观测到布设有4～6个反射体，并且在边坡上固定6个反射体(这六个反射体中心点坐标不进行点云拼接)，其位置必须稳固，不易移动；并用全站仪测出每一个测站每一个反射片的中心坐标。步骤2、一个月后对尾矿边坡进行第二次扫描，再一次通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标；步骤3、对三维激光扫描仪获取的边坡三维点云数据进行去噪和拼接预处理；步骤3-1、由于尾矿边坡点云既包含坡上杂草、小灌木、电力设施和岩石早点，又包含三维激光扫描系统自身存在着一些系统误差，所以点云情况较复杂，本专利技术实施例中提出了基本文档来自技高网...

【技术保护点】
一种边坡点云曲面重建及变形信息提取方法，其特征在于，包括以下步骤：步骤1、通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标；步骤2、根据实际需求间隔一段时间，再一次通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标；步骤3、对三维激光扫描仪获取的边坡三维点云数据进行去噪和拼接预处理；步骤4、基于最小包围盒的RBF点云曲面重建，构建边坡点云曲面；步骤5、对三维激光扫描仪和全站仪测得的两次点云坐标进行比较，并根据预处理过的两期三维点云数据生成DEM模型，将两个DEM模型进行比较，获得边坡变形信息。

【技术特征摘要】
1.一种边坡点云曲面重建及变形信息提取方法，其特征在于，包括以下步骤：步骤1、通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标；步骤2、根据实际需求间隔一段时间，再一次通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标；步骤3、对三维激光扫描仪获取的边坡三维点云数据进行去噪和拼接预处理；步骤4、基于最小包围盒的RBF点云曲面重建，构建边坡点云曲面；步骤5、对三维激光扫描仪和全站仪测得的两次点云坐标进行比较，并根据预处理过的两期三维点云数据生成DEM模型，将两个DEM模型进行比较，获得边坡变形信息。2.根据权利要求1所述的边坡点云曲面重建及变形信息提取方法，其特征在于，步骤1所述的通过三维激光扫描仪获取边坡三维点云数据，并通过全站仪测出每个反射片的中心坐标，具体为：三维激光扫描仪扫描时每站需布设有4～6个反射片。3.根据权利要求1所述的边坡点云曲面重建及变形信息提取方法，其特征在于，步骤4所述的基于最小包围盒的RBF点云曲面重建，构建边坡点云曲面，具体步骤如下：步骤4-1、输入预处理后的点云数据；步骤4-2、用坐标系旋转确定最小包围盒，具体如下：设三个坐标轴的初始转角均为0，按设定的角度步长，用三重循环分别使三个坐标由0°至360°进行旋转；在每一个旋转位置，确定包围盒的边界数据和体积，求出体积最小包围盒及相对应的转角；步骤4-3、沿坐标系x轴方向分解包围盒，若点云数据集的点数大于设定的包围盒内最多点数，并且包围盒的边长度大于设定的最小盒子边长，则获取点云数据集所有点的重心；步骤4-4、采用重心x坐标所在平面对点云数据集和包围盒进行切分，若切分后包围盒内的数据...

【专利技术属性】
技术研发人员：杨帆，邵阳，常俊飞，李龙飞，
申请(专利权)人：辽宁工程技术大学，
类型：发明
国别省市：辽宁;21