Graph Fourier transform

In mathematics, the graph Fourier transform is a mathematical transform which eigendecomposes the Laplacian matrix of a graph into eigenvalues and eigenvectors. Analogously to the classical Fourier transform, the eigenvalues represent frequencies and eigenvectors form what is known as a graph Fourier basis.

The Graph Fourier transform is important in spectral graph theory. It is widely applied in the recent study of graph structured learning algorithms, such as the widely employed convolutional networks.

Definition

Given an undirected weighted graph [math]\displaystyle{ G = (V,E) }[/math], where [math]\displaystyle{ V }[/math] is the set of nodes with [math]\displaystyle{ |V| = N }[/math] ([math]\displaystyle{ N }[/math] being the number of nodes) and [math]\displaystyle{ E }[/math] is the set of edges, a graph signal [math]\displaystyle{ f: V \rightarrow \mathbb{R} }[/math] is a function defined on the vertices of the graph [math]\displaystyle{ G }[/math]. The signal [math]\displaystyle{ f }[/math] maps every vertex [math]\displaystyle{ \{v_i\}_{i=1,\ldots,N} }[/math] to a real number [math]\displaystyle{ f(i) }[/math]. Any graph signal can be projected on the eigenvectors of the Laplacian matrix [math]\displaystyle{ L }[/math]. Let [math]\displaystyle{ \lambda_l }[/math] and [math]\displaystyle{ \mu_l }[/math] be the [math]\displaystyle{ l_\text{th} }[/math] eigenvalue and eigenvector of the Laplacian matrix [math]\displaystyle{ L }[/math] (the eigenvalues are sorted in an increasing order, i.e., [math]\displaystyle{ 0= \lambda_0\leq\lambda_1\leq\cdots\leq\lambda_{N-1} }[/math]), the graph Fourier transform (GFT) [math]\displaystyle{ \hat{f} }[/math] of a graph signal [math]\displaystyle{ f }[/math] on the vertices of [math]\displaystyle{ G }[/math] is the expansion of [math]\displaystyle{ f }[/math] in terms of the eigenfunctions of [math]\displaystyle{ L }[/math]. It is defined as:

[math]\displaystyle{ \mathcal{G F}[f](\lambda_{l})=\hat{f}\left(\lambda_{l}\right)=\langle f, \mu_{l}\rangle=\sum_{i=1}^{N} f(i) \mu_{l}^*(i), }[/math]

where [math]\displaystyle{ \mu_l^* = \mu_l^\text{T} }[/math].

Since [math]\displaystyle{ L }[/math] is a real symmetric matrix, its eigenvectors [math]\displaystyle{ \{\mu_l\}_{l=0,\cdots, N-1} }[/math] form an orthogonal basis. Hence an inverse graph Fourier transform (IGFT) exists, and it is written as:

[math]\displaystyle{ \mathcal{I} \mathcal{G} \mathcal{F}[\hat{f}](i)=f(i)=\sum_{l=0}^{N-1} \hat{f}(\lambda_l) \mu_l(i) }[/math]

Analogously to the classical Fourier transform, graph Fourier transform provides a way to represent a signal in two different domains: the vertex domain and the graph spectral domain. Note that the definition of the graph Fourier transform and its inverse depend on the choice of Laplacian eigenvectors, which are not necessarily unique. The eigenvectors of the normalized Laplacian matrix are also a possible base to define the forward and inverse graph Fourier transform.

Properties

Parseval's identity

The Parseval relation holds for the graph Fourier transform, that is, for any [math]\displaystyle{ f,h\in \mathbb{R}^N }[/math]

[math]\displaystyle{ \langle f, h\rangle=\langle\hat{f}, \hat{h}\rangle. }[/math]

This gives us Parseval's identity:

[math]\displaystyle{ \sum_{i=1}^N |f(i)|^{2}=\|f\|_2^2=\langle f, f\rangle=\langle\hat{f}, \hat{f}\rangle=\|\hat{f}\|_2^2=\sum_{l=0}^{N-1} \left|\hat{f}\left(\lambda_l\right)\right|^2. }[/math]

Generalized convolution operator

The definition of convolution between two functions [math]\displaystyle{ f }[/math] and [math]\displaystyle{ g }[/math] cannot be directly applied to graph signals, because the signal translation is not defined in the context of graphs. However, by replacing the complex exponential shift in classical Fourier transform with the graph Laplacian eigenvectors, convolution of two graph signals can be defined as:

[math]\displaystyle{ (f * g)=\mathcal{I} \mathcal{G} \mathcal{F}[\mathcal{G} \mathcal{F}[f] \cdot \mathcal{G} \mathcal{F}[g]], }[/math]

[math]\displaystyle{ (f * g)(i)=\sum_{l=0}^{N-1} \hat{f}(\lambda_l) \hat{g}(\lambda_l) \mu_l(i). }[/math]

Properties of the convolution operator

The generalized convolution operator satisfies the following properties:

• Generalized convolution in the vertex domain is multiplication in the graph spectral domain: [math]\displaystyle{ \widehat{f * g}=\hat{f} \hat{g}. }[/math]
• Commutativity: [math]\displaystyle{ f * g=g * f }[/math]
• Distributivity: [math]\displaystyle{ f * (g+h)=f*g + f*h }[/math]
• Associativity: [math]\displaystyle{ (f*g)*h = f*(g*h) }[/math]
• Associativity with scalar multiplication: [math]\displaystyle{ \alpha(f * g)=(\alpha f) * g=f *(\alpha g) }[/math], for any [math]\displaystyle{ \alpha \in \mathbb{R} }[/math].
• Multiplicative identity: [math]\displaystyle{ f * g_0=f }[/math], where [math]\displaystyle{ g_0(i)=\sum_{l=0}^{N-1} \mu_{l}(i) }[/math] is an identity for the generalized convolution operator.
• The sum of the generalized convolution of two signals is a constant times the product of the sums of the two signals:

[math]\displaystyle{ \sum_{i=1}^N(f * g)(i)=\sqrt{N} \hat{f}(0) \hat{g}(0)=\frac{1}{\sqrt{N}}\left[\sum_{i=1}^N f(i)\right] \left[\sum_{i=1}^N g(i)\right]. }[/math]

Generalized translation operator

As previously stated, the classical translation operator [math]\displaystyle{ T_v }[/math] cannot be generalized to the graph setting. One way to define a generalized translation operator is through generalized convolution with a delta function centered at vertex [math]\displaystyle{ n }[/math]:[math]\displaystyle{ \left(T_{n} f\right)(i)=\sqrt{N}\left(f * \delta_{n}\right)(i){=} \sqrt{N} \sum_{l=0}^{N-1} \hat{f}\left(\lambda_{l}\right) u_{l}^{*}(n) u_{l}(i), }[/math]

where [math]\displaystyle{ \delta_i(n)=\begin{cases} 1, & \text {if } i=n, \\ 0, & \text {otherwise.} \end{cases} }[/math]

The normalization constant [math]\displaystyle{ \sqrt{N} }[/math] ensures that the translation operator preserves the signal mean, i.e.,

[math]\displaystyle{ \sum_{i=1}^N(T_n f)(i)=\sum_{i=1}^N f(i). }[/math]

Properties of the translation operator

The generalized convolution operator satisfies the following properties:

For any [math]\displaystyle{ f, g \in \mathbb{R}^N }[/math], and [math]\displaystyle{ j,k\in \{1,2,\dots,N\} }[/math],

• [math]\displaystyle{ T_j(f * g)=(T_j f) * g=f *(T_{j} g) }[/math]
• [math]\displaystyle{ T_j T_k f=T_k T_j f }[/math]
• [math]\displaystyle{ \sum_{i=1}^N (T_j f)(i)=\sqrt{N} \hat{f}(0)=\sum_{i=1}^N f(i) }[/math]
• [math]\displaystyle{ \left\|T_j f\right\| \neq \| f \| }[/math]

Applications

Image compression

Representing signals in frequency domain is a common approach to data compression. As graph signals can be sparse in their graph spectral domain, the graph Fourier transform can also be used for image compression.

Graph noise reduction

Similar to classical noise reduction of signals based on Fourier transform, graph filters based on the graph Fourier transform can be designed for graph signal denoising.

Data classification

As the graph Fourier transform enables the definition of convolution on graphs, it makes possible to adapt the conventional convolutional neural networks (CNN) to work on graphs. Graph structured semi-supervised learning algorithms such as graph convolutional network (GCN), are able to propagate the labels of a graph signal throughout the graph with a small subset of labelled nodes, theoretically operating as a first order approximation of spectral graph convolutions without computing the graph Laplacian and its eigendecomposition.

Toolbox

GSPBOX is a toolbox for signal processing of graphs, including the graph Fourier transform. It supports both Python and MATLAB languages.

External links

• DeepGraphLibrary A free Python package built for easy implementation of graph neural networks.