Nonuniform DFT Applications in MRI: A Novel Algorithm and Parallel Implementation

Typically, the discrete Fourier transform (DFT) is performed on a sequence of data points that are equidistant on a given domain; however, it may be desirable to relax this constraint. In certain types of magnetic resonance imaging (MRI), there is the opportunity to take several unequally spaced samples simultaneously. In one MRI scheme known as “pure phase encoding”, each time a point is sampled in Fourier space it is possible to sample several nearby points very cheaply, which can improve the resulting image. Reconstructing an image from these points requires an inverse non-uniform DFT.