球面上不规则网格的球面谐波变换

时间:2017-07-26 06:23:38

标签: python fft data-processing

假设我在球体表面上有许多数据点。这些数据点随机散布在球体表面上,即它们不会分布在规则网格上,例如高斯网格。

我的问题是:是否有可用的Python库可以计算这个数据集的球谐函数变换?否则,在应用标准球面谐波变换之前,我首先必须自己将数据点插入到规则网格中。

对于非均匀傅立叶变换,我发现了几个选项,例如: this one.但是,对于非均匀球面谐波变换(当然与傅里叶变换有关),到目前为止我没有找到。

2 个答案:

答案 0 :(得分:1)

我认为这是一个老问题,但无论如何我都会发布这个问题。

python或Fortran中可用于球谐函数的最佳库之一是shtools。具体来说,检查函数SHExpandLSQ,它使用最小二乘法在不规则采样网格上找到球谐调节展开系数。

答案 1 :(得分:0)

这是一种方法,您首先将数据点插入到常规网格中吗?如果是,numpy / scipy有这样做的方法。看一下手册:

https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.interpolate.griddata.html

https://docs.scipy.org/doc/scipy-0.16.1/reference/generated/scipy.interpolate.RegularGridInterpolator.html

和一个简单的例子解释了here(方法2)