An automated iterative algorithm for water and fat decomposition in three-point Dixon magnetic resonance imaging

An iterative, outlier exclusion, second-order surface fitting algorithm has been developed to solve the well-known phase wraparound problem associated with in vivo applications of the three-point Dixon magnetic resonance imaging method. The technique was optimized for speed by reducing the problem to a pair of planar fits. The spatial misalignment between water and fat components due to the chemical shift was handled on a subpixel level by invoking the shift theorem of Fourier transformation. From the chemical shift corrected water and fat images, high quality recombined MR images were generated. The algorithm was validated in both phantom and patient studies. In vivo breast images and pelvic images are provided as a demonstration of the method.