C[(i), (j), :, :] = np.array([a[i - 1, j - 1], a[i - 1, j], a[i - 1, j + 1],
    a[i, j - 1], a[i, j], a[i, j + 1], a[i + 1, j - 1], a[i + 1, j], a[i + 
    1, j + 1]])
Разделяемый фильтр в массиве numpy