Fast image reconstruction with L2-regularization

Purpose We introduce L2-regularized reconstruction algorithms with closed-form solutions that achieve dramatic computational speed-up relative to state of the art L1- and L2-based iterative algorithms while maintaining similar image quality for various applications in MRI reconstruction. Materials and Methods We compare fast L2-based methods to state of the art algorithms employing iterative L1- and L2-regularization in numerical phantom and in vivo data in three applications; (i) Fast Quantitative Susceptibility Mapping (QSM), (ii) Lipid artifact suppression in Magnetic Resonance Spectroscopic Imaging (MRSI), and (iii) Diffusion Spectrum Imaging (DSI). In all cases, proposed L2-based methods are compared with the state of the art algorithms, and two to three orders of magnitude speed up is demonstrated with similar reconstruction quality. Results The closed-form solution developed for regularized QSM allows processing of a three-dimensional volume under 5 s, the proposed lipid suppression algorithm takes under 1 s to reconstruct single-slice MRSI data, while the PCA based DSI algorithm estimates diffusion propagators from undersampled q-space for a single slice under 30 s, all running in Matlab using a standard workstation. Conclusion For the applications considered herein, closed-form L2-regularization can be a faster alternative to its iterative counterpart or L1-based iterative algorithms, without compromising image quality.National Institute for Biomedical Imaging and Bioengineering (U.S.) (Grant NIBIB K99EB012107)National Institutes of Health (U.S.) (Grant NIH R01 EB007942)National Institute for Biomedical Imaging and Bioengineering (U.S.) (Grant NIBIB R01EB006847)Grant K99/R00 EB008129National Center for Research Resources (U.S.) (Grant NCRR P41RR14075)National Institutes of Health (U.S.) (Blueprint for Neuroscience Research U01MH093765)Siemens CorporationSiemens-MIT AllianceMIT-Center for Integration of Medicine and Innovative Technology (Medical Engineering Fellowship

Denoising and fast diffusion imaging with physically constrained sparse dictionary learning

International audienceDiffusion-weighted imaging (DWI) allows imaging the geometry of water diffusion in biological tissues. However, DW images are noisy at high b-values and acquisitions are slow when using a large number of measurements, such as in Diffusion Spectrum Imaging (DSI). This work aims to denoise DWI and reduce the number of required measurements, while maintaining data quality. To capture the structure of DWI data, we use sparse dictionary learning constrained by the physical properties of the signal: symmetry and positivity. The method learns a dictionary of diffusion profiles on all the DW images at the same time and then scales to full brain data. Its performance is investigated with simulations and two real DSI datasets. We obtain better signal estimates from noisy measurements than by applying mirror symmetry through the q-space origin, Gaussian denoising or state-of- the-art non-local means denoising. Using a high-resolution dictionary learnt on another subject, we show that we can reduce the number of images acquired while still generating high resolution DSI data. Using dictionary learning, one can denoise DW images effectively and perform faster acquisitions. Higher b-value acquisitions and DSI techniques are possible with approximately 40 measurements. This opens important perspectives for the connectomics community using DSI

Probing holography in pp-adic CFT

We holographically calculate the partition functions of CFTs dual to Bruhat-Tits trees and $p$-adic BTZ black holes. Along the way, we propose new spectral decompositions of the Laplacian operator other than the plane-wave basis on these two types of backgrounds, with both analytical and numerical evidence. We extract the density of states and hence entropy from BTZ partition function via inverse Laplace transform. Then the one-loop Witten diagram is computed in the $p$-adic BTZ black hole background, yielding constraints on the heavy-heavy-light averaged three-point coefficient of its boundary $p$-adic CFT. Finally, for general $p$-adic CFTs (not necessarily holographic), we analyze the representation theory of their global conformal group $PGL\left(2,\mathbb{Q}_p\right)$, and discuss the suitability of different representations as Hilbert spaces of $p$-adic CFT.Comment: 52 pages, 13 figure

Compendio de métodos para caracterizar la geometría de los tejidos cerebrales a partir de imágenes de resonancia magnética por difusión del agua.

221 p.FIDMAG Hermanas Hospitalarias Research Foundation; CIBERSAM:Centro de Investigación Biomédica en Re

Large N Fields and Holography

We study (nearly) AdS/CFT holography within the context of the Sachdev-Ye- Kitaev (SYK) model. We present a systematic procedure to extract the dynamics of the low energy Schwarzian mode in SYK type models with a single energy scale J and emergent reparametrization symmetry in the infrared within the framework of perturbation theory. We develop a systematic approach using Feynman diagrams in bilocal theory to obtain a formal expression for the enhanced and O(1) corrections to the bilocal propagator and apply this general technique to large q SYK. We show that the Schwarzian theory describes a sector of a general class of two dimensional CFTs using conformal bootstrap techniques. We also provide a gravitational interpretation of this fact by showing that the dynamics in the near horizon throat of a specific class of BTZ black holes is described very well by Jackiw-Teitelboim (JT) gravity which gives rise to the Schwarzian action. We study the highly nontrivial conformal matter spectrum of the SYK model at arbitrary q. We provide a three dimensional bulk interpretation and carry out Kaluza-Klein (KK) reduction to reproduce the SYK spectrum and more significantly, the conformal bilocal propagator. We provide an interpretation of the two dimensional dual spacetime of the conformal sector of the SYK model using techniques from bilocal holography. We present a resolution of a conundrum about the signature of the dual spacetime using nonlocal integral transformations

Modélisation locale en imagerie par résonance magnétique de diffusion : de l'acquisition comprimée au connectome

L’imagerie par résonance magnétique pondérée en diffusion est une modalité d’imagerie médicale non invasive qui permet de mesurer les déplacements microscopiques des molécules d’eau dans les tissus biologiques. Il est possible d’utiliser cette information pour inférer la structure du cerveau. Les techniques de modélisation locale de la diffusion permettent de calculer l’orientation et la géométrie des tissus de la matière blanche. Cette thèse s’intéresse à l’optimisation des métaparamètres utilisés par les modèles locaux. Nous dérivons des paramètres optimaux qui améliorent la qualité des métriques de diffusion locale, de la tractographie de la matière blanche et de la connectivité globale. L’échantillonnage de l’espace-q est un des paramètres principaux qui limitent les types de modèle et d’inférence applicable sur des données acquises en clinique. Dans cette thèse, nous développons une technique d’échantillonnage de l’espace-q permettant d’utiliser l’acquisition comprimée pour réduire le temps d’acquisition nécessaire

Joint Spatial-Angular Sparse Coding, Compressed Sensing, and Dictionary Learning for Diffusion MRI

Neuroimaging provides a window into the inner workings of the human brain to diagnose and prevent neurological diseases and understand biological brain function, anatomy, and psychology. Diffusion Magnetic Resonance Imaging (dMRI) is an emerging medical imaging modality used to study the anatomical network of neurons in the brain, which form cohesive bundles, or fiber tracts, that connect various parts of the brain. Since about 73% of the brain is water, measuring the flow, or diffusion of water molecules in the presence of fiber bundles, allows researchers to estimate the orientation of fiber tracts and reconstruct the internal wiring of the brain, in vivo. Diffusion MRI signals can be modeled within two domains: the spatial domain consisting of voxels in a brain volume and the diffusion or angular domain, where fiber orientation is estimated in each voxel. Researchers aim to estimate the probability distribution of fiber orientation in every voxel of a brain volume in order to trace paths of fiber tracts from voxel to voxel over the entire brain. Therefore, the traditional framework for dMRI processing and analysis has been from a voxel-wise vantage point with added spatial regularization considered post-hoc. In contrast, we propose a new joint spatial-angular representation of dMRI data which pairs signals in each voxel with the global spatial environment, jointly. This has the ability to improve many aspects of dMRI processing and analysis and re-envision the core representation of dMRI data from a local perspective to a global one. In this thesis, we propose three main contributions which take advantage of such joint spatial-angular representations to improve major machine learning tasks applied to dMRI: sparse coding, compressed sensing, and dictionary learning. First, we will show that we can achieve sparser representations of dMRI by utilizing a global spatial-angular dictionary instead of a purely voxel-wise angular dictionary. As dMRI data is very large in size, we provide a number of novel extensions to popular spare coding algorithms that perform efficient optimization on a global-scale by exploiting the separability of our dictionaries over the spatial and angular domains. Next, compressed sensing is used to accelerate signal acquisition based on an underlying sparse representation of the data. We will show that our proposed representation has the potential to push the limits of the current state of scanner acceleration within a new compressed sensing model for dMRI. Finally, sparsity can be further increased by learning dictionaries directly from datasets of interest. Prior dictionary learning for dMRI learn angular dictionaries alone. Our third contribution is to learn spatial-angular dictionaries jointly from dMRI data directly to better represent the global structure. Traditionally, the problem of dictionary learning is non-convex with no guarantees of finding a globally optimal solution. We derive the first theoretical results of global optimality for this class of dictionary learning problems. We hope the core foundation of a joint spatial-angular representation will open a new perspective on dMRI with respect to many other processing tasks and analyses. In addition, our contributions are applicable to any general signal types that can benefit from separable dictionaries. We hope the contributions in this thesis may be adopted in the larger signal processing, computer vision, and machine learning communities. dMRI signals can be modeled within two domains: the spatial domain consisting of voxels in a brain volume and the diffusion or angular domain, where fiber orientation is estimated in each voxel. Computationally speaking, researchers aim to estimate the probability distribution of fiber orientation in every voxel of a brain volume in order to trace paths of fiber tracts from voxel to voxel over the entire brain. Therefore, the traditional framework for dMRI processing and analysis is from a voxel-wise, or angular, vantage point with post-hoc consideration of their local spatial neighborhoods. In contrast, we propose a new global spatial-angular representation of dMRI data which pairs signals in each voxel with the global spatial environment, jointly, to improve many aspects of dMRI processing and analysis, including the important need for accelerating the otherwise time-consuming acquisition of advanced dMRI protocols. In this thesis, we propose three main contributions which utilize our joint spatial-angular representation to improve major machine learning tasks applied to dMRI: sparse coding, compressed sensing, and dictionary learning. We will show that sparser codes are possible by utilizing a global dictionary instead of a voxel-wise angular dictionary. This allows for a reduction of the number of measurements needed to reconstruct a dMRI signal to increase acceleration using compressed sensing. Finally, instead of learning angular dictionaries alone, we learn spatial-angular dictionaries jointly from dMRI data directly to better represent the global structure. In addition, this problem is non-convex and so we derive the first theories to guarantee convergence to a global minimum. As dMRI data is very large in size, we provide a number of novel extensions to popular algorithms that perform efficient optimization on a global-scale by exploiting the separability of our global dictionaries over the spatial and angular domains. We hope the core foundation of a joint spatial-angular representation will open a new perspective on dMRI with respect to many other processing tasks and analyses. In addition, our contributions are applicable to any separable dictionary setting which we hope may be adopted in the larger image processing, computer vision, and machine learning communities