The optic radiation (OR) is one of the major components of the visual system and a key structure at risk in white matter diseases such as multiple sclerosis (MS). However, it is challenging to perform track reconstruction of the OR using diffusion MRI due to a sharp change of direction in the Meyer's loop and the presence of kissing and crossing fibers along the pathway. As such, we aimed to provide a highly precise and reproducible framework for tracking the OR from thalamic and visual cortex masks. The framework combined the generation of probabilistic streamlines by high order fiber orientation distributions estimated with constrained spherical deconvolution and an automatic post-processing based on anatomical exclusion criteria (AEC) to compensate for the presence of anatomically implausible streamlines. Specifically, those ending in the contralateral hemisphere, cerebrospinal fluid or grey matter outside the visual cortex were automatically excluded. We applied the framework to two distinct high angular resolution diffusion weighted imaging (HARDI) acquisition protocols on one cohort, comprised of ten healthy volunteers and five MS patients. The OR was successfully delineated in both HARDI acquisitions in the healthy volunteers and MS patients. Quantitative evaluation of the OR position was done by comparing the results with histological reference data. Compared with histological mask, the OR reconstruction into a template (OR-TCT) was highly precise (percentage of voxels within the OR-TCT correctly defined as OR), ranging from 0.71 to 0.83. The sensitivity (percentage of voxels in histological reference mask correctly defined as OR in OR-TCT) ranged from 0.65 to 0.81 and the accuracy (measured by F1 score) was 0.73 to 0.77 in healthy volunteers. When AEC was not applied the precision and accuracy decreased. The absolute agreement between both HARDI datasets measured by the intraclass correlation coefficient was 0.73. This improved framework allowed us to reconstruct the OR with high reliability and accuracy independently of the acquisition parameters. Moreover, the reconstruction was possible even in