此函数采用大地水准面波动网格,并使用带状 2D-FFT 方法将它们转换为重力异常值,以评估逆斯托克斯积分。
该函数最适合用于移除-计算-恢复类型的过程(即 N 应该是剩余大地水准面波动的网格是长波长全局模型的减法)。 用 100% 零填充输入 N 以避免循环卷积也是有益的。
Gravity_Anom=Inv_Stokes(Longm,Latm,N,res)
输入如下:
-N是一个nxm大地水准面波动网格- Longm 是 N 中每个网格节点的经度值网格(以度为单位) - Latm 是 N 中每个网格节点的纬度值网格(以度为单位) - res 是 N 网格的分辨率(即一弧分网格将具有 res=1/60)(以度为单位)
该代码是提供的公式的实现; Hirt, C.、Featherstone, WE 和 Claessens, SJ,2011 年。关于大地测量卷积积分的准确数值评估。
2021-12-22 14:03:45
2KB
matlab
1