劉 彬，薄 華

（上海海事大學(xué) 信息工程學(xué)院，上海 201306）

艦船在航海中拍攝的視頻圖像可以廣泛應(yīng)用于船員的人臉識(shí)別、表情識(shí)別、疲勞檢測(cè)、艦船安檢工作等。圖像質(zhì)量的好壞關(guān)系到整個(gè)圖像處理系統(tǒng)的性能，一幅清晰的圖像是這些應(yīng)用的前提。但在傳感器成像過程中，記錄介質(zhì)積分時(shí)間內(nèi)拍攝目標(biāo)和攝像機(jī)之間的相對(duì)運(yùn)動(dòng)會(huì)造成圖像的模糊，給后續(xù)的圖像處理和分析帶來一定的困難。

運(yùn)動(dòng)模糊的恢復(fù)算法研究大部分是針對(duì)全局模糊圖像的，而有一些時(shí)候圖像會(huì)出現(xiàn)局部模糊的情況，這時(shí)需要針對(duì)局部模糊區(qū)域進(jìn)行恢復(fù)，并將復(fù)原的區(qū)域無縫合成到整幅圖中，因此，對(duì)由運(yùn)動(dòng)產(chǎn)生的局部模糊區(qū)域恢復(fù)和合成是一項(xiàng)很有意義的研究。目前針對(duì)局部模糊圖像的恢復(fù)有：非負(fù)支持域遞歸逆濾波算法[1]、維納濾波法、L-R 算法[2]、基于Z 變換的恢復(fù)算法[3]等，圖像合成有：基于小波變換的圖像融合、基于a 分量的圖像合成、基于梯度場(chǎng)的圖像合成、基于多分辨率模型的圖像合成[4]。為使合成圖像在紋理和噪聲上具有一致性，Sunkavalli[5]提出了平滑直方圖匹配的過采樣的拉普拉斯金字塔構(gòu)建方法，該金字塔中進(jìn)行子帶的平滑直方圖匹配能有效的減小邊界及周圍存在明顯的光暈和合成圖像的人造痕跡。

在文獻(xiàn)[2，5]的基礎(chǔ)上提出一種新的局部運(yùn)動(dòng)模糊圖像恢復(fù)與合成算法，把模糊區(qū)域從圖像中分割出來，把復(fù)雜的局部模糊恢復(fù)問題轉(zhuǎn)化為前景模糊恢復(fù)、前景和背景融合兩個(gè)子問題來解決，對(duì)局部模糊區(qū)域進(jìn)行恢復(fù)，最后用基于直方圖匹配的小波變換對(duì)恢復(fù)圖像與背景圖像無縫的合成。實(shí)驗(yàn)結(jié)果表明，采用基于直方圖匹配的小波變換能較好地將恢復(fù)圖像合成到原圖中去，同時(shí)還能改善運(yùn)動(dòng)模糊圖像恢復(fù)效果。

1 局部運(yùn)動(dòng)模糊圖像退化模型

運(yùn)動(dòng)模糊的模型中，水平方向的勻速直線運(yùn)動(dòng)模糊具有代表性和簡(jiǎn)單性，其它方向上的運(yùn)動(dòng)模糊可以由水平方向上的運(yùn)動(dòng)模糊方法推廣得到。一幅理想的清晰圖像為f（x，y），其中前景圖為q（x，y），背景圖像為b（x，y）[3]，則：

f（x，y）經(jīng)過退化過程h（x，y），假定在退化過程中背景保持不變，只有前景區(qū)域發(fā)生運(yùn)動(dòng)模糊，疊加一個(gè)噪聲后，得到一幅局部退化的圖像，這一過程可描述為：

其中，* 表示卷積，g（x，y）表示退化圖像，h（x，y）為退化函數(shù)，n（x，y）又稱為點(diǎn)擴(kuò)散函數(shù)，為加性噪聲，為了簡(jiǎn)化運(yùn)動(dòng)模糊恢復(fù)問題，本文暫不考慮噪聲的影響。

2 模糊區(qū)域恢復(fù)

運(yùn)動(dòng)模糊圖像復(fù)原一般分為兩步：首先通過系統(tǒng)辨識(shí)求解h（x，y），然后采取相應(yīng)算法由模糊圖像g（x，y）和點(diǎn)擴(kuò)散函數(shù)h（x，y）來恢復(fù)出原圖像f（x，y）?？紤]到分割模糊區(qū)域會(huì)引入背景干擾，文中提出L-R 算法與亮度調(diào)整相結(jié)合的恢復(fù)算法。Lucy-Richardson（L-R）算法是一種非線性迭代方法，是從最大似然公式引出來的，在這種方程中，圖像是用泊松統(tǒng)計(jì)加以模型化的。當(dāng)這個(gè)迭代收斂時(shí)，模型的最大似然函數(shù)可以得到一個(gè)方程：

對(duì)L-R 算法恢復(fù)后的圖像采用亮度調(diào)整的方法，可以調(diào)整圖像的亮度，得到一幅清晰的圖像。為了調(diào)整每個(gè)像素的亮度，在每個(gè)點(diǎn)上設(shè)置一個(gè)虛擬的曝光時(shí)間作為參數(shù)，通過延長(zhǎng)或縮短曝光時(shí)間來增加或降低圖像的亮度，做法是在每個(gè)點(diǎn)乘上一個(gè)亮度縮放因子。

這里，I 表示輸入的圖像亮度，I′是調(diào)整后的亮度；GMap表示每個(gè)點(diǎn)的曝光調(diào)整因子，當(dāng)GMap（x，y）＞1.0 時(shí)提升（x，y）的亮度，GMap（x，y）＜1.0 時(shí)降低（x，y）的亮度。

2.1 模糊區(qū)域的提取

物體相對(duì)于像機(jī)是運(yùn)動(dòng)的，用背景差值法可以快速的從背景中提取模糊區(qū)域。當(dāng)前相機(jī)獲取的圖像為C（x，y），背景圖像B（x，y），兩者相減得到差值圖像D（x，y），對(duì)差值圖像D（x，y）進(jìn)行閾值化處理，如下式：

其中（x，y）是B（x，y）和C（x，y）重疊區(qū)域中的點(diǎn)，Th是給定的閾值。差值圖像不為0 的像素點(diǎn)是物體的位置區(qū)域，對(duì)閾值化后的差值圖像進(jìn)行開閉運(yùn)算之后再進(jìn)行水平和垂直投影到水平和垂直軸上[1]，通過對(duì)兩軸上的灰度進(jìn)行零值判斷即可得出模糊物體所在的區(qū)域，利用所得到的矩形區(qū)域，就可以把模糊物體從圖像中提取出來。

2.2 基于Radon 變換的模糊角度估計(jì)

對(duì)運(yùn)動(dòng)模糊圖像來說，點(diǎn)擴(kuò)展函數(shù)可以描述為：

其中d為模糊長(zhǎng)度，θ為模糊角度。

Radon 變換是計(jì)算圖像在某一特定角度射線上投影（線積分）的方法：

R0（x′）表示積分結(jié)果，其中

二值圖像而言，某方向上積分（這時(shí)即點(diǎn)個(gè)數(shù)）很大表明該方向上有較強(qiáng)的線性，也就是存在線段；對(duì)模糊圖像的二維頻譜函數(shù)沿條紋方向計(jì)算線積分，所得結(jié)果為投影變換，其中心位置附近的積分值最大[7]。用Radon 變換來獲取參數(shù)，在運(yùn)動(dòng)模糊圖像的傅立葉頻譜上能看到明顯的平行暗線，暗線與軸正方向夾角是模糊角度逆時(shí)針旋轉(zhuǎn)90°。算法流程如下[8]：

1）對(duì)提取出的模糊區(qū)域進(jìn)行傅里葉變換，將圖像變換到頻域，對(duì)頻域圖像進(jìn)行LOG 變換；

2）將頻譜圖二值化，對(duì)二值圖像進(jìn)行邊緣檢測(cè)，得到頻譜圖的邊緣檢測(cè)圖；

3）對(duì)邊緣圖像進(jìn)行Radon 變換，并畫出直線，在中心區(qū)域得到多條直線，根據(jù)直線公式可以求出模糊角度θ。

2.3 基于自相關(guān)函數(shù)的運(yùn)動(dòng)模糊尺度的估計(jì)

在Radon 變換中，頻譜上暗線的個(gè)數(shù)是模糊尺度，在模糊尺度較大的情況下，直接讀取暗線條數(shù)會(huì)不準(zhǔn)確，因此在得到模糊角度θ 后，采用自相關(guān)函數(shù)來估計(jì)模糊尺度。

模糊帶內(nèi)相近的像素點(diǎn)灰度值形成了低頻區(qū)域，所以模糊像素軌跡內(nèi)部的像素值更加相關(guān)；而沿運(yùn)動(dòng)方向求微分后圖像的自相關(guān)函數(shù)可以得到一個(gè)對(duì)稱圖，圖中有一個(gè)中心峰值和對(duì)稱分布在峰值兩邊的最大負(fù)峰值，負(fù)峰與中心峰的距離就是模糊尺度[7]。算法流程如下：

2）計(jì)算g′（i，j）在行方向上自相關(guān)s（i，j），s（i，j）的每行都包含一對(duì)共軛相關(guān)峰，對(duì)稱分布在中心峰兩側(cè)；

3）s（i，j）在列方向上求和，得到一行數(shù)據(jù)sadd（·），低通濾波后可以突出最大負(fù)峰，有效抑制噪聲和局部尖峰點(diǎn)的干擾；

3 復(fù)原圖像與背景圖像合成

對(duì)模糊區(qū)域進(jìn)行了恢復(fù)后，需要將其無縫地合成到背景圖像中去。本文在模糊區(qū)域分割過程中使用背景差值法，對(duì)物體的提取會(huì)得到一個(gè)大于物體面積的矩形區(qū)域，因此背景可能會(huì)產(chǎn)生干擾，直接合成將造成圖像間分層和明顯跳躍過渡現(xiàn)象。

為使合成圖像在紋理和噪聲上具有一致性，分析并利用圖像在不同尺度下的特性是一類有效的方法?；诙喾直媛屎铣傻挠行Х椒ㄊ窍葘?duì)拉普拉斯金字塔分解后的子帶圖像進(jìn)行直方圖匹配，然后由這些子帶重構(gòu)出合成圖像；但直接對(duì)子帶進(jìn)行直方圖匹配會(huì)導(dǎo)致邊界周圍有明顯的光暈，合成痕跡十分明顯[4]。本文提出一種基于直方圖匹配的非下采樣的小波變換法，靜態(tài)二維離散小波變換[9]，與離散小波變換不同的是靜態(tài)二維離散小波分解得到的子帶圖像與原圖像具有相同尺寸，它不對(duì)分解系數(shù)進(jìn)行下采樣，所以單層和多層分解的結(jié)果一樣；重構(gòu)時(shí)，對(duì)經(jīng)過重構(gòu)濾波后的信號(hào)不做上采樣。對(duì)分解后的子帶進(jìn)行直方圖匹配，匹配后的子帶可以通過求解與分解濾波器相同的線性方程重構(gòu)為合成圖像。

4 實(shí)驗(yàn)結(jié)果與分析

為了驗(yàn)證該方法的有效性，在MATLAB 平臺(tái)下，本文對(duì)一張含局部運(yùn)動(dòng)模糊區(qū)域的圖像進(jìn)行模糊區(qū)域的恢復(fù)與合成的實(shí)驗(yàn)。

圖1 原始圖像Fig.1 Original image

圖2 模糊區(qū)域Fig.2 The blurred region

圖3 維納濾波（左）與本文方法（右）恢復(fù)效果對(duì)比Fig.3 Contrast of the restoration effect between wiener filter（left）and the proposed method（right）

圖4 直接合成法（左）與直方圖匹配法（右）合成效果對(duì)比Fig.4 Contrast of the compositing effect between direct compositing（left）and histogram matching method（right）

圖5 直接合成細(xì)節(jié)部分（左）與直方圖匹配法細(xì)節(jié)部分（右）對(duì)比Fig.5 Contrast of the detail section between direct compositing（left）and histogram matching method（right）

從維納濾波和本文算法（迭代次數(shù)為50 次）恢復(fù)效果對(duì)比圖3 來看，本文提出的算法更適合分割的模糊區(qū)域中存在背景干擾的情況，得到的恢復(fù)效果比維納濾波恢復(fù)的效果要更好一些；圖4 表示直接將恢復(fù)區(qū)域合成到背景中與本文算法合成效果對(duì)比圖，本文提出直方圖匹配的算法，能更好地做到了無縫合成，但是在細(xì)節(jié)上還存在著人工合成的痕跡；圖5 是對(duì)圖4 中合成效果圖的細(xì)節(jié)部分放大，結(jié)果顯示，在直方圖匹配的過程中，原圖的紋理細(xì)節(jié)被部分加入到恢復(fù)圖像，從而達(dá)到了解決恢復(fù)圖像中的震蕩問題。

5 結(jié)論

本文提出一種針對(duì)航海過程中拍攝的圖像出現(xiàn)局部模糊現(xiàn)象進(jìn)行恢復(fù)與合成的算法，該方法在部分背景的干擾下，仍可以達(dá)到較好的圖像的去模糊效果，是一種有效可行的方法；直方圖匹配的小波變換算法能將恢復(fù)圖無縫地合成到背景圖像中；不足之處在于沒有討論復(fù)雜背景下運(yùn)動(dòng)模糊物體的提取和噪聲對(duì)參數(shù)求解的影響，算法的實(shí)時(shí)性還有待改善并且分層現(xiàn)象尚未完全消除，下一步將就這幾個(gè)問題對(duì)方法進(jìn)行優(yōu)化與改進(jìn)。

[1]李晉武，張娟，鄒麗暉.含運(yùn)動(dòng)物體的運(yùn)動(dòng)場(chǎng)景模糊圖像恢復(fù)方法[C]//Proceedings of the 30th Chinese Conference，Yantai：3287-3291.

[2]季冰川.局部運(yùn)動(dòng)模糊圖像恢復(fù)系統(tǒng)研究和實(shí)現(xiàn)[D].贛州：江西理工大學(xué)，2008.

[3]黎明和，何斌，岳繼光，等.基于Z變換的局部勻速運(yùn)動(dòng)模糊圖像恢復(fù)算法[J].光學(xué)學(xué)報(bào)，2009，29（5）：1193-1197.LI Ming-he，HE Bin，YUE Ji-guang，et al.Blurred image resration of local uniform motion based on Z transrorm[J].Acta Optica Sinica，2009，29（5）：1193-1197.

[4]吳昊，徐丹.數(shù)字圖像合成技術(shù)綜述[J].中國(guó)圖象圖形學(xué)報(bào)，2012，17（11）：1333-1346.WU Hao，XU Dan.Survey of digital image compositing[J].Journal of Image and Graphics，2012，17（11）：1333-1346.

[5]Sunkavalli K，Johnson M，Matusik W，et al.Multi-scale image harmonization [C]//Proceedings of ACM SIGGRAPH 2010.New York，USA：ACM，2010（125）：1-10.

[6]Gonzalez R C，Woods R E.數(shù)字圖像處理（M）.MATLAB版.阮秋琦等譯.北京：電子工業(yè)出版社，2005.

[7]周玉，周序生，戴圣偉，等.一種改進(jìn)的運(yùn)動(dòng)模糊圖像復(fù)原算法[J].計(jì)算機(jī)工程與科學(xué)，2009，31（10）：49-52.ZHOU Yu，ZHOU Xu-sheng，DAI Sheng-wei，et al.An improved motion-blurred image restoration algorithm [J].Computer Engineering &Science，2009，31（10）：49-52.

[8]LIU Yuan-min，TIAN Lian-fang.Approach to optimizing restoration of motion blur images with robust blind deconvolution based on new cepstrum and total variation[C]//2012 International Conference on Computer Science and Electronics Engineering，2012：554-558.

[9]張德豐.MATLAB小波分析[M].北京：機(jī)械工業(yè)出版社，2009.