Recovery of singularities for the weighted cone transform appearing in the Compton camera imaging
Yang Zhang

TL;DR
This paper investigates the mathematical properties of the weighted cone transform in Compton camera imaging, demonstrating stable recovery of singularities and establishing microlocal stability under certain conditions.
Contribution
It provides a microlocal analysis showing the ellipticity of the normal operator and stable recovery of singularities for the weighted cone transform with specific weights and geometric conditions.
Findings
Normal operator is elliptic at accessible singularities.
Accessible singularities are stably recoverable from local data.
The analysis applies to both full and restricted cone transforms.
Abstract
We study the weighted cone transform of distributions with compact support in a domain of , over cone surfaces whose vertexes are located on a smooth surface away from and opening angles are limited to an open interval of . We show that when the weight function has compact support and satisfies certain nonvanishing assumptions, the normal operator is an elliptic DO at the accessible singularities. Then the accessible singularities are stably recoverable from local data. We prove a microlocal stability estimate for . Moreover, we show the same analysis can be applied to the restricted cone transform.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Recovery of Singularities for the weighted cone transform appearing in the Compton camera imaging
Yang Zhang
Purdue University
Department of Mathematics
Abstract.
We study the weighted cone transform of distributions with compact support in a domain of , over cone surfaces whose vertexes are located on a smooth surface away from and opening angles are limited to an open interval of . We show that when the weight function has compact support and satisfies certain nonvanishing assumptions, the normal operator is an elliptic DO at the accessible singularities. Then the accessible singularities are stably recoverable from local data. We prove a microlocal stability estimate for . Moreover, we show the same analysis can be applied to the restricted cone transform.
Partly supported by NSF Grant DMS-1600327
1. Introduction
Let be a circular cone in with vertex , central axis , and opening angle , as shown in Figure 1. We study the weighted cone transform
[TABLE]
of distributions supported in a domain in over cones of which the vertexes are restricted to a smooth surface , where is a smooth weight and is the Euclidean measure on the cone. The goal of this work is to study the microlocal invertibility of this transform.
The cone transform arises in Compton camera imaging dating back to [31, 4, 22]. A Compton camera is composed of two detectors: a scatter and an absorber. Both detectors are position and energy sensitive. When incoming gamma photons hit the camera, they have Compton scattering at various angles in the first detector and are completely absorbed in the second one. Photons can be traced back to the surface of cones.
From the setting of Compton camera, there are three important points about the transform to note here. Firstly, it is natural and unavoidable to have a weight function . Since the integral over the cone is a superposition of the line integrals, even without attenuation, there should be weight with depending on the cone. Secondly, the probability of scattering angles (i.e. the opening angle ) governed by the Klein-Nishina distribution excludes angles that are too close to [math] or . For example, in [22] the scattering angles range from to at keV incident energy and one obtains angular resolution by requiring the energy resolution to be 3.8-18.9 keV. On the other hand, when the scattering angle equals zero, the cone transform reduces to the weighted X-ray transform over a set of rays and when the scattering angle equals , it reduces to the weighted Radon transform. In both cases the inversion are easier. The limits of scattering angles can be modeled by choosing supported in such intervals w.r.t. . Thirdly, the detectors are located outside the region of interest so it is reasonable to require the vertexes are restricted to a smooth surface that does not intersect .
A lot of work has been done on the inversion of the cone transform and some of them are [2, 20, 29, 23, 17, 6, 25, 15, 14, 21, 18, 27, 19, 28, 26]. For a more complete and detailed list of previous works, see [28]. Some of the inversion formula are for special geometries, or consider the opening angle , or constant weight, or use transform over cones whose vertexes can be everywhere. Most recently, [26] considers a polynomial weight function and computes the normal operator of the cone transform over all cones whose vertexes are in to show it is a DO. It gives certain integral formula for the amplitude of the normal operator but it is not so clear how to resolve the singularities in the denominator of the integrand to obtain a smooth amplitude.
For the cone transform that we define, it is harder to analyze the Schwartz kernel of the normal operator and the normal operator may fail to be a DO in some microlocal region, if Tuy’s condition is not satisfied. Instead, we use the clean intersection calculus of FIOs in [3, 12] to show under which conditions the microlocalized normal operator is a DO of order and it is elliptic with certain nonvanishing assumptions of the weight, see Proposition 4 and 5. This approach was proposed by Guillemin in [7] for the generalized Radon transform . One difficulty to apply Hörmander’s clean intersection composition to our transform is that the composition fails to be proper. For this purpose, we modify the clean composition theorem slightly by assuming the microsupport of composed FIO is conically compact instead of the properness condition following András Vasy’s suggestion. Our first result describes the singularities of that we can recover from the data in a stable way. We note that this is the intrinsic property of the transform itself no matter what inversion algorithm is used. More specifically, let
[TABLE]
be the set of all cones that are conormal to fixed a , of which the dimension equals . We can only expect to recover singularities conormal to the cones of which the weight is nonvanishing there. The definition of accessible singularities and Tuy’s condition can be found in Section 2.
Theorem 1**.**
Suppose is accessible. If for some , then is recoverable. In other words, if is smooth near , then is smooth near .
This theorem shows that accessible singularities are recoverable with some nonvanishing assumptions on . The recovery comes from the fact that with the assumptions the normal operator in certain microlocal region is an elliptic DO and therefore it is microlocally invertible. We also shows the mapping properties of and a microlocal stability estimate when all singularities are accessible.
Theorem 2**.**
Suppose satisfy Tuy’s condition w.r.t. . Then for any , we have are continuous. For , there exists constant such that
[TABLE]
For surfaces satisfying Tuy’s condition, see Example 2. We should mention that assuming are analytic, one might be able to show the injectivity of by applying the analytic microlocal analysis, see [5]. Once we have the transform is injective on some closed subspace of , we can show the microlocal stability actually implies the stability estimate, i.e. we do not have the term in above inequality.
This work is structured as follows. Section 2 introduces the notations and definitions. Section 3 proves that is an FIO and show in which cases the projection is an injective immersion. Section 4 presents a slightly modified version of the clean composition theorem. Section 5 shows with certain positive assumptions on the weight , the normal operator is an elliptic DO at accessible singularities, which implies Theorem 1. Section 6 contains the proof of Theorem 2. Section 7 studies the weighted cone transform over cones whose vertexes are restricted to a smooth curve and the opening angle is fixed.
Acknowledgments
The author would like to thank Plamen Stefanov for suggesting this problem and for numerous helpful discussion with him throughout this project, and to thank András Vasy for helpful suggestions on the clean intersection calculus part.
2. Preliminaries
Here we introduce some notations to be used in the following sections. As mentioned before, let be an open domain in and we always assume . Suppose is a smooth regular surface without boundary where the vertexes of cones are located. Let . Let be the family of cones that we consider. Notice is a smooth manifold without boundary.
Since both are in smooth surface, we consider them in local coordinates. Suppose has the regular parameterization locally given by . Let be the Jacobian matrix. Notice form a basis for the Tangent space at of . Suppose is locally parameterized by for . Let , where and . Notice form an orthonormal basis in .
Definition 1**.**
We say is accessible by if the hyperplane conormal to has a non-tangential intersection with .
If we donate the set of all in that are accessible by , then is an open set. Observe that is accessible implies that exists a cone that it is conormal to, and for any qualified cone , the covector cannot be perpendicular to the tangent space at in .
Definition 2**.**
If every is accessible, then we say the surface satisfying the Tuy’s condition with respect to .
This coincides with the definition in [27] that Tuy’s condition is satisfied if any hyperplane intersecting has a non-tangential intersection with .
3. as an FIO
The weighted cone transform can be written as
[TABLE]
where the distribution has a nonzero factor but we can regard the factor as a part of the weight.
Proposition 1**.**
The weighted cone transform is an FIO of order associated with the canonical relation
[TABLE]
where and ; the vertex locally and is the Jacobian matrix; the unit vector is locally parameterized in the spherical coordinates and is the Jacobian matrix.
Proof.
We rewrite as the oscillatory integral
[TABLE]
where is defined as above. Notice that . Its Schwartz kernel is a Lagrangian distribution conormal to the characteristics manifold
[TABLE]
Then the conormal bundle is
[TABLE]
where we abuse the notation , to denote the differential w.r.t. the parameterization of , and we have
[TABLE]
Let
[TABLE]
One can show this is a closed conic Lagrangian submanifold of .
Additionally, since we always have is away from , from (1)(2) we get that
[TABLE]
Since , therefore the Lagrangian satisfies
[TABLE]
Then the canonical relation is given by the twisted Lagrangian. The order is given by by [10, Definition 3.2.2]. ∎
Note that we have , , and . Let , be the natural projection of to respectively. Notice in cone transform, neither of these projections can be local diffeomorphism. The following proposition describes the mapping properties of them.
Proposition 2**.**
Recall is the set of all accessible covectors in . Notice . Let be the hyperplane conormal to . We have the following statements.
- (a)
*For each , the set is a surface that can be parametrized by , where , see Figure 2 as examples. *
- (b)
*For each , there is one unique solution for the equation , which is given by (7) and below. *
- (c)
The projection restricted to is an injective immersion. In particular, if satisfies Tuy’s condition, then itself is an injective immersion.
Before the proof, the following are two examples to illustrate its first statement.
Example 1**.**
Suppose is a plane. W.L.O.G., assume . Let . If , then there are no cones that are conormal to . Here we assume . We solve
[TABLE]
to get . Thus, for , we have
[TABLE]
Example 2**.**
Suppose is a unit sphere. Note in this case satisfies Tuy’s condition w.r.t. any inside domain away from it. It can be covered by six coordinates charts. We consider the special case that and . Then
[TABLE]
In a small neighborhood of , the vertex can be parameterized by one of the following
[TABLE]
[TABLE]
Notice that and therefore can be parameterized by or from the proof above. This can also be seen from itself.
Proof of Proposition 2.
For (a), given , we are going to find out all possible solutions of from the canonical relation in Proposition 1. We have some freedom to choose , but with , the vector must be conormal to . This coincides with the fact that the singularity can only be possibly detected by the cones that it is conormal to. In other words, the possible choice of is the set . Indeed, the vertex should satisfy the equation , where . The Jacobian matrix is listed in the following,
with , . Here form a basis of the tangent space . Since is accessible, the inner products and cannot vanish at the same time; otherwise, the covector is normal to at . We simply assume in a small neighborhood of fixed . In this neighborhood, the derivative . Applying implicit function theorem, we get is a smooth function of near . Locally can be parameterized by .
Next, with given, by choosing , the axis can be determined by and . From Proposition 1,
[TABLE]
These are all smooth functions of and therefore smooth functions of . The map is an immersion. Thus, is a local parameterization of . This proves statement (a).
For (b), to recover from given , we are solving the following system of equations.
[TABLE]
Recall . We have is always nonzero and is nonzero with the assumption that . Indeed, vanishes if and only if . As we stated before, this contradicts with the assumption.
We are solving the system. Divided by and respectively, the equations (LABEL:one) and (LABEL:two) give us the projection of along , , and . These vectors form a orthogonal basis in and we can get from the projection. Plugging back into equation (LABEL:four), we can solve . Thus, when or
[TABLE]
where
[TABLE]
When or , this argument still holds if we choose another regular parameterization of there.
For (c), first we prove is an immersion. It suffices to prove that \text{\operatorname{rank}}(d\pi_{\mathcal{M}})=\dim\mathcal{C}_{I}=8. The canonical relation has a parametrization , which implies . Indeed, we can solve directly from these parameters and therefore can be represented by them. We list the Jacobian matrix with respect to these variables in the following. Here we write with and , which form a basis of ; the matrix , where and are unit vectors orthogonal to . The matrix has \text{\operatorname{rank}}=3; otherwise we have for ,
[TABLE]
which implies is normal to at . By the same arguments as before, this contradicts with the assumption. Observe that \text{\operatorname{rank}}J_{1}=\text{\operatorname{rank}}J_{2}=2. This proves that \text{\operatorname{rank}}(d\pi_{\mathcal{M}})=8.
The injectivity of comes from the statement (b). In particular, if satisfies Tuy’s condition, then . And for every , there is a unique such that
[TABLE]
∎
4. Clean Intersection Calculus
In this section, we present the clean intersection calculus stated in [12, Theorem 25.2.3]. We modify it slightly to better suit our problem.
Proposition 3**.**
Let , , be three manifolds. Suppose from to is a homogeneous canonical relation closed in and another from to . Let
[TABLE]
be properly supported FIOs with principal symbol and respectively. Assume we have
- (A1)
* intersects cleanly with excess . We call the intersection ,*
- (A2)
*there exist conically compact subsets of respectively, such that the microsupport of is in some open set of respectively. *
- (A3)
the inverse image under the projection of any is connected.
Then there exists an open set in such that is a conic Lagrangian submanifold and
[TABLE]
and for the principal symbol of we have
[TABLE]
where is the density on as is defined in [12, Theorem 25.2.3].
Remark**.**
If we compare this proposition with Theorem 25.2.3, the difference is that Hörmander assumes the restricted projection is proper instead of condition (2) above. The properness of restricted has the following implications in the original proof. First, since by [11, Theorem 21.2.14] the restricted has constant rank, if is proper and has connected fibers (i.e. the preimage of a single point is connected, that is condition (3)), then the image is an embedded submanifold, see Lemma 1 and its remark. Second, since a continuous map is proper if and only if it is closed and has compact fibers (i.e. the preimage of a single point is compact), we have is a closed submanifold and is compact. The later implies the integral is well-defined.
Proof of Proposition 3: part 1.
We prove in the following that there exists an open set such that is an embedded submanifold of .
Let be the projection
[TABLE]
Since is an open map, we can choose by choosing an open subset of by condition (2). More specifically, recall and we define
[TABLE]
Let be an open set such that . Then we have is an open subset.
We can show is an embedded submanifold in two ways. On the one hand, since restricted on is proper and , by the proof of Lemma 1, each point in has a submanifold coordinate chart. On other other hand, we can prove it by the following claim.
Claim 1**.**
The restricted projection is a closed map.
Also by Lemma 1, we have is an embedded submanifold of and therefore that of .
Then, by Claim 2, the wave front set of is contained in and thus is contained in . Now we can define the Fourier integral distributions on the embedding Lagrangian submanifold by [12, Lemma 25.1.2] and its remark. Although is not necessarily closed, by the remark in [9, p. 147], we can require the symbols (the amplitude in any specific representation) vanishing outside a closed conic subset and have the same conclusions about the principal symbols. ∎
Proof of Proposition 3: part 2.
We follow the proof of Hörmander and skip most of it here. One can first show can be written as a sum of an FIO associated with and a smoothing operator, and then compute the principal symbol. The only difference is that we still need to verify that the integration in (8) is well-defined. We claim that the principal symbol and have conically compact support and therefore it is well defined. Indeed, consider the local representation of
[TABLE]
where are non-degenerate phase functions and are amplitudes. For more details see [12, Theorem 25.2.3]. Then the principal symbols of are
[TABLE]
according to [12, p. 14]. Since have conically compact microsupport respectively, it suffices to show that the principal symbol vanishes outside the microsupport. Indeed, for fixed (x_{0},y_{0},\xi^{0},\eta^{0})\notin\text{\operatorname{WF}}(A_{1}) in the conical relation, there exists a small conic neighborhood of such that the local representation of above is smooth. By [30, Lemma 4.1], there is such that on . On the other hand, the principal symbol defined in (9) only depends the amplitude restricted to the the manifold . Therefore, it vanishes in this neighborhood.
∎
Lemma 1**.**
Let be a smooth map with constant rank. If is closed and the preimage of any is connected, then is an embedded submanifold of .
Remark**.**
In particular, this lemma holds if we assume is proper instead of simply closed. A slightly different version of this lemma can be found in [8]. It claims the proof can be found in Appendix C.3 in [11]. The proof we write in the following is based on the outlines of Hörmander and we borrow most of it from an online answer that proves the case when is proper.
Remark**.**
It is necessary to assume that is closed. Consider projection from to the plane. Let be a smooth curve as is shown in Figure LABEL:counterexamples (a), whose image under projection is the figure . Here is the projection restricted on . The preimage for each point in under is connected but is not closed. The image is not an embedded submanifold.
It is also necessary to assume that each fiber of is connected. A counterexample shown in Figure LABEL:counterexamples (b) is constructed from the immersed manifold ”figure ”. By lifting it into and smoothly extended the two ends, we get a smooth curve in and is defined as before. Notice is a closed map but the preimage of certain point is not connected.
Proof of Lemma 1.
We followed the suggestions posted in [13] to prove this lemma. The proof can be divided into three steps.
Step 1. By [11, C 3.3], the constant rank of shows is locally a submanifold. That is, for every and , there exists an open neighborhood of and an open neighborhood of such that is a submanifold of . It suffices to prove that there is an open such that , thus we have is a submanifold of .
Step 2. Since is connected, for any , we can show maps the neighborhood of to the ”same” neighborhood of by defining an equivalence relation. More specifically, we say if for any open neighborhood of and of , there exists open neighborhood and such that . Observe that each equivalence class is an open set, and different classes are disjoint. Therefore, if there are more than one equivalence class, then they form partition of , which contradicts with the connectedness.
Step 3. Back to what we want to prove, assume there is no such exists. Then there exists a sequence converging to with each point distinct from but . Pick arbitrary preimage of . When is closed instead of proper, we have two cases. If has limit points in , then we can choose a subsequence , which is forced to be in . By Step 2, there should be and such that . For large enough , we have and therefore , which contradicts the assumption. If has no limit points, then is a closed subset . The set has a limit point which is not contained in . This contradicts with the assumption that is closed. ∎
Claim 2**.**
If the microsupport \text{\operatorname{WF}}(A_{1}) and \text{\operatorname{WF}}(A_{2}) are conically compact, then the wave front set \text{\operatorname{WF}}(A) is conically compact and is contained in .
Proof.
By [10, Theorem 8.2.14], since \text{\operatorname{WF}}(A_{1}),\text{\operatorname{WF}}(A_{1}) are away from the zero sections, then we have
[TABLE]
where \text{\operatorname{WF}}^{\prime}(\cdot) is the twisted relation. Hence, \text{\operatorname{WF}}^{\prime}(A) is a closed conic set contained in the image of the projection from the intersection (\text{\operatorname{WF}}^{\prime}(A_{1})\times\text{\operatorname{WF}}^{\prime}(A_{2}))\cap(T^{*}X\times T^{*}\Delta(Y)\times T^{*}Z) to . The intersection is conically compact so it is closed. Then the projection restricted there is continuous and it maps compact set to compact set. Moreover, it commutes with the multiplication by positive scalars in the covariant variables, and therefore the image of the intersection is conically compact. Thus, we have \text{\operatorname{WF}}^{\prime}(A) is conically compact. ∎
Proof of Claim 1.
Since is an open set in , any closed subset is the intersection of with some closed subset of . Notice is conically compact and preserves the fiber, which implies restricted to is proper (therefore, is closed). It follows that is closed and therefore is a closed subset of . ∎
Clean Composition
A special case is when we compose the operator with its adjoint. As the following lemma shows, with certain condition the composition is clean.
Lemma 2**.**
Suppose from to is a homogeneous canonical relation closed in . Let be a properly supported FIO associated with . If the projection is an injective immersion, then the composition is clean. In particular, the canonical relation of the composition map is the identity.
We follow the same arguments in [16].
Proof.
From [12, Thm.25.2.2], we have . By definition, the composition is clean if is a smooth manifold and its tangent space equals to the intersection of the tangent space of the intersecting manifolds. Indeed, let , and be the interchanging map. We have
[TABLE]
Since the projection is injective, then . This implies
[TABLE]
On the other hand, we have is contained in the tangent space of , if and only if
[TABLE]
Since is an injective immersion, it follows that and , which proves the lemma. ∎
Remark**.**
In the setting of this lemma, the connectedness condition (A3) in Proposition 3 is not needed. Observe that from the proof of Proposition 3 and Lemma 1, the connectedness of is required to guarantee that the composition does not intersect itself. However, in this case we have the composition is the diagonal of , which is automatically not self-intersecting.
5. The normal operator as a DO
In order to apply the clean composition theorem in Proposition 3 to , we need the composition satisfying three assumptions (A1), (A2), (A3).
For (A1), by Proposition 2 and Lemma 2, if is accessible, then the composition is clean. For (A3), see the remark after Lemma 2. In most case in application, the surface is a plane, which makes the situation simpler. As for (A2), we would like to show that with certain assumptions on the support of (or by choosing proper smooth cutoff functions), we can find a compact subset of such that the microsupport of (or multiplied by cutoffs) is contained in some open subset of .
Lemma 3**.**
Let be smooth cutoff functions with compact supports. Then is a Lagrangian distribution with conically compact microsupport support in . Additionally, there exists a compact set such that \text{\operatorname{WF}}(\chi_{2}I_{\kappa}\chi_{1}) is contained in some open subset of . In particular, these statements are true when has compact support in .
Proof.
As is shown in the proof of Propostion 2 (c), we have is a parameterization of . Thus, the map
[TABLE]
is a continuous submersion. Since we have continuous functions map compact sets to compact sets and submersions are open maps that map open sets to open sets, this implies maps compact (or open) set in to compact (or open) set in .
In fact, we have maps conically compact (or open) set to conically compact (or open) set. Indeed, suppose we have a conically compact neighborhood in , we can modify it to a compact neighborhood by restricting . Then the image of the compact set is compact in and therefore is compact in restricted to the cosphere bundle. Since is homogeneous of order 1 w.r.t. , we have the image of the conically compact (or open) set is conically compact (or open).
For the first statement of this lemma, notice that is compact and we have the compact supports w.r.t . Since might not be the whole range of , we additionally assume we have compact support w.r.t . To show the existence of the compact set , observe that for any compact set in or , we can find a larger compact set such that is in some open set of this larger compact set. By the arguments above, it proves the second statement. ∎
With the lemma above, assuming has compact support, we can apply Proposition 3 to to show it is a DO. Moreover, with additional assumptions it is an elliptic DO, according to the formula 8 for composed principal symbols.
Proposition 4**.**
Assume is accessible. Suppose the weight function has compact support in . If we have for some , then is a DO of order elliptic at .
Remark**.**
The order is calculated by , where the excess equals to the dimension of . By Proposition 2, we have , of which the dimension is .
If it is not true that the whole cotangent bundle is accessible, then we can have a microlocal version Proposition 4 by choosing proper cutoff DOs. More specifically, for a fixed accessible covector , there exists a conic neighborhood of such that each covector in this neighborhood is accessible. By choosing a cutoff DO which is supported in this neighborhood, then we can prove the following proposition.
Proposition 5**.**
Suppose is accessible. Suppose for some covector in . Then there exists a DO of order zero elliptic at with microsupport in a conically compact neighborhood that is accessible, such that for any that is a DO of order zero elliptic at with conically compact support, the microlocalized normal operator is a DO of order elliptic at .
Theorem 1 is the direct result of this proposition and the following example is a special case when we have unrecoverable singularities.
Example 3**.**
Let the surface of vertices be an open disk with radius in the plane. Let be an open domain above plane.
Consider a covector with or nonzero. Otherwise, it is not accessible. There is a cone with vertex in conormal to only if
[TABLE]
[TABLE]
[TABLE]
If we have , we can solve from the first equation, and plug it in the second one to have
[TABLE]
For simplification, we denote by to get
[TABLE]
This is a parabola opening to the top. There exists a solution of if and only if
[TABLE]
which implies
[TABLE]
If , then we must have . By symmetry we should get the same result. Notice, if we can find such in that , then we can construct a cone conormal to by properly choosing and . Thus, the set of the unrecoverable singularities is
[TABLE]
6. proof of theorem 2
For the estimates, we have the following restatement of [Hoermander1791, Theorem 4.1.9 and Thmeorem 4.3.2], which indicates the mapping properties in the general case.
Proposition 6**.**
Let be a homogeneous canonical relation from to such that
- (1)
the maps and have subjective differentials,
- (2)
the projections and have constant rank.
Let and respectively. Then . Suppose is properly supported. Then
[TABLE]
are continuous, for any and .
Proof.
This corollary is a direct result of the theorems above. Let be the DO of order such that . We apply [9, Thm.4.3.2] to . Then continuously, for any . Notice if satisfies the conditions (1)(2) in above corollary, then satisfies them as well. Since , we can apply the same argument. ∎
Now we consider the conical Radon transform and its canonical relation . We have the following claim by Proposition 6.
Claim**.**
When is accessible, the canonical relation satisfies condition (1) and (2). In particular, we have and the inequality in Corollary 6 is .
Proof.
By Proposition 2, the projection is an injective immersion, which implies is a submersion. Thus, . From the proof of Proposition 2, we have is parameterized by . It is obvious that has subjective differential. For the projection , the Jacobian is in the following,
where . Since is solved from , differentiating w.r.t. we have
[TABLE]
Thus, the vector is nonzero and the differential of projection has rank equals . ∎
By the claim above, we can prove Theorem 2 in the following based on similar arguments in [24, 1]
Proof of Theorem 2.
By Corollary 6 and the claim above, we have the second inequality and for some constant . From Proposition 4, we have Combining these two we have desired result. ∎
Another proof of the estimate for .
First we prove is bounded, for , where is the space of all supported in . Indeed, we have
[TABLE]
Here we have
[TABLE]
where are multi-indexes and . Since we always have , locally we can write as , where is the index such that . Therefore by integration by parts and induction, we get the similar integral transform of derivatives of up to order with a new weight function. This implies is a DO of order . Thus, we have the estimates
[TABLE]
When , the proof is simplified. For non-integer , we can use interpolation methods. Then by duality, we have is bounded, for and . Since we always assume is compactly supported, we have has support in . Thus,
[TABLE]
Combining these two inequality and the ellipticity of , we have for ,
[TABLE]
We abuse the notation to denote different constants. ∎
7. Restricted Cone Transform
In this section, suppose is a smooth regular curve that is parameterized by . For fixed , we define the restricted cone transform as
[TABLE]
We have the following corollaries.
Corollary 1**.**
The restricted cone transform is an FIO of order associated with the canonical relation
[TABLE]
where and ; the vertex ; the unit vector is parameterized in the spherical coordinates and is the Jacobian matrix defined as before.
Corollary 2**.**
Let be the set of all in that are accessible w.r.t. . Then is an open set. Let . We have the following properties.
- (a)
*For every , there is one unique solution for the equation , which is given by (11) and (12) . *
- (b)
The projection restricted to is an injective immersion. In particular, if satisfies Tuy’s condition, then itself is an injective immersion.
Proof.
For (a), notice we have . By (7) and its proof, when or ,
[TABLE]
where
[TABLE]
We still need to solve from . Note that , since is accessible. Therefore, we have
[TABLE]
For (b), to prove that restricted to is an immersion, it suffices to show that is parameterized by . Indeed, from (b), can be represented by . By writing and performing a same argument as in the proof of Proposition 2, we can show it is a parameterization and therefore the rank of the differential of equals . ∎
Remark**.**
Based on similar arguments, we can also show the same results hold if we only restrict the vertexes on a smooth curve without fixed. Additionally, one can show an analog of recovery of singularities in both cases as in Theorem 1.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Yernat Assylbekov and Plamen Stefanov. A sharp stability estimate for the geodesic ray transform.
- 2[2] Roman Basko, Gengsheng L Zeng, and Grant T Gullberg. Application of spherical harmonics to image reconstruction for the Compton camera. Physics in Medicine & Biology , 43(4):887, 1998.
- 3[3] J. J. Duistermaat and V. W. Guillemin. The spectrum of positive elliptic operators and periodic bicharacteristics. In Mathematics Past and Present Fourier Integral Operators , pages 243–283. Springer Berlin Heidelberg, 1994.
- 4[4] D.B. Everett, J.S. Fleming, R.W. Todd, and J.M. Nightingale. Gamma-radiation imaging system based on the Compton effect. Proceedings of the Institution of Electrical Engineers , 124(11):995, 1977.
- 5[5] Bela Frigyik, Plamen Stefanov, and Gunther Uhlmann. The x-ray transform for a generic family of curves and weights. Journal of Geometric Analysis , 18(1):89–108, nov 2007.
- 6[6] Rim Gouia-Zarrad and Gaik Ambartsoumian. Exact inversion of the conical Radon transform with a fixed opening angle. Inverse Problems , 30(4):045007, mar 2014.
- 7[7] Victor Guillemin. On some results of Gel’fand in integral geometry, 1985.
- 8[8] Victor Guillemin and Shlomo Sternberg. The moment map revisited. Journal of Differential Geometry , 69(1):137–162, jan 2005.
