Probabilistic tractography using Q-ball imaging and particle filtering: application to adult and in-utero fetal brain studies.

By assuming that orientation information of brain white matter fibers can be inferred from Diffusion-Weighted Magnetic Resonance Imaging (DW-MRI) measurements, tractography algorithms provide an estimation of the brain connectivity in vivo. The two key ingredients of tractography are the diffusion model (tensor, high-order tensor, Q-ball, etc.) and the means to deal with uncertainty during the tracking process (deterministic vs probabilistic mathematical framework). In this paper, we investigate the use of an analytical Q-ball model for the diffusion data within a well-formalized particle filtering framework. The proposed method is validated and compared to other tracking algorithms on the MICCAI'09 contest Fiber Cup phantom. Tractographies of in vivo adult and fetal brain Diffusion-Weighted Images (DWIs) are also shown to illustrate the robustness of the algorithm.

Cardiac motion estimation by joint alignment of tagged MRI sequences.

Image registration has been proposed as an automatic method for recovering cardiac displacement fields from tagged Magnetic Resonance Imaging (tMRI) sequences. Initially performed as a set of pairwise registrations, these techniques have evolved to the use of 3D+t deformation models, requiring metrics of joint image alignment (JA). However, only linear combinations of cost functions defined with respect to the first frame have been used. In this paper, we have applied k-Nearest Neighbors Graphs (kNNG) estimators of the α-entropy (H(α)) to measure the joint similarity between frames, and to combine the information provided by different cardiac views in an unified metric. Experiments performed on six subjects showed a significantly higher accuracy (p<0.05) with respect to a standard pairwise alignment (PA) approach in terms of mean positional error and variance with respect to manually placed landmarks. The developed method was used to study strains in patients with myocardial infarction, showing a consistency between strain, infarction location, and coronary occlusion. This paper also presents an interesting clinical application of graph-based metric estimators, showing their value for solving practical problems found in medical imaging.

Wall motion estimation in intracranial aneurysms.

The quantification of wall motion in cerebral aneurysms is becoming important owing to its potential connection to rupture, and as a way to incorporate the effects of vascular compliance in computational fluid dynamics simulations. Most of papers report values obtained with experimental phantoms, simulated images or animal models, but the information for real patients is limited. In this paper, we have combined non-rigid registration with signal processing techniques to measure pulsation in real patients from high frame rate digital subtraction angiography. We have obtained physiological meaningful waveforms with amplitudes in the range 0 mm-0.3 mm for a population of 18 patients including ruptured and unruptured aneurysms. Statistically significant differences in pulsation were found according to the rupture status, in agreement with differences in biomechanical properties reported in the literature.

Efficient computational fluid dynamics mesh generation by image registration.

Most implementations of computational fluid dynamics (CFD) solutions require a discretisation or meshing of the solution domain. The production from a medical image of a computationally efficient mesh representing the structures of interest can be time consuming and labour-intensive, and remains a major bottleneck in the clinical application of CFD. This paper presents a method for deriving a patient-specific mesh from a medical image. The method uses volumetric registration of a pseudo-image, produced from an idealised template mesh, with the medical image. The registration algorithm used is robust and computationally efficient. The accuracy of the new algorithm is measured in terms of the distance between a registered surface and a known surface, for image data derived from casts of the lumen of two different vessels. The true surface is identified by laser profiling. The average distance between the surface points measured by the laser profiler and the surface of the mapped mesh is better than 0.2 mm. For the images analysed, the new algorithm is shown to be 2-3 times more accurate than a standard published algorithm based on maximising normalised mutual information. Computation times are approximately 18 times faster for the new algorithm than the standard algorithm. Examples of the use of the algorithm on two clinical examples are also given. The registration methodology lends itself immediately to the construction of dynamic mesh models in which vessel wall motion is obtained directly using registration.

Myocardial motion estimation in tagged MR sequences by using alphaMI-based non rigid registration.

Tagged Magnetic Resonance Imaging (MRI) is currently the reference MR modality for myocardial motion and strain analysis. NMI-based non rigid registration has proven to be an accurate method to retrieve cardiac deformation fields. The use of alphaMI permits higher dimensional features to be implemented in myocardial deformation estimation through image registration. This paper demonstrates that this is feasible with a set of Haar wavelet features of high dimension. While we do not demonstrate performance improvement for this set of features, there is no significant degradation as compared to implementing the registration method with the traditional NMI metric. We use Entropic Spanning Graphs (ESGs) to estimate the alphaMI of the wavelet feature vectors WFVs since this is not possible with histograms. To the best of our knowledge, this is the first time that ESGs are used for non rigid registration.