example.py 文件源码

def kernel(s):
    """
    Constructs a normalized discrete 1D gaussian kernel
    """
    size_grid = int(s*4)
    x= scipy.mgrid[-size_grid:size_grid+1]
    g = scipy.exp(-(x**2/float(s**2)/2.))
    return g / np.sum(g)