def np_gaussian_kernel_1d(sigma, window_radius=None): assert type(sigma) in (float, int) if window_radius is None: window_radius = int(np.ceil(3*sigma)) index = np.arange(2*window_radius + 1) - window_radius kernel = np.exp(-0.5*index**2/sigma**2) kernel = kernel / kernel.sum() return cast_to_floatx(kernel)
def np_gaussian_kernel_1d(sigma, window_radius=None): assert type(sigma) in (float, int) if window_radius is None: window_radius = int(np.ceil(3 * sigma)) index = np.arange(2 * window_radius + 1) - window_radius kernel = np.exp(-0.5 * index**2 / sigma**2) kernel = kernel / kernel.sum() return cast_to_floatx(kernel)
def astronaut(): astronaut = skimage.data.astronaut() / 255. astronaut = skimage.color.rgb2gray(astronaut) return cast_to_floatx(astronaut[::4, ::4])