These notes are based on Gunther Uhlmann’s lectures for MATH 581 taught at the University of Washington in Autumn 2009.
An index to all of the notes is available here.
Now we are prepared to study the range of the X-Ray Transform. We will will look at the image of in detail — once we understand this we can easily use our understanding of the Radon transform to characterize the image of for smooth, compactly supported functions.
So let . Recall that
It will be convenient to extend to all of as follows
So is said to be positive homogeneous of degree -1 in .
Claim 1 For , The partial differential equations
hold for all .
To show this, check the following
which is symmetric in and .
These equations are called John’s Equations, named for Fritz John who was the first to point this fact out. The main result is that these are the only conditions we need to characterize the range of
Theorem 1 (Fritz John) If , define its extension to as
Then is in the range of if and only if and
This is nice: we characterize the range with only a finite number of conditions. Contrast this with the range characterization of the Radon Transform, which required an infinite number of conditions and could not be effectively checked.
Proof: We want to solve
By the Fourier Slice Theorem for the X-Ray transform, this is solved if we take such that
But for any fixed , there are many such that , so at first glance this definition looks problematic. To make it work we must show that, for fixed , for all orthogonal to .
We will do this by showing that the directional derivatives of are zero in all directions orthogonal to . It is easier to work on than on the manifold , so we will study an extension of to and show that the -gradient of this extension is parallel to , thus the directional derivatives in all directions orthogonal to will be zero.
Defining this extension, we need to make sure that it does not vary as we scale , or the -gradient will have a component in the direction and will not be parallel to . With that in mind, define the extension as follows
using the homogeneity in the definition of . Applying the inverse Fourier Transform gives
When . For general we have
Because , and hence , is well behaved, we can differentiate this expression under the integral sign.
at points where . The general expression is quite complicated, but many of the terms are multiples of and . It is a worthwhile exercise to verify this expression.
Now we can use John’s equation to say that, when
But this is an (inverse) Fourier Transform of a function and can only vanish if the transformed function vanishes ( is injective). So
The expression on the left hand side is nothing but the coordinate of the wedge product . So this wedge product vanishes. That can only happen when the two (co)vectors are parallel, so there is a scalar valued function such that
If you aren’t comfortable with wedge products, it is just a short hand way of making the following calcuation.
which clearly implies that and are parallel.
This is exactly what we were looking for: is an extension of , and its derivatives vanish in all directions orthogonal to . So when we look at a small change in , as long as that change remains perpendicular to :
So for fixed , is locally constant on . This set is connected, so is constant on and our definition of is sound.
The rest of the proof is straightforward.
(A pdf version of these notes is available here.)