If y want to convlute your function with a gaussian in multidimensional space,
scipy.ndimage.gaussian_filter is useful:
from scipy.ndimage import gaussian_filter
data3_blur = gaussian_filter(data3, sigma=[1, 0.5, 0])
sigma is a scalar or list of scalars for respective axes.
