1、精選優質文檔-傾情為你奉上圖像復原方法綜述1、摘要圖像是人類視覺的基礎，給人具體而直觀的作用。圖像的數字化包括取樣和量化兩個步驟。數字圖像處理就是將圖像信號轉換成數字格式，并利用計算機進行加工和處理的過程。圖像復原是圖像處理中的一個重要問題，對于改善圖像質量具有重要的意義。解決該問題的關鍵是對圖像的退化過程建立相應的數學模型，然后通過求解該逆問題獲得圖像的復原模型并對原始圖像進行合理估計。本文主要介紹了圖像退化的原因、圖像復原技術的分類和目前常用的幾種圖像復原方法，詳細的介紹了維納濾波、正則濾波、LR算法和盲區卷積，并通過實驗證明了該方法的可行性和有效性。關鍵詞：圖像退化、圖像復原、維納濾波、

2、正則濾波、LR算法、盲區卷積、2、圖像復原概述在圖像的獲取、傳輸以及保存過程中，由于各種因素，如大氣的湍流效應、攝像設備中光學系統的衍射、傳感器特性的非線性、光學系統的像差、成像設備與物體之間的相對運動、感光膠卷的非線性及膠片顆粒噪聲以及電視攝像掃描的非線性等所引起的幾何失真，都難免會造成圖像的畸變和失真。通常，稱由于這些因素引起的質量下降為圖像退化。圖像退化的典型表現是圖像出現模糊、失真，出現附加噪聲等。由于圖像的退化，在圖像接受端顯示的圖像已不再是傳輸的原始圖像，圖像效果明顯變差。為此，必須對退化的圖像進行處理，才能恢復出真實的原始圖像，這一過程就稱為圖像復原1。圖像復原技術是圖像處理領域

3、中一類非常重要的處理技術，與圖像增強等其他基本圖像處理技術類似，也是以獲取視覺質量某種程度的改善為目的，所不同的是圖像復原過程實際上是一個估計過程，需要根據某些特定的圖像退化模型，對退化圖像進行復原。簡言之，圖像復原的處理過程就是對退化圖像品質的提升，并通過圖像品質的提升來達到圖像在視覺上的改善。由于引起圖像退化的因素眾多，且性質各不相同，目前沒有統一的復原方法，眾多研究人員根據不同的應用物理環境，采用了不同的退化模型、處理技巧和估計準則，從而得到了不同的復原方法。圖像復原算法是整個技術的核心部分。目前，國內在這方面的研究才剛剛起步，而國外卻已經取得了較好的成果。早期的圖像復原是利用光學的方法

4、對失真的觀測圖像進行校正，而數字圖像復原技術最早則是從對天文觀測圖像的后期處理中逐步發展起來的。其中一個成功例子是NASA的噴氣推進實驗室在1964年用計算機處理有關月球的照片。照片是在空間飛行器上用電視攝像機拍攝的，圖像的復原包括消除干擾和噪聲，校正幾何失真和對比度損失以及反卷積。另一個典型的例子是對肯尼迪遇刺事件現場照片的處理。由于事發突然，照片是在相機移動過程中拍攝的，圖像復原的主要目的就是消除移動造成的失真2。早期的復原方法有：非鄰域濾波法，最近鄰域濾波法以及效果較好的維納濾波和最小二乘濾波等。隨著數字信號處理和圖像處理的發展，新的復原算法不斷出現，在應用中可以根據具體情況加以選擇。目

5、前國內外圖像復原技術的研究和應用主要集中于諸如空間探索、天文觀測、物質研究、遙感遙測、軍事科學、生物科學、醫學影象、交通監控、刑事偵察等領域。如生物方面，主要是用于生物活體細胞內部組織的三維再現和重構，通過復原熒光顯微鏡所采集的細胞內部逐層切片圖，來重現細胞內部構成；醫學方面，如對腫瘤周圍組織進行顯微觀察，以獲取腫瘤安全切緣與癌腫原發部位之間關系的定量數據；天文方面，如采用迭代盲反卷積進行氣動光學效應圖像復原研究等。3、圖像退化模型圖像復原問題的有效性關鍵之一取決于描述圖像退化過程模型的精確性。要建立圖像的退化模型，則首先必須了解、分析圖像退化的機理并用數學模型表現出來。在實際的圖像處理過程中

6、，圖像均需以數字離散函數表示，所以必須將退化模型離散化3。對于退化圖像： (1) 如果上式中，按相同間隔采樣，產生相應的陣列、，然后將這些陣列補零增廣得到大小為的周期延拓陣列，為了避免重疊誤差，這里，。由此，當k=0,1,L,M-1;l=0,1,L,N-1時，即可得到二維離散退化模型形式： (2)如果用矩陣表示上式，則可寫為： (3)其中，,為一個行堆疊形成的列向量，H為階的塊循環矩陣?，F實中造成圖像降質的種類很多，常見的圖像退化模型及點擴展函數有如下情景15：(1) 線性移動降質在拍照時，成像系統與目標之間有相對直線移動會造成圖像的降質。水平方向線性移動可以用以下降質函數來描述： (4) 式

7、中，d是降質函數的長度。在應用中如果線性移動降質函數不在水平方向，則可類似地定義移動降質函數。(2) 散焦降質當鏡頭散焦時，光學系統造成的圖像降質相應的點擴展函數是一個均勻分布的圓形光斑。此時，降質函數可表示為： (5)式中，R是散焦半徑。(3) 高斯(Gauss)降質Gauss降質函數是許多光學測量系統和成像系統最常見的降質函數。對于這些系統，決定系統點擴展函數的因素比較多。眾多因素綜合的結果總是使點擴展函數趨于Gauss型。典型的系統可以舉出光學相機和CCD攝像機、相機、CT相機、成像雷達、顯微光學系統等。Gauss降質函數可以表達為： (6)式中，K是歸一化常數，是一個正常數，C是的圓形

8、支持域。4、幾種較經典的復原方法介紹圖像復原算法有線性和非線性兩類。線性算法通過對圖像進行逆濾波來實現反卷積，這類方法方便快捷，無需循環或迭代，直接可以得到反卷積結果，然而，它有一些局限性，比如無法保證圖像的非負性。而非線性方法通過連續的迭代過程不斷提高復原質量，直到滿足預先設定的終止條件，結果往往令人滿意。但是迭代程序導致計算量很大，圖像復原時耗較長，有時甚至需要幾個小時。所以實際應用中還需要對兩種處理方法綜合考慮，進行選擇4。(1) 維納濾波法維納濾波法是由Wiener首先提出的，應用于一維信號處理，取得了很好的效果。之后，維納濾波法被用于二維信號處理，也取得了不錯的效果，尤其在圖像復原領

9、域，由于維納濾波計算量小，復原效果好，從而得到了廣泛的應用和發展。維納濾波器尋找一個使統計誤差函數 (7)最小的估計。E是期望值操作符，是未退化的圖像。該表達式在頻域可表示為 (8)其中， 表示退化函數 表示的復共軛表示噪聲的功率譜表示未退化圖像的功率譜比率稱為信噪功率比。在IPT中維納濾波使用函數deconvwnr來實現的。 模擬實驗結果如下：(2) 正則濾波法另一個容易實現線性復原的方法稱為約束的最小二乘方濾波，在IPT中稱為正則濾波，并且通過函數deconvreg來實現。在最小二乘復原處理中，常常需要附加某種約束條件。例如令Q為f的線性算子，那么最小二乘方復原的問題可以看成使形式為的函數

10、，服從約束條件的最小化問題，這種有附加條件的極值問題可以用拉格朗日乘數法來處理。尋找一個，使下述準則函數為最?。?(9)式中叫拉格朗日系數。通過指定不同的Q，可以得到不同的復原目標。模擬實驗結果如下：(3)Lucy-Richardson算法L-R算法是一種迭代非線性復原算法，它是從最大似然公式印出來的，圖像用泊松分布加以模型化的。當下面這個迭代收斂時模型的最大似然函數就可以得到一個令人滿意的方程： (10)*代表卷積，代表未退化圖像的估計，g和h和以前定義一樣。在IPT中，L-R算法由名為deconvlucy的函數完成的。模擬實驗結果如下：(4)盲去卷積在圖像復原過程中，最困難的問題之一是，如

11、何獲得PSF的恰當估計。那些不以PSF為基礎的圖像復原方法統稱為盲區卷積。它以MLE為基礎的，即一種用被隨機噪聲所干擾的量進行估計的最優化策略。工具箱通過函數deconvblind來執行盲區卷積。模擬實驗結果如下：圖像分別迭代5次，10次，20次的結果。5、參考文獻1 羅軍輝馮平等。MATLAB7.0在圖像處理中的應用M，北京：機械工業出版社2006，202204，257278。2 劉維一，于德月，王肇圻等。用迭代法消除數字圖像放大后的模糊J。光電子.激光，2002，13（4）：398400。3 陶洪。數字共焦顯微技術及其圖像復原算法研究四川大學碩士學位論文CNKI:CDMD:10610.2.

12、2003.6632。4 岡薩雷斯等著。數字圖像處理（MATLAB版）M，電子工業出版社2009。6、程序附錄清單 (1) 維納濾波clcclear allclose allI=imread('C:Documents and SettingsAdministrator桌面原始圖.jpg');ubplot(231);Imshow(I);I=rgb2gray(I);subplot(232);imshow(I);noise=0.1*randn(size(I);PSF=fspecial('motion',21,11);Blurred=imfilter(I,PSF,'

13、;circular');BlurredNoisy=im2uint8(Blurred);NP=abs(fftn(noise).2;NPOW=sum(NP(:)/prod(size(noise);NCORR=fftshift(real(ifftn(NP);IP=abs(fftn(I).2;IPOW=sum(IP(:)/prod(size(noise);ICORR=fftshift(real(ifftn(IP);ICORR1=ICORR(:,ceil(size(I,1)/2);NSR=NPOW/IPOW;subplot(233);imshow(BlurredNoisy,);title(

14、9;A=Blurred and Noisy');subplot(234);imshow(deconvwnr(BlurredNoisy,PSF),);title('deconbwnr(A,PSF,NSR)');subplot(235);imshow(deconvwnr(BlurredNoisy,PSF,NCORR,ICORR),);title('deconbwnr(A,PSF,NCORR,ICORR)');subplot(236);imshow(deconvwnr(BlurredNoisy,PSF,NPOW,ICORR1),);title('dec

15、onbwnr(A,PSF,NPOW,ICORR_1_D)');(2) 正則濾波clcclear allclose allI=imread('C:Documents and SettingsAdministrator桌面原始圖.jpg');subplot(231);imshow(I);I=rgb2gray(I);subplot(232);imshow(I);PSF=fspecial('gaussian',7,10);V=.01;BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian',0,V);NOIS

16、EPOWER=V*prod(size(I);J LAGRA=deconvreg(BlurredNoisy,PSF,NOISEPOWER);subplot(233);imshow(BlurredNoisy);title('A=Blurred and Noisy');subplot(234);imshow(J);title('J LAGRA=deconvreg(A,PSF,NP)');subplot(235);imshow(deconvreg(BlurredNoisy,PSF,LAGRA/10);title('deconvreg(A,PSF,0.1*LAGR

17、A)');subplot(236);imshow(deconvreg(BlurredNoisy,PSF,LAGRA*10);title('deconvreg(A,PSF,10*LAGRA');(3) L-R算法clcclear allclose allI=imread('C:Documents and SettingsAdministrator桌面原始圖.jpg');subplot(231);imshow(I);I=rgb2gray(I);subplot(232);imshow(I);PSF=fspecial('gaussian',7,1

18、0);V=.0001;BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian',0,V);BlurredNoisy=double(BlurredNoisy);WT=zeros(size(I);WT(5:end-4,5:end-4)=1;J1=deconvlucy(BlurredNoisy,PSF);J2=deconvlucy(BlurredNoisy,PSF,20,sqrt(V);J3=deconvlucy(BlurredNoisy,PSF,20,sqrt(V),WT);subplot(233);imshow(mat2gray(Blurre

19、dNoisy);title('A=Blurred and Noisy');subplot(234);imshow(mat2gray(J1);title('deconvlucy(A,PSF)');subplot(235);imshow(mat2gray(J2);title('deconvlucy(A,PSF,NI,DP)');subplot(236);imshow(mat2gray(J3);title('deconvlucy(A,PSF,NI,DP,WT)');(4) 盲去卷積clcclear allclose allI=imread('C:Documents and SettingsAdministrator桌面原始圖.jpg');subplot(231);imshow(I);I=rgb2gray(I);subplot(232);imshow(I);PSF=fspecial('gaussian',7,10);V=.0001;BlurredNoisy=imnoise(imfilter(I,