TECNOLOGIE PER LA RIABILITAZIONE. lezione # 5 la cattura del movimento. Both movement and morphology are represented. t k. t 3. t 2.

Similar documents
TECNOLOGIE PER LA RIABILITAZIONE. lezione # 4 il movimento digitale. The representation of a bone using a cloud of points. René Descartes ( )

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

5200/7200 Fall 2007 Concurrence theorems for triangles

If the center of the sphere is the origin the the equation is. x y z 2ux 2vy 2wz d 0 -(2)

Minima, Maxima, Saddle points

Vector Calculus: Understanding the Cross Product

Dipartimento di Elettronica, Informazione e Bioingegneria Robotics

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

the remaiis c p i j fi r a lent'tli of time, can havetspeclii and case any affi :er shall refuse" or neglect <to «J3,b.01Jllis p. m.

Absolute Scale Structure from Motion Using a Refractive Plate

CHICAGO, A P R I L

-t:ma. u V:. 1! *VV; jf >*4: t i f n - v i ' ONABI E :VISITd:Rs ; 8 A EASTBOURNE, SAT. 5? -tesiste- COLLEGE. D ressm aker.

Section 4.2 selected answers Math 131 Multivariate Calculus D Joyce, Spring 2014

Animation Curves and Splines 2

3.1 Iterated Partial Derivatives

Sensor fusion for motion processing and visualization

In-plane principal stress output in DIANA

Relative Positioning from Model Indexing

Partial Derivatives. Partial Derivatives. Partial Derivatives. Partial Derivatives. Partial Derivatives. Partial Derivatives

Interactive Computer Graphics. Hearn & Baker, chapter D transforms Hearn & Baker, chapter 5. Aliasing and Anti-Aliasing

Hw 4 Due Feb 22. D(fg) x y z (

DAY,.JANUARY 5, 1918 s.:i. A a T B O D B H « gg P li. SO UTH ST R E ET. BOY S A R B T A K E N FR O M T H E A G E O F E IG H T. i I..

Structure from Motion. Prof. Marco Marcon

Boolean Algebra. P1. The OR operation is closed for all x, y B x + y B

Measurement Equipment. I. Introduction

(1) Given the following system of linear equations, which depends on a parameter a R, 3x y + 5z = 2 4x + y + (a 2 14)z = a + 2

OPENRAVE TUTORIAL. Robot Autonomy Spring 2014

.ps11 INTERPRETATION OF IMAGE FLOW: A SPATIO-TEMPORAL APPROACH. Muralidhara Subbarao. Abstract

LOWELL WEEKLY JOURNAL

Lecture outline Graphics and Interaction Scan Converting Polygons and Lines. Inside or outside a polygon? Scan conversion.

CS 231. Inverse Kinematics Intro to Motion Capture. 3D characters. Representation. 1) Skeleton Origin (root) Joint centers/ bones lengths

-LOWELL WEEI^LY JOURNAL.

Factoring. Factor: Change an addition expression into a multiplication expression.

Introduction to Computer Graphics. Modeling (3) April 27, 2017 Kenshi Takayama

1.5 Equations of Lines and Planes in 3-D

21-256: Lagrange multipliers

(X 2,Y 2 ) (X 1,Y 1 ) (X 0,Y 0 ) (X c,y c ) (X 3,Y 3 )

Section 1.8. Simplifying Expressions

7.3 3-D Notes Honors Precalculus Date: Adapted from 11.1 & 11.4

f (Pijk ) V. may form the Riemann sum: . Definition. The triple integral of f over the rectangular box B is defined to f (x, y, z) dv = lim

Winter 2012 Math 255 Section 006. Problem Set 7

Graphics and Interaction Transformation geometry and homogeneous coordinates

COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates

extracted occurring from the spatial and temporal changes in an image sequence. An image sequence

Data Association for SLAM

3.3 Optimizing Functions of Several Variables 3.4 Lagrange Multipliers

You may know these...

Kinematics and dynamics analysis of micro-robot for surgical applications

Dubna 2018: lines on cubic surfaces

The Distributive Property and Expressions Understand how to use the Distributive Property to Clear Parenthesis

and A L T ( ) SOLO "It means a lot to the iner-i IS FOR FARM H O M E S v ^ r. ' 0 Iak, "" h i s ruujie which was begun in July of last ) a i,, Mr

Available online at ScienceDirect. 6th CIRP International Conference on High Performance Cutting, HPC2014

Problem Grade Total

Software correction of geometric errors for multiaxis

2 Second Derivatives. As we have seen, a function f (x, y) of two variables has four different partial derivatives: f xx. f yx. f x y.

Final Exam Review. Name: Class: Date: Short Answer

. 1. Chain rules. Directional derivative. Gradient Vector Field. Most Rapid Increase. Implicit Function Theorem, Implicit Differentiation

DOWN PLUNGE CROSS SECTIONS

(1) Tangent Lines on Surfaces, (2) Partial Derivatives, (3) Notation and Higher Order Derivatives.

[Bharathi, 3(2): February, 2014] ISSN: Impact Factor: 1.852

Backing-up Fuzzy Control of a Truck-trailer Equipped with a Kingpin Sliding Mechanism

HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008

Section 13.5: Equations of Lines and Planes. 1 Objectives. 2 Assignments. 3 Lecture Notes

Transformations. Write three rules based on what you figured out above: To reflect across the y-axis. (x,y) To reflect across y=x.

Chapter 6. Curves and Surfaces. 6.1 Graphs as Surfaces

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

A few multilinear algebraic definitions I Inner product A, B = i,j,k a ijk b ijk Frobenius norm Contracted tensor products A F = A, A C = A, B 2,3 = A

EXTRA-CREDIT PROBLEMS ON SURFACES, MULTIVARIABLE FUNCTIONS AND PARTIAL DERIVATIVES

ME5286 Robotics Spring 2014 Quiz 1 Solution. Total Points: 30

7.3 FIRST INVERSE ANTECEDENT BAF(BAF) OPERATION

GEOMETRIC TRANSFORMATIONS AND VIEWING

Kinematic Model of Robot Manipulators

= w. w u. u ; u + w. x x. z z. y y. v + w. . Remark. The formula stated above is very important in the theory of. surface integral.

LOGIC CIRCUITS. Kirti P_Didital Design 1

Stereo. CS 510 May 2 nd, 2014

13.1. Functions of Several Variables. Introduction to Functions of Several Variables. Functions of Several Variables. Objectives. Example 1 Solution

Research on error detection technology of numerical control machine tool. Cao YongJie

Visual Recognition: Image Formation

Definition 3.1 The partial derivatives of a function f(x, y) defined on an open domain containing (x, y) are denoted by f

521493S Computer Graphics Exercise 2 Solution (Chapters 4-5)

Working with Algebraic Expressions

Robot Inverse Kinematics Asanga Ratnaweera Department of Mechanical Engieering

Chapter 15: Functions of Several Variables

Ray casting. Ray casting/ray tracing

Linear First-Order PDEs

Response Surface Model Updating for Nonlinear Structures

Solution 2. ((3)(1) (2)(1), (4 3), (4)(2) (3)(3)) = (1, 1, 1) D u (f) = (6x + 2yz, 2y + 2xz, 2xy) (0,1,1) = = 4 14

CS348B Lecture 2 Pat Hanrahan, Spring Greeks: Do light rays proceed from the eye to the light, or from the light to the eye?

Math 240 Practice Problems

Jacobian Range Space

Research of Classification based on Deep Neural Network

Z+z 1 X2 Y2. or y, Graph / 4 25 jj y=±x. x2+y 2=

Motion Control (wheeled robots)

Fundamentals of Structural Geology Exercise: concepts from chapter 2

(c) 0 (d) (a) 27 (b) (e) x 2 3x2

Dynamic IBVS Control of an Underactuated UAV

CSE 167: Introduction to Computer Graphics Lecture #2: Coordinate Transformations

Sensitivity Analysis of Hopfield Neural Network in Classifying Natural RGB Color Space

Computation of the gravity gradient tensor due to topographic masses using tesseroids

1.5 Equations of Lines and Planes in 3-D

Transcription:

TENOLOGIE PER L RIILITZIONE 2005-2006 Università dei Studi di Napoi "Federico II Facotà di Ineneria orso di Laurea in Ineneria iomedica oth movement and morphooy are represented the movement requires time variant information about the pose the morphooy may be time-invariant t 3 t t t 2 1 ezione # 5 a cattura de movimento... Prof. ureio appozzo Dipartimento di Scienze de Movimento Umano e deo Sport Laboratorio di ioineneria Istituto Universitario di Scienze Motorie - Roma dismus a cattura de movimento 1 Let s separate the two probems!! dismus a cattura de movimento 2 Data coection movement data acquisition/estimation (variabes) anatomica data acquisition/estimation (parameters) Data coection movement data acquisition/estimation (variabes) t 1 t 2 t 3 t t 1 t 2 t 3 t...... dismus a cattura de movimento 3 dismus a cattura de movimento 4 1

The description of the pose The term pose audes to the ocation in space of a body The description of the pose In order to describe the pose of the bone, we substitute a compex morphooy with a simpe and time invariant morphooy (riid body hypothesis) dismus a cattura de movimento 5 dismus a cattura de movimento 6 The description of the pose Motion capture In order to describe the pose of the bone, we substitute a compex morphooy with a simpe and time invariant morphooy (riid body hypothesis) Goba system of reference (Goba frame) z For each bone invoved in the anaysis, and in each samped instant of time, motion capture provides two vectors, i.e., six scaar quantities t O t j [ t t t ]; j 1,..., xj position vector yj zj z Loca system of reference (Loca frame) dismus a cattura de movimento 7 orientation vector [ θ θ ]; j 1,..., dismus a cattura de movimento 8 j θ This is, normay, done in an indirect fashion xj yj zj 2

Estimate of the instantaneous bone pose usin Stereophotorammetry the reconstructed instantaneous position of marers Given a point in motion in the 3-D aboratory space (marer), stereophotorammetry provides its position vector (three artesian coordinates - X, Y, Z) in each samped instant of time. Minima set up: three non-ained marers [ Y position vectors pij pxij pyij ] pzij ; i 1,...,c ; j 1,..., position vector [ t j t xj j [ t yj ] t zj ; j 1,..., θ xj θ yj X Z orientation vector θzj ]; j 1,..., dismus a cattura de movimento 9 dismus a cattura de movimento 10 amera Stereophotorammetry Given a point in motion in the 3-D aboratory space (marer), stereophotorammetry provides its position vector Principa pane (three artesian coordinates - X, Y, Z) in each samped instant of time. Noda point (N) Optica axis dismus a cattura de movimento 11 dismus a cattura de movimento 12 3

amera and object point amera and object point Principa pane Principa pane Noda point (N) Noda point (N) Optica axis Optica axis Object point Object point dismus a cattura de movimento 13 dismus a cattura de movimento 14 amera, object and imae points Imae point Imae point Principa pane Noda point (N) Optica axis Object point chip D Imae pane dismus a cattura de movimento 15 dismus a cattura de movimento 16 4

Stereophotorammetry Stereophotorammetry recordin phase N 1 P N 1 P N 2 N 2 dismus a cattura de movimento 17 dismus a cattura de movimento 18 Stereophotorammetry recordin phase Stereophotorammetry recordin phase N 1 P N 1 P N 2 N 2 dismus a cattura de movimento 19 dismus a cattura de movimento 20 5

Stereophotorammetry reconstruction phase Stereophotorammetry reconstruction phase N 1 P N 1 P N 2 N 2 dismus a cattura de movimento 21 dismus a cattura de movimento 22 Stereophotorammetry reconstruction phase (with errors) Stereophotorammetry reconstruction phase N 1 P N 1 P N 2 N 2 Object point position reconstruction has been achieved by usin the foowin information: Goba camera position and orientation Loca position of the noda points caibration parameters (time invariant) Loca position of the imae points measured variabe (time variant) dismus a cattura de movimento 23 dismus a cattura de movimento 24 6

naytica stereophotorammetry Position and orientation of a camera y aibration parameters x N z Measured variabes Mathematica mode Object point oba coordinates dismus a cattura de movimento 25 dismus a cattura de movimento 26 Position and orientation of a camera Position and orientation of a camera y θ p parameters y θ p parameters x p p (3 numbers) x p p (3 numbers) N θ p (3 numbers) N θ p (3 numbers) z d z d (1 number) p p p p dismus a cattura de movimento 27 dismus a cattura de movimento 28 7

Imae coordinates naytica stereophotorammetry y data x x p x p y p N y p p p1, θ p1, d 1, p p2, q p2, d 2 z x p1, y p1, x p2, y p2 Mathematica mode X, Y, Z dismus a cattura de movimento 29 dismus a cattura de movimento 30 System caibration aibration object The contro points (marers) are ocated in nown positions. The system of reference with respect to which their ocation is iven becomes the stereophotorammetric oba frame. p p1, θ p1, d 1, p p2, q p2, d 2 Y x p1, y p1, x p2, y p2 Mathematica mode X, Y, Z X dismus a cattura de movimento 31 Z dismus a cattura de movimento 32 8

aibration object So named dynamic caibration fter a first approximation caibration usin a simpe and stationary caibration object, the marers mounted on a riid wand are traced whie movin within the measurement voume. aibration parameters are iterativey modified whie optimizin an objective function. Y ourtesy of NIH Z X dismus a cattura de movimento 33 dismus a cattura de movimento 34 So named dynamic caibration System caibration p p1, θ p1, d 1, p p2, q p2, d 2 x p1, y p1, x p2, y p2 Mathematica mode X, Y, Z ourtesy of NIH dismus a cattura de movimento 35 dismus a cattura de movimento 36 9

The experiment ctive marers In each samped instant of time p p1, θ p1, d 1, p p2, q p2, d 2 x p1, y p1, x p2, y p2 Mathematica mode X, Y, Z dismus a cattura de movimento 37 dismus a cattura de movimento 38 http://www.charndyn.com/products/products_intro.htm The movement anaysis aboratory Marer trajectory reconstruction Determination of the instantaneous bone pose usin the reconstructed instantaneous position of marers? p dismus a cattura de movimento 39 p ij [ p p p ]; i 1,..., c; j 1,..., xij position vectors yij zij orientation vector [ θ θ ]; j 1,..., dismus a cattura de movimento 40 j t j position vector [ t t t ]; j 1,..., xj θ xj yj yj zj zj 10

simpe exampe Determination of the pose of a oca set of axes the position vectors or three non-ained marers are iven dismus a cattura de movimento 41 dismus a cattura de movimento 42 Determination of the pose of a oca set of axes Determination of the pose of a oca set of axes z 1. Definition of a pane 2. Definition of two orthoona axes on that pane 3. The third axis is orthoona to the former two axes dismus a cattura de movimento 43 1. Definition of a pane 2. Definition of two orthoona axes in that pane 3. The third axis is orthoona to the former two axes dismus a cattura de movimento 44 11

Determination of the pose of a oca set of axes Determination of the instantaneous bone pose usin the reconstructed instantaneous position of marers y c ppendix z 1. Definition of a pane 2. Definition of two orthoona axes on that pane 3. The third axis is orthoona to the former two axes dismus a cattura de movimento 45 This oca frame is referred to as marer-custer technica frame dismus a cattura de movimento 46 Determination of the poses of the marer-custer technica frames y c Determination of the poses of the marer-custer technica frames y c Minima set up: three non-ained marers y c Redundant set up y c y 3 mathematica operator yc y >3 mathematica operator yc x x z z p ij [ p p p ]; i 1,..., c; j 1,..., xij position vectors yij zij cj orientation vector [ θ θ ]; j 1,..., dismus a cattura de movimento 47 t cj position vector [ t t t ]; j 1,..., cxj θ cxj cyj cyj czj czj p ij [ p p p ]; i 1,..., c; j 1,..., xij position vectors yij zij cj orientation vector [ θ θ ]; j 1,..., dismus a cattura de movimento 48 t cj position vector [ t t t ]; j 1,..., cxj θ cxj cyj cyj czj czj 12

Determination of the poses of the marer-custer technica frames y c Redundant set up y >3 mathematica operator y c yc fine dea ezione # 5 x z p ij [ p p p ]; i 1,..., c; j 1,..., xij position vectors yij zij cj [ t t t ]; j 1,..., orientation matrix R ; j 1,..., dismus a cattura de movimento 49 t position vector cxj cj cyj czj dismus a cattura de movimento 50 TENOLOGIE PER L RIILITZIONE The riid body 2005-2006 Università dei Studi di Napoi "Federico II Facotà di Ineneria orso di Laurea in Ineneria iomedica ppendix riid-body pose determination y Prof. ureio appozzo Dipartimento di Scienze de Movimento Umano e deo Sport Laboratorio di ioineneria x Istituto Universitario di Scienze Motorie - Roma dismus a cattura de movimento 51 z dismus a cattura de movimento 52 13

Point position vectors Frame unity vectors determination y p p p y p p j j ( p p ) p p x x z dismus a cattura de movimento 53 z dismus a cattura de movimento 54 Frame unity vectors determination Frame unity vectors determination y p p p j i i ( p p ) j ( p p ) j y z j i i j x x z dismus a cattura de movimento 55 z dismus a cattura de movimento 56 14

Frame orientation-matrix and position vector determination Frame orientation-matrix and direction cosines R i j y p z x j i R t i j p c z y z j x i i xx yx zx x R y z j x x x xy yy zy xy yy zy xz yz zz xz yz zz z dismus a cattura de movimento 57 dismus a cattura de movimento 58 Frame orientation-vector determination Reationship between orientation matrix and orientation vector x y z x zy zz x xy xz R x yy yz [ ] θx θy θz xa ya z a xb zayb zazb xb xayb xazb a a a a a R b x b yay b yaz [ ] b b θbx θby θbz Given two frames a and b y z j i x z dismus a cattura de movimento 59 a b tn θ θ 2sinθ θ θ 2sinθ θ θ 2sinθ a bx a by a bz 2 2 2 \ [( ) + ( ) + ( ) ] zayb ( ) zayb ( ) xazb ( ) yaxb yazb yazb zaxb xayb xaxb xazb + dismus a cattura de movimento 60 yayb zaxb + zazb 1 yaxb xayb 1 2 15

In summary In summary Given the marer position vectors in the oba frame: p p p The orientation of a riid body can therefore be described by a position vector y p z p p We were abe to estimate the position and orientation vectors of the oca frame (marer custer frame): t c c [ tcx tcy tcz] [ θ θ θ ] cx cy cz y z c or by an orientation matrix R z x In addition, the orientation matrix has been presented and its reationship with the orientation vector found. dismus a cattura de movimento 61 z x dismus a cattura de movimento 62 The end of ppendix dismus a cattura de movimento 63 16