[MSU] Lomonosov Moscow State University
Research Computer Centre

Laboratory of Automated Image Processing Systems
tel: +7(495)939-2759, fax: +7(495)939-2768
English English
 

 Inverse problems of Diagnostics

  X-ray tomography     Ultrasonic tomography     Synthetic aperture  

 

Ultrasonic tomography

Inverse problems of wave tomography arise in various applications. These include some of the inverse problems of electromagnetic sounding, seismology, civil engineering, sonar, industrial diagnosis via acoustic radiation, etc.

One of the major applications of wave tomography is ultrasound tomography in medicine. The primary task of ultrasound tomography is to perform differential diagnosis of breast cancer one of the most urgent problems in modern medicine. More than 40% of the patients have advanced breast cancer diagnosed at the first examination.

Prototypes of ultrasound tomography devices, which have a number of advantages over magnetic-resonance and x-ray tomography, are now developed extensively in the USA, Germany, Japan, and Russia [1, 2, 3], Mathematically, the inverse problems arising in ultrasound tomography are much more challenging compared to those of x-ray tomography. Ray optics model works well in x-ray tomography where inverse problems are linear and can be solved on a common PC.

Unlike x-ray tomography, ultrasound tomography uses both transmitted and reflected radiation emitted by broadband sounding sources. Detectors and sources are placed on all sides of the object examined and the signal is recorded as a time series, which makes it possible to process a very large volume of data, allowing high resolution to be achieved.

Inverse problems of ultrasound tomography can be treated as coefficient inverse problems. Let us now describe the formulation of the inverse problem. Figure 1(a) shows the schemes of tomographic experiments in the case of traditional layer-by-layer ("2.5D") tomography, where two-dimensional problems are solved on layers of a three-dimensional object. Figure 1(b) shows the scheme of 3D tomography, where the acquired experimental data are used to solve the inverse problem right away for the three-dimensional unknown function c(r), r∊R3. In Figure 1(a) the sources are labeled by digit "1" and the detectors arranged along the perimeter of the computational domain are labeled by digit "2". The sources S and detectors of ultrasound pulses in Figure 1(b) are located on the faces of the cube.

 

(a) (b)
Figure 1. Ultrasound tomography schemes:   (a) layer-by-layer 2D scheme,   (b) - 3D scheme.

 

Consider now the wave equation that describes the acoustic pressure field u(r,t) in the domain Ω⊂RN, N=2,3, confined by surface S over the time interval (0,T) with a point source located at a point r0 .

c(r)utt(r,t)-Δu(r,t)=δ(r-r0)˙f(t)

u(r,t=0)=ut(r,t=0)=0,   ∂nu|ST=p(r,t).

Here c-0.5(r)=v(r) is the unknown wave velocity in the diagnosed region. The pulse generated by the source is described by the given function f(t), nu|ST is the derivative along the normal to surface S in region S×(0,T), and p(r,t) is a known function. We further assume that the inhomogeneity of the medium is due solely to velocity variations, and outside the inhomogeneity region c(r)≡c0=const, where c0 is known. The inverse problem consists in finding a function c(r) that describes the inhomogeneity given the measurements of acoustic pressure U(r,t) at the boundary S of the domain Ω with various source positions r0 . The inverse problem is nonlinear. The number of unknowns to be determined when solving the inverse problem of ultrasound tomography on a 320×320×320 grid exceeds 10 million.

Mathematically the inverse problem is formulated as the problem of minimization of the residual functional Φ(u(c)) in the norm of L2 space  (S×(0,T)):

T
Φ(u(c))=
1
2
(u(s,t)-U(s,t))2ds dt = ||u(s,t)-U(s,t)||2
0S

The breakthrough results in the solution of ultrasound tomography problems are associated with the possibility of explicitly, analytically formulating the gradient of the residual functional [4, 5, 6]. The approach that we propose differs fundamentally from those currently employed [1-3], which use substantially higher sounding frequencies. The iterative algorithms that we developed to solve inverse problems of ultrasound tomography and the simulations performed over a wide range of parameter values allow the sounding frequencies to be reduced several times. This point is of fundamental importance, because attenuation in soft tissues depends heavily on sounding frequency. The lower the frequency, the smaller the attenuation, and hence the smaller the error of experimental data. Simulations show that about 5-mm wavelength ultrasound allows a resolution of about 2-3 mm to be achieved. In our paper [7] we propose efficient methods to solve inverse problems of ultrasound tomography in terms of models that take into account attenuation in the medium.

Figure 2 shows the results of the reconstruction in the case of layer-by-layer tomography. Figure 2(a) shows the image of a 2D phantom. Figures 2(b) and 2(c) show the reconstructed velocity structure in the case of zero- and 1.5%-error input data, respectively. The simulations were performed for eight sources with the detectors spaced 2.5 mm, i.e., half the wavelength, apart. The size of the sounding region is 200×200 mm.

 

(a) (b) (c)
Figure 2. Layer-by-layer scheme:
(a) phantom,   (b) reconstruction from exact input data;   (c) reconstruction from noisy input data with 1.5% error.

 

The minimum size of inhomogeneities in the phantom is about 2 mm. As is evident from a comparison of the initial and reconstructed images, the resulting resolution is 2-3 mm or better.

Figures 3 and 4 illustrate the capabilities of 3D tomographic schemes involving the reconstruction of the three-dimensional function (x,,z). Figures 3 (a-c) shows the cross sections of the three-dimensional phantom in several planes Z=const and Fig.3 (d-f), the cross sections of the reconstructed 3D velocity structure (x,,z) in the same planes. In our simulations we used 17 sources in a 3D tomography scheme without elements of rotation. The computations were performed on a 320×320×320 grid. As is evident from the figures, the quality of reconstruction is high, allowing inhomogeneities as small as 2mm to be reconstructed fairly well.

 

(a) (b) (c)
(d) (e) (f)
Figure 4. Cross-sections of the 3D velocity structure in Z=const planes:
(a)-(c) images of the phantom (exact sound speed);   (d)-(f) reconstructed images

 

Figure 4 shows the vertical cross sections of the phantom (a) and of the reconstructed image (b) in the X=const plane.

 

a) b)
Figure 4. Cross-section of the 3D velocity structure in X=const plane:
(a) phantom;   (b) reconstructed image

 

Figure 5 shows the profile of the reconstructed sound velocity structure along one of the A-A cross sections (shown in Figure 3(e)). As is evident from the figure, the absolute sound velocity values are reconstructed very accurately. The reconstructed and exact (phantom) sound velocity values are shown by the solid and dashed lines, respectively. We started all our simulations with the initial approximation c0(r)=const. Note that common ultrasound facilities operating in the reflection mode are in principle incapable of reconstructing the absolute value of sound velocity in an inhomogeneous medium.

 

Figure 5. Plot of the reconstructed velocity structure along the A-A line

 

Three-dimensional inverse problems of ultrasound tomography can be solved only with supercomputers. The algorithms that we developed to this end are highly scalable, up to petaflops-level supercomputers based on general-purpose processors, and were tested on "Lomonosov" supercomputer of Moscow State University Supercomputing Center (Moscow, Russia). However, because of their technical parameters, general-purpose supercomputers cannot be incorporated into ultrasound tomographic facilities. We developed efficient algorithms that allow inverse problems of ultrasound tomography both in the layer-by-layer 2.5D and explicit 3D formulation to be solved on GPU-based supercomputers. The number of parallel operating graphical processors is equal to that of the sources employed. The technical parameters of GPU-based supercomputers allow them to be incorporated into tomographic facilities as autonomous computing units.

References

1. J. Wiskin, D. T. Borup, S. A. Johnson, and M. Berggren. Non-linear inverse scattering: High resolution quantitative breast tissue tomography, - J. Acoust. Soc. Am. Volume 131, Issue 5, pp. 3802-3813 (2012);

2. Duric N, Littrup P, Poulo L, Babkin A, Pevzner R, Holsapple E, Rama 0, Glide C. Detection of breast cancer with ultrasound tomography: First results with the Computed Ultrasound Risk Evaluation (CURE) prototype. Medical Physics 2007;34:773-785.

3. H. Gemmeke, A. Menshikov, D. Tchernikovski, L. Berger, G. G?bel, M. Birk, M. Zapf and N.V. Ruiter: Hardware Setup for the Next Generation of 3D Ultrasound Computer Tomography. IEEE NSS MIC 2010.

4. Natterer F., Wubbeling F., A propagation-backpropagation method for uzltrasound tomography. Inverse Problems 11, 1225-1232 (1995) IOP Publishing Ltd.

[PDF] 5. Goncharsky AV, Romanov SY. Supercomputer technologies in inverse problems of ultrasound tomography. Inverse Probl. 2013;29:075004.

[PDF] 6. Goncharsky AV, Romanov SY, Seryozhnikov SY. Inverse problems of 3D ultrasonic tomography with complete and incomplete range data. Wave motion 2014;51:389-404.

[PDF] 7. Goncharsky AV, Romanov SY. Inverse problems of ultrasound tomography in models with attenuation. Phys Med Biol. 2014 Apr 21;59(8):1979-2004. doi: 10.1088/0031-9155/59/8/1979. Epub 2014 Apr 2.

8. Goncharsky A.V, Romanov S.Y, Seryozhnikov S.Y. A computer simulation study of soft tissue characterization using low-frequency ultrasonic tomography. Ultrasonics 2016;67:136-150.

[PDF] 9. Goncharsky A.V, Romanov S.Y, Seryozhnikov S.Y. SUPERCOMPUTER TECHNOLOGIES IN DEVELOPMENT OF DIAGNOSTIC TOMOGRAPHY FACILITIES (in Russian) - SPb.:Politechnical University publishing, 2016. - 230 p.