Are Fast (or the so-called quasi-Fast) Hankel Transforms implemented or
implementable in Mathematica?
Two-dimensional Fourier Transforms of functions f[r] where r=Sqrt[x^2+y^2]
are actually one-dimensional Hankel Transforms of order zero. So I was
wondering if one could use somehow 1D FFTs to optimize 2D FTs of f[r].
What I do now is rather lame:
Fourier[ Table[ f[ Sqrt[i^2+j^2] ],{j,-n/2+1,n/2},{i,-n/2+1,n/2}] ]
and I would like something similar to
FHT[ Table[ f[i],{i,n} ] ]
to have directly the radially symmetric 2D spectrum.
Thanks,
Bye
Hyper