convol4.py

N = 15
kernel = np.ones((N, N))/N**2
fig, axes = plt.subplots(2, 2)
p_reflect  = PeriodicArray(scipy.ndimage.convolve(img, kernel, mode='reflect'))
p_wrap     = PeriodicArray(scipy.ndimage.convolve(img, kernel, mode='wrap'))
p_nearest  = PeriodicArray(scipy.ndimage.convolve(img, kernel, mode='nearest'))
p_constant = PeriodicArray(scipy.ndimage.convolve(img, kernel, mode='constant'))

p_reflect.plot(axes[0,0]); axes[0,0].set_title("mode='reflect'")
p_wrap.plot(axes[0,1]); axes[0,1].set_title("mode='wrap'")
p_nearest.plot(axes[1,0]); axes[1,0].set_title("mode='nearest'")
p_constant.plot(axes[1,1]); axes[1,1].set_title("mode='constant'")
fig.tight_layout()