# Inversion Formula for a Radon-Type Transform Arising in Photoacoustic Tomography with Circular Integrating Detectors

Authors: Sunghwan Moon
Journal: Advances in Mathematical Physics
 (2018)

---

## Abstract

This paper is devoted to a Radon-type transform arising in a version of Photoacoustic Tomography that uses integrating circular detectors. The Radon-type transform that arises can be decomposed into the known Radon-type transforms: the spherical Radon transform and the sectional Radon transform. An inversion formula is obtained by combining existing inversion formulas for the above two Radon-type transforms.


---

## Body

## 1. Introduction

Hybrid biomedical imaging is an emerging technology using the combination of different physical signals in order to utilize their advantages to enhance image. Among hybrid biomedical imaging, thermoacoustic tomography (TAT) and its sibling Photoacoustic Tomography (PAT) are the most successful examples. Many computed tomography researchers have been considering the mathematical modeling of PAT/TAT and many breakthroughs have been achieved. Also, mathematical problems arising in PAT/TAT also relate to sonar and radar research (e.g., [1, 2]).TAT/PAT utilizes ultrasound and optical or radiofrequency electromagnetic waves. While pure ultrasound imaging typically offers high resolution images with low contrast, optical or radiofrequency electromagnetic imaging offers an enormous contrast between cancerous and healthy tissues and low resolution. After the photoacoustic effect discovered by A.G. Bell [3], one can take advantage of the strengths of the pure optical and ultrasound imaging and form the basis of TAT/PAT.In TAT, a very short radiofrequency (RF) pulse is sent through a biological object, more or less uniformly irradiated. RF energy, absorbed throughout the object, causes thermoelastic expansion of the tissue and emergence and propagation of a pressure wave (an ultrasound signal). This pressure wave can be measured by transducers placed outside the object. The spatial distribution of the absorption of RF energy will have a very good contrast. This is due to the fact that cancerous cells absorb several times more energy than healthy tissue in the RF range. Since the intensity of the source of the acoustic wave describes the spatial distribution of the absorption of RF energy, one of the classical mathematics problems of TAT is the recovery of this initial acoustic pressure field from the ultrasound measurements made outside the object.PAT differs from TAT only in the way the thermoacoustic signal is triggered. PAT uses a laser pulse rather than an RF. The rest of the procedure is the same as that of TAT. From the mathematical view, there is no difference between TAT and PAT [4]. Hence, from now on, we will just mention PAT.As the classic detector, small piezoelectric transducers were used in PAT. Since these transducers mimic point-like measurements, reconstruction algorithms produce images with a spatial resolution limited by the size of the transducers. Another drawback of such detectors is the difficulty in manufacturing small transducers with high sensitivity. To overcome these weak points, various other types of acoustic detectors have been introduced, e.g., linear, planar, and cylindrical detectors [5–9]. These detectors are modeled as measuring the integrals of the pressure over the shape of the detector. Linear, planar, and cylindrical detectors do not provide a simple and compact experimental buildup (because their ideal size is infinite).Zangerl et al. started to consider PAT with circular integrating detectors [10–12] (which can be realized with line-profile transducer) [13]. Some works [11, 12, 14] have considered PAT with circular integrating detectors on cylindrical geometry: a stack of parallel circular detectors whose centers lie on a cylinder. The cases where the centers of the detector circles are located on a plane or a sphere have been discussed in [14]. In [10–12] they showed that the data from PAT with circular detectors was a solution of a certain initial value problem and converted this problem into a circular Radon transform using this fact. In [10] the spherical geometry, which can be useful in investigating small object, was discussed. Here we also study PAT with the circular integrating detectors on a spherical geometry which is already discussed in [10]. Our approach is to define a new Radon-type transform arising in this version of PAT, and we show that this transform can be decomposed into the spherical Radon transform and the sectional Radon transform introduced by Rubin [15, 16], both of which are already studied.This paper is organized as follows. Section2 is devoted to the construction of a Radon-type transform arising in PAT with circular detectors. In Section 3, we show that the Radon-type transform is the composition of the sectional Radon transform and the spherical Radon transform with centers on the boundary of a sphere, and we present the inversion formula using this fact.


## 2. Photoacoustic Tomography with Circular Integrating Detectors

In PAT, the acoustic pressurep(x,t) satisfies the following initial value problem:(1)∂t2px,t=Δxpx,tx,t=x1,x2,x3,t∈R3×0,∞,px,0=fxx∈R3,∂tpx,0=0x∈R3.(We assume that the sound speed everywhere including the interior of the object is equal to one.) The goal of PAT is to recover the initial pressure f from measurements of p outside the support of f.Throughout this article, it is assumed that the initial pressure fieldf is smooth and supported in the unit ball B3. We also assume that the acoustic signals are measured by a set of circular detectors on the sphere in the plane perpendicular to (θ,0)∈S2 with distance s∈[-1,1] from the origin, where θ∈S1 (see Figure 1).Figure 1
The circle on the sphere in the plane perpendicular to(θ,0)∈S2 with distance s from the origin.Basically, the spherical geometry could be very useful for the investigation of small objects. Here we studied the spherical geometry Zangerl et al. suggested [10].The measured dataP(θ,s,t) for (θ,s,t)∈S1×[-1,1]×[0,∞) can be written as (2)Pθ,s,t=12π∫θ¯·α=s,α∈S2pα,tdSα,where θ¯=(θ,0)∈S2 and dS (everywhere in the text dS denotes the measure either on the circle or on the sphere; the dimension will be clear from the context) is the measure on the unit circle S1. The shape of detectors is the circle which lies in a plane perpendicular to θ¯ with distance s from the origin (see Figure 1). Also, it is a well-known fact that (3)px,t=∂t14πt∫∂B3x,tfβdSβis a solution of PDE (1) (see [17]). Here Bt3(x)=B3(x,t) is a ball in R3 centered at x∈R3 with radius t and dS is the standard measure on the sphere. Hence P(θ,s,t) becomes(4)Pθ,s,t=12π∫θ¯·α=s,α∈S2∂t14πt∫∂B3α,tfβdSβdSα=18π2∂tt∫θ¯·α=s,α∈S2∫S2fα+tβdSβdSα.Let us define, for (θ,s,t)∈S1×[-1,1]×[0,∞), (5)RPfθ,s,t=∫θ¯·α=s,α∈S2∫S2fα+tβdSβdSα.Then, (4) reads as(6)Pθ,s,t=18π2∂ttRPfθ,s,t.If f is odd in x3, i.e., f(x1,x2,x3)=-f(x1,x2,-x3), then RPf is equal to zero. We thus assume that f is even in x3, i.e., f(x1,x2,x3)=f(x1,x2,-x3). Indeed, any function compactly supported in the upper hemisphere of B3 can be extended to an even function with support in B3. This means that, for any object located on only the hemisphere of B3, we can reconstruct its image through even extension. We will demonstrate that RP is the composition of the spherical Radon transform centered at S2 and the sectional Radon transform. This will allow us to recover f from P.


## 3. Reconstruction

We first show thatRP can be decomposed into the spherical Radon transform centered at S2 and the sectional Radon transform.Definition 1.
(1)
The sectional Radon transformQ maps a function ϕ on S1×[-1,1]×[0,∞) with ϕ(α1,α2,α3,t)=ϕ(α1,α2,-α3,t) into (7)Qϕθ,s,t=∫θ¯·α=sϕα,tdSα,forθ,s,t∈S1×-1,1×0,∞. (2)
The spherical Radon transformRS maps a locally integrable function f on R3 into (8)RSfα,t=∫S2fα+tβdSβ,forα,t∈S2×0,∞.Now it is easily seen that the following representation ofRPf holds.Theorem 2.
For anyf∈Cc∞(R3), we have RPf(θ,s,t)=Q(RSf)(θ,s,t).The spherical Radon transform is a classic object in computed tomography. Many inversion formulas for this transform have been derived in [18–20]. The sectional Radon transform was introduced in [15, 16, 21, 22] and the relation between the sectional Radon transform and the regular Radon transform has been derived in [21, 22].Combining the inversion formula for the spherical Radon transform derived in [18] and the relation discussed in [22], we present the inversion formula for RPf as follows.Theorem 3.
Letf∈C∞(R3) have compact support in B3 and be even in x3. Then we have (9)fx=-164π4∫S2∫S1∫Rα3∂tt∂tt∂sRPfθ,s,tt=α-xθ·α1,α2-sdsdSθα-xdSα.Proof.
From [21, 22], we know that, for θ=(θ1,θ2)∈S1, (10)Qϕθ,s,t=∫α1θ1+α2θ2=sϕα,tdSα=RΦα1,α2,tθ,s,t,where Φ(α1,α2,t)=2ϕ(α,t)/α3 and R is the 2D regular Radon transform of Φ with respect to the first two variables; that is, (11)RΦα1,α2,tω,s,t=∫RΦsω+νω⊥,tdν,forω,s∈S1×R.Applying the inversion formula for the 2D regular Radon transform, we have(12)ϕα,t=2-1α3Φα1,α2,t=α38π2∫S1∫R∂sQϕω,s,tθ·α1,α2-sdsdSθ.Also, it is well known (see, e.g., [18]) that(13)fx=-18π2∫S2∂tt∂ttRSfα,ttt=α-xdSα.Combining (12) and (13) completes the proof.Combining Theorem3 and (6) shows how to recover f from P.Corollary 4.
Letf∈C∞(R3) have compact support in B3 and be even in x3. Then we have (14)fx=-164π4∫S2∫S1∫Rα3∂tt∂sPθ,s,tt=α-xθ·α1,α2-sdsdSθα-xdSα.


## 4. Conclusion

Here we studied a Radon-type transform arising in PAT with circular detectors. We show that this transform can be decomposed into the spherical Radon transformRS centered at S2 and the sectional Radon transform Q. To recover f from RPf, we reconstruct RSf first. The formula to recover RSf from RPf is essentially (12) which is equivalent to the reconstruction of 2D Radon transform. This formula is well studied, and Matlab has “iradon” as a built-in function. The second part about the reconstruction f of RSf is also well studied in [18–20]. Some works present numerical simulations for the inversion formula. Therefore, the numerical simulations are not within the scope of this work.



---

*Source: 1727582-2018-07-03.xml*