謝君利,郭魏源,張春生
(1. 黃河科技學院應用技術學院,河南 濟源 459000;2. 河南理工大學土木工程學院,河南 焦作 454003)
我國交通運輸業的迅速發展使公路橋梁被大量建造并投入使用[1],鋼橋是公路橋梁的典型形式之一,具有受力明確、跨越能力大等優點,但在長期降雨和風沙等自然環境侵蝕下,鋼橋難免出現劣變情況,橋面裂縫是最常見的鋼橋劣變問題,裂縫不但會降低鋼橋的路用性能,而且對鋼橋受力也存在很大的負面影響,導致鋼橋整體性能下降[2]。因此,建立科學有效的鋼橋橋面圖像細小裂縫動態識別方法至關重要。
胡成雪等人[3]利用灰度校正和形態學重建降低外界干擾和噪聲影響,通過顯著性差異提取圖像鄰域顯著特征,引入方向可調濾波器獲取方向各異的梯度顯著特征,融合鄰域顯著特征和梯度顯著特征生成特征融合顯著圖并采用閾值分割,得到可能的裂縫聚集區域,采用聚類算法依據區域內不同幾何特征篩選裂縫候選區域,定位候選區域端點,剔除無端點區域子集,將不同區域端點連接,實現鋼橋橋面圖像細小裂縫動態識別。?;莸热薣4]混合傳統卷積塊和空洞卷積塊,構建金字塔空洞卷積模塊作為殘差網絡底部模塊,通過改進的殘差網絡實現鋼橋橋面圖像細小裂縫動態識別。謝雄耀等人[5]采用基于深度學習算法的Faster R-CNN網絡識別原始鋼橋橋面圖像裂縫目標,判定裂縫存在與否,若存在,標記裂縫區域并裁剪裂縫區域,進而得到全部裂縫圖像,去除圖像干擾背景后通過U-Net語義分割網絡實現鋼橋橋面圖像細小裂縫動態識別。以上方法沒有在識別到細小裂縫后對裂縫優化處理,導致誤差較大。
為了解決上述方法中存在的問題,提出鋼橋橋面圖像細小裂縫最優化動態識別方法。通過Mask勻光法和灰度級腐蝕法,提高圖像質量。采用改進Canny算子,提取細小裂縫邊緣。通過形態學方法實現裂縫最優化動態識別。
鋼橋橋面圖像中往往存在由橋面顆粒紋理、積水水漬、拋灑物、光照不均、對比度較低等因素造成的圖像質量低問題,在識別鋼橋橋面圖像細小裂縫前務必對圖像預處理。
2.1.1 圖像勻光處理
采用Mask勻光法平衡鋼橋圖像整體光強。
用Ig(x,y)表示原始圖像,Ig(x,y)表示高斯濾波法濾波后得到的背景圖像[6],Ib(x,y)表示線性拉伸前圖像,(x,y)表示圖像中任意像素,φ表示灰度偏移量,則存在以下關系式:
Ir(x,y)=Ig(x,y)-Ib(x,y)+φ
(1)
在[0,255]范圍內對Ib(x,y)的像素值拉伸處理,用ε表示線性拉伸因子,得到勻光后圖像Im(x,y)如下所示:
Im(x,y)=255×[Ir(x,y)-ε]/(255-2×ε)
(2)
2.1.2 圖像對比度增強處理
由于拍攝環境和拍攝設備等因素影響,鋼橋橋面圖像對比度較低,因此采用灰度級腐蝕法增強圖像對比度。
用H(x′,y′)表示結構元素[7],H(x′,y′)表示結構元素中任意像素,Gs表示H(x′,y′)的定義域,得到經灰度腐蝕增強對比度后圖像Ie(x,y)如下所示:
Ie(x,y)=min[Im(x+x′,y+y′)-H(x′,y′)]
(3)
在經過圖像預處理得到高質量鋼橋橋面圖像后,引入模糊C均值聚類算法(FCM)對圖像細小裂縫粗提取[8]。
(4)
FCM算法主要流程如下:
①選取大于0的一個小數作為停止閾值θ,初始化隸屬度矩陣記作U0,迭代次數t設置為0;
因為在鋼橋橋面圖像中光線被細小裂縫所吸收,所以裂縫所處區域反射作用較弱,即對應圖像區域像素灰度值較低。將灰度值作為聚類特征,記作xi,將聚類中心zj的最小值對應索引記作j=J,若第i個像素點對于第J個聚類中心的隸屬度μiJ在{μij|1≤j≤c}中為最大值,則將xi判定為裂縫。完成對全部像素點的判斷,即可實現鋼橋橋面圖像細小裂縫的粗提取。
鋼橋橋面圖像細小裂縫可視為不規則曲線段組合,存在一定的邊緣性質,基于粗識別結果提取細小裂縫邊緣有助于后續對細小裂縫的分析和識別,采用改進Canny算子提取鋼橋橋面圖像細小裂縫邊緣[9]。
傳統Canny算子的濾波對非線性噪聲沒有濾波效果,濾波能力受標準差σ的直接影響,其一階差分會造成部分圖像信息的丟失[10],并且其采用人工定義閾值,會導致結果出現一定偏差。因此所提方法對傳統Canny算子加以改進。
采用非線性中值濾波替代高斯濾波作為濾波算法,中值濾波依據升序排列鋼橋橋面圖像像素值并選取中間值實現濾波[11]。
引入Sobel算子的8方向模板計算圖像梯度[12],假設存在鋼橋橋面圖像f(x,y),(x,y)表示圖像中任意像素值,則經Sobel算子計算得到方向梯度Sx和Sy如下所示:
(5)
通過Otsu法替代人工閾值選取[13],假設鋼橋橋面圖像中共存在Y個灰度級,每級灰度出現次數記作mi,對于一幅尺寸為M×N的鋼橋橋面圖像,每級灰度出現概率pi計算方式如下所示:
pi=mi/(M×N),i∈[0,1,…,Y-1]
(6)
設置閾值K用于劃分圖像背景和邊緣區域,基于該閾值計算背景和邊緣區域在圖像中的概率Pb和Po,如下所示:
(7)
(8)
(9)
在[0,Y-1]范圍內不斷調節閾值K,重復計算不同閾值下的類間方差,其中,類間方差最大值對應的閾值即為所求閾值,通過該閾值即可區分鋼橋橋面圖像細小裂縫邊緣區域與背景區域,實現細小裂縫邊緣提取。
區域生長法通過合并相鄰且具有近似屬性的像素點為一個區域實現細小裂縫識別[14],依據灰度閾值制定生長規則,將生長點灰度值記作p,生長點8鄰域內像素點灰度值記作qi,i∈[1,2,…,8],生長閾值記作Td,則有:
|p-qi|
(10)
通過區域生長法動態識別鋼橋橋面圖像細小裂縫的主要流程如下所示:
①在粗識別區域內,將通過改進Canny算子提取到的鋼橋橋面圖像細小裂縫邊緣點設定為種子點集合;
②設定種子點生長規則,根據式(10)確定生長閾值;
③以種子點為起始生長點,標記種子點所在區域為目標區域Q;
④對生長點8鄰域內全部像素點加以搜索,保留不屬于Q的像素點和未超出粗識別區域像素點;
⑤對比保留像素點和相應生長點,將滿足生長規則的像素點記作新生長點,添加至Q中;
⑥若新生長點個數0,則算法結束,最終得到的區域Q即為鋼橋橋面圖像細小裂縫區域;反之返回步驟④重復執行直到滿足終止要求。
由于鋼橋橋面圖像中存在無法避免的噪聲,識別到的細小裂縫可能會出現不連續、空洞或毛刺等問題,因此對識別到的細小裂縫進一步處理,通過形態學方法去除細小裂縫圖像毛刺和空洞[15],依據像素位置信息拼接細小裂縫。
引入形態學方法中的開操作和閉操作處理細小裂縫圖像,前者負責去除比結構元素更小的毛刺,后者負責填補比結構元素更小的空洞或缺口,進而使得細小裂縫邊緣更為平滑。
識別到的細小裂縫圖像中可能存在不連續的問題,同一條細小裂縫斷開區域間距離較近,所提方法基于該特征拼接細小裂縫,主要流程如下所示:
①將通過區域生長法得到的細小裂縫均標記為候選裂縫,由左上至右下排列全部裂縫,生成候選裂縫集合,選取首個裂縫作為當前待拼接裂縫;
②由該待拼接裂縫端點開始判斷裂縫走向,遍歷該裂縫端點方向附近像素點,如果在端點指定像素點距離內存在其它候選裂縫,則將該裂縫與待拼接裂縫拼接;
③完成當前裂縫拼接后選取下一個候選裂縫作為待拼接裂縫執行步驟②,直到完成全部候選裂縫拼接。
所提方法具有計算過程簡單、效率高等特點,采用所提方法對不斷更新的鋼橋橋面圖像持續識別、監測,即可實現對鋼橋橋面圖像細小裂縫的實時最優化動態識別。
為了驗證鋼橋橋面圖像細小裂縫最優化動態識別方法(所提方法)的整體有效性,需要對其測試。實驗選取362幅鋼橋橋面細小裂縫圖像,構建裂縫數據庫。采用鄰域與梯度顯著特征融合的裂縫檢測方法(參考文獻[3]方法)和基于改進殘差網絡的鐵路隧道裂縫檢測算法(參考文獻[4]方法)為對照檢測所提方法的有效性。任意選取4幅圖像對比所提方法、參考文獻[3]方法和參考文獻[4]方法的細小裂縫識別結果,測試結果如圖1所示。
圖1 三種方法細小裂縫識別結果視覺對比圖
由圖1可以看出,所提方法的細小裂縫識別結果更為清晰,邊緣更為明確,而參考文獻[3]方法和參考文獻[4]方法雖然能夠識別到細小裂縫,但與背景區別較小,難以區分,識別效果不理想。
為了更加客觀的評估所提方法、參考文獻[3]方法和參考文獻[4]方法的檢測結果,采用ROC曲線、AUC值、平均絕對誤差量化362幅鋼橋橋面圖像細小裂縫檢測結果,測試結果如圖3所示。
1)ROC曲線和AUC值
ROC曲線圖能夠用來衡量裂縫識別的準確性,ROC曲線與橫縱坐標圍成的面積越大,即AUC值越大,則對應方法的裂縫識別準確性越高,不同方法的檢測結果如圖2所示:
圖2 ROC曲線和AUC值檢測結果
由圖2可以看出,在相同假陽性率處,所提方法真陽性率高于參考文獻[3]方法和參考文獻[4]方法,說明所提方法檢測出的點是真實裂縫點占比最高,同時所提方法的AUC值高于參考文獻[3]方法和參考文獻[4]方法,因為所提方法在識別到鋼橋橋面細小裂縫后采用形態學方法和像素位置信息拼接方法解決細小裂縫不連續、空洞或毛刺等問題,優化識別效果。
2)平均絕對誤差
ROC曲線主要衡量的是潛在識別目標,為了全面評估裂縫識別性能,引入平均絕對誤差進一步檢測所提方法、參考文獻[3]方法和參考文獻[4]方法的可靠性與準確性。
用M和N表示圖像長和寬,(x,y)表示任意像素值,G(x,y)表示真實像素值,H(x,y)表示不同方法檢測下像素值,則平均絕對誤差MAE計算方法如下所示:
(11)
通過式(11)計算三種方法在不同鋼橋橋面圖像細小裂縫類別中識別結果的平均絕對誤差,如下表所示。
表1 平均絕對誤差檢測結果
由圖3可以看出,所提方法在不同裂縫類別中的平均絕對誤差始終保持在0.015以下,參考文獻[3]方法的平均絕對誤差均已超過在0.02,參考文獻[4]方法的平均絕對誤差總體趨勢雖略低于參考文獻[3]方法,但明顯高于所提方法,說明所提方法具有一定的有效性和優越性。
圖3 不同方法的運行時間檢測結果
為了進一步檢測三種方法在細小裂縫識別中的性能,統計三種方法在識別362幅鋼橋橋面圖像細小裂縫時的平均運行時間,結果如圖3所示。
由圖3可以看出,所提方法的運行時間在0.5s以下,參考文獻[3]方法運行時間在1.5s以上,參考文獻[4]方法運行時間接近3s,所提方法的運行時間明顯少于參考文獻[3]方法和參考文獻[4]方法,說明在相同時間內所提方法能夠完成更多的裂縫識別任務,更符合動態識別需求。
為了保障鋼橋使用壽命和運行安全,對鋼橋橋面裂縫及時檢測識別并維護至關重要。為此,提出鋼橋橋面圖像細小裂縫最優化動態識別方法。采用Mask勻光法和灰度級腐蝕法預處理圖像,利用改進的Canny算子提取細小裂縫邊緣,通過區域生長法識別細小裂縫并結合形態學方法和像素位置信息拼接方法緩解細小裂縫不連續、空洞或毛刺等問題,不斷分析、識別持續更新的鋼橋橋面圖像細小裂縫信息,完成細小裂縫最優化動態識別。該方法能夠有效地增大AUC值、減小平均絕對誤差、縮短運行時間,為鋼橋的安全運行奠定基礎。