Quantitative MRI techniques as well as methods such as blood oxygen level-dependent (BOLD) imaging and in vivo spectroscopy require stringent optimization of magnetic field homogeneity, particularly when using high main magnetic fields. Automated shimming approaches require a method of measuring the main magnetic field, B 0 , followed by adjusting the currents in resistive shim coils to maximize homogeneity. A robust automated shimming technique using arbitrary mapping acquisition parameters (RASTAMAP) using a 3D multiecho gradient echo sequence that measures B 0 with high precision was developed. Inherent compensation and postprocessing methods enable removal of artifacts due to hardware timing errors, gradient propagation delays, gradient amplifier asymmetry, and eddy currents. This allows field maps to be generated for any field of view, bandwidth, resolution, or acquisition orientation without custom tuning of sequence parameters. The homogeneity of the static magnetic field, B 0 , is critical for many fast, quantitative, and spectroscopic imaging techniques. The trend towards higher static magnetic fields, both clinically and experimentally, is tempered by the increased magnetic field distortions caused by susceptibility differences. One method to compensate for these distortions is using higher-order resistive shim coils. However, the practical optimization of numerous higher-order resistive shim currents is only possible with automated techniques. Manual shimming by monitoring integrated signal magnitude and line shape during repeated adjustment of individual shim currents according to operator judgment or some mathematical algorithm (1) requires prohibitive amounts of time. Iterative readjustment of shim currents is required due to shim coil cross-terms resulting from either imperfect coils or off-center localization. In addition, signal weighting due to variations in longitudinal and transverse relaxation rates can cause inappropriate spatial biasing of these FID envelope methods.Automated methods of improving homogeneity developed to alleviate these inherent difficulties fall into two classes, projection mapping (2,3) and volume mapping (4,5). Projection mapping is advantageous for its short acquisition time. However, it relies on the incorrect assumption that shim coil fields are always fully characterized by a minimal set of spherical harmonics. In general, projection mapping also involves localization techniques that make it ill-suited for disjoint regions such as in multivoxel spectroscopy. Volume field mapping eliminates the reliance on spherical harmonics by using full 3D maps of the magnetic field generated by the shim coils to determine the optimum current settings. Regions of signal void do not affect the shim current calculations, and arbitrary, potentially disjoint regions over which to shim can be specified. This additional flexibility comes at the expense of increased time, since chemical shift imaging (6,7) and phase mapping (8,9) techniques that have been proposed to create volume B 0 ma...