【能源人都在看，點擊右上角加“關注”】
北極星環境監測網訊:摘要：圖像在形成、傳輸、接收和處理的過程中，不可避免地受到噪聲的影響，從而影響了圖像的視覺效果及后續的圖像處理過程。本文詳細論述了基于噪聲檢測的圖像去噪算法。
圖像去噪是圖像處理領域中一項基本，而又十分關鍵的技術，一直是圖像處理領域的一個難題。在圖像的獲取、傳輸和存貯的過程中不可避免地受到各種噪聲源的干擾。為了從圖像中獲取更準確的信息，圖像去噪預處理算法的好壞成為后續處理的關鍵。
一、噪音簡介
噪聲一般被認為是多余的干擾，所以圖像中的噪聲可定義為圖像中不需要的干擾部分。若從信號的角度來描述圖像，那么圖像是一種信號。同樣噪聲從廣義的角度也可被看作是一種信號。只不過這種信號所攜帶的信息是無意義的，或者說沒有攜帶“信息”。因此噪聲和圖像中的有用信號之間既對立又統一。在研究噪聲的路徑上，既可利用噪聲與信號的統一關系，通過信號的特性來對噪聲的特性進行推理和預測；也可從二者的對立出發，將噪聲看作與信號無關的統計量，進而單獨進行研究。
圖像中的噪聲主要來源于兩個過程，一是圖像的獲取過程，二是圖像的傳輸過程。圖像獲取過程中，由于成像傳感器自身質量的影響及環境條件的作用，產生圖像退化，從而使圖像質量下降。而圖像在傳輸過程中，傳輸信道中的干擾是造成圖像污染主要因素。另外，噪聲對圖像的破壞作用主要表現在以下方面：①影響主觀視覺感受，降低了人們觀看、理解和欣賞圖像的能力。②使圖像的后續中高層研究及應用無法進行或得出錯誤結論。噪聲降低了數據源（原始圖像）的質量和精度，產生局部二義性，最直接的影響就是對邊緣檢測和圖像分割產生誤導。
二、圖像中去噪方法概述
圖像去噪的目的是為了減少圖像噪聲。圖像噪聲來自于多方面，它由來自于系統外部干擾,如電磁波和經電源串進系統內部而引起的外部噪聲，也有來自于系統內部的干擾,如攝像機的熱噪聲,電器機械運動而產生的抖動噪聲等內部噪聲。減少噪聲的方法可在圖像空間域或在圖像變換域完成。
圖像空間域去噪方法較多，鄰域平均法是一種局部空間域處理的算法。圖像經過鄰域平均法處理后會變得相對模糊,這是因平均本來是以圖像的模糊為代價來換取噪聲的減少。空間域低通濾波方法也可平滑圖像的噪聲,它實際上是通過一個低通卷積模板在圖像空間域進行二維卷積來達到去除圖像噪聲的目的。多幅圖像平均法是利用對同一景物的多幅圖像取平均來消除噪聲的。中值濾波是一種空間域非線性濾波方法,由于它在實際運算過程中并不需要圖像的統計特性,所以比較方便。在一定的條件下,它能克服線性濾波器所帶來的圖像細節模糊,而且對濾除脈沖干擾及圖像噪聲最為有效。
圖像變換域去噪方法是對圖像進行某種變換,將圖像從空間域轉換到變換域,再對變換域中的變換系數進行處理,再進行反變換將圖像從變換域轉換到空間域來達到去除圖像噪聲的目的。將圖像從空間轉換到變換域的變換方法很多,如付利葉變換、小波變換等。每種變換它的變換域得到的系數都有不同的特點,合理地處理變換系數再通過反變換將圖像還原到空間域，能有效達到去除噪聲的目的。
三、新的濾波器設計
設濾波器窗口大小為L(L為奇數)，將濾波器窗口內的原始圖像的樣本像素從上到下、從左到右排列，寫成
(1)
為了方便，記這4個方向核為Kh(h=1，2，3，4)，由于沖擊噪聲一般都是基于下面兩個假設的：1)沒有受到沖擊噪聲干擾的圖像區域是局部平滑的，這些局部平滑區域由邊緣分隔開；2)受到沖擊噪聲污染的像素的灰度值與它臨近沒有受到噪聲污染的像素的平均灰度值的差異較大。所以，利用上面的4個方向算子，設計如下的公式去判斷是否為噪聲像素并實現本文提出的濾波算法，即
基于相關特征約束的分水嶺路面裂縫圖像分割算法
摘要：針對路面裂縫圖像的分割問題，提出一種相關特征約束的分水嶺分割算法。該算法基于形態學分水嶺思想，為了解決分水嶺過分割問題，分析了裂縫圖像狹長山谷的特點，選擇從梯度圖像的高頻成分中提取與裂縫特征信息相關的局部極小值，從而約束偽極小值的出現，再對原始梯度圖像強制標記，消除了過分割現象。同時解決了在消除過分割時邊緣信息被濾除的問題，保留了細微有用的邊緣。實驗表明，該算法分割出的路面裂縫圖像準確清晰，有效降低了過分割現象。
0引言
路面裂縫檢測是高速公路及時維護的有效工作手段，利用計算機進行路面裂縫圖像處理技術近年來取得了很大的進展。由于拍攝到的路面裂縫圖像受自然光照陰影及各種噪聲的影響，再加上裂縫信息雜亂微弱的特點，給分割提取工作帶來極大困難［1］。如何準確分割出裂縫信息，保證能夠及時維護路面是目前路政系統亟需解決的一個難點問題。利用各種改進的邊緣檢測算子對裂縫信息進行檢測，因常常伴隨有大量的噪聲信息，結果很不理想。而隨后出現的結合各種工具進行裂縫信息處理的方法，包括小波變換、人工種群、模糊邏輯等，這些算法對檢測結果有一定改觀，但實現復雜、運算量過大，尤其對一些網狀或輕微路面裂縫圖像仍然不能得到理想效果［2］。
分水嶺算法是數學形態學中的經典分割算法，具有直觀、簡便的特點，而且該算法以分割的邊緣輪廓單像素寬、連續封閉的優點，得到很多學者的青睞［3］。多種實驗證明，分水嶺分割對微弱的信息能夠有良好的響應，非常有利于路面裂縫圖像的分割。但分水嶺的過分割問題一直是分割中要處理的一個難點［4］。
因此，本文在利用分水嶺分割路面裂縫圖像時，為了解決分水嶺的過分割問題，分析了傳統方法在防止偽極小值過多的同時要保留邊緣細節信息的矛盾過程，提出了一種相關特征約束的分水嶺分割方法，該方法有效地解決了過分割現象，不僅檢測出了細微的裂縫信息，也避免了大量雜亂噪聲分支。
1分水嶺算法
1.1分水嶺算法原理
分水嶺算法原理是受地形學的啟發，模擬了地貌浸水的過程。把灰度圖像看作是一片地形，如圖1所示［5］，圖像中的每個像素點的值代表海拔高度，局部最小值表示“谷底”，最大值表示“山峰”。假設將地形放入水中浸沒，每個谷底被打洞，水位將快速上升，淹沒整個地形，處在每個谷底的水將因山峰的邊界分割開，形成分水嶺，從而達到分割圖像的目的。Vincent［6］提出的浸沒法分水嶺分割，是最經典的分水嶺算法，該算法由排序過程和淹沒過程完成，使分割速率大大提高。
1.2解決過分割現象
分水嶺算法可以得到單像素的邊緣輪廓，而且大多數都是封閉的，但其對細節信息非常敏感，雖然能識別出細微邊緣，但同時受噪聲影響較大，檢測到的區域極小值過多導致分割區域過多，從而出現過分割問題，這對后續的研究工作很不利。因此，很多研究人員從各種角度提出了若干解決過分割問題的方法，主要的解決方法大致有兩種：
（1）預處理
分割之前對原圖進行去噪濾波，將不必要的信息從原圖中先去除掉，使后面的分割過程中不再產生大量的偽極小區域，從而消弱過分割現象。但存在的問題是去噪的同時細小的有用邊緣信息也被濾除掉了［7］。
（2）分割后區域合并
這種方法用在分水嶺分割之后，采用合適的準則對分割后有相似性質的區域合并達到消除過分割的目的［8］。區域合并法往往因為在分割后出現的區域過多，導致合并處理運算量很大，時間復雜度較高。而且合并準則的選擇是否最優很難確定，目前在過分割問題上還需要繼續深入研究。
2改進的分水嶺分割算法
2.1改進算法的基本思想
根據以上分析，從大量的實驗結果來看，濾波后的圖像大量的目標邊緣信息被破壞。本文提出的改進的分水嶺分割算法不進行預處理，直接在原始梯度圖像上進行分水嶺分割，從而保留了有用的細節信息。
同時，依據經典的分水嶺思想，在淹沒過程中采用一種新的標記方法。首先根據圖像光照模型可知，圖像中的高頻成分對應著灰度變化激烈的地方，比如邊緣、噪聲；而低頻成分對應著灰度變化不大的平坦區域。從所要研究的路面裂縫圖像來看，要分割出的裂縫目標是灰度變化劇烈的地方，其對應的是圖像的高頻部分，因此，僅在高頻成分中提取與裂縫目標有強烈相關性的局部極小值，即可以保證極小值范圍被縮小，減少過分割區域。且因為沒有濾噪，其細微的裂縫信息不會丟失。然后，使用相關特征約束提取的局部極小值圖像對原始梯度圖像強制標記，屏蔽掉原梯度圖像所有局部極小值。最后，用分水嶺算法分割新標記的梯度圖像，完成路面裂縫圖像的分割過程。
算法總流程圖如圖2所示。
2.2梯度圖像
由于分水嶺算法與圖像的梯度聯系很大，所以，在梯度圖像上進行分水嶺分割，其結果更加準確。利用數學形態學梯度算子可以得到較好的梯度圖像。形態學梯度計算過程如下：
定義I為原始圖像，B為結構元素，那么數學形態學梯度定義有3種：
（1）腐蝕型：
Ie(x,y)=I－IΘB（1）
（2）膨脹型：
Id(x,y)=I⊕B－I（2）
（3）膨脹腐蝕型：
Ide(x,y)=I⊕B－IΘB（3）
腐蝕型和膨脹型梯度圖像都發生了檢測位置偏移，而式（3）同時檢測內邊緣和外邊緣，其定位準確，因此，本文選擇式（3）計算梯度圖像I。
2.3梯度圖像高通濾波
巴特沃斯濾波器（）頻帶內的頻率響應曲線是最平滑的，是比較理想、經典的濾波器。利用濾波器對上節得到的梯度圖像I進行高通濾波，得到梯度圖像的高頻成分，步驟如下：
（1）將梯度圖像I進行FFT變換到頻率域：
F(I)(u,v)=f［(－1)x+yI(x,y)］（4）
（2）在頻率域中對F(I)(u,v)進行二階高通濾波(BHPF)，得到高頻成分：
GBHPF(I)(u,v)=HBHPF(u,v)×F(I)(u,v)={1+［D0/［(u－M/2)2+(v－N/2)2］1/2］2}×F(I)(u,v)（5）
（3）再經過IFFT變換回到空間域，得到空間域的高頻成分：
IBHPF(x,y)=(－1)x+y×I－1×［R［GBHPF(I)(u,v)］］（6）
其中，R()代表取實部計算。
2.4相關特征約束提取標記
采用標記來約束分水嶺的過分割，是通過搜尋出與圖像中不同目標區域具有相關性的一系列局部極小值進行標記，每個極小值作為一個種子點，將其周圍具有相似性的像素點合并到該區域中，不斷擴展生長，得到像素具有相關性的聚水盆。而不符合每個局部極小值相關性的像素將不被合并，最終形成了聚水盆之間的分水嶺。
從獲取的高頻成分IBHPF(x,y)中提取與路面裂縫信息高度相關的局部極小值作為標記，縮小了局部極小值的范圍，防止了偽極小區域的形成，構成了標記圖像Imark。
本文利用形態學擴展最小變換技術（Hminima）提取標記，根據梯度圖像中目標物體對應的極小值一般小于噪聲對應的極小值的思想，預先設定一個集水盆深度閾值h，深度大于h的極小值即為與目標物體相關的、要提取的標記；深度小于h的極小值為噪聲對應的，不被標記，是需要被填平的淺水盆，這樣極大地消弱了過分割。參數h的確定是提取的關鍵，過大的h易造成真實的極小值丟失，出現欠分割的情況；h過小又不能達到消除過分割的目的。
本文在對路面裂縫圖像進行分析的基礎上，根據H-minima思想，對大量圖像局部極小值統計分析，得到獲取h的自適應方法：
其中，M0為梯度圖像的均值，M1為梯度圖像極小值的均值，M2為局部極大值的均值，M2－M1為聚水盆平均深度。圖3為示意圖。
圖3h相關參數示意圖對得到的高頻成分IBHPF(x,y)進行Hminima變換，得到標記圖像Imark：
Imark=Hmin(IBHPF,h)（8）
其中，Hmin(*)是 Hminima變換。
3改進的分水嶺算法分割路面裂縫圖像
通過Hminima變換技術大大消弱了過分割，因此，得到與裂縫信息相關的極小值標記圖像Imark后，采用Solid提出的形態學強制極小值標定操作lm min(*)強制將提取的標記作為原梯度圖像的極小值，從而屏蔽掉所有原梯度圖像的極小值，達到解決過分割的目的。
Imark對原梯度圖像I采用強制極小值標定運算，得到修改極小值后的梯度圖像Iws：
Iws=lm min(I,Imark)(9)
最后，利用形態學分水嶺變換分割標定后的梯度圖像Iws，得到最終分割結果圖像Iws:
Iws=(I)(10)
其中，(*)表示分水嶺變換操作。
4仿真研究與結果分析
圖4對裂縫信息較弱的網狀裂縫圖像分割結果比較為驗證本文算法，在MATLAB7.0環境下進行實驗，實驗采用連霍高速公路三門峽段拍攝的50幅路面裂縫圖像。選擇了三種比較有代表性的路面裂縫圖像，分別采用各種經典算法和本文算法進行分割，對最終分割出的效果進行比較，如圖4~圖6所示。其中，原圖1為裂縫信息比較弱的網狀裂縫圖像，原圖2為白噪聲較多的圖像，原圖3為顆粒紋理比較大的瀝青路面。
傳統的閾值分割法由于閾值難以確定，容易出現兩種情況：閾值過小產生的噪聲比較多，閾值過大微弱裂縫信息被去除掉。使用經典的分水嶺法分割出的圖像過分割現象比較嚴重，出現過多雜亂的小分支，甚至出現一些偽裂縫。使用本文算法，不僅縮小了分割過程中標定的極小區域，防止了過分割，還保留了微小的裂縫信息，非常利于路面裂縫圖像中裂縫信息的提取，從實驗結果圖中可以看出，其分割結果準確、清晰，消除了偽裂縫。
為進一步驗證本文算法的優勢，采用信噪比評價各種算法分割的結果，PSNR值越大，說明提取的裂縫信息越準確，反之，效果越不好。
其中，f(x,y)為原圖灰度像素值，I(x,y)為分割后圖像像素值，M和N分別為圖像的長和寬。對以上三幅圖像進行信噪比計算，結果如表1。
從表1數據可以看出，本文算法信噪比明顯高于其他兩種算法。
最后，將50幅圖像用以上3種方法進行綜合性能比較，結果如表2。
實驗證明，改進后的分割算法結果準確、噪聲少、效率高，在保證細微裂縫能被檢測出的同時極大縮減了過分割噪聲信息。
5結論
高速公路路面裂縫圖像復雜多樣，而且裂縫信息弱小，使用傳統方法很難檢測出準確清晰的裂縫信息。本文在分析了路面裂縫圖像本身的特點后，選擇基于三維地形模型的形態學分水嶺思想，利用裂縫圖像背景區域為低頻、裂縫信息對應高頻的思想，選擇與裂縫信息相關的局部極小值進行標記圖像，使得最終的分割結果避免了過分割現象。但本文算法在提取到的裂縫信息中有少量斷裂，對后續的計算測量有局部影響，還需要進一步的研究。
參考文獻
［1］ 唐磊,趙春霞,王鴻南,等.基于圖像三維地形模型的路面裂縫自動檢測［J］.計算機工程,2008,34(5):20-22.
［2］ 伯紹波.瀝青路面裂縫圖像檢測算法研究［D］.西安：長安大學，2008.
［3］ 孫惠杰，鄧廷權，李艷超.改進的分水嶺圖像分割算法［J］.哈爾濱工程大學學報，2014,35（7）：857-864.
［4］ 余旺盛,候志強,宋建軍.基于標記分水嶺和區域合并的彩色圖像分割［J］.電子學報，2011,39（5）：1007-1011.
［5］ 黎蔚，高璐.基于改進的分水嶺算法的路面裂縫檢測［J］.計算機工程與應用，2013，49（20）：263-266.
［6］ 伊力哈木.基于改進的自適應分水嶺圖像分割方法研究［J］.計算機仿真，2013,30（2）：373-377.
［7］ 夏倩倩，宋余慶，劉哲，等.結合分水嶺和互信息的醫學圖像自適應分割方法［J］ .計算機應用與軟件，2015,32（1）：221-224.
［8］ 高麗,楊樹元,李海強.一種基于標記的分水嶺圖像分割新算法［J］.中國圖象圖形學報，2007,12（6）：1025-1032.