[发明专利]平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法

摘要

本发明提供的是一种平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法。先进行平面地磁异常场测量数据预处理；利用径向基函数插值法将观测平面上不规则离散分布的地磁异常场测量数据进行网格化，获得网格化的地磁异常测量数据。再利用递归余弦变换将网格化的地磁异常测量数据转换为余弦变换波谱；利用一步波数域转换因子，由地磁异常余弦变换波谱计算地磁模量梯度场各分量的波谱。最后分别对水平x分量、水平y分量、垂直z分量的波谱进行行向移相π/2的递归反余弦变换得到延拓平面上的地磁模量梯度场的分量。本发明的一步法的计算量约为两步法计算量的一半；余弦变换法减小Gibbs边界效应，提高了算法精度；可应用于任意长度的数据序列。

主权项

一种平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法，其特征是：步骤1、选取网格化参数，对平面地磁异常场测量数据进行预处理；步骤2、利用径向基函数插值法，将观测平面上的离散地磁异常场测量数据进行网格化，获得观测平面上的网格化地磁异常测量数据；步骤3、先后对观测平面上的地磁异常场进行行向和列向的一维递归余弦变换，得到地磁异常场的二维余弦变换波谱；步骤4、由观测平面上的地磁异常场波谱计算延拓平面上的地磁模量梯度场的水平x、水平y和垂直z波谱；步骤5、先后对延拓平面上的地磁模量梯度场的水平x波谱进行行向的移相π/2的一维递归反余弦变换和列向的一维递归反余弦变换，定义为行向正弦、列向余弦反变换，得到延拓平面上的地磁模量梯度场的水平x分量；步骤6、先后对延拓平面上的地磁模量梯度场的水平y和波谱进行行向的一维递归反余弦变换和列向的移相π/2的一维递归反正弦变换，定义为行向余弦、列向正弦反变换，得到延拓平面上的地磁模量梯度场的水平y分量；步骤7、先后对延拓平面上的地磁模量梯度场的垂直z波谱进行行向和列向的一维递归反余弦变换，得到延拓平面上的地磁模量梯度场的垂直z分量。