  Wavelength-resolved (WR) neutron transmission tomography is an emerging
technique to characterize engineering materials. While tomographic
reconstruction for amorphous samples is straightforward, it is challenging to
reconstruct samples with single-crystal domains because the attenuation of the
sample varies as a function of its orientation with respect to the incident
beam due to Bragg scattering. In this paper, we present an algorithm that can
reconstruct samples with single-crystal domains from WR neutron tomographic
measurements. In particular, we use a model-based iterative reconstruction
(MBIR) technique that reconstructs the volume by identifying and leaving out
the regions of the measurement that are affected by Bragg scatter. We combine
the output of the MBIR method with an algorithm that matches the reconstruction
to the identified Bragg scatter to reconstruct a feature that corresponds to
the local crystallography of the sample being measured. Using simulated data,
we demonstrate how our algorithm can reconstruct materials with single-crystal
domains, thereby adding a powerful new capability for WR neutron imaging
instruments.
