matlab高斯滤波器

来源：互联网发布：珠海金控知乎编辑：程序博客网时间：2024/05/13 19:11

1、使用matlab自带的函数

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear;img = imread('lena.tif');w = fspecial('gaussian',[5,5],1);%replicate:图像大小通过赋值外边界的值来扩展%symmetric 图像大小通过沿自身的边界进行镜像映射扩展I11 = imfilter(img,w,'replicate');figure(1);imshow(img);title('原图像');figure(2);imshow(I11);title('matlab高斯滤波后的图像');

结果：可以看到帽子和头发都非常显著的平滑了

2、自己编写的高斯滤波函数

%own gaussian function I = imread('lena.tif');I = double(I);%convert unit8 to double myimg = I;[rows, cols] = size(I);size = 5;conv = zeros(size,size);%高斯核segma = 1;sum = 0;center  = 3;for i=1:sizefor j=1:size%%利用高斯函数计算模板的各个参数temp = exp ( -((i-center)*(i-center) + (j-center)*(j-center)) / 2*segma*segma );conv(i,j) = temp / (2*3.14*segma*segma);sum =sum +  conv(i,j);endendconv = conv ./ sum;disp('kernel');display(conv);%use gaussian kernel to smooth imagefor k=1:rowsfor m= 1:cols    sum =0;%%卷积：模板与图像对应部分相乘再累加for i=1:sizefor j=1:sizeif( ( (k-center+i) > 0 && (k-center+i) <= rows ) && ((m-center+j)>0 && (m-center+j)<=cols) )sum = sum + conv(i,j) * I(k-center+i, m-center+j);endendend%%替换掉原像素myimg(k,m) = sum;endendfigure(3);imshow(myimg,[]);title('Gaussion Image');

结论：使用自己编写的滤波，虽然也起到了平滑效果，但是图像亮度偏高些，还不知道原因，希望大神指点指点。

阅读全文

0 0