Curl-constrained Gradient Estimation for Image Recovery from Highly Incomplete Spectral Data

C. Ravazzi, G. Coluccia, E. Magli

IEEE Transactions on Image Processing, vol. 26, no. 6, pp. 2656-2668, June 2017

Abstract

In this paper, we introduce new gradient-based methods for image recovery from a small collection of spectral coefficients of the Fourier transform, which is of particular interest for several scanning technologies, such as magnetic resonance imaging.

Since gradients of a medical image are much more sparse or compressible than the corresponding image, classical ℓ_{1}-minimization methods have been used to recover these relative differences. The image values can then be obtained by integration algorithms imposing boundary constraints.

Compared to classical gradient recovery methods, we propose two new techniques that improve reconstruction. First, we cast the gradient recovery problem as a compressed sensing problem taking into account that the curl of the gradient field should be zero. Second, inspired by the emerging field of signal processing on graphs, we formulate the gradient recovery problem as an inverse problem on graphs. Iteratively reweighted ℓ_{1} recovery methods are proposed to recover these relative differences and the structure of the similarity graph.

Once the gradient field is estimated, the image is recovered from the compressed Fourier measurements using least squares estimation.

Numerical experiments shows that the proposed approach outperforms the state-of-the-art image recovery methods.

Additional material

Click on an item to open a preview, then on (top-left) to download it.