Phase imaging is especially important at high to ultra high field strengths, both to map local field inhomogeneities and to highlight important anatomical structures (e.g. with Susceptibility Weighted Imaging, SWI). Although nowadays multiple coils are commonly used for MR image acquisition, accurate phase combination is not yet generally offered by the generic reconstruction pipeline of commercial scanners. The reason is that the various receiver coils experience different phase sensitivities due to their geometry, and the complex data need to be manipulated before combination in order to avoid destructive interferences. Last year, a generic framework was presented by Santini et al.¹ that allowed phase reconstruction for arbitrary sequences without any user interaction or prior knowledge, termed Generic Referenceless Phase Combination (GRPC). While the solution was optimal in a mathematical sense, the most robust embodiment was implementing a constant phase correction, therefore not taking into account spatial variations. In this work, we use the zeroth-order GRPC method as a first step of the correction, and then apply an additional correction as proposed by Parker et al.². The method is demonstrated in the brain and in the abdomen using two different acquisition sequences and the quality of the reconstruction of the two steps is compared by using a quality of phase matching metric³.The two-step GRPC correction was implemented as a generic functor of the reconstruction pipeline of a commercial 3T scanner by separately performing two correction step. The first step was a zeroth-order correction (removal of average coil phase) for preliminary coil rephasing. Subsequently, the rephased images were processed through a modified Parker virtual coil approach² by filtering the difference between the true coils and the rephased complex sum with a 21x21 2D gaussian low-pass filter. The reconstruction was applied to two different sequences, a gradient echo and a balanced steady-state free precession. Images were acquired on a healthy volunteer in the abdomen with a 18-channel body array together with an integrated spine coil (matrix size 256x154x1, resolution 1.6x1.6x5mm³, GRE TR/TE 15/3.6ms, bSSFP TR/TE 4.2/2.1ms), and in the head with 20- and a 64-channel coils (matrix size 192x192x1, resolution 1.0x1.0x5mm³, GRE TR/TE 45/4.7ms, bSSFP TR/TE 4.6/2.3ms). In the reconstructed complex images, the quality of phase matching (Q) map was calculated pixelwise according to the formula: $$Q(x,y)=\frac{||\sum I_c(x,y)||}{\sum||I_c(x,y)||}$$, where I_c is the complex signal intensity detected by each coil c. The median value of Q was calculated on each map (only in the portion of the images where the signal level is above the noise).The correction produced artifact-free phase images for all location and coils (Fig. 1). The phase images produced by GRPC and 2-step correction appear similar, however a lower noise level is recognizable in the two-pass images (Fig. 2). Quantitatively, the Quality of phase matching maps (Q-maps) showed important differences between the two methods (Fig. 3). The median Q values across the different coils and position averaged 0.56±0.19 for GRPC and 0.95±0.06 for the two-step correction (Fig. 4).For the combination of multi-element coils, both the GRPC and the proposed method deliver consistent phase images. However, the addition of the second step greatly increases the phase matching with which the coil elements sum, meaning that the coil sum constructively over the whole field of view. This increases the signal level in the combined image and has a noticeable impact on the signal-to-noise of the phase images. This method was used to implement a very generic reconstruction functor, which could be applied to any combination of protocol parameters and delivered fast reconstruction seamlessly at the scanner console, without need for external postprocessing. We showed that the method works equally well, and without operator interaction, with different coils and type of acquisition sequences.The addition of a second processing step to the GRPC method is beneficial to delivering accurate phase images from multi-element coils.