python实现高斯滤波_如何在python中获得高斯滤波器

一般来说,如果您真的关心如何获得与MATLAB完全相同的结果,最简单的方法就是直接查看MATLAB函数的源代码.

在这种情况下,编辑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('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]])