二階矩&灰度差分的橋梁裂縫快速識別方法

李 娜 許元飛 賈澎濤

(西安科技大學(xué)計算機科學(xué)與技術(shù)學(xué)院 陜西 西安 710054)

0 引 言

橋梁裂縫識別是數(shù)字圖像處理技術(shù)與橋梁監(jiān)測的主要結(jié)合點，圖像質(zhì)量評價成為該問題的有效解決方法。盡管已有裂縫識別算法眾多，但裂縫圖像的連通性和完整性方面均表現(xiàn)欠佳，快速、準(zhǔn)確地識別完整的裂縫圖像依然是一大挑戰(zhàn)[1]。近年來研究者趨向于人眼視覺特征的圖像質(zhì)量評價方法[2-3]研究，具體是將圖像在顏色空間中進(jìn)行分解，在已分解的子分量上利用已有人眼視覺模型獲取質(zhì)量測度，但應(yīng)用于復(fù)雜紋理圖像時，存在計算復(fù)雜的問題。對此，研究者提出了各種可計算的質(zhì)量測度方法。HDR[4]方法利用HDR-VDP-2.2中的基于視覺感知的模型得到關(guān)于亮度與對比度的視覺保真度特征，通過支持向量回歸對特征進(jìn)行融合，得到高動態(tài)圖像的客觀評價值。文獻(xiàn)[5]結(jié)合峰值信噪比(PSNR)的計算優(yōu)勢與人眼視覺特征信息熵的結(jié)構(gòu)特性，提出一種基于視覺特征信息熵的彩色圖像質(zhì)量評價方法。VIF[6]方法通過計算失真和參考圖像間的互信息相關(guān)性實現(xiàn)質(zhì)量評價，具有較好的性能，缺點是小波分解帶來復(fù)雜計算。

圖像質(zhì)量主觀評價和客觀評價的利弊體現(xiàn)在：主觀評價方法[7]準(zhǔn)確、合理，但受觀察者知識水平等影響，導(dǎo)致結(jié)果不穩(wěn)定，且難以滿足實時處理的要求[8-9]?？陀^評價具有實時穩(wěn)定、簡單高效等優(yōu)點?？陀^評價從自然圖像統(tǒng)計和空域變換角度出發(fā)，提出無參客觀質(zhì)量評價方法[10-12]。基于梯度信息與HVS濾波器[13]的方法建立了一種新的清晰度評價算法，該算法計算每個像素點的梯度值，得到每點的清晰度描述。Lab色彩空間更接近于人眼視覺，通過提取Lab色彩空間的信息熵特征來表征圖像質(zhì)量信息[14]。文獻(xiàn)[15]將圖像的模糊度評價轉(zhuǎn)化為待測圖像與參考圖像對應(yīng)分量的相似度測量。

為了提高識別的實時性與準(zhǔn)確性，本文在無災(zāi)害參考圖像和待檢測圖像間建立安全評價函數(shù)，利用評價值對橋梁裂縫的病害情況進(jìn)行自動評價并快速識別。

1 基于圖像質(zhì)量評價的橋梁裂縫識別算法

1.1 橋梁裂縫安全評價時間序列

選擇橋梁初建成時且通過相關(guān)質(zhì)量認(rèn)定即質(zhì)量最好時的圖像作為基礎(chǔ)參考圖像R。設(shè)定一定時間間隔，采集橋梁表面圖像。設(shè)當(dāng)前n(n>1)時刻的圖像為待評價圖像En，將其與參考圖像R利用評價函數(shù)進(jìn)行比較，得出質(zhì)量評價參數(shù)Dn，如圖1所示。

圖1 基于參考圖像的橋梁輔助安全評價

圖1產(chǎn)生兩組時間序列，分別是待評價圖像序列[E1,E2,…,En]、評價參數(shù)序列[D1,D2,…,Dn]。評價圖像En和參考圖像R分別通過相關(guān)度比較算法即安全評價函數(shù)F計算評價參數(shù)Dn，如式(1)所示。若Dn有較大異常，提示人工檢測。

Dn=F(En)-F(R)

(1)

1.2 基于顏色矩與灰度差分的安全評價

HSV顏色模型與其他顏色空間相比較，能夠更好地適應(yīng)人類對顏色的主觀理解，本文選擇HSV模型實現(xiàn)裂縫安全評價。

通過對圖像進(jìn)行H、S、V分量觀察，發(fā)現(xiàn)V分量的圖像與灰度圖視覺效果相同，式(1)變成：

Dn=F(Env)-F(Rv)

(2)

式(2)在HSV模型中的V分量上進(jìn)行裂縫安全性計算。

通過二階矩計算待測區(qū)域的顏色方差，即顏色的不均勻性。En與R的第i像素在V分量上的二階矩分別表示為：

(3)

式中：N為像素的個數(shù)，M1Eiv和M1Riv分別為Ei和Ri的一階矩。

(4)

V(Eiv)和V(Riv)分別代表Ei與Ri上第i個像素在V分量上的值。

裂縫在顏色特征變化的同時，紋理特征也有改變，通過熵計算紋理圖像的灰度差分統(tǒng)計特征。Ei與Ri的熵分別表示為：

(5)

利用顏色二階矩與灰度差分表示橋梁裂縫圖像在顏色和紋理上的變化，并進(jìn)行裂縫安全性比較，則F表示為：

F=M2+Ent

(6)

評價參數(shù)Dn表示為：

Dn=(M2Eiv-M2Riv)+(EntEiv-EntRiv)

(7)

基于顏色二階矩和灰度差分的橋梁裂縫病害評價流程如圖2所示。

圖2 安全評價算法流程

在圖2中，通過基于顏色二階矩與熵的人眼視覺特征相關(guān)度方法評價裂縫安全性，評價參數(shù)D值越大，表示相關(guān)度越低，代表安全質(zhì)量較差，病害程度越嚴(yán)重。

對公路監(jiān)控設(shè)備采集到的橋梁表面圖像信息，利用數(shù)字圖像的質(zhì)量評價，提出道路安全自動輔助評價方法。通過待評價圖像En與參考圖像R相關(guān)程度的客觀比較，利用得到的評價參數(shù)輔助說明橋梁表面裂縫的質(zhì)量安全。

2 實驗仿真

選取6組病害程度不同的橋梁裂縫圖像作為待檢測對象，資料搜集建立初期的圖像作為參考圖像R。

2.1 裂縫提取結(jié)果

為了與安全評價算法進(jìn)行對比，分別選取邊緣提取sobel算子、prewitt算子、roberts算子、log算子和canny算子對6組裂縫實現(xiàn)特征提取，觀察病害程度，具體見圖3-圖8。

圖3 橋梁裂縫1

圖4 橋梁裂縫2

圖5 橋梁裂縫3

圖6 橋梁裂縫4

圖7 橋梁裂縫5

圖8 橋梁裂縫6

對于6組病害圖，5種算子中前三種基本能夠提取出裂縫特征，后2種算子噪聲影響較大。通過算子提取的裂縫特征，觀察得到6組裂縫寬度由小到大變化。

2.2 安全評價參數(shù)比較

對圖3-圖8的6組病害圖用本文方法進(jìn)行安全評價。按照圖2的算法流程，通過式(3)、式(5)、式(7)計算質(zhì)量評價參數(shù)，具體數(shù)據(jù)見表1。

表1 6組病害圖的安全評價參數(shù)

表1中，均值、方差不符合裂縫提取6組病害由小到大的變化趨勢，數(shù)值有前后波動，無法獲得變化規(guī)律?；叶炔罘衷诘?組圖處有波動，評價參數(shù)D的值由小到大變化，結(jié)合1.2節(jié)對D的說明，評價參數(shù)檢測到6組病害逐級嚴(yán)重，與2.1節(jié)的特征提取結(jié)果一致。

將本文方法與峰值信噪比方法和文獻(xiàn)[7]中彩色圖像質(zhì)量評價方法計算的圖像質(zhì)量相關(guān)度參數(shù)進(jìn)行比較，計算結(jié)果見表2，值越大代表相關(guān)度越低。

表2 幾種方法的質(zhì)量評價參數(shù)

2.3 實驗分析

根據(jù)2.1節(jié)中的裂縫提取結(jié)果，將病害圖按照安全質(zhì)量優(yōu)劣順序排序，依次是橋梁裂縫1、橋梁裂縫2、橋梁裂縫3、橋梁裂縫4、橋梁裂縫5、橋梁裂縫6。評價參數(shù)序列依照值越小相關(guān)度越高，安全質(zhì)量較好的原則。本文方法的客觀質(zhì)量評價從小到大依次是橋梁裂縫1、橋梁裂縫2、橋梁裂縫3、橋梁裂縫4、橋梁裂縫5、橋梁裂縫6，結(jié)果與裂縫提取判斷一致，驗證了方法的有效性。

而對于表2中峰值信噪比的評價參數(shù)值，橋梁裂縫1、橋梁裂縫3、橋梁裂縫4由于受到噪聲的影響，值較小，不易按照最小值判斷其裂縫情況，故不適合于橋梁裂縫的自動安全評價。對于視覺感知的評價參數(shù)值，能夠說明橋梁裂縫6的相關(guān)度最低，但是該方法的時間復(fù)雜度高于本文方法。

3 結(jié) 語

針對橋梁監(jiān)控圖像的計算機自動安全輔助，本文提出基于參考圖像的輔助安全識別方法。建立具體的顏色二階矩與人眼視覺信息熵特征客觀評價函數(shù)及算法流程。利用參考圖像序列、評價圖像序列，計算評價參數(shù)序列，本文方法自動化程度較高，且與人眼主觀感受一致，對于橋梁的自動安全監(jiān)控有重要的現(xiàn)實應(yīng)用價值。由于監(jiān)控圖像的持續(xù)性難度和不穩(wěn)定性特點，本文方法在特殊環(huán)境下有可能造成評價失誤。下一步工作主要是對評價參數(shù)矩陣的特點進(jìn)行分析，研究其中蘊含的質(zhì)量評價規(guī)律。