Matlab中进行高斯滤波-学习笔记
最近在学习的时候遇到了问题,稍稍总结一下。
比如说我现在要对一幅图像进行高斯滤波,要求SIGMA=1,卷积核为5×5。网上查找了一下,总共分为三种办法。
fspecial + filter2
B = imread('jessica.bmp');
core = fspecial('gaussian',[5,5],1);
C = filter2(core,B);
fspecial + imfilter
B = imread('jessica.bmp');
core = fspecial('gaussian',[5,5],1);
D = imfilter(B,core);
直接使用imgaussfilt
B = imread('jessica.bmp');
E = imgaussfilt(B,1);
从结果中发现,这三种滤波方式中,前两种的结果相同,而第三种和前两种不同。
查询帮助文档,发现原因是imgaussfilt的默认边界填充为replicate,而并非其他两种方法中边界的默认值zero。
故若将第三种方法改成
B = imread('jessica.bmp');
E = imgaussfilt(B,1,‘Padding’,0);
就可以得到和上面相同的结果。