The present disclosure provides a method and system for quantifying and making images of tissue anisotropy property based on magnetic resonance imaging (MRI). The systems and methods provided herein utilize orientation distribution of magnetic susceptibility to characterize magnetic susceptibility anisotropy (MSA) inside biological tissues. This MSA may be intrinsic property of the tissue or may be induced by the presence of external agents. In certain embodiments, the MSA is displayed as an orientation distribution function of susceptibility and/or may be described by mathematical quantities such as tensors (e.g., symmetric or asymmetric second order or higher order tensors) and spherical harmonics. In other embodiments, the MSA is characterized using a second order tensor named apparent susceptibility tensor (AST).