Affine curvature

From Wikipedia, the free encyclopedia
Jump to: navigation, search
This article is about the curvature of affine plane curves, not to be confused with the curvature of an affine connection.

Special affine curvature, also known as the equi-affine curvature or affine curvature, is a particular type of curvature that is defined on a plane curve that remains unchanged under a special affine transformation (an affine transformation that preserves area). The curves of constant equi-affine curvature k are precisely all non-singular plane conics. Those with k > 0 are ellipses, those with k = 0 are parabolas, and those with k < 0 are hyperbolas.

The usual Euclidean curvature of a curve at a point is the curvature of its osculating circle, the unique circle making second order contact (having three point contact) with the curve at the point. In the same way, the special affine curvature of a curve at a point P is the special affine curvature of its hyperosculating conic, which is the unique conic making fourth order contact (having five point contact) with the curve at P. In other words it is the limiting position of the (unique) conic through P and four points P1, P2, P3, P4 on the curve, as each of the points approaches P:

P_1,P_2,P_3,P_4\to P.

In some contexts, the affine curvature refers to a differential invariant κ of the general affine group, which may readily obtained from the special affine curvature k by κ = k−3/2dk/ds, where s is the special affine arc length. Where the general affine group is not used, the special affine curvature k is sometimes also called the affine curvature (Shirokov 2001b).

Formal definition[edit]

Special affine arclength[edit]

To define the special affine curvature, it is necessary first to define the special affine arclength (also called the equi-affine arclength). Consider an affine plane curve \beta (t). Choose co-ordinates for the affine plane such that the area of the parallelogram spanned by two vectors a = (a_1, \; a_2) and b = (b_1, \; b_2) is given by the determinant

\det\left[ a\; b \right] = a_{1} b_{2} - a_{2} b_{1}.

In particular, the determinant

\det\begin{bmatrix}\frac{d\beta}{dt} & \frac{d^2\beta}{dt^2}\end{bmatrix}

is a well-defined invariant of the special affine group, and gives the signed area of the parallelogram spanned by the velocity and acceleration of the curve β. Consider a reparameterization of the curve β, say with a new parameter s related to t by means of a regular reparameterization s = s(t). This determinant undergoes then a transformation of the following sort, by the chain rule:

\begin{align}
\det\begin{bmatrix}\frac{d\beta}{dt} & \frac{d^2\beta}{dt^2}\end{bmatrix} &= \det\begin{bmatrix}\frac{d\beta}{ds}\frac{ds}{dt} & \left(\frac{d^2\beta}{ds^2}\left(\frac{ds}{dt}\right)^2+\frac{d\beta}{ds}\frac{d^2s}{dt^2}\right)\end{bmatrix}\\
&=\left(\frac{ds}{dt}\right)^3\det\begin{bmatrix}\frac{d\beta}{ds} & \frac{d^2\beta}{ds^2}\end{bmatrix}.
\end{align}

The reparameterization can be chosen so that

\det\begin{bmatrix}\frac{d\beta}{ds} & \frac{d^2\beta}{ds^2}\end{bmatrix} = 1

provided the velocity and acceleration, dβ/dt and d2β/dt2 are linearly independent. Existence and uniqueness of such a parameterization follows by integration:

s(t) = \int_a^t\sqrt[3]{\det\begin{bmatrix}\frac{d\beta}{dt} & \frac{d^2\beta}{dt^2}\end{bmatrix}}\,\,dt.

This integral is called the special affine arclength, and a curve carrying this parameterization is said to be parameterized with respect to its special affine arclength.

Special affine curvature[edit]

Suppose that β(s) is a curve parameterized with its special affine arclength. Then the special affine curvature (or equi-affine curvature) is given by

k(s) = \det\begin{bmatrix}\beta''(s) & \beta'''(s) \end{bmatrix}.

Here β′ denotes the derivative of β with respect to s.

More generally (Guggenheimer 1977, §7.3; Blaschke 1923, §5), for a plane curve with arbitrary parameterization

t \mapsto (x(t), y(t)),

the special affine curvature is:


\begin{align}
k(t)&=\frac{x''y'''-x'''y''}{(x'y''-x''y')^{5/3}}-\frac{1}{2}\left[\frac{1}{(x'y''-x''y')^{2/3}}\right]''\\
&= \frac{4(x''y'''-x'''y'')+(x'y''''-x''''y')}{3(x'y''-x''y')^{5/3}} -\frac{5}{9}\frac{(x'y'''-x'''y')^2}{(x'y''-x''y')^{8/3}}
\end{align}

provided the first and second derivatives of the curve are linearly independent. In the special case of a graph y = y(x), these formulas reduce to

k=-\frac{1}{2}\left(\frac{1}{(y'')^{2/3}}\right)''=\frac{1}{3}\frac{y''''}{(y'')^{5/3}}-\frac{5}{9}\frac{(y''')^2}{(y'')^{8/3}}

where the prime denotes differentiation with respect to x (Blaschke 1923, §5; Shirokov 2001a).

Affine curvature[edit]

Suppose as above that β(s) is a curve parameterized by special affine arclength. There are a pair of invariants of the curve that are invariant under the full general affine group (Shirokov 2001b) — the group of all affine motions of the plane, not just those that are area-preserving. The first of these is

\sigma = \int \sqrt{k(s)}\, ds,

sometimes called the affine arclength (although this risks confusion with the special affine arclength described above). The second is referred to as the affine curvature:

\kappa = \frac{1}{k^{3/2}} \frac{dk}{ds}.

Conics[edit]

Suppose that β(s) is a curve parameterized by special affine arclength with constant affine curvature k. Let

C_\beta(s) = \begin{bmatrix}\beta'(s) & \beta''(s)\end{bmatrix}.

Note that det Cβ, since β is assumed to carry the special affine arclength parameterization, and that

k = \det(C_\beta').\,

It follows from the form of Cβ that

C_\beta' = C_\beta\begin{bmatrix}0&-k\\1&0\end{bmatrix}.

By applying a suitable special affine transformation, we can arrange that Cβ(0) = I is the identity matrix. Since k is constant, it follows that Cβ is given by the matrix exponential

\begin{align}
C_\beta(s) &= \exp\left\{s\cdot\begin{bmatrix}0&-k\\1&0\end{bmatrix}\right\}\\
&=\begin{bmatrix}\cos\sqrt{k}s&\sqrt{k}\sin\sqrt{k}s\\ -\frac{1}{\sqrt{k}}\sin\sqrt{k}s&\cos\sqrt{k}s\end{bmatrix}.
\end{align}

The three cases are now as follows.

k = 0

If the curvature vanishes identically, then upon passing to a limit,

C_\beta(s) = \begin{bmatrix}1&0\\s&1\end{bmatrix}

so β'(s) = (1,s), and so integration gives

\beta(s)=(s,s^2/2)\,

up to an overall constant translation, which is the special affine parameterization of the parabola y = x2/2.

k > 0

If the special affine curvature is positive, then it follows that

\beta'(s) = \left(\cos\sqrt{k}s,\frac{1}{\sqrt{k}}\sin\sqrt{k}s\right)

so that

\beta(s) = \left(\frac{1}{\sqrt{k}}\sin\sqrt{k}s, -\frac{1}{k}\cos\sqrt{k}s\right)

up to a translation, which is the special affine parameterization of the ellipse kx2 + k2y2 = 1.

k < 0

If k is negative, then the trigonometric functions in Cβ give way to hyperbolic functions:

C_\beta(s) =\begin{bmatrix}\cosh\sqrt{|k|}s&\sqrt{|k|}\sinh\sqrt{|k|}s\\ \frac{1}{\sqrt{|k|}}\sinh\sqrt{|k|}s&\cosh\sqrt{|k|}s\end{bmatrix}.

Thus

\beta(s) = \left(\frac{1}{\sqrt{|k|}}\sinh\sqrt{|k|}s,\frac{1}{|k|}\cosh\sqrt{|k|}s\right)

up to a translation, which is the special affine parameterization of the hyperbola

-|k|x^2 + |k|^2y^2 = 1.

Characterization up to affine congruence[edit]

The special affine curvature of an immersed curve is the only (local) invariant of the curve in the following sense:

  • If two curves have the same special affine curvature at every point, then one curve is obtained from the other by means of a special affine transformation.

In fact, a slightly stronger statement holds:

  • Given any continuous function k : [a,b] → R, there exists a curve β whose first and second derivatives are linearly independent, such that the special affine curvature of β relative to the special affine parameterization is equal to the given function k. The curve β is uniquely determined up to a special affine transformation.

This is analogous to the fundamental theorem of curves in the classical Euclidean differential geometry of curves, in which the complete classification of plane curves up to Euclidean motion depends on a single function κ, the curvature of the curve. It follows essentially by applying the Picard–Lindelöf theorem to the system

C_\beta' = C_\beta\begin{bmatrix}0&-k\\1&0\end{bmatrix}

where Cβ = [β′ β′′]. An alternative approach, rooted in the theory of moving frames, is to apply the existence of a primitive for the Darboux derivative.

Derivation of the curvature by affine invariance[edit]

The special affine curvature can be derived explicitly by techniques of invariant theory. For simplicity, suppose that an affine plane curve is given in the form of a graph y = y(x). The special affine group acts on the Cartesian plane via transformations of the form

\begin{align}
x&\mapsto ax+by + \alpha\\
y&\mapsto cx+dy + \beta,
\end{align}

with ad − bc = 1. The following vector fields span the Lie algebra of infinitesimal generators of the special affine group:

T_1 = \partial_x, \quad T_2 = \partial y
X_1 = x\partial_y, \quad X_2 = y\partial_x, \quad H=x\partial_x - y\partial_y.

An affine transformation not only acts on points, but also on the tangent lines to graphs of the form y = y(x). That is, there is an action of the special affine group on triples of coordinates

(x,y,y').\,

The group action is generated by vector fields

T_1^{(1)},T_2^{(1)},X_1^{(1)},X_2^{(1)},H^{(1)}

defined on the space of three variables (x,y,y′). These vector fields can be determined by the following two requirements:

  • Under the projection onto the xy-plane, they must to project to the corresponding original generators of the action T_1,T_2,X_1,X_2,H, respectively.
  • The vectors must preserve up to scale the contact structure of the jet space
\theta_1 = dy - y'dx.
Concretely, this means that the generators X(1) must satisfy
L_{X^{(1)}}\theta_1 \equiv 0 \pmod{\theta_1}
where L is the Lie derivative.

Similarly, the action of the group can be extended to the space of any number of derivatives

(x,y,y',y'',\dots,y^{(k)}).

The prolonged vector fields generating the action of the special affine group must then inductively satisfy, for each generator X ∈ {T1,T2,X1,X2,H}:

  • The projection of X(k) onto the space of variables (x,y,y′,…,y(k−1)) is X(k−1).
  • X(k) preserves the contact ideal:
L_{X^{(k)}}\theta_k \equiv 0 \pmod{\theta_1,\dots, \theta_k}
where
\theta_i = dy^{(i-1)} - y^{(i)}dx.

Carrying out the inductive construction up to order 4 gives

T_1^{(4)} = \partial_x, \quad T_2^{(4)} = \partial_y
X_1^{(4)} = x\partial_y + \partial_{y'}
\begin{align}X_2^{(4)} = y\partial_x&-y'^2\partial_{y'}-3y'y''\partial_{y''}-(3y''^2+4y'y''')\partial_{y'''}\\
&-(10y''y'''+5y'y'''')\partial_{y''''}
\end{align}
H^{(4)} = x\partial_x - y\partial_y - 2y'\partial_{y'} - 3y''\partial_{y''}-4y'''\partial_{y'''}-5y''''\partial_{y''''}.

The special affine curvature

k=\frac{1}{3}\frac{y''''}{(y'')^{5/3}}-\frac{5}{9}\frac{(y''')^2}{(y'')^{8/3}}

does not depend explicitly on x, y, or y′, and so satisfies

T_1^{(4)}k=T_2^{(4)}k=X_1^{(4)}k=0.

The vector field H acts diagonally as a modified homogeneity operator, and it is readily verified that H(4)k = 0. Finally,

X_2^{(4)}k = \frac{1}{2}[H,X 1]^{(4)}k = \frac{1}{2}[H^{(4)},X 1^{(4)}]k = 0.

The five vector fields

T_1^{(4)},T_2^{(4)},X_1^{(4)},X_2^{(4)},H^{(4)}

form an involutive distribution on (an open subset of) R6 so that, by the Frobenius integration theorem, they integrate locally to give a foliation of R6 by five-dimensional leaves. Concretely, each leaf is a local orbit of the special affine group. The function k parameterizes these leaves.

Human Motor System[edit]

Human curvilinear 2-dimensional drawing movements tend to follow the equi-affine parametrization.[1] This is more commonly known as the two thirds power law, according to which the hand's speed is proportional to the Euclidean curvature raised to the minus third power.[2] Namely,

 v = \gamma \kappa^{-\frac{1}{3}},

where v is the speed of the hand, \kappa is the Euclidean curvature and \gamma is a constant termed the velocity gain factor.

See also[edit]


References[edit]

  1. ^ Flash,Tamar; Handzel,Amir A (2007). "Affine differential geometry analysis of human arm movements". Biological cybernetics 96: 577–601. doi:10.1007/s00422-007-0145-5. 
  2. ^ Lacquaniti, Francesco and Terzuolo, Carlo and Viviani, Paolo (1983). "The law relating the kinematic and figural aspects of drawing movements". Acta psychologica 54: 115–130. doi:10.1016/0001-6918(83)90027-6.