劉 莉,靳 鴻,唐 波,陸 真,李志敏
(1.中北大學 電子測試技術國家重點實驗室,山西 太原 030051;2.中北大學 儀器科學與動態(tài)測試教育部重點實驗室,山西 太原 030051;3.北京航天自動控制研究所,北京 100854;4.宇航智能控制技術國防科技重點實驗室,北京 100854)
圖像增強是圖像處理的一種重要技術,已被廣泛應用在生物醫(yī)學工程、指紋識別、天文圖像處理等不同領域[1]。圖像增強的主要思想是通過增大圖像中明亮區(qū)域和暗淡區(qū)域的對比度,以獲得更好的視覺效果。紅外圖像受紅外成像傳感器、背景輻射和成像環(huán)境操作的影響,具有較低的對比度[2]。為了提取紅外小目標,需要對雜波背景里的暗淡目標進行增強。
數(shù)學形態(tài)學的基本運算是腐蝕和膨脹,以及由兩種運算組合的開和閉[11]。先腐蝕后膨脹的過程稱為開運算,先膨脹后腐蝕的過程稱為閉運算。灰度數(shù)學形態(tài)學是本文研究的基礎,其基本運算如式(1)-式(4)所示。
設f 是輸入圖像,B 為結構元素,Df和DB分別是f 和B 的定義域,則:
1)灰度膨脹
2)灰度腐蝕
3)灰度開運算
4)灰度閉運算
式中,a=x-u,b=y-v,c=x+u,d=y+v。
晚上9時許,專家義診車剛抵達一師醫(yī)院,轉運患兒的救護車也隨后趕到。由于提前獲得了阿拉爾醫(yī)院的術前相關檢查結果,經過簡短的術前談話和準備,患兒被直接送進了手術室,專家們立即投身到緊張的急診手術中……手術由錢云忠主任和方家杰副主任主刀,援疆專家張德林主任和曹振剛副主任實施麻醉管理,陳正副院長進行全程監(jiān)護,專家們一起為患兒手術保駕護航……最終,經過一個多小時的緊張手術,雖然患兒一側睪丸已經完全壞死,但另一側睪丸得以保全,術后的感染炎癥指標也均處于正常水平。
對灰度圖像的膨脹(或腐蝕)操作有兩類效果:如果結構元素的值都為正的,則輸出圖像會比輸入圖像亮(暗);根據(jù)輸入圖像細節(jié)中亮(暗)灰度值與它們的形狀相對于結構元素的關系,它們在運算中或被消減或被除掉;開運算可以平滑對應于所用結構元素的尺寸的圖像的明亮區(qū)域;閉運算可以平滑對應于所用結構元素的尺寸的圖像的暗淡區(qū)域。運用開運算和閉運算的算子,Meyer提出了高帽變換算子,高帽變換可以分為白高帽變換和黑高帽變換,分別表示為WTH 和BTH,定義如式(5)和式(6)所示:
5)白高帽變換
6)黑高帽變換
WTH 用來提取明亮圖像區(qū)域,BTH 用來提取暗淡圖像區(qū)域[11]。
WTH 和BTH 的定義表明,在傳統(tǒng)的高帽變換里只使用了一個結構元素,其僅僅能增強對應于所用結構元素尺寸的圖像特征,目標區(qū)域和其周圍臨近區(qū)域間的不同信息不能有效地得到處理。因此,高帽變換的性能有待提升和改進。
在圖像增強里,多尺度數(shù)學形態(tài)學發(fā)揮著重要作用。具有正方形狀的結構元素是一種廣泛使用的結構元素的形狀。nL 是正方形結構元素Bb的大小尺寸,nW 是正方形邊緣結構元素ΔB 的大小尺寸,nM 是在ΔB 里的邊緣的大小尺寸,如圖1所示。利用ΔB 和Bb,改進的高帽變換可以提取大小尺寸不大于nW -2×nM 的明亮或暗淡圖像區(qū)域。
圖1 使用的正方形結構元素Fig.1 A square structure element in use
設n個尺度的結構元素被使用,結構元素在每個尺度i(1≤i≤n)上的大小尺寸逐漸改變,如式(7)和式(8)所示。
式中,nS 是對于每個尺度的大小尺寸增長步長。
運用結構元Bi(u,v)(0≤i≤n-1)的多尺度的膨脹、腐蝕、開運算、閉運算、高帽變換可以表示如式(9)-式(14)所示:
式中,a=x-u,b=y-v,c=x+u,d=y+v。
改進的高帽變換利用兩個具有相同形狀的平坦結構元素Bin和Bout,Bin尺寸小于等于目標區(qū)域,Bout尺寸大于目標區(qū)域。令ΔB=Bout-Bin為邊緣結構元素,Bb為目標區(qū)域中的結構元素,Bin≤Bb≤Bout。基于ΔB 和Bb,提出了兩種新算子,分別為fBoi和f□Boi,其定義如式(15)和(16)所示:
式中,Bbi的大小為nLi;ΔBi的大小是nWi;從1到n改變i,多尺度明亮和暗淡圖像區(qū)域可以被提取出來。
在紅外圖像中,目標區(qū)域通常是明亮或暗淡的圖像。NWTH 可以用來提取明亮圖像區(qū)域,NBTH可以用來提取暗淡圖像區(qū)域。因為含有ΔB 和Bb,結構元素里包含了目標區(qū)域的不同信息,這提升了高帽變換對區(qū)域提取的性能。此外,負灰度值的抑制減少了一些噪聲,使目標增強更加有效。
低質量紅外圖像通常具有極少的灰度級并且目標區(qū)域和背景間的對比度低。為了有效地增強圖像的對比度,本文建立了一種基于加權的對比度增強算法,其定義如式(19)所示:
式(19)中,f 是原始圖像;fEn是最終增強圖像;w1,w2和w3是調整增強效果的權重。通過調整w1,w2和w3,最終圖像的對比度可以針對不同的應用而改變。
首先,提取多尺度明亮和暗淡圖像區(qū)域;然后,通過對提取的多尺度明亮和暗淡圖像區(qū)域上的最大值進行操作,建立了明亮和暗淡的圖像區(qū)域;最后,通過加權策略,使目標圖像的對比度得到了有效地增強。改進算法的執(zhí)行框圖如圖2所示。
圖2 改進算法的執(zhí)行框圖Fig.2 The implementation block diagram of improved algorithm
為了顯示提出算法的效力,選用了天空海洋等一系列紅外圖像進行增強,并用直方圖均衡化(HE),多尺度形態(tài)學方法(MSM)與提出的算法作比較。本實驗選用的尺度數(shù)為n=9,ΔB 的邊緣大小nM =2,nW =nL=5,增長步長nS=11,w1=1,w2=3并且w3=3,試驗在Matlab7.1上運行[12],其實驗結果如圖3和圖4所示。
圖3是海上航行紅外圖像的一個示例,目標區(qū)域是模糊的并且背景比較暗淡。圖像增強結果表明,HE可以增強目標區(qū)域,但是同時生成了很多噪聲,特別是在天空上;相對HE,MSM 對目標區(qū)域增強效果較好,但是很多其他圖像區(qū)域也被損壞,影響了視覺效果。而提出的算法很好地增強了目標區(qū)域并且生成了極少的噪聲,得到了一個良好的視覺效果,有效地提取了目標。
圖4是被大氣中的霧模糊的紅外圖像的一個示例。在圖像中有三個暗淡船的目標并且圖像具有天和海的背景。從圖4中可以看出,HE 不能很好地增強圖像,并且三個船目標區(qū)域均被平滑。MSM增強了圖像目標區(qū)域,但是很多其他非目標區(qū)域也被過度增強并且產生了很多噪聲。與其他算法相比,提出的算法性能相對最好,不僅增強了原始圖像的目標區(qū)域,也抑制了噪聲。
圖3 海上航行的紅外圖像Fig.3 The infrared images of the sea
圖4 被大氣中的霧模糊的紅外圖像Fig.4 The infrared images blurred by atmospheric fog
本文選用線性模糊性指數(shù)量度[8-10]對實驗中的算法作一個定量比較。線性模糊性指數(shù)基于空間域分析,廣泛用來定量比較不同算法的增強性能,其表示為γ,定義如式(20)、式(21)所示。
式中,pxy為圖像的像素;f(x,y)為像素pxy的灰度值;fmax是M×N 的圖像f 的最大灰度值。γ值越小,表明圖像增強算法的性能越好。
通過計算HE、MSM 和提出算法的γ 值,將每種算法的γ平均值作比較,如表1所示。從表1中可以看出,提出算法的γ值最小,這表明提出算法對于紅外圖像增強的性能最好。
表1 利用γ值的圖像增強定量比較Tab.1 The quantitative comparison of image enhancement byγ
本文提出了多尺度高帽變換的紅外圖像增強的優(yōu)化算法,該方法利用多尺度理論,對傳統(tǒng)的高帽變換進行了改進,尋找有效的結構元素,提取不同紅外圖像上的明亮或暗淡圖像信息。實驗結果表明,同其他算法相比,提出的算法對紅外圖像的對比度增強效果相對最好。通過計算幾種算法的線性模糊性指數(shù)量度,對幾種算法的增強效果進行定量比較得知,提出算法對于紅外圖像增強的性能最好。因此,可將其應用在各種低對比度紅外圖像的目標增強中,應用前景廣泛。
[1]孫慧賢,羅飛路,張玉華.基于小波變換和同態(tài)濾波的內窺圖像增強算法[J].探測與控制學報,2008,30(5):68-72.
[2]田鵬輝,隋立春,燕莎.紅外運動小目標檢測方法綜述[J].探測與控制學報,2013,35(2):76-80.
[3]徐衛(wèi)昌,黃威,李永峰,等.低通濾波與灰度值調整在圖像增強中的應用[J].激光與紅外,2012,42(4):458-462.
[4]尹士暢,喻松林.基于小波變換和直方圖均衡的紅外圖像增強[J].激光與紅外,2013,43(2):225-228.
[5]Gonzalez R C,Wintz P.Digital Image Processing[M].2nd Edition.US:Addison-Wesley,Reading,MA,2006.
[6]阮秋琦.數(shù)字圖像處理學[M].北京:電子工業(yè)出版,2001.
[7]崔屹.圖像處理及分析—數(shù)學形態(tài)學方法及應用[M].北京:科學出版社,2009.
[8]Serra J.Image Analysis Using Mathematical Morphology[M].London:Academic Press,1982.
[9]Serra J.Image Mathematical Morphology and its Applications to Image Processing[M].UK:Kluwer Academic Publishers,1994.
[10]袁俊.數(shù)學形態(tài)學理論及其在圖像處理中的應用[D].武漢:武漢理工大學,2007.
[11]朱士虎.形態(tài)學高帽變換與低帽變換功能擴展及應用[J].計算機工程與應用,2011,47(34):190-192.
[12]Stephen J Chapman.MATLAB編程[M].北京:科學出版社出版,2007.