We present a novel method for high-order phase reduction in networks of weakly coupled oscillators and, more generally, perturbations of reducible normally hyperbolic (quasi-)periodic tori. Our method works by computing an asymptotic expansion for an embedding of the perturbed invariant torus, as well as for the reduced phase dynamics in local coordinates. Both can be determined to arbitrary degrees of accuracy, and we show that the phase dynamics may directly be obtained in normal form. We apply the method to predict remote synchronisation in a chain of coupled Stuart-Landau oscillators.