基于椭球立体表征的裂缝三维检测方法技术

基于椭球立体表征的裂缝三维检测方法。本发明专利技术公开了基于视频流的路面裂缝动态检测方法，包括以下步骤：步骤1，读取路面三维图像数据矩阵；步骤2，对路面三维图像数据矩阵进行预处理，得到预处理后三维图像数据；步骤3，对预处理后的三维图像数据矩阵进行滤波处理，得到去噪后的三维图像数据矩阵；步骤4：提取去噪后的三维图像数据矩阵的每行数据，提取出完整的路面裂缝；步骤5，将路面裂缝划分为多个裂缝区域，求取每个附加深度的椭圆模型的表征参数；步骤6，根据表征参数对所有附加深度的椭圆模型进行聚类，得到聚类后的裂缝；步骤7，计算聚类后的裂缝的裂缝特征值，将裂缝分为网状裂缝和线性裂缝。本发明专利技术降低了计算量，复杂度，对细节信息丢失、虚假深度裂缝等原因造成的裂缝段，能够准确还原为完整裂缝。

Three dimensional crack detection method based on ellipsoid stereoscopic representation

Three dimensional crack detection method based on ellipsoid stereoscopic representation. The invention discloses a dynamic pavement crack detection method based on video stream, which comprises the following steps: 1, read the road 3D image data matrix; step 2 on the road, the 3D image data matrix preprocessing preprocessing after three-dimensional image data; step 3, the 3D image data matrix after pretreatment filtering go, get the 3D image data matrix after noise; step 4: each row of data to extract the 3D image data matrix after noise, extract the pavement cracks complete; step 5, the pavement crack is divided into a plurality of crack area, characterization parameters of solving elliptic model of each additional depth; step 6. According to the characterization parameters for clustering elliptical model all the additional depth of the fractures were obtained after clustering; step 7, the fracture characteristics calculation of crack after clustering value, the cracks are divided into Reticulate cracks and linear cracks. The invention reduces the computation complexity and the crack segment caused by details missing and false depth cracks, and can accurately restore to complete cracks.

【技术实现步骤摘要】
基于椭球立体表征的裂缝三维检测方法
本专利技术属于道路工程领域，涉及一种图像处理技术，具体提出了基于椭球立体表征的裂缝三维检测方法。
技术介绍
随着我国公路事业的迅猛发展，特别是高速公路里程的不断增长，公路在国民日常生活中发挥的作用越来越大，公路事业的发展极大地促进了我国经济的发展。由于公路及相关设施大都裸露于自然界，除受到交通荷载等各种力系的作用外，还直接或者间接地遭受光照、雨、雪、风、温度变化等各种自然因素的作用，由此会产生各种路面病害。路面裂缝类病害是路面各类病害中最常见、最早期发生并且最容易发生的病害之一，在恰当的时间对裂缝病害进行处理，不仅可以大大节约路面养护的成本，而且可以最大化地延长道路的使用周期。传统的裂缝检测主要依靠人工目视检测，这种检测方法不仅具有极高的危险性，而且检测结果有很强的主观性，同时成本高昂、耗时耗力。因此，研究高效且准确的路面自动化裂缝检测技术十分必要，对路面破损状况评估、道路养护有着很大的现实意义和重要价值。路面是一个三维实体，由于受到目前检测技术的限制，用于裂缝病害的检测技术大多是二维检测技术。二维检测技术依赖于路面图像的色彩或强度信息来识别路面裂缝病害，容易受到轮胎印迹、光照、油污、阴影、水渍等的影响。随着近年来公路事业的迅速发展，公路路面的评价体系也在不断完善，智能交通系统与路面管理系统对路面裂缝检测技术的要求也在不断提高，传统的二维检测技术由于存在的一些缺陷无法满足这些要求。三维检测系统中的成像系统可以利用高度信息对高度信息存在明显差别的路面裂缝进行全面的评价，克服了二维检测系统的缺陷，因而受到广泛的应用。目前现有的路面裂缝三维检测技术大都采用基于光学三维信息的获取方法来获取裂缝的深度信息。这种方法有一定的局限性，其局限性主要表现：在采集裂缝深度信息的过程中，存在相当一部分的虚假深度裂缝信息(所谓虚假深度裂缝是指实际路面裂缝中广泛存在的部分裂缝段被尘土等异物填充后的裂缝)，即采集到的路面裂缝深度信息中存在虚假深度现象。
技术实现思路
针对上述问题，本专利技术针对三维裂缝检测中虚假深度裂缝进行研究，根据三维路面裂缝的检测要求，针对虚假深度裂缝的数据特点，结合数字图像处理技术，提出了一种基于椭球立体表征的裂缝三维检测方法。本专利技术所采用的技术方案是，基于椭球立体表征的裂缝三维检测方法，包括：步骤1，读取路面三维图像数据矩阵；步骤2，对路面三维图像数据矩阵进行预处理，得到预处理后三维图像数据；步骤3，对预处理后的三维图像数据矩阵进行滤波处理，得到去噪后的三维图像数据矩阵；步骤4：提取去噪后的三维图像数据矩阵的每行数据，分别进行smooth平滑处理和曲线拟合处理，得到平滑曲线和拟合曲线；根据得到的平滑曲线和拟合曲线提取基准线，即为路面的横断面轮廓线，将每一行的数据处理后，最终提取出完整的路面裂缝；步骤5，将路面裂缝划分为多个裂缝区域，再将每个裂缝区域转化为附加深度的椭圆模型，每个附加深度的椭圆模型与其对应的裂缝区域具有等价归一化二阶中心距，求取每个附加深度的椭圆模型的表征参数，所述表征参数包括：椭圆中心点坐标、偏转角、长轴的长度、短轴的长度和深度；步骤6，根据表征参数对所有附加深度的椭圆模型进行聚类，得到聚类后的裂缝；所述根据表征参数对所有附加深度的椭圆模型进行聚类，包括：步骤61，对所有附加深度的椭圆模型进行编号，将所有附加深度的椭圆模型标记为未归类模型；步骤62，将未归类模型中长轴最长的附加深度的椭圆模型标记为当前起始聚类模型，当前起始聚类模型的编号为m，m＝1,2,...,M，M为所有附加深度的椭圆模型的数量；步骤621，将除过当前起始聚类模型后未归类模型中编号最小的附加深度的椭圆模型标记为待聚类模型，待聚类模型的编号为m′，m′≠m，m′＝1,2,...,M；步骤622，若当前起始聚类模型的表征参数与待聚类模型的表征参数满足聚类准则条件，则执行步骤623；否则，执行步骤624；所述聚类准则条件为：(条件1∪条件2)∩条件3∩条件4；其中，条件1为：(|Oa-Ob|≤30)∪(|Oa-Ob|≥150)；条件2为：条件3为：条件4为：XCa为中心点Ca的横坐标，YCa为中心点Ca的纵坐标，XCb为中心点Cb的横坐标，YCb为中心点Cb的纵坐标；其中，Ca为当前起始聚类模型的中心点、Oa为当前起始聚类模型的偏转角、La为当前起始聚类模型的长轴、Sa为当前起始聚类模型的短轴、Da为当前起始聚类模型的深度；Cb为待聚类模型的中心点、Ob为待聚类模型的偏转角、Lb为待聚类模型的长轴、Sb为待聚类模型的短轴、Db为待聚类模型的深度；步骤623，将待聚类模型标记为已归类模型，令待聚类模型的编号m′＝m，执行步骤624；步骤624，重复步骤621至623，直至未归类模型中的所有模型都被作为待聚类模型；步骤63，重复步骤62，直至所有未归类模型中的所有模型都被标记为已归类模型；步骤7，计算聚类后的裂缝的裂缝特征值，所述裂缝特征值包括裂缝长宽比R、裂缝块度H，根据裂缝特征值将裂缝分为网状裂缝和线性裂缝。进一步地，步骤5中求取每个附加深度的椭圆模型的表征参数，包括：步骤51，将路面裂缝划分成多个连通域，每个连通域即为一个裂缝区域；步骤52，任选一裂缝区域作为当前裂缝区域，计算当前裂缝区域对应的附加深度的椭圆模型的表征参数，包括：通过式(1)计算椭圆中心点坐标(μx,μy)：其中，xi、yi分别为当前裂缝区域中第i个像素点的坐标值，N为当前裂缝区域中像素点的总数；通过式(2)计算偏转角度O：其中，偏转角度O是指椭圆长轴与水平方向的夹角，Δ为当前裂缝区域中数值的二阶中心距，通过式(3)得到长轴的长度L：通过式(4)得到短轴的长度S：将当前裂缝区域中所有裂缝像素点处的裂缝深度值按照[depmin,depmin+0.01),[depmin+0.01,depmin+0.02),…,[depmax-0.01,depmax进行区间划分，找到每一裂缝深度值区间内所对应的裂缝像素点个数最多的区间d，将区间d的中值作为当前裂缝区域对应的附加深度的椭圆模型的深度D；其中，depmax为所有裂缝像素点处的裂缝深度值中最大的裂缝深度值，depmin为所有裂缝像素点处的裂缝深度值中最小的裂缝深度值；步骤53，重复步骤52，直至得到所有裂缝区域所对应的附加深度的椭圆模型的表征参数。进一步地，所述步骤7中根据裂缝特征值将裂缝分为网状裂缝和线性裂缝，包括：步骤71，若裂缝块度H满足：则该裂缝为网状裂缝；否则，执行步骤72；其中，hx1,hx2,…,hxm为水平方向上的裂缝块度，m为水平方向裂缝块度的数量，hy1,hy2,…,hyn为垂直方向上的裂缝块度，n为垂直方向裂缝块度的数量，I是裂缝区域的行数，J是裂缝区域的列数，h0是块度的阈值，h1是块度数量的阈值；步骤72，若裂缝长宽比R满足：R＜R0∩|r|＜r0，则该裂缝为网状裂缝；否则该裂缝为线性裂缝；其中，R0为长宽比的阈值，r0为相关系数的阈值，T为裂缝像素点的数量，xt为第t个裂缝像素点的行号，yt为第t个裂缝像素点的列号，为裂缝像素点行号的平均值，为裂缝像素点行号与列号乘积的平均值。本专利技术的有益效果是:(1)本专利技术在分析采集到的路面三维数据的结构特征后，提出本文档来自技高网...

【技术保护点】
基于椭球立体表征的裂缝三维检测方法，包括：步骤1，读取路面三维图像数据矩阵；步骤2，对路面三维图像数据矩阵进行预处理，得到预处理后三维图像数据；步骤3，对预处理后的三维图像数据矩阵进行滤波处理，得到去噪后的三维图像数据矩阵；步骤4：提取去噪后的三维图像数据矩阵的每行数据，分别进行smooth平滑处理和曲线拟合处理，得到平滑曲线和拟合曲线；根据得到的平滑曲线和拟合曲线提取基准线，即为路面的横断面轮廓线，将每一行的数据处理后，最终提取出完整的路面裂缝；其特征在于，还包括以下步骤：步骤5，将路面裂缝划分为多个裂缝区域，再将每个裂缝区域转化为附加深度的椭圆模型，每个附加深度的椭圆模型与其对应的裂缝区域具有等价归一化二阶中心距，求取每个附加深度的椭圆模型的表征参数，所述表征参数包括：椭圆中心点坐标、偏转角、长轴的长度、短轴的长度和深度；步骤6，根据表征参数对所有附加深度的椭圆模型进行聚类，得到聚类后的裂缝；所述根据表征参数对所有附加深度的椭圆模型进行聚类，包括：步骤61，对所有附加深度的椭圆模型进行编号，将所有附加深度的椭圆模型标记为未归类模型；步骤62，将未归类模型中长轴最长的附加深度的椭圆模型标记为当前起始聚类模型，当前起始聚类模型的编号为m，m＝1,2,...,M，M为所有附加深度的椭圆模型的数量；步骤621，将除过当前起始聚类模型后未归类模型中编号最小的附加深度的椭圆模型标记为待聚类模型，待聚类模型的编号为m′，m′≠m，m′＝1,2,...,M；步骤622，若当前起始聚类模型的表征参数与待聚类模型的表征参数满足聚类准则条件，则执行步骤623；否则，执行步骤624；所述聚类准则条件为：(条件1∪条件2)∩条件3∩条件4；其中，条件1为：(|Oa‑Ob|≤30)∪(|Oa‑Ob|≥150)；条件2为：...

【技术特征摘要】
1.基于椭球立体表征的裂缝三维检测方法，包括：步骤1，读取路面三维图像数据矩阵；步骤2，对路面三维图像数据矩阵进行预处理，得到预处理后三维图像数据；步骤3，对预处理后的三维图像数据矩阵进行滤波处理，得到去噪后的三维图像数据矩阵；步骤4：提取去噪后的三维图像数据矩阵的每行数据，分别进行smooth平滑处理和曲线拟合处理，得到平滑曲线和拟合曲线；根据得到的平滑曲线和拟合曲线提取基准线，即为路面的横断面轮廓线，将每一行的数据处理后，最终提取出完整的路面裂缝；其特征在于，还包括以下步骤：步骤5，将路面裂缝划分为多个裂缝区域，再将每个裂缝区域转化为附加深度的椭圆模型，每个附加深度的椭圆模型与其对应的裂缝区域具有等价归一化二阶中心距，求取每个附加深度的椭圆模型的表征参数，所述表征参数包括：椭圆中心点坐标、偏转角、长轴的长度、短轴的长度和深度；步骤6，根据表征参数对所有附加深度的椭圆模型进行聚类，得到聚类后的裂缝；所述根据表征参数对所有附加深度的椭圆模型进行聚类，包括：步骤61，对所有附加深度的椭圆模型进行编号，将所有附加深度的椭圆模型标记为未归类模型；步骤62，将未归类模型中长轴最长的附加深度的椭圆模型标记为当前起始聚类模型，当前起始聚类模型的编号为m，m＝1,2,...,M，M为所有附加深度的椭圆模型的数量；步骤621，将除过当前起始聚类模型后未归类模型中编号最小的附加深度的椭圆模型标记为待聚类模型，待聚类模型的编号为m′，m′≠m，m′＝1,2,...,M；步骤622，若当前起始聚类模型的表征参数与待聚类模型的表征参数满足聚类准则条件，则执行步骤623；否则，执行步骤624；所述聚类准则条件为：(条件1∪条件2)∩条件3∩条件4；其中，条件1为：(|Oa-Ob|≤30)∪(|Oa-Ob|≥150)；条件2为：条件3为：条件4为：XCa为中心点Ca的横坐标，YCa为中心点Ca的纵坐标，XCb为中心点Cb的横坐标，YCb为中心点Cb的纵坐标；其中，Ca为当前起始聚类模型的中心点、Oa为当前起始聚类模型的偏转角、La为当前起始聚类模型的长轴、Sa为当前起始聚类模型的短轴、Da为当前起始聚类模型的深度；Cb为待聚类模型的中心点、Ob为待聚类模型的偏转角、Lb为待聚类模型的长轴、Sb为待聚类模型的短轴、Db为待聚类模型的深度；步骤623，将待聚类模型标记为已归类模型，令待聚类模型的编号m′＝m，执行步骤624；步骤624，重复步骤621至623，直至未归类模型中的所有模型都被作为待聚类模型；步骤63，重复步骤62，直至所有未归类模型中的所有模型都被标记为已归类模型；步骤7，计算聚类后的裂缝的裂缝特征值，所述裂缝特征值包括裂缝长宽比R、裂缝块度H，根据裂缝特征值将裂缝分为网状裂缝和线性裂缝。2.如权利要求1所述的基于椭球立体表征的裂缝三维检测方法，其特征在于，步骤5中求取每个附加深度的椭圆模型的表征参数，包括：步骤51，将路面裂缝划分成多个连通域，每个连通域即为一个裂缝区域；步骤52，任选一裂缝区域作为当前裂缝区域，计算当前裂缝区域对应的附加深度...

【专利技术属性】
技术研发人员：沙爱民，李伟，孙朝云，邵娜娜，郝雪丽，张欣，马志丹，
申请(专利权)人：长安大学，
类型：发明
国别省市：陕西,61