The classic paradigm for MRI requires a homogeneous B 0 field in combination with linear encoding gradients. Distortions are produced when the B 0 is not homogeneous, and several postprocessing techniques have been developed to correct them. Field homogeneity is difficult to achieve, particularly for shortbore magnets and higher B 0 fields. Nonlinear magnetic components can also arise from concomitant fields, particularly in low-field imaging, or intentionally used for nonlinear encoding. In any of these situations, the second-order component is key, because it constitutes the first step to approximate higher-order fields. We propose to use the fractional Fourier transform for analyzing and reconstructing the object's magnetization under the presence of quadratic fields. The fractional fourier transform provides a precise theoretical framework for this. We show how it can be used for reconstruction and for gaining a better understanding of the quadratic field-induced distortions The standard paradigm for MRI requires a strong magnetic field with uniform intensity and time-varying linear encoding gradients across the entire field of view. However, deviations in the main field are common as uniform fields are physically difficult to achieve and also because of off-resonance effects from susceptibility changes. Such frequency variations introduce an accumulating phase over time, which cannot be demodulated easily as it varies spatially. This problem is worse for stronger fields and for sequences with long acquisition times. Great efforts are put in building the system with the highest possible homogeneity, for example, with passive or active shimming which partially correct first-and second-order field variations.Additionally, there has been an increasing interest in spatial encoding by nonhomogeneous, nonbijective spatial encoding magnetic fields (SEMs; Ref. 1). Starting from a general nonlinear field concept, novel techniques including PatLoc (1), O-Space (2), Null Space (3), and PhaseScrambled imaging (4-6) have all introduced the use of second-order SEMs as the first and simplest approach in simulations, custom-built hardware and experiments (7-12). One of the challenges of these approaches is to have an appropriate reconstruction technique.Higher-order fields also appear as a natural component of linear gradients. These concomitant fields can be well approximated by quadratic functions (13,14). In several applications, the magnitude of these fields is not negligible and artifacts appear in image reconstructions (15,16), most notably in very low-field or microtesla imaging (17,18).Several image reconstruction methods have been proposed to correct distortions, or to reconstruct an image produced by nonhomogeneous fields, being an active field of research (19-28). There is a well-known theoretical background for the linear correction approaches, in which an exact analytical solution is provided (19,23). For secondorder and arbitrary field maps, there are several correction techniques that typically trade-off corr...