A variational formulation of the fast marching eikonal solver

Next: Variational principles on a Up: Fomel: Fast marching Previous: A brief description of

# The theoretic grounds of variational principles

This section serves as a brief reminder of the well-known theoretical connection between Fermat's principle and the eikonal equation. The reader, familiar with this theory, can skip safely to the next section.

fermat
Figure 2.
Illustration of the connection between Fermat's principle and the eikonal equation. The shortest distance between a wavefront and a neighboring point is along the wavefront normal.

Both Fermat's principle and the eikonal equation can serve as the foundation of traveltime calculations. In fact, either one can be rigorously derived from the other. A simplified derivation of this fact is illustrated in Figure 2. Following the notation of this figure, let us consider a point in the immediate neighborhood of a wavefront . Assuming that the source is on the other side of the wavefront, we can express the traveltime at the point as the sum

 (3)

where is a point on the front, is the length of the ray segment between and , and is the local slowness. As follows directly from equation (3),
 (4)

Here denotes the angle between the traveltime gradient (normal to the wavefront surface) and the line from to , and is the directional traveltime derivative along that line.

If we accept the local Fermat's principle, which says that the ray from the source to corresponds to the minimum-arrival time, then, as we can see geometrically from Figure 2, the angle in formula (4) should be set to zero to achieve the minimum. This conclusion leads directly to the eikonal equation (2). On the other hand, if we start from the eikonal equation, then it also follows that , which corresponds to the minimum traveltime and constitutes the local Fermat's principle. The idea of that simplified proof is taken from Lanczos (1966), though it has obviously appeared in many other publications. The situations in which the wavefront surface has a discontinuous normal (given raise to multiple-arrival traveltimes) require a more elaborate argument, but the above proof does work for first-arrival traveltimes and the corresponding viscosity solutions of the eikonal equation (Lions, 1982).

The connection between variational principles and first-order partial-differential equations has a very general meaning, explained by the classic Hamilton-Jacobi theory. One generalization of the eikonal equation is

 (5)

where represents the vector of space coordinates, and the coefficients form a positive-definite matrix . Equation (5) defines the characteristic surfaces for a linear hyperbolic second-order differential equation of the form
 (6)

where F is an arbitrary function.

A known theorem (Smirnov, 1964) states that the propagation rays [characteristics of equation (5) and, correspondingly, bi-characteristics of equation (6)] are geodesic (extreme-length) curves in the Riemannian metric

 (7)

where are the components of the matrix . This means that a ray path between two points and has to correspond to the extreme value of the curvilinear integral

For the isotropic eikonal equation (2), , and metric (7) reduces to the familiar traveltime measure
 (8)

where is the usual Euclidean distance metric. In this case, the geodesic curves are exactly Fermat's extreme-time rays.

From equation (7), we see that Fermat's principle in the general variational formulation applies to a much wider class of situations if we interpret it with the help of non-Euclidean geometries.

 A variational formulation of the fast marching eikonal solver

Next: Variational principles on a Up: Fomel: Fast marching Previous: A brief description of

2013-03-03