反卷积复原算法.pdfVIP

几种图像复原方法的对比 一、 Richardson-Lucy 算法 R-L 算法是目前世界上应用最广泛的函数恢复技术之一，它是一 种迭代方法。 MATLAB提供的 deconvlucy （）函数还能够用于实现复 杂图像重建的多种算法中， 这些算法都基于 Lucy-Richardson 最大化 可能性算法。 R-L 算法是一种迭代非线性复原算法，它是从最大似然公式推导 出来的，图像用泊松分布加以模型化的。 当下面这个迭代收敛时模型 的最大似然函数就可以得到一个令人满意的方程： g (x , y ) fk 1(x , y ) f k (x, y)[ h(x, y)] h( x, y) f k (x, y) 其中， * 代表卷积， 代表相关， f 代表未退化图像的估计， g 和 h 和以前定义一样。在 IPT 中， L-R 算法由名为 deconvlucy 的函 数完成的。 deconvlucy() 函数的调用格式： J=deconvlucy(I ，PSF，NUMIT， 。 DAMPAR，WEIGHT) 其中， I 表示输入图像， PSF表示点扩散函数。其 他参数都是可选参数： NUMIT表示算法的迭代次数，默认为 10 次； DAMPAR是一个标量，它指定了结果图像与原图像 I 之间的偏离阈值 表，默认值为 0( 无衰减）；WEIGHT是一个与 I 同样大小的数组，它为 每一个像素分配一个权重来反映其重量， 表示像素加权值， 默认值为 原始图像的数值。 图像复原源代码： %%Deblurring Gray Images Using the Lucy-Richardson Algorithm clc clear close all I=imread(E:\lena512color.tif); % 彩色图像的像素为 512*512 I1=rgb2gray(I); % 灰度图像的像素为 512*512 % figure,imshow(I),title(Original color image); % figure,imshow(I1),title(Original gray image); I2=I1(1:2:end,1:2:end); % 图像的像素设置为 256*256 figure,imshow(I2),title(Gray Image 256*256); PSF = fspecial(gaussian,5,5); % 点扩散函数 Blurred = imfilter(I2,PSF,symmetric,conv); figure; imshow(Blurred); title(Gaussian Blurred); V = 0.0001; BlurredNoisy = imnoise(Blurred,gaussian,0,V); figure; imshow(BlurredNoisy); title(Blurred Noisy); K=size(I2); WT=zeros(K); WT(5:end-4,5:end-4)=1; J1 = deconvlucy(BlurredNoisy,PSF); % H1 = deconvlucy(BlurredNoisy,PSF,5);