Create a 3x3 matrix M = [1/3, -1/3, 0; 0, 1/3, -1/3; -1/3, 0, 1/3]. Create a column vector of line currents I = [IA; IB; IC]. Multiply M*I to get vector of phase currents [IAB; IBC; ICA].
The derivation of M was a little painful; it's the pseudoinverse of the non-invertible transformation...