Abstract:Magnetic Resonance Imaging (MRI) is a widely used non-invasive clinical imaging modality. Unlike other medical imaging tools, such as X-rays or computed tomography (CT), the advantage of MRI is that it uses non-ionizing radiation. In addition, MRI can provide images with multiple contrast by using different pulse sequences and protocols. However, acquisition speed, which remains the main challenge for MRI, limits its clinical application. Clinicians have to compromise between spatial resolution, SNR, and scan time, which leads to sub-optimal performance. The acquisition speed of MRI can be improved by collecting fewer data samples. However, according to the Nyquist sampling theory, undersampling in k-space will lead to aliasing artifacts in the recovered image. The recent mathematical theory of compressed sensing has been developed to exploit the property of sparsity for signals/images. It states that if an image is sparse, it can be accurately reconstructed using a subset of the k-space data under certain conditions. Generally, the reconstruction is formulated as an optimization problem. The sparsity of the image is enforced by using a sparsifying transform. Total variation (TV) is one of the commonly used methods, which enforces the sparsity of the image gradients and provides good image quality. However, TV introduces patchy or painting-like artifacts in the reconstructed images. We introduce novel regularization penalties involving higher degree image derivatives to overcome the practical problems associated with the classical TV scheme. Motivated by novel reinterpretations of the classical TV regularizer, we derive two families of functionals, which we term as isotropic and anisotropic higher degree total variation (HDTV) penalties, respectively. The numerical comparisons of the proposed scheme with classical TV penalty, current second order methods, and wavelet algorithms demonstrate the performance improvement. Specically, the proposed algorithms minimize the staircase and ringing artifacts that are common with TV schemes and wavelet algorithms, while better preserving the singularities. Higher dimensional MRI is also challenging due to the above mentioned trade-offs. We propose a three-dimensional (3D) version of HDTV (3D-HDTV) to recover 3D datasets. One of the challenges associated with the HDTV framework is the high computational complexity of the algorithm. We introduce a novel computationally efficient algorithm for HDTV regularized image recovery problems. We find that this new algorithm improves the convergence rate by a factor of ten compared to the previously used method. We demonstrate the utility of 3D-HDTV regularization in the context of compressed sensing, denoising, and deblurring of 3D MR dataset and fluorescence microscope images. We show that 3D-HDTV outperforms 3D-TV schemes in terms of the signal to noise ratio (SNR) of the reconstructed images and its ability to preserve ridge-like details in the 3D datasets. To address speed limitations…

Directional TV algorithm for image reconstruction from sparse-view projections in EPR imaging

Directional TV algorithm for fast EPR imaging

Fast and Quantitative Reconstruction of 4D SPECT for Fan-beam Geometry Based on Novikov's Inverse Formula

An edge-preserving total nuclear variation minimization algorithm in EPR image reconstruction

Optimization-based Image Reconstruction from Sparsely Sampled Data in Electron Paramagnetic Resonance Imaging.

Low Dose PET Image Reconstruction with Total Variation Using Alternating Direction Method.

4D-image reconstruction directly from limited-angular-range data in continuous-wave electron paramagnetic resonance imaging

Directional-TV algorithm for image reconstruction from limited-angular-range data

Accurate Reconstruction of 4D Spectral-Spatial Images from Sparse-View Data in Continuous-Wave EPRI

A Fast Barzilai-Borwein Gradient Projection for Sparse Reconstruction Algorithm Based on 3D Modeling: Application to ERT Imaging

3D Tensor Based Nonlocal Low Rank Approximation in Dynamic PET Reconstruction.

Reconstruction for 3D PET Based on Total Variation Constrained Direct Fourier Method

3-D Weighted CBCT Iterative Reconstruction with TV Minimization from Angularly Under-Sampled Projection Data

Adaptive-weighted High Order TV Algorithm for Sparse-View CT Reconstruction

Preconditioned alternating projection algorithm for solving the penalized-likelihood SPECT reconstruction problem

Image Reconstruction From Sparse Data In Echo-Planar Imaging

Prototyping optimization-based image reconstructions from limited-angular-range data in dual-energy CT

Accelerated exponential parameterization of T2 relaxation with model-driven low rank and sparsity priors (MORASA).

Pet Image Reconstruction From Under-Sampled Data

Novel Compressed Sensing Algorithms with Applications to Magnetic Resonance Imaging

Spherical Polar Fourier EAP and Odf Reconstruction Via Compressed Sensing in Diffusion Mri