Title: Nonlinear LeastSquares
1Non-linear Least-Squares
2Why non-linear?
- The model is non-linear (e.g. joints, position,
..) - The error function is non-linear
3Setup
- Let f be a function such that
- where x is a vector of parameters
4Setup
- Let f be a function such that
- where x is a vector of parameters
- Let ak,bk be a set of measurements/constraints.
- We fit f to the data by solving
5Setup
- Let f be a function such that
- where x is a vector of parameters
- Let ak,bk be a set of measurements/constraints.
- We fit f to the data by solving
6Example
7Overview
- Existence and uniqueness of minimum
- Steepest-descent
- Newtons method
- Gauss-Newtons method
- Levenberg-Marquardt method
8A non-linear functionthe Rosenbrock function
Global minimum at (1,1)
9Existence of minimum
- A local minima is characterized by
-
-
10Existence of minimum
11Descent algorithm
- Start at an initial position x0
- Until convergence
- Find minimizing step dxk
- xk1xk dxk
- Produce a sequence x0, x1, , xn such that
-
f(x0) gt f(x1) gt gt f(xn)
12Descent algorithm
- Start at an initial position x0
- Until convergence
- Find minimizing step dxk
- using a local approximation of f
- xk1xk dxk
- Produce a sequence x0, x1, , xn such that
-
f(x0) gt f(x1) gt gt f(xn)
13Approximation using Taylor series
14Approximation using Taylor series
15Approximation using Taylor series
16Steepest descent
- Step
- where is chosen such that
- using a line search algorithm
17(No Transcript)
18(No Transcript)
19In the plane of the steepest descent direction
20Rosenbrock function (1000 iterations)
21Newtons method
- Step second order approximation
- At the minimum of N
22(No Transcript)
23Problem
- If is not positive semi-definite, then
- is not a descent direction the step increases
the error - function
- Uses positive semi-definite approximation of
Hessian based on the jacobian - (quasi-Newton methods)
24Gauss-Newton method
- Step use
- with the approximate hessian
- Advantages
- No second order derivatives
- is positive semi-definite
25Rosenbrock function (48 evaluations)
26Levenberg-Marquardt algorithm
- Blends Steepest descent and Gauss-Newton
- At each step solve, for the descent direction
27Managing the damping parameter
- General approach
- If step fails, increase damping until step is
successful - If step succeeds, decrease damping to take larger
step - Improved damping
28Rosenbrock function (90 evaluations)
29Resynthesizing Facial Animation through 3D
Model-Based Tracking, Pighin etal.,
??
??
??
??
??
Video frame
Model
30Model fitting
??
2
??
-
??
??
??
Video frame
Model
31Model fitting
??
2
??
-
??
??
??
Video frame
Fitted model
Model
32Resynthesizing Facial Animation through 3D
Model-Based Tracking, Pighin etal.,
- Facial tracking by fitting face model
33Resynthesizing Facial Animation through 3D
Model-Based Tracking, Pighin etal.,
- Facial tracking by fitting face model
- Problem
- Given an image I
- A 3D face model
34Resynthesizing Facial Animation through 3D
Model-Based Tracking, Pighin etal.,
- Facial tracking by fitting face model
- Problem
- Given an image I
- A 3D face model
- Minimize
- Where is the image produced by the face
model - And is a regularization term.
35Problem setup
36Solving
- Using a gaussian pyramid
- of the input image
- For level in coarserLevel finerLever
- Initialize using solution at level-1
- Minimize at level using Levenberg-Marquardt
- Solution is solution at finerLevel
37Tracking results
38Conclusion
- Get a good initial guess
- Prediction (temporal/spatial coherency)
- Partial solve
- Prior knowledge
- Levenberg-Marquardt is your friend