Ray tracing is an integral part of modern imaging and inversion techniques. Several methods have been proposed depending on the requirements of the application. Some algorithms are best applied to fast changing material properties, like an interface between two differing media, while others are well suited to media with gradually changing properties, like composite materials. In this paper, an enhanced numerical algorithm for ray tracing is presented. Focus is given to solutions involving ordinary differential equations with initial-value conditions. The proposed algorithm is the result of a combination of two classical implementations and the authors show that it is well suited for media with both sharp and gradual changes in the index of refraction. Additionally, the authors present an application of ray path computation by using the technique known as the shooting method.

1.
Andersen
,
A.
, and
Kak
,
A.
(
1982
). “
Digital ray tracing in two-dimensional refractive fields
,”
J. Acoust. Soc. Am.
72
(
5
),
1593
1606
.
2.
Barry
,
G.
(
1963
). “
Ray tracings of acoustic waves in the upper atmosphere
,”
J. Atmos. Terr. Phys.
25
(
11
),
621
629
.
3.
Bold
,
G. E.
, and
Birdsall
,
T. G.
(
1986
). “
A top-down philosophy for accurate numerical ray tracing
,”
J. Acoust. Soc. Am.
80
(
2
),
656
660
.
4.
Brath
,
A. J.
, and
Simonetti
,
F.
(
2017
). “
Phased array imaging of complex-geometry composite components
,”
IEEE Trans. Ultrason. Ferroelectr. Freq. Control
64
(
10
),
1573
1582
.
5.
Camacho
,
J.
,
Cruza
,
J. F.
,
Brizuela
,
J.
, and
Fritsch
,
C.
(
2014
). “
Automatic dynamic depth focusing for NDT
,”
IEEE Trans. Ultrason. Ferroelectr. Freq. Control
61
(
4
),
673
684
.
6.
de Kool
,
M.
,
Rawlinson
,
N.
,
Hauser
,
J.
, and
Sambridge
,
M.
(
2006
). “
A practical grid-based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous media
,”
Geophys. J. Int.
167
,
253
270
.
7.
Denis
,
F.
,
Basset
,
O.
, and
Gimenez
,
G.
(
1995
). “
Ultrasonic transmission tomography in refracting media: Reduction of refraction artifacts by curved-ray techniques
,”
IEEE Trans. Med. Imag.
14
(
1
),
173
188
.
8.
Engdahl
,
E. R.
,
Sleep
,
N. H.
, and
Lin
,
M.-T.
(
1977
). “
Plate effects in North Pacific subduction zones
,”
Tectonophysics
37
(
1-3
),
95
116
.
9.
Holmes
,
C.
,
Drinkwater
,
B.
, and
Wilcox
,
P.
(
2004
). “
The post-processing of ultrasonic array data using the total focusing method
,”
Insight—Non-Destr. Test. Cond. Monit.
46
(
11
),
677
680
.
10.
Jakowatz
,
C. V.
, and
Kak
,
A. C.
(
1976
). “
Computerized tomographic imaging using x-rays and ultrasound, research report
,” TR-EE 76-26.
11.
Jeune
,
L.
,
Robert
,
S.
,
Dumas
,
P.
,
Membre
,
A.
, and
Prada
,
C.
(
2014
). “
Adaptive ultrasonic imaging with the total focusing method for inspection of complex components immersed in water
,” in
41st Annual Review of Progress in Quantitative Nondestructive Evaluation
, Boise, ID, 20 – 25 July, pp.
1037
1046
.
12.
Johnson
,
S.
,
Greenleaf
,
J.
,
Samayoa
,
W.
,
Duck
,
F.
, and
Sjostrand
,
J.
(
1975
). “
Reconstruction of three-dimensional velocity fields and other parameters by acoustic ray tracing
,” in
1975 Ultrasonics Symposium
, Los Angeles, CA, 22 – 24 September, pp.
46
51
.
13.
Julian
,
B. R.
, and
Gubbins
,
D.
(
1977
). “
Three-dimensional seismic ray tracing
,”
J. Geophys.
43
(
1
),
95
113
.
14.
Junghyun
Kwon
, and
Sang-Jin
Choi
, and
Song
,
S. M.-H.
(
2005
). “
Implementation and evaluation of the ultrasonic TOF tomography for the NDT of concrete structures
,” in
Computational Imaging III
, San Jose, CA, Vol.
5674
.
15.
Kak
,
A. C.
(
1979
). “
Computerized tomography with x-ray, emission, and ultrasound sources
,”
Proc. IEEE
67
(
9
),
1245
1272
.
16.
Keller
,
H. B.
(
1968
).
Numerical Methods for Two-Point Boundary Value Problems
, 1st ed. (
Blaisdell Publishing Co
.,
Waltham, MA
), Chap. 2.
17.
Kelso
,
J. M.
(
1967
). “
Ray tracing in the ionosphere
,” External Technical Memorandum No. 74.
18.
Langan
,
R. T.
,
Lerche
,
I.
, and
Cutler
,
R. T.
(
1985
). “
Tracing rays through heterogeneous media: An accurate and efficient procedure
,”
Geophysics
50
(
9
),
1456
1465
.
19.
Moser
,
T. J.
(
1991
). “
Shortest path calculation of seismic rays
,”
Geophysics
56
(
1
),
59
67
.
20.
Press
,
W. H.
,
Teukolsky
,
S. A.
,
Vetterling
,
W. T.
, and
Flannery
,
B. P.
(
2007
).
Numerical Recipes
, 3rd ed. (
Cambridge University Press
,
London
), Chap. 18.
21.
Rawlinson
,
N.
,
Hauser
,
J.
, and
Sambridge
,
M.
(
2008
). “
Seismic ray tracing and wavefront tracking in laterally heterogeneous media
,”
Adv. Geophys.
49
,
203
273
.
22.
Schomberg
,
H.
(
1978
). “
An improved approach to reconstructive ultrasound tomography
,”
J. Phys. D: Appl. Phys.
11
(
15
),
L181
.
23.
Smith
,
A.
,
Goldberg
,
M.
, and
Liu
,
E.
(
1980
). “
Numerical ray tracing in media involving continuous and discrete refractive boundaries
,”
Ultrason. Imaging
2
(
4
),
291
301
.
24.
Stoer
,
J.
, and
Bulirsch
,
R.
(
1980
).
Introduction to Numerical Analyisis
, 2nd ed. (
Springer-Verlag
,
Harrisonburg, VA
), Chap. 7.
25.
Turnbull
,
A.
, and
Garton
,
M.
(
1994
). “
Ultrasound ray tracing in arbitrary complex geometries
,” in
Review of Progress in Quantitative Nondestructive Evaluation
, Snowmass Village, CO, July 31–August 4, pp.
1105
1110
.
You do not currently have access to this content.