115 research outputs found
Recommended from our members
Nonconvex compressive sensing reconstruction for tensor using structures in modes
This paper focuses on the reconstruction of a tensor captured using Compressive Sensing (CS). Instead of processing the signals via vectorization as is done in conventional CS, in tensor CS high dimensional signals are kept in their original formats, which benefits hardware implementation and eases memory requirements. In addition, more structures exist in a tensor along its various dimensions than in its vectorized format. Utilizing these various structures, this paper proposes a general reconstruction approach for tensor CS. Employing the proximity operator of a nonconvex norm function, a special case for a tensor with low rank and sparse structures is elaborated, which is shown to outperform the state-of-art tensor CS reconstruction methods when applied to magnetic resonance imaging and hyper-spectral imaging.This work is supported by EPSRC Research Grant EP/K033700/1; the Natural Science Foundation of China (61401018, U1334202); the Fundamental Research Funds for the Central Universities (No. 2014JBM149); the State Key Laboratory of Rail Traffic Control and Safety (RCS2016ZT014) of Beijing Jiaotong University; the Key Grant Project of Chinese Ministry of Education (313006).This is the author accepted manuscript. It is currently embargoed pending publication
Low-rank and sparse reconstruction in dynamic magnetic resonance imaging via proximal splitting methods
Dynamic magnetic resonance imaging (MRI) consists of collecting multiple MR images in time, resulting in a spatio-temporal signal. However, MRI intrinsically suffers from long acquisition times due to various constraints. This limits the full potential of dynamic MR imaging, such as obtaining high spatial and temporal resolutions which are crucial to observe dynamic phenomena. This dissertation addresses the problem of the reconstruction of dynamic MR images from a limited amount of samples arising from a nuclear magnetic resonance experiment. The term limited can be explained by the approach taken in this thesis to speed up scan time, which is based on violating the Nyquist criterion by skipping measurements that would be normally acquired in a standard MRI procedure. The resulting problem can be classified in the general framework of linear ill-posed inverse problems. This thesis shows how low-dimensional signal models, specifically lowrank and sparsity, can help in the reconstruction of dynamic images from partial measurements. The use of these models are justified by significant developments in signal recovery techniques from partial data that have emerged in recent years in signal processing. The major contributions of this thesis are the development and characterisation of fast and efficient computational tools using convex low-rank and sparse constraints via proximal gradient methods, the development and characterisation of a novel joint reconstruction–separation method via the low-rank plus sparse matrix decomposition technique, and the development and characterisation of low-rank based recovery methods in the context of dynamic parallel MRI. Finally, an additional contribution of this thesis is to formulate the various MR image reconstruction problems in the context of convex optimisation to develop algorithms based on proximal splitting methods
Joint Sensing Matrix and Sparsifying Dictionary Optimization for Tensor Compressive Sensing.
Tensor compressive sensing (TCS) is a multidimensional framework of compressive sensing (CS), and it is
advantageous in terms of reducing the amount of storage, easing
hardware implementations, and preserving multidimensional
structures of signals in comparison to a conventional CS system.
In a TCS system, instead of using a random sensing matrix and
a predefined dictionary, the average-case performance can be
further improved by employing an optimized multidimensional
sensing matrix and a learned multilinear sparsifying dictionary.
In this paper, we propose an approach that jointly optimizes
the sensing matrix and dictionary for a TCS system. For the
sensing matrix design in TCS, an extended separable approach
with a closed form solution and a novel iterative nonseparable
method are proposed when the multilinear dictionary is fixed.
In addition, a multidimensional dictionary learning method that
takes advantages of the multidimensional structure is derived,
and the influence of sensing matrices is taken into account in the
learning process. A joint optimization is achieved via alternately
iterating the optimization of the sensing matrix and dictionary.
Numerical experiments using both synthetic data and real images
demonstrate the superiority of the proposed approache
Hyperspectral Image Restoration via Total Variation Regularized Low-rank Tensor Decomposition
Hyperspectral images (HSIs) are often corrupted by a mixture of several types
of noise during the acquisition process, e.g., Gaussian noise, impulse noise,
dead lines, stripes, and many others. Such complex noise could degrade the
quality of the acquired HSIs, limiting the precision of the subsequent
processing. In this paper, we present a novel tensor-based HSI restoration
approach by fully identifying the intrinsic structures of the clean HSI part
and the mixed noise part respectively. Specifically, for the clean HSI part, we
use tensor Tucker decomposition to describe the global correlation among all
bands, and an anisotropic spatial-spectral total variation (SSTV)
regularization to characterize the piecewise smooth structure in both spatial
and spectral domains. For the mixed noise part, we adopt the norm
regularization to detect the sparse noise, including stripes, impulse noise,
and dead pixels. Despite that TV regulariztion has the ability of removing
Gaussian noise, the Frobenius norm term is further used to model heavy Gaussian
noise for some real-world scenarios. Then, we develop an efficient algorithm
for solving the resulting optimization problem by using the augmented Lagrange
multiplier (ALM) method. Finally, extensive experiments on simulated and
real-world noise HSIs are carried out to demonstrate the superiority of the
proposed method over the existing state-of-the-art ones.Comment: 15 pages, 20 figure
Recommended from our members
Structured Tensor Recovery and Decomposition
Tensors, a.k.a. multi-dimensional arrays, arise naturally when modeling higher-order objects and relations. Among ubiquitous applications including image processing, collaborative filtering, demand forecasting and higher-order statistics, there are two recurring themes in general: tensor recovery and tensor decomposition. The first one aims to recover the underlying tensor from incomplete information; the second one is to study a variety of tensor decompositions to represent the array more concisely and moreover to capture the salient characteristics of the underlying data. Both topics are respectively addressed in this thesis.
Chapter 2 and Chapter 3 focus on low-rank tensor recovery (LRTR) from both theoretical and algorithmic perspectives. In Chapter 2, we first provide a negative result to the sum of nuclear norms (SNN) model---an existing convex model widely used for LRTR; then we propose a novel convex model and prove this new model is better than the SNN model in terms of the number of measurements required to recover the underlying low-rank tensor. In Chapter 3, we first build up the connection between robust low-rank tensor recovery and the compressive principle component pursuit (CPCP), a convex model for robust low-rank matrix recovery. Then we focus on developing convergent and scalable optimization methods to solve the CPCP problem. In specific, our convergent method, proposed by combining classical ideas from Frank-Wolfe and proximal methods, achieves scalability with linear per-iteration cost.
Chapter 4 generalizes the successive rank-one approximation (SROA) scheme for matrix eigen-decomposition to a special class of tensors called symmetric and orthogonally decomposable (SOD) tensor. We prove that the SROA scheme can robustly recover the symmetric canonical decomposition of the underlying SOD tensor even in the presence of noise. Perturbation bounds, which can be regarded as a higher-order generalization of the Davis-Kahan theorem, are provided in terms of the noise magnitude
- …