This paper proposes a method for aligning image volumes acquired from different imaging modalities (e.g. MR, CT) based on 3D scale-invariant image features. A novel method for encoding invariant feature geometry and appearance is developed, based on the assumption of locally linear intensity relationships, providing a solution to poor repeatability of feature detection in different image modalities. The encoding method is incorporated into a probabilistic feature-based model for multi-modal image alignment. The model parameters are estimated via a group-wise alignment algorithm, that iteratively alternates between estimating a feature-based model from feature data, then realigning feature data to the model, converging to a stable alignment solution with few pre-processing or pre-alignment requirements. The resulting model can be used to align multi-modal image data with the benefits of invariant feature correspondence: globally optimal solutions, high efficiency and low memory usage. The method is tested on the difficult RIRE data set of CT, T1, T2, PD and MP-RAGE brain images of subjects exhibiting significant inter-subject variability due to pathology.
Publications by Year: 2013
Diffusion magnetic resonance imaging (dMRI) is an important tool that allows non-invasive investigation of the neural architecture of the brain. Advanced dMRI protocols typically require a large number of measurements for accurately tracing the fiber bundles and estimating the diffusion properties (such as, FA). However, the acquisition time of these sequences is prohibitively large for pediatric as well as patients with certain types of brain disorders (such as, dementia). Thus, fast echo-planar imaging (EPI) acquisition sequences were proposed by the authors in [6, 16], which acquired multiple slices simultaneously to reduce scan time. The scan time in such cases drops proportionately to the number of simultaneous slice acquisitions (which we denote by R). While preliminary results in [6, 16] showed good reproducibility, yet the effect of simultaneous acquisitions on long range fiber connectivity and diffusion measures such as FA, is not known. In this work, we use multi-tensor based fiber connectivity to compare data acquired on two subjects with different acceleration factors (R = 1, 2, 3). We investigate and report the reproducibility of fiber bundles and diffusion measures between these scans on two subjects with different spatial resolutions, which is quite useful while designing neuroimaging studies.
This paper reports on a new computational methodology, inter-slice correspondence (ISC), for robustly aligning sets of 2D ultrasound (US) slices during image-guided medical procedures. Correspondences are derived from distinctive, local scale-invariant features, which are used in one-to-many matching of US slices in near real-time despite out-of-plane rotation, in addition to global in-plane similarity transforms, occlusion, missing tissue, US plane mirroring, changes in US probe depth settings. Experiments demonstrate that ISC can align manually-acquired US slices without probe tracking information in the context of image-guided neurosurgery, with an accuracy of 1.3mm. A novel reconstruction-without-calibration application based on ISC is proposed, where 3D US reconstruction results are very similar to those obtained via traditional phantom-based calibration.
The average diffusion propagator (ADP) obtained from diffusion MRI (dMRI) data encapsulates important structural properties of the underlying tissue. Measures derived from the ADP can be potentially used as markers of tissue integrity in characterizing several mental disorders. Thus, accurate estimation of the ADP is imperative for its use in neuroimaging studies. In this work, we propose a simple method for estimating the ADP by representing the acquired diffusion signal in the entire q-space using radial basis functions (RBF). We demonstrate our technique using two different RBF’s (generalized inverse multiquadric and Gaussian) and derive analytical expressions for the corresponding ADP’s. We also derive expressions for computing the solid angle orientation distribution function (ODF) for each of the RBF’s. Estimation of the weights of the RBF’s is done by enforcing positivity constraint on the estimated ADP or ODF. Finally, we validate our method on data obtained from a physical phantom with known fiber crossing of 45 degrees and also show comparison with the solid spherical harmonics method of . We also demonstrate our method on in-vivo human brain data.