Proceedings of the second "international Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST'14)

The implicit objective of the biennial "international - Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST) is to foster collaboration between international scientific teams by disseminating ideas through both specific oral/poster presentations and free discussions. For its second edition, the iTWIST workshop took place in the medieval and picturesque town of Namur in Belgium, from Wednesday August 27th till Friday August 29th, 2014. The workshop was conveniently located in "The Arsenal" building within walking distance of both hotels and town center. iTWIST'14 has gathered about 70 international participants and has featured 9 invited talks, 10 oral presentations, and 14 posters on the following themes, all related to the theory, application and generalization of the "sparsity paradigm": Sparsity-driven data sensing and processing; Union of low dimensional subspaces; Beyond linear and convex inverse problem; Matrix/manifold/graph sensing/processing; Blind inverse problems and dictionary learning; Sparsity and computational neuroscience; Information theory, geometry and randomness; Complexity/accuracy tradeoffs in numerical methods; Sparsity? What's next?; Sparse machine learning and inference.Comment: 69 pages, 24 extended abstracts, iTWIST'14 website: http://sites.google.com/site/itwist1

3D single breath-hold MR methodology for measuring cardiac parametric mapping at 3T

Mención Internacional en el título de doctorOne of the foremost and challenging subfields of MRI is cardiac magnetic resonance imaging (CMR). CMR is becoming an indispensable tool in cardiovascular medicine by acquiring data about anatomy and function simultaneously. For instance, it allows the non-invasive characterization of myocardial tissues via parametric mapping techniques. These mapping techniques provide a spatial visualization of quantitative changes in the myocardial parameters. Inspired by the need to develop novel high-quality parametric sequences for 3T, this thesis's primary goal is to introduce an accurate and efficient 3D single breath-hold MR methodology for measuring cardiac parametric mapping at 3T. This thesis is divided into two main parts: i) research and development of a new 3D T1 saturation recovery mapping technique (3D SACORA), together with a feasibility study regarding the possibility of adding a T2 mapping feature to 3D SACORA concepts, and ii) research and implementation of a deep learning-based post-processing method to improve the T1 maps obtained with 3D SACORA. In the first part of the thesis, 3D SACORA was developed as a new 3D T1 mapping sequence to speed up T1 mapping acquisition of the whole heart. The proposed sequence was validated in phantoms against the gold standard technique IR-SE and in-vivo against the reference sequence 3D SASHA. The 3D SACORA pulse sequence design was focused on acquiring the entire left ventricle in a single breath-hold while achieving good quality T1 mapping and stability over a wide range of heart rates (HRs). The precision and accuracy of 3D SACORA were assessed in phantom experiments. Reference T1 values were obtained using IR-SE. In order to further validate 3D SACORA T1 estimation accuracy and precision, T1 values were also estimated using an in-house version of 3D SASHA. For in-vivo validation, seven large healthy pigs were scanned with 3D SACORA and 3D SASHA. In all pigs, images were acquired before and after administration of MR contrast agent. The phantom results showed good agreement and no significant bias between methods. In the in-vivo experiments, all T1-weighted images showed good contrast and quality, and the T1 maps correctly represented the information contained in the T1-weighted images. Septal T1s and coefficients of variation did not considerably differ between the two sequences, confirming good accuracy and precision. 3D SACORA images showed good contrast, homogeneity and were comparable to corresponding 3D SASHA images, despite the shorter acquisition time (15s vs. 188s, for a heart rate of 60 bpm). In conclusion, the proposed 3D SACORA successfully acquired a whole-heart 3D T1 map in a single breath-hold at 3T, estimating T1 values in agreement with those obtained with the IR-SE and 3D SASHA sequences. Following the successful validation of 3D SACORA, a feasibility study was performed to assess the potential of modifying the acquisition scheme of 3D SACORA in order to obtain T1 and T2 maps simultaneously in a single breath-hold. This 3D T1/T2 sequence was named 3D dual saturation-recovery compressed SENSE rapid acquisition (3D dual-SACORA). A phantom of eight tubes was built to validate the proposed sequence. The phantom was scanned with 3D dual-SACORA with a simulated heart rate of 60 bpm. Reference T1 and T2 values were estimated using IR-SE and GraSE sequences, respectively. An in-vivo study was performed with a healthy volunteer to evaluate the parametric maps' image quality obtained with the 3D dual-SACORA sequence. T1 and T2 maps of the phantom were successfully obtained with the 3D dual-SACORA sequence. The results show that the proposed sequence achieved good precision and accuracy for most values. A volunteer was successfully scanned with the proposed sequence (acquisition duration of approximately 20s) in a single breath-hold. The saturation time images and the parametric maps obtained with the 3D dual-SACORA sequence showed good contrast and homogeneity. The septal T1 and T2 values are in good agreement with reference sequences and published work. In conclusion, this feasibility study's findings open the door to the possibility of using 3D SACORA concepts to develop a successful 3D T1/T2 sequence. In the second part of the thesis, a deep learning-based super-resolution model was implemented to improve the image quality of the T1 maps of 3D SACORA, and a comprehensive study of the performance of the model in different MR image datasets and sequences was performed. After careful consideration, the selected convolutional neural network to improve the image quality of the T1 maps was the Residual Dense Network (RDN). This network has shown outstanding performance against state-of-the-art methods on benchmark datasets; however, it has not been validated on MR datasets. In this way, the RDN model was initially validated on cardiac and brain benchmark datasets. After this validation, the model was validated on a self-acquired cardiac dataset and on improving T1 maps. The RDN model improved the images successfully for the two benchmark datasets, achieving better performance with the brain dataset than with the cardiac dataset. This result was expected as the brain images have more well-defined edges than the cardiac images, making the resolution enhancement more evident. On the self-acquired cardiac dataset, the model also obtained an enhanced performance on image quality assessment metrics and improved visual assessment, particularly on well-defined edges. Regarding the T1 mapping sequences, the model improved the image quality of the saturation time images and the T1 maps. The model was able to enhance the T1 maps analytically and visually. Analytically, the model did not considerably modify the T1 values while improving the standard deviation in both myocardium and blood. Visually, the model improved the T1 maps by removing noise and motion artifacts without losing resolution on the edges. In conclusion, the RDN model was validated on three different MR datasets and used to improve the image quality of the T1 maps obtained with 3D SACORA and 3D SASHA. In summary, a 3D single breath-hold MR methodology was introduced, including a ready to-go 3D single breath-hold T1 mapping sequence for 3T (3D SACORA), together with the ideas for a new 3D T1/T2 mapping sequence (3D dual-SACORA); and a deep learning-based post-processing implementation capable of improving the image quality of 3D SACORA T1 maps.This thesis has received funding from the European Union Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement N722427.Programa de Doctorado en Multimedia y Comunicaciones por la Universidad Carlos III de Madrid y la Universidad Rey Juan CarlosPresidente: Carlos Alberola López.- Secretario: María Jesús Ledesma Carbayo.- Vocal: Nathan Mewto

Improving low-dose blood-brain barrier permeability quantification using sparse high-dose induced prior for Patlak model

Blood-brain barrier permeability (BBBP) measurements extracted from the perfusion computed tomography (PCT) using the Patlak model can be a valuable indicator to predict hemorrhagic transformation in patients with acute stroke. Unfortunately, the standard Patlak model based PCT requires excessive radiation exposure, which raised attention on radiation safety. Minimizing radiation dose is of high value in clinical practice but can degrade the image quality due to the introduced severe noise. The purpose of this work is to construct high quality BBBP maps from low-dose PCT data by using the brain structural similarity between different individuals and the relations between the high- and low-dose maps. The proposed sparse high-dose induced (shd-Patlak) model performs by building a high-dose induced prior for the Patlak model with a set of location adaptive dictionaries, followed by an optimized estimation of BBBP map with the prior regularized Patlak model. Evaluation with the simulated low-dose clinical brain PCT datasets clearly demonstrate that the shd-Patlak model can achieve more significant gains than the standard Patlak model with improved visual quality, higher fidelity to the gold standard and more accurate details for clinical analysis. Copyright 2013 Elsevier B.V. All rights reserved

CINENet: deep learning-based 3D cardiac CINE MRI reconstruction with multi-coil complex-valued 4D spatio-temporal convolutions

Cardiac CINE magnetic resonance imaging is the gold-standard for the assessment of cardiac function. Imaging accelerations have shown to enable 3D CINE with left ventricular (LV) coverage in a single breath-hold. However, 3D imaging remains limited to anisotropic resolution and long reconstruction times. Recently deep learning has shown promising results for computationally efficient reconstructions of highly accelerated 2D CINE imaging. In this work, we propose a novel 4D (3D + time) deep learning-based reconstruction network, termed 4D CINENet, for prospectively undersampled 3D Cartesian CINE imaging. CINENet is based on (3 + 1)D complex-valued spatio-temporal convolutions and multi-coil data processing. We trained and evaluated the proposed CINENet on in-house acquired 3D CINE data of 20 healthy subjects and 15 patients with suspected cardiovascular disease. The proposed CINENet network outperforms iterative reconstructions in visual image quality and contrast (+ 67% improvement). We found good agreement in LV function (bias ± 95% confidence) in terms of end-systolic volume (0 ± 3.3 ml), end-diastolic volume (- 0.4 ± 2.0 ml) and ejection fraction (0.1 ± 3.2%) compared to clinical gold-standard 2D CINE, enabling single breath-hold isotropic 3D CINE in less than 10 s scan and ~ 5 s reconstruction time

MRI reconstruction using Markov random field and total variation as composite prior

Reconstruction of magnetic resonance images (MRI) benefits from incorporating a priori knowledge about statistical dependencies among the representation coefficients. Recent results demonstrate that modeling intraband dependencies with Markov Random Field (MRF) models enable superior reconstructions compared to inter-scale models. In this paper, we develop a novel reconstruction method, which includes a composite prior based on an MRF model and Total Variation (TV). We use an anisotropic MRF model and propose an original data-driven method for the adaptive estimation of its parameters. From a Bayesian perspective, we define a new position-dependent type of regularization and derive a compact reconstruction algorithm with a novel soft-thresholding rule. Experimental results show the effectiveness of this method compared to the state of the art in the field

Advanced VLBI Imaging

Very Long Baseline Interferometry (VLBI) is an observational technique developed in astronomy for combining multiple radio telescopes into a single virtual instrument with an effective aperture reaching up to many thousand kilometers and enabling measurements at highest angular resolutions. The celebrated examples of applying VLBI to astrophysical studies include detailed, high-resolution images of the innermost parts of relativistic outflows (jets) in active galactic nuclei (AGN) and recent pioneering observations of the shadows of supermassive black holes (SMBH) in the center of our Galaxy and in the galaxy M87. Despite these and many other proven successes of VLBI, analysis and imaging of VLBI data still remain difficult, owing in part to the fact that VLBI imaging inherently constitutes an ill-posed inverse problem. Historically, this problem has been addressed in radio interferometry by the CLEAN algorithm, a matching-pursuit inverse modeling method developed in the early 1970-s and since then established as a de-facto standard approach for imaging VLBI data. In recent years, the constantly increasing demand for improving quality and fidelity of interferometric image reconstruction has resulted in several attempts to employ new approaches, such as forward modeling and Bayesian estimation, for application to VLBI imaging. While the current state-of-the-art forward modeling and Bayesian techniques may outperform CLEAN in terms of accuracy, resolution, robustness, and adaptability, they also tend to require more complex structure and longer computation times, and rely on extensive finetuning of a larger number of non-trivial hyperparameters. This leaves an ample room for further searches for potentially more effective imaging approaches and provides the main motivation for this dissertation and its particular focusing on the need to unify algorithmic frameworks and to study VLBI imaging from the perspective of inverse problems in general. In pursuit of this goal, and based on an extensive qualitative comparison of the existing methods, this dissertation comprises the development, testing, and first implementations of two novel concepts for improved interferometric image reconstruction. The concepts combine the known benefits of current forward modeling techniques, develop more automatic and less supervised algorithms for image reconstruction, and realize them within two different frameworks. The first framework unites multiscale imaging algorithms in the spirit of compressive sensing with a dictionary adapted to the uv-coverage and its defects (DoG-HiT, DoB-CLEAN). We extend this approach to dynamical imaging and polarimetric imaging. The core components of this framework are realized in a multidisciplinary and multipurpose software MrBeam, developed as part of this dissertation. The second framework employs a multiobjective genetic evolutionary algorithm (MOEA/D) for the purpose of achieving fully unsupervised image reconstruction and hyperparameter optimization. These new methods are shown to outperform the existing methods in various metrics such as angular resolution, structural sensitivity, and degree of supervision. We demonstrate the great potential of these new techniques with selected applications to frontline VLBI observations of AGN jets and SMBH. In addition to improving the quality and robustness of image reconstruction, DoG-HiT, DoB-CLEAN and MOEA/D also provide such novel capabilities as dynamic reconstruction of polarimetric images on minute time-scales, or near-real time and unsupervised data analysis (useful in particular for application to large imaging surveys). The techniques and software developed in this dissertation are of interest for a wider range of inverse problems as well. This includes such versatile fields such as Ly-alpha tomography (where we improve estimates of the thermal state of the intergalactic medium), the cosmographic search for dark matter (where we improve forecasted bounds on ultralight dilatons), medical imaging, and solar spectroscopy