Finite element-based invariant-domain preserving approximation of hyperbolic systems: Beyond second-order accuracy in space