python高斯滤波和降噪_如何在python中获得高斯滤波器

一般来说，如果你真的想得到与MATLAB完全相同的结果，最简单的方法就是直接查看MATLAB函数的源代码。

在这种情况下，edit fspecial：...

case 'gaussian' % Gaussian filter

siz = (p2-1)/2;

std = p3;

[x,y] = meshgrid(-siz(2):siz(2),-siz(1):siz(1));

arg = -(x.*x + y.*y)/(2*std*std);

h = exp(arg);

h(h

sumh = sum(h(:));

if sumh ~= 0,

h = h/sumh;

end;

...

很简单，嗯？将它移植到Python需要10分钟：import numpy as np

def matlab_style_gauss2D(shape=(3,3),sigma=0.5):

"""

2D gaussian mask - should give the same result as MATLAB's

fspecial('gaussian',[shape],[sigma])

"""

m,n = [(ss-1.)/2. for ss in shape]

y,x = np.ogrid[-m:m+1,-n:n+1]

h = np.exp( -(x*x + y*y) / (2.*sigma*sigma) )

h[ h < np.finfo(h.dtype).eps*h.max() ] = 0

sumh = h.sum()

if sumh != 0:

h /= sumh

return h

这与fspecial对舍入误差的回答相同：>> fspecial('gaussian',5,1)

0.002969 0.013306 0.021938 0.013306 0.002969

0.013306 0.059634 0.09832 0.059634 0.013306

0.021938 0.09832 0.1621 0.09832 0.021938

0.013306 0.059634 0.09832 0.059634 0.013306

0.002969 0.013306 0.021938 0.013306 0.002969

: matlab_style_gauss2D((5,5),1)

array([[ 0.002969, 0.013306, 0.021938, 0.013306, 0.002969],

[ 0.013306, 0.059634, 0.09832 , 0.059634, 0.013306],

[ 0.021938, 0.09832 , 0.162103, 0.09832 , 0.021938],

[ 0.013306, 0.059634, 0.09832 , 0.059634, 0.013306],

[ 0.002969, 0.013306, 0.021938, 0.013306, 0.002969]])