Boundary coding for geodesic flows in higher dimensions