Reproducible Research in Signal Processing

Reproducible research results become more and more an important issue as systems under investigation are growing permanently in complexity, and it becomes thus almost impossible to judge the accuracy of research results merely on the bare paper presentation.Peer ReviewedPreprin

Variability and anatomical specificity of the orbitofrontothalamic fibers of passage in the ventral capsule/ventral striatum (VC/VS): precision care for patient-specific tractography-guided targeting of deep brain stimulation (DBS) in obsessive compulsive disorder (OCD)

Deep Brain Stimulation (DBS) is a neurosurgical procedure that can reduce symptoms in medically intractable obsessive-compulsive disorder (OCD). Conceptually, DBS of the ventral capsule/ventral striatum (VC/VS) region targets reciprocal excitatory connections between the orbitofrontal cortex (OFC) and thalamus, decreasing abnormal reverberant activity within the OFC-caudate-pallidal-thalamic circuit. In this study, we investigated these connections using diffusion magnetic resonance imaging (dMRI) on human connectome datasets of twenty-nine healthy young-adult volunteers with two-tensor unscented Kalman filter based tractography. We studied the morphology of the lateral and medial orbitofrontothalamic connections and estimated their topographic variability within the VC/VS region. Our results showed that the morphology of the individual orbitofrontothalamic fibers of passage in the VC/VS region is complex and inter-individual variability in their topography is high. We applied this method to an example OCD patient case who underwent DBS surgery, formulating an initial proof of concept for a tractography-guided patient-specific approach in DBS for medically intractable OCD. This may improve on current surgical practice, which involves implanting all patients at identical stereotactic coordinates within the VC/VS region

Imagerie de diffusion en temps-réel (correction du bruit et inférence de la connectivité cérébrale)

La plupart des constructeurs de systèmes d'imagerie par résonance magnétique (IRM) proposent un large choix d'applications de post-traitement sur les données IRM reconstruites a posteriori, mais très peu de ces applications peuvent être exécutées en temps réel pendant l'examen. Mises à part certaines solutions dédiées à l'IRM fonctionnelle permettant des expériences relativement simples ainsi que d'autres solutions pour l'IRM interventionnelle produisant des scans anatomiques pendant un acte de chirurgie, aucun outil n'a été développé pour l'IRM pondérée en diffusion (IRMd). Cependant, comme les examens d'IRMd sont extrêmement sensibles à des perturbations du système hardware ou à des perturbations provoquées par le sujet et qui induisent des données corrompues, il peut être intéressant d'investiguer la possibilité de reconstruire les données d'IRMd directement lors de l'examen. Cette thèse est dédiée à ce projet innovant. La contribution majeure de cette thèse a consisté en des solutions de débruitage des données d'IRMd en temps réel. En effet, le signal pondéré en diffusion peut être corrompu par un niveau élevé de bruit qui n'est plus gaussien, mais ricien ou chi non centré. Après avoir réalisé un état de l'art détaillé de la littérature sur le bruit en IRM, nous avons étendu l'estimateur linéaire qui minimise l'erreur quadratique moyenne (LMMSE) et nous l'avons adapté à notre cadre de temps réel réalisé avec un filtre de Kalman. Nous avons comparé les performances de cette solution à celles d'un filtrage gaussien standard, difficile à implémenter car il nécessite une modification de la chaîne de reconstruction pour y être inséré immédiatement après la démodulation du signal acquis dans l'espace de Fourier. Nous avons aussi développé un filtre de Kalman parallèle qui permet d'appréhender toute distribution de bruit et nous avons montré que ses performances étaient comparables à celles de notre méthode précédente utilisant un filtre de Kalman non parallèle. Enfin, nous avons investigué la faisabilité de réaliser une tractographie en temps-réel pour déterminer la connectivité structurelle en direct, pendant l'examen. Nous espérons que ce panel de développements méthodologiques permettra d'améliorer et d'accélérer le diagnostic en cas d'urgence pour vérifier l'état des faisceaux de fibres de la substance blanche.Most magnetic resonance imaging (MRI) system manufacturers propose a huge set of software applications to post-process the reconstructed MRI data a posteriori, but few of them can run in real-time during the ongoing scan. To our knowledge, apart from solutions dedicated to functional MRI allowing relatively simple experiments or for interventional MRI to perform anatomical scans during surgery, no tool has been developed in the field of diffusion-weighted MRI (dMRI). However, because dMRI scans are extremely sensitive to lots of hardware or subject-based perturbations inducing corrupted data, it can be interesting to investigate the possibility of processing dMRI data directly during the ongoing scan and this thesis is dedicated to this challenging topic. The major contribution of this thesis aimed at providing solutions to denoise dMRI data in real-time. Indeed, the diffusion-weighted signal may be corrupted by a significant level of noise which is not Gaussian anymore, but Rician or noncentral chi. After making a detailed review of the literature, we extended the linear minimum mean square error (LMMSE) estimator and adapted it to our real-time framework with a Kalman filter. We compared its efficiency to the standard Gaussian filtering, difficult to implement, as it requires a modification of the reconstruction pipeline to insert the filter immediately after the demodulation of the acquired signal in the Fourier space. We also developed a parallel Kalman filter to deal with any noise distribution and we showed that its efficiency was quite comparable to the non parallel Kalman filter approach. Last, we addressed the feasibility of performing tractography in real-time in order to infer the structural connectivity online. We hope that this set of methodological developments will help improving and accelerating a diagnosis in case of emergency to check the integrity of white matter fiber bundles.PARIS11-SCD-Bib. électronique (914719901) / SudocSudocFranceF

Fiber consistency measures on brain tracts from digital streamline, stochastic and global tractography

La tractografía es el proceso que se emplea para estimar la estructura de las fibras nerviosas del interior del cerebro in vivo a partir de datos de Resonancia Magnética (MR). Existen varios métodos de tractografía, que generalmente se dividen en locales y globales. Los primeros intentan reconstruir cada fibra por separado, mientras que los segundos intentan reconstruir todas las estructuras neuronales a la vez, buscando una configuración que mejor se ajusta a los datos proporcionados. Dichos métodos globales han demostrado ser más precisos y fiables que los métodos de tractografía local, para datos sintéticos. Sin embargo hasta la fecha no hay estudios que definan la relación entre los parámetros de adquisición de la MR y los resultados de tractografía estocástica o global con datos reales. Esta tésis de Master pretende mostrar la influencia de ciertos parámetros de adquisición como el factor de difusión de las secuencias de adquisición, el espaciado entre voxels o el número de gradientes en la variabilidad de las tractografías obtenidas.Teoría de la Señal, Comunicaciones e Ingeniería TelemáticaMáster en Investigación en Tecnologías de la Información y las Comunicacione

Higher-Order Tensors and Differential Topology in Diffusion MRI Modeling and Visualization

Diffusion Weighted Magnetic Resonance Imaging (DW-MRI) is a noninvasive method for creating three-dimensional scans of the human brain. It originated mostly in the 1970s and started its use in clinical applications in the 1980s. Due to its low risk and relatively high image quality it proved to be an indispensable tool for studying medical conditions as well as for general scientific research. For example, it allows to map fiber bundles, the major neuronal pathways through the brain. But all evaluation of scanned data depends on mathematical signal models that describe the raw signal output and map it to biologically more meaningful values. And here we find the most potential for improvement. In this thesis we first present a new multi-tensor kurtosis signal model for DW-MRI. That means it can detect multiple overlapping fiber bundles and map them to a set of tensors. Compared to other already widely used multi-tensor models, we also add higher order kurtosis terms to each fiber. This gives a more detailed quantification of fibers. These additional values can also be estimated by the Diffusion Kurtosis Imaging (DKI) method, but we show that these values are drastically affected by fiber crossings in DKI, whereas our model handles them as intrinsic properties of fiber bundles. This reduces the effects of fiber crossings and allows a more direct examination of fibers. Next, we take a closer look at spherical deconvolution. It can be seen as a generalization of multi-fiber signal models to a continuous distribution of fiber directions. To this approach we introduce a novel mathematical constraint. We show, that state-of-the-art methods for estimating the fiber distribution become more robust and gain accuracy when enforcing our constraint. Additionally, in the context of our own deconvolution scheme, it is algebraically equivalent to enforcing that the signal can be decomposed into fibers. This means, tractography and other methods that depend on identifying a discrete set of fiber directions greatly benefit from our constraint. Our third major contribution to DW-MRI deals with macroscopic structures of fiber bundle geometry. In recent years the question emerged, whether or not, crossing bundles form two-dimensional surfaces inside the brain. Although not completely obvious, there is a mathematical obstacle coming from differential topology, that prevents general tangential planes spanned by fiber directions at each point to be connected into consistent surfaces. Research into how well this constraint is fulfilled in our brain is hindered by the high precision and complexity needed by previous evaluation methods. This is why we present a drastically simpler method that negates the need for precisely finding fiber directions and instead only depends on the simple diffusion tensor method (DTI). We then use our new method to explore and improve streamsurface visualization.<br /

Kalman-filter-based EEG source localization

This thesis uses the Kalman filter (KF) to solve the electroencephalographic (EEG) inverse problem to image its neuronal sources. Chapter 1 introduces EEG source localization and the KF and discusses how it can solve the inverse problem. Chapter 2 introduces an EEG inverse solution using a spatially whitened KF (SWKF) to reduce the computational burden. Likelihood maximization is used to fit spatially uniform neural model parameters to simulated and clinical EEGs. The SWKF accurately reconstructs source dynamics. Filter performance is analyzed by computing the innovations’ statistical properties and identifying spatial variations in performance that could be improved by use of spatially varying parameters. Chapter 3 investigates the SWKF via one-dimensional (1D) simulations. Motivated by Chapter 2, two model parameters are given Gaussian spatial profiles to better reflect brain dynamics. Constrained optimization ensures estimated parameters have clear biophysical interpretations. Inverse solutions are also computed using the optimal linear KF. Both filters produce accurate state estimates. Spatially varying parameters are correctly identified from datasets with transient dynamics, but estimates for driven datasets are degraded by the unmodeled drive term. Chapter 4 treats the whole-brain EEG inverse problem and applies features of the 1D simulations to the SWKF of Chapter 2. Spatially varying parameters are used to model spatial variation of the alpha rhythm. The simulated EEG here exhibits wave-like patterns and spatially varying dynamics. As in Chapter 3, optimization constrains model parameters to appropriate ranges. State estimation is again reliable for simulated and clinical EEG, although spatially varying parameters do not improve accuracy and parameter estimation is unreliable, with wave velocity underestimated. Contributing factors are identified and approaches to overcome them are discussed. Chapter 5 summarizes the main findings and outlines future work

Kalman-filter-based EEG source localization

This thesis uses the Kalman filter (KF) to solve the electroencephalographic (EEG) inverse problem to image its neuronal sources. Chapter 1 introduces EEG source localization and the KF and discusses how it can solve the inverse problem. Chapter 2 introduces an EEG inverse solution using a spatially whitened KF (SWKF) to reduce the computational burden. Likelihood maximization is used to fit spatially uniform neural model parameters to simulated and clinical EEGs. The SWKF accurately reconstructs source dynamics. Filter performance is analyzed by computing the innovations’ statistical properties and identifying spatial variations in performance that could be improved by use of spatially varying parameters. Chapter 3 investigates the SWKF via one-dimensional (1D) simulations. Motivated by Chapter 2, two model parameters are given Gaussian spatial profiles to better reflect brain dynamics. Constrained optimization ensures estimated parameters have clear biophysical interpretations. Inverse solutions are also computed using the optimal linear KF. Both filters produce accurate state estimates. Spatially varying parameters are correctly identified from datasets with transient dynamics, but estimates for driven datasets are degraded by the unmodeled drive term. Chapter 4 treats the whole-brain EEG inverse problem and applies features of the 1D simulations to the SWKF of Chapter 2. Spatially varying parameters are used to model spatial variation of the alpha rhythm. The simulated EEG here exhibits wave-like patterns and spatially varying dynamics. As in Chapter 3, optimization constrains model parameters to appropriate ranges. State estimation is again reliable for simulated and clinical EEG, although spatially varying parameters do not improve accuracy and parameter estimation is unreliable, with wave velocity underestimated. Contributing factors are identified and approaches to overcome them are discussed. Chapter 5 summarizes the main findings and outlines future work

Generative models of brain connectivity for population studies

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2012.Cataloged from PDF version of thesis.Includes bibliographical references (p. 131-139).Connectivity analysis focuses on the interaction between brain regions. Such relationships inform us about patterns of neural communication and may enhance our understanding of neurological disorders. This thesis proposes a generative framework that uses anatomical and functional connectivity information to find impairments within a clinical population. Anatomical connectivity is measured via Diffusion Weighted Imaging (DWI), and functional connectivity is assessed using resting-state functional Magnetic Resonance Imaging (fMRI). We first develop a probabilistic model to merge information from DWI tractography and resting-state fMRI correlations. Our formulation captures the interaction between hidden templates of anatomical and functional connectivity within the brain. We also present an intuitive extension to population studies and demonstrate that our model learns predictive differences between a control and a schizophrenia population. Furthermore, combining the two modalities yields better results than considering each one in isolation. Although our joint model identifies widespread connectivity patterns influenced by a neurological disorder, the results are difficult to interpret and integrate with our regioncentric knowledge of the brain. To alleviate this problem, we present a novel approach to identify regions associated with the disorder based on connectivity information. Specifically, we assume that impairments of the disorder localize to a small subset of brain regions, which we call disease foci, and affect neural communication to/from these regions. This allows us to aggregate pairwise connectivity changes into a region-based representation of the disease. Once again, we use a probabilistic formulation: latent variables specify a template organization of the brain, which we indirectly observe through resting-state fMRI correlations and DWI tractography. Our inference algorithm simultaneously identifies both the afflicted regions and the network of aberrant functional connectivity. Finally, we extend the region-based model to include multiple collections of foci, which we call disease clusters. Preliminary results suggest that as the number of clusters increases, the refined model explains progressively more of the functional differences between the populations.by Archana Venkataraman.Ph.D