def hessian(array): (dy, dx) = np.gradient(array) (dydy, dxdy) = np.gradient(dy) (dydx, dxdx) = np.gradient(dx) return np.dstack((dxdx, dydx, dxdy, dydy))