Microstructure evolution-based constitutive models in macro-scale simulation tools require statistical microstructure input at each integration point. This paper reports several robust procedures for interpolating orientation distribution functions (ODFs) from coarsely spaced experimental measurement grids to finely spaced modeling grids. The procedures are based on representing ODFs using generalized spherical harmonics (GSH) functions. Relying on the linearity of the expansion space, the core of the procedures involves weighting of the GSH expansion coefficients over a given variable such as space, strain, or temperature, thus providing ODFs as a function of location, deformation, or thermal treatment. Interpolation approaches including regression and linear/bilinear/barycentric weighting are explored and discussed in terms of their accuracy and suitability. Upon establishing the expansion coefficients of the interpolated ODF, a linear programing problem is solved within the expansion space to construct the interpolated ODF. The procedures are applied to a comprehensive data set obtained by neutron diffraction for a hemispherical part made of depleted α-uranium. Utilizing the ODF interpolation procedures, texture is interpolated over the part revealing significant anisotropy in the texture dependent thermal expansion and elastic stiffness coefficients. The results of this work facilitate introducing thermal expansion and elastic anisotropy into numerical tools for simulating thermo-mechanical loading of microstructurally heterogeneous components using the ODF informed crystal mechanics-based models to more accurately estimate the effective properties required by such simulations.