宋蘭平，田長(zhǎng)超，魏 敏

（山東省交通科學(xué)研究院，山東 濟(jì)南 250031）

引言

路面損壞狀況指數(shù)的大小可以直觀地反映道路表面的破損情況，其中路面裂縫的檢測(cè)是其計(jì)算的一項(xiàng)重要依據(jù)[1-2]。傳統(tǒng)的人工檢測(cè)方式檢測(cè)效率低、投入成本大。SUBIRATS等[3]將小波變換方法引入裂縫檢測(cè)中，同時(shí)檢測(cè)提取裂縫信息。李剛和賀昱曜[4]提出使用Sobel算子，結(jié)合最大熵圖像分割算法處理不同光照下的裂縫圖像，并提取裂縫信息。張磊[5]利用傳統(tǒng)數(shù)字圖像處理方法檢測(cè)裂縫。

1 基本方法

1.1 圖像預(yù)處理

市政道路現(xiàn)場(chǎng)采集路面裂縫樣本一般由紅色、綠色、藍(lán)色三個(gè)通道組成，即RGB模型，是一個(gè)三維矩陣。為了使數(shù)字圖像在操作過程中更加方便準(zhǔn)確，將RGB圖像轉(zhuǎn)化為只有一個(gè)采樣顏色圖像，即灰度圖像?；叶葓D像通常顯示為從最暗的黑色到最亮的白色灰度，理論上這個(gè)采樣可以是任何顏色的不同深淺，甚至可以是不同亮度上的不同顏色。

1.2 中值反差濾波方法

不同的濾波窗在用于圖像處理時(shí)有著不同的作用，中值濾波窗可以去除一些獨(dú)立的散點(diǎn)，均值濾波窗能將圖像中的空缺點(diǎn)進(jìn)行填補(bǔ)。從路面病害圖樣中可以看出，裂縫在整幅圖像中具有異于路面特征的形態(tài)，同時(shí)由于一些光照及路面不清潔的影響，會(huì)產(chǎn)生離散的噪點(diǎn)信息，這些信息會(huì)對(duì)裂縫信息的整體提取工作起到誤導(dǎo)作用。為準(zhǔn)確提取路面裂縫的分布信息，以中值濾波窗為基礎(chǔ)，設(shè)計(jì)中值反差濾波器：

式中：X—路面病害采集圖樣數(shù)組；Xsize—經(jīng)過中值反差濾波后的圖樣數(shù)組；in,out,dc—濾波器的尺寸變量，其定義見圖1。

圖1 濾波器輸入、輸出和忽略區(qū)域分布

濾波器只會(huì)保留幅度值在2in與2dc之間對(duì)數(shù)據(jù)，處于以內(nèi)或者以外的這部分?jǐn)?shù)據(jù)都將被濾波器濾除。

1.3 連通域方法

路面裂縫在空間分布上具有一定的連續(xù)性和連通性，但經(jīng)過濾波器處理的圖像噪點(diǎn)分布噪聲往往是孤立的[6]，在空間上是不連續(xù)分布的、雜亂無章的。因此可采用連通域標(biāo)記法，利用連通域的面積特征約束來去除噪聲。同時(shí)采用形態(tài)學(xué)的膨脹和腐蝕[7]運(yùn)算操作來修復(fù)斷裂的裂縫。

在一幅M×N的灰度圖像中，符號(hào)b（x，y）表示坐標(biāo)值為（x，y），（0≤x≤M-1，0≤y≤N-1）的像素，像素b（x，y）在水平和垂直方向上分別有一個(gè)鄰接像素，分布基本關(guān)系見圖2。像素b（x-1，y），b（x，y-1），b（x+1，y），b（x，y+1）稱為像素b（x，y）的4-鄰域像素，見圖3。

圖2 像素分布的基本關(guān)系

圖3 像素4-鄰域像素分布

圖5 現(xiàn)場(chǎng)采集到的路面裂縫原始

圖6 轉(zhuǎn)換為灰度圖像的路面裂縫

2 路面裂縫識(shí)別處理

2.1 圖像色彩模型轉(zhuǎn)化

利用圖像預(yù)處理中介紹的RGB模型與灰度圖像之間的轉(zhuǎn)換關(guān)系，將讀取的RGB圖像轉(zhuǎn)化為灰度圖像，見圖5、圖6。

2.2 圖像的濾波處理

獲得灰度圖像后，對(duì)圖像中的像素信息進(jìn)行降噪處理，主要是因?yàn)椴杉降牧芽p圖像在接收和處理過程中會(huì)受到不同程度的噪聲干擾，噪聲強(qiáng)度大的圖像質(zhì)量會(huì)下降，有的甚至?xí)⒂杏玫男畔⑻卣鞲采w掉。設(shè)置濾波器的尺寸為in=3，dc=7，out=12，對(duì)圖像從左到右，從上到下進(jìn)行過濾處理。

目前對(duì)圖像降噪的方法主要有均值模塊濾波、中值濾波法、高斯濾波法以及Sobel梯度算子濾波法。圖7為利用不同的濾波方法對(duì)轉(zhuǎn)化為灰度圖像的裂縫圖像進(jìn)行濾波處理。

可以看出：（1）均值濾波處理裂縫邊界較模糊；（2）中值濾波平滑效果好一些，但不能區(qū)分處理裂縫部分和樣本其余部分特征；（3）高斯濾波為線性濾波器，能較好地抑制噪聲保留圖像的細(xì)節(jié)部分，但在圖像處理過程中無法對(duì)裂縫部分做專門識(shí)別處理；（4）梯度算子濾波主要是根據(jù)子模的分布形式進(jìn)行遍歷，Sobel梯度算子濾波后，許多存在梯度差值的敵方都會(huì)將邊界信息保留，會(huì)給后期處理帶來復(fù)雜的處理過程；（5） 中值反差濾波處理后的圖像，保留了裂縫的主要信息，并且其余地方的背景分布被濾波器弱化，從而增強(qiáng)了裂縫位置的信息，使得裂縫更加明顯。

圖7 經(jīng)過濾波處理后的圖像

中值反差濾波器不僅很好地保留了裂縫信息，且能弱化其余的噪點(diǎn)信息，濾波器在圖像遍歷過程中，除了對(duì)裂縫形成的區(qū)域增強(qiáng)外，還會(huì)將其余的噪點(diǎn)（與裂縫像素相近）形成孤立區(qū)域，連通域標(biāo)記方法可以將此類噪點(diǎn)標(biāo)識(shí)出來。

2.3 圖像的連通域標(biāo)記

結(jié)合預(yù)設(shè)好的4-鄰域連通域方法對(duì)濾波處理后的圖像進(jìn)行處理，并利用形態(tài)學(xué)的膨脹、腐蝕[8]操作進(jìn)行圖像內(nèi)部的填充及邊緣的細(xì)化。

連通域標(biāo)記法是首先逐行掃描一幅圖像，以單素為標(biāo)記對(duì)象，逐一處理圖片中的每個(gè)像素。當(dāng)一個(gè)待標(biāo)記像素b（x，y）是背景時(shí)，無需做任何操作，遍歷下一個(gè)像素。當(dāng)一個(gè)待標(biāo)記像素b（x，y）是前景時(shí)，則按識(shí)別處理步驟進(jìn)行：（1）步驟1：如果b（x，y）是一個(gè)未被標(biāo)記的像素，并且它上行鄰接像素b（x，y-1）是一個(gè)背景像素，則b（x，y）一定是個(gè)新的未標(biāo)記連通域外部邊界像素，如圖8中P位置，需用新的標(biāo)號(hào)值標(biāo)記像素，然后利用邊界跟蹤找到所有的像素，并用-1標(biāo)記圍繞外部邊界的所有背景像素。（2）步驟2：如果b（x，y）的正下方的鄰接像素b（x，y+1）是一個(gè)背景像素，如圖9中的像素q，而左邊相近的像素b（x-1，y）是一個(gè)已標(biāo)記的像素，如圖9中的p像素，則b（x，y）是內(nèi)部邊界像素。用像素b（x，y）的標(biāo)號(hào)標(biāo)記當(dāng)前像素b（x，y）。由像素b（x，y）出發(fā)查找內(nèi)部邊界像素，并且用b（x，y）標(biāo)號(hào)標(biāo)記，用-1標(biāo)記圍繞邊界的所有背景像素。（3）步驟3：如果像素b（x，y）不屬于上述步驟1、2所描述的情況，那么像素點(diǎn)的左側(cè)像素b（x-1，y）會(huì)是個(gè)已標(biāo)記的像素。則用b（x-1，y）的標(biāo)記號(hào)來標(biāo)記像素b（x，y），并終止當(dāng)前像素的標(biāo)記處理。

圖8 連通域標(biāo)記過程中的外部邊

圖9 連通域標(biāo)記過程中的內(nèi)部邊緣追蹤

由連通域標(biāo)記后的圖像，會(huì)得到許多獨(dú)立的區(qū)域，共計(jì)3 782組獨(dú)立的分布區(qū)域，其中裂縫的區(qū)域同一標(biāo)號(hào)的像素點(diǎn)最多，所以只要將標(biāo)記號(hào)相同、數(shù)量最多的像素位置分布找到，就可確定路面裂縫的分布。圖10、圖11分別為提取后的裂縫骨架和將提取的骨架信息與原始圖疊合后的示意圖，從疊合后的圖像可以看出，本文算法所提取的裂縫骨架描跡信息與原始圖中的裂縫信息吻合較好。

圖10 提取后的裂縫結(jié)果

圖11 裂縫結(jié)果畫在原始圖上

2.4 試驗(yàn)處理

利用市政道路現(xiàn)場(chǎng)采集的106幅路面裂縫圖片，以相同的樣本集作試驗(yàn)對(duì)比，分別采用較流行的裂縫提取方法：動(dòng)態(tài)閾值+小波變換[8](算法1)和算法（算法2）進(jìn)行比較，對(duì)比結(jié)果見表1。

表1 算法對(duì)比數(shù)據(jù)

可以看出，算法2在各指標(biāo)上都具有一定的優(yōu)勢(shì)，提高了裂縫的提取精度，可以滿足施工檢測(cè)的實(shí)際需求；降低了裂縫處理時(shí)間，基本滿足工程實(shí)時(shí)性要求。

3 結(jié)語(yǔ)

（1）闡述了一種基于新型濾波器，融合連通域標(biāo)記等圖像處理的方法進(jìn)行路面裂縫識(shí)別，并且詳細(xì)地描述了濾波器的構(gòu)建過程。（2）結(jié)合路面裂縫的連通性特點(diǎn)，采用形態(tài)學(xué)中的連通域標(biāo)記法、膨脹和腐蝕運(yùn)算等操作，完成裂縫的修復(fù)和提取。（3）相較于傳統(tǒng)處理算法，無論在精度上還是在處理時(shí)間上，本文算法都具有一定的優(yōu)勢(shì)。

現(xiàn)場(chǎng)道路病害樣本的采集，受采集現(xiàn)場(chǎng)的光照、油污及雜物等外在因素影響，采集到的路面裂縫樣本具有數(shù)據(jù)量大、類型多樣的特點(diǎn)，所以對(duì)算法的魯棒性提出更大的要求，因此提高算法的魯棒性是下一步的工作方向。