matlab 高斯滤波(原创)

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;