Toward 0-norm Toward 0-norm Reconstruction, and A Reconstruction, and A Nullspace Technique Nullspace Technique for Compressive for Compressive Sampling Sampling Christine Law Christine Law Gary Glover Gary Glover Dept. of EE, Dept. of Dept. of EE, Dept. of Radiology Radiology Stanford University Stanford University
49
Embed
Toward 0-norm Reconstruction, and A Nullspace Technique for Compressive Sampling Christine Law Gary Glover Dept. of EE, Dept. of Radiology Stanford University.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Toward 0-norm Toward 0-norm Reconstruction, and A Reconstruction, and A Nullspace Technique Nullspace Technique
for Compressive for Compressive SamplingSamplingChristine LawChristine LawGary GloverGary Glover
Dept. of EE, Dept. of Dept. of EE, Dept. of RadiologyRadiology
Stanford UniversityStanford University
OutlineOutline
0-norm Magnetic Resonance 0-norm Magnetic Resonance Imaging (MRI) reconstruction.Imaging (MRI) reconstruction. HomotopicHomotopic Convex iterationConvex iteration
Signal separation exampleSignal separation example
1-norm deconvolution in fMRI. 1-norm deconvolution in fMRI. Improvement in cardinality Improvement in cardinality
constraint problem with nullspace constraint problem with nullspace technique.technique.
Suppose Suppose xx in in RRnn is is kk-sparse and we are given -sparse and we are given mm Fourier coefficients with frequencies selected Fourier coefficients with frequencies selected uniformly at random. If uniformly at random. If
mm ≥ ≥ kk log log22 (1 + (1 + n n / / kk))
then then
reconstructs reconstructs xx exactly with overwhelming exactly with overwhelming probability. probability.
1arg min s.t. =x x x y
n
k = 2
Rat dies 1 week after drinking poisoned wine
Example by Anna Gilbert
x1 x2 x3 x4 x5 x6 x7
y1
y2
y3
10011 01
0101011
0010111
x1 x2 x3 x4 x5 x6 x7
y1
y2
y3
1001 01
0101011
00
1
110 11
n=7m=3
k=1
Reconstruction by Reconstruction by OptimizationOptimization
Compressed Sensing theory (2004 Donoho, Candes): under certain conditions,
Candes et al. IEEE Trans. Information Theory 2006 52(2):489 Donoho. IEEE Trans. Information Theory 2006 52(4):1289
1arg min s.t. =x x x y
0arg min s.t. =x x x y
y are measurements (rats)
x are sensors (wine)
0-norm reconstruction0-norm reconstruction
Try to solve 0-norm directly.Try to solve 0-norm directly.
For For pp-norm, where 0 < -norm, where 0 < pp < 1 < 1
Chartrand (2006) demonstrated fewer samples Chartrand (2006) demonstrated fewer samples y y required than 1-norm formulation.required than 1-norm formulation.
Chartrand. IEEE Signal Processing Letters. 2007: 14(10) 707-710.Chartrand. IEEE Signal Processing Letters. 2007: 14(10) 707-710.
0arg min s.t. =u u u y
Trzasko (2007): Rewrite the Trzasko (2007): Rewrite the problem problem
0, .min lim s.t
i iu u u y
00, lim
i iu u
Trzasko Trzasko et alet al. IEEE SP 14th workshop on statistical signal processing. 2007. 176-180.. IEEE SP 14th workshop on statistical signal processing. 2007. 176-180.
0arg min s.t. =u u u y
where where is tanh, laplace, log, is tanh, laplace, log, etcetc. . such thatsuch that
Homotopic function in Homotopic function in 1D1D
Laplace function: 1 xe
Start as 1-norm problem, then reduce slowly and approach 0-norm function.
1000
1
Homotopic methodHomotopic method
2
2min ,
2iiE
uu u y
T H,
,
Find minimum of by zeroing Lagrangian gradient:
= 0
isdiag
u
E
E u u u y
d d uu
u
T
11 T
1
,
,
It's a contraction
Iterate Conjugate Gradient method until .
.
t t
t t
u u y
u u y
u u
DemonstrationDemonstration when is big (1when is big (1stst iteration), solving 1-norm problem. iteration), solving 1-norm problem. reduce reduce to approach 0-norm solution. to approach 0-norm solution.
061%x
02%x
original x∆
: image to solve
: k-space measurements (undersampled)
: Total Variation (TV) operator
: Fourier matrix (masked)
u
y
Example 1Example 1
F
F
original phantom - reconstructionerror 20 log
original phantom
original subsampled
Zero-filled Reconstruction Fourier sample mask
HomotopicHomotopic result: use 4% Fourier data error: -66.2 dB 85 seconds
1-norm result: use 4% Fourier data error: -11.4 dB 542 seconds
In cardinality constraint problem, In cardinality constraint problem, convex iteration and nullspace convex iteration and nullspace technique success more often than 1-technique success more often than 1-norm.norm.