---
title: 'MGFRFT: Multi-dimensional Graph FRFT'
url: https://www.emergentmind.com/topics/multi-dimensional-gfrft-mgfrft
type: topic
---

# MGFRFT: Multi-dimensional Graph FRFT

The Multi-dimensional Graph Fractional Fourier Transform (MGFRFT) generalizes the classical fractional Fourier transform (FRFT) to signals defined on Cartesian products of graphs, allowing independent fractional spectral analysis along each axis or domain of a multi-dimensional graph signal. This framework addresses the limitations of one-dimensional or single-parameter GFRFTs, enabling flexible adaptation to heterogeneous signal behaviors. MGFRFT encompasses both undirected and directed graphs via eigendecomposition or singular value decomposition (SVD), supports arbitrary order parameterization per dimension, and admits both closed-form and learnable filtering schemes. Notable specializations include the two-dimensional Graph Bi-Fractional Fourier Transform (2D-GBFRFT), multi-parameter GFRFTs, and geodesic-coupled variants that interpolate between standard time and graph-based fractional bases. MGFRFT has demonstrated superior denoising, compression, and restoration performance on synthetic, spatiotemporal, and real-world datasets.

## 1. Mathematical Definition and Transform Structure

Let \( G_n = (\mathcal{V}_n, \mathbf{A}_n) \) for \( n=1,\dots,M \) denote \( M \) simple, undirected weighted graphs with adjacency matrices \( \mathbf{A}_n \), each admitting the eigendecomposition
\[
\mathbf{A}_n = \boldsymbol{\chi}_n \boldsymbol{\Lambda}_n \boldsymbol{\chi}_n^T
\]
where \( \boldsymbol{\chi}_n \in \mathbb{R}^{N_n \times N_n} \) is orthonormal and \( \boldsymbol{\Lambda}_n = \mathrm{diag}(\lambda^{(n)}_1, \dots, \lambda^{(n)}_{N_n}) \). The fractional-order transform on \( G_n \) with continuous order \( \alpha_n \in \mathbb{R} \) is
\[
\mathbf{F}_{G_n}^{\alpha_n} = \boldsymbol{\chi}_n \boldsymbol{\Lambda}_n^{\alpha_n} \boldsymbol{\chi}_n^T,
\quad
\boldsymbol{\Lambda}_n^{\alpha_n} = \mathrm{diag}((\lambda_1^{(n)})^{\alpha_n}, \dots, (\lambda_{N_n}^{(n)})^{\alpha_n})
\]
For a signal \( \mathbf{x} \in \mathbb{R}^{N_1 N_2 \dots N_M} \) defined on the vertex set of a Cartesian product graph, the MGFRFT applies these transforms along each axis and stacks them via Kronecker products:
\[
\mathbf{x}_f^{(\alpha_1, \dots, \alpha_M)} = (\mathbf{F}_{G_M}^{\alpha_M} \otimes \cdots \otimes \mathbf{F}_{G_1}^{\alpha_1})\,\mathbf{x}
\]
In the 2D-GBFRFT case (\( M=2 \)), for \(\alpha_1=\alpha\), \(\alpha_2=\beta\):
\[
\mathbf{F}_{2D}^{(\alpha, \beta)} = \mathbf{F}_{G_2}^\beta \otimes \mathbf{F}_{G_1}^\alpha
\]
For an \(\mathbf{X} \in \mathbb{R}^{N_1 \times N_2}\), the transform is
\[
\mathbf{X}_f^{(\alpha, \beta)} = \mathbf{F}_{G_1}^\alpha\,\mathbf{X}\,(\mathbf{F}_{G_2}^\beta)^T
\]
The framework accommodates both Laplacian- and adjacency-based constructions [2510.11279, 2109.04358, 2603.01484].

For directed graphs, the Laplacian \( L_\ell \) is decomposed via SVD:
\[
L_\ell = U_\ell\,\Sigma_\ell\,V_\ell^T,\quad
L_\ell^\alpha = P_\ell\,R_\ell\,Q_\ell^T
\]
and multi-dimensional transforms can use Kronecker products or Kronecker sums of the left/right singular vectors [2506.03925].

## 2. Theoretical Properties

Key properties of MGFRFT and 2D-GBFRFT:

- **Identity:** \( \mathbf{F}_{2D}^{(0,0)} = \mathbf{I}_{N_1 N_2} \)
- **Reduction:** For \( \alpha = \beta \), recovers single-order 2D-GFRFT: \( \mathbf{F}_{2D}^{\alpha} \)
- **Invertibility and Index Additivity:** 
  \[
  (\mathbf{F}_{G_n}^{\alpha_n})^{-1} = \mathbf{F}_{G_n}^{-\alpha_n},\ 
  \mathbf{F}_{2D}^{(\alpha,\beta)}\,\mathbf{F}_{2D}^{(-\alpha,-\beta)} = \mathbf{I}
  \]
  and more generally, 
  \[
  \mathbf{F}_{2D}^{(\alpha,\beta)}\,\mathbf{F}_{2D}^{(\gamma,\delta)} = \mathbf{F}_{2D}^{(\alpha+\gamma,\,\beta+\delta)}
  \]
- **Unitarity / Parseval:** If graph eigenvectors are orthonormal (e.g., symmetric adjacency/Laplacian), the MGFRFT is unitary: 
  \[
  \|f\|_2^2 = \|\widehat{f}_\alpha\|_2^2
  \]
- **Linearity:** The transform is linear in the input.
- **Bandlimitedness:** MGFRFT provides a framework for defining multi-dimensional or multi-parameter bandlimitedness.

These properties guarantee well-posedness of inversion, stable energy representation, and support for efficient spectral manipulation [2510.11279, 2109.04358, 2506.03925, 2507.23570].

## 3. Filtering, Parameter Estimation, and Learning

MGFRFT domains admit specialized filtering strategies that exploit their multi-dimensional separability and parametrization:

- **Wiener-style Filtering (Grid Search):** For denoising with observation model \( \mathbf{Y} = \mathbf{G}_{G_1} \mathbf{X} \mathbf{G}_{G_2} + \mathbf{N} \), filtering proceeds by
  \[
  \hat{\mathbf{x}} = \mathbf{F}_{2D}^{(-\alpha, -\beta)}\,\mathbf{H}_{2D}\,\mathbf{F}_{2D}^{(\alpha,\beta)}\,\mathbf{y}
  \]
  and minimizing expected MSE over a grid of (\(\alpha, \beta\)), optimizing \(\mathbf{h}\) via a Wiener system [2510.11279].
- **Differentiable Joint Optimization:** Orders (\(\alpha, \beta\)) and the diagonal filter vector (\(\mathbf{h}\)) can be optimized jointly by gradient descent:
  \[
  \mathcal{L}(\alpha, \beta, \mathbf{h}) = \mathbb{E}\left\|\mathbf{F}_{2D}^{(-\alpha,-\beta)}\,\mathrm{diag}(\mathbf{h})\,\mathbf{F}_{2D}^{(\alpha,\beta)}\,\mathbf{y} - \mathbf{x}\right\|_2^2
  \]
  exploiting differentiability with respect to fractional orders [2510.11279].
- **Learnable Order Vectors:** Multi-parameter GFRFTs (MPGFRFT) generalize this by learning a vector of orders over spectral components, with closed-form gradients supporting end-to-end adaptation for compression or denoising tasks [2507.23570].
- **Hybrid Time-Vertex/Graph Interpolation:** Hybridization with the joint time-vertex FRFT (JFRFT) leverages a convex combination parameter (\(\lambda\)) to interpolate between classical time-based FRFT and graph-based FRFT bases.

These frameworks enable adaptation to both structural (graph) and temporal variability in complex signals [2510.11279, 2507.23570, 2603.01484].

## 4. Hybrid and Coupled Transform Extensions

A notable generalization is the geodesic-coupled GFRFT (GC-GFRFT), which interpolates between two unitary bases (e.g., graph-induced and classical FRFT) along the principal geodesic of the unitary group. For spatial order \(\alpha\), temporal order \(\beta\), and coupling \(\lambda \in [0,1]\), the transform is
\[
X_{\mathrm{GC}}^{(\lambda;\alpha,\beta)} = F_{G_1}^\alpha\,X\, [F_{t,\mathrm{GC}}^{(\lambda;\beta)}]^T
\]
with
\[
F_{t,\mathrm{GC}}^{(\lambda;\beta)} = F_{G_2}^\beta\,\exp(\lambda\,\log(W_t)),\quad W_t = (F_{G_2}^\beta)^H F^\beta
\]
yielding a unitary transform that at \(\lambda=0\) is the 2D-GBFRFT and at \(\lambda=1\) is the JFRFT [2603.01484, 2510.11279].

This coupling enables unified control over spectral bases, allowing the transform to adaptively select the optimal blend of graph and classical temporal analysis.

## 5. Computational Complexity and Implementation

MGFRFTs offer significant computational advantages when acting on product graphs:

- **Eigendecomposition:** \(O(\sum N_i^3)\) for the factors, versus \(O(\prod N_i^3)\) for a full product graph.
- **Forward/Inverse Transform:** Applied sequentially or via tensor contractions, yields \(O(M N \max_i N_i)\) total complexity for \(M\) factors and \(N = \prod N_i\) [2109.04358, 2506.03925].
- **Fast Implementation:** For regular graphs (e.g., paths, grids), fast transforms analogous to the FFT can be used. For large sparse graphs, Lanczos or partial EVD is effective.
- **Storage:** No need to instantiate large Kronecker matrices explicitly; transforms can be applied dimension-wise via mode-multiplication [2506.03925].

These efficiencies make MGFRFT practical for high-dimensional or large-scale graph signals.

## 6. Practical Applications and Experimental Findings

MGFRFT, and in particular 2D-GBFRFT and its variants, have been evaluated in a range of graph signal processing tasks:

- **Compression:** On spatiotemporal weather data, MGFRFT yields highly sparse representations, with only a few percent of coefficients sufficient for low relative error and high PSNR [2109.04358].
- **Denoising:** Bandlimiting in the MGFRFT domain effectively denoises temperature and air-quality datasets, outperforming classical GFRFT and single-parameter methods; energy concentration is substantially higher for multi-dimensional transforms [2506.03925].
- **Image Restoration and Deblurring:** For dynamic image deblurring (e.g., REDS’A/B datasets), 2D-GBFRFT produces lower MSE, higher PSNR, and sharper visual detail relative to 2D-GFRFT, especially around edges; hybrid transforms provide further improvement [2510.11279, 2603.01484].
- **Spectral Adaptivity:** Hybrid time-vertex/graph transforms and geodesic coupling allow selection of the spectral basis most suited to the temporal or spatial inhomogeneity of the data, with the coupling parameter tuned structurally [2510.11279, 2603.01484].
- **Compression and Encryption:** Multi-parameter MGFRFTs (Type I and II) show greatly improved performance for ultra-low-ratio compression and show key-sensitivity advantages for image encryption [2507.23570].

Consistently, MGFRFT-based approaches surpass one-order or non-fractional methods in representational efficiency, denoising, and restoration tasks across synthetic and real datasets.

## 7. Extensions and Relation to Other Transforms

MGFRFT frameworks encompass and generalize numerous previously proposed transforms:

- **Single-order GFRFT:** MGFRFT reduces to the standard GFRFT when all orders are identical.
- **Joint Time-Vertex FRFT (JFRFT):** A special case of the hybrid/interpolated scheme, acting on time–graph signals with coupled fractional orders [2510.11279].
- **Multi-parameter GFRFTs (MPGFRFT-I/II):** Permit non-constant order-vectors, further increasing adaptability at the expense of additional storage or computation, and supporting gradient-based learning schemes [2507.23570].
- **SVD-based and Kronecker-sum/product MGFRFTs:** For directed graphs, transforms constructed via SVD of fractional Laplacians extend the framework to non-symmetric structures [2506.03925].

These connections position MGFRFT as a unifying scheme for multi-dimensional and flexible spectral analysis of graph-supported signals. The framework’s adaptability and efficiency have solidified its foundational role in modern graph signal processing [2510.11279, 2109.04358, 2603.01484, 2506.03925, 2507.23570].

Source: https://www.emergentmind.com/topics/multi-dimensional-gfrft-mgfrft