3D space-varying coefficient models with application to diffusion tensor imaging

  • Authors:

  • S. Heim;L. Fahrmeir;P. H. C. Eilers;B. D. Marx

  • Affiliations:

  • Department of Statistics, Ludwig-Maximilians-University, Ludwigstrasse 33/II, 80539 Munich, Germany;Department of Statistics, Ludwig-Maximilians-University, Ludwigstrasse 33/II, 80539 Munich, Germany;Department of Medical Statistics, Leiden University Medical Center, 2300 RA, Leiden, The Netherlands;Department of Experimental Statistics, Louisiana State University, Baton Rouge, LA 70803 USA

  • Venue:
  • Computational Statistics & Data Analysis
  • Year:
  • 2007

Quantified Score

Hi-index 0.03

Visualization

Abstract

The present methodological development and the primary application field originate from diffusion tensor imaging (DTI), a powerful nuclear magnetic resonance technique which enables the quantification of microscopical tissue properties. The current analysis framework of separate voxelwise regressions is reformulated as a 3D space-varying coefficient model (SVCM) for the entire set of diffusion tensor images recorded on a 3D voxel grid. The SVCM unifies the three-step cascade of standard data processing (voxelwise regression, smoothing, interpolation) into one framework based on B-spline basis functions. Thereby strength is borrowed from spatially correlated voxels to gain a regularization effect right at the estimation stage. Two SVCM variants are conceptualized: a full tensor product approach and a sequential approximation, rendering the SVCM numerically and computationally feasible even for the huge dimension of the joint model in a realistic setup. A simulation study shows that both approaches outperform the standard method of voxelwise regression with subsequent regularization. Application of the fast sequential method to real DTI data demonstrates the inherent ability to increase the grid resolution by evaluating the incorporated basis functions at intermediate points. The resulting continuous regularized tensor field may serve as basis for multiple applications, yet, ameloriation of local adaptivity is desirable.