FAST EXPANSION INTO HARMONICS ON THE DISK: A STEERABLE BASIS WITH FAST RADIAL CONVOLUTIONS

Nicholas F. Marshall, Oscar Mickelin, Amit Singer

Research output: Contribution to journalArticlepeer-review

Abstract

We present a fast and numerically accurate method for expanding digitized L× L images representing functions on [-1, 1]2 supported on the disk \{x ∊ ℝ2 : |x| < 1} in the harmonics (Dirichlet Laplacian eigenfunctions) on the disk. Our method, which we refer to as the Fast Disk Harmonics Transform (FDHT), runs in O(L2 log L) operations. This basis is also known as the Fourier-Bessel basis, and it has several computational advantages: it is orthogonal, ordered by frequency, and steerable in the sense that images expanded in the basis can be rotated by applying a diagonal transform to the coefficients. Moreover, we show that convolution with radial functions can also be efficiently computed by applying a diagonal transform to the coefficients.

Original languageEnglish (US)
Pages (from-to)A2431-A2457
JournalSIAM Journal on Scientific Computing
Volume45
Issue number5
DOIs
StatePublished - Oct 2023

All Science Journal Classification (ASJC) codes

  • Computational Mathematics
  • Applied Mathematics

Keywords

  • Fourier-Bessel basis
  • Laplacian eigenfunctions
  • steerable basis

Fingerprint

Dive into the research topics of 'FAST EXPANSION INTO HARMONICS ON THE DISK: A STEERABLE BASIS WITH FAST RADIAL CONVOLUTIONS'. Together they form a unique fingerprint.

Cite this