常用图像去模糊算法分析与对比.doc

2012-2013 学年度第二学期数字图像处理课程设计题目: 常用图像去模糊算法分析与对比院系: 自动化学院专业: 自动化班级: 控实 1001 班姓名: 朱延俊学号: U201014240 指导教师: 谭山报告时间: 由于光学系统的装配、拍摄对焦不准或拍摄时的移动等均会使图像模糊,对诸如红外成像系统等离散成像系统,探测单元的非点元性质是图像模糊的重要原因之一,这些模糊图像一般可以用卷积过程描述。图像复原的基本问题是去模糊,即依据某种模糊和噪声降质的观测来估计原来的图像。反卷积是使模糊的图像复原的基本方法,如果成像系统的模糊函数已知,则去模糊成为常规的反卷积问题,否则,是盲解卷积问题。一般的,模糊函数是一个低通滤波器,使输入图像的高频成分收到抑制甚至丧失。反卷积是逆过程,需要恢复低频信息并找回丢失的高频成分。噪声的存在将可能使反卷积的结果偏离真实的解,因此需要在图像信号复原和噪声放大之间做出适当的折中。图像退化/复原可以用图 0所示的模型来表示( g= Hf+n): 图0 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