维纳滤波和约束最小二乘滤波图像复原自编matlab代码

维纳滤波和约束最小二乘滤波图像复原自编matlab代码，共有两个文件CLSFilter.m，WienerFilter.m和一张测试图，可直接在R2013b上可以运行，有详细注释，注释里还有参考资料的网页链接，可帮助理解代码。

clear;clc
% 参考资料：http://blog.csdn.net/bluecol/article/details/47359421
% 参考数字图像处理（第三版） 冈萨雷斯 英文版公式（5.9-4）
% 约束最小二乘方滤波约束最小二乘方复原方法主要是针对有约束退化模型而言的，
% 这种图像复原方法只需有关噪声均值和方差的知识就能够对每一幅给定的图像进行
% 复原就以得到最优的效果。但有时该方程的解振荡的非常厉害。
I = im2double（imread（‘testimg.jpg‘））;
tic
[mn~] = size（I）;
subplot（131）imshow（I）;
title（‘原图‘）;


% Simulate a motion blur and additive noise.
LEN = 20;
THETA = 30;
PSF = fspecial（‘motion‘ LEN THETA）;
blurred = imfilter（I PSF ‘conv‘ ‘circular‘）;% 生成运动模糊图像

% psf2otf是把一个空间点扩散函数转换为频谱面的光学传递函数，其转换过程也是
% 一个二维的傅里叶变换。不过与fft2有区别
% psf2otf和fft2的区别：http://blog.csdn.net/bluecol/article/details/48288739
H = psf2otf（PSF[mn]）;% 退化函数H（uv）

noise_mean = 0;% 噪声均值
noise_var = 0.00001;% 噪声方差
blurred_noisy = imnoise（blurred‘gaussian‘...
    noise_mean noise_var）;% 加噪的运动模糊图像
subplot（132） imshow（blurred_noisy）
title（‘模糊退化且添加高斯噪声的图像‘）

% restoration using  Home Made Constrained Least Squares Filtering.
p = [0 -1 0;-1 4 -1;0 -1 0];% 拉普拉斯算子
P = psf2otf（p[mn]）;% 对拉普拉斯算子做傅立叶变换

gama = 0.001;% 公式里左边分式分母的一个参数
G = fft2（blurred_noisy）;% 对模糊图像做傅立叶变换

numerator = conj（H）;%公式里左边分式的分子，也就是H（uv）的共轭
denominator = H.^2 + gama*（P.^2）;%公式里左边分式的分母

deblurred2 = ifft2（numerator.*G./denominator）;% 再从频域变换回时域注意此处是点除
subplot（133） imshow（deblurred2）
title（‘约束最小二乘复原图‘）;
toc

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1646  2015-11-12 23:32  code\CLSFilter.m

     文件     263222  2013-11-16 21:11  code\testimg.jpg

     文件       1589  2015-11-12 23:33  code\WienerFilter.m

     目录          0  2015-11-12 23:39  code

----------- ---------  ---------- -----  ----

               266457                    4