Injectivity and Stability for a Generic Class of Generalized Radon Transforms

Let (M, g) be an analytic, compact, Riemannian manifold with boundary, of dimension \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n \ge 2$$\end{document}n≥2. We study a class of generalized Radon transforms, integrating over a family of hypersurfaces embedded in M, satisfying the Bolker condition (in: Quinto, Proceedings of conference “Seventy-five Years of Radon Transforms”, Hong Kong, 1994). Using analytic microlocal analysis, we prove a microlocal regularity theorem for generalized Radon transforms on analytic manifolds defined on an analytic family of hypersurfaces. We then show injectivity and stability for an open, dense subset of smooth generalized Radon transforms satisfying the Bolker condition, including the analytic ones.


Introduction
Let (M, g) be an analytic, compact, Riemannian manifold with boundary, of dimension n ≥ 2, with volume form denoted by dVol. Let be a family of embedded hypersurfaces. A generalized Radon transform takes each f ∈ C ∞ (M) to the set of its integrals over the hypersurfaces of , with respect to the surface measure induced by the volume form. Often is itself a smooth manifold; for example, the Euclidean Radon transform is defined over the set of affine hypersurfaces in R n , which has R × S n−1 as a double cover. In the current paper, we assume that is a smooth manifold of dimension n. Such transforms are found in applications to many other fields, including harmonic analysis, scattering theory, seismology and medical imaging.
The main questions regarding these transforms include determining conditions under which they are injective, finding when the transform has a stable inversion, and characterizing the range. We concentrate on the first two questions here. These questions are also important in the partial data case, where integrals are known only for a subset of . In the case of the Euclidean Radon transform, we refer the reader to [8,16,17,23] and references therein for the resolution of these problems, and their generalization to the Radon transform over symmetric spaces and other related contexts.
The geometric data of a generalized Radon transform can be encoded by an incidence relation between points on M and the hypersurfaces in that contain them. Let ⊂ M × be this relation, i.e., the set of ordered pairs (x, σ ) such that x ∈ σ . One says is a double fibration when it is a smooth, embedded submanifold of M × such that both canonical projections are smooth and their restrictions to form a fiber bundle over M and , respectively [10].
Guillemin and Sternberg [12][13][14] showed that given a Radon transform R defined by a double fibration , both R and its adjoint R * (often called the generalized backprojection by analogy with the Euclidean case) are Fourier integral operators, and the canonical relation of R is the conormal bundle N * of the incidence relation. If in addition satisfies the Bolker condition, which says that the induced projection π * : N * → T * is an embedding, then the normal operator R * R is an elliptic pseudodifferential operator, which yields invertibility up to smoothing error.
A stronger result than invertibility is a Helgason-type support theorem, by analogy with that of the Euclidean Radon transform [15]. Such a support theorem implies that if f is a priori of compact support, and R f = 0 for all hypersurfaces intersecting the support, then f = 0. In the analytic category, there has been much work in this area (for example, [3,4,25,26]) using analytic microlocal analysis and the Bolker condition to prove support theorems for analytic generalized Radon transforms (i.e., with M, analytic manifolds) with nonvanishing analytic weight. On the other hand, in n = 2 there is a counterexample in the smooth category due to Boman [2] of a function supported in the disk such that some weighted Radon transform over lines vanishes. For the weighted X-ray transform over curves, there is an analogous Bolker condition [11] and support theorems are known for a class of such transforms in n ≥ 3 [24,31,34] including the geodesic ray transform on an analytic, simple manifold over functions [21] and over symmetric tensor fields [22].
Our first result considers such analytic generalized Radon transforms and shows their analytic microlocal regularity. This builds upon similar results for the weighted X-ray transform over a generic class of curves [7,9] and in particular, geodesics [29,30]. To avoid complications at the boundary of M, we embed it isometrically in a slightly larger, open analytic manifold M 1 . In Sect. 2, we show how to extend the definition of R w to a transform on M 1 in a stable way. We show: Theorem 1 Let R w be an analytic generalized Radon transform satisfying the Bolker condition, with w an analytic, nonvanishing weight. Let f ∈ E (M 1 ) be such that R w f (σ ) = 0 in a neighborhood of some hypersurface σ 0 ∈ . Then the analytic wavefront set WF A ( f ) does not intersect the conormal bundle N * σ 0 .
The main tool is a complex stationary phase lemma of Sjöstrand [28, Theorems 2.8, 2.10 ff.] and related techniques, which suffice in lieu of a hypothetical analytic calculus of Fourier integral operators. The proof of this theorem is given in Sect. 3.
Theorem 1 also implies a local support theorem [25,Prop. 2.3] and in particular the injectivity of R w : L 2 (M) → H (n−1)/2 ( ). In fact, the proof of the theorem shows that it suffices for R w f (σ ) to be analytic in a neighborhood of σ 0 . One then obtains a unique continuation result of the following type: if f ∈ E (M 1 ) is analytic on one side of σ 0 , and R w f (σ ) is analytic in a neighborhood of σ 0 , then there is a neighborhood of σ 0 on which f is analytic.
Our second result is a stability estimate for a generic class of smooth generalized Radon transforms satisfying the Bolker condition. We restrict ourselves to those generalized Radon transforms studied by Beylkin [1], which have parameterized globally by the level sets of a smooth defining function ϕ satisfying some conditions to be made explicit later.
Theorem 2 Let (M, g) be an analytic Riemannian manifold with boundary. Take R w : L 2 (M) → H (n−1)/2 ( ) to be an injective generalized Radon transform defined by ϕ with weight w, satisfying the Bolker condition. Then there exists K n and a neighborhood of (ϕ, w) ∈ C K such that the generalized Radon transformRw defined on (M, g) by a defining function and weight in this neighborhood is injective and for all f ∈ L 2 (M) there exists C > 0 such that This follows from an analysis of the symbol of the normal operator R * w R w . As mentioned above, under the Bolker condition it is an elliptic pseudodifferential operator. We show that perturbing the defining function and weight slightly in C K perturbs the operator slightly, preserving the stability estimate.
While we work entirely on an analytic Riemannian manifold (M, g) and do not perturb the metric in this result, we use the metric only to provide a convenient choice of surface measure, and to ensure the existence of a dense set of injective, generalized Radon transforms. We may then conclude: Corollary 1 On each analytic, compact Riemannian manifold with boundary, there is a generic set of generalized Radon transforms satisfying the Bolker condition that are both injective and stable.
We defer the proof of Theorem 2 to Sect. 4.

Generalized Radon Transforms
In this section we fix notation and establish some basic facts about the generalized Radon transform, including a statement of the Bolker condition. For concreteness we consider the space of hypersurfaces as parameterized by a defining function, following Beylkin [1], though we only consider oriented hypersurfaces. To avoid difficulties occurring at the boundary of M, we assume M is isometrically embedded in a slightly larger open manifold M 1 with compact closure (i.e., M 1 is compact), whose metric we also refer to by g. If we are considering the analytic category of Radon transforms, we will also assume M 1 is analytic. In the sequel, we will always consider L 2 (M) to be functions on M 1 , extended by zero.
). ϕ is a defining function when it satisfies the following conditions: The level sets of ϕ will be denoted by Note that by homogeneity, H s,θ = H λs,λθ for λ > 0. Therefore we can consider as globally parameterized by (s, θ) ∈ R × S n−1 , which is of dimension n. Often we will also implicitly consider ϕ as a function on M 1 × S n−1 . Moreover, the map (x, θ) → (x, ϕ(x, θ), θ) ∈ M 1 × provides local coordinates on the incidence relation , which implies that the projection π : → M 1 has compact fiber S n−1 , and π is proper.
The third condition imposed on a defining function is a local form of Bolker's condition. This allows us to locally identify x M 1 . We will assume in addition a stronger, global Bolker condition.
The first condition is roughly analogous to the "no conjugate points" condition assumed by [9,21] for similar results regarding the geodesic ray transform, and the second ensures that every singularity is observable from some hypersurface in . Note that generalized Radon transforms defined by a double fibration satisfying the Bolker condition as stated by Guillemin et al. also satisfy this Bolker condition; see [25,Lemma 3.5].
Given h ∈ C ∞ (M), let the generalized Radon transform R w determined by (M, g, ϕ, w) be defined by is the extension of h to M 1 by zero, and dμ s,θ is the volume form on H s,θ induced by dVol. There exists a smooth, nonvanishing function J (x, θ) such that We now calculate a formula for the adjoint R * w . Given f ∈ C ∞ (M) and g ∈ C ∞ c (R × S n−1 ), we derive the following formula for R w f, g L 2 (R×S n−1 ) : Therefore, we may define This is simply a generalized backprojection with weight w J , and by differentiating under the integral sign we see that it maps

Microlocal Regularity
In this section, we take (M 1 , g) to be an analytic Riemannian manifold, ϕ to be an analytic defining function, and w to be an analytic nowhere vanishing weight. Given f ∈ E (M 1 ), we are interested in the microlocal analyticity of f given that of R w f . (We extend R w to E (M 1 ) by duality.) We will use the following definition of the analytic wavefront set, following Sjöstrand. There are alternative approaches to analytic wavefront set by Sato et al. [27] and also Bros and Iagolnitzer [6], which were shown to be equivalent by Bony [5].
Definition ([28, Def. 6.1]). Let (x 0 , ξ 0 ) ∈ T * R n \ 0 and let ψ(x, y, ξ) be an analytic function defined in a neighborhood U of (x 0 , x 0 , ξ 0 ) ∈ C 3n such that (2) There exists C > 0 such that for all (x, y, ξ) ∈ U , we have Let a(x, y, ξ) be an elliptic classical analytic symbol defined on U ; see, e.g., [28,Theorem 1.5]. We for some C > 0, uniformly in a conic neighborhood of (x 0 , ξ 0 ). The analytic wavefront set is the closed conic set WF A (u) ⊂ T * R n \ 0, which is the complement of the set of covectors near which u is microlocally analytic.
We note that this definition is microlocal and invariantly defined, and therefore can be extended to distributions on analytic manifolds (see [20,Theorem 8.5.1] and the remarks following). In this case, for u ∈ D (M 1 ), WF A (u) is a closed conic subset of T * M 1 \ 0.
Proof Let us fix a coordinate system. We already have local coordinates (x, θ) on T * M \ 0. Without loss of generality we can take s 0 = 0 and |θ 0 | = 1. To simplify the coordinates on , we perform a stereographic projection onto the tangent plane of the sphere at θ 0 , which is an analytic diffeomorphism mapping a neighborhood of θ 0 ∈ S n−1 to a neighborhood of the origin in R n−1 . We refer to the coordinates on this tangent plane by ξ , and pass to a perhaps smaller neighborhood of with |s| < 2 and |ξ | < δ, with , δ > 0 being small parameters.
Much of the complexity of analytic microlocal calculus is due to the difficulty of localizing in the analytic category, as there are no suitable cut-off functions. Instead one often uses a sequence of quasianalytic cut-off functions χ N ∈ C ∞ c (R), depending on , for whose construction we refer to [18,33]. We will only use the following properties of this sequence: (1) supp χ N ⊂ (−2 , 2 ) and χ N (− , ) = 1.
(2) For all N ∈ N and k ≤ N , the estimate holds for a constant C > 0 independent of N .
By assumption R w f (s, ξ) = 0 for |s| < 2 and |ξ | < δ. Let λ 1 be a large parameter, to be fixed later. This implies that Recall that ξ are analytic coordinates for the neighborhood of θ 0 in S n−1 that we are concerned with, and so here and in the sequel we write for brevity, e.g., w(x, ξ) = w(x, θ (ξ )) and dμ s,ξ = dμ s,θ(ξ ) . It follows from Beylkin's construction that where J (x, ξ) is an analytic, nonvanishing Jacobian and dVol is the volume form on M 1 associated with the metric. Hence (1) reduces to the oscillating integral Here a N (x, ξ) is a sequence of classical analytic symbols on the same neighborhood of (x 0 , 0) ∈ M 1 × R n−1 . The coordinates on x and ξ are real-analytic, and so we may extend their domain of definition slightly by analytic continuation to a Grauert tube of a small neighborhood of H 0,0 ⊂ M 1 (for x) and a small neighborhood of the origin in C n−1 for ξ . This continuation in principle depends on the choice of analytic coordinates, but as the analytic wavefront set is invariantly defined the final result does not depend on this choice. We choose a perhaps smaller δ such that {ξ ∈ C n−1 : |ξ | < δ/2} is contained in this neighborhood. We denote the local complex coordinate patch of x 0 as U ⊂ C n . Let y ∈ U and η ∈ C n−1 , with |η| < δ/2. Let ρ(ξ ) = 1 when |ξ | ≤ δ and zero otherwise. Then we multiply (2) by and integrate with respect to ξ . The resulting integral is of the form Here b N is a sequence of classical analytic symbols defined on a complex neighborhood of H 0,0 × {0} × {0} and is the augmented phase function given by To estimate the left-hand side of (3), we intend to use the method of complex stationary phase. Therefore, we are interested in the critical points of the function ξ → (x, y, ξ, η). Note that There are clearly real critical points ξ when ξ = η and x = y. These critical points are non-degenerate, and therefore induce complex critical points ξ c (x, y, Consider the situation when y = 0. Then for x = 0, the only real critical points are where ∂ ξ ϕ(x, ξ) = ∂ ξ ϕ(y, ξ). However, this cannot happen by the global Bolker condition that we imposed on the defining function. By non-degeneracy again we see there are no real or complex critical points other than ξ c (x, y, η) for (x, y, ξ, η) where |y| < δ and |ξ − η| < δ. Now we apply the complex stationary phase lemma [28, Theorems 2.8, 2.10] to (3). As a preparatory step divide the integral into two regions; one over the region y, ξ, η) : |x − y| ≤ δ/C 0 , |ξ − η| < δ} and one over the region Here C 0 > 0 is a constant chosen so that the critical points ξ c (x, y, η) lie within I + and none lie in I − .
In I − , we may define the usual operator L such that Le iλ = e iλ via This is well defined as there are no critical points in I − , so we may repeatedly integrate by parts: The terms B k are boundary terms that decay exponentially, due to the fact that Im > 0 for |ξ − η| = O(δ). As for the integral on the right-hand side, we recall that b N is defined by The worst possible growth of (L * ) N b N in terms of N occurs when all derivatives are applied to χ N (ϕ(x, ξ)), and in this case we may apply the estimate which follows from the construction of the sequence of quasianalytic cut-off functions. Therefore, As for the integral over I + , the cut-off functions χ N (ϕ(x, ξ)) are all equal to one. Therefore the amplitude on I + does not depend on N ; we remove this dependence and refer to the amplitude restricted to this region as b. We know all of the critical points of ξ → and can therefore apply the complex stationary phase lemma. This yields an estimate of the form Here ψ (x, y, η) = (x, y, ξ c (x, y, η), η) and B(x, y, η) = b(x, ξ c (x, y, η), η). We may now fix N such that N ≤ (λ/Ce) ≤ N + 1 to ensure the error is exponentially small. (x, y, η) = 0 and ψ(x, y, η) is a non-degenerate phase function. To show this implies (x 0 , θ 0 ) / ∈ WF A ( f ), we check the details of the characterization of the analytic wave front set given above. Recall ξ c (x, y, η)).

Now B(x, y, η) is an elliptic analytic symbol near
Note that ξ c (x, x, η) = η for x real, and therefore ψ(x, x, η) = 0. In addition By the global Bolker condition we can make a change of variables η so that η = d x ϕ(x, η). Finally, it is clear that Im ψ(x, y, ξ) ≥ C|x − y| 2 for x, y real. Therefore, Theorem 1 follows from applying the proposition to all conormals of a fixed hypersurface σ 0 .
Remark From the proof we see that it suffices for R w f (σ ) to be analytic in a neighborhood of (s 0 , θ 0 ). After microlocalization, the right-hand side of (2) will be O(e −λ/C ) instead of zero, but this poses no problem.

Stability
We now return to generalized Radon transforms with smooth defining function ϕ : M 1 × S n−1 and smooth, nonvanishing weight w : M 1 × S n−1 . The object of interest in this section is the normal operator N w = R * w R w . Note that as mentioned in Sect. 2, the projection π : → M 1 is proper, thus R * w and R w can be composed. It is known that the global Bolker condition implies that N w is a pseudodifferential operator [13,Prop. 8.2], at least when the weight w(x, θ) satisfies symbol estimates of order zero in x. Here we have considered R w as acting on functions in L 2 (M) after the latter have been extended by zero; therefore, we are free to modify the weight outside a neighborhood of M (e.g., by cutoff function) to satisfy the symbol estimates, without affecting the action of R w on L 2 (M).
We require detailed knowledge of the symbol of N w for the kind of stability estimates we prove later. First we obtain a representation of the Schwartz kernel of R w .
Proof We perform a partial Fourier transform of R w f (s, θ) in the s variable, taking s to be the dual variable of s. The change of variables then yields Therefore −ϕ(y,θ))s w(y, θ)J (y, θ) f (y) dVol(y) ds Similarly, the kernel of the generalized backprojection R * w is From this we see that the kernel of N w is We can now use this representation to find the principal symbol of the normal operator N w .
where W is the auxiliary function Proof Beginning from (6), we split the integration over R into {s > 0} and {s < 0}.
Using the positive homogeneity of the defining function, we rewrite the integral as Here K + N w and K − N w are the Schwartz kernels of the operators N + w and N − w respectively, so that N w = N + w + N − w . We work with each term separately. Let ξ = s θ be polar coordinates for R n . This change of variables is justified when the kernel is applied to a test function in C ∞ c (M 1 ); using the proof of [20,Theorem 7.8.2] it can be shown that it is justified for the kernel itself. Then we obtain By the global Bolker condition, ∂ ξ ϕ(x, ξ) = ∂ ξ ϕ(y, ξ) implies x = y. A stationary phase argument implies that K + N w is a smooth function away from the diagonal of Fix x 0 ∈ M 1 . There exists a neighborhood U of x 0 on which we have normal coordinates, which we refer to again with (x i ), such that x(x 0 ) = 0. We then use (x i , y i ) as coordinates on U × U , with x i = y i . We consider the localized kernel In these local coordinates, we can expand the phase function near the diagonal x = y.
Near the diagonal, this map is smooth, and and so the map (x, y, ξ) → (x, y, ξ ) is a diffeomorphism of a neighborhood of the diagonal onto another neighborhood of the diagonal. It may be necessary to shrink the support of χ slightly for the change of coordinates to be well defined. In these variables, Both sides are positive homogeneous of degree one, which implies that |ξ | = c(x, y)|ξ | with c(x, y) a strictly positive, smooth function defined near the diagonal. Clearly c(x, x) = |∂ x ϕ(x, ξ)|. This reduces the cut-off kernel of N + w to an honest pseudodifferential operator To evaluate the principal symbol of K + N w , we restrict the amplitude to the diagonal x = y. The principal symbol of N w is the sum of those for N + w and N − w .
This reconfirms that N w is an elliptic pseudodifferential operator of order 1 − n, provided the weight is nonvanishing and the global Bolker condition is satisfied.
We now consider the stability of reconstructing f ∈ L 2 c (M 1 ) from R w f , using the analysis of the normal operator N w f = R * w R w f from the previous section. The basic estimate follows from elliptic regularity; the stability estimate follows from [32, Prop. V.3.1].

Lemma 3
Let w ∈ C ∞ (M 1 × S n−1 ) be a nonvanishing weight and let ϕ ∈ C ∞ (M × (R n \ 0)) be a defining function. Let f ∈ L 2 (M) be extended by zero tof ∈ L 2 (M 1 ). Then for all and s > 0 there exists C > 0 and C s > 0 depending on s such that If, in addition, N w : L 2 (M) → H n−1 (M 1 ) is injective, then we have a stability estimate with a loss of n − 1 derivatives, In particular, by Proposition 1, the latter stability estimate holds when the geometric data (i.e., M, ϕ and w) are analytic. This also follows directly from [3]. Our main contribution is to extend this stability estimate by perturbation to a generic set of smooth geometric data. We begin by using the standard pseudodifferential calculus to show that the normal operator depends continuously on finitely many derivatives of the data. (M 1 , g) be an open Riemannian manifold with an embedded compact manifold M with boundary. Let ϕ 1 , ϕ 2 be two defining functions and w 1 , w 2 be two nonvanishing weights. Let N 1 = R * w 1 R w 1 and N 2 = R * w 2 R w 2 . Let f ∈ L 2 (M) and let f ∈ L 2 (M 1 ) denote its extension by zero. There exists a K n such that if

Lemma 4 Let
then there exists C > 0 depending a priori on the C K (M 1 × S n−1 ) norm of ϕ 1 and w 1 such that Proof We have seen in the previous lemmas that N 1 and N 2 are both elliptic pseudodifferential operators with symbols depending on ϕ 1 , ϕ 2 and w 1 , w 2 respectively. Let K be an arbitrary, large natural number to be fixed later. If the defining functions and weights are δ-close in C K (M 1 ), then it follows from Lemma 2 that the amplitudes are O(δ) in C K −2 (M 1 ). By the continuity of pseudodifferential operators [19, Theorem 18.1.6] the operator norm of N ± 1 − N ± 2 is bounded by a constant multiplied by some C ∞ (M 1 )-seminorm of the difference of the amplitudes. Take K large enough so that The lemma follows from adding the positive and negative parts of the estimate together. Notice that K , the necessary number of derivatives, does not depend on the defining functions or the weights themselves.
It is of interest to determine the minimal regularity necessary for the perturbation result of the previous lemma. For the related geodesic ray transform, it is known that the geometric data need only be δ-close in C 2 [9]. One would then expect the above to hold for data δ-close in C n .
We now able to prove with our main result, Theorem 2.
Proof of Theorem 2 Recall R w is an injective generalized Radon transform. Lemma 3 yields the following stability estimate: wheref is the extension of f ∈ L 2 (M) to M 1 by zero. Then Lemma 4 allows us to perturb this estimate using Therefore, For δ < min{(2C 1 C 2 ) −1 , 1/2}, the second term on the right-hand side may be absorbed into the left. The resulting stability estimate for the perturbed normal operator implies injectivity ofRw.