Diffusion wavelets

Diffusion wavelets are a fast multiscale framework for the analysis of functions on discrete (or discretized continuous) structures like graphs, manifolds, and point clouds in Euclidean space. Diffusion wavelets are an extension of classical wavelet theory from harmonic analysis. Unlike classical wavelets whose basis functions are predetermined, diffusion wavelets are adapted to the geometry of a given diffusion operator ${\displaystyle T}$ (e.g., a heat kernel or a random walk). Moreover, the diffusion wavelet basis functions are constructed by dilation using the dyadic powers (powers of two) of ${\displaystyle T}$. These dyadic powers of ${\displaystyle T}$ diffusion over the space and propagate local relationships in the function throughout the space until they become global. And if the rank of higher powers of ${\displaystyle T}$ decrease (i.e., its spectrum decays), then these higher powers become compressible. From these decaying dyadic powers of ${\displaystyle T}$ comes a chain of decreasing subspaces. These subspaces are the scaling function approximation subspaces, and the differences in the subspace chain are the wavelet subspaces.

Diffusion wavelets were first introduced in 2004 by Ronald Coifman and Mauro Maggioni at Yale University. ^[1]

Algorithm

This algorithm constructs the scaling basis functions and the wavelet basis functions along with the representations of the diffusion operator ${\displaystyle T}$ at these scales.

In the algorithm below, the subscript notation ${\displaystyle \Phi _{a}}$ and ${\displaystyle \Psi _{b}}$ represents the scaling basis functions at scale ${\displaystyle a}$ and the wavelet basis functions at scale ${\displaystyle b}$ respectively. The notation ${\displaystyle [\Phi _{b}]_{\Phi _{a}}}$ denotes the matrix representation of the scaling basis ${\displaystyle \Phi _{b}}$ represented with respect to the basis ${\displaystyle \Phi _{a}}$. Lastly, the notation ${\displaystyle [T]_{\Phi _{a}}^{\Phi _{b}}}$ denotes the matrix represents of the operator ${\displaystyle T}$, where the row space of ${\displaystyle T}$ is represented with respect to the basis ${\displaystyle \Phi _{a}}$, and the column space of ${\displaystyle T}$ is represented with respect to the basis ${\displaystyle \Phi _{b}}$. Otherwise put, the domain of operator ${\displaystyle T}$ is represented with respect to the basis ${\displaystyle \Phi _{a}}$ and the range is represented with respect to the basis ${\displaystyle \Phi _{b}}$. The function ${\displaystyle QR}$ is a sparse QR decomposition with ${\displaystyle \epsilon }$ precision. ^[2]

// Input: //    ${\displaystyle T}$ is the matrix representation of the diffusion operator. //    ${\displaystyle \epsilon }$ is the precision of the QR decomposition, e.g., 1e-6. //    ${\displaystyle J}$ is the maximum number of scale levels (note: this is an optional upper bound, it may converge sooner.) // Output: //    ${\displaystyle \lbrace \Phi _{j}\rbrace }$ is the set of scaling basis functions indexed by scale ${\displaystyle j}$. //    ${\displaystyle \lbrace \Psi _{j}\rbrace }$ is the set of wavelet basis functions indexed by scale ${\displaystyle j}$.    ${\displaystyle \lbrace \Phi _{j}\rbrace ,\lbrace \Psi _{j}\rbrace \leftarrow {\text{function DiffusionWaveletTree}}(T,\epsilon ,J):}$${\displaystyle {\textbf {for}}j\leftarrow 0{\text{ to }}J-1}$:         ${\displaystyle [\Phi _{j+1}]_{\Phi _{j}},[T^{2^{j}}]_{\Phi _{j}}^{\Phi _{j+1}}\leftarrow QR\left([T^{2^{j}}]_{\Phi _{j}}^{\Phi _{j}},\epsilon \right)}$${\displaystyle [T^{2^{j+1}}]_{\Phi _{j+1}}^{\Phi _{j+1}}\leftarrow \left([T^{2^{j}}]_{\Phi _{j}}^{\Phi _{j+1}}[\Phi _{j+1}]_{\Phi _{j}}\right)^{2}}$${\displaystyle [\Psi _{j}]_{\Phi _{j}}\leftarrow QR\left(I_{\langle \Phi _{j}\rangle }-[\Phi _{j+1}]_{\Phi _{j}}\left([\Phi _{j+1}]_{\Phi _{j}}\right)^{*},\epsilon \right)}$${\displaystyle {\textbf {endfor}}}$

Applications

Mathematics

Diffusion wavelets are of general interest in mathematics. Specifically, they allow for the direct calculation of the Green′s function and the inverse graph Laplacian.

Computer science

Diffusion wavelets have been used extensively in computer science, especially in machine learning. They have been applied to the following fields:

• solving Markov decision processes and Markov chains for machine learning, ^[3]
• transfer learning, ^[4]
• value function approximation in reinforcement learning, ^[5]
• dimensionality reduction, ^[6]
• mesh compression for 3D graphics, ^[7]
• topic model analysis of document corpora. ^[8]
• relation extraction. ^[9]

See also

• Wavelets

References

1. Coifman, Ronald; Mauro Maggioni (May 2008). "Diffusion Wavelets" (PDF). Applied and Computational Harmonic Analysis. 24 (3): 329–353. Archived from the original (PDF) on 2012-04-22.
2. Maggioni, Mauro; Mahadevan, Sridhar (2006). Fast Direct Policy Evaluation using Multiscale Analysis of Markov Diffusion Processes (PDF). The 23rd International Conference on Machine Learning.
3. Mahadevan, Sridhar (2008). "Learning Representation and Control in Markov Decision Processes". Foundations and Trends in Machine Learning. 1 (4).
4. Wang, Chang; Mahadevan, Sridhar (2010). "Multiscale Manifold Alignment" (PDF). Univ. of Massachusetts Technical Report (UM-CS-2010-049).
5. Mahadevan, Sridhar; Maggioni, Mauro (2006). "Value Function Approximation using Diffusion Wavelets and Laplacian Eigenfunctions" (PDF). Advances in Neural Information Processing Systems.
6. Wang, Chang; Mahadevan, Sridhar (2009). "Multiscale Dimensionality Reduction with Diffusion Wavelets" (PDF). Univ. of Massachusetts Technical Report (UM-CS-2009-030).
7. Mahadevan, Sridhar (2007). Adaptive Mesh Compression in 3D Computer Graphics using Multiresolution Manifold Learning (PDF). The 24th International Conference on Machine Learning.
8. Wang, Chang; Mahadevan, Sridhar (2009). Multiscale Analysis of Document Corpora Based on Diffusion Models (PDF). The 21st International Joint Conference on Artificial Intelligence.^{[ permanent dead link ]}
9. Wang, Chang; James Fan; Aditya A. Kalyanpur; David Gondek (2011). Relation Extraction with Relation Topics (PDF). The 2011 Conference on Empirical Methods in Natural Language Processing.^{[ permanent dead link ]}

External links

• Mauro Maggioni's MATLAB code implementation
• Chang Wang's diffusion wavelet page