Journal article icon

Journal article

Approximate message passing with a colored aliasing model for variable density Fourier sampled images

Abstract:
The Approximate Message Passing (AMP) algorithm eciently reconstructs signals which have been sampled with large i.i.d. sub-Gaussian sensing matrices. However, when Fourier coecients of a signal with non-uniform spectral density are sampled, such as in Magnetic Resonance Imaging (MRI), the aliasing is intrinsically colored. Consequently, AMP’s i.i.d. state evolution is no longer accurate and the algorithm encounters convergence problems. In response, we propose an algorithm based on Orthogonal Approximate Message Passing (OAMP) that uses the wavelet domain to model the colored aliasing. We present empirical evidence that a structured state evolution occurs, where the e↵ective noise covariance matrix is diagonal with one unique entry per subband. A benefit of state evolution is that Stein’s Unbiased Risk Estimate (SURE) can be e↵ectively implemented, yielding an algorithm with no free parameters. We empirically evaluate the e↵ectiveness of the parameterfree algorithm on a synthetic image with three variable density sampling schemes and find that it converges in over 20x fewer iterations than optimally tuned Fast Iterative Shrinkage-Thresholding (FISTA).
Publication status:
Published
Peer review status:
Peer reviewed

Actions


Access Document


Files:
Publisher copy:
10.1109/OJSP.2020.3025228

Authors


More by this author
Institution:
University of Oxford
Division:
MPLS
Department:
Mathematical Institute
Sub department:
Mathematical Institute
Role:
Author


Publisher:
Institute of Electrical and Electronics Engineers
Journal:
IEEE Open Journal of Signal Processing More from this journal
Volume:
1
Pages:
146-158
Publication date:
2020-09-18
Acceptance date:
2020-09-06
DOI:
EISSN:
2644-1322


Language:
English
Keywords:
Pubs id:
1130741
Local pid:
pubs:1130741
Deposit date:
2020-09-06

Terms of use



Views and Downloads






If you are the owner of this record, you can report an update to it here: Report update to this record

TO TOP