Introduction to quaternions. Mathematics. Operations

Size: px
Start display at page:

Download "Introduction to quaternions. Mathematics. Operations"

Transcription

1 Introduction to quaternions Topics: Definition Mathematics Operations Euler Angles (optional) intro to quaternions 1 noel.h.hughes@gmail.com

2 Euler's Theorem y y Angle! rotation follows right hand rule x Euler Rotation Axis e = e 1 e 2 e 3 e 2 x e 3 z e 1 z intro to quaternions 2 Euler's Theorem: (paraphrased) The rotational relationship between any two coordinate frames can be described by a unit vector about which rotation takes place and a total rotation angle.

3 Quaternion Definition Discovered by William Rowan Hamilton in 1843 while walking with Lady Hamilton, crossing the Broom Bridge on the way to Dublin. He scratched the fundamental formulation i 2 = j 2 = k 2 = ijk = 1 into a stone on the side of the bridge. Q = q 1 q 2 q 3 q 4 = e sin(θ/2) cos(θ/2) = (e 1 e 2 e 3 )sin(θ/2) cos(θ/2) = q 1 i q 2 j q 3 k q 4 e 1 e 2 e 3 must be unit vector q 1 q 2 q 3 called the "vector elements" q 4 called the "scalar element" intro to quaternions 3 There is not an accepted standard on the order. In some instances, the scalar element is first, sometimes denoted q 0

4 Example 90 degree rotation about z axis e = Q = 0.0*sin(45), 0.0*sin(45), 1.0*sin(45), cos(45) intro to quaternions 4 =

5 Coordinate Frames Earth Centered Inertial ECI ECI has: Z axis along Earth rotation axis, X axis along the Vernal equinox (intersection of Equatorial and Ecliptic planes) Fixed in inertial space. Sometimes ECI is defined as aligned with ECR at a particular time intro to quaternions 5 Earth Centered Relative ECR or ECEF (Earth Centered, Earth Fixed) Mars Centered Inertial ECR has: Z axis along Earth rotation axis, X axis through intersection of Prime Meridian and Equator Rotates with Earth.

6 More Coordinate Frames Body Fixed Sensor Boresight An infinite number of coordinate frames can be defined Bottom Line: Always know what coordinate frame(s) you are working in! intro to quaternions 6

7 What Do We Do with Quaternions? ECR Describe vehicle orientation: e.g. where, rotationally, is the vehicle relative to a given reference frame? Body Fixed Transform vectors: e.g. I know where a target is in ECR, where is it relative to the body (in the body fixed frame)? intro to quaternions 7 Rotate vectors: e.g. I know a vector in my vehicle coordinate frame, where is it pointing in ECR?

8 Quaternion Mathematics Basic Functions: Normalization Multiplication Inverse, Identity Equality Successive Rotations Two Coordinate frames Vector cross and dot products Operations with Quaternions: Vector Rotation, Transformation Attitude/Orientation Propagation Vehicle maneuver/re-orientation Attitude/orientation knowledge update Appendage Pointing intro to quaternions 8

9 Quaternion Normalization To normalize quaternion, make it's magnitude = 1.0 Q n = q 1 q 2 q 3 q 4 q 1 2 q 2 2 q 3 2 q 4 2 Any quaternion that represents a rotation MUST be normal. Normalize before and after all operations. intro to quaternions 9

10 Quaternion Multiplication i i k j k j multiplication rules i 2 = j 2 = k 2 = ijk = -1 ij = k jk = i ki = j ji = -k ik = -j kj = -i Q A Q B = (q a1 i q a2 j q a3 k q a4 ) (q b1 i q b2 j q b3 k q b4 ) intro to quaternions 10 multiplication performed term by term

11 Commutivity, Associativity Q A Q B Q B Q A Not Commutive (Q A Q B ) Q C = Q A (Q B Q C ) Associative intro to quaternions 11

12 Quaternion Inverse, Identity inverse: Q -1 = Q* = -q 1 i -q 2 j -q 3 k q 4 (also called conjugate) identity: Q I = (a zero angle rotation) Multiplication by conjugate Q* Q = Q Q* = Q I = Multiplication by identity Q Q I = Q I Q = Q These are the only multiplicative operations that are commutative intro to quaternions 12 Q A Q B Q B Q A in general

13 Quaternion Equality Reversing signs of all four elements of a quaternion yields the same quaternion, mathematically Interpretation: Negative rotation about the opposite rotation axis Warning: Be very careful about software that may encounter an equivalent quaternion with reversed signs. A 1 o error in one direction may be interpreted as a 359 o error in the opposite direction IT HAS HAPPENED!!!!!!!!!!!!!!!! intro to quaternions 13

14 Successive Rotations Given: Q 1 Q 2 Q 3... Q n = sequence of rotations Q total = Q 1 Q 2 Q 3... Q n successive rotations are described by successive multiplications of the quaternions intro to quaternions 14

15 Two Coordinate Frames Given: Two coordinate frames, described by Q A Q B both referenced to the same coordinate frame To determine the quaternion going from A to B: start with: Q B = Q A Q AB multiply both sides by inverse of Q A Q A * Q B = Q A * Q A Q AB or Q AB = Q A * Q B intro to quaternions 15

16 Vector Cross, Dot Products V C V 2 V 1 Angle! Cross Product V C = V 1 X V 2 Dot Product P = V 1 V 2 = V 1 V 2 cos(θ) V C perpendicular to both V 1 and V 2 V C = V 1 V 2 sin(θ) intro to quaternions 16 If V1 and V2 are unit vectors, Θ = cos -1 (P)

17 Vector Transformation y y frame A frame B V A Q AB x z z Vector, V A, is known in coordinate frame A What is it in frame B? - V B x intro to quaternions 17 V B = Q* AB V A Q AB a zero (0) is appended to V to make it compatible with quaternion multiplication

18 Vector Transformation, Example: Vector from Star Catalog J2000 Earth Centered Inertial Star Vector, referenced to ECI V ECI = What is star vector in Spacecraft coordinate frame? V SC Spacecraft Coordinate Frame Q sc = Relative to ECI intro to quaternions 18

19 Vector Transformation, Example: Vector from Star Catalog (in ECI) to Spacecraft Coordinate Frame V SC = Q * SC V ECI Q SC from two charts back Q SC = V ECI = V SC = intro to quaternions 19

20 Vector Rotation y y frame A frame B Q AB V B x z z Vector, V B, is known in coordinate frame B What is it in frame A? - V A x intro to quaternions 20 V A = Q AB V B Q* AB a zero (0) is appended to V to make it compatible with quaternion multiplication

21 Vector Rotation Example: Star Camera Boresight J2000 Earth Centered Inertial Star Camera Boresight Star Vector Camera in boresight Body Frame In body frame V STb = V sb = What is star camera boresight in J2000 coordinate frame? V 2000 intro to quaternions 21

22 Vector Rotation Example: Star Camera Boresight from SC Body Frame to ECI V A = Q AB V B Q* AB from two charts back Q sc = V sb = V 2000 = intro to quaternions 22

23 Another View V 2000 = intro to quaternions 23

24 Attitude/Orientation Propagation We know orientation at time t 1 Q(t 1 ) and rotational rates ω What is orientation at time t 2? t 2 = t 1 + Δt Q(t 2 ) ω 2 ω = ω 1 ω 2 ω 3 0 = body rotation vector with zero (0) appended ω 3 ω 1 Calculate quaternion derivative Integrate forward in time intro to quaternions 24 Q = ½ Qω Q(t 2 ) = t 2!. Q dt + Q(t 1 ) t 1 (do not normalize)

25 Attitude/Orientation Knowledge Update Attitude/Orientation Propagation Process Inertial Measurement Unit Body Rotational Rates [ω 0.0] ½ Qω Q est Q est is "estimated" attitude Errors accumulate over time, requires update from external reference Star Camera, Sun sensor, planetary horizon sensor are common external reference Updating using quaternions: point sensor at known star compare observed with anticipated star position intro to quaternions 25 generate delta quaternion update Q est

26 Updating Attitude V o - Observed star vector V a - Anticipated star vector, based on initial Q est sensor V d = V a X V o delta rotation axis (normalize) Θ d = cos -1 (V a V o ) delta rotation angle Q d = (V d )sin(θ d /2) cos(θ d /2) delta quaternion Update Q est Q upd = Q est Q d updated attitude estimate intro to quaternions 26

27 Attitude Change/Reorientation V a V w Where we are commanded attitude = Q a Where we want to be commanded attitude = Q w rotation axis V R = V a X V w (normalize) rotation angle Θ R = cos -1 (V a V w ) rotation quaternion Q R = V R sin(θ R /2) cos(θ R /2) intro to quaternions 27 final attitude Q w = Q a Q R

28 Appendage Pointing Azimuth Azimuth Axis Gimbal Base Frame intro to quaternions 28

29 Appendage Pointing Azimuth - Elevation Azimuth Frame (after -30 deg rotation) Gimbal Base Frame Boresight 30 deg intro to quaternions 29

30 intro to quaternions 30 Appendage Pointing Example Newest Top Secret Satellite

31 Appendage Pointing Example Azimuth and Elevation Gimbal Vehicle frame Z Y Antenna Boresight is Z axis El = 30 deg Relative to Az Frame Azimuth Frame Az = -30 deg Relative to ABF X Antenna Base Frame (ABF) Pitch 90, yaw 180 Relative to vehicle frame Azimuth axis fixed in Vehicle Frame Elevation axis rotates with azimuth and Is fixed to the antenna intro to quaternions 31

32 Appendage Pointing Example Quaternions What is the orientation of the antenna relative to the universe (reference - ECI)? Vehicle attitude Q v relative to ECI Vehicle to ABF Q P90 = 0.0 sin(45) 0.0 cos(45) Q Y180 = Q VABF = Q P90 Q Y180 Azimuth Q az = sin(θ AZ /2) cos(θ AZ /2) Elevation Q El = 0.0 sin(θ EL /2) 0.0 cos(θ EL /2) Antenna boresight relative to ECI Q ANT = Q v Q VABF Q az Q El intro to quaternions 32

33 Appendage Pointing: The Problem Ø We have: Target Vector o Offboard Sensor In ECI, ECR, MCI, etc. o Onboard sensor In sensor frame, vehicle frame, etc Gimbal frame o Relative to Vehicle Many geometries Ø What we need to find Gimbal Angles To Point Payload At Target intro to quaternions 33

34 Gimbal Geometries Azimuth, Elevation Roll, Pitch Pitch, Roll Many Payload boresight may be offset Gimbal Axis Vehicle may be one or two axes of pointing Nadir, yaw for azimuth Spinner with spin axis orbit normal: despun platform, into, out of spin for azimuth elevation gimbal Single axis gimbal relative to nadir for push broom Single axis scanning gimbal Each Has Unique Math to Point intro to quaternions 34

35 Example: Find Azimuth, Elevation Gimbal Angles to Point at Target Vector Vehicle coordinate frame Z Gimbal Frame: Antenna Base Frame Az, 0.0, 0.0 Azimuth Gimbal Axis Fixed to Vehicle Target Vector Offboard sensor, => Given in ECI X Y Azimuth Gimbal Axis In ABF: V AZ = [1 0 0] Elevation Gimbal Axis Rotates with Azimuth intro to quaternions 35

36 Step 1: Transform Target Vector into Gimbal Frame (from earlier slide) Vehicle attitude Q v relative to reference (ECI) Vehicle to ABF Q P90 = 0.0 sin(45) 0.0 cos(45) Q Y180 = Q VABF = Q P90 Q Y180 ECI to ABF Q E2ABF = Q v Q VABF Target Vector in ECI = V ECI Target Vector in ABF V TABF = Q * E2ABF V ECI Q E2ABF intro to quaternions 36

37 Step 2: Find Azimuth Angle Cross Product of Target Vector with Azimuth Axis V 1 = V TABF X V AZ (normalize) Azimuth Axis V AZ = [1 0 0 ] V 1 is in ABF Y-Z Plane, Perpendicular to Target Vector V 1 V TABF intro to quaternions 37 Elevation Axis

38 Step 2: Find Azimuth Angle (cont d) Cross Product Az axis with V1 V 2 = V AZ X V 1 (normalize) V AZ V 1 is in ABF Y-Z Plane, Projection of V TABF On Y-Z plane V TABF V 1 V 2 intro to quaternions 38

39 Step 2: Find Azimuth Angle (more) Azimuth Magnitude is Angle between ABF boresight vector (Z ABF ) and V 2 Θ AZ = cos -1 (V 2 Z ABF ) (dot product) Right Hand Rule Determines Sign Of Azimuth Angle V AZ V TABF Θ AZ V 1 Z ABF = [0 0 1] V 2 intro to quaternions 39

40 Step 3: Find Elevation Angle Elevation is Angle between V TABF and V 2 Θ EL = cos -1 (V 2 V TABF ) Right Hand Rule Determines Sign Of Elevation Angle V AZ Elevation axis Rotated by Az angle Θ EL V 1 Θ AZ V 2 Z ABF intro to quaternions 40

41 Conclusion Quaternions simple, elegant to use unambiguous non-singular straightforward operations Questions? intro to quaternions 41

42 Euler Angles Three successive rotations about independent axes. Sometimes called "yaw, pitch, roll" angles or "heading, attitude and bank". Describe rotational relationship between two coordinate frames Quaternions and Euler angles contain the same information. Conversion from one to the other: one-to-one correspondence between a quaternion and a set of Euler angles intro to quaternions 42

43 Advantages, Disadvantages Euler Angles Advantages: Familiarity - nose up/down, left/right, roll left/right Disadvantages: Ambiguous - 12 rotation sequences - ypr (321), rpy (123), ypy (323),... Vector tranformations and rotations require calculating corresponding matrix or quaternion Singularities - e.g. 90 pitch causes instantaneous 180 deg change in euler angles Equations of motion are horrible intro to quaternions 43

44 Advantages, Disadvantages (cont'd) Quaternions: Advantages: No ambiguity No singularities Transformations: two quaternion multiplications Equations of motion simple Successive rotations: successive multiplications Disadvantages: Unfamiliar intro to quaternions 44 let's look at some maneuvers

45 1) pick correct Euler sequence 2) calculate matrix 3) perform vector/matrix multiplication Transformation with Euler Angles intro to quaternions 45

46 intro to quaternions 46 Propagation with Euler Angles

Quaternions and Rotations

Quaternions and Rotations CSCI 520 Computer Animation and Simulation Quaternions and Rotations Jernej Barbic University of Southern California 1 Rotations Very important in computer animation and robotics Joint angles, rigid body

More information

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

ME 597: AUTONOMOUS MOBILE ROBOTICS SECTION 2 COORDINATE TRANSFORMS. Prof. Steven Waslander ME 597: AUTONOMOUS MOILE ROOTICS SECTION 2 COORDINATE TRANSFORMS Prof. Steven Waslander OUTLINE Coordinate Frames and Transforms Rotation Matrices Euler Angles Quaternions Homogeneous Transforms 2 COORDINATE

More information

12.1 Quaternions and Rotations

12.1 Quaternions and Rotations Fall 2015 CSCI 420 Computer Graphics 12.1 Quaternions and Rotations Hao Li http://cs420.hao-li.com 1 Rotations Very important in computer animation and robotics Joint angles, rigid body orientations, camera

More information

Quaternions and Rotations

Quaternions and Rotations CSCI 420 Computer Graphics Lecture 20 and Rotations Rotations Motion Capture [Angel Ch. 3.14] Rotations Very important in computer animation and robotics Joint angles, rigid body orientations, camera parameters

More information

Quaternions & Rotation in 3D Space

Quaternions & Rotation in 3D Space Quaternions & Rotation in 3D Space 1 Overview Quaternions: definition Quaternion properties Quaternions and rotation matrices Quaternion-rotation matrices relationship Spherical linear interpolation Concluding

More information

Quaternions and Rotations

Quaternions and Rotations CSCI 480 Computer Graphics Lecture 20 and Rotations April 6, 2011 Jernej Barbic Rotations Motion Capture [Ch. 4.12] University of Southern California http://www-bcf.usc.edu/~jbarbic/cs480-s11/ 1 Rotations

More information

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

Quaternion to Euler Angle Conversion for Arbitrary Rotation Sequence Using Geometric Methods uaternion to Euler Angle Conversion for Arbitrary Rotation Sequence Using Geometric Methods ê = normalized Euler ation axis i Noel H. Hughes Nomenclature = indices of first, second and third Euler

More information

Quaternions and Rotations

Quaternions and Rotations CSCI 520 Computer Animation and Simulation Quaternions and Rotations Jernej Barbic University of Southern California 1 Rotations Very important in computer animation and robotics Joint angles, rigid body

More information

1 Historical Notes. Kinematics 5: Quaternions

1 Historical Notes. Kinematics 5: Quaternions 1 Historical Notes Quaternions were invented by the Irish mathematician William Rowan Hamilton in the late 1890s. The story goes 1 that Hamilton has pondered the problem of dividing one vector by another

More information

EE 570: Location and Navigation: Theory & Practice

EE 570: Location and Navigation: Theory & Practice EE 570: Location and Navigation: Theory & Practice Navigation Mathematics Tuesday 15 Jan 2013 NMT EE 570: Location and Navigation: Theory & Practice Slide 1 of 14 Coordinate Frames - ECI The Earth-Centered

More information

Satellite Attitude Determination

Satellite Attitude Determination Satellite Attitude Determination AERO4701 Space Engineering 3 Week 5 Last Week Looked at GPS signals and pseudorange error terms Looked at GPS positioning from pseudorange data Looked at GPS error sources,

More information

3D Kinematics. Consists of two parts

3D Kinematics. Consists of two parts D Kinematics Consists of two parts D rotation D translation The same as D D rotation is more complicated than D rotation (restricted to z-ais) Net, we will discuss the treatment for spatial (D) rotation

More information

CS354 Computer Graphics Rotations and Quaternions

CS354 Computer Graphics Rotations and Quaternions Slide Credit: Don Fussell CS354 Computer Graphics Rotations and Quaternions Qixing Huang April 4th 2018 Orientation Position and Orientation The position of an object can be represented as a translation

More information

Quaternion Rotations AUI Course Denbigh Starkey

Quaternion Rotations AUI Course Denbigh Starkey Major points of these notes: Quaternion Rotations AUI Course Denbigh Starkey. What I will and won t be doing. Definition of a quaternion and notation 3 3. Using quaternions to rotate any point around an

More information

Quaternions and Euler Angles

Quaternions and Euler Angles Quaternions and Euler Angles Revision #1 This document describes how the VR Audio Kit handles the orientation of the 3D Sound One headset using quaternions and Euler angles, and how to convert between

More information

Inertial Measurement Units II!

Inertial Measurement Units II! ! Inertial Measurement Units II! Gordon Wetzstein! Stanford University! EE 267 Virtual Reality! Lecture 10! stanford.edu/class/ee267/!! wikipedia! Polynesian Migration! Lecture Overview! short review of

More information

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

a a= a a =a a 1 =1 Division turned out to be equivalent to multiplication: a b= a b =a 1 b MATH 245 Extra Effort ( points) My assistant read through my first draft, got half a page in, and skipped to the end. So I will save you the flipping. Here is the assignment. Do just one of them. All the

More information

Transformations: 2D Transforms

Transformations: 2D Transforms 1. Translation Transformations: 2D Transforms Relocation of point WRT frame Given P = (x, y), translation T (dx, dy) Then P (x, y ) = T (dx, dy) P, where x = x + dx, y = y + dy Using matrix representation

More information

Orientation & Quaternions

Orientation & Quaternions Orientation & Quaternions Orientation Position and Orientation The position of an object can be represented as a translation of the object from the origin The orientation of an object can be represented

More information

CS612 - Algorithms in Bioinformatics

CS612 - Algorithms in Bioinformatics Fall 2017 Structural Manipulation November 22, 2017 Rapid Structural Analysis Methods Emergence of large structural databases which do not allow manual (visual) analysis and require efficient 3-D search

More information

Lecture Note 3: Rotational Motion

Lecture Note 3: Rotational Motion ECE5463: Introduction to Robotics Lecture Note 3: Rotational Motion Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio, USA Spring 2018 Lecture 3 (ECE5463

More information

Astromechanics. 12. Satellite Look Angle

Astromechanics. 12. Satellite Look Angle Astromechanics 12. Satellite Look Angle The satellite look angle refers to the angle that one would look for a satellite at a given time from a specified position on the Earth. For example, if you had

More information

Visual Recognition: Image Formation

Visual Recognition: Image Formation Visual Recognition: Image Formation Raquel Urtasun TTI Chicago Jan 5, 2012 Raquel Urtasun (TTI-C) Visual Recognition Jan 5, 2012 1 / 61 Today s lecture... Fundamentals of image formation You should know

More information

CMSC 425: Lecture 6 Affine Transformations and Rotations

CMSC 425: Lecture 6 Affine Transformations and Rotations CMSC 45: Lecture 6 Affine Transformations and Rotations Affine Transformations: So far we have been stepping through the basic elements of geometric programming. We have discussed points, vectors, and

More information

MTRX4700 Experimental Robotics

MTRX4700 Experimental Robotics MTRX 4700 : Experimental Robotics Lecture 2 Stefan B. Williams Slide 1 Course Outline Week Date Content Labs Due Dates 1 5 Mar Introduction, history & philosophy of robotics 2 12 Mar Robot kinematics &

More information

3D Transformations. CS 4620 Lecture Kavita Bala w/ prior instructor Steve Marschner. Cornell CS4620 Fall 2015 Lecture 11

3D Transformations. CS 4620 Lecture Kavita Bala w/ prior instructor Steve Marschner. Cornell CS4620 Fall 2015 Lecture 11 3D Transformations CS 4620 Lecture 11 1 Announcements A2 due tomorrow Demos on Monday Please sign up for a slot Post on piazza 2 Translation 3 Scaling 4 Rotation about z axis 5 Rotation about x axis 6

More information

CS770/870 Spring 2017 Quaternions

CS770/870 Spring 2017 Quaternions CS770/870 Spring 2017 Quaternions Primary resources used in preparing these notes: 1. van Osten, 3D Game Engine Programming: Understanding Quaternions, https://www.3dgep.com/understanding-quaternions 2.

More information

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

CS 445 / 645 Introduction to Computer Graphics. Lecture 21 Representing Rotations CS 445 / 645 Introduction to Computer Graphics Lecture 21 Representing Rotations Parameterizing Rotations Straightforward in 2D A scalar, θ, represents rotation in plane More complicated in 3D Three scalars

More information

3D Game Engine Programming. Understanding Quaternions. Helping you build your dream game engine. Posted on June 25, 2012 by Jeremiah van Oosten

3D Game Engine Programming. Understanding Quaternions. Helping you build your dream game engine. Posted on June 25, 2012 by Jeremiah van Oosten 3D Game Engine Programming Helping you build your dream game engine. Understanding Quaternions Posted on June 25, 2012 by Jeremiah van Oosten Understanding Quaternions In this article I will attempt to

More information

Transformation. Jane Li Assistant Professor Mechanical Engineering & Robotics Engineering

Transformation. Jane Li Assistant Professor Mechanical Engineering & Robotics Engineering RBE 550 MOTION PLANNING BASED ON DR. DMITRY BERENSON S RBE 550 Transformation Jane Li Assistant Professor Mechanical Engineering & Robotics Engineering http://users.wpi.edu/~zli11 Announcement Project

More information

ξ ν ecliptic sun m λ s equator

ξ ν ecliptic sun m λ s equator Attitude parameterization for GAIA L. Lindegren (1 July ) SAG LL Abstract. The GAIA attaitude may be described by four continuous functions of time, q1(t), q(t), q(t), q4(t), which form a quaternion of

More information

Geometric Transformations

Geometric Transformations Geometric Transformations CS 4620 Lecture 9 2017 Steve Marschner 1 A little quick math background Notation for sets, functions, mappings Linear and affine transformations Matrices Matrix-vector multiplication

More information

Visualizing Quaternions

Visualizing Quaternions Visualizing Quaternions Andrew J. Hanson Computer Science Department Indiana University Siggraph 25 Tutorial OUTLINE I: (45min) Twisting Belts, Rolling Balls, and Locking Gimbals: Explaining Rotation Sequences

More information

3D Transformations. CS 4620 Lecture 10. Cornell CS4620 Fall 2014 Lecture Steve Marschner (with previous instructors James/Bala)

3D Transformations. CS 4620 Lecture 10. Cornell CS4620 Fall 2014 Lecture Steve Marschner (with previous instructors James/Bala) 3D Transformations CS 4620 Lecture 10 1 Translation 2 Scaling 3 Rotation about z axis 4 Rotation about x axis 5 Rotation about y axis 6 Properties of Matrices Translations: linear part is the identity

More information

CALCULATING TRANSFORMATIONS OF KINEMATIC CHAINS USING HOMOGENEOUS COORDINATES

CALCULATING TRANSFORMATIONS OF KINEMATIC CHAINS USING HOMOGENEOUS COORDINATES CALCULATING TRANSFORMATIONS OF KINEMATIC CHAINS USING HOMOGENEOUS COORDINATES YINGYING REN Abstract. In this paper, the applications of homogeneous coordinates are discussed to obtain an efficient model

More information

Chapter 13. Vision Based Guidance. Beard & McLain, Small Unmanned Aircraft, Princeton University Press, 2012,

Chapter 13. Vision Based Guidance. Beard & McLain, Small Unmanned Aircraft, Princeton University Press, 2012, Chapter 3 Vision Based Guidance Beard & McLain, Small Unmanned Aircraft, Princeton University Press, 22, Chapter 3: Slide Architecture w/ Camera targets to track/avoid vision-based guidance waypoints status

More information

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

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

More information

Quaternions. Mike Bailey. A Useful Concept: Spherical Linear Interpolation. sin(1 t) sin t

Quaternions. Mike Bailey. A Useful Concept: Spherical Linear Interpolation. sin(1 t) sin t 1 Quaternions This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License Mike Bailey mjb@cs.oregonstate.edu Quaternions.pptx A Useful Concept: Spherical

More information

Rotation with Quaternions

Rotation with Quaternions Rotation with Quaternions Contents 1 Introduction 1.1 Translation................... 1. Rotation..................... 3 Quaternions 5 3 Rotations Represented as Quaternions 6 3.1 Dynamics....................

More information

Autonomous Navigation for Flying Robots

Autonomous Navigation for Flying Robots Computer Vision Group Prof. Daniel Cremers Autonomous Navigation for Flying Robots Lecture 3.1: 3D Geometry Jürgen Sturm Technische Universität München Points in 3D 3D point Augmented vector Homogeneous

More information

CHAPTER 2 SENSOR DATA SIMULATION: A KINEMATIC APPROACH

CHAPTER 2 SENSOR DATA SIMULATION: A KINEMATIC APPROACH 27 CHAPTER 2 SENSOR DATA SIMULATION: A KINEMATIC APPROACH 2.1 INTRODUCTION The standard technique of generating sensor data for navigation is the dynamic approach. As revealed in the literature (John Blakelock

More information

QUATERNIONS AND ROTATIONS

QUATERNIONS AND ROTATIONS 1 CHAPTER 6. QUATERNIONS AND ROTATIONS 1 INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW QUATERNIONS AND ROTATIONS 1 Quaternions and Rotations 1.1 Introduction William Rowan Hamilton

More information

Basics of Computational Geometry

Basics of Computational Geometry Basics of Computational Geometry Nadeem Mohsin October 12, 2013 1 Contents This handout covers the basic concepts of computational geometry. Rather than exhaustively covering all the algorithms, it deals

More information

Answers. Chapter 2. 1) Give the coordinates of the following points:

Answers. Chapter 2. 1) Give the coordinates of the following points: Answers Chapter 2 1) Give the coordinates of the following points: a (-2.5, 3) b (1, 2) c (2.5, 2) d (-1, 1) e (0, 0) f (2, -0.5) g (-0.5, -1.5) h (0, -2) j (-3, -2) 1 2) List the 48 different possible

More information

Computer Animation II

Computer Animation II Computer Animation II Orientation interpolation Dynamics Some slides courtesy of Leonard McMillan and Jovan Popovic Lecture 13 6.837 Fall 2002 Interpolation Review from Thursday Splines Articulated bodies

More information

Rotations in 3D Graphics and the Gimbal Lock

Rotations in 3D Graphics and the Gimbal Lock Rotations in 3D Graphics and the Gimbal Lock Valentin Koch Autodesk Inc. January 27, 2016 Valentin Koch (ADSK) IEEE Okanagan January 27, 2016 1 / 37 Presentation Road Map 1 Introduction 2 Rotation Matrices

More information

Robotics (Kinematics) Winter 1393 Bonab University

Robotics (Kinematics) Winter 1393 Bonab University Robotics () Winter 1393 Bonab University : most basic study of how mechanical systems behave Introduction Need to understand the mechanical behavior for: Design Control Both: Manipulators, Mobile Robots

More information

Vector Algebra Transformations. Lecture 4

Vector Algebra Transformations. Lecture 4 Vector Algebra Transformations Lecture 4 Cornell CS4620 Fall 2008 Lecture 4 2008 Steve Marschner 1 Geometry A part of mathematics concerned with questions of size, shape, and relative positions of figures

More information

Smartphone Video Guidance Sensor for Small Satellites

Smartphone Video Guidance Sensor for Small Satellites SSC13-I-7 Smartphone Video Guidance Sensor for Small Satellites Christopher Becker, Richard Howard, John Rakoczy NASA Marshall Space Flight Center Mail Stop EV42, Huntsville, AL 35812; 256-544-0114 christophermbecker@nasagov

More information

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

Animating orientation. CS 448D: Character Animation Prof. Vladlen Koltun Stanford University Animating orientation CS 448D: Character Animation Prof. Vladlen Koltun Stanford University Orientation in the plane θ (cos θ, sin θ) ) R θ ( x y = sin θ ( cos θ sin θ )( x y ) cos θ Refresher: Homogenous

More information

Satellite Attitude Determination II

Satellite Attitude Determination II Satellite Attitude Determination II AERO4701 Space Engineering 3 Week 7 Last Week Looked at the problem of attitude determination for satellites Examined several common methods such as inertial navigation,

More information

Jorg s Graphics Lecture Notes Coordinate Spaces 1

Jorg s Graphics Lecture Notes Coordinate Spaces 1 Jorg s Graphics Lecture Notes Coordinate Spaces Coordinate Spaces Computer Graphics: Objects are rendered in the Euclidean Plane. However, the computational space is better viewed as one of Affine Space

More information

DEVELOPMENT OF CAMERA MODEL AND GEOMETRIC CALIBRATION/VALIDATION OF XSAT IRIS IMAGERY

DEVELOPMENT OF CAMERA MODEL AND GEOMETRIC CALIBRATION/VALIDATION OF XSAT IRIS IMAGERY DEVELOPMENT OF CAMERA MODEL AND GEOMETRIC CALIBRATION/VALIDATION OF XSAT IRIS IMAGERY Leong Keong Kwoh, Xiaojing Huang, Wee Juan Tan Centre for Remote, Imaging Sensing and Processing (CRISP), National

More information

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

Analysis of Euler Angles in a Simple Two-Axis Gimbals Set Vol:5, No:9, 2 Analysis of Euler Angles in a Simple Two-Axis Gimbals Set Ma Myint Myint Aye International Science Index, Mechanical and Mechatronics Engineering Vol:5, No:9, 2 waset.org/publication/358

More information

Rotation and Orientation: Fundamentals. Perelyaev Sergei VARNA, 2011

Rotation and Orientation: Fundamentals. Perelyaev Sergei VARNA, 2011 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

More information

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

To do this the end effector of the robot must be correctly positioned relative to the work piece. Spatial Descriptions and Transformations typical robotic task is to grasp a work piece supplied by a conveyer belt or similar mechanism in an automated manufacturing environment, transfer it to a new position

More information

GYRO MISALIGNMENT AND SCALE FACTOR ERROR DETERMINATION IN MARS ORBITER MISSION

GYRO MISALIGNMENT AND SCALE FACTOR ERROR DETERMINATION IN MARS ORBITER MISSION IAA-AAS-DCoSS2-14-07-11 GYRO MISALIGNMENT AND SCALE FACTOR ERROR DETERMINATION IN MARS ORBITER MISSION Naga Manjusha, * M. Srikanth, Ritu Karidhal, $ and V. Kesava Raju INTRODUCTION This paper deals with

More information

Quaternions and Dual Coupled Orthogonal Rotations in Four-Space

Quaternions and Dual Coupled Orthogonal Rotations in Four-Space Quaternions and Dual Coupled Orthogonal Rotations in Four-Space Kurt Nalty January 8, 204 Abstract Quaternion multiplication causes tensor stretching) and versor turning) operations. Multiplying by unit

More information

LOCAL GEODETIC HORIZON COORDINATES

LOCAL GEODETIC HORIZON COORDINATES LOCAL GEODETIC HOIZON COODINATES In many surveying applications it is necessary to convert geocentric Cartesian coordinates X,,Z to local geodetic horizon Cartesian coordinates E,N,U (East,North,Up). Figure

More information

3D Rotation: more than just a hobby

3D Rotation: more than just a hobby 3D Rotation: more than just a hobby Eric Yang s special talent is the mental rotation of threedimensional objects. Mental rotation is Yang s hobby. I can see textures and imperfections and the play of

More information

Visualizing Quaternion Multiplication

Visualizing Quaternion Multiplication Received April 12, 2017, accepted May 8, 2017, date of publication May 17, 2017, date of current version June 28, 2017. Digital Object Identifier 10.1109/ACCESS.2017.2705196 Visualizing Quaternion Multiplication

More information

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

Quaternion properties: addition. Introduction to quaternions. Quaternion properties: multiplication. Derivation of multiplication Introduction to quaternions Definition: A quaternion q consists of a scalar part s, s, and a vector part v ( xyz,,, v 3 : q where, [ s, v q [ s, ( xyz,, q s+ ix + jy + kz i 2 j 2 k 2 1 ij ji k k Quaternion

More information

Geometric transformations in 3D and coordinate frames. Computer Graphics CSE 167 Lecture 3

Geometric transformations in 3D and coordinate frames. Computer Graphics CSE 167 Lecture 3 Geometric transformations in 3D and coordinate frames Computer Graphics CSE 167 Lecture 3 CSE 167: Computer Graphics 3D points as vectors Geometric transformations in 3D Coordinate frames CSE 167, Winter

More information

Mathematics in Orbit

Mathematics in Orbit Mathematics in Orbit Dan Kalman American University Slides and refs at www.dankalman.net Outline Basics: 3D geospacial models Keyhole Problem: Related Rates! GPS: space-time triangulation Sensor Diagnosis:

More information

1 Transformations. Chapter 1. Transformations. Department of Computer Science and Engineering 1-1

1 Transformations. Chapter 1. Transformations. Department of Computer Science and Engineering 1-1 Transformations 1-1 Transformations are used within the entire viewing pipeline: Projection from world to view coordinate system View modifications: Panning Zooming Rotation 1-2 Transformations can also

More information

Homogeneous Coordinates

Homogeneous Coordinates COMS W4172 3D Math 2 Steven Feiner Department of Computer Science Columbia Universit New York, NY 127 www.cs.columbia.edu/graphics/courses/csw4172 Februar 1, 218 1 Homogeneous Coordinates w X W Y X W Y

More information

METR 4202: Advanced Control & Robotics

METR 4202: Advanced Control & Robotics Position & Orientation & State t home with Homogenous Transformations METR 4202: dvanced Control & Robotics Drs Surya Singh, Paul Pounds, and Hanna Kurniawati Lecture # 2 July 30, 2012 metr4202@itee.uq.edu.au

More information

CRYOSAT-2: STAR TRACKER ATTITUDE CALCULATION TOOL FILE TRANSFER DOCUMENT

CRYOSAT-2: STAR TRACKER ATTITUDE CALCULATION TOOL FILE TRANSFER DOCUMENT Page: 1 / 7 CRYOSAT-: STAR TRACKER ATTITUDE CALCULATION TOOL FILE TRANSFER DOCUMENT 1. INTRODUCTION This is the File Transfer Document for the executable routine that computes attitude data (quaternions,

More information

Animation and Quaternions

Animation and Quaternions Animation and Quaternions Partially based on slides by Justin Solomon: http://graphics.stanford.edu/courses/cs148-1-summer/assets/lecture_slides/lecture14_animation_techniques.pdf 1 Luxo Jr. Pixar 1986

More information

CS184: Using Quaternions to Represent Rotation

CS184: Using Quaternions to Represent Rotation Page 1 of 5 CS 184 home page A note on these notes: These notes on quaternions were created as a resource for students taking CS184 at UC Berkeley. I am not doing any research related to quaternions and

More information

( ) Derivation of Polar Reduction Formula for a Calculator Robert Bernecky April, 2018 ( )

( ) Derivation of Polar Reduction Formula for a Calculator Robert Bernecky April, 2018 ( ) Derivation of Polar Reduction Formula for a Calculator Robert Bernecky April, 2018 1 Problem Statement The polar reduction formula takes an observer's assumed position (lat, lon), and a body's celestial

More information

Attitude Representation

Attitude Representation Attitude Representation Basilio Bona DAUIN Politecnico di Torino Semester 1, 015-16 B. Bona (DAUIN) Attitude Representation Semester 1, 015-16 1 / 3 Mathematical preliminaries A 3D rotation matrix R =

More information

θ x Week Date Lecture (M: 2:05p-3:50, 50-N202) 1 23-Jul Introduction + Representing Position & Orientation & State 2 30-Jul

θ x Week Date Lecture (M: 2:05p-3:50, 50-N202) 1 23-Jul Introduction + Representing Position & Orientation & State 2 30-Jul θ x 2018 School of Information Technology and Electrical Engineering at the University of Queensland Lecture Schedule Week Date Lecture (M: 2:05p-3:50, 50-N202) 1 23-Jul Introduction + Representing Position

More information

Quaternions & Octonions Made Easy. A Series Illustrating Innovative Forms of the Organization & Exposition of Mathematics by Walter Gottschalk

Quaternions & Octonions Made Easy. A Series Illustrating Innovative Forms of the Organization & Exposition of Mathematics by Walter Gottschalk Quaternions & Octonions Made Easy #20 of Gottschalk's Gestalts A Series Illustrating Innovative Forms of the Organization & Exposition of Mathematics by Walter Gottschalk Infinite Vistas Press PVD RI 2001

More information

Transforms. COMP 575/770 Spring 2013

Transforms. COMP 575/770 Spring 2013 Transforms COMP 575/770 Spring 2013 Transforming Geometry Given any set of points S Could be a 2D shape, a 3D object A transform is a function T that modifies all points in S: T S S T v v S Different transforms

More information

An Observatory for Ocean, Climate and Environment. SAC-D/Aquarius. MWR Geometric Correction Algorithms. Felipe Madero

An Observatory for Ocean, Climate and Environment. SAC-D/Aquarius. MWR Geometric Correction Algorithms. Felipe Madero An Observatory for Ocean, Climate and Environment SAC-D/Aquarius MWR Geometric Correction Algorithms Felipe Madero 7th Aquarius SAC-D Science Meeting Buenos Aires 1 April 11-13, 2012 Geometric Processing

More information

Homogeneous coordinates, lines, screws and twists

Homogeneous coordinates, lines, screws and twists Homogeneous coordinates, lines, screws and twists In lecture 1 of module 2, a brief mention was made of homogeneous coordinates, lines in R 3, screws and twists to describe the general motion of a rigid

More information

BASIC ELEMENTS. Geometry is the study of the relationships among objects in an n-dimensional space

BASIC ELEMENTS. Geometry is the study of the relationships among objects in an n-dimensional space GEOMETRY 1 OBJECTIVES Introduce the elements of geometry Scalars Vectors Points Look at the mathematical operations among them Define basic primitives Line segments Polygons Look at some uses for these

More information

Robot Inverse Kinematics Asanga Ratnaweera Department of Mechanical Engieering

Robot Inverse Kinematics Asanga Ratnaweera Department of Mechanical Engieering PR 5 Robot Dynamics & Control /8/7 PR 5: Robot Dynamics & Control Robot Inverse Kinematics Asanga Ratnaweera Department of Mechanical Engieering The Inverse Kinematics The determination of all possible

More information

Working With 3D Rotations. Stan Melax Graphics Software Engineer, Intel

Working With 3D Rotations. Stan Melax Graphics Software Engineer, Intel Working With 3D Rotations Stan Melax Graphics Software Engineer, Intel Human Brain is wired for Spatial Computation Which shape is the same: a) b) c) I don t need to ask for directions Translations A childhood

More information

Motion Control (wheeled robots)

Motion Control (wheeled robots) Motion Control (wheeled robots) Requirements for Motion Control Kinematic / dynamic model of the robot Model of the interaction between the wheel and the ground Definition of required motion -> speed control,

More information

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

Transformations. Ed Angel Professor of Computer Science, Electrical and Computer Engineering, and Media Arts University of New Mexico Transformations Ed Angel Professor of Computer Science, Electrical and Computer Engineering, and Media Arts University of New Mexico Angel: Interactive Computer Graphics 4E Addison-Wesley 25 1 Objectives

More information

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

Rotation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J. Fowler, WesternGeco. otation parameters for model building and stable parameter inversion in orthorhombic media Cintia Lapilli* and Paul J Fowler, WesternGeco Summary Symmetry groups commonly used to describe seismic anisotropy

More information

JUST THE MATHS SLIDES NUMBER 5.2. GEOMETRY 2 (The straight line) A.J.Hobson

JUST THE MATHS SLIDES NUMBER 5.2. GEOMETRY 2 (The straight line) A.J.Hobson JUST THE MATHS SLIDES NUMBER 5.2 GEOMETRY 2 (The straight line) by A.J.Hobson 5.2.1 Preamble 5.2.2 Standard equations of a straight line 5.2.3 Perpendicular straight lines 5.2.4 Change of origin UNIT 5.2

More information

Zero Robotics Autonomous Space Capture Challenge Manual

Zero Robotics Autonomous Space Capture Challenge Manual Zero Robotics Autonomous Space Capture Challenge Manual v1.3 1 Introduction 1.1 Conventions Vectors All vectors in this document are denoted with a bold face font. Of special note is the position vector

More information

Use of n-vector for Radar Applications

Use of n-vector for Radar Applications Use of n-vector for Radar Applications Nina Ødegaard, Kenneth Gade Norwegian Defence Research Establishment Kjeller, NORWAY email: Nina.Odegaard@ffi.no Kenneth.Gade@ffi.no Abstract: This paper aims to

More information

Quaternions: From Classical Mechanics to Computer Graphics, and Beyond

Quaternions: From Classical Mechanics to Computer Graphics, and Beyond Proceedings of the 7 th Asian Technology Conference in Mathematics 2002. Invited Paper Quaternions: From Classical Mechanics to Computer Graphics, and Beyond R. Mukundan Department of Computer Science

More information

Mobile Robotics. Mathematics, Models, and Methods. HI Cambridge. Alonzo Kelly. Carnegie Mellon University UNIVERSITY PRESS

Mobile Robotics. Mathematics, Models, and Methods. HI Cambridge. Alonzo Kelly. Carnegie Mellon University UNIVERSITY PRESS Mobile Robotics Mathematics, Models, and Methods Alonzo Kelly Carnegie Mellon University HI Cambridge UNIVERSITY PRESS Contents Preface page xiii 1 Introduction 1 1.1 Applications of Mobile Robots 2 1.2

More information

Perspective Projection Describes Image Formation Berthold K.P. Horn

Perspective Projection Describes Image Formation Berthold K.P. Horn Perspective Projection Describes Image Formation Berthold K.P. Horn Wheel Alignment: Camber, Caster, Toe-In, SAI, Camber: angle between axle and horizontal plane. Toe: angle between projection of axle

More information

Computer Graphics. P04 Transformations. Aleksandra Pizurica Ghent University

Computer Graphics. P04 Transformations. Aleksandra Pizurica Ghent University Computer Graphics P4 Transformations Aleksandra Pizurica Ghent Universit Telecommunications and Information Processing Image Processing and Interpretation Group Transformations in computer graphics Goal:

More information

Math for Geometric Optics

Math for Geometric Optics Algebra skills Math for Geometric Optics general rules some common types of equations units problems with several variables (substitution) Geometry basics Trigonometry Pythagorean theorem definitions,

More information

Object Representation Affine Transforms. Polygonal Representation. Polygonal Representation. Polygonal Representation of Objects

Object Representation Affine Transforms. Polygonal Representation. Polygonal Representation. Polygonal Representation of Objects Object Representation Affine Transforms Polygonal Representation of Objects Although perceivable the simplest form of representation they can also be the most problematic. To represent an object polygonally,

More information

Lab 2A Finding Position and Interpolation with Quaternions

Lab 2A Finding Position and Interpolation with Quaternions Lab 2A Finding Position and Interpolation with Quaternions In this Lab we will learn how to use the RVIZ Robot Simulator, Python Programming Interpreter and ROS tf library to study Quaternion math. There

More information

Animation. Animation

Animation. Animation CS475m - Computer Graphics Lecture 5 : Interpolation for Selected (key) frames are specified. Interpolation of intermediate frames. Simple and popular approach. May give incorrect (inconsistent) results.

More information

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

3D Transformations and Complex Representations. Computer Graphics CMU /15-662, Fall 2016 3D Transformations and Complex Representations Computer Graphics CMU 15-462/15-662, Fall 2016 Quiz 4: Trees and Transformations Student solutions (beautiful!): Rotations in 3D What is a rotation, intuitively?

More information

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

+ i a y )( cosφ + isinφ) ( ) + i( a x. cosφ a y. = a x Rotation Matrices and Rotated Coordinate Systems Robert Bernecky April, 2018 Rotated Coordinate Systems is a confusing topic, and there is no one standard or approach 1. The following provides a simplified

More information

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

Vector Addition. Qty Item Part Number 1 Force Table ME-9447B 1 Mass and Hanger Set ME Carpenter s level 1 String rev 05/2018 Vector Addition Equipment List Qty Item Part Number 1 Force Table ME-9447B 1 Mass and Hanger Set ME-8979 1 Carpenter s level 1 String Purpose The purpose of this lab is for the student to gain

More information

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

CS 475 / CS 675 Computer Graphics. Lecture 16 : Interpolation for Animation CS 475 / CS 675 Computer Graphics Lecture 16 : Interpolation for Keyframing Selected (key) frames are specified. Interpolation of intermediate frames. Simple and popular approach. May give incorrect (inconsistent)

More information

Visualizing Quaternions

Visualizing Quaternions Visualizing Quaternions Andrew J. Hanson Computer Science Department Indiana University Siggraph 1 Tutorial 1 GRAND PLAN I: Fundamentals of Quaternions II: Visualizing Quaternion Geometry III: Quaternion

More information

EECE 478. Learning Objectives. Learning Objectives. Linear Algebra and 3D Geometry. Linear algebra in 3D. Coordinate systems

EECE 478. Learning Objectives. Learning Objectives. Linear Algebra and 3D Geometry. Linear algebra in 3D. Coordinate systems EECE 478 Linear Algebra and 3D Geometry Learning Objectives Linear algebra in 3D Define scalars, points, vectors, lines, planes Manipulate to test geometric properties Coordinate systems Use homogeneous

More information