Discrete Shearlets as a Sparsifying Transform in Low-Rank Plus Sparse Decomposition for Undersampled (k, t)-Space MR
Nicholas E Protonotarios1, Evangelia Tzampazidou2,3, George A Kastis2
1Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 0WA, UK.
Journal of Imaging
|February 24, 2022
Summary
Discrete shearlets effectively decompose magnetic resonance imaging (MRI) data, separating motion from enhancement. This method accurately captures dynamic contrast-enhanced (DCE) and small bowel data, even with undersampling.
Area of Science:
- Medical Imaging
- Signal Processing
- Biomedical Engineering
Background:
- Magnetic resonance imaging (MRI) requires robust methods for analyzing complex dynamic data.
- Traditional sparsifying transforms face challenges in accurately representing discontinuities and edges in MRI.
- Low-rank plus sparse (L+S) decomposition is a promising technique for separating signal components.
Purpose of the Study:
- To evaluate discrete shearlets as a sparsifying transform within an L+S decomposition framework for undersampled MRI data.
- To compare the performance of discrete shearlets against other sparsifying transforms, specifically in dynamic contrast-enhanced (DCE) and small bowel imaging.
- To assess the accuracy of motion estimation and the isolation of physiological signals like bowel motility.
Main Methods:
- Implementation of an L+S decomposition algorithm utilizing discrete shearlets.
- Evaluation on simulated dynamic contrast-enhanced (DCE) and small bowel MRI datasets.
- Comparison of reconstruction performance against the k-t FOCUSS algorithm.
- Analysis of motion estimation accuracy and motility metrics derived from the sparse component.
Main Results:
- Discrete shearlets successfully separated low-rank (background, periodic motion) from sparse components (enhancement, bowel motility) in both DCE and small bowel data.
- Motion estimated from the low-rank component of DCE data showed higher fidelity to ground truth deformations compared to other methods.
- Bowel motility metrics derived from the sparse component of free-breathing data were comparable to breath-holding data, even with significant undersampling.
- The discrete shearlet-based L+S decomposition demonstrated robust performance in isolating rapid/random bowel motility.
Conclusions:
- Discrete shearlets are a highly effective sparsifying transform for L+S decomposition in undersampled MRI.
- The proposed method accurately separates physiological motion and enhancement, improving quantitative analysis.
- This approach holds significant potential for improving the quality and interpretability of dynamic MRI studies, particularly for gastrointestinal applications.
Related Concept Videos
Singularity Functions for Shear
227
In structural analysis, singularity functions are crucial in simplifying the representation of shear forces in beams under discontinuous loading. These functions describe discontinuous variations in shear force across a beam with varying loads by using a single mathematical expression, regardless of the complexity of the loading conditions. The singularity functions are derived from creating a free-body diagram of the beam and then making conceptual cuts at specific points to examine the...
227
Extraction: Partition and Distribution Coefficients
3.5K
The distribution law or Nernst's distribution law is the law that governs the distribution of a solute between two immiscible solvents. This law, also known as the partition law, states that if a solute is added to the mixture of two immiscible solvents at a constant temperature, the solute is distributed between the two solvents in such a way that the ratio of solute concentrations in the solvents remains constant at equilibrium.
For extracting a solute from an aqueous phase into an...
For extracting a solute from an aqueous phase into an...
3.5K
Upsampling
346
Managing signal sampling rates is essential in digital signal processing to maintain signal integrity. A decimated signal, characterized by a reduced frequency range due to its lower sampling rate, can be upsampled by inserting zeros between each sample. This upsampling process expands the original spectrum and introduces repeated spectral replicas at intervals dictated by the new Nyquist frequency. To refine this zero-inserted sequence, it is passed through a lowpass filter with a cutoff...
346
Discrete Fourier Transform
453
The Discrete Fourier Transform (DFT) is a fundamental tool in signal processing, extending the discrete-time Fourier transform by evaluating discrete signals at uniformly spaced frequency intervals. This transformation converts a finite sequence of time-domain samples into frequency components, each representing complex sinusoids ordered by frequency. The DFT translates these sequences into the frequency domain, effectively indicating the magnitude and phase of each frequency component present...
453
Downsampling
291
When considering a sampled sequence with zero values between sampling instants, one can replace it by taking every N-th value of the sequence. At these integer multiples of N, the original and sampled sequences coincide. This process, known as decimation, involves extracting every N-th sample from a sequence, thereby creating a more efficient sequence.
The Fourier transform of the decimated sequence reveals a combination of scaled and shifted versions of the original spectrum. This...
The Fourier transform of the decimated sequence reveals a combination of scaled and shifted versions of the original spectrum. This...
291
Three-Dimensional Analysis of Strain
318
Three-dimensional strain analysis is crucial for understanding how materials deform under stress, particularly in elastic, homogeneous materials. This method employs principal stress axes to simplify complex stress states into more understandable forms. Subjected to stress, a small cubic element within a material either expands or contracts along these axes, transforming into a rectangular parallelepiped. This transformation effectively illustrates the material's deformation. The principal...
318


