We present a novel approach to shape-based interpolation of gray-level volume data. In contrast to the segmentation-based techniques our method directly processes the scalar volume requiring no user interaction. The key idea is to perform the interpolation in the directions given by analysis of the eigensystem of the structure tensor. Our method processes a 256