kovac.maths

Class EllipsoidAlgorithm

java.lang.Object

kovac.maths.EllipsoidAlgorithm

public class EllipsoidAlgorithm
extends java.lang.Object

This class handles the generation of the ellipsoid fitting a set of given points. It uses Singular Value Decomposition to generate the characteristic Matrix of the ellipsoid, and center of mass to calculate its center.

Author:

bastien.kovac & pierre.weiss

Constructor Summary

Constructors

Constructor and Description
EllipsoidAlgorithm()
EllipsoidAlgorithm(java.util.List<icy.type.point.Point3D> basePoints) Builds a new EllipsoidAlgorithm from a given list of three dimensional points

Method Summary

Modifier and Type	Method and Description
static Jama.Matrix	diag(Jama.Matrix q) Return a vector column (nx1 Matrix) corresponding to the diagonal of the input Matrix.
kovac.shapes.EllipsoidOverlay	generateEllipsoid() Builds and return the ellipsoid calculated by the algorithm as an Overlay
kovac.shapes.EllipsoidOverlay	generateEllipsoid(kovac.res.quadric.QuadricExpression quadric, kovac.res.quadric.QuadricExpression quadricMicro)
kovac.res.quadric.QuadricExpression	getFinalQuadric()
static Jama.Matrix	getK(Jama.Matrix basePoints)
static Jama.Matrix	getM(Jama.Matrix K)
static Jama.Matrix	getQ0(Jama.Matrix baseMatrix)
static java.util.List<java.lang.Double>	projsplx(java.util.List<java.lang.Double> y) Original Matlab/C code from Xiaojing Ye Algorithm explained here : http://arxiv.org/abs/1101.6081 (VALIDATED)
static Jama.Matrix	proxf1(Jama.Matrix M, Jama.Matrix q) Returns the mlDivide result for the baseSphereMatrix and the input matrix
static Jama.Matrix	proxf2(Jama.Matrix q0) Function Project_On_B from matlab (VALIDATED)
static void	WriteMatrix(java.lang.String name, Jama.Matrix M)

Methods inherited from class java.lang.Object

equals, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

EllipsoidAlgorithm

public EllipsoidAlgorithm(java.util.List<icy.type.point.Point3D> basePoints)

Builds a new EllipsoidAlgorithm from a given list of three dimensional points

Parameters:

basePoints - The base points

EllipsoidAlgorithm

public EllipsoidAlgorithm()

Method Detail

getQ0

public static Jama.Matrix getQ0(Jama.Matrix baseMatrix)

WriteMatrix

public static void WriteMatrix(java.lang.String name,
                               Jama.Matrix M)

proxf1

public static Jama.Matrix proxf1(Jama.Matrix M,
                                 Jama.Matrix q)

Returns the mlDivide result for the baseSphereMatrix and the input matrix

Parameters:

q - The input matrix

m - The other input matrix

Returns:

A result matrix A as A = M\q

proxf2

public static Jama.Matrix proxf2(Jama.Matrix q0)

Function Project_On_B from matlab (VALIDATED)

Parameters:

q - Input Matrix

Returns:

a vector column

diag

public static Jama.Matrix diag(Jama.Matrix q)

Return a vector column (nx1 Matrix) corresponding to the diagonal of the input Matrix. If the input matrix is a vector column, it returns a square Matrix with these diagonal values instead (VALIDATED)

Parameters:

q - A square input matrix of n x n dimension

Returns:

The diagonal values

getK

public static Jama.Matrix getK(Jama.Matrix basePoints)

Returns:

The K matrix corresponding to the points given to the algorithm as parameters (VALIDATED)

projsplx

public static java.util.List<java.lang.Double> projsplx(java.util.List<java.lang.Double> y)

Original Matlab/C code from Xiaojing Ye Algorithm explained here : http://arxiv.org/abs/1101.6081 (VALIDATED)

Parameters:

y - An input n-dimension vector

Returns:

?

getM

public static Jama.Matrix getM(Jama.Matrix K)

getFinalQuadric

public kovac.res.quadric.QuadricExpression getFinalQuadric()

generateEllipsoid

public kovac.shapes.EllipsoidOverlay generateEllipsoid()

Builds and return the ellipsoid calculated by the algorithm as an Overlay

Returns:

An overlay displaying the ellipsoid

generateEllipsoid

public kovac.shapes.EllipsoidOverlay generateEllipsoid(kovac.res.quadric.QuadricExpression quadric,
                                                       kovac.res.quadric.QuadricExpression quadricMicro)