这个函数使用大地水准面波动网格,并采用带状2D-FFT方法将其转换为重力异常值,以进行逆斯托克斯积分评估。这种方法适用于移除和恢复计算过程中的长波长全局模型。为了避免循环卷积,建议用100%零填充输入N。输入参数包括:N是一个nxm的大地水准面波动网格;Longm和Latm分别是N中每个网格节点的经度和纬度值网格(以度为单位);res是N网格的分辨率(以度为单位)。该代码实现了Hirt, C.、Featherstone, WE和Claessens, SJ(2011年)提出的公式。