维纳滤波图像去噪

clear all; close all
I = imread（‘eight.tif‘）;
figure（1） imshow（I） title（‘Original Image‘）;

% 加入方差0.005的Gaussian 噪声
J1 = imnoise（I‘gaussian‘00.005）;
% 加入方差0.01的Gaussian 噪声
J2 = imnoise（I‘gaussian‘00.01）;
% 加入方差0.02的Gaussian 噪声
J3 = imnoise（I‘gaussian‘00.02）;
% 加入0.02的椒盐噪声
F1 = imnoise（I‘salt & pepper‘0.02）;
% 加入0.05的椒盐噪声
F2 = imnoise（I‘salt & pepper‘0.05）;
% 加入0.1的椒盐噪声
F3 = imnoise（I‘salt & pepper‘0.1）;

% 对J1、J2、J3进行模版为3的维纳滤波
K1 = wiener2（J1[3 3]）;
K2 = wiener2（J2[3 3]）;
K3 = wiener2（J3[3 3]）;
figure（2）
subplot（131） imshow（K1） title（‘J1 Filtered Image With 3*3 ‘）
subplot（132） imshow（K2） title（‘J2 Filtered Image With 3*3 ‘）
subplot（133） imshow（K3） title（‘J3 Filtered Image With 3*3 ‘）

% 对F1、F2、F3进行模版为3的维纳滤波
S1 = wiener2（F1[3 3]）;
S2 = wiener2（F2[3 3]）;
S3 = wiener2（F3[3 3]）;
figure（3）
subplot（131） imshow（S1） title（‘F1 Filtered Image With 3*3 ‘）
subplot（132） imshow（S2） title（‘F2 Filtered Image With 3*3 ‘）
subplot（133） imshow（S3） title（‘F3 Filtered Image With 3*3 ‘）

% 对J1、J2、J3进行模版为5的维纳滤波
K1_5 = wiener2（J1[5 5]）;
K2_5 = wiener2（J2[5 5]）;
K3_5 = wiener2（J3[5 5]）;
figure（4）
subplot（1