Use phase info in EPI images

Once EPI scan done, there are 2 DICOM folders. One with amplitude, another with phase values.
I heard that one can incorporate the phase information, in calculating average EPI image over many time points.

Is this correct? Anyone knows a way to do that? With AFNI or nipype, etc.