我想采用以下形式的 FFT： 球坐标中的傅立叶
所以实际上，myf(theta,phi)是一个 3D 图像体素，即它有 f(radius,theta,phi)。如果我必须对 dx 进行积分，那相当于只选择我的 IMG(ylocation,:,zlocation)=currentvoxel 并在其上运行 fft(currentvoxel)。但是，现在我只需要沿着 phi 进行 fft。这很难，因为它需要我找到所有具有 theta=currenttheta 的位置，并且对于每个 r=currentr，我必须找到所有的 phi 来获取 fft。
我什至不知道如何展开 phi 以获取 fft。当它是笛卡尔坐标时，这很容易，只需从 x=0,..x=end 开始并运行 fft。phi大小写一样吗？获取某个 r,theta 上的所有 phi 并将它们排序为 0 到 2pi并取 fft？
或者可能像这样的正弦/余弦变换： 余弦/正弦变换
如果这是 af(x,y)，我认为在 MATLAB 中执行此操作的方法是采用 dct(f*(x.^2+y.^2))。换句话说，运行离散余弦变换。我对吗？
第二部分的额外说明：我感兴趣的积分是：bessel where j_l(kr):扩展 j_l 并 最终进一步扩展 F_(l)(t) 由前面给出：余弦/正弦变换
