CN104537651B - 一种路面图像中裂缝比例检测方法及系统 - Google Patents

Abstract

本发明公开了一种路面图像中裂缝比例检测算法及系统，读取路面三维数据矩阵；对路面三维数据矩阵进行去噪处理；由去噪后的路面三维数据矩阵识别路面裂缝，得到最终裂缝图像；从最终裂缝图像中提取裂缝种子点并进行区域生长，得到裂缝二值化图像；将裂缝二值化图像等分为多个正方形；计算路面图像中裂缝所占比例。本发明只需输入采集到的三维数据矩阵，即可快速而准请地计算出裂缝在整幅图像中的比例。

Description

一种路面图像中裂缝比例检测方法及系统

技术领域

本发明属于道路工程领域，具体涉及一种路面图像中裂缝比例检测算法及系统。

背景技术

根据我国的《公路养护技术规范》，我国高速公路路面评价包括四部分内容，即道路路面的平整度(行驶舒适度)、路面破损(路面破损状况指数)、弯沉(路面的结构强度)和抗滑性能(安全性)。其中，路面破损状况指数(PCI)是决策养护方案的最重要的数据，它不仅反应了路面结构的完好程度，又直接影响道路的服务寿命，为了了解和掌握路面使用性能的衰变情况，以便及时采用相应的养护和改善措施，延缓其衰变或恢复其使用性能，就必须对路面破损状况进行正确的评价，这是科学预测路面使用性能、合理制定养护维修计划、进行投资决策的重要依据之一，是路面养护中最重要的环节。而路面裂缝类病害作为路面破损检测的重要内容，其自动化检测一直是公路路面破损检测的热点和难点。

目前，国内外的裂缝检测技术大都是采用CCD摄像机获取路面图像，然后对采集到的二维图像进行后续的处理来识别裂缝，进而对裂缝的类别进行判断和提取裂缝的特征信息。但是，采集到的二维图像往往会受到路面光照、油污、建筑物以及树等的影子的影响，这样会使得基于二维图像的裂缝分析受到很大的干扰，大大影响裂缝占路面比例检测的准确度，因此研究一种效率高、检测精确高的路面图像中裂缝比例检测技术是十分有必要的。

发明内容

针对上述现有技术中存在的缺陷或不足，本发明的一个目的在于，提供一种路面图像中裂缝比例检测算法。

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

一种路面图像中裂缝所占比例检测算法，具体包括如下步骤：

步骤1：读取路面三维数据矩阵；

步骤2：对路面三维数据矩阵进行去噪处理，得到去噪后的路面三维数据矩阵；

步骤3：由去噪后的路面三维数据矩阵识别路面裂缝，得到最终裂缝图像；

步骤4：从最终裂缝图像中提取裂缝种子点并进行区域生长，得到裂缝二值化图像。

步骤5：将步骤4得到的裂缝二值化图像等分为多个正方形；

步骤6：计算路面图像中裂缝所占比例：统计包含有裂缝像素点的正方形个数和网格划分后的整幅图像中的正方形块个数，计算前者在后者中所占的比例即得到路面图像中裂缝所占比例。

进一步的，所述步骤1中所述的路面三维数据矩阵O_m×n如下：

(i＝1,2…m,j＝1,2…n)

z_ij表示行号为i，列号为j所对应的路面高度。

进一步的，所述步骤2对路面三维数据矩阵进行去噪处理的操作如下：

画出路面高度直方图，图中横坐标为路面三维数据矩阵中的高度数据，纵坐标为路面三维数据矩阵中的对应于每个高度数据段的元素个数；将路面高度直方图中两个高度数据段所对应的元素分别标记；其他高度数据段对应的元素标记为噪声点；对步骤22中标记的噪声点进行滤波处理。

进一步的，所述步骤3具体包括如下步骤：

步骤31：对步骤2得到的去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的水平单相扫描，得到裂缝图像I1；具体是：将矩阵中每一行的数据个数N等分得到N个数据段，N为能被每行元素的个数整除的数；对每个数据段中的元素进行曲线拟合，用拟合值减去对应元素的值得到差值Δ；然后取多个阈值将所有的Δ值分为多个段，将每一段中对应的元素的值用对应的阈值替换，得到裂缝图像I1；

步骤32：将步骤2得到的去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的垂直单相扫描，得到裂缝图像I2；具体是：将矩阵中每列的数据个数M等分得到M个数据段，M是能被每列元素的个数整除的数；对每个数据段中的元素进行曲线拟合，用拟合值减去其对应的元素的值Δ；然后将所有的Δ值分为多个段，将每一段中对应的元素的值用对应的阈值替换；

步骤33：将裂缝图像I1与裂缝图像I2取交集，得到交集裂缝图像I3；

步骤34：将裂缝图像I1与裂缝图像I2取并集，得到并集裂缝图像I4；

步骤35：将裂缝图像I1、裂缝图像I2、交集裂缝图像I3、并集裂缝图像I4分别划分为m*n个小块，m是能被每行元素的个数整除的数，n是能被每列元素的个数整除的数；

对于交集裂缝图像I3，从左到右从上而下扫描每个小块，对于每个小块，分别计算该小块在裂缝图像I1、裂缝图像I2以及交集裂缝图像I3中对应位置的小块中黑色像素点所占的比例R1、R2和R3；计算R3与R1的相对误差R31＝|R1-R3|/R1，若R1＝0，则R31＝0；计算R3与R2的相对误差R32＝|R2-R3|/R2，若R2＝0，则R32＝0；计算交集裂缝图像I3中该小块与裂缝图像I1、I2中对应小块的相似度R312＝0.5*R31+0.5*R32。

对于并集裂缝图像I4，从左到右从上而下扫描每一个小块，对于每个小块，分别计算其在裂缝图像I1、裂缝图像I2、并集裂缝图像I4中对应的小块中黑色像素点所占的比例R1、R2和R4，计算R4与R1的相对误差记为R41＝|R1-R4|/R1，若R1＝0，则R41＝0；计算R4与R2的相对误差R42＝|R2-R4|/R2，若R2＝0，则R42＝0；R412＝0.5*R41+0.5*R42。

步骤36：对于每个小块，比较其对应的R312与R412的大小，若R312大于等于R412，则取I4中的相应小块作为最终裂缝图像中的相应位置的小块。否则取I3中的相应小块作为最终裂缝图像中的相应位置的小块，得到最终裂缝图像I5。

本发明的另一个目的是，提供一种路面图像中裂缝所占比例检测系统，包括如下依次相连接的模块：

矩阵读取模块，是用于实现路面三维数据矩阵读取的模块；

去噪模块，是用于对路面三维数据矩阵进行去噪处理得到去噪后的路面三维数据矩阵的模块；

路面裂缝识别模块，是用于从去噪后的路面三维数据矩阵得到最终裂缝图像的模块。

裂缝二值化模块，是用于从最终裂缝图像中提取裂缝种子点并进行区域生长，得到裂缝二值化图像的模块；

图像网格划分模块，是用于将裂缝二值化图像等分为多个正方形的模块；

裂缝比例计算模块：是用于实现以下功能的模块：统计包含有裂缝像素点的正方形个数和网格划分后的整幅图像中的正方形块个数，计算前者在后者中所占的比例得到路面图像中裂缝所占比例。

进一步的，所述矩阵读取模块中所述的路面三维数据矩阵O_m×n如下：

(i＝1,2…m,j＝1,2…n)

z_ij表示行号为i，列号为j所对应的路面高度。

进一步的，所述路面裂缝识别模块是用于实现以下功能的模块：画出路面高度直方图，图中横坐标为路面三维数据矩阵中的高度数据，纵坐标为路面三维数据矩阵中的对应于每个高度数据段的元素个数；将路面高度直方图中两个高度数据段所对应的元素分别标记；其他高度数据段对应的元素标记为噪声点；对步骤22中标记的噪声点进行滤波处理，得到去噪后的路面三维数据矩阵。

进一步的，所述裂缝二值化模块是用于实现以下功能的模块：

对去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的水平单相扫描，得到裂缝图像I1；具体是：将矩阵中每一行的数据个数N等分得到N个数据段，N为能被每行元素的个数整除的数；对每个数据段中的元素进行曲线拟合，用拟合值减去对应元素的值得到差值Δ；然后取多个阈值将所有的Δ值分为多个段，将每一段中对应的元素的值用对应的阈值替换，得到裂缝图像I1；

将去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的垂直单相扫描，得到裂缝图像I2；具体是：将矩阵中每列的数据个数M等分得到M个数据段，M是能被每列元素的个数整除的数；对每个数据段中的元素进行曲线拟合，用拟合值减去其对应的元素的值Δ；然后将所有的Δ值分为多个段，将每一段中对应的元素的值用对应的阈值替换；

将裂缝图像I1与裂缝图像I2取交集，得到交集裂缝图像I3；

将裂缝图像I1与裂缝图像I2取并集，得到并集裂缝图像I4；

将裂缝图像I1、裂缝图像I2、交集裂缝图像I3、并集裂缝图像I4分别划分为m*n个小块，m是能被每行元素的个数整除的数，n是能被每列元素的个数整除的数；

对于交集裂缝图像I3，从左到右从上而下扫描每个小块，对于每个小块，分别计算该小块在裂缝图像I1、裂缝图像I2以及交集裂缝图像I3中对应位置的小块中黑色像素点所占的比例R1、R2和R3；计算R3与R1的相对误差R31＝|R1-R3|/R1，若R1＝0，则R31＝0；计算R3与R2的相对误差R32＝|R2-R3|/R2，若R2＝0，则R32＝0；计算交集裂缝图像I3中该小块与裂缝图像I1、I2中对应小块的相似度R312＝0.5*R31+0.5*R32。

对于并集裂缝图像I4，从左到右从上而下扫描每一个小块，对于每个小块，分别计算其在裂缝图像I1、裂缝图像I2、并集裂缝图像I4中对应的小块中黑色像素点所占的比例R1、R2和R4，计算R4与R1的相对误差记为R41＝|R1-R4|/R1，若R1＝0，则R41＝0；计算R4与R2的相对误差R42＝|R2-R4|/R2，若R2＝0，则R42＝0；R412＝0.5*R41+0.5*R42。

对于每个小块，比较其对应的R312与R412的大小，若R312大于等于R412，则取I4中的相应小块作为最终裂缝图像中的相应位置的小块。否则取I3中的相应小块作为最终裂缝图像中的相应位置的小块，得到最终裂缝图像I5。

与现有技术相比，本发明提出的方法具有以下优点：

1、采用面测量，只需输入采集到的路面三维数据矩阵，即可完成裂缝的种类判别以及特征提取，其效率高、检测精确，适合在实时系统中采用。

2、对于所需要的裂缝数据以及正常的路面数据保持原值不予进行滤波处理，而只对毛刺点数据进行滤波去噪，从而更好的保证了去噪处理后的三维信息数据的准确性。

3、由去噪处理后的三维数据矩阵直接识别裂缝，避免了由二维图像识别裂缝时由于路面光照、油污、建筑物以及树等的影子的影响而造成的裂缝识别率低的问题。本发明采用基于双相扫描的交并拼接方法进行裂缝识别，兼顾了单独采用交集裂缝图像以及单独采用并集裂缝图像的优缺点，在交集裂缝图像与并集裂缝图像中择优取块而进行拼接，使得最终得到的裂缝图像在尽可能的保留裂缝的全部信息的同时，也使得噪声点减到最小，这样得到的裂缝图像对于后续的特征提取更有实用意义。

4、采用合适数量的正方形将裂缝二值图像进行网格化，并将裂缝区域标出，使得计算操作简单、直观、高效，便于在实时系统中采用。

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

附图说明

图1是本发明的路面图像中裂缝比例检测算法的流程图。

图2是原始路面三维数据的三维显示图。

图3是去噪后的路面三维数据的三维显示图。

图4是由三维数据直接识别出来的最终裂缝图像。

图5是提取种子点算法流程图。

图6是提取出来的种子点图像。

图7是区域生长后得到的裂缝二值化图像。

图8是4*4标定后的裂缝图像。

图9是5*5标定后的裂缝图像。

图10是8*8标定后的裂缝图像。

图11是10*10标定后的裂缝图像。

图12是16*16标定后的裂缝图像。

图13是20*20标定后的裂缝图像。

图14是25*25标定后的裂缝图像。

图15是原始路面三维数据高度直方图。

图16是本发明的路面图像中裂缝比例检测系统的功能结构框图。

具体实施方式

以下是发明人提供的具体实施例，需要说明的是，所提供的实施例是对本发明的进一步解释说明，本发明的保护范围并不限于所给实施例。

参见图1，遵循本发明的技术方案，本实施例的路面裂缝在路面图像中所占比例的检测算法包括如下步骤：

步骤2：对路面三维数据矩阵进行去噪处理，得到去除噪声后的路面三维数据矩阵；包括如下步骤：

步骤21：画出路面高度直方图：即对路面三维数据矩阵的处于各个高度数据段的元素个数进行统计，在该路面高度直方图中，横坐标为路面三维数据矩阵中的高度数据，纵坐标为路面三维数据矩阵中的对应于每个高度数据段的元素个数；

步骤22：进行元素标记：路面高度直方图存在两个波峰，一个波峰所在的高度数据段代表正常路面的元素，另一个波峰所在的高度数据段代表裂缝的元素，将这两个高度数据段所对应的元素分别在路面三维数据矩阵中进行标记；其他高度数据段对应的元素在路面三维数据矩阵中标记为噪声点。

步骤23：对步骤22中标记的噪声点进行滤波处理。1)逐行处理：对噪声点元素所在行的所有数据求算术平均值，然后将该噪声点元素用该算术平均值代替。2)逐列 处理：在逐行处理后的路面三维数据矩阵中，对噪声点元素所在列的所有数据求算术平均值，然后将该噪声点元素用该列的算术平均值代替，得到去噪后的路面三维数据矩阵。

步骤3：由去噪后的三维数据矩阵直接识别出路面裂缝；

对步骤2得到的去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的水平单相扫描，得到裂缝图像I1。将每一行的数据8等分得到8个数据段。对每个数据段中的元素进行曲线拟合，用元素的拟合值减去其对应的值得到差值Δ(亦即该元素在去噪后的路面三维数据矩阵中的值)。然后进行阈值分割，即选取不同的阈值对Δ的值进行判断，根据Δ所处的阈值范围赋予矩阵中的元素新的数值。借鉴裂缝的破损等级划分为轻、中、重的思想，本实施例的阈值设置为3个，采用基于行标准值与行平均值组合的自适应阈值，即对某一行的数据进行处理时，该行数据的阈值是通过自动计算该行数据的平均值与标准差，然后由该平均值与该标准差组合成三个不同的阈值。这三个阈值将数据分为四个深度范围，根据Δ所处的深度范围可将矩阵中的元素分别赋予四个不同的数值(即四值化处理)：0,64,128,255(0-255之间平均分为四个值)。然后以这些新值作为图像的灰度值，显示图像，便可得到水平单相扫描检测后的四值化裂缝图像I1。由于去噪后的三维高度矩阵中的高度数据有大小之分，亦即采集到的裂缝有深浅之分，因此I1中初步识别出的裂缝颜色也有深浅之分。借鉴裂缝破损程度的轻、中、重等级划分思想，将I1中裂缝点灰度值处于[0,64)认为是“重”裂缝点，处于[64,128)认为是“中”裂缝点，处于[128,255]认为是“轻”裂缝点。对步骤2得到的去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的垂直单相扫描，得到裂缝图像I2，与水平单相扫描检测方法相似，不再赘述。

分别将裂缝图像I1与裂缝图像I2取交集和并集，得到交集裂缝图像I3和并集裂缝图像I4，然后将I1、I2、I3、I4分别划分为20*20个小块。然后对于每个小块，比较其对应的R312与R412的大小，若R312大于等于R412，则说明该小块在I4中的相应小块比在I3中的相应小块更接近I1与I2中的相应小块，因此，对于该小块，取I4中的相应小块作为最终裂缝图像中的相应位置的小块。否则，对于该小块，取I3中的相应小块作为最终裂缝图像中的相应位置的小块；得到最终裂缝图像I5。

步骤3中，交集裂缝图像I3中噪声点减弱，确保噪声点之外的肯定是裂缝信息；并集裂缝图像I4中噪声点加重，但确保含有所有的裂缝信息。因此，单独采用交集裂缝图像或者单独采用并集裂缝图像都不能很好地识别裂缝。本发明采用基于双相扫描的 交并拼接方法进行裂缝识别，兼顾了单独采用交集裂缝图像以及单独采用并集裂缝图像的优缺点，在交集裂缝图像I3与并集裂缝图像I4中择优取块而进行拼接，使得最终裂缝图像在尽可能的保留裂缝的全部信息的同时，也使得噪声点减到最小，这样得到的裂缝图像对于后续的特征提取更有实用意义。

步骤4：提取裂缝种子点并进行区域生长以得到裂缝二值化图像；

为了将步骤3得到最终裂缝图像中存在的一些非裂缝点去除，需要进一步分割图像。在分割图像时，存在不可避免的过度分割，即在去除非裂缝点的同时也会去掉一些浅裂缝点。由于过度分割后保留下来的深裂缝点与被分割掉的浅裂缝点存在一定的灰度相似性以及空间连续性，而与分割掉的非裂缝点虽有一定的灰度相似性但是空间连续性极差，因此本发明的步骤4中，将过度分割后的深裂缝点当作种子点来进行区域生长，以恢复分割掉的裂缝点同时避免非裂缝点的再次出现。在步骤3得到的裂缝图像中，颜色越深(即灰度值越小)说明所在部位裂缝越深，因此可以确定该裂缝图像中一些颜色较深的像素点一定是位于裂缝上的像素点。这样只需对步骤3得到的最终裂缝图像选取合适的阈值进行二值化处理便可以得到种子点。经试验，本实施例中提取种子点时的阈值取为64。种子点提取流程图如图5所示，提取得到的种子点如图6所示，将提取得到的种子点放到步骤3提取出来的裂缝图像中进行区域生长，生长准则是将种子点与其8领域的像素点逐个比较，判断待生长像素与种子点像素的灰度值的差的绝对值是否在设定的阈值范围内，若在，则认为该待生长像素包含在种子区域内，否则认为该待生长像素是背景点。对新的种子点继续相同的区域生长，直到没有新的像素包括到种子里，得到如图7所示的裂缝二值化图像。

步骤5：将步骤4得到的裂缝二值化图像等分为多个正方形；

为了选取合适的正方形进行了试验，将步骤4得到的裂缝二值化图像分别划分为4*4、5*5、8*8、10*10、16*16、20*20、25*25的正方形块。4*4中裂缝如图8所示，裂缝所占比例为43.75％；5*5中裂缝如图9所示，裂缝所占比例为28％；8*8中裂缝如图10所示，裂缝所占比例为20.31％；10*10中裂缝如图11所示，裂缝所占比例为15％；16*16中裂缝标定结果图如图12所示，裂缝所占比例为9.77％；20*20中裂缝如图13所示，裂缝所占比例为7.50％；25*25中裂缝如图14所示，裂缝所占比例为7.20％。分析上述不同的划分方法的结果，考虑到结果的准确性以及程序的运行速度等因素，本系统中采用将裂缝图像划分为20*20，亦即将裂缝图像划分为400个小正方形块，每个小正方形块的大小 为50*50(以像素为单位)，然后通过对划分后的裂缝图像中的小正方形块从上到下从左到右的扫描，将有裂缝点的小正方形块标定为红色，这样便可标定出裂缝所在的区域，裂缝如图14所示。

步骤6：计算路面图像中裂缝所占比例。

由图13可知，含有裂缝的正方形块的个数为num＝30(个)，而整幅裂缝图像中小正方形块的个数为sum＝400(个)，得到裂缝占图像比例r＝num/sum＝30/400＝0.075＝7.50％。

遵循本发明的技术方案，进行了本发明的路面图像中裂缝所占比例检测系统的开发，以下是运行该系统的结果与传统的水准仪法的结果进行比较。如图16所示，为本发明的系统的功能结构框图。

Claims (6)

1.一种路面图像中裂缝比例检测方法，其特征在于，具体包括如下步骤：

步骤1：读取路面三维数据矩阵；

步骤2：对路面三维数据矩阵进行去噪处理，得到去噪后的路面三维数据矩阵；

步骤3：由去噪后的路面三维数据矩阵识别路面裂缝，得到最终裂缝图像；所述步骤3具体包括如下步骤：

步骤31：对步骤2得到的去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的水平单相扫描，得到裂缝图像I1；具体是：将矩阵中每一行的数据个数N等分得到N个数据段，N为能被每行元素的个数整除的数；对每个数据段中的元素进行曲线拟合，用拟合值减去对应元素的值得到水平单相扫描中的差值Δ1；然后取多个阈值将所有的Δ1值分为多个段，将每一段中对应的元素的值用对应的阈值替换，得到裂缝图像I1；

步骤32：将步骤2得到的去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的垂直单相扫描，得到裂缝图像I2；具体是：将矩阵中每列的数据个数M等分得到M个数据段，M是能被每列元素的个数整除的数；对每个数据段中的元素进行曲线拟合，用拟合值减去其对应的元素得到垂直单相扫描中的差值Δ2；然后将所有的Δ2值分为多个段，将每一段中对应的元素的值用对应的阈值替换；

步骤33：将裂缝图像I1与裂缝图像I2取交集，得到交集裂缝图像I3；

步骤34：将裂缝图像I1与裂缝图像I2取并集，得到并集裂缝图像I4；

步骤35：将裂缝图像I1、裂缝图像I2、交集裂缝图像I3、并集裂缝图像I4分别划分为p*q个小块，p是能被每行元素的个数整除的数，q是能被每列元素的个数整除的数；

对于交集裂缝图像I3，从左到右从上而下扫描每个小块，对于每个小块，分别计算该小块在裂缝图像I1、裂缝图像I2以及交集裂缝图像I3中对应位置的小块中黑色像素点所占的比例R1、R2和R3；计算R3与R1的相对误差R31＝|R1-R3|/R1，若R1＝0，则R31＝0；计算R3与R2的相对误差R32＝|R2-R3|/R2，若R2＝0，则R32＝0；计算交集裂缝图像I3中该小块与裂缝图像I1、I2中对应小块的相似度R312＝0.5*R31+0.5*R32；

对于并集裂缝图像I4，从左到右从上而下扫描每一个小块，对于每个小块，分别计算其在裂缝图像I1、裂缝图像I2、并集裂缝图像I4中对应的小块中黑色像素点所占的比例R1、R2和R4，计算R4与R1的相对误差记为R41＝|R1-R4|/R1，若R1＝0，则R41＝0；计算R4与R2的相对误差R42＝|R2-R4|/R2，若R2＝0，则R42＝0；R412＝0.5*R41+0.5*R42；

步骤36：对于每个小块，比较其对应的R312与R412的大小，若R312大于等于R412，则取I4中的相应小块作为最终裂缝图像中的相应位置的小块；否则取I3中的相应小块作为最终裂缝图像中的相应位置的小块，得到最终裂缝图像I5；

步骤4：从最终裂缝图像中提取裂缝种子点并进行区域生长，得到裂缝二值化图像；

步骤5：将步骤4得到的裂缝二值化图像等分为多个正方形；

步骤6：计算路面图像中裂缝所占比例：统计包含有裂缝像素点的正方形个数和网格划分后的整幅图像中的正方形块个数，计算前者在后者中所占的比例即得到路面图像中裂缝所占比例。

2.如权利要求1所述的路面图像中裂缝比例检测方法，其特征在于，所述步骤1中所述的路面三维数据矩阵O_m×n如下：

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

z_ij表示行号为i，列号为j所对应的路面高度。

3.如权利要求1所述的路面图像中裂缝比例检测方法，其特征在于，所述步骤2对路面三维数据矩阵进行去噪处理的操作如下：

画出路面高度直方图，图中横坐标为路面三维数据矩阵中的高度数据，纵坐标为路面三维数据矩阵中的对应于每个高度数据段的元素个数；将路面高度直方图中两个高度数据段所对应的元素分别标记；其他高度数据段对应的元素标记为噪声点；对标记的噪声点进行滤波处理，得到去噪后的路面三维数据矩阵。

4.一种路面图像中裂缝比例检测系统，其特征在于，包括如下依次相连接的模块：

矩阵读取模块，是用于实现路面三维数据矩阵读取的模块；

去噪模块，是用于对路面三维数据矩阵进行去噪处理得到去噪后的路面三维数据矩阵的模块；

路面裂缝识别模块，是用于从去噪后的路面三维数据矩阵得到最终裂缝图像的模块；所述路面裂缝识别模块用于实现以下功能的模块：

对去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的水平单相扫描，得到裂缝图像I1；具体是：将矩阵中每一行的数据个数N等分得到N个数据段，N为能被每行元素的个数整除的数；对每个数据段中的元素进行曲线拟合，用拟合值减去对应元素的值得到水平单相扫描中的差值Δ1；然后取多个阈值将所有的Δ1值分为多个段，将每一段中对应的元素的值用对应的阈值替换，得到裂缝图像I1；

将去噪后的路面三维数据矩阵进行基于平均值法曲线拟合的垂直单相扫描，得到裂缝图像I2；具体是：将矩阵中每列的数据个数M等分得到M个数据段，M是能被每列元素的个数整除的数；对每个数据段中的元素进行曲线拟合，用拟合值减去其对应的元素得到垂直单相扫描中的差值Δ2；然后将所有的Δ2值分为多个段，将每一段中对应的元素的值用对应的阈值替换；

将裂缝图像I1与裂缝图像I2取交集，得到交集裂缝图像I3；

将裂缝图像I1与裂缝图像I2取并集，得到并集裂缝图像I4；

将裂缝图像I1、裂缝图像I2、交集裂缝图像I3、并集裂缝图像I4分别划分为p*q个小块，p是能被每行元素的个数整除的数，q是能被每列元素的个数整除的数；

对于交集裂缝图像I3，从左到右从上而下扫描每个小块，对于每个小块，分别计算该小块在裂缝图像I1、裂缝图像I2以及交集裂缝图像I3中对应位置的小块中黑色像素点所占的比例R1、R2和R3；计算R3与R1的相对误差R31＝|R1-R3|/R1，若R1＝0，则R31＝0；计算R3与R2的相对误差R32＝|R2-R3|/R2，若R2＝0，则R32＝0；计算交集裂缝图像I3中该小块与裂缝图像I1、I2中对应小块的相似度R312＝0.5*R31+0.5*R32；

对于并集裂缝图像I4，从左到右从上而下扫描每一个小块，对于每个小块，分别计算其在裂缝图像I1、裂缝图像I2、并集裂缝图像I4中对应的小块中黑色像素点所占的比例R1、R2和R4，计算R4与R1的相对误差记为R41＝|R1-R4|/R1，若R1＝0，则R41＝0；计算R4与R2的相对误差R42＝|R2-R4|/R2，若R2＝0，则R42＝0；R412＝0.5*R41+0.5*R42；

对于每个小块，比较其对应的R312与R412的大小，若R312大于等于R412，则取I4中的相应小块作为最终裂缝图像中的相应位置的小块；否则取I3中的相应小块作为最终裂缝图像中的相应位置的小块，得到最终裂缝图像I5；

裂缝二值化模块，是用于从最终裂缝图像中提取裂缝种子点并进行区域生长，得到裂缝二值化图像的模块；

图像网格划分模块，是用于将裂缝二值化图像等分为多个正方形的模块；

裂缝比例计算模块：是用于实现以下功能的模块：统计包含有裂缝像素点的正方形个数和网格划分后的整幅图像中的正方形块个数，计算前者在后者中所占的比例得到路面图像中裂缝所占比例。

5.如权利要求4所述的路面图像中裂缝比例检测系统，其特征在于，所述矩阵读取模块中所述的路面三维数据矩阵O_m×n如下：

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

z_ij表示行号为i，列号为j所对应的路面高度。

6.如权利要求4所述的路面图像中裂缝比例检测系统，其特征在于，所述路面裂缝识别模块是用于实现以下功能的模块：画出路面高度直方图，图中横坐标为路面三维数据矩阵中的高度数据，纵坐标为路面三维数据矩阵中的对应于每个高度数据段的元素个数；将路面高度直方图中两个高度数据段所对应的元素分别标记；其他高度数据段对应的元素标记为噪声点；对标记的噪声点进行滤波处理。