Correction of susceptibility artifacts in diffusion tensor data using non-linear registration

Diffusion tensor imaging can be used to localize major white matter tracts within the human brain. For surgery of tumors near eloquent brain areas such as the pyramidal tract this information is of importance to achieve an optimal resection while avoiding post-operative neurological deficits. However, due to the small bandwidth of echo planar imaging, diffusion tensor images suffer from susceptibility artifacts resulting in positional shifts and distortion. As a consequence, the fiber tracts computed from echo planar imaging data are spatially distorted. We present an approach based on non-linear registration using Bézier functions to efficiently correct distortions due to susceptibility artifacts. The approach makes extensive use of graphics hardware to accelerate the non-linear registration procedure. An improvement presented in this paper is a more robust and efficient optimization strategy based on simultaneous perturbation stochastic approximation (SPSA). Since the accuracy of non-linear registration is crucial for the value of the presented correction method, two techniques were applied in order to prove the quality of the proposed framework. First, the registration accuracy was evaluated by recovering a known transformation with non-linear registration. Second, landmark-based evaluation of the registration method for anatomical and diffusion tensor data was performed. The registration was then applied to patients with lesions adjacent to the pyramidal tract in order to compensate for susceptibility artifacts. The effect of the correction on the pyramidal tract was then quantified by measuring the position of the tract before and after registration. As a result, the distortions observed in phase encoding direction were most prominent at the cortex and the brainstem. The presented approach allows correcting fiber tract distortions which is an important prerequisite when tractography data are integrated into a stereotactic setup for intra-operative guidance.