EFFECTIVE COMPRESSIVE SENSING VIA REWEIGHTED TOTAL VARIATION
AND WEIGHTED NUCLEAR NORM REGULARIZATION
Mingli Zhang
1
, Christian Desrosiers
1
, Caiming Zhang
2,3
1
Software and IT Engineering department,
´
Ecole de technologie sup
´
erieure, Montreal, Canada, H3C 1K3
2
Shandong Provincial Key Laboratory of Digital Media Technology, Jinan 250014, China
3
School of Computer Science and Technology, Shandong University, Jinan 250100, China
ABSTRACT
Total variation (TV) and non-local patch similarity have been
used successfully to enhance the performance of compressive
sensing (CS) approaches. However, such techniques can of-
ten remove important details in the image or introduce recon-
struction artifacts. This paper presents a novel CS method,
which uses an adaptive reweighted TV strategy to better pre-
serve image edges. Our method also leverages the redundancy
of non-local image patches through the use of weighted low
rank regularization. An optimization strategy based on the
ADMM algorithm is used to reconstruct images efficiently.
Experimental results show our method to outperform state-
of-the-art CS approaches, for various sampling ratios.
Index Terms— Compressive sensing, Total variation,
Reweighted TV, Non-local self similarity, Low rank, ADMM.
1. INTRODUCTION
Compressive (or compressed) sensing (CS) has been widely
exploited in image processing, with numerous applications in
photography [1], video [2], spectral imaging [3] and medical
imaging [4, 5]. The key idea of this technique is that, if the
sampling matrix satisfies a condition known as the restricted
isometry property (RIP), a sparse signal can be accurately
reconstructed from undersampled measurements [6–8]. For-
mally, the model for reconstructing an image x ∈ R
n
from
measurements y ∈ R
m
can be formulated as
y = Φx + ν, (1)
where Φ ∈ R
m×n
is a known measurement operator, and ν is
additive noise (e.g., Gaussian, Rice, etc.).
In most CS methods, the task of recovering x is defined
as an inverse problem:
ˆx = argmin
x
1
2
kΦx − yk
2
2
+ λΨ(x), (2)
This work is supported by
´
ETS, FRQNT (210057) and National Natural
Science Foundation of China (61272431,61373080).
where Ψ(x) is a regularization prior, and λ is a parameter con-
trolling the trade-off between data fidelity and the regulariza-
tion of x. Generally, the regularization prior is based on the
principle that x is sparse under some suitable transform such
as wavelets [9]. Another popular regularization approach is
total variation (TV) [9,10], which is represented as:
TV(∇x) =
n
X
i=1
p
(∇
1
x
i
)
2
+ (∇
2
x
i
)
2
. (3)
Because this approach penalizes gradients uniformly, it may
result in the loss of image details like edges. To overcome this
problem, a reweighted TV model was proposed in [7], where
the gradient of a pixel i is penalized according to a weight w
i
,
i.e.
TV(∇x, w) =
n
X
i=1
w
i
p
(∇
1
x
i
)
2
+ (∇
2
x
i
)
2
. (4)
Weights w are updated iteratively from x in such a way that
regions with sharp gradients (e.g., edges) are less penalized
than uniform ones:
w
t+1
i
=
1
k∇x
t
i
k
2
+
, (5)
where is a small positive constant.
Another strategy to improve reconstruction is to exploit
the redundancy of local patterns in the image, represented
as small patches of pixels [11–14]. In this strategy, groups
of similar patches are regularized by applying a sparsifying
transform, for instance based on wavelets [13] or dictionary
learning [14]. Similar patches can also be regularized us-
ing the fact that matrices having these patches as columns or
rows have a low-rank [12, 15]. While CS methods based on
patch similarity can significantly improve the reconstruction
when few samples are available, they also have the tendency
to over-smooth images by “averaging” similar patches [16].
Motivated by the aforementioned observations, we present
a novel CS method based on reweighted TV and weighted
low-rank regularization of similar patches. The main contri-
butions of our work are as follows: