Interpolation is an essential step in the visualization process. While most data from simulations or experiments are discrete many visualization methods are based on smooth, continuous data approximation or interpolation methods. We introduce a new interpolation method for symmetrical tensor fields given on a triangulated domain. Differently from standard tensor field interpolation, which is based on the tensor components, we use tensor invariants, eigenvectors and eigenvalues, for the interpolation. This interpolation minimizes the number of eigenvectors and eigenvalues computations by restricting it to mesh vertices and makes an exact integration of the tensor lines possible. The tensor field topology is qualitatively the same as for the component wise-interpolation. Since the interpolation decouples the ``shape'' and ``direction'' interpolation it is shape-preserving, what is especially important for tracing fibers in diffusion MRI data.