Robust and highly scalable estimation of directional couplings from time-shifted signals