280 likes | 408 Views
Wavelets on manifolds. Mingzhen Tan National University of Singapore. Overview. W avelets on Euclidean spaces Continuation of the work from [Hammond, 2011] Defining wavelet bases on closed manifolds Using eigenfunctions of Laplace-Beltrami
E N D
Wavelets on manifolds Mingzhen Tan National University of Singapore
Overview • Wavelets on Euclidean spaces • Continuation of the work from [Hammond, 2011] • Defining wavelet bases on closed manifolds • Using eigenfunctions of Laplace-Beltrami • Construction of wavelet transforms of functions defined on closed manifolds • Inverse wavelet transforms • Application on Alzheimer’s disease data
Wavelets in Euclidean Spaces Admissible Function: Mother Wavelet: Wavelet coefficients: E.g. Haar Wavelets in 1-dimension E.g. Stereographic dilations for spheres
Bases on Closed Manifolds Inner Product on 2-manifolds: We consider surfaces to be discretized meshes: Decomposition of functions defined on the surface:
Wavelet Bases on Closed Manifolds Definition of wavelet bases on closed manifolds:
Wavelet Bases on Closed Manifolds ‘Fourier’ transform: We consider dilation in the Fourier domain: Inverse ‘Fourier’ transform: Wavelet coefficients in the spatial domain: Wavelet Basis:
Weight functions Weight function for wavelet bases: Weight function for scaling bases:
Localization in both frequency and spatial domains Lemma (spatial localization):
Scaling Effects Simulation with varying Simulation with varying
Eigenfunctions of Laplace-Beltrami Operator References: Meyer et. al. Discrete differential-geometry operators for triangulated 2-manifolds
Plot of the spectrum of the Laplace-Beltrami (c) (b) (a)
Wavelet Transform Wavelet coefficients, Example 1
Wavelet Transform Wavelet and scaling coefficients, Example 2
Wavelet Frames Question: Is this set of bases well behaved for representing functions on the surface? To examine this, we consider the wavelets at discretized scales as a frame, and check the frame bounds. Definition: Frames Bounds:
Inverse Wavelet Transform Reconstruction formula: Inverse Wavelet Transform of wavelet coefficients of Wavelet Transform of
Inverse Wavelet Transform Reconstructing Reconstruction formula:
Wavelet Transform (Fast Approximation) Fast Approximation Scheme using Chebyshev polynomials Chebyshevpolynomials: Chebyshev expansion: Approximation of weighting functions:
Wavelet Transform (Fast Approximation) Fast Approximation Scheme using Chebyshev polynomials Approximation of weighting functions: Approximation of wavelet coefficients:
Classification of subjects as AD or control • Subjects selection: • Source of controls: community and clinics • Patient groups with dementia were recruited from the stroke service and memory clinics in Singapore • Normal controls were defined as subjects without any cognitive complaints or functional loss & • MMSE scores of at least 23 if they had secondary/tertiary education • MMSE scores of at least 21 if they had primary/no education on initial screening and had no significant cognitive impairments on formal neuropsychological testing • AD was diagnosed in accordance with the NINCDS-ADRDA critieria • All normal and AD subjects were required to have no history of stroke, and no evidence of severe cerebrovascular disease on MRI (no infarcts) and/or presence of significant white matter lesions, defined as a grade of at least 2 in more than 4 regions • From August 2010 to November 2012, a total of 172 subjects were recruited, out of which 25 were normal controls and 20 are AD subjects.
Classification of subjects as AD or control • Experiment Objectives • A classifier constructed from a mix of hippocampi shapes could serve as an important biomarker to differentiate between the diseased and normal subjects • In particular, we want to improve the classification performance by using information on the possible hippocampal variations across the different resolutions brought about by the wavelet transform • Inputs • Jacobian determinants of transformations from a common template to the individual hippocampal shapes • Wavelet transform of the Jacobian determinants into 30 scales • Data reduction (each done separately for the Jacobian det. and their wavelet transforms) using PCA from d=1184 to d=4
Problems • Over-complete • Larger number of wavelet coefficients are used • No multi-resolution analysis (MRA)