The past decade has seen rapid advances in the development of diffusion MRI techniques that permit for the first time the noninvasive characterization of neural architecture. While the exact biophysical determinants of the diffusion signal have yet to be completely elucidated, it is now generally accepted that microscopic boundaries to diffusion in the brain coincide with the local orientations of white matter (WM) fiber tracts (1). The advent of diffusion tensor imaging (DTI) as a tool for modeling intravoxel diffusion has inspired a number of promising applications in which WM connectivity can be evaluated in both health and disease (2,3). Measures derived from the tensor are now widely used to characterize regional anisotropy and orientation of WM throughout the brain, and more recently have been incorporated into tractography algorithms to allow 3D delineation of fiber pathways (4 -6).A major shortcoming of DTI lies in its inability to accurately characterize diffusion in complex WM, where fiber tracts with different orientations intersect or are otherwise partial volume averaged within a voxel (7-9). This limitation presents a significant obstacle to routine clinical application of DTI. For example, in brain tumor patients, presurgical DTI tractography fails to delineate the course of the pyramidal tract in regions of fiber crossing (10). Although intravoxel diffusion in these areas of complex WM can be more accurately characterized using the qspace formalism (11), long experiment times and heavy gradient demands preclude routine in vivo application of this technique in humans. High angular resolution diffusion imaging (HARDI) represents an alternative approach that strives to improve imaging efficiency by recovering the angular structure of diffusion in lieu of the 3D spindisplacement probability function. Shorter imaging times make HARDI particularly promising for clinical applications, especially when combined with the greater signalto-noise ratios (SNRs) afforded by high-field MR systems (12) and the reduction in single-shot echo-planar image artifacts provided by parallel imaging (12)(13)(14).Several approaches for reconstructing fiber orientations from HARDI data have been proposed, including spherical harmonic modeling of the apparent diffusion coefficient (ADC) profile (15,16), multitensor modeling (17), generalized tensor representations (18,19), circular spectrum mapping (20), q-ball imaging (QBI) (21-24), persistent angular structure (25), and spherical deconvolution (26). Among these techniques, QBI and spherical deconvolution have generated considerable interest due to their linearity and sensitivity to multimodal diffusion. QBI has the additional desirable property of model independence because it does not require ad hoc measurement of a "response function" in order to reconstruct fiber orientations. While the two methods differ significantly in their approach to recovering angular information from the measured diffusion data, both define a continuous function over the sphere that encodes the ...