# A General Total Variation Minimization Theorem for Compressed Sensing Based Interior Tomography

Authors: Weimin Han; Hengyong Yu; Ge Wang
Journal: International Journal of Biomedical Imaging
 (2009)

---

## Abstract

Recently, in the compressed sensing framework we found that a two-dimensional interior region-of-interest (ROI) can be exactly reconstructed via the total variation minimization if the ROI is piecewise constant (Yu and Wang, 2009). Here we present a general theorem charactering a minimization property for a piecewise constant function defined on a domain in any dimension. Our major mathematical tool to prove this result is functional analysis without involving the Dirac delta function, which was heuristically used by Yu and Wang (2009).


---

## Body

## 1. Introduction

While in general an interior region-of-interest (ROI) cannot be uniquely reconstructed from projection data only associated with lines through the ROI [1, 2],  in the compressed sensing framework,  we recently found that a two-dimensional interior ROI can be exactly reconstructed via the total variation minimization if the function on the ROI is piecewise constant [3, 4].  The major idea behind our analysis is that the total variations of a piecewise constant function and a smooth artifact function are separable. The main mathematical tool is the expression of the two-dimensional gradient in terms of the Dirac delta function. In our analysis [3], the Delta function was instrumental but applied heuristically without mathematical rigor.  In this note, we will prove rigorously a more general theorem, as an extension of the total variation minimization property presented in [3], to characterize the total variation minimization property for a piecewise constant function defined on a domain in any dimension. Such a theorem may serve as a theoretical basis for further development of interior tomography algorithms.


## 2. Theoretical Result

For piecewise constant or piecewise smooth functions, it is natural to use the space of functions of bounded variation [5] to capture the discontinuities. For an integer d>0, let Ω⊂ℝd be a d-dimensional open bounded set and denote its boundary by ∂Ω. Then the space of functions of bounded variation is
(1)BV(Ω)={υ∈L1(Ω)∣∥υ∥BV(Ω)<∞}.
It is a Banach space with the norm
(2)∥υ∥BV(Ω)=∥υ∥L1(Ω)+∫Ω|Dυ|,
where ∥υ∥L1(Ω) is the integral of |υ(x)| over Ω, and
(3)∫Ω|Dυ|=sup{∫Ωυdivϕdx:ϕ∈C01(Ω)d,|ϕ|≤1inΩ}
is the total variation of the function υ∈BV(Ω). Here C01(Ω) is the space of continuously differentiable functions that vanish on ∂Ω, and for ϕ=(ϕ1,…,ϕd)T∈C01(Ω)d, divϕ=∑i=1d∂ϕi/∂xi. The Sobolev space
(4)W1,1(Ω)={υ∈L1(Ω):|∇υ|∈L1(Ω)}
is a subspace ofBV(Ω) and
(5)∫Ω|Dυ|=∫Ω|∇υ|dx,∀υ∈W1,1(Ω).We that assumeΩ has a piecewise C1 boundary, and it is decomposed into a union of a finite number of subsets with disjoint interiors
(6)Ω̅=⋃m=1MΩm¯
such that each subset Ωm has a piecewise C1 boundary. The unit outward normal vector on ∂Ωmis denoted by υm. Denote Γij=Ωi¯⋂Ωj¯, which may be empty for some pairs of i and j between 1 and M. We write meas(Γij) for the (d-1)-dimensional measure of Γij; it is the area of Γij for d=3, and the length of Γijfor d=2. The symbol ∑i<j will refer to a summation for those  i and j with a nonempty Γij in the range 1≤i<j≤M. The main result of this note is the following.Theorem 1.
Letf be a piecewise constant function corresponding to the decomposition (6): f(x)=cm∈ℝ for x∈Ωm, 1≤m≤M. Then we have
(7)∫Ω|D(f+g)|=∫Ω|Df|+∫Ω|∇g|dx,∀g∈W1,1(Ω),(8)∫Ω|Df|=∑i<j|ci-cj|meas(Γij).
Consequently, we have the minimization property
(9)∫Ω|Df|≤∫Ω|D(f+g)|,∀g∈W1,1(Ω).Proof.
After an integration by parts and some rearrangement, we have, for anyϕ∈C01(Ω)d,
(10)∫Ω(f+g)divϕdx=∑m=1Mcm∫∂Ωmϕ(x)·υm(x)ds-∫Ω∇g(x)·ϕ(x)dx=∑i<j(ci-cj)∫Γijϕ(x)·υi(x)ds-∫Ω∇g(x)·ϕ(x)dx.
By the definition (3),
(11)∫Ω|D(f+g)|=sup{∑i<j(ci-cj)∫Γijϕ(x)·υi(x)ds-∫Ω∇g(x)·ϕ(x)dx:ϕ∈C01(Ω)d,|ϕ|≤1inΩ{∑i<j(ci-cj)∫Γijϕ(x)}.
Taking g(x)=0 in (11), the formula (8) follows (cf. the argument in the next paragraph for a more general situation). Moreover, from (11) again,
(12)∫Ω|D(f+g)|≤∫Ω|Df|+∫Ω|∇g|dx.
For the opposite inequality, we first consider the case whereg∈C1(Ω̅). For any ε>0, define two open subsets
(13)Ωε∂={x∈Ω:dist(x,∂Ω)<ε},ΩεΓ={x∈Ω∖Ωε∂:mini<jdist(x,Γij)<ε}.
Here, dist(x,D)=min{|x-y|:y∈D} is the distance between x and a closed set D. Obviously, for some constant c>0,
(14)meas(Ωε∂)+meas(ΩεΓ)<cε.
We start with a function  ψε(x)∈C(Ω̅)d satisfying
(15)|ψε(x)|≤1forx∈Ω,|ψε(x)|=0forx∈Ωε/4∂,ψε(x)=sgn(ci-cj)υi(x)|υi(x)|forx∈Γij⋂(Ω∖Ωε/2∂),ψε(x)·∇g(x)=-|∇g(x)|forx∈Ω∖(Ωε/2∂⋃Ωε/2Γ)
and then apply the well-known mollification technique in the theory of Sobolev space [6] to define
(16)ϕε,δ(x)=∫Bδηδ(x-y)ψε(y)dy,
where Bδ is the ball of radius δ centered at the origin, ηδ(x)=η(x/δ)/δd, and
(17)η(x)={c0e1/(|x|2-1)if|x|<1,0if|x|≥1,c0-1=∫|x|<1e1/(|x|2-1)dx.
Then |ϕε,δ(x)|≤1 for x∈Ω, and for δ sufficiently small, ϕε,δ(x)∈C0∞(Ω)d. Moreover, as δ→0, ϕε,δ(x) converges uniformly to ψε(x) for x∈Ω∖(Ωε/2∂⋃Ωε/2Γ). Thus from (11),
(18)∫Ω|D(f+g)|≥∑i<j(ci-cj)∫Γijϕε,δ(x)·υi(x)ds-∫Ω∇g(x)·ϕε,δ(x)dx
and as  δ→0, we obtain, with some constant c1>0(19)∫Ω|D(f+g)|≥∑i<j(ci-cj)∫Γijψε(x)·υi(x)ds-∫Ω∇g(x)·ψε(x)dx-c1ε.
Using the defining properties of ψε, we further have
(20)∫Ω|D(f+g)|≥∫Ω|Df|+∫Ω|∇g|dx-c2ε
for some other constant  c2>0.  Since ε>0 is arbitrary, we obtain from the above relation that
(21)∫Ω|D(f+g)|≥∫Ω|Df|+∫Ω|∇g|dx.
Combining (12) and (21), we conclude (7) for g∈C1(Ω̅).
Forg∈W1,1(Ω), we use the density ofC1(Ω̅) in W1,1(Ω) [6] and choose {gn}⊂C1(Ω̅) such that
(22)gn→ginW1,1(Ω),asn→∞.
Since (3) defines a seminorm on BV(Ω), we have
(23)|∫Ω|D(f+g)|-∫Ω|D(f+gn)||≤∫Ω|∇(g-gn)|dx.
Thus, taking this limit n→∞ in (7) for gn⊂C1(Ω̅), we obtain (7) for g∈W1,1(Ω).As an example of (8), let Ω=Br0⊂ℝ2 be a disk of radius r0 centered at the origin. Consider a piecewise constant, radial function f(r) defined on Br0 such that it has a jump jm∈ℝ at rm, 1≤m≤M, where 0<r1<⋯<rM<r0. Then by (8), we have
(24)∫Br0|Df|=2π∑m=1M|jm|rm
(cf.  [3, Theorem 2.2]).


## 3. Discussions and Conclusion

Some comments on the name of our approach “CS-based interior tomography” are in order.  In the strict sense, compressed sensing refers to situations where the sampling scheme is built (often with random techniques) to achieve specific properties for satisfactory recovery of an underlying signal, rather than imposed by a specific detector arrangement as in limited data tomography. However, in a broad sense, compressed sensing can be interpreted as achieving better reconstruction from less data relative to the common practice.  Hence, while the current name is not far off, an alternative phrase for our approach can be “total variation minimization-based interior tomography.”In conclusion, we have extended the total variation minimization property of a piecewise function from two-dimensions to any dimensionality in the Sobolev space, which can be used for exact reconstruction of any piecewise function on an ROI by minimizing its total variation under the constraint of the truncated projection data through the ROI.  Previously, we implemented an alternating iterative reconstruction algorithm to minimize the total variation, which is time-consuming and needs improvement. Under the guidance of the theoretical finding presented here, we are working to develop a multidimensional ROI reconstruction algorithm for better performance.  Clearly, major efforts are still needed in this direction.



---

*Source: 125871-2009-11-17.xml*