Kronnecker product on just two dimensions of 3 dimensional arrays

Question

Suppose two matrices:

a = np.array([[1,2],
              [3,4]])
b = np.array([[5,6],
              [7,8]])

which would give Kronecker produc