Fast Optimal Transport Averaging of Neuroimaging Data

Knowing how the Human brain is anatomically and functionally organized at the level of a group of healthy individuals or patients is the primary goal of neuroimaging research. Yet computing an average of brain imaging data defined over a voxel grid or a triangulation remains a challenge. Data are large, the geometry of the brain is complex and the between subjects variability leads to spatially or temporally non-overlapping effects of interest. To address the problem of variability, data are commonly smoothed before group linear averaging. In this work we build on ideas originally introduced by Kantorovich to propose a new algorithm that can average efficiently non-normalized data defined over arbitrary discrete domains using transportation metrics. We show how Kantorovich means can be linked to Wasserstein barycenters in order to take advantage of an entropic smoothing approach. It leads to a smooth convex optimization problem and an algorithm with strong convergence guarantees. We illustrate the versatility of this tool and its empirical behavior on functional neuroimaging data, functional MRI and magnetoencephalography (MEG) source estimates, defined on voxel grids and triangulations of the folded cortical surface.Comment: Information Processing in Medical Imaging (IPMI), Jun 2015, Isle of Skye, United Kingdom. Springer, 201

Photon-assisted tunneling in a Fe8 Single-Molecule Magnet

The low temperature spin dynamics of a Fe8 Single-Molecule Magnet was studied under circularly polarized electromagnetic radiation allowing us to establish clearly photon-assisted tunneling. This effect, while linear at low power, becomes highly non-linear above a relatively low power threshold. This non-linearity is attributed to the nature of the coupling of the sample to the thermostat.These results are of great importance if such systems are to be used as quantum computers.Comment: 4 pages, 4 figure

Bayesian Joint Detection-Estimation of cerebral vasoreactivity from ASL fMRI data

International audienceAlthough the study of cerebral vasoreactivity using fMRI is mainly conducted through the BOLD fMRI modality, owing to its relatively high signal-to-noise ratio (SNR), ASL fMRI provides a more interpretable measure of cerebral vasoreactivity than BOLD fMRI. Still, ASL suffers from a low SNR and is hampered by a large amount of physiological noise. The current contribution aims at improving the re- covery of the vasoreactive component from the ASL signal. To this end, a Bayesian hierarchical model is proposed, enabling the recovery of per- fusion levels as well as fitting their dynamics. On a single-subject ASL real data set involving perfusion changes induced by hypercapnia, the approach is compared with a classical GLM-based analysis. A better goodness-of-fit is achieved, especially in the transitions between baseline and hypercapnia periods. Also, perfusion levels are recovered with higher sensitivity and show a better contrast between gray- and white matter

Unsupervised Fiber Bundles Registration using Weighted Measures Geometric Demons

International audienceBrain image registration aims at reducing anatomical variability across subjects to create a common space for group analysis. Multi-modal approaches intend to minimize cortex shape variations along with internal structures, such as fiber bundles. A di ficulty is that it requires a prior identi fication of these structures, which remains a challenging task in the absence of a complete reference atlas. We propose an extension of the log-Geometric Demons for jointly registering images and fi ber bundles without the need of point or ber correspondences. By representing fi ber bundles as Weighted Measures we can register subjects with di fferent numbers of fiber bundles. The ef ficacy of our algorithm is demonstrated by registering simultaneously T1 images and between 37 and 88 ber bundles depending on each of the ten subject used. We compare results with a multi-modal T1 + Fractional Anisotropy (FA) and a tensor-based registration algorithms and obtain superior performance with our approach

Magnetic Anisotropy of a Single Cobalt Nanoparticle

Using a new microSQUID set-up, we investigate magnetic anisotropy in a single 1000-atoms cobalt cluster. This system opens new fields in the characterization and the understanding of the origin of magnetic anisotropy in such nanoparticles. For this purpose, we report three-dimensional switching field measurements performed on a 3 nm cobalt cluster embedded in a niobium matrix. We are able to separate the different magnetic anisotropy contributions and evidence the dominating role of the cluster surface.Comment: 4 pages, 8 figure

Joint T1 and Brain Fiber Log-Demons Registration Using Currents to Model Geometry

International audienceWe present an extension of the diffeomorphic Geometric Demons algorithm which combines the iconic registration with geometric constraints. Our algorithm works in the log-domain space, so that one can efficiently compute the deformation field of the geometry. We represent the shape of objects of interest in the space of currents which is sensitive to both location and geometric structure of objects. Currents provides a distance between geometric structures that can be defined without specifying explicit point-to-point correspondences. We demonstrate this framework by registering simultaneously T1 images and 65 fiber bundles consistently extracted in 12 subjects and compare it against non-linear T1, tensor, and multi-modal T1+ Fractional Anisotropy (FA) registration algorithms. Results show the superiority of the Log-domain Geometric Demons over their purely iconic counterparts