This paper is concerned with a general condition for the solvability of integral geometry problems along the plane curves of given curvatures. As two important results, the solvabilities of integral geometry problems along the family of circles with fixed radius and along the family of circles of varying radius centered on a fixed circle are given. By using some extension of the class of unknown functions, the proofs are based on the solvabilities of equivalent inverse problems for transport-like equation.
MSC: 35R30, 53C65, 65N30.
Keywords:integral geometry problem; inverse problem; Galerkin method; transport-like equation
The problems of integral geometry are to determine a function, given (weighted) integrals of this function over a family of manifolds, and there has been significant progress in the classical Radon problem when manifolds are hyperplanes and the weight function is unity, there are interesting results in the plane case when a family of curves is regular or in the case of a family of straight lines with arbitrary regular attenuation [, Chapter 7]. It is assumed that the basis of the integral geometry problems is the Radon transform . The Radon transform R integrates a function f on over hyperplanes. Let be the hyperplane perpendicular to (unit sphere) with signed distance from the origin, and the Radon transform is defined as the integral of f over , i.e.,
(see [, Chapter 2]).
The problems of integral geometry have important applications in imaging and provide the mathematical background of tomography, where the main goal is to recover the interior structure of a nontransparent object using external measurements. The object under investigation is exposed to radiation at different angles, and the radiation parameters are measured at the points of observation. The basic problem in computerized tomography is the reconstruction of a function from its line or plane integrals, and there are many applications related with computerized tomography: medical imaging, geophysics, diagnostic radiology, astronomy, seismology, radar and many other fields (see, e.g., ).
From the applied point of view, the importance of integral geometry problem over a family of straight lines in the plane is indicated in , where the problem models X-rays, and applicable to the problems of radiology and radiotherapy. Because of their many practical applications, a considerable attention has been devoted to other family of curves in the plane as well as straight lines. Invertibility of the Radon transforms on some families of curves in the plane is given with explicit inversion formulas via circular harmonic decomposition in  and for the explicit inversion formulas of the attenuated Radon transform, see, e.g., [7,8]. Note that the circle is the simplest non-trivial curve in the plane next to the straight line, and the representation of a function by its circular Radon transform also arises in applications. In , invertibility of the Radon transforms over all translations of a circle of fixed radius and circles of varying radius centered on a fixed circle is considered, where the proofs require microlocal analysis of the Radon transforms and a microlocal Holmgren theorem. In , some existence and uniqueness results on recovering a function from its circular Radon transform with partial data are presented and the relations to applications in medical imaging are described. There are several other ways related to the selection of a family of curves, such as circles of varying radius centered on a straight line or a fixed curve, circles passing through a fixed point, along paths that are not on the zero sets of harmonic polynomials, circular arcs having a chord of fixed length rotating around its middle point etc., which are meaningful in applications on thermo-acoustic and photoacoustic tomography, synthetic aperture radar, Compton scattering tomography, ultrasound tomography etc. (see, e.g., [10-12] and the references therein).
In fact, since the seminal work of Radon , the various integral geometry problems with numerous applications have been considered in several important aspects which are not mentioned here, but for a comprehensive list, see, e.g., [13-16] and the references therein. Furthermore, the problems of integral geometry and inverse problems for transport equations are interrelated and the latter are also of great importance in theory and applications; see, e.g., [17-20] and for the derivation and applications of transport equations, see, e.g., [21-23].
In this paper, a general condition for the solvability of integral geometry problems along plane curves of given curvatures is presented and its relation with some previous results is indicated. Moreover, as two important results, the solvabilities of integral geometry problems along the family of circles with fixed radius and along the family of circles of varying radius centered on a fixed circle are given. Since the curvature of a circle is defined to be the reciprocal of the radius of the circle, in the former case the curvature is a constant, while in the latter one the curvature depends on the point and the direction. To investigate the solvability of the integral geometry problem (IGP) given in Section 2.1, which is overdetermined since the underlying operator of the IGP is compact and its inverse operator is unbounded (see Section 2.2), it is reduced to an equivalent overdetermined inverse problem for a transport-like equation, and then, with the use of a similar method which was proposed in  (see also [, Chapter 1]), on using some extension of the class of unknown functions, this inverse problem is replaced by a determined one. Thus, the solvability of IGP is proved via the solvability of an inverse problem for a transport-like equation. The above mentioned method on the solvability was also previously utilized in [25-27] for IGP along some family of plane curves of given curvatures and straight lines, in  for IGP along geodesics and in  for IGP along the family of curves whose curvatures are given by the Christoffel symbols. Here, the presented general condition for the solvability covers those of [25-28], in the manner indicated by Remark 1 in Section 3.1. Moreover, these previous solvability conditions do not hold for IGP along the above given two families of circles, and this is the main importance and motivation of this study.
In Section 2, IGP and its reduction to the equivalent inverse problem for a transport-like equation, the method to overcome the difficulty on investigating the solvability arising from overdeterminacy of these problems are presented and some definitions and notations which will be used throughout the paper are introduced. Section 3 is devoted to the statements of main results, and finally in Section 4, proofs of the results are given.
The investigation of approximate solutions of the concerned integral geometry problems is beyond the scope of this paper, but similar procedures as in  can be carried out by using the Galerkin method or the finite difference method.
2 Statement of the problem and overdeterminacy
2.1 Statement of the problem
Let D be a bounded domain in . It is assumed that in D, a family of regular curves is given by curvature which is the curvature of the curve passing from the point in the direction , and there exists a unique sufficiently smooth curve of this family which is passing from any point in the arbitrary direction ν, with the endpoints on the boundary of D. Suppose that the lengths of these curves in D are bounded above by the same constant. Let us denote the family of these curves by . IGP is stated below.
IGPDetermine a function in the domainDfrom the integrals ofλalong the curves of a given family of curves .
Suppose that vanishes outside D, and let us introduce an auxiliary function
where is a part of the curve that belongs to , with one end of it being the point x and the other one on ∂D, and dσ is the arc length element along .
Investigating the uniqueness of a solution of a problem of integral geometry by reducing it to the equivalent inverse problem for a differential equation was first carried out in . Similar reduction is demonstrated for IGP formulated below.
Differentiating (1) in the direction ν at x, we obtain the following transport-like equation:
From (1), u is 2π-periodic with respect to φ, and since the integrals of λ along the curves of are known, u is known on , i.e.,
Problem 1 Determine a pair of functions from the transport-like equation (2) provided that the function K is known and u satisfies conditions (3).
Generally, in the theory of integral geometry, reconstruction of a function of n variables from a function of variables is said to be an overdetermined problem of integral geometry (see, e.g., [, Chapter 5]), and in the theory of inverse problems, overdetermination usually means that the number m of independent variables in the data exceeds the number n of independent variables in the unknown target function (see, e.g., [, Section 1.3]). However, since the data given on the two-dimensional surface and is a function of two variables, these are not the cases for IGP or Problem 1. Here, since the operator given in (1) (for , where ) is compact, its inverse is unbounded, and therefore it is not possible to prove a general existence result. So, IGP and Problem 1 are called overdetermined in this sense. Hence, because of the overdeterminacy, the initial data for these problems should not be arbitrary and satisfy some ‘solvability conditions’ (see [, p.4] and [, p.18, Theorem 1.4]) which are difficult to establish. It should be noted that the set of functions for which IGP is solvable is not everywhere dense in any of the spaces , and . Moreover, the data in problems of integral geometry are of quasianalytic character, i.e., their values specified in a domain of the Lebesgue measure can be as small as desired, determine their values in an essentially larger domain (see [, Chapter 6, Section 17] and [, Chapter 6, Section 1]). In particular, this implies that it is impossible to avoid overdeterminacy of the problem by specifying the data on a part of the boundary rather than on the whole boundary. Even if it were possible to find the solvability conditions for the mentioned overdetermined problems, since the real data usually have some errors in practice, and thus fall out of the data class for which the existence of a solution is established, it appears that these conditions would not always be satisfactory in applications. Therefore, to prove the existence results, such special conditions on the data have to be posed.
Let us propose the procedure for establishing the solvability of IGP. Assume that the unknown function λ in IGP depends not only upon the space variables x, but also upon the direction φ in some special manner, i.e., consider , where this dependence upon φ is impossible to be arbitrary, for in the opposite case the problem would be underdetermined and the examples on the nonuniqueness of a solution can be easily constructed. Herein the special dependence of upon the direction means that satisfies a certain differential equation ( , where the expression of is given in Section 2.3) with the following properties:
(1) The IGP or Problem 1 with the function becomes a determined one.
(2) The sufficiently smooth functions λ depending only on x satisfy this equation.
Suppose that a differential equation for satisfying properties (1) and (2) has been found and that a priori the function , which represents the exact data of IGP related to a function λ depending only on x, is known. Then, utilizing , a solution to IGP can be constructed. By uniqueness of a solution, and coincide. At the same time, knowing the approximate data with , an approximate solution can be constructed such that . Recall that if λ depends only on x and does not satisfy the ‘solvability conditions’, the solution depending only x does not exist. Here the data are specified on and is independent of and . In other words, a regularizing procedure is constructed for the IGP.
In general, the equation with the properties (1) and (2) for the same problem is not uniquely defined. Hence, the class of unknown functions λ extends so that IGP for this class becomes a determined problem and all sufficiently smooth functions in x belong to it. On using some extension of the class of functions λ, the overdetermined Problem 1 is replaced by a determined one (Problem 2 in Section 3.1).
The above method of solvability of the IGP or Problem 1 leads to the Dirichlet-type problem with conditions (3) for the third-order equation of the form . In investigating and proving the solvability of IGP over any regular family of curves with curvature K, since the quadratic form in (14) is required to be positive definite, the construction of is important and to be able to this, in Theorem 1 condition (6) is given.
2.3 Definitions and notations
In this section, some notations are given based on . Let with the boundary , and be the closure of Ω. By we denote a scalar product of functions u and v in and by the set of all functions defined in Ω which have continuous partial derivatives of order up to all , whose supports are compact subsets of Ω. For a differential expression A, by we denote the conjugate of A in the sense of Lagrange. For ,
Let denote the set of real-valued functions that are 2π-periodic with respect to φ in the domain Ω, i.e., , where are nonnegative integers such that .
The proof of Theorem 1 involves energy-like estimates and the Galerkin method (see, e.g., [, Chapter 5, Section 2.3], [, Chapter 7]), and therefore some class of functions are introduced below. In , let us introduce the scalar product
Let , and and be the completions of with respect to the norms and ( ). Let us take a set which is complete and orthonormal in , then we may assume that the linear span of this set is everywhere dense in . Since is separable, there exists a countable set which is everywhere dense in this space and this set up can be extended to a set which is everywhere dense in . Orthonormalizing the latter in , we obtain . We denote by the orthogonal projector of onto which is the linear span of .
where and it can be easily verified that
The existence of the function g in the expression of leads to a generalization on conditions for the solvability of integral geometry problems. In Theorem 1, it is shown that if there exists a function g satisfying condition (6) which depends on the curvature K and the domain D, the solvability holds.
Let and be the set of all functions such that for any there exists such that in the generalized functions sense, i.e., holds for every . Take a subset such that for any there exists a sequence such that weakly in and as . If we denote the closure of with respect to the norm by , then we have and it can be shown that the inclusions hold.
3 Statements of results
3.1 Solvability of IGP along plane curves
Since Problem 1 is overdetermined, as indicated in Section 2.2, we consider the following determined problem.
Problem 2 Determine a pair of functions defined in Ω that satisfies
provided that , u is 2π-periodic with respect to φ, and K are known.
In (4), it is assumed that the unknown function λ depends also on φ and the condition holds in the generalized functions sense.
If and , then there exists a function (see [, p.130, Theorem 2]) such that , and we can consider the new unknown function . Hence from (4) we obtain the equation , where and . Let us denote again by u for simplicity, then Problem 2 can be reduced to Problem 3 given below (see [, p.20]) and the solvability of the former follows from that of the latter and does not depend on the choice of G.
Problem 3 Determine a pair of functions defined in Ω that satisfies
provided that , u is 2π-periodic with respect to φ, , K and F are known.
The existence, uniqueness and the stability of the solution of Problem 3 are given by the following theorem.
Theorem 1If and a function exists such that
for all , then Problem 3 has a unique solution such that , and
holds, where depends onKand the Lebesgue measure ofD.
Remark 1 In fact, without being aware of (6), the function g, with the appropriate choices of it, was used previously in [25,26,28]. The convenience of (6) with those of previous solvability results is indicated below.
(i) In , when the curvature K is sufficiently smooth and 2π-periodic, the condition for the solvability is for all , where this condition holds for in (6).
(ii) In , when , where and are sufficiently smooth functions, the condition for the solvability is for all , which holds for in (6).
(iii) In , when , where and are sufficiently smooth 2π-periodic functions, the condition for the solvability is for all , and this condition holds for in (6).
Note that in none of the above cases, the solvability conditions hold for the IGP along the family of circles of varying radius centered on a fixed circle and the family of curves of constant curvatures, i.e., the family of circles with fixed radius where the curvature is a nonzero constant or the family of straight lines where the curvature is zero. The former cases are investigated in Section 3.2 below and the latter case is considered in , where the term in the expression of Lu will not be present since and the proof of solvability is given for . In fact, the strict inequality in (6) can be written as a non-strict inequality, with the equality only for .
3.2 Solvability of IGP along the family of circles
Since Theorem 1 was given for IGP along a regular family of plane curves for the general case, in this section the given results on the solvabilities of IGP along the family of circles depend on finding an appropriate function g satisfying (6) for the given curvature K and the domain D under the assumptions of Theorem 1.
3.2.1 Solvability of IGP along the family of circles with fixed radius
Let and . Let us take the family of curves in IGP as the family of circular arcs (the segments of circles inside D with the endpoints on ∂D) with fixed radius r, passing from the point in the direction and denote this family of circles by . Since the curvature of a circle is defined to be the reciprocal of the radius of this circle, the curvature of the elements of the family becomes .
The solvability of IGP along the curves of follows from the following lemma.
Lemma 1Let us definegon as
If , then (6) holds for .
Remark 2 The above choice of g is not unique, and since the curvature is constant, and hence , condition (6) reduces to
for all , which depends on r and the domain D.
3.2.2 Solvability of IGP along the family of circles of varying radius centered on a fixed circle
Let and . Let be a fixed disk of radius and, without any loss of generality, centered at the origin such that , and take the family of curves in IGP as the family of circular arcs (the segments of circles inside D with the endpoints on ∂D) of varying radius centered on , passing from the point in the direction and denote this family of circular arcs by . It can be shown that the radius of the circles of the family is defined by the function
and hence the curvature of the elements of the family is defined by
on . Note that since , we have and on .
The solvability of IGP along the curves of follows from Lemma 2 given below.
Lemma 2If the functiongis defined on as
then (6) holds forKdefined in (11).
4 Proof of results
Proof of Theorem 1 First, we will prove the uniqueness of the solution of Problem 3 under the assumptions of the theorem. To this end, it is sufficient to show that the corresponding homogeneous problem has only a trivial solution. Then, by taking into account that and , from (5) we obtain .
Since , there exists a sequence such that weakly in and as . Now we want to decompose the product into the sum of a positive definite quadratic form and a divergence form. For this purpose, we have the following identities:
If (6) holds, then the quadratic form in , , is positive definite, where
Indeed, we can estimate the terms and as follows:
where , and we obtain
Moreover, whenever (6) holds, for sufficiently close value of ε to 1, there exists an such that in Ω, and we obtain
Since the domain D is bounded, on and is positive definite, we have
where and depends on the Lebesgue measure of D and does not depend on k.
Thus, since , K and g are 2π-periodic with respect to φ, after integrating (13) over Ω, the divergent terms disappear and we obtain
Since , from (15) we have
which implies that , and since , from (5) we get . So, the uniqueness part of the proof is completed.
Now we will prove that there exists a solution of Problem 3 in by means of the following auxiliary problem.
Determineudefined in Ω that satisfies
The solution u of problem (17)-(18) will be approximated by
construction of which is based on finding the vector from the system of linear algebraic equations
where the system of functions is taken as indicated in Section 2.3 and on .
We must show that the solution of system (19) exists and is unique for any . To demonstrate this, let us assume that the homogeneous version of (19), i.e., the system
has a nonzero solution . Substituting for , multiplying the jth equation of the above homogenous system by and summing with respect to j from 1 to N, we obtain
where . So, from (15) and (20) we obtain
and since the quadratic form defined in (14) is positive definite and on , we have in Ω. But is linearly independent and this implies that , , which contradicts with the assumption . So, it is shown that system (19) has a unique solution for any .
Now we estimate the solution of system (19) in terms of F. For this purpose, we multiply both sides of the jth equation of (19) by and sum the obtained equations with respect to j from 1 to N to obtain
Since , applying integration by parts, the right-hand side of (21) can be estimated as
for , and from (15) and (21) we have
It can be verified that for sufficiently large , from (22) we obtain
where the constant C is independent of N. This implies that is bounded in and , and since and are Hilbert spaces, it is weakly compact in and . Therefore, there exists a subsequence, which we again denote by , such that weakly in and as and
holds. Since and weakly in , we have . From (23) we have also that , and are bounded and there exists a subsequence of , which is again denoted by , such that , and converge weakly in to , and , respectively. Taking into account that , and applying integration by parts in (19), for , we obtain
Since the linear span of is dense on the space , passing to the limit as , we get
for every . If we set , since , from (24) we have in the generalized functions sense. Moreover,
holds and, by using , it can be seen that (7) holds. In the above expressions, by C we denote generic constants which depend only on the given functions and Lebesgue measure of the domain D.
Now it remains to show that . Since and , for , from (24) we have
for any , which implies that in the generalized functions sense, i.e., .
Moreover, from (19) we have , where is the orthogonal projector. Since the system is orthogonal and complete in , strongly, i.e., strongly in as . Then, since weakly in as , we have as . On the other hand, since the projection operator is self-adjoint in ,
thus as , and so . The proof of Theorem 1 is complete. □
Proof of Lemma 1 Since is constant, and hence
to prove that condition (6) holds, we only need to show that g satisfies
By taking into account that and , for the function
given in (8) defined on , we obtain
and the proof is complete. □
Proof of Lemma 2 For
Since and , for
given in (12) defined on , we have
Hence, to prove that condition (6) holds, we only need to show that
holds. If we take into account again that and , then we obtain
The proof is complete. □
The author declares that he has no competing interests.
Cormack, AM: Representation of a function by its line integrals, with some radiological applications. J. Appl. Phys.. 34, 2722–2727 (1963). Publisher Full Text
Cormack, AM: The Radon transform on a family of curves in the plane. Proc. Am. Math. Soc.. 83, 325–330 (1981). Publisher Full Text
Natterer, F: Inversion of attenuated Radon transform. Inverse Probl.. 17, 113–119 (2001). Publisher Full Text
Novikov, RG: An inversion formula for the attenuated X-ray transformation. Ark. Mat.. 40(1), 145–167 (2002). Publisher Full Text
Rigaud, G, Nguyen, MK, Louis, AK: Novel numerical inversions of two circular-arc Radon transforms in Compton scattering tomography. Inverse Probl. Sci. Eng.. 20(6), 809–839 (2012). Publisher Full Text
Arridge, SR: Optical tomography in medical imaging. Inverse Probl.. 15, R41–R93 (1999). Publisher Full Text
Tamasan, A: An inverse boundary value problem in two-dimensional transport. Inverse Probl.. 18, 209–219 (2002). Publisher Full Text
Klibanov, MV, Yamamoto, M: Exact controllability for the time dependent transport equation. SIAM J. Control Optim.. 46(6), 2071–2095 (2007). Publisher Full Text
Amirov, A, Yildiz, M, Ustaoglu, Z: Solvability of a problem of integral geometry via an inverse problem for a transport-like equation and a numerical method. Inverse Probl.. 25, Article ID 095002 (2009)
Ustaoglu, Z, Heydarov, B, Amirov, S: On the solvability and approximate solution of a two dimensional coefficient inverse problem for a transport-like equation. Inverse Probl.. 26, Article ID 115019 (2010)
Amirov, A, Ustaoglu, Z, Heydarov, B: Solvability of a two dimensional coefficient inverse problem for transport equation and a numerical method. Transp. Theory Stat. Phys.. 40(1), 1–22 (2011). Publisher Full Text
Golgeleyen, I: An integral geometry problem along geodesics and a computational approach. An. Univ. “Ovidius” Constanţa, Ser. Mat.. 18(2), 91–112 (2010). PubMed Abstract