non-euclidean example: the unit sphere
DESCRIPTION
Non-Euclidean Example: The Unit Sphere. Differential Geometry. Formal mathematical theory Work with small ‘patches’ the ‘patches’ look Euclidean Do calculus over the patches. Manifolds. Open Sets Coordinate neighbourhood Compatible neighbourhoods. Tangents. Tangent Vectors - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/1.jpg)
Non-Euclidean Example: The Unit Sphere
![Page 2: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/2.jpg)
Differential Geometry
• Formal mathematical theory
• Work with small ‘patches’– the ‘patches’ look Euclidean
• Do calculus over the patches
![Page 3: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/3.jpg)
Manifolds
• Open Sets
• Coordinate neighbourhood
• Compatible neighbourhoods
![Page 4: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/4.jpg)
Tangents
• Tangent Vectors
• Tangent Space,
• Inner Product– Norm:
– depends on – varies smoothly
![Page 5: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/5.jpg)
Geodesics and Metrics
• The shortest path between two points is the geodesic
• The length of the geodesic is the distance between the points
![Page 6: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/6.jpg)
Exponential and Logarithm Maps
• : Maps tangents to the manifold
• : Maps points on the manifold to
•
• Both maps are locally well defined
![Page 7: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/7.jpg)
Gradient
• In Euclidean space: direction of fastest increase
• On a manifold: tangent of fastest increase
• Definition: is a real valued function. The gradient at is satisfying
directional derivative along delta
![Page 8: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/8.jpg)
The Conversion Table
1. X. Pennec, P. Fillard and N. Ayache , “A Riemannian Framework for Tensor Computing,” International. Journal of Computer Vision., 66(1), 41–66, 2006.
![Page 9: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/9.jpg)
Matrix Lie Groups
• Sets of matrices which– form a group under matrix multiplication– are Riemannian manifolds
• Examples– Rigid body transformations SE(n)– Rotations SO(n)– Affine motions A(n)
1. W. Rossman, “Lie Groups: An Introduction through Linear Groups,” Oxford University Press, 2003.
![Page 10: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/10.jpg)
Grassmann Manifolds, .
• Each point on the Grassmann manifold, , represents a dimensional subspace of .
• Numerically, represented by an orthonormal basis – matrix such that
• Representation is not unique– computation should account for this
1. A. Edelman, T. A. Arias and S. T. Smith, “The Geometry of Algorithms with Orthogonality Constraints,” SIAM Journal on Matrix Analysis and Applications, 20(2), 303–353, 1998.
![Page 11: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/11.jpg)
The Essential Manifold
• Set of matrices with– two equal and one zero singular value– let the two equal singular values be 1
• Equivalent to SO(3)xSO(3)– two-time covering of the essential manifold
• Can also be expressed as a homogeneous space
1. S. Soatto, R. Frezza and P. Perona , “Recursive Estimation on the Essential Manifold,” 3rd Europan Conference on Computer Vision, Stockholm, Sweden, May 1994, vol.II, p.61-72.
![Page 12: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/12.jpg)
Twisted Pairs
![Page 13: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/13.jpg)
The Symmetric Manifold
• contains symmetric positive definite matrices. e.g. diffusion tensor MRI– it has two different metrics
• The Affine Invariant metric• The Log-Euclidean Metric
– practically similar to the affine invariant metric– computationally easier to work with
1. V. Arsigny , P. Fillard, X. Pennec and N. Ayache , “Geometric Means in a Novel Vector Space Structure on Symmetric Positive-Definite Matrices,” SIAM Journal of Matrix Analysis and Applications, 29(1), 328–347, 2006.
![Page 14: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/14.jpg)
Mean Shift for Euclidean Spaces
• The kernel density estimate
• Mean shift as normalized gradient of
where
• The iteration
![Page 15: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/15.jpg)
Mean Shift
• Gradient ascent on kernel density– but, no line search
• Equivalent, to expectation-maximization
• Nonparametric Clustering
1. D. Comaniciu and P. Meer , “Mean Shift: A Robust Approach Towards Feature Space Analysis,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol.24, 603–619, 2002.
2. D. Comaniciu, V. Ramesh and P. Meer , “Kernel-based Object Tracking,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol.25, 564–577, 2003.
![Page 16: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/16.jpg)
Mean Shift for Manifolds
• The kernel density estimate
• Mean shift as normalized gradient of
• The iteration
![Page 17: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/17.jpg)
Mean Shift for Riemannian Manifolds
• Map points to tangent space
• Get weighted average of tangent vectors– this is the mean shift vector
• Map the mean shift vector back to the manifold
![Page 18: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/18.jpg)
Theoretical Properties
• Gradient ascent on kernel density
• Nonlinear mean shift is provably convergent– upper limit on allowed bandwidth
• nonlinear mean shift is equivalent to EM– for homogeneous spaces
![Page 19: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/19.jpg)
Motion Segmentation
• Hypothesis Generation Lie groups– Randomly pick elemental subset– Generate parameter hypothesis
• Clustering– Cluster parameters on the manifold
• Return– Number of dominant modes– Positions of dominant modes
![Page 20: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/20.jpg)
Affine Motion: Results A(2)
![Page 21: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/21.jpg)
Camera Pose Estimation SE(3)
![Page 22: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/22.jpg)
Camera Pose Estimation: Results
![Page 23: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/23.jpg)
Factorization: Results G10,3
![Page 24: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/24.jpg)
Essential Matrix: Results SO(3)xSO(3)
![Page 25: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/25.jpg)
Discontinuity Preserving Filtering
• An image is a mapping from a lattice in to data lying on a manifold
• Filtering: Run mean shift in the space
Iterations update spatial and parameter values.
• If the iteration from converges to , set in the filtered image
![Page 26: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/26.jpg)
Chromatic Noise Filtering M = G3,1
![Page 27: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/27.jpg)
Chromatic Noise Filtering
![Page 28: Non-Euclidean Example: The Unit Sphere](https://reader036.vdocument.in/reader036/viewer/2022062316/56816786550346895ddc9fcb/html5/thumbnails/28.jpg)
DTI Images R3xSym+3