Multifrequency interpolation for fast off-resonance correction

Field inhomogeneities or susceptibility variations produce blurring in images acquired using non-2DFT k-space readout trajectories. This problem is more pronounced for sequences with long readout times such as spiral imaging. Theoretical and practical correction methods based on an acquired field map have been reported in the past. This paper introduces a new correction method based on the existing concept of frequency segmented correction but which is faster and the-oretically more accurate. It consists of reconstructing the data at several frequencies to form a set of base images that are then added together with spatially varying linear coefficients derived from the field map. The new algorithm is applied to phantom and in vivo images acquired with projection reconstruction and spiral sequences, yielding sharply focused images.