Abstract:This paper develops a constructive numerical scheme for Fourier-Bessel approximations on disks compatible with convolutions supported on disks. We address accurate finite Fourier-Bessel transforms (FFBT) and inverse finite Fourier-Bessel transforms (iFFBT) of functions on disks using the discrete Fourier Transform (DFT) on Cartesian grids. Whereas the DFT and its fast implementation (FFT) are ubiquitous and are powerful for computing convolutions, they are not exactly steerable under rotations. In contrast, Fourier-Bessel expansions are steerable, but lose both this property and the preservation of band limits under convolution. This work captures the best features of both as the band limit is allowed to increase. The convergence/error analysis and asymptotic steerability of FFBT/ iFFBT are investigated. Conditions are established for the FFBT to converge to the Fourier-Bessel coefficient and for the iFFBT to uniformly approximate the Fourier-Bessel partial sums. The matrix form of the finite transforms is discussed. The implementation of the discrete method to compute numerical approximation of convolutions of compactly supported functions on disks is considered as well.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop a constructive numerical scheme to achieve the Fourier - Bessel approximation of functions on the disk and ensure that these approximations are compatible with the convolution operation on the disk. Specifically, the authors focus on how to accurately calculate the finite Fourier - Bessel transform (FFBT) and its inverse transform (iFFBT), and study their performance in numerical implementation. ### Specific description of the main problem 1. **Numerical implementation of the Fourier - Bessel transform**: - Although the traditional discrete Fourier transform (DFT) and fast Fourier transform (FFT) are widely used in convolution calculations, they are not exactly steerable under rotation. - In contrast, the Fourier - Bessel expansion is steerable, but it will lose the band - limited property and steerability under convolution. - This paper aims to combine the advantages of both and develop a numerical method that has steerability and can maintain good performance in convolution while allowing the band - limit to increase. 2. **Convergence and error analysis**: - Study the convergence of the finite Fourier - Bessel transform (FFBT) and the inverse finite Fourier - Bessel transform (iFFBT). - Establish conditions so that FFBT can converge to the Fourier - Bessel coefficients and iFFBT can uniformly approximate the Fourier - Bessel partial sum. 3. **Matrix form and numerical experiments**: - Explore the matrix form of the finite transform and use MATLAB to conduct numerical experiments to verify the effectiveness of the proposed numerical method. - In particular, study how to implement FFBT and iFFBT in MATLAB through discrete Fourier transform, matrix multiplication, and matrix trace function. ### Key points of the solution - **Asymptotic steerability**: Prove that when the band - limit is large enough, FFBT and iFFBT have asymptotic steerability. - **Numerical approximation**: Through appropriate truncation and sampling, transform the infinite series into a finite sum, thereby achieving numerical approximation. - **Matrix representation**: Use the matrix form to simplify the calculation and accelerate the calculation process through efficient algorithms such as fast Fourier transform (FFT). ### Application background The application fields of this research include but are not limited to: - **Image processing**: Process image signals with rotational symmetry. - **Cryo - electron microscopy (cryo - EM)**: Used for the structural analysis of biological macromolecules. - **Template matching technology**: Used for pattern recognition in digital image analysis. - **Computer vision**: Improve the performance of convolutional neural networks (CNNs). Through these methods, the authors hope to provide an efficient and accurate numerical tool for processing two - dimensional functions with rotational symmetry and convolution requirements.

Asymptotically Steerable Finite Fourier-Bessel Transforms and Closure under Convolution

Generalized Fourier-Bessel operator and almost-periodic interpolation and approximation

Method Based on Fast Fourier Transform for Calculating Conical Refraction of Beams with Noncircular Symmetry

Efficient polar convolution based on the discrete Fourier-Bessel transform for application in computational biophotonics

Convergence of fractional Fourier series on the torus and applications

Underdetermined Fourier Extensions for Partial Differential Equations on Surfaces

Circular Convolution and Product Theorem for Affine Discrete Fractional Fourier Transform

Computation of highly oscillatory integrals using a Fourier extension approximation

Asymptotic Expansions of Finite Hankel Transforms and the Surjectivity of Convolution Operators

A Fast Fourier-Galerkin Method for Solving Boundary Integral Equations on Torus-Shaped Surfaces

High-order FDTD schemes for Maxwell's interface problems with discontinuous coefficients and complex interfaces based on the Correction Function Method

Convergence analysis of variable steps BDF2 method for the space fractional Cahn-Hilliard model

A new construction of the Clifford-Fourier kernel

Development and comparison of spectral algorithms for numerical modeling of the quasi-static mechanical behavior of inhomogeneous materials

Accurate, high-order representation of complex three-dimensional surfaces via Fourier continuation analysis

Conforming Finite Element Function Spaces in Four Dimensions, Part 1: Foundational Principles and the Tesseract

Slice Fourier transform and convolutions

On the Fourier coefficients of powers of a Blaschke factor and strongly annular fonctions

Fourier Galerkin approximation of mean field control problems