clear all;
close all;
tic;
P=input('input picture data:');
d=size(P);
if(d(3)>1)
P=rgb2gray(P);
end
subplot(131);
imshow(P);
title('原图像');
P=imnoise(P,'gaussian',0.1,0.002);
subplot(132);imshow(P);
title('加入高斯噪声后的图片');
P=double(P);
[m n]=size(P);
k=floor([(m+1)/2 (n+1)/2]);%计算图象中心
sigma=input('input smoothing parameter:');
b=zeros(m,n);
for i=1:m
for j=1:n
% b(i,j) =exp(-((i-k(1))^2+(j-k(2))^2/(2*sigma^2)));
b(i,j) =exp(-((i-k(1))^2+(j-k(2))^2)/(4*sigma))/(4*pi*sigma);
% b(i,j)=exp(-((i^2+j^2)/(2*sigma^2)));
end
end
Img1=conv2(P,b,'same'); %用生成的高斯序列卷积运算，进行高斯滤波
d=uint8(Img1);
subplot(133);imshow(d);
title('高斯滤波去噪结果');
toc;
因篇幅问题不能全部显示，请点此查看更多更全内容