一种基于k均值聚类分析的图像落石体积测定方法
【技术领域】
[0001] 本发明涉及数字图像处理技术领域，具体来说，涉及一种基于K均值聚类分析的 图像落石体积测定方法。
【背景技术】
[0002] 山体滑坡预警系统目前都采用监测仪器进行监测，传统的用于山体滑坡监测的方 法大概有如下几种：常规大地测量法；液体静力水准测量、重力测量法；地下水位监测法； 电测法、地下钻孔倾斜等。尽管这些方法都在山体滑坡预警方面起到了积极作用，但它们具 有一些弊端，如受地形影响，不能连续观测，自动化程度不高，人力投入过大，数据不能实时 处理等。
[0003] 采用基于数字图像处理技术的监测方法则可以消除上述部分弊端，如无接触式测 量，智能化程度高，成本低，可以实时预警等。在以落石作为预警信息的视频监测系统中，监 测落石有无的同时，常常希望给出落石的具体测量参数，如体积。不同数量和不同大小的落 石表征了不同程度的预警，甚至极小的落石可能只是虚警。因此，监测落石的同时，给出相 应的体积测量参数意义及其重要。
[0004] 基于图像处理算法的体积计算技术都是使用多个相机从不同角度拍摄，增大了成 本，且实现困难；而基于单相机的技术则大部分是预先建立待测目标参数的模型库，然后通 过不同算法进行匹配和拟合，不适用于落石的监测。
[0005] 针对相关技术中的问题，目前尚未提出有效的解决方案。
【发明内容】
[0006] 针对相关技术中的问题，本发明提出一种基于K均值聚类分析的图像落石体积测 定方法，其基于落石在运动过程中随机滚动的考虑，利用不同图像序列的面积和名义高度 参数建立落石目标不同姿态下体积数据的集合，并采用K均值聚类算法对集合进行分簇， 数量最大的簇即为真实体积的随机分布，由此估算出落石的体积。
[0007] 本发明的技术方案是这样实现的：
[0008] 一种基于K均值聚类分析的图像落石体积测定方法，包括：
[0009] 通过预先搭建的视频采集系统采集落石视频图像，并对采集的落石视频图像的每 帧图像进行落石检测，确定每帧检测后的落石二值图像；
[0010] 根据确定出的每帧落石二值图像，确定所述每帧落石二值图像所对应的面积和名 义高度；
[0011] 根据所述面积和名义高度，计算落石目标的体积，并根据K均值聚类分析算法，对 落石目标的体积集合进行聚类；
[0012] 对聚类后的若干簇的元数数据进行统计，并根据统计结果，确定出数目最多的簇， 并根据该簇，估算落石目标的体积，其中，估算的方法为取簇中所有元数的均值。
[0013] 进一步的，上述基于K均值聚类分析的图像落石体积测定方法还包括：
[0014] 在通过预先搭建的视频采集系统采集落石视频图像之前，对所述视频采集系统进 行相机标定。
[0015] 进一步的，在通过预先搭建的视频采集系统采集落石视频图像之前，对所述视频 采集系统进行相机标定包括：
[0016] 将预先配置的两标志物置于山体上，促使两标志物连线位于视场水平方向，并测 定两标志物的实际距离；
[0017] 更换两标志物的实际位置，促使两标志物连线位于视场垂直方向，并测定两标志 物的实际距离；
[0018] 根据上述水平方向上的实际距离和垂直方向上得实际距离，确定视频采集系统 中，每像素所表示的实际长度，并根据该实际长度对视频采集系统进行相机标定。
[0019] 进一步的，所述落石检测的方法包括基于背景差分法和目标特征检测算法。
[0020] 进一步的，根据确定出的每帧落石二值图像，确定所述每帧落石二值图像所对应 的面积包括：
[0021] 对单个落石目标进行像素的遍历，其中，像素的数据即为当前帧落石的面积。
[0022] 进一步的，根据确定出的每帧落石二值图像，确定所述每帧落石二值图像所对应 的名义高度包括：
[0023] 根据确定出的每帧落石二值图像，计算落石目标图像的惯性矩；
[0024] 根据所述落石目标图像的惯性矩，确定落石目标图像的长轴值和短轴值；
[0025] 根据所述长轴值和所述短轴值，确定所述长轴值和所述短轴值的均值，并将该均 值作为所述名义高度。
[0026] 本发明的有益效果：本发明以数字图像处理技术和模式识别技术为基础，采用单 相机、无接触式测量技术，利用不同图像序列的面积和名义高度参数建立落石目标不同姿 态下体积数据的集合，并采用K均值聚类算法对集合进行分簇最终估算落石的体积。该算 法无需提前建立待测目标参数的模型库，图像处理算法简单高效。
【具体实施方式】
[0027] 下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通 技术人员所获得的所有其他实施例，都属于本发明保护的范围。
[0028] 根据本发明的实施例，提供了一种基于K均值聚类分析的图像落石体积测定方 法，包括通过预先搭建的视频采集系统采集落石视频图像，并对采集的落石视频图像的每 帧图像进行落石检测，确定每帧检测后的落石二值图像；根据确定出的每帧落石二值图像， 确定所述每帧落石二值图像所对应的面积和名义高度；根据所述面积和名义高度，计算落 石目标的体积，并根据K均值聚类分析算法，对落石目标的体积集合进行聚类；对聚类后的 若干簇的元数数据进行统计，并根据统计结果，确定出数目最多的簇，并根据该簇，估算落 石目标的体积，其中，估算的方法为取簇中所有元数的均值。
[0029] 进一步的，上述基于K均值聚类分析的图像落石体积测定方法还包括：在通过预 先搭建的视频采集系统采集落石视频图像之前，对所述视频采集系统进行相机标定。
[0030] 进一步的，在通过预先搭建的视频采集系统采集落石视频图像之前，对所述视频 采集系统进行相机标定包括：将预先配置的两标志物置于山体上，促使两标志物连线位于 视场水平方向，并测定两标志物的实际距离；更换两标志物的实际位置，促使两标志物连线 位于视场垂直方向，并测定两标志物的实际距离；根据上述水平方向上的实际距离和垂直 方向上得实际距离，确定视频采集系统中，每像素所表示的实际长度，并根据该实际长度对 视频米集系统进行相机标定。
[0031] 进一步的，所述落石检测的方法包括基于背景差分法和目标特征检测算法。
[0032] 进一步的，根据确定出的每帧落石二值图像，确定所述每帧落石二值图像所对应 的面积包括：对单个落石目标进行像素的遍历，其中，像素的数据即为当前帧落石的面积。
[0033] 进一步的，根据确定出的每帧落石二值图像，确定所述每帧落石二值图像所对应 的名义高度包括：根据确定出的每帧落石二值图像，计算落石目标图像的惯性矩；根据所 述落石目标图像的惯性矩，确定落石目标图像的长轴值和短轴值；根据所述长轴值和所述 短轴值，确定所述长轴值和所述短轴值的均值，并将该均值作为所述名义高度。
[0034] 为了方便理解本发明的上述技术方案，以下通过具体工作原理和工作方式，对本 发明的上述技术方案进行详细说明。
[0035] 第一步相机标定