Coupled Dictionary Learning for Multi-Contrast MRI Reconstruction.

Magnetic resonance (MR) imaging tasks often involve multiple contrasts, such as T1-weighted, T2-weighted and fluid-attenuated inversion recovery (FLAIR) data. These contrasts capture information associated with the same underlying anatomy and thus exhibit similarities in either structure level or gray level. In this paper, we propose a coupled dictionary learning based multi-contrast MRI reconstruction (CDLMRI) approach to leverage the dependency correlation between different contrasts for guided or joint reconstruction from their under-sampled k -space data. Our approach iterates between three stages: coupled dictionary learning, coupled sparse denoising, and enforcing k -space consistency. The first stage learns a set of dictionaries that not only are adaptive to the contrasts, but also capture correlations among multiple contrasts in a sparse transform domain. By capitalizing on the learned dictionaries, the second stage performs coupled sparse coding to remove the aliasing and noise in the corrupted contrasts. The third stage enforces consistency between the denoised contrasts and the measurements in the k -space domain. Numerical experiments, consisting of retrospective under-sampling of various MRI contrasts with a variety of sampling schemes, demonstrate that CDLMRI is capable of capturing structural dependencies between different contrasts. The learned priors indicate notable advantages in multi-contrast MR imaging and promising applications in quantitative MR imaging such as MR fingerprinting.

Low-rank magnetic resonance fingerprinting.

PURPOSE: Magnetic resonance fingerprinting (MRF) is a relatively new approach that provides quantitative MRI measures using randomized acquisition. Extraction of physical quantitative tissue parameters is performed offline, without the need of patient presence, based on acquisition with varying parameters and a dictionary generated according to the Bloch equation simulations. MRF uses hundreds of radio frequency (RF) excitation pulses for acquisition, and therefore, a high undersampling ratio in the sampling domain (k-space) is required for reasonable scanning time. This undersampling causes spatial artifacts that hamper the ability to accurately estimate the tissue's quantitative values. In this work, we introduce a new approach for quantitative MRI using MRF, called magnetic resonance fingerprinting with low rank (FLOR). METHODS: We exploit the low-rank property of the concatenated temporal imaging contrasts, on top of the fact that the MRF signal is sparsely represented in the generated dictionary domain. We present an iterative recovery scheme that consists of a gradient step followed by a low-rank projection using the singular value decomposition. RESULTS: Experimental results consist of retrospective sampling that allows comparison to a well defined reference, and prospective sampling that shows the performance of FLOR for a real-data sampling scenario. Both experiments demonstrate improved parameter accuracy compared to other compressed-sensing and low-rank based methods for MRF at 5% and 9% sampling ratios for the retrospective and prospective experiments, respectively. CONCLUSIONS: We have shown through retrospective and prospective experiments that by exploiting the low-rank nature of the MRF signal, FLOR recovers the MRF temporal undersampled images and provides more accurate parameter maps compared to previous iterative approaches.

PEAR: PEriodic And fixed Rank separation for fast fMRI.

PURPOSE: In functional MRI (fMRI), faster acquisition via undersampling of data can improve the spatial-temporal resolution trade-off and increase statistical robustness through increased degrees-of-freedom. High-quality reconstruction of fMRI data from undersampled measurements requires proper modeling of the data. We present an fMRI reconstruction approach based on modeling the fMRI signal as a sum of periodic and fixed rank components, for improved reconstruction from undersampled measurements. METHODS: The proposed approach decomposes the fMRI signal into a component which has a fixed rank and a component consisting of a sum of periodic signals which is sparse in the temporal Fourier domain. Data reconstruction is performed by solving a constrained problem that enforces a fixed, moderate rank on one of the components, and a limited number of temporal frequencies on the other. Our approach is coined PEAR - PEriodic And fixed Rank separation for fast fMRI. RESULTS: Experimental results include purely synthetic simulation, a simulation with real timecourses and retrospective undersampling of a real fMRI dataset. Evaluation was performed both quantitatively and visually versus ground truth, comparing PEAR to two additional recent methods for fMRI reconstruction from undersampled measurements. Results demonstrate PEAR's improvement in estimating the timecourses and activation maps versus the methods compared against at acceleration ratios of R = 8,10.66 (for simulated data) and R = 6.66,10 (for real data). CONCLUSIONS: This paper presents PEAR, an undersampled fMRI reconstruction approach based on decomposing the fMRI signal to periodic and fixed rank components. PEAR results in reconstruction with higher fidelity than when using a fixed-rank based model or a conventional Low-rank + Sparse algorithm. We have shown that splitting the functional information between the components leads to better modeling of fMRI, over state-of-the-art methods.

Reference-based MRI.

PURPOSE: In many clinical MRI scenarios, existing imaging information can be used to significantly shorten acquisition time or to improve Signal to Noise Ratio (SNR). In this paper the authors present a framework, referred to as FASTMER, for fast MRI by exploiting a reference image. METHODS: The proposed approach utilizes the possible similarity of the reference image to the acquired image, which exists in many clinical MRI imaging scenarios. Examples include similarity between adjacent slices in high resolution MRI, similarity between various contrasts in the same scan and similarity between different scans of the same patient. To account for the fact that the reference image may exhibit low similarity with the acquired image the authors develop an iterative weighted reconstruction approach, which tunes the weights according to the degree of similarity. RESULTS: Experimental results demonstrate the performance of the method in three different clinical MRI scenarios: The first example demonstrates SNR improvement in high resolution brain MRI, the second scenario exploits similarity between T2-weighted and fluid-attenuated inversion recovery (FLAIR) for fast FLAIR scanning and the last application utilizes similarity between baseline and follow-up scans for fast follow-up acquisition. The results show that FASTMER outperforms image reconstruction of existing state-of-the-art methods. CONCLUSIONS: The authors present a framework for fast MRI by exploiting a reference image. Recovery is based on an iterative algorithm that supports cases in which similarity to the reference scan is not guaranteed. This extends the applicability of the FASTMER to different MRI scanning scenarios. Thanks to the existence of reference images in various clinical imaging tasks, the proposed framework can play a major role in improving reconstruction in many MR applications.

Low rank magnetic resonance fingerprinting.

Magnetic Resonance Fingerprinting (MRF) is a relatively new approach that provides quantitative MRI using randomized acquisition. Extraction of physical quantitative tissue values is preformed off-line, based on acquisition with varying parameters and a dictionary generated according to the Bloch equations. MRF uses hundreds of radio frequency (RF) excitation pulses for acquisition, and therefore high under-sampling ratio in the sampling domain (k-space) is required. This under-sampling causes spatial artifacts that hamper the ability to accurately estimate the quantitative tissue values. In this work, we introduce a new approach for quantitative MRI using MRF, called Low Rank MRF. We exploit the low rank property of the temporal domain, on top of the well-known sparsity of the MRF signal in the generated dictionary domain. We present an iterative scheme that consists of a gradient step followed by a low rank projection using the singular value decomposition. Experiments on real MRI data demonstrate superior results compared to conventional implementation of compressed sensing for MRF at 15% sampling ratio.

Compressed sensing for longitudinal MRI: An adaptive-weighted approach.

PURPOSE: Repeated brain MRI scans are performed in many clinical scenarios, such as follow up of patients with tumors and therapy response assessment. In this paper, the authors show an approach to utilize former scans of the patient for the acceleration of repeated MRI scans. METHODS: The proposed approach utilizes the possible similarity of the repeated scans in longitudinal MRI studies. Since similarity is not guaranteed, sampling and reconstruction are adjusted during acquisition to match the actual similarity between the scans. The baseline MR scan is utilized both in the sampling stage, via adaptive sampling, and in the reconstruction stage, with weighted reconstruction. In adaptive sampling, k-space sampling locations are optimized during acquisition. Weighted reconstruction uses the locations of the nonzero coefficients in the sparse domains as a prior in the recovery process. The approach was tested on 2D and 3D MRI scans of patients with brain tumors. RESULTS: The longitudinal adaptive compressed sensing MRI (LACS-MRI) scheme provides reconstruction quality which outperforms other CS-based approaches for rapid MRI. Examples are shown on patients with brain tumors and demonstrate improved spatial resolution. Compared with data sampled at the Nyquist rate, LACS-MRI exhibits signal-to-error ratio (SER) of 24.8 dB with undersampling factor of 16.6 in 3D MRI. CONCLUSIONS: The authors presented an adaptive method for image reconstruction utilizing similarity of scans in longitudinal MRI studies, where possible. The proposed approach can significantly reduce scanning time in many applications that consist of disease follow-up and monitoring of longitudinal changes in brain MRI.

The effect of chemotherapy on optic pathway gliomas and their sub-components: A volumetric MR analysis study.

BACKGROUND: Optic pathway gliomas (OPG) represent 5% of pediatric brain tumors and compose a major therapeutic dilemma to the treating physicians. While chemotherapy is widely used for these tumors, our ability to predict radiological response is still lacking. In this study, we use volumetric imaging to examine in detail the long-term effect of chemotherapy on the tumor as well as its various sub-components. PROCEDURE: The tumors of 15 patients with OPG, treated with chemotherapy, were longitudinally measured using our novel, previously described volumetric method. Patients were treated with up to five lines of chemotherapy. Sufficient follow-up imaging data, and patient's numbers, allowed for analysis of two treatment lines. Volumetric measurements of the tumors were segmented into solid-non-enhancing, solid-enhancing, and cystic components. Outcome analysis was done per specific treatment line and for the overall follow-up period. RESULTS: An average reduction of 9.7% (±23%) in the gross-total-solid volume (GTSV) was noted following treatment with vincristine and carboplatin. The cystic component grew under therapy by an average of 12.6% (±39%). When measured over the course of the whole study period, the cystic component grew by an average of 35% (±100%) and the GTSV increased by 12% (±35%). CONCLUSION: Initial treatment with vincristine and carboplatin seems to have a minimal initial effect, mostly on the solid components. The cystic component in itself seems to be unaffected by chemotherapy, and contributes to the subsequent growth of the total volume. During the overall treatment period, both solid and cystic components grew regardless of combined treatment methods.

Semiautomatic segmentation and follow-up of multicomponent low-grade tumors in longitudinal brain MRI studies.

PURPOSE: Tracking the progression of low grade tumors (LGTs) is a challenging task, due to their slow growth rate and associated complex internal tumor components, such as heterogeneous enhancement, hemorrhage, and cysts. In this paper, the authors show a semiautomatic method to reliably track the volume of LGTs and the evolution of their internal components in longitudinal MRI scans. METHODS: The authors' method utilizes a spatiotemporal evolution modeling of the tumor and its internal components. Tumor components gray level parameters are estimated from the follow-up scan itself, obviating temporal normalization of gray levels. The tumor delineation procedure effectively incorporates internal classification of the baseline scan in the time-series as prior data to segment and classify a series of follow-up scans. The authors applied their method to 40 MRI scans of ten patients, acquired at two different institutions. Two types of LGTs were included: Optic pathway gliomas and thalamic astrocytomas. For each scan, a "gold standard" was obtained manually by experienced radiologists. The method is evaluated versus the gold standard with three measures: gross total volume error, total surface distance, and reliability of tracking tumor components evolution. RESULTS: Compared to the gold standard the authors' method exhibits a mean Dice similarity volumetric measure of 86.58% and a mean surface distance error of 0.25 mm. In terms of its reliability in tracking the evolution of the internal components, the method exhibits strong positive correlation with the gold standard. CONCLUSIONS: The authors' method provides accurate and repeatable delineation of the tumor and its internal components, which is essential for therapy assessment of LGTs. Reliable tracking of internal tumor components over time is novel and potentially will be useful to streamline and improve follow-up of brain tumors, with indolent growth and behavior.

PNist: interactive volumetric measurements of plexiform neurofibromas in MRI scans.

PURPOSE: Volumetric measurements of plexiform neurofibromas (PNs) are time consuming and error prone, as they require the delineation of the PN boundaries, which is mostly impractical in the daily clinical setup. Accurate volumetric measurements are seldom performed for these tumors mainly due to their great dispersion, size and multiple locations. This paper presents a semiautomatic method for segmentation of PN from STIR MRI scans. METHODS: Plexiform neurofibroma interactive segmentation tool (PNist) is a new tool to segment PNs in STIR MRI scans. The method is based on histogram tumor models computed from a training set. RESULTS: Experimental results from 28 datasets show an average absolute volume difference of 6.8 % with an average user time of approximately 7 min versus more than 13 min with manual delineation. In complex cases, the PNist user time is less than half in compared to state-of-the-art tools. CONCLUSIONS: PNist is a new method for the semiautomatic segmentation of PN lesions. Its simplicity and reliability make it unique among other state-of-the-art methods. It has the potential to become a clinical tool that allows the reliable evaluation of PN burden and progression.

Exploiting similarity in adjacent slices for compressed sensing MRI.

Due to fundamental characteristics of MRI that limit scan speedup, sub-sampling techniques such as compressed sensing (CS) have been developed for rapid MRI. Current CS MRI approaches utilize sparsity of the image in the wavelet or other transform domains to speed-up acquisition. Another drawback of MRI is its poor signal-to-noise ratio (SNR), which is proportional to the image slice thickness. In this paper, we use the difference between adjacent slices as the sparse domain for CS MRI. We propose to acquire thick MRI slices and to reconstruct the thin slices from the thick slices' data, utilizing the similarity between adjacent thin slices. The acquisition of thick slices, instead of thin ones, improves the total SNR of the reconstructed image. Experimental results show that the image reconstruction quality of the proposed method outperforms existing CS MRI methods using the same number of measurements.

Interactive segmentation of plexiform neurofibroma tissue: method and preliminary performance evaluation.

Plexiform neurofibromas (PNs) are a major manifestation of neurofibromatosis-1 (NF1), a common genetic disease involving the nervous system. Treatment decisions are mostly based on a gross assessment of changes in tumor using MRI. Accurate volumetric measurements are rarely performed in this kind of tumors mainly due to its great dispersion, size, and multiple locations. This paper presents a semi-automatic method for segmentation of PN from STIR MRI scans. The method starts with a user-based delineation of the tumor area in a single slice and automatically segments the PN lesions in the entire image based on the tumor connectivity. Experimental results on seven datasets, with lesion volumes in the range of 75-690 ml, yielded a mean absolute volume error of 10 % (after manual adjustment) as compared to manual segmentation by an expert radiologist. The mean computation and interaction time was 13 versus 63 min for manual annotation.

Prediction of brain MR scans in longitudinal tumor follow-up studies.

We present a new method for the estimation of the next brain MR scan in a longitudinal tumor follow-up study. Our method effectively incorporates information of the past scans in the time series to predict the future scan of the patient. Its advantages are that it requires no user intervention and does not assume any particular tumor growth model. Instead, the patient-specific tumor growth parameters are estimated individually from the past patient scans. To validate our method, we conducted an experimental study on four patients with Optic Path Gliomas (OPGs) and four patients with glioblastomas multiforma (GBM), each scanned at five time points. The tumor volumes in the predicted and actual future scans, both segmented by an expert radiologist, yield a mean volume overlap difference of 13.65% for the OPG patients, and 34.23% for the GBM patients.

MRI internal segmentation of optic pathway gliomas: clinical implementation of a novel algorithm.

PURPOSE: Optic pathway gliomas (OPGs) are diagnosed based on typical MR features and require careful monitoring with serial MRI. Reliable, serial radiological comparison of OPGs is a difficult task, where accuracy becomes very important for clinical decisions on treatment initiation and results. Current radiological methodology usually includes linear measurements that are limited in terms of precision and reproducibility. METHOD: We present a method that enables semiautomated segmentation and internal classification of OPGs using a novel algorithm. Our method begins with co-registration of the different sequences of an MR study so that T1 and T2 slices are realigned. The follow-up studies are then re-sliced according to the baseline study. The baseline tumor is segmented, with internal components classified into solid non-enhancing, solid-enhancing, and cystic components, and the volume is calculated. Tumor demarcation is then transferred onto the next study and the process repeated. Numerical values are correlated with clinical data such as treatment and visual ability. RESULTS: We have retrospectively implemented our method on 24 MR studies of three OPG patients. Clinical case reviews are presented here. The volumetric results have been correlated with clinical data and their implications are also discussed. CONCLUSIONS: The heterogeneity of OPGs, the long course, and the young age of the patients are all driving the demand for more efficient and accurate means of tumor follow-up. This method may allow better understanding of the natural history of the tumor and provide a more advanced means of treatment evaluation.

Automatic segmentation and components classification of optic pathway gliomas in MRI.

We present a new method for the automatic segmentation and components classification of brain Optic Pathway Gliomas (OPGs) from multi-spectral MRI datasets. Our method accurately identifies the sharp OPG boundaries and consistently delineates the missing contours by effectively incorporating prior location, shape, and intensity information. It then classifies the segmented OPG volume into its three main components--solid, enhancing, and cyst--with a probabilistic tumor tissue model generated from training datasets that accounts for the datasets grey-level differences. Experimental results on 25 datasets yield a mean OPG boundary surface distance error of 0.73mm and mean volume overlap difference of 30.6% as compared to manual segmentation by an expert radiologist. A follow-up patient study shows high correlation between the clinical tumor progression evaluation and the component classification results. To the best of our knowledge, ours is the first method for automatic OPG segmentation and component classification that may support quantitative disease progression and treatment efficacy evaluation.

Detection of one-lung intubation incidents.

Lung sounds are very common source for monitoring and diagnosis of pulmonary function. This approach can be used for detecting one lung intubation (OLI) during anesthesia or intensive care. In this paper, an algorithm for detecting OLI from lung sounds is presented. The algorithm assumes a multiple-input-multiple-output system, in which a multi-dimensional auto-regressive model relates the input (lungs) and the output (recorded sounds). An OLI detector is developed based on the generalized likelihood ratio test (GLRT), assuming coherent distributed sources for each lung. This method exhibited reliable results also when the lungs were modeled by incoherent distributed sources, which is a more accurate model for lung sources. The algorithm was tested using real breathing sounds recorded in an operating room, and it achieved an OLI detection rate of more than 95%, for each breathing cycle.