Rotation and Orientation: Fundamentals. Perelyaev Sergei VARNA, 2011

Similar documents
Quaternions and Rotations

Quaternions and Rotations

3D Kinematics. Consists of two parts

CS770/870 Spring 2017 Quaternions

Visual Recognition: Image Formation

Lecture Note 3: Rotational Motion

Reading. Topics in Articulated Animation. Character Representation. Animation. q i. t 1 t 2. Articulated models: Character Models are rich, complex

CS612 - Algorithms in Bioinformatics

Computer Graphics. P04 Transformations. Aleksandra Pizurica Ghent University

Quaternion to Euler Angle Conversion for Arbitrary Rotation Sequence Using Geometric Methods

ME 597: AUTONOMOUS MOBILE ROBOTICS SECTION 2 COORDINATE TRANSFORMS. Prof. Steven Waslander

Quaternion Rotations AUI Course Denbigh Starkey

CS5620 Intro to Computer Graphics

12.1 Quaternions and Rotations

CS354 Computer Graphics Rotations and Quaternions

Quaternions and Rotations

Quaternions and Rotations

Quaternions & Rotation in 3D Space

Analysis of Euler Angles in a Simple Two-Axis Gimbals Set

CMSC 425: Lecture 6 Affine Transformations and Rotations

Animation. Animation

3D Transformations and Complex Representations. Computer Graphics CMU /15-662, Fall 2016

Lecture «Robot Dynamics»: Kinematics 2

Chapter 3 : Computer Animation

Animating orientation. CS 448D: Character Animation Prof. Vladlen Koltun Stanford University

Jacobian: Velocities and Static Forces 1/4

Visualizing Quaternions

Orientation & Quaternions

3D Rotations and Complex Representations. Computer Graphics CMU /15-662, Fall 2017

Jorg s Graphics Lecture Notes Coordinate Spaces 1

Transformations: 2D Transforms

3-D D Euclidean Space - Vectors

Transformation. Jane Li Assistant Professor Mechanical Engineering & Robotics Engineering

Quaternions and Euler Angles

By taking products of coordinate charts, we obtain charts for the Cartesian product of manifolds. Hence the Cartesian product is a manifold.

Rotation with Quaternions

+ i a y )( cosφ + isinφ) ( ) + i( a x. cosφ a y. = a x

Kinematics of Closed Chains

CS 475 / CS 675 Computer Graphics. Lecture 16 : Interpolation for Animation

Inertial Measurement Units II!

MTRX4700 Experimental Robotics

Straight Line motion with rigid sets

Quaternions and Exponentials

[ ] [ ] Orthogonal Transformation of Cartesian Coordinates in 2D & 3D. φ = cos 1 1/ φ = tan 1 [ 2 /1]

Rotational Joint Limits in Quaternion Space. Gino van den Bergen Dtecta

A Detailed Look into Forward and Inverse Kinematics

MEM380 Applied Autonomous Robots Winter Robot Kinematics

A GENTLE INTRODUCTION TO THE BASIC CONCEPTS OF SHAPE SPACE AND SHAPE STATISTICS

CS 445 / 645 Introduction to Computer Graphics. Lecture 21 Representing Rotations

Computer Animation II

CALCULATING TRANSFORMATIONS OF KINEMATIC CHAINS USING HOMOGENEOUS COORDINATES

Vector Addition. Qty Item Part Number 1 Force Table ME-9447B 1 Mass and Hanger Set ME Carpenter s level 1 String

Introduction to quaternions. Mathematics. Operations

Matrix Transformations. Affine Transformations

Kinematics of Wheeled Robots

Graphics and Interaction Transformation geometry and homogeneous coordinates

COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates

Rotation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J. Fowler, WesternGeco.

Polynomials. Math 4800/6080 Project Course

From isothermic triangulated surfaces to discrete holomorphicity

Why animate humans? Why is this hard? Aspects of the Problem. These lectures. Animation Apreciation 101

Anatomical Descriptions That Compute Functional Attributes

To do this the end effector of the robot must be correctly positioned relative to the work piece.

Lecture 4: Transforms. Computer Graphics CMU /15-662, Fall 2016

Geometric Algebra for Computer Graphics

Exponential Maps for Computer Vision

Rotations in 3D Graphics and the Gimbal Lock

Comparison of practically applicable mathematical descriptions of orientation and rotation in the three-dimensional Euclidean space

Today s class. Geometric objects and transformations. Informationsteknologi. Wednesday, November 7, 2007 Computer Graphics - Class 5 1

Transformations. Ed Angel Professor of Computer Science, Electrical and Computer Engineering, and Media Arts University of New Mexico

1 Historical Notes. Kinematics 5: Quaternions

Jacobian: Velocities and Static Forces 1/4

a a= a a =a a 1 =1 Division turned out to be equivalent to multiplication: a b= a b =a 1 b

Aalto CS-C3100 Computer Graphics, Fall 2016

CS-184: Computer Graphics. Today. Lecture #5: 3D Transformations and Rotations. Transformations in 3D Rotations

Quaternion properties: addition. Introduction to quaternions. Quaternion properties: multiplication. Derivation of multiplication

Curve Subdivision in SE(2)

QUATERNIONS AND ROTATIONS

Two Dimensional Viewing

ORTHOGONAL FAMILIES OF CURVES

10. Cartesian Trajectory Planning for Robot Manipulators

Computer Graphics. 2D transformations. Transforma3ons in computer graphics. Overview. Basic classes of geometric transforma3ons

Image Warping : Computational Photography Alexei Efros, CMU, Fall Some slides from Steve Seitz

Math 734 Aug 22, Differential Geometry Fall 2002, USC

Homogeneous coordinates, lines, screws and twists

Measuring Lengths The First Fundamental Form

Animation Curves and Splines 2

Linear algebra deals with matrixes: two-dimensional arrays of values. Here s a matrix: [ x + 5y + 7z 9x + 3y + 11z

Kinematic Model of Robot Manipulators

5/27/12. Objectives. Plane Curves and Parametric Equations. Sketch the graph of a curve given by a set of parametric equations.

Kinematics, Kinematics Chains CS 685

Position and Orientation Control of Robot Manipulators Using Dual Quaternion Feedback

3.3 Interpolation of rotations represented by quaternions

Animation and Quaternions

Animation. CS 4620 Lecture 32. Cornell CS4620 Fall Kavita Bala

Surface Parameterization

1. We ll look at: Types of geometrical transformation. Vector and matrix representations

PSE Game Physics. Session (3) Springs, Ropes, Linear Momentum and Rotations. Oliver Meister, Roland Wittmann

Centre for Autonomous Systems

Polar Coordinates. 2, π and ( )

Transcription:

Rotation and Orientation: Fundamentals Perelyaev Sergei VARNA, 0

What is Rotation? Not intuitive Formal definitions are also confusing Many different ways to describe Rotation (direction cosine) matri Euler angles Ais-angle Rotation vector Helical angles Unit uaternions

Orientation vs. Rotation Rotation Circular movement Orientation The state of being oriented Given a coordinate system, the orientation of an object can be represented as a rotation from a reference pose Analogy (point : vector) is similar to (orientation : rotation) Both represent a sort of (state : movement)

D Orientation π π or π θ 0 Polar Coordinates π

D Orientation π θ (t) π θ π or π 0 π Time π Although the motion is continuous, its representation could be discontinuous

D Orientation π θ (t) π θ π or π 0 π Time π Many-to-one correspondences between D orientations and their representations

Etra Parameter Y y (, y) θ X

Etra Parameter Rotation matri is yet another method of using etra parameters Y y cosθ sinθ sinθ cosθ θ (, y) X

Comple Number Imaginary iy θ Real

Comple Eponentiation Imaginary iy cosθ i sinθ e iθ iy θ Real

Rotation Composition Imaginary i( θ φ ) e e iθ e iφ θ φ Real

D Rotation Comple numbers are good for representing D orientations, but inadeuate for D rotations A comple number cannot distinguish different rotational movements that result in the same final orientation Turn 0 degree counter-clockwise Turn -40 degree clockwise Turn 480 degree counter-clockwise Imaginary θ π θ Real

D Rotation and Orientation D Rotation The conseuence of any D rotational movement can be uniuely represented by a turning angle A turning angle is independent of the choice of the reference orientation D Orientation The non-singular parameterization of D orientations reuires etra parameters Eg) Comple numbers, rotation matrices The parameterization is dependent on the choice of the reference orientation

3D Rotation Given two arbitrary orientations of a rigid object, X X Z Z Y Y How many rotations do we need to transform one orientation to the other?

3D Rotation Given two arbitrary orientations of a rigid object, vˆ θ we can always find a fied ais of rotation and a rotation angle about the ais

Euler s Rotation Theorem The general displacement of a rigid body with one point fied is a rotation about some ais Leonhard Euler (707-783) In other words, Arbitrary 3D rotation euals to one rotation around an ais Any 3D rotation leaves one vector unchanged

Euler Angles Rotation about three orthogonal aes combinations XYZ, XYX, XZY, XZX YZX, YZY, YXZ, YXY ZXY, ZXZ, ZYX, ZYZ Gimble lock Coincidence of inner most and outmost gimbles rotation aes Loss of degree of freedom

Rotation Vector vˆ vˆ : unit vector θ : scalar angle θ Rotation vector (3 parameters) v θ vˆ (, y, z) Ais-Angle ( parameters) ( θ, vˆ )

3D Orientation Unhappy with three parameters Euler angles Discontinuity (or many-to-one correspondences) Gimblelock Rotation vector (a.k.a Ais/Angle) Discontinuity (or many-to-one correspondences) π π

Using an Etra Parameter Euler parameters sin ˆ cos 3 0 θ θ v e e e e : rotation ais ˆ : rotation angle v θ

Quaternions William Rowan Hamilton (805-865) Algebraic couples (comple number) 833 iy where i

Quaternions William Rowan Hamilton (805-865) Algebraic couples (comple number) 833 Quaternions 843 iy where i w i jy kz where i j k ijk ij k, jk i, ki j ji k, kj i, ik j

Quaternions William Thomson though beautifully ingenious, have been an unmied evil to those who have touched them in any way. Arthur Cayley which contained everything but had to be unfolded into another form before it could be understood.

Unit Quaternions Unit uaternions represent 3D rotations w i jy kz 3 S ( w,, y, z) ( w, v) w y z

Rotation about an Arbitrary Ais Rotation about ais vˆ by angle (, y, z ) (, y, z) vˆ θ θ θ θ cos, v ˆ sin p p where p ( 0,, y, z) Purely Imaginary Quaternion

Unit Quaternion Algebra Identity Multiplication Inverse Unit uaternion space is closed under multiplication and inverse, but not closed under addition and subtraction ), ( ), )(, ( v v v v v v v v w w w w w w ) ) /(,,, ( ) ) /(,,, ( z y w z y w z y w z y w (,0,0,0)

Tangent Vector (Infinitesimal Rotation) 3 T S

Tangent Vector (Infinitesimal Rotation) 3 T S

Tangent Vector (Infinitesimal Rotation) 3 T I S I (, 0, 0, 0) ( 0,, y, z) Angular Velocity ω &

Rotation Vector p 3 R 3 S u p p p p p p u ( p p )

Rotation Vector 3 R 3 S I p u p p p p p p u ( p p )

Rotation Vector 3 R v log ( ) I p u p p p p p p u ( p p ) ep( v) ep ( ( log )

Rotation Vector Finite rotation Eg) Angular displacement Be careful when you add two rotation vectors e Infinitesimal rotation u e v e u v Eg) Instantaneous angular velocity Addition of angular velocity vectors are meaningful

Coordinate-Invariant Operations

Analogy (point : vector) is similar to (orientation : rotation)

Rotation Conversions In theory, conversion between any representations is possible In practice, conversion is not simple because of different conventions Quaternion to Matri 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 z y z y y z z y z y z y y z z y z y R

Method for Mapping the Four-Dimensional Space onto the Oriented Three-Dimensional Space For further presentation, we recall the notion of three-dimensional sphere S3 R4. Such a sphere defined as a subspace of the four-dimensional vector space R4 is determined by the well-known epression

The sphere S3 has the structure of a three-dimensional analytic connected closed oriented manifold, just as the three-dimensional rotation group SO(3). Therefore, such a sphere S3 can in a standard way be embedded in a four-dimensio nal arithmetic space R4 euipped with the structure of uaternion algebra. In this case, the four-dimensional vector (,, 3, 4)T whose coordinates are λ0, λ, 3 λ, 4 λ3, respectively, can be represented in the well-known algebraic form (.) of the classical Hamiltonian uaternion Λ. The sphere projection S3 RP3 associates each s uch a uaternion Λ S3 R4 with a pair of uaternions (Λ, Λ), which correspond to two identified opposite points on the surface of the three-dimensional sphere S3. If the four real parameters λ0, λ, λ, λ3 R of the classical Hamiltonian uaternions Λ R4 are used, the mapping of the sphere S3 R4 onto the space SO(3) of all possible configurations of a rigid body with a single immovable (fied) point is two-sheeted.

METHOD OF LOCAL THREE-DIMENSIONAL PARAMETRIZATION Consider the stereographic projection of the above-introduced three-dimensional sphere S3 R4 onto the oriented three-dimensional vector subspace R3 (the hyperplane Γ3 R3) in more detail. For the standard sphere S3 of unit radius r, we have the well-known relation (.6). Inturn, the sphere S3 itself as a subspace of the space R4 has the structure of an analytic connected oriented manifold, which is a submanifold of the space R4. In the case of stereographic projection (mapping) S3 R3, any point on S3 opposite to the hyperplane Γ3 R3 can be t he center of the projection. Note that, in addition, the mapping considered here is also a conformal mapping. Indeed, the stereographic projection of the sphere S3 canbe considered as part of the conformal mapping of the finite four-dimensional R4 R4 (into itself), because the stereographic projection can be continued to such a mapping. An eception is the projection center α, which corresponds to the single point at infinity in R4. Under the stereographic projection, the point at infinity of the hyperplane Γ3 R3 is associated with a single point of the sphere S3, i.e., the pole point α. Because of the above property and the fact that the mapping itself is confor mal, we use the method of the stereographic projection S3 R3. The mapping considered here associates the four co-ordinates (,, 3, 4) of a global vector R4 with the three coordinates (y, y, y3) of a local vector y R3. Usually, the operation of such projection can be written symbolic-ally as the chain S3\{α} R3. We prescribe the center of the stereographic projection α, n amely, the pole of the three-dimensional sphere S3, for which we take the chosen

CONTINUE Then the straight line passing through the given pole α(0, 0, 0, ) and an arbitrary p oint S3 on the surface of the sphere S3 intersects the oriented vector subspa-c e R3 at some point, which we denote by ϕ(). Just themapping taking such a point R4 to the oriented subspace R3 ( ϕ() R3) homeomorphism between the sphere S3 (with a single punctured p oint α) and the space R3. In this case, there eists a stereographic projection of the four-dimensional vector M3 R4 onto the oriented subspace R3. Therefore, the point of intersection of the straight line drawn from the pole α M3 through an arbitrary point R4 on the surface of the sphere S3 corresponding to the vector (,, 3, 4) with the oriented space R3 gives a single point of intersection ϕ() on the hyperplane Γ3 R3, i.e., the desired three-dimensional vector y R3. Here we present the three coordinates of this point in the form ϕ ( ) 4, 4, 3 4

Eplaining slide For the subseuent calculations, we introduce a rectangular 3 4 matri V o f the projective transformation satisfying the identities VV T V α 0 E where E is the unit 3 3 matri and α ( 0, 0, 0, )T is a 4 column vector. Under the mapping considered here, i.e., under the stereographic projection, the intersection point ϕ() R3 coincides with the desired three-dimensional v ector of local parameters y R3. Then, changing the notation ϕ() y and us ing identities (3.) and (3.), we have the coupling euation for the two vecto-rs R4 and y R3 introduced above: y V α where M3 R4 and V is the rectangular 3 4 matri of projection written as two matri ces: V E3 3 03. Thus, E. (3.3) obtained above is the point of intersection ϕ() y R3 of the straight line connecting the point α of the center (pole) of the stereographic proj ection and an arbitrary point M3 R4 on the sphere S3 itself with the oriented sub space R3. Note that E. (3.3) relating three- and four-dimensional vectors is defined for all M3 R4 ecept α. The latter can readily be proved, because the point α of the projection center (pole) does not belong to the set M3. Then, prescribing the four line ar coordinates,, 3, 4 of a point M3 R4 and using (3.3), one can readily obtain the three desired local parameters, i.e., the coordinates y, y, y3 of the point of inter T

Eplaining slide