Constrained Spherical Deconvolution (CSD) estimates the Fibre Orientation Distribution (FOD) in each voxel in the Spherical Harmonic (SH) basis. Although features of the underlying fibre configuration are captured within the FOD, the SH basis is not wholly appropriate for their analysis. We propose a re-parametrization of the FOD, fitting a sparse set of SH delta functions to each FOD, that additionally captures the anisotropic dispersion of each discrete fibre bundle.