This paper presents a novel method for phase unwrapping for phase sensitive reconstruction in MR imaging. The unwrapped phase is obtained by integrating the phase gradient by solving a Poisson equation. An efficient solver, which has been made publicly available, is used to solve the equation. The proposed method is demonstrated on a fat quantification MRI task that is a part of a prospective study of fat accumulation. The method is compared to a phase unwrapping method based on region growing. Results indicate that the proposed method provides more robust unwrapping. Unlike region growing methods, the proposed method is also straight-forward to implement in 3D.