CN103061237A - 一种水泥混凝土路面错台量的检测算法 - Google Patents

Abstract

本发明公开了一种水泥混凝土路面错台量的检测算法，具体包括如下步骤：输入图像三维数据矩阵，并对数据进行滤波处理；依次取滤波处理之后的三维图像数据矩阵O″的每一行，逐行求每一行所对应的错台量，即找出该行的所有拐点后求该行的错台量：在拐点矩阵G中，取横坐标相差为M的拐点两两作差；对所有行处理完毕，即得高度差矩阵，该高度差矩阵的最大值即为该行相应的错台量；逐行求得错台量h₁,h₂,…,h_m，对此m个错台量求算术平均，即得该图像采集区域路面错台量。本发明计算简单、运行时间短，无需人工参与。采用面测量，只需输入采集到水泥混凝土路面的三维图像数据，即可完成对路面错台量的检测，因此，该检测算法效率高、检测精确。

Description

一种水泥混凝土路面错台量的检测算法

技术领域

本发明属于道路工程领域，特别是指一种水泥混凝土路面错台量的检测算法。

背景技术

水泥混凝土路面错台量是指水泥混凝土路面中横向接缝或裂缝处相邻两块水泥板的高程差。错台产生的原因主要有以下几种：路基基层碾压不密实，强度不足，致使基层在行车荷载作用下发生塑性累计位移；局部地基不均匀沉降；相邻板间的传荷能力下降；水浸入基层，行车荷载使路面板产生泵吸现象，动水将面板与基层间的碎屑抛向后方，把后方的板抬起等等。错台是影响水泥混凝土路面平整度、寿命以及加铺后反射裂缝的重要因素，也是水泥混凝土路面技术状况评定、养护、加铺设计时必须考虑的一个重要指标。因此，如何对错台量进行实时、准确和高效的检测，成为公路管理养护部门十分关注的问题。

目前，常用对错台量的检测方法有人工法和自动检测法。人工法利用直尺、游标卡尺或水准仪等实现，速度慢，交通干扰大，精度低；自动检测法利用错台的存在会显著影响国际平整度指标（IRI）的特性，有激光断面仪、超声波断面仪等，断面仪价格昂贵，并且不能对错台量准确测量。综上，现有检测技术存在误差大、效率低、造价昂贵等问题，研究一种效率高、检测精确的水泥混凝土路面错台量的检测算法是十分有必要的。

发明内容

有鉴于此，本发明的主要目的在于，提供一种水泥混凝土路面错台量的检测算法，从而实现对水泥混凝土路面错台量的实时、高效和精确的检测。

为了达到上述目的，本发明采用如下的技术方案：

一种水泥混凝土路面错台量的检测算法，具体包括如下步骤：

步骤1：输入图像三维数据矩阵，并对数据进行滤波处理；

步骤2：依次取滤波处理之后的三维图像数据矩阵O″的每一行，逐行求每一行所对应的错台量，包括以下步骤：

（1）找出该行的所有拐点：从每一行的第二个数据点开始，到倒数第二个数据点为止，依次判断，如果该数据点减去前一个数据点的差值与后一个数据点减去该数据点的差值异号，则说明该数据点是拐点，保存拐点的高度值，对所有行处理完毕，得到拐点矩阵G；

（2）求该行的错台量：在拐点矩阵G中，取横坐标相差为M（M取5~10）的拐点两两作差；对所有行处理完毕，即得高度差矩阵，该高度差矩阵的最大值即为该行相应的错台量；

步骤3：经过步骤2逐行求得错台量h₁,h₂,…,h_m，对此m个错台量求算术平均，即得该图像采集区域路面错台量；

步骤4：结束。

进一步的，所述步骤1的输入图像三维数据矩阵，并对数据进行滤波处理具体包括如下步骤：

将图像采集设备采集到的三维图像数据矩阵输入计算机，计算机读取三维图像数据矩阵

$O_{m\×n}=\left[\begin{array}{ccccccccccc}{z}_{11}& z_{12}& z_{13}& .& .& .& z_{1j}& .& .& .& z_{1n}\\ z_{21}& z_{22}& z_{23}& .& .& .& z_{2j}& .& .& .& z_{2n}\\ .& .& .& & & & .& & & & .\\ .& .& .& & & & .& & & & .\\ .& .& .& & & & .& & & & .\\ z_{i1}& z_{i2}& z_{i3}& .& .& .& z_{\mathrm{ij}}& .& .& .& z_{\mathrm{in}}\\ .& .& .& & & & .& & & & .\\ .& .& .& & & & .& & & & .\\ .& .& .& & & & .& & & & .\\ z_{m1}& z_{m2}& z_{m3}& .& .& .& z_{\mathrm{mj}}& .& .& .& z_{\mathrm{mn}}\end{array}\right],(i=\mathrm{1,2}\·\·\·m,j=\mathrm{1,2}\·\·\·n)$

z_ij表示行号为i，列号为j所对应的图像高度数据。

对数据的滤波处理包括粗滤波和细滤波两个部分：

（1）对数据的粗滤波部分包括逐行消除图像的噪声点和逐列消除图像的噪声点两个过程，步骤如下：

1）逐行取三维图像数据矩阵O的数据，如第i行R_i＝(z_i1,z_i2,…z_in),i＝(1,2,…m),求出第i行数据R_i的算术平均值

和标准差S_i；

2）依次取该行中的每一个数据z_ij，如果满足

则用算术平均值

代替该点数据值z_ij，k是行滤波系数，k取3~8；

3）逐行消除图像噪声点后，用同样的方法逐列处理，再次消除噪声点，对所有行处理完毕，得到粗滤波之后的三维图像数据矩阵O′。

（2）对数据的细滤波部分是针对错台三维图像数据中特有的暗点所进行的处理，具体步骤如下：

1）依次取粗滤波之后的三维图像数据矩阵O′的每一行，并定义一个可变长度N，宽度为1的滑动窗口H_N，H_N＝[1,1,…,1]_N，其中，N取11~21之间的奇数；

2）让滑动窗口H_N的中心点从每行的第一个数据点滑动到最后一个数据点，如果在滑动窗口范围内该数据点的前

个数据点与该数据点的差值的绝对值或后

个数据点与该数据点的差值的绝对值有一端均不超过设定阈值T（T取3~5），其中，对于每一行的第一个数据点只需判断其后面数据点，对于最后一个数据点只需判断其前面数据点，则该数据点具有连通性，否则该数据点即为错台三维图像数据中的暗点，暗点需要被删除，并将该暗点后面的数据依次向前移动，对所有行处理完毕，得到细滤波之后的三维图像数据矩阵O″。

本发明提出的方法具有以下优点：

1、计算简单、运行时间短，适合在实时的系统中采用。

2、无需人工参与，克服了人工检测方法具有的劳动强度大、安全性低、行车受干扰、工作效率低和检测精确度较低的缺点。

3、采用面测量，只需输入采集到水泥混凝土路面的三维图像数据，即可完成对路面错台量的检测，因此，该检测算法效率高、检测精确。

4、消除了暗点的影响，适用范围广。

5、可为水泥混凝土路面的养护管理提供有力的信息支持，提高了公路养护和管理水平，同时，为进一步开发公路检测设备，改变我国公路工程检测设备绝大部分依赖进口的现状，节省资源，培养自己的研究和开发技术力量打下了人力与技术基础。

附图说明

图1是本发明的算法的总流程图。

图2是粗滤波算法流程图。

图3是细滤波算法流程图。

图4计算一行数据错台量算法流程图。

以下结合附图和具体实施方式对本发明作进一步解释说明。

具体实施方式

参见图1-图4，本发明的一种水泥混凝土路面错台量的检测算法，具体包括如下步骤：

步骤1：输入图像三维数据矩阵，并对数据进行滤波处理；

将图像采集设备采集到的三维图像数据矩阵输入计算机，计算机读取三维图像数据矩阵O_m×n，本实施例中取m＝1000，n＝1536。

（1）对数据进行粗滤波处理，包括逐行消除图像的噪声点和逐列消除图像的噪声点两个过程，步骤如下：

1）逐行取三维图像数据矩阵O的每一行数据，如第i行R_i＝(z_i1,z_i2,…z_i，1536),i＝(1,2,…1000),求出R_i的算术平均值

和标准差S_i；

2）依次取该行中的每一个数据z_ij，如果满足则用算术平均值代替该点数据值z_ij，这里k＝3；

3）逐行消除图像噪声点后，用同样的方法逐列处理，再次消除噪声点，对所有行处理完毕，得到滤波之后的三维图像数据矩阵O′；

（2）对数据的细滤波部分主要是针对错台三维图像数据中特有的暗点所进行的处理，步骤如下：

1）依次取粗滤波之后的三维图像数据矩阵O′的每一行，并定义一个滑动窗口H_N，H_N＝[1,1,…,1]_N，这里N取11；

2）让滑动窗口H_N的中心点从每行的第一个数据点滑动到最后一个数据点，如果在滑动窗口范围内该数据点的前

个数据点与该数据点的差值的绝对值或后

个数据点与该数据点的差值的绝对值有一端均不超过设定阈值T（T＝5），其中，对于每一行的第一个数据点只需判断其后面数据点，对于最后一个数据点只需判断其前面数据点，则该数据点具有连通性，否则该数据点即为错台三维图像数据中的暗点，若发现暗点，则删除暗点，并将该暗点后面的数据依次向前移动，对所有行处理完毕，得到细滤波之后的三维图像数据矩阵O″。

步骤2：依次取滤波处理之后的三维图像数据矩阵O″的每一行，逐行求每一行所对应的错台量，包括以下步骤：

（1）找出该行的所有拐点：从每一行的第二个数据点开始，到倒数第二个数据点为止，依次判断，如果该数据点减去前一个数据点的差值和后一个数据点减去该数据点的差值异号，则说明该数据点是拐点，保存拐点的高度值，对所有行处理完毕，得到拐点矩阵G；

（2）求该行的错台量：在拐点矩阵G中，取横坐标相差为M（M＝5）的拐点两两作差，对所有行处理完毕，即得高度差矩阵，该高度差矩阵的最大值，即为该行相应的错台量；

步骤3：经过步骤2逐行求得错台量h₁,h₂,…,h₁₀₀₀，对此1000个错台量求算术平均，即得该图像采集区域路面错台量。

对需要进行计算错台量的水泥混泥土路面，抽样采集该路面的三维图像数据，对于待评估路面上的每幅图像数据，根据上述算法求出采集区域水泥混凝土路面错台的错台量。

Claims (2)

1.一种水泥混凝土路面错台量的检测算法，其特征在于，具体包括如下步骤：

步骤1：输入图像三维数据矩阵，并对数据进行滤波处理；

步骤2：依次取滤波处理之后的三维图像数据矩阵O″的每一行，逐行求每一行所对应的错台量，包括以下步骤：

（1）找出该行的所有拐点：从每一行的第二个数据点开始，到倒数第二个数据点为止，依次判断，如果该数据点减去前一个数据点的差值与后一个数据点减去该数据点的差值异号，则说明该数据点是拐点，保存拐点的高度值；对所有行处理完毕，得到拐点矩阵G；

（2）求该行的错台量：在拐点矩阵G中，取横坐标相差为M（M取5~10）的拐点两两作差，对所有行处理完毕，即得高度差矩阵，该高度差矩阵的最大值即为该行相应的错台量；

步骤3：经过步骤2逐行求得错台量h₁,h₂,…,h_m，对此m个错台量求算术平均，即得该图像采集区域路面错台量；

步骤4：结束。

2.如权利要求1所述的方法，其特征在于，所述步骤1的输入图像三维数据矩阵，并对数据进行滤波处理具体包括如下步骤：

将图像采集设备采集到的三维图像数据矩阵输入计算机，计算机读取三维图像数据矩阵

$O_{m\×n}=\left[\begin{array}{ccccccccccc}{z}_{11}& z_{12}& z_{13}& .& .& .& z_{1j}& .& .& .& z_{1n}\\ z_{21}& z_{22}& z_{23}& .& .& .& z_{2j}& .& .& .& z_{2n}\\ .& .& .& & & & .& & & & .\\ .& .& .& & & & .& & & & .\\ .& .& .& & & & .& & & & .\\ z_{i1}& z_{i2}& z_{i3}& .& .& .& z_{\mathrm{ij}}& .& .& .& z_{\mathrm{in}}\\ .& .& .& & & & .& & & & .\\ .& .& .& & & & .& & & & .\\ .& .& .& & & & .& & & & .\\ z_{m1}& z_{m2}& z_{m3}& .& .& .& z_{\mathrm{mj}}& .& .& .& z_{\mathrm{mn}}\end{array}\right],(i=\mathrm{1,2}\·\·\·m,j=\mathrm{1,2}\·\·\·n)$

z_ij表示行号为i，列号为j所对应的图像高度数据。

对数据的滤波处理包括粗滤波和细滤波两个部分：

（1）对数据的粗滤波部分包括逐行消除图像的噪声点和逐列消除图像的噪声点两个过程，步骤如下：

1）逐行取三维图像数据矩阵O的数据，如第i行R_i＝(z_i1,z_i2,…z_in),i＝(1,2,…m),求出第i行数据R_i的算术平均值

和标准差S_i；

2）依次取该行中的每一个数据z_ij，如果满足

则用算术平均值代替该点数据值z_ij，k是行滤波系数，k取3~8；

3）逐行消除图像噪声点后，用同样的方法逐列处理，再次消除噪声点，对所有行处理完毕，得到粗滤波之后的三维图像数据矩阵O′。

（2）对数据的细滤波部分是针对错台三维图像数据中特有的暗点所进行的处理，具体步骤如下：

1）依次取粗滤波之后的三维图像数据矩阵O′的每一行，并定义一个可变长度N，宽度为1的滑动窗口H_N，H_N＝[1,1,…,1]_N，其中，N取11~21之间的奇数；

2）让滑动窗口H_N的中心点从每行的第一个数据点滑动到最后一个数据点，如果在滑动窗口范围内该数据点的前

个数据点与该数据点的差值的绝对值或后

个数据点与该数据点的差值的绝对值有一端均不超过设定阈值T（T取3~5），其中，对于每一行的第一个数据点只需判断其后面数据点，对于最后一个数据点只需判断其前面数据点，则该数据点具有连通性，否则该数据点即为错台三维图像数据中的暗点，暗点需要被删除，并将该暗点后面的数据依次向前移动，对所有行处理完毕，得到细滤波之后的三维图像数据矩阵O″。

Images