常用图像去模糊算法分析与对比综述

2012-2013学年度第二学期

数字图像处理课程设计

题目： 常用图像去模糊算法分析与对比

院 系： 自动化学院 专 业： 自动化 班 级： 控实1001班 姓 名： 朱延俊 学 号： U201014240 指导教师： 谭山 报告时间： 2013.6.24

由于光学系统的装配、拍摄对焦不准或拍摄时的移动等均会使图像模糊，对诸如红外成像系统等离散成像系统，探测单元的非点元性质是图像模糊的重要原因之一，这些模糊图像一般可以用卷积过程描述。

图像复原的基本问题是去模糊，即依据某种模糊和噪声降质的观测来估计原来的图像。反卷积是使模糊的图像复原的基本方法，如果成像系统的模糊函数已知，则去模糊成为常规的反卷积问题，否则，是盲解卷积问题。一般的，模糊函数是一个低通滤波器，使输入图像的高频成分收到抑制甚至丧失。反卷积是逆过程，需要恢复低频信息并找回丢失的高频成分。噪声的存在将可能使反卷积的结果偏离真实的解，因此需要在图像信号复原和噪声放大之间做出适当的折中。

图像退化/复原可以用图0所示的模型来表示（g = Hf + n）：

图0

f(x,y)原图像H(x,y)退化函数(PSF)g(x,y)退化函数n(x,y)附加噪声f’(x,y)复原的图像M(x,y)复原函数 1.盲目解卷积算法（Blind Deconvolution Algorithm）

图像恢复是一种改善图像质量的处理技术, 将降质了的图像恢复成原来的图像。目前, 图像恢复的方法很多, 然而在图像恢复过程中, 最难解决的问题之一是如何获得恢复算法中PSF的恰当估计, 那些不以PSF 知识为基础的图像恢复方法统称为盲去卷积算法。盲去卷积的方法已经受到了人们的极大重视, 对于给定的原图像, 使其退化, 得到退化图像, 再利用盲去卷积的方法使其恢复, 得到视觉质量上更好的图像。

盲解卷积的方法是以最大似然估计(MLE )为基础的，即一种用被随机噪声所干扰的量进行估计的最优化策略。似然函数用g ( x, y )、f ( x, y ) 和h( x, y )来加以表达, 然后问题就变成了寻求最大似然函数。在盲解卷积中, 最优化问题用规定的约束条件并假定收敛时通过迭代来求解, 得到的最大f( x, y )和h (x, y )就是还原的图像和PSF。

【函数】deconvblind

【功能】使用盲解卷积算法对图像进行去模糊

[J,PSF] = DECONVBLIND(I,INITPSF) deconvolves image I using maximum likelihood algorithm, returning both deblurred image J and a restored point-spread function PSF. The resulting PSF is a positive array of the same size as the INITPSF, normalized so its sum adds to 1. The PSF

restoration is affected strongly by the size of its initial guess, INITPSF, and less by its values (an array of ones is a safer guess).

使用盲解卷积对图像I进行去模糊，得到去模糊后的图像J和重建点扩散函数矩阵PSF。参量INITPSF为矩阵，表示重建点扩散函数矩阵的初始值。

[J,PSF] = DECONVBLIND(I,INITPSF,NUMIT)

参量NUMIT为迭代次数，默认值为10。

[J,PSF] = DECONVBLIND(I,INITPSF,NUMIT,DAMPAR)

参量DAMPAR表示输出图像与输入图像的偏离阈值，该函数对于偏离阈值的像素不再进行迭代计算，这抑制了像素上的噪声，又保存了图像的细节。

[J,PSF] = DECONVBLIND(I,INITPSF,NUMIT,DAMPAR,WEIGHT)

参量WEIGHT为矩阵，其元素为图像每个像素的权值，默认值为与输入图像相同维数的单位矩阵。

[J,PSF] = DECONVBLIND(I,INITPSF,NUMIT,DAMPAR,WEIGHT,READOUT)

参量READOUT制定噪声类型，默认值为0。 【编程实现】

I = checkerboard(8);%创建棋盘：为8或者30 PSF = fspecial('gaussian',7,10);%设置滤波器 V = .0001;

BlurredNoisy = imnoise(imfilter(I,PSF),'gaussian',0,V);%添加高斯噪声 %迭代次数为20，输出图像与输入图像的偏离阈值为10*sqrt(V)，阈值为zeros(size(I) WT = zeros(size(I)); WT(5:end-4,5:end-4) = 1; INITPSF = ones(size(PSF));

[J P] = deconvblind(BlurredNoisy,INITPSF,20,10*sqrt(V),WT); %对图像去模糊，显示原始图像和PSF图像 subplot(221);imshow(BlurredNoisy); title('Blurred and Noisy'); subplot(222);imshow(PSF,[]); title('True PSF'); %复原图像和PSF图像 subplot(223);imshow(J);

title('Deblurred Image'); subplot(224);imshow(P,[]);

title('Recovered PSF');

得到的图形如图1所示，可以看出同时恢复了图像和点扩张函数,在对失真情况毫无先验知识的情况下，仍能实现对模糊图像的恢复操作。利用MATLAB 实现的图像恢复, 并对恢复图像的失真情况做了改善。在进行图像恢复时，重建PSF，对图像进行重建, 得到恢复的图像。若I = checkerboard(30);则效果如图2所示。

图1

图2

2.L-R算法（Lucy-Richardson Algorithm）

L-R算法假设图像服从Possion分布，采用最大似然估计进行估算，其迭代方程为

f(x,y)n?1?f(x,y)n[(g(x,y))?h(x,y)] nh(x,y)*f(x,y)该算法考虑到了物体的先验信息，无需模糊图像的噪声信息，适用于线性和非线性成像模型，解具有唯一性。在MATLAB中，L-R算法附加了一些改进，可以在多次迭代过于逼近噪声易产生错误信息斑点的情况下，使用参数来控制收敛次数，使每次过程中结果图像和原始图像的每个相应局部的背离程度的灰度值在阈值内。

【函数】deconvlucy

【功能】使用L-R算法对图像进行去模糊处理

J = DECONVLUCY(I,PSF)

deconvolves image I using Lucy-Richardson algorithm, returning deblurred image J. The assumption is that the image I was created by convolving a true image with a point-spread function PSF and possibly by adding noise. J = DECONVLUCY(I,PSF,NUMIT) J = DECONVLUCY(I,PSF,NUMIT,DAMPAR) J = DECONVLUCY(I,PSF,NUMIT,DAMPAR,WEIGHT) J = DECONVLUCY(I,PSF,NUMIT,DAMPAR,WEIGHT,READOUT)

J = DECONVLUCY(I,PSF,NUMIT,DAMPAR,WEIGHT,READOUT,SUBSMPL), where

【编程实现】

I = checkerboard(8);%创建棋盘，8或30 PSF = fspecial('gaussian',7,10);%设置滤波器 V = .0001;

BlurredNoisy = imnoise(imfilter(I,PSF),'gaussian',0,V);%添加高斯噪声 WT = zeros(size(I));%阈值为zeros(size(I)

WT(5:end-4,5:end-4) = 1;%重建点扩散函数矩阵的初始值 J1 = deconvlucy(BlurredNoisy,PSF);

%设置输出图像与输入图像的偏离阈值为sqrt(V)

J2 = deconvlucy(BlurredNoisy,PSF,20,sqrt(V)); %设置阈值为zeros(size(I)

J3 = deconvlucy(BlurredNoisy,PSF,20,sqrt(V),WT);

subplot(221);imshow(BlurredNoisy); title('Blurred and Noisy'); subplot(222);imshow(J1);

title('deconvlucy(A,PSF)'); subplot(223);imshow(J2);

title('deconvlucy(A,PSF,NI,DP …… 此处隐藏：2252字，全部文档内容请下载后查看。喜欢就下载吧 ……