P1: JQX Mathematical Geology [mg] PP1377-matg-496105 November 3, 2004 7:56 Style file version June 25th, 2002
Mathematical Geology, Vol. 36, No. 8, November 2004 (C 2004)
The Bingham Distribution of Quaternions and Its Spherical Radon Transform in Texture Analysis1 Karsten Kunze2 and Helmut Schaeben3
Spherical geometry of quaternions is employed to characterize the Bingham distribution on the 3-dimensional sphere S3 ⊂ R4 as being uniquely composed of a bipolar, a circular and a spherical component. A new parametrization of its dispersion parameters provides a classification of patterns of crystallographic preferred orientations (CPO, or textures). It is shown that the Bingham distribution can represent most types of ideal CPO patterns; in particular single component, fiber and surface textures are represented by rotationally invariant bipolar, circular and spherical distributions, respec- tively. Pole figures of given crystal directions are derived by the spherical Radon transform of the Bingham probability density function of rotations, which are displayed for general and special cases.
KEY WORDS: rotations, crystallographic preferred orientations, texture components, fiber textures, pole figures, spherical statistics.
INTRODUCTION
Crystallographic preferred orientations (CPO) of minerals in rocks are commonly used to infer constraints on their tectono-metamorphic history (e.g., Wenk, 1985). They are usually presented and discussed as pole figures (fabric diagrams) of some specific crystal planes or directions. The actual frequency distribution of crystal orientations typically receives less attention, even though it provides the basic link between the diverse pole figures, because it is harder to visualize in an intuitive fashion. Model orientation distributions aid in interpreting CPO patterns, as they are characterized by a rather small number of meaningful parameters, which may be the centers of concentration or depletion, the degree and anisotropy of scattering, or certain patterns of rotational symmetry and invariance. They provide hypotheses
1Received 17 December 2003; accepted 27 July 2004. 2Department of Earth Sciences, ETH Zurich, CH-8092 Zurich, Switzerland; e-mail: karsten.kunze@ erdw.ethz.ch 3Geoscience Mathematics and Informatics, University of Mining and Technology, D-09596 Freiberg, Germany; e-mail: [email protected]
917
0882-8121/04/1100-0917/1 C 2004 International Association for Mathematical Geology P1: JQX Mathematical Geology [mg] PP1377-matg-496105 November 3, 2004 7:56 Style file version June 25th, 2002
918 Kunze and Schaeben
for measured data to be tested against by statistical methods. The Bingham distri- bution is well suited for this purpose, as its parameter domain allows to simulate a wide range of different CPO; preliminary results were given by Schaeben (1996). Neglecting crystal symmetry for the sake of simplicity, a crystallographic orientation is a rotation between the individual crystal and the specimen. The frequency of crystal orientations in a representative polycrystalline sample volume 1 is given by the orientation (probability) density function (ODF) f : SO(3) → R+. If rotations are represented by unit quaternions, the ODF is an even function defined on S3 ⊂ R4. Pole figures are derived by the spherical Radon transform R f of the ODF f , representing the probability distribution of a given unit vector h ∈ S2 ⊂ R3, i.e., of a specific direction in a crystal, which has been subjected to a random rotation with a given—say Bingham—distribution f . In this contribution we will show that the Bingham distribution of unit quater- nions is uniquely composed of distinct amounts to a bipolar, a circular and a spherical component, which is made transparent by a new parametrization of the dispersion parameters. We will particularly discuss the limiting cases, for which the Bingham distribution degenerates to rotationally invariant distributions, and which are related to known basic CPO types such as single orientation components, fiber textures and surface textures. It is further demonstrated, how transitional cases of the Bingham distribution relate to intermittent texture types. Pole figures for general and special cases are derived by the spherical Radon transform. It is shown by example that very different, even complementary, orientation distributions may give rise to hardly distinguishable pole figures, at least for certain crystal directions, stressing the point to rely on the entire orientation distribution when characterizing the CPO of a rock.
ROTATIONS AND QUATERNIONS
The orientation g of an individual crystal in a polycrystalline sample is the active rotation g ∈ SO(3) : KS → KC that maps a right-handed orthonormal co- ordinate system KS fixed to the specimen onto another right-handed orthonormal coordinate system KC fixed to the crystal,
gKS = KC, g ∈ SO(3). (1)
If a unique direction is represented by unit vector h with respect to the crystal frame KC, and by unit vector r with respect to the specimen frame KS , then the coordinates of the unique direction transform according to
= . rKS ghKC (2) P1: JQX Mathematical Geology [mg] PP1377-matg-496105 November 3, 2004 7:56 Style file version June 25th, 2002
Bingham Distribution and Spherical Radon Transform in Texture Analysis 919
The rotation g may be represented by a (3 × 3) orthogonal matrix X(g) such that
r = X(g)h,
where matrix multiplication applies, or by a unit quaternion q(g) ∈ S3 such that
r = q(g)h q∗(g) (3)
where quaternion multiplication applies. The quaternion representation will be pursued here as it proved particularly appropriate to treat statistics of rotations by spherical statistics. Basic properties of quaternions are compiled in the Appendix.
Quaternion Representation of Rotations in R3
Equations (2) and (3) explicitly read (e.g., Altmann, 1986)
r = g h = q h q∗ = h cos ω + (n × h) sin ω + (n · h) n(1 − cos ω)
using the parametrization of unit quaternions by rotation axis n and rotation angle ω [Appendix, Eq. (33)]. Both unit quaternions ±q ∈ S3 represent one and the same rotation in R3, i.e., S3 stands for a double covering of the group SO(3). The unit sphere S2 ⊂ VecH consists of all quaternions representing rotations by the angle π about arbitrary axes in S2, as every unit vector q can be considered as the pure quaternion q = cos(π/2) + q sin(π/2). The unit quaternions ±q∗ represent the inverse rotation g−1r = h by angle ω and axis −n. Henceforward, no distinction will be made between the rotation g and its (up to the sign) unique quaternion representation q. Similarly, a unit vector q ∈ S2 ⊂ VecH is identified with its corresponding pure unit quaternion q = (0 + q) ∈ S3 ⊂ H, for which the rules of quaternion multiplication apply. Both will be denoted here by q; the distinction is made by context. In terms of passive rotations, Equation (2) provides the coordinate transfor-
mation of a unique vector given by hKC or rKS with respect to the crystal (KC)or the specimen (KS ) coordinate system, respectively, that are related to each other by the rotation g: KS → KC. The commonly applied convention in texture analysis (Bunge, 1982) is — in a strict mathematical sense incorrectly — to use one and the same unique symbol g to refer to both the (active) rotation g of coordinate systems = −1 = −1 g KS KC and the (passive) transformation g of unit vectors hKC g rKS . Formal consistency based on the vector transformation rule [Eq. (2)] is maintained P1: JQX Mathematical Geology [mg] PP1377-matg-496105 November 3, 2004 7:56 Style file version June 25th, 2002
920 Kunze and Schaeben
between previous references and the formulae given here by identifying g = g−1, and vice versa. Additional details concerning the geometry of rotations when represented by unit quaternions are given by Meister and Schaeben (2004).
Parameterizations of Unit Quaternions
3 Any set of four orthonormal quaternions ai ∈ S , i = 1,...,4, spans the sphere S3 of all unit quaternions, which may be parameterized in either way by
q(r, s, t) = ((a1 cos t + a2 sin t) sin s + a3 cos s) sin r + a4 cos r, or (4)
q(r, s, t) = (a1 cos t + a2 sin t) sin s + (a3 sin r + a4 cos r) cos s. (5) S3 ={q(r, s, t) ∈ H | r, s ∈ [0,π], t ∈ [0, 2π)}
If the quaternions ai are identified with the basis in H by ai = ei , i = 1,...,3, and a4 = 1, then the variables are easily interpreted. Since any rotation is represented twice in S3 (by q and −q), the parameter domains for the following rotations cover half of S3. • In the first case [Eq. (4)] they correspond to rotation angle ω ∈ [0,π] and axis n ∈ S2 as introduced by Equation (33) according to
r = ω/2, s = θ, t = ϕ, (6)
where θ ∈ [0,π] and ϕ ∈ [0, 2π) are the spherical angles of the rotation axis n = (cos ϕ sin θ, sin ϕ sin θ, cos θ)t . • In the second case [Eq. (5)] they correspond to linear combinations of the Euler angles (φ1, ,φ2) conventionally used in texture analysis (Bunge, 1982) with (φ1,φ2 ∈ [0, 2π); ∈ [0,π])
r = (φ1 + φ2)/2, s = /2, t = (φ1 − φ2)/2. (7)
• The variables in Equation (5) may alternatively be interpreted in terms of a product pq of two rotations represented by p, q ∈ S3 about mutually per- pendicular axes. Let p and q have rotation axes parallel and perpendicular to e3, respectively, i.e., ω ω p = cos p + sin p e , 2 2 3 ω ω q = cos q + sin q (cos ϕ e + sin ϕ e ), 2 2 q 1 q 2 P1: JQX Mathematical Geology [mg] PP1377-matg-496105 November 3, 2004 7:56 Style file version June 25th, 2002
Bingham Distribution and Spherical Radon Transform in Texture Analysis 921
with ωp,ϕq ∈ [0, 2π); ωq ∈ [0,π]. Then the parameters are given by
r = ωp/2, s = ωq /2, t = ωp/2 + ϕq .
Points, Circles and Spheres in S3
3 Let ai ∈ S , i = 1,...,4, be any set of four orthonormal quaternions. Let L(a4) ⊂ H be a line parallel to a4, E(a3, a4) ⊂ H be a plane spanned by a3, a4, and F(a2, a3, a4) ⊂ H be a hyperplane spanned by a2, a3, a4. Their intersections 3 with S constitute a single bipolar unit quaternion B(a4), a unit circle C(a3, a4), or a unit sphere S(a2, a3, a4), which are of dimension zero, one or two, respectively, and given by
3 B(a4):= L(a4) ∩ S ={−a4, +a4}, 3 C(a3, a4):= E(a3, a4) ∩ S 3 ={q(t) ∈ S | q(t) = a4 cos t + a3 sin t, t ∈ [0, 2π)}, 3 (8) S(a2, a3, a4):= F(a2, a3, a4) ∩ S 3 ={q(s, t) ∈ S | q(s, t) = (a4 cos t + a3 sin t) sin s + a2 cos s, × s ∈ [0,π], t ∈ [0, 2π)}.
, , = ⊥ Obviously, B(a1) and S(a2 a3 a4) : a1 are mutually orthogonal complements 3 with respect to S , and so are C(a3, a4) and C(a1, a2). These circles and spheres are centered at O, the origin of the coordinate system; thus they are geodetic lines or surfaces, respectively. The orientation distances between x ∈ S3 and a bimodal point, a circle or a sphere are derived from orthogonal projections of x onto L, E,orF, respectively, as ω cos2 b = (Sc(a∗x))2 = 1 − (Sc(a∗x))2 − (Sc(a∗x))2 − (Sc(a∗x))2, 2 4 3 2 1 ω cos2 c = (Sc(a∗x))2 + (Sc(a∗x))2 = 1 − (Sc(a∗x))2 − (Sc(a∗x))2, (9) 2 4 3 2 1 ω cos2 s = (Sc(a∗x))2 + (Sc(a∗x))2 + (Sc(a∗x))2 = 1 − (Sc(a∗x))2, 2 4 3 2 1
where ωb denotes the orientation distance between x and the bimodal point B(a4), ωc the one between x and the circle C(a3, a4), and ωs the one between x and the sphere S(a2, a3, a4). 2 For vectors hc, rc ∈ S defined by
= ∗ = ∗ =− ∗ = ∗ = ∗ =− ∗ , hc : Vec(a4 a3) a4 a3 a3 a4 Vec(a2 a1) a2 a1 a1 a2 = ∗ = ∗ =− ∗ =− ∗ =− ∗ = ∗, (10) rc : Vec(a3a4 ) a3a4 a4a3 Vec(a1a2 ) a1a2 a2a1 P1: JQX Mathematical Geology [mg] PP1377-matg-496105 November 3, 2004 7:56 Style file version June 25th, 2002
922 Kunze and Schaeben
the circle C(a3, a4) represents the set of all rotations mapping hc on rc, i.e., ∗ q hc q = rc for all q ∈ C(a3, a4). The complementary circle C(a1, a2) is the set of all rotations mapping hc on −rc (Kunze, 1991; Meister and Schaeben, 2004). ω ω = · ∗ The orientation distance c is also given by cos c rc x hc x as follows using = 4 i i = ∗ the decomposition x i=1 x ai with x Sc(ai x)
· ∗ =− ∗ =− ∗ ∗ ∗ rc xhc x Sc(rc x hc x ) Sc(a3a4 xa4 a3 x ) 4 4 =− i j ∗ ∗ ∗ =− 1 2− 2 2+ 3 2+ 4 2 x x Sc(a3a4 ai a4 a3 a j ) (x ) (x ) (x ) (x ) i=1 j=1 ω = 2 [(Sc(a∗x))2 + (Sc(a∗x))2] − 1 = 2 cos2 c − 1 = cos ω . 3 4 2 c
, ∈ S2 ∗ = , , = ⊥ = For any pair h r with a1ha1 r, the sphere S(a2 a3 a4) a1 2 3 S a1 ⊂ S ⊂ H contains the circle representing all rotations mapping h on −r, and it does not contain any other circle mapping h on a direction different of −r (Meister and Schaeben, 2004).
Invariant Measures in S3
Let S3 be embedded in R4 so that p = cq ∈ R4, q ∈ S3, c ∈ R. Then the in- variant measure on S3 is obtained from the (unsigned) Jacobian for the transforma- tion from (p0, p1, p2, p3)to(c, r, s, t) for c = 1, which reads from Equations (4) and (6)
1 ω ds = sin2 r sin sdrdsdt = sin2 sin θ dω dθ dϕ, 3 2 2
or from Equations (5) and (7) it is
1 ds = sin s | cos s| dr ds dt = sin d dφ dφ . (11) 3 8 1 2
/ = φ φ = 2 ω θ ω θ ϕ Notice that ds3 is 1 8ofdg sin d d 1 d 2 4 sin 2 sin d d d , = π 2 which is commonly defined in texture analysis yielding SO(3) dg 8 . Further recall that SO(3) is covered twice by S3. From Equations (4) and (6), the invariant measure on the sphere is analogously obtained for c = 1,ω = 2r = π
ds2 = sin sdsdt = sin θ dθ dϕ, (12) P1: JQX Mathematical Geology [mg] PP1377-matg-496105 November 3, 2004 7:56 Style file version June 25th, 2002
Bingham Distribution and Spherical Radon Transform in Texture Analysis 923
and the invariant measure on the circle is found for c = 1,ω= 2r = π, θ = s = π/2as
ds1 = dt = dϕ.
Integration results in