A new Compton scattered tomography modality and its

International Refereed Journal of Engineering and Science (IRJES)
ISSN (Online) 2319-183X, (Print) 2319-1821
Volume 4, Issue 1 (January 2015), PP.13-19
A new Compton scattered tomography modality and its
application to material non-destructive evaluation
E.R. Lopez-Villaverde 1, S. -N. Vu1, M. K. Nguyen 1, A. Noumowe2
Equipes Traitement de l'Information et Systèmes (ETIS) / ENSA/Université de Cergy-Pontoise/UMR CNRS
8051/, F-95302 Cergy-Pontoise, France)
Laboratoire de Mécanique et Matériaux du Génie Civil(L2MGC) EA 4114/ Université de Cergy-Pontoise,
F-95000 Cergy-Pontoise, France)
Abstract: - Imaging modalities exploiting the use of Compton scattering are currently under active
investigation. However, despite many innovative contributions, this topic still poses a formidable mathematical
and technical challenge. Due to the very particular nature of the Compton effect, the main problem consists of
obtaining the reconstruction of the object electron density. Investigations on Compton scatter imaging for
biological tissues, organs and the like have been performed and studied widely over the years. However in
material sciences, in particular in non-destructive evaluation and control, this type of imaging procedure is just
at its beginning. In this paper, we present a new scanning process which collects scattered radiation to
reconstruct the internal electronic distribution of industrial materials. As an illustration, we shall look at one of
the most widely used construction material: concrete and its variants in civil engineering. The Compton
scattered radiation approach is particularly efficient in imaging steel frame and voids imbedded in bulk concrete
We present numerical simulation results to demonstrate the viability and performances of this imaging modality.
Keywords: - Compton scattering , Gamma-ray imaging , Non-destructive testing/evaluation (NDT/NDE),
Concrete: structure and defects, Radon transform
In medicine, hidden structures of objects are routinely revealed by ultrasound imaging and radiation
tomography, which comprises three modalities: transmission, emission and scattering. In material sciences, in
addition of these, many more technological imaging processes are available, such as magnetic flux leakage,
eddy current, infrared and thermal testing, Radar technology, etc.). However these methods although less costly
and easy to operate lack accuracy and cannot deliver key physical quantities necessary for material evaluation.
In this respect, radiation imaging, despite stringent safety restrictions, is more reliable in bringing accurate data
on relevant physical quantities. This is the reason why many medical tomographic imaging modalities have been
in transferred to non-destructive evaluation systematically.
Although transmission radiation imaging (tomography) has been in operation for a long time, it has
been realized that there is a need to determine other physical parameters than the ones deduced from the
absorption of radiation in matter. Moreover working with the propagation of radiation through matter,
conventional transmission imaging cannot handle some large objects or objects that are improper to be inserted
in an industrial tomograph. Transmission imaging provides line-integrated information along the radiation
propagation path, which masks the position of an anomaly present along the line. Therefore, it is difficult,
among other things, to determine the position of this anomaly from the transmission data.
Thus the need to remedy to that situation has led workers to propose the use of radiation scattering
imaging. By recording the amount of deflected radiation from the object, any anomaly will not be hidden in the
integrated data. Scattering also eliminates the need of accessing the target object from two opposing sides.
Therefore, since both source and detector can be placed on one side of the object, examination of extended
structures becomes possible. Pioneering works, done in recent years, have confirmed these advantages in nondestructive evaluation of industrial objects and materials.
One of the most widely used industrial material is concrete. It enters a large variety of objects in civil
engineering. Therefore the monitoring of its status in time becomes of primary concern for safety and
maintenance. What is needed is a non-destructive evaluation (NDE) or a non-destructive testing (NDT) of its
status which could provide images of its interior, so as to see and determine anomalies ad their time evolution
13 | Page
A new Compton scattered tomography modality and its application to material non-destructive…
[1,2]. The most common defects are internal cracks, voids, shallow delamination, honeycombing and surface
opening cracks. They are the factors that affect most the durability of concrete. In the case of reinforced
concrete structures, other defects such as thickness variations, deteriorated zones, corroded reinforcing bars, see
e.g. [3,4] and moisture may occur. Imaging of concrete structures is very important because it reveals aging and
deterioration of infrastructures, as pointed out by [5]. Hence their detection is of utmost importance in order to
prevent material failure which leads to accidents. Common detection techniques using ultra-sounds or
radiography generally provide only qualitative information. Infrared thermography, radar and acoustic imaging
are also applied but with coarser resolution. Radiography can be used to locate internal defects and reinforced
re-bars position in concrete and general condition through differences in radiation intensity passing through the
structure that is captured on photographic film placed on the opposite side of the structure from the source, see
[6,7,8]. Computerized Tomography (CT) is a method that permits the development of three-dimensional
radiograph views of an object. The degree of damage in concrete structures can be assessed using high
resolution three-dimensional images obtained from this CT technique [9]. However the performances of CT
have limitations as pointed out earlier. They may be overcome by scattered radiation imaging which also brings
about useful complementary information.
Since gamma rays from radioactive isotopes have sharply defined energies, the analysis of the detected scattered
signal is easier. Moreover, gamma-ray isotopic sources are readily portable, self-contained an usable in hostile
environments. The most common modality in scattered radiation imaging proposed so far is a step by step
procedure. It consists in determining the electron density at a site by measuring the deflected amount of
radiation by Compton effect by a fixed scattering angle, given an incoming calibrated amount of incident
radiation. This procedure has been tested and shown to be viable over the years. However it is time consuming
and not very practical, although it enables the detection of local defects and the discrimination between
materials of different density such as voids and steel in concrete. We propose an alternative modality whereby
data acquisition is continuous: the detector is set to absorb only scattered radiation of given energy and the pair
source-detector, set at a constant distance from each other, rotate around its middle point in a plane. Thus one
gets an integrated data along various circular arcs of the electron density. The reconstruction of this density is
shown to be possible thanks to an analytic inversion method, derived recently from that of A M Cormack for the
transmission scanning modality of computerized tomography (CT) [24].
In this work, we describe this Compton scatter imaging modality. This gamma scattering tomography
resolves some of the problems encountered in radiography and CT. A detailed comparison of these two
techniques can be found in [9]. Since the sought quantity is the electron density of matter, Compton scatter
imaging is particularly efficient in the detection and imaging of voids and steel rabars in reinforced concrete
because of large discontinuities in traversing these defects. Moreover in concrete bulk a change in water content,
which is due to higher porosity, may be also seen in variations of measured electron density. This may indicate a
change in mechanical properties of concrete via a change of its compression Young modulus. Thus the
knowledge of an electron density map for concrete may be extremely relevant for nondestructive evaluation.
This quantity is becoming the focus of attention in recent years [10,11].
Up to now the most currently used method is the point-wise scanning, by which information can be
obtained by focusing the field of view of the source and detector around a scattering volume. Then the measured
scattered signal allows the determination of the density of the material. The disadvantages are relatively poor
resolution, slow scanning speed and difficulties associated with numerical reconstruction of object’s internal
characteristics. The scattered signal is also affected by the attenuation of both incident and scattered photons
which makes the reconstruction difficult. Collimation of the initial and scattered radiation allows the detection
of relatively few scattered photons, and this is the reason why the statistics relative to this process are very poor,
see [12,13]. Finally by scanning through a plane of interest within an object using raster motion, one may obtain
density distribution in a whole plane.
Here we propose a different imaging approach by scattering of ionizing radiation. Previous work on the
use of scattered radiation has originated from medical imaging but this idea has made quick headway to NDT
and NDE through many proposals such as [14]. To generate scattering tomography at high resolution one
follows the idea of Kondic by installing wide-angle collimated source and detector, see [15]. This allows
counting almost all of the scattered photons leaving the object in the direction of the detector. Moreover to get
the necessary complete data we must allow for the rotation of the pair source-detector around the object. In fact
this was already foreseen long ago by [16] and [17]. But at this time no analytic inversion formula was available
and numerical inversion of data has led to poor reconstruction results. Our proposal differs from that of [18,19],
as well as [20,21,22,23], since in their scanning procedure, the detector runs along a line passing through the
source which is not very convenient for practical operations except in the case of large flat structures. Of course
their image reconstruction is based on a different approach. The recent derivation of the reconstruction formula
14 | Page
A new Compton scattered tomography modality and its application to material non-destructive…
for the electron density has prompted applications in medical imaging [26]. We propose now an application in
non-destructive evaluation of concrete as an alternative to existing NDE methods. In section 2 we describe the
new Compton scatter tomography modality, with image formation as well as image reconstruction. Then in the
next section 3 we discuss the corresponding simulation results. Conclusion and possible future perspectives are
given in the last section.
The principle of this imaging modality is simple: it makes use of the Compton effect or scattering of
ionizing radiation by an electron. Fig. 1 briefly recalls the mechanism of Compton scattering. An incident X or
gamma radiation of energy E0 hits an electron at rest and is deflected from the initial direction of motion by an
angle ω.
The scattered radiation (or photon) carries the energy E0 given by the famous Compton relation
where mc2 is the rest energy of the electron or 0,511 MeV.
Image formation and circular arc Radon transform.
Fig. 2 shows how data acquisition is performed in this new Compton scattering tomography modality. Test
objects under study are put inside a circular frame and are illuminated by an isotropic radiation source S of
definite energy E0 situated on the circular frame. A radiation detector at site D, diametrically opposite to S via a
multichannel analyzer registers scattered radiation at scattered energy E. From the known kinematics of
Compton scattering, the radiation amount collected by a unit surface of detector per unit time is proportional to
the integral of the object (concrete) electron density along a circular arc joining the point source site S to the
detector site D, the circular arc subtending the angle (π-ω), for details see [24]. We denote the sought electron
density by f(r,θ) which describes the inner state of the object. The line SD has a length 2p and rotates around its
center O during scanning, the rotation state of the apparatus being conveniently labeled by the angle φ. The
detected radiation flux density at D is Cf(φ,ω) and given by the integral of f(r,θ) on the circular arc C(φ,ω) up to
a factor representing the Compton scattering phenomena (this factor is absorbed in the definition of f(r,θ) for
where ds is the arc integration element, computed from the arc equation
with τ = cotan ω and γ = π – ω. Note that for ω = 0, the circular arc reduces to a line, and the detector receives
non scattered radiation, which is to be excluded
15 | Page
A new Compton scattered tomography modality and its application to material non-destructive…
The set of all integrals Cf(φ,ω) represents what is called the circular arc Radon transform of the object electron
density f(r,θ) . It is a function of two angular variables (φ,ω)
Remark The electronic density of a material is in fact given by the formula
where ρ is the mass density , NA=6.022x1023 mol-1 is the Avogadro number, A the mass number and Z the
atomic number.
Reconstruction of the electronic density
The inverse formula of equation (4) has been worked out in [24]. Hereafter we give a short account of its
derivation. Let us introduce the angular Fourier components of f(r,θ)
The angular Fourier components of Cf(φ,ω) are similarly defined. Then equation (4) is transformed into the
following integral equation for the Fourier angular components Cfl(τ)) and fl(r)
where τ = cotan ω. To extract the sought components fl(r) we use the method of [25] and get
where q= 1/τ. Another change of variable in equation (8) leads to the final result
Finally f(r,θ) is reconstructed by summing over its angular components fl(r) according to equation (6).
16 | Page
A new Compton scattered tomography modality and its application to material non-destructive…
Concrete is an interesting test material for this new investigation method. Let us recall how it behaves
under standard non-destructive evaluation using standard radiation attenuation method. Here the mathematical
tool is of course the classical Radon transform. What is reconstructed is the attenuation map of the material. By
using a careful pixel thresholding of the attenuation image, it is possible to deduce the porosity of the concrete
This would yield, through the Féret formula, its compression limit and elasticity modulus. So long as
the found numerical values are acceptable, the state of the investigated material can be declared safe for further
use. It is expected that the reconstructed electron density would lead to similar conclusions and thereby
complement the information obtained from the attenuation map.
Simulations are performed using a synthetic image of the reinforced concrete with iron grid. The
original concrete electron density is given in Fig. 3. For simplicity, the pixel intensity is normalized. It is to be
noted that iron grid holds the highest electron density and the lowest electron density is located in the empty
spaces. The acquired data is given in Fig. 4 in the coordinate system (φ,ω). Using equation (9) and the previous
data the reconstructed concrete electron density is given in Fig. 5. This reconstruction procedure is recent and
need to be perfected. The results are however convincing. The Mean Quadratic Error is 0.5266%. The image is
that of a reinforced concrete with iron grid.
Another series of simulations are performed on a crack inside a concrete piece. The results are
presented in Figs. 6, 7 and 8 with the same image quality.
An analysis of concrete properties can be carried out in the same way as in X-ray CT. Here a notable
advantage resides in the fact that large doses of radiation are not needed due to a better confinement of radiation.
Moreover this scanning modality can be used on small objects in civil engineering. For large objects the
modality advocated by [18] is appropriate since the scanning by scattered radiation is done only on one side of
the object. However it presents an inconvenient displacement of the detector along a line passing through the
source. There is a possibility of using Compton scattering at angles larger than /2 with the present modality
to scan the half space of a large object. This has been already pointed out in [26].
We have presented in this work a new NDT/NDE modality based the phenomena of Compton
scattering by gamma rays produced by radioactive isotopes. The use of gamma radiation allows to get an
excellent spatial resolution as it is able to penetrate deeply in bulk material, hence perfectly adapted to concrete.
As there exists no universal NDT/NDE method, capable of providing all interesting parameters, our proposal
provides complementary information to the existing ones via the reconstruction of matter electron density. In
this way, the porosity of concrete can be located and evaluated, so that resistance to compression, elasticity
modulus and permeability of concrete can be deduced. This method has the advantage that it need not to be in
physical contact with the object to be investigated such as the geo-radar. It gives very convincing simulation
results and stands as a valid competitor for current methods of NDT/NDE. It is evident that not all the properties
of this Compton scatter tomography have been exploited. A major challenge is this research is the problem
posed by radiation attenuation. The way attenuation comes into the problem differs from that of the known
classical Radon transform. It leads to a very cumbersome inversion problem for which there is no solution at
present. However one may proceed by performing attenuation corrections, which have turned out to be very
successful in recent studies, see [27]. Efforts in this direction will be dispatched in the future.
The authors would like to thank the “Fondation de l’Université de Cergy-Pontoise”- France for financial support
to this work.
Carino N J. Nondestructive testing of concrete: history and challenges. ACI SP-144, Concrete
Technology - Past, Present and Future, Mehta P K, Ed., American Concrete Institute, Detroit, MI,
1994, 623-678.
la-Rucki G F and Eberhard M O, Imaging of reinforced concrete: State-of-the-Art Review, Journal of
Infrastructure Systems, 1995, 134-141.
Tuzi S, Sato O. Locating the positions of reinforcing bars in reinforced concrete using backscattered
gamma rays. Applied Radiation and Isotopes, 41,1990, 1013-1018.
17 | Page
A new Compton scattered tomography modality and its application to material non-destructive…
Tuzi S, Sato O, Determination of positions of reinforcing bars in reinforced concrete using
backscattered gamma rays II. Experimental and Monte Carlo methods, Applied Radiation and Isotopes,
44, 1993, 1279-1284.
Buyukozturk O, Imaging of concrete structures. NDT&E International, 31, 1998, 233-243.
Anghaie S, Humphries L L, Diaz N J, Material characterization and flaw detection, sizing, and location
by the differential gamma scattering spectroscopy technique. Part1: Development of theoretical basis,
NuclearTechnology, 91, 1990, 361-375.
Babot D, Le Floch C, Peix G, Contrôle et caractérisation des matériaux composites par
tomodensitométrie Compton, Revue Pratique du Contrôle Industriel 173, 1992, 64-67.
Berodias M G,Peix M G. Nondestructive measurement of density and effective atomic number by
photon scattering. Material Evaluation, 46(8), 1988, 1209-1213.
Gorshkov V A, Kroening M, Anosov Y V, Dorjgochoo O, X-ray scattering tomography, Nondestructive Test Evaluation, 20, 2005, 147-157.
Akkurt I and El-Khayatt A M, Effective atomic number and electron density of marble concrete,
Journal of Radioanalytical and Nuclear Chemistry, 295, 2013, 633-638.
Manohara S R, Hanagodimath S M, Thind K S and Gerward L, On the effective atomic number and
electron density: A comprehensive set of formulas for all types of materials and energies above 1 keV,
Nuclear Instruments and Methods in Physics Research, B266, 2008, 3906-3912.
de Oliveira Jr J M, Martins A C G, De Milito J A, Analysis of concrete material through gamma ray
computerized tomography, Brazilian Journal of Physics, 34(3A), 2004, 1020-1023.
Priyada P, Ramar R, Shivaramu, Application of gamma ray scattering technique for non-destructive
evaluation of voids in concrete, Applied Radiation and Isotopes, 74, 2013, 13-22.
Hussein E M A, Whynot T M, A Compton scattering method for inspecting concrete structures,
Nuclear Instrumentation Methods, A283, 1989, 100-106.
Kondic N N, Density field determination by an external stationary radiation source using a kernel
technique, Measurements in polyphase flows, ASME Winter Annual Meeting, San Francisco
California USA, December 10-15, 1978. Stock DE Ed., pp. 37-51.
Prettyman T H, Gardner R P, Russ J C, Verghese K, A combined transmission and scattering
tomographic approach to composition and density imaging, Appl. Radiat. Isot,.44(10/11), 1993,13271341.
Bodette D E, Jacobs A M, Tomographic two-phase flow distribution measurement using gamma-ray
scattering, ISA 1985, Paper 85-0362, pp. 23-34.
Norton S J, Compton scattering tomography, J. Appl. Phys., 76(4), 1994, 2007-2015.
Gautam S R, Hopkins F F, Klinksiek R, Morgan I L, Compton interaction tomography, IEEE
Transactions on Nuclear Science NS-30(2), 1983, 1680-1684.
Evans B L, Martin J B, Burggraf L W, Roggemann M C, Nondestructive inspection using Compton
scatter tomography. IEEE Transactions on Nuclear Science 45(3), 1998,950-956.
Harding G, X-ray scatter imaging in nondestructive testing. International Advances in Nondestructive
Testing 11, 1985, 271-295.
Holt R S, Cooper M J., Non-destructive examination with a Compton scanner. British Journal of
Nondestructive Testing, 30(2), 1988, 75-80.
Hussein E M A, Radiation scattering methods for nondestructive testing and imaging. Int. Advances
in Nondestructive Testing 14, 1989, 301-321.
Nguyen M K, Truong T T, Inversion of a new circular-arc Radon transform for Compton scattering
tomography. Inverse Problems 26, 2010, 065005.
Cormack A M, The Radon transform on a family of curves in the plane, Proceedings of the American
Mathematical Society, 83(2), 1981, 325-330.
Truong T T, Nguyen M K., Recent Developments on Compton Scatter Tomography: Theory and
Numerical Simulations, in Andriychuk M. (Ed) Numerical Simulation - From Theory to Industry,
Rijeka Craotia, INTECH; 2012, 101-128.
Rigaud G, Régnier R, Nguyen M K, New bimodal scattered radiation tomographic imaging with
attenuation and electron density correction algorithm. 38th International Conference on Acoustics,
Speech, and Signal Processing (ICASSP), Vancouver Canada, May 2013.
18 | Page
A new Compton scattered tomography modality and its application to material non-destructive…
19 | Page