1. A method of nuclear magnetic resonance imaging of an object, the method comprising: receiving MR data corresponding to n complex MR images generated using an MR scan having a series of n different echo times, where n is an integer equal to two or more, each image corresponding to a respective echo time; and processing the images to correct phase wrap aliasing, the processing comprising: i) complex dividing a second image corresponding to the second echo time with a first image corresponding to the first echo time to create a divided complex image; ii) using a scaling of the phase from the divided complex image to determine a predicted phase at one or more pixels in the second image; iii) for one or more pixels in the second image, comparing the predicted phase with a measured phase, and if the difference between the predicted phase and the measured phase lies within a first threshold range, leaving the measured phase in the second image unchanged; if the difference between the phases lies outside the first threshold range, changing the measured phase value by adding multiples of 2π until the difference falls within the threshold range.