約束最小二乘濾波的高斯模糊圖像復原研究

張艷紅 覃鳳清 姜 麗 文興東 何文杰 萬浩飛

約束最小二乘濾波的高斯模糊圖像復原研究

張艷紅 覃鳳清 姜 麗 文興東 何文杰 萬浩飛

（宜賓學院，四川 宜賓 644000）

文章分析了兩種常見的圖像復原濾波的復原原理及其局限后，在已知圖像退化函數的情況下，主要對高斯模糊圖像分別進行維納濾波和約束最小二乘法濾波復原處理并比較這兩種處理效果圖像的優劣之處，在上述基礎上，重點是分析不同高斯點擴展函數的尺寸及標準差的參數對高斯模糊圖像復原效果的影響。最后再綜合視覺效果和PSNR值得出高斯模糊圖像使用約束最小二乘濾波復原的效果較佳，并且在真實高斯點擴展函數參數附近的參數復原效果會更好。

高斯模糊；圖像復原；約束最小二乘濾波；維納濾波；點擴展函數

引言

圖像復原算法有很多種，主要包含非盲復原算法和盲復原算法兩類。最早的是Harris提出使用點擴散函數模型來處理由于大氣干擾引起的模糊圖片，隨后相繼提出了逆濾波、維納濾波、約束最小二乘濾波等非盲復原算法[1]。圖像復原首先要有一個明確的圖像復原質量準則，處理時需要對圖像退化過程先進行估算，并且補償退化過程造成的失真，以便獲得未經干擾退化的原始圖像或原始圖像的最優估值[2]，衡量與原始場景圖像的接近程度。上述問題的關鍵在于復原模型的建立，復原模型必須根據引起圖像退化的數學模型對退化圖像進行處理，然后通過學習建立它們之間的映射關系，從而對高分辨率圖像進行引導和估計，復原丟失的信息，進而達到恢復圖像的目的。實際運用可以使用不同濾波的手段對同一個圖像進行圖像復原處理，使復原后圖像符合一定的視覺標準，進而達到復原圖像的目的。

本文主要解決高斯模糊圖像復原的問題，高斯模糊是生活中常見的模糊圖像，解決此類圖像的模糊問題可以使后續應用更加方便和高速。文中首先介紹了高斯模糊圖像的退化模型，然后研究了維納濾波和約束最小二乘濾波的復原原理和局限性，然后通過上述兩種不同復原算法對高斯模糊圖像進行處理，得出處理高斯模糊圖像較好的復原算法。最后，分析了不同高斯點擴散函數的大小和標準差參數對高斯模糊圖像恢復的影響，并得出了相應的結論。

1 高斯模糊

1.1 高斯圖像退化模型

圖像模糊是由于像素值的局部平均化而導致圖像內容的平滑[3]，在圖像的成像過程中，目標會高速運動、散射、成像系統畸變和噪聲等干擾，使最終的圖像出現很多問題，稱為“退化”。圖像復原的關鍵是建立退化模型[4]，通常，使用一致的數學模型來描述圖像退化過程，然后建立退化模型。對模型進行精確處理，然后反向推導，確定復原圖像復原的準確性。然后，將圖像降晰過程抽象化，粗略概括為將高清圖像與未知點擴展函數進行卷積，并加入噪聲，最終得到降晰圖像。

圖1 高斯圖像退化模型

2 圖像復原處理算法

2.1 維納濾波

常見的圖像復原處理包括維納濾波處理、逆濾波處理、約束最小二乘濾波處理等復原算法，在這些復原方法中，維納濾波的圖像復原處理是一種常見的濾波器，是去除退化圖像的模糊和噪聲的理想濾波器[6]，被廣泛應用于各種圖像檢測處理。在Wiener（維納）濾波圖像復原中，中和了退化函數和噪聲統計這兩個方面來進行圖像的復原處理。

維納濾波復原的主要原理是假設圖像噪聲是一個隨機過程，找到一個沒有失真圖像的估值，使得它們之間的均方誤差達到最小，從而實現圖像復原的效果。在圖像復原過程之中，實際使用到的是非因果Wiener（維納）濾波器，因為它的頻域形式比較簡單，計算效率比較高[7]。

維納濾波是一種線性圖像復原方法，尋找一個使統計誤差函數：

在頻域可以表示為：

2.2 約束最小二乘濾波

其約束條件為：

從公式中可以看到，約束最小二乘濾波在高頻區域中信噪比較低，但是由于通常噪聲是出現在高頻部分，所以濾波器在抑制噪聲的同時也會過濾掉一些圖像本身就有用的高頻信息，這說明約束最小二乘濾波在濾波復原過程中，可以減少對噪聲的放大作用[10]，從而達到圖像復原的效果。對于不同原始圖像的先驗知識，采取克服模糊的不同方案，約束最小二乘解具有不同的具體形式。約束最小二乘算法對噪聲較為敏感，在無噪聲或可以忽略的情況下可以獲得比較理想的復原效果[11]。

3 高斯點擴展函數

圖像模糊的成因有成像系統的像素差、有限寬帶造成的圖像失真、模擬圖像數字化過程中，由于損失掉部分細節，因而造成圖像質量下降等問題。在現實生活中，很多光學測量和成像系統的降晰函數均趨近于高斯點擴展函數（PSF），如顯微鏡成像系統、光學相機都屬于高斯型。

通常，高斯點擴展函數表示為如下形式：

其中Sup表示支持域，通常被表示為一個K*K的矩形區域、K被稱為高斯點擴展函數的尺寸，通常是3～15的奇數、σ為高斯點擴展函數的標準差。因此，高斯模糊圖像主要改變的是高斯點擴展函數的尺寸和標準差。

4 實驗及對比

為了更好地驗證約束最小二乘濾波對高斯模糊圖像的復原效果，實驗是通過對Lena原始圖像如圖2（尺寸為256×256）加以高斯點擴展函數的尺寸為9、標準差為3時生成的高斯模糊圖像如圖3所示，高斯模糊圖像的PSNR值為31.6277 dB，通過對高斯模糊圖像分別使用維納濾波和約束最小二乘濾波復原處理，進而比較出約束最小二乘濾波與維納濾波的復原效果的差異，接著再研究不同高斯點擴展函數的參數對復原效果的影響。

圖2 原始圖像

圖3 高斯模糊

4.1 不同復原算法的復原效果比較

為了評價不同復原算法的復原效果，在原始高斯點擴展函數參數情況下進行復原。此時高斯點擴展函數估計尺寸為9、標準為3，通過維納濾波及約束最小二乘濾波算法進行圖像復原，所得到復原圖像及圖像PSNR值如圖4、圖5所示。

（PSNR：27.1937 dB）

（PSNR：33.4747 dB）

通過這組對比實驗可以看出對于高斯模糊圖像，進行維納濾波和約束最小二乘法濾波復原處理后，從復原圖像質量以及客觀評價標準來看，約束最小二乘濾波復原的效果比維納濾波復原效果好，細節保持能力較優，因此對于高斯模糊圖像，應更多使用約束最小二乘濾波復原方法。

4.2 高斯點擴展函數參數對復原圖像質量的影響

通過改變高斯點擴展函數的參數（即高斯點擴展函數的尺寸及標準差）對同一模糊圖像進行約束最小二乘法濾波復原處理，并對這幾組數據進行比較，分析高斯點擴展函數參數對復原圖像質量的影響。

4.2.1 高斯點擴展函數尺寸對復原圖像質量的影響

為了驗證高斯點擴展函數尺寸參數對復原圖像質量的影響，當高斯點擴展函數的標準差為3時，在尺寸分別為5、7、11、13時通過約束最小二乘濾波進行圖像復原，得到的圖像如圖6所示，復原圖像的PSNR值見表1。

圖6 不同尺寸的復原圖像

表1 不同尺寸復原圖像的PSNR值

高斯點擴展函數尺寸PSNR（dB） 530.8912 729.8729 933.4747 1128.8632 1328.7225

圖7 不同尺寸復原圖像的PSNR值

通過圖7及PSNR值可以看出在高斯點擴展函數的標準差一定的情況下，改變高斯點擴展函數的尺寸對圖像復原質量的影響是明顯的。當高斯點擴展函數的尺寸逐漸接近原始尺寸時復原效果會越來越好，PSNR值逐漸在增大，一旦到達原始尺寸后，接著增加數值圖像復原效果會逐漸下降。

4.2.2 高斯點擴展函數標準差對復原圖像質量的影響

為了驗證高斯點擴展函數標準差參數對復原圖像質量的影響，當高斯點擴展函數的尺寸為9時，在標準差分別取2、2.5、3.5、4時通過約束最小二乘濾波進行圖像復原，得到的圖像如圖8所示，復原圖像的PSNR值見表2。

圖8 不同標準差的復原圖像

表2 不同標準差復原圖像的PSNR值

高斯點擴展函數標準差PSNR（dB） 230.3736 2.529.2014 333.4747 3.529.8002 428.7415

圖9 不同標準差復原圖像的PSNR值

通過圖9及PSNR值可以看出在高斯點擴展函數的尺寸一定的情況下，改變高斯點擴展函數的標準差得到不同的圖像。當高斯點擴展函數的標準差逐漸接近原始標準差時效果會越來越好，當達到原始標準差時，再增加標準差復原效果不會很明顯。因此，要想得到更精確的復原結果，高斯點擴展函數的參數取值還需進一步深入研究。

5 結束語

圖像復原技術與其他基本圖像處理技術相類似，都是一種以獲取視覺效果在某個方面提升為目的，是圖像處理領域中一種非常重要的技術。文中通過對上述兩種不同算法的研究分析及對計算機模擬出的高斯模糊退化圖像進行復原處理，從復原圖像質量以及客觀評價標準來看，得到約束最小二乘濾波復原處理能更好的復原高斯模糊圖像，在此基礎上改變高斯點擴展函數參數時越靠近真實值圖像復原效果越好這一結論。此外對于不同的模糊圖像對應的最佳復原方法均不同，后續需要在復原方法的參數估值方面研究，更重要的是對算法的深入研究。

[1] 李亞成. 高斯模糊圖像復原與匹配算法研究[D]. 武漢: 華中科技大學，2019.

[2] 黃力宇，黃登鳳，張園. “醫學影像處理”課程實驗軟件的開發[J]. 電氣電子教學學報，2012，34(6): 82-84，109.

[3] 張墨華. 基于先驗學習的圖像復原技術研究[D]. 鄭州: 戰略支援部隊信息工程大學，2019.

[4] 王鵬飛. 運動模糊圖像復原算法研究[D]. 重慶: 重慶大學，2013.

[5] 謝盛華，張啟衡，宿丁. 基于APEX方法的改進圖像復原算法[J]. 激光與紅外，2007(2): 185-188.

[6] 陳樹越，朱雙雙，蔣星，等. 基于高斯型點擴展函數的紅外圖像熱源復原[J]. 激光技術，2016，40(2): 270-273.

[7] 李慎之. 毫米波無源成像超分辨算法研究及DSP實現[D]. 成都: 電子科技大學，2011.

[8] 劉書君. 運動模糊圖像復原技術研究[J]. 信息通信，2017(4): 4-5.

[9] 萬宇. 空間變化模糊的圖像復原算法[D]. 杭州: 浙江工業大學，2018.

[10] 李國立，段汕. 基于約束最小二乘數字圖像恢復[J]. 電腦知識與技術，2008(14): 914，925.

[11] 石明珠. 運動模糊圖像全變分復原理論及關鍵技術研究[D]. 北京: 北京理工大學，2014.

Research on Gaussian Blur Image Restoration Based on Constrained Least Squares Filtering

After analyzing the restoration principles and limitations of two common image restoration filters, this paper mainly performs Wiener filtering and constrained least square filtering restoration on Gaussian blur image respectively when the image degradation function is known, and compares the advantages and disadvantages of these two processing effects. On the basis of the above, the emphasis is to analyze the influence of the size of different Gaussian point spread functions and the parameters of standard deviation on the restoration effect of Gaussian blur image. Finally, combining the visual effect and PSNR value, it is found that the effect of Gaussian blur image restoration using constrained least square filter is better, and the parameter restoration effect near the parameters of the real Gaussian point spread function is better.

Gaussian blur; image restoration; constrained least squares filtering; Wiener filter; point spread function

TN911

A

1008-1151(2022)09-0010-04

2022-07-10

四川省教育廳項目（18ZA0538）；宜賓市科技局項目（2017ZSF009-9）；宜賓學院科研項目（2020YY02）；國家級大學生創新項目（202010641021、201910641016）；四川省大學生創新項目（S201910641096、S202010641126）；宜賓學院大學生創新項目（X106412020336）。

張艷紅（1999－），女，安徽阜陽人，宜賓學院學生，研究方向為數字圖像處理。

覃鳳清（1976－），女，四川南充人，宜賓學院研究員，博士，研究方向為數字圖像處理。