A method for sparse iterative phase correction for Magnetic Resonance (MR) partial Fourier reconstruction includes acquiring a partial Fourier k-space dataset using an MR scanner and estimating, by a control computer, a coil sensitivity map associated with the MR scanner from fully sampled k-space center. The control computer extracts a symmetrically sampled k-space center dataset from the partial Fourier k-space dataset and determines a low-resolution image based on the symmetrically sampled k-space center dataset and the coil sensitivity map. The control computer also determines phase corresponding to the low-resolution image. An iterative reconstruction process may then be applied to generate an image based on the partial Fourier k-space dataset. This iterative reconstruction process applies a Fast Iterative Shrinkage Thresholding Algorithm (FISTA) with phase correction based on the phase corresponding to the low-resolution image.