A magnetic resonance imaging method includes generating spatially resolved fiber orientation distributions (FODs) from magnetic resonance signals acquired from a patient tissue using a plurality of diffusion encodings, each acquired magnetic resonance signal corresponding to one of the diffusion encodings and being representative of a three-dimensional distribution of displacement of magnetic spins of gyromagnetic nuclei present in each imaging voxel. Generating the spatially resolved FODs includes performing generalized spherical deconvolution using the acquired magnetic resonance signals and a modeled tissue response matrix (TRM) to reconstruct the spatially resolved FODs. The method also includes using the spatially resolved FODs to generate a representation of fibrous tissue within the patient tissue.