|Home | About | Journals | Submit | Contact Us | Français|
Recently, in the compressed sensing framework we proved that an interior ROI can be exactly reconstructed via the total variation minimization if the ROI is piecewise constant. In the proofs, we implicitly utilized the property that if an artifact image assumes a constant value within the ROI then this constant must be zero. Here we prove this property in the space of square integrable functions.
It is well known that an interior region-of-interest (ROI) cannot be uniquely reconstructed from projection data only associated with lines through the ROI. Interestingly, our group and others independently established that this interior problem can be exactly and stably solved if a sub-region in the ROI is precisely known (Ye et al., 2007; Yu et al., 2008; Courdurier et al., 2008; Kudo et al., 2008). Recently, in the compressed sensing framework we proved that an interior ROI can be exactly reconstructed via the total variation minimization if the imaging object is piecewise constant (Yu and Wang, 2009; Han et al., 2009). In our proofs of Theorems 2.2 and 2.3, after we made the conclusion that the minimum was achieved at λ = 0 in (Yu and Wang, 2009) we implicitly utilized the property that if an artifact image assumes a constant value within the ROI then this constant must be zero. Here we provide a proof for this property so that our compressed sensing based interior tomography theory is expressed rigorously.
Let us consider a 2D smooth image f(ρ,θ), ρ[0,1], θ[0,2π) on a compact support within a unit disk Ω . The Radon transform of f(ρ,θ) can be written as R(s,), s [–1,1], [0,π) . Suppose that we are only interested in reconstructing its interior part defined as ρ<a (0<a<1) from the corresponding local Radon transform data R(s,),|s|<a , which is also referred to as purely local projections. Based on the classic analysis, given these local data there is in general no unique solution (Natterer, 2001). Thus, any reconstructed image from such a local dataset can be viewed as a superposition of the exact reconstruction from a complete dataset R(s,),s[–1,1],[0,π) and an artificial reconstruction, which is an artifact image corresponding to another projection dataset (s,), a <|s|≤1,[0,π) . Although (s,)=0 for |s|<a, it can still produce a non-zero 2D local image g(ρ,θ), ρ[0,a), θ[0,2π) , inside the ROI. That is why we have the non-uniqueness of this interior problem (Natterer, 2001). It is well known that g(ρ,θ) is smooth and bounded inside the ROI if (s,) is continuous and bounded. Previously, we proved that both and are smooth and bounded inside the ROI (Yu and Wang, 2009). In the this note we will prove that g(ρ,θ)C for ρ[0,a) if and only if C=0.
The organization of this report is as follows. In the next section, we first prove that the constant C must be zero in the symmetric case using Taylor's expansion technique and Weierstrass's approximation theorem, and then extend it to the general non-symmetric case. In the third section, we numerically verify that a non-zero constant cannot satisfy the so-called Picard condition in the framework of singular value decomposition. In the last section, we discuss the relationship between this report and our previous work as well as future directions
Assuming that a circularly symmetric artifact image g(ρ,θ)=g(ρ) is reconstructed from a projection dataset (s,) = \(s) (a <|s|≤1,[0,π)) and (s)0 for s [–a,a] (Figure 1) and g(ρ)C (ρ[0,a) ), where C is a non-zero constant. Then, by the definition of the line integral we must have the following relationship for g(ρ) (ρ[a,1]):
Our main theoretical results are as follows.
Theorem 2.1. Assuming that a circularly symmetric artifact image g(ρ,θ) = g(ρ) is reconstructed from a projection dataset (s,)=(s) (|s|≤1,[0,π)) and (s) 0 for s [–a,a]. If g(ρ) is a square integrable function on [0,1] with g(ρ)C (ρ[0,a) ), then C=0.
Proof: Eq. (1) can be changed to
Denote (2k–1)!!=1×3×5---× (2k–1) and k!=1×2×3---×k, we have the following Taylor's expansions
Making the variable transformation t=a/ρ, Eq. (5) becomes
Eq. (6) implies
Considering the even function ĝ(t)=((|t|) on the interval [–1,1], we have
By Weierstrass's approximation theorem and the fact that the continuous function space is dense in the square integrable function space on [–1,1], there exist a polynomial series on interval [–1,1] such that
This leads to
Hence, ĝ(t)=0 almost everywhere on [–1,1], which implies C= 0.
The above theorem can be extended to the general non-symmetric case.
Theorem 2.2. Assuming that a circularly non-symmetric artifact image g(ρ,θ) is reconstructed from a projection dataset (s,) (|s|≤1,[0,π) ) and (s,)0 for s[–a,a]. If g(ρ,θ) is a square integrable function on ρ[0,1] and θ[0,2π] with g(ρ,θ)C (ρ[0,a) ), then C=0.
Proof: Let us denote g(ρ,θ) in the Cartesian coordinates as gd(x, y) with x = ρcosθ and y = ρsinθ. Define , correspondingly we have
Clearly, (ρ,θ) is a circularly symmetric square integrable function and satisfies (ρ,θ)C (ρ[0,a) ). For s [–a,a] the Radon transform (s,) of (ρ,θ) can be expressed as
where Eq.(12) has been used. By Theorem 2.1., we have C=0.
Remark: For real-world applications, all the reconstructed CT artifact images are bounded function and hence square integrable. Therefore, our assumption of square integrable functions can cover all the practical cases.
Assuming that an integral kernel K(s,ρ) defined on Is×Iρ is continuous, square integrable and non-degenerate, by means of the singular value expansion (see Chapter 1.2.1 in (Hansen, 1998) ), we have
where σ1 ≥σ2 ≥---> 0 , and ui(s) and vi(ρ) are orthonormal basis functions with respect to the conventional inner product
For any continuous and square integral function q(s), sIs, the Fredholm integral equation of the first kind
can be expressed as
If there exists solutions for Eq. (17), one of them would be in the following form:
Eqs. (18) and (19) hold for uniform convergence (see Chapter 2.4 and Theorem 8.3.2 in (Smithies, 1958)). Although K(s,ρ) and q(s) are continuous and square integrable, g(ρ) expressed by Eq. (19) is not necessarily continuous or square integrable. In order that there exists a square integrable solution g(ρ) to Eq. (17), the right side of Eq. (19) must satisfy the Picard condition (see Chapter 1.2.3 in (Hansen, 1998)):
This means that the (ui, q) must decay faster than the singular value σi. This requirement is identical to that q(s) must belong to the range of the kernel K(s,ρ).
Let , Iρ = [a, 1], Is = [0, a] and , Eq. (1) can be considered as a a Fredholm integral equation of the first kind as defined in Eq. (17). If there exists a square integrable solution g(ρ) to Eq. (1), it should satisfy the Picard condition. However, it is not easy to perform a thorough theoretical analysis to verify the Picard condition for our problem. Nevertheless, the discrete Picard condition is often used in practical applications to study the existence of a solution to Eq. (17) using the singular value decomposition (SVD) approach. As pointed out by Hansen (Hansen, 1988), although a discrete Picard condition does not exist in a strict mathematical sense, it is certainly informative when numerically solving a Fredholm integral equation of the first kind Eq.(17).
In our case, let us discrete s as , n = 1,2,---,N and ρ as , m = 1,2,---,M. Accordingly, K(s,ρ), g(ρ) and q(s) are discretized as K(sn,ρm), g(ρm), q(sn). Hence, we have
Let M=N, the matrix K has the following SVD form (Hansen, 1988):
where U=(u1 --- um --- uM) and V = (v1 --- vm --- vM) are orthogonal matrices, and Λ is a diagonal matrix whose diagonal elements satisfying σ1≥σ2≥---≥ 0. The solution to Eq. (20) can be written as
According to the Picard condition, we can compute to test the existence of a solution to Eq. (1).
Without loss of generality, we set a=0.5 and . For different combinations of sampling point numbers M and N, we computed and plotted in Figure 2, which clearly show that the singular value σm decays faster than , suggesting that there is no square integrable solution for C ≠ 0 because any non-zero constant does not lead to a sufficiently fast decay of . These numerical results are consistent to our theoretical analysis in the preceding section that if g(ρ) is square integrable constrained by an interior scan dataset then we must have C = 0 in Eq. (1).
Although we assumed that an object to be reconstructed was piecewise constant both inside and outside an ROI, the piecewise constancy is only needed inside the ROI in the proofs in (Yu and Wang, 2009; Han et al., 2009). In other words, the piecewise constancy outside the ROI is not necessary. Note that a generalized total variation minimization theorem was proved for compressed sensing based interior tomography in (Han et al., 2009), without involving the Dirac delta function. Because most objects in CT applications can be approximately modeled as piecewise constant, the piecewise constancy assumption is quite reasonable (Wang et al., 2004), and our compressed sensing based interior tomography approach is practically useful. Most importantly, our finding suggests a new direction of interior tomography, and may lead to powerful methods for better interior reconstruction. Currently, we are working to relax the piecewise constancy assumption so that real-world images can be modeled with higher accuracy for new interior reconstruction results.
As discussed in (Yu and Wang, 2009), the key idea of compressed sensing based interior tomography is to define an appropriate sparsifying transform and an associated objective function. Then, the minimization of the objective function will lead to the true ROI image without any ambiguity. For that purpose, the most important step is to prove that the artifact image inside the ROI must be zero when its total variation is zero. Actually, it might appear in (Yu and Wang, 2009) that the artifact image inside the ROI can be any constant when the total variation is zero. Hence, it is the major contribution of this report that has proved the constant being necessarily zero, providing the uniqueness of our interior reconstruction. Using this methodology, we are working to extend compressed sensing based interior reconstruction under different conditions and for other modalities.
In conclusion, we have studied the characteristic of the artifact image inside an interior ROI subject to the local projection data constraint. Our theoretical and numerical analysis has demonstrated that the artifact image inside the ROI can be expressed as a constant if and only if the constant is zero. This analysis has consolidated our proofs of Theorems 2.2 and 2.3 in (Yu and Wang, 2009).
This work is partially supported by NIH/NIBIB grants (EB002667, EB004287, EB007288), a grant from Toshiba Medical Research Institute USA, Inc., NSFC (60325101, 60532080, 60628102), Key Laboratory of Machine Perception (Ministry of Education) of Peking University, and Microsoft Research of Asia. The first author thanks Dr. Per Christian Hansen for offering a hardcopy of his paper (Hansen, 1988).