Method of MR imaging of an object placed in a main magnetic field within an examination volume of a MR device, the method comprising the steps of: - subjecting the object to an imaging sequence to generate two sets of echo signals at two different echo times using bipolar pairs of readout magnetic field gradients, one set of echo signals being generated at a first echo time and the other set of echo signals being generated at a second echo time, - acquiring the echo signals from the object, - reconstructing a first image from the echo signals attributed to the first echo time and a second image from the echo signals attributed to the second echo time, - computing modified first and second images by compensating for a fat shift in the reconstructed first and second images respectively, which compensation of the fat shift involves spatially shifting the first and second images in the positive and negative readout direction, depending on the respective polarity of the readout magnetic field gradient, whereby a spatial displacement of signal contributions from fat in the reconstructed first and second images is corrected for in the modified first and second images, - estimating phase errors in the acquired echo signals on the basis of the first and second images and the modified first and second images by a two-point Dixon method, wherein the phase errors are estimated in two individual computation steps using the first and second images in one computation step and using the modified first and second images in another computation step, wherein, for each voxel position, the phase error estimated on the basis of the first and second images is selected if the contribution from water is found to dominate at the respective voxel position and the phase error estimated on the basis of the modified first and second images is selected if the contribution from fat is found to dominate at the respective voxel position, and - reconstructing a water image and/or a fat image by separating the signal contributions of fat and water to the acquired echo signals using the estimated phase errors. - subjecting the object to an imaging sequence to generate two sets of echo signals at two different echo times using bipolar pairs of readout magnetic field gradients, one set of echo signals being generated at a first echo time and the other set of echo signals being generated at a second echo time, - acquiring the echo signals from the object, - reconstructing a first image from the echo signals attributed to the first echo time and a second image from the echo signals attributed to the second echo time, - computing modified first and second images by compensating for a fat shift in the reconstructed first and second images respectively, which compensation of the fat shift involves spatially shifting the first and second images in the positive and negative readout direction, depending on the respective polarity of the readout magnetic field gradient, whereby a spatial displacement of signal contributions from fat in the reconstructed first and second images is corrected for in the modified first and second images, - estimating phase errors in the acquired echo signals on the basis of the first and second images and the modified first and second images by a two-point Dixon method, wherein the phase errors are estimated in two individual computation steps using the first and second images in one computation step and using the modified first and second images in another computation step, wherein, for each voxel position, the phase error estimated on the basis of the first and second images is selected if the contribution from water is found to dominate at the respective voxel position and the phase error estimated on the basis of the modified first and second images is selected if the contribution from fat is found to dominate at the respective voxel position, and - reconstructing a water image and/or a fat image by separating the signal contributions of fat and water to the acquired echo signals using the estimated phase errors.