Affine curvature

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 P₁, P₂, P₃, P₄ 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

Special affine arclength

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{aligned}\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^{2}s}{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{aligned}}

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 d²β/dt² 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

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{aligned}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{aligned}}

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

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

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{aligned}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{aligned}}

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 = x²/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 kx² + k²y² = 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|^{2}y^{2}=1.

Characterization up to affine congruence

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

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{aligned}x&\mapsto ax+by+\alpha \\y&\mapsto cx+dy+\beta ,\end{aligned}}

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⁽¹⁾ 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 ∈ {T₁,T₂,X₁,X₂,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{aligned}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{aligned}}

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⁽⁴⁾k = 0. Finally,

X_{2}^{(4)}k={\frac {1}{2}}[H,X1]^{(4)}k={\frac {1}{2}}[H^{(4)},X1^{(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) R⁶ so that, by the Frobenius integration theorem, they integrate locally to give a foliation of R⁶ 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

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.

References

Blaschke, Wilhelm (1923), Affine Differentialgeometrie, Vorlesungen über Differentialgeometrie und geometrische Grundlagen von Einsteins Relativitätstheorie, II, Berlin: Springer-Verlag OHG
Guggenheimer, Heinrich (1977), Differential Geometry, New York: Dover Publications, ISBN 978-0-486-63433-3
Shirokov, A.P. (2001a), "Affine curvature", in Hazewinkel, Michiel, Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
Shirokov, A.P. (2001b), "Affine differential geometry", in Hazewinkel, Michiel, Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
Spivak, Michael (1999), A Comprehensive introduction to differential geometry (Volume 2), Houston, TX: Publish or Perish, ISBN 978-0-914098-71-3

↑ 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.
↑ 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. PMID 6666647.

This article is issued from Wikipedia - version of the 3/1/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.

Affine curvature

Formal definition

Special affine arclength

Special affine curvature

Affine curvature

Conics

Characterization up to affine congruence

Derivation of the curvature by affine invariance

Human Motor System

See also

References