# Calculate Hankel Transform by Fast Fourier Transform

I am currently working on computing convolution by Fast Fourier Transform (FFT). The functions are in 2D however with circular symmetry $f(\mathbf{r}) = f(r)$. I have already implemented 2D FFT to solve this problem.

If I change the real space and Fourier space variables into cylindrical coordinates, I will end up with a Hankel transform in 1D.

$$F(\mathbf{k})=\iint f(\mathbf{r})e^{i\mathbf{k}\cdot\mathbf{r}}\operatorname{d}\!\mathbf{r}$$

$$F(\mathbf{k})=\int_{r=0}^\infty \int_{\theta=0}^{2\pi}f(r)e^{ikr\cos(\theta)}\,r\operatorname{d}\!\theta\operatorname{d}\!r$$

$$F(\mathbf{k})=F(k)= 2\pi\int_0^\infty f(r) J_0(kr) r\operatorname{d}\!r$$

I am wondering if I could carry out this 1D Hankel Transform by FFT?