The “primary” or “first-order relay” nuclei of the thalamus feed the cerebral cortex with information about ongoing activity in the environment or the subcortical motor systems. Because of the small size of these nuclei and the high specificity of their input and output pathways, new imaging protocols are required to investigate thalamocortical interactions in human perception, cognition and language. The goal of the present study was twofold: I) to develop a reconstruction protocol based on in vivo diffusion MRI to extract and measure the axonal fiber tracts that originate or terminate specifically in individual first-order relay nuclei; and, II) to test the reliability of this reconstruction protocol. In left and right hemispheres, we investigated the thalamocortical/corticothalamic axon bundles linking each of the first-order relay nuclei and their main cortical target areas, namely, the lateral geniculate nucleus (optic radiation), the medial geniculate nucleus (acoustic radiation), the ventral posterior nucleus (somatosensory radiation) and the ventral lateral nucleus (motor radiation). In addition, we examined the main subcortical input pathway to the ventral lateral posterior nucleus, which originates in the dentate nucleus of the cerebellum. Our protocol comprised three components: defining regions-of-interest; preprocessing diffusion data; and modeling white-matter tracts and tractometry. We then used computation and test-retest methods to check whether our protocol could reliably reconstruct these tracts of interest and their profiles. Our results demonstrated that the protocol had nearly perfect computational reproducibility and good-to-excellent test-retest reproducibility. This new protocol may be of interest for both basic human brain neuroscience and clinical studies and has been made publicly available to the scientific community.