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

权利要求书

1.一种平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法，其特征是：

步骤1、选取网格化参数，对平面地磁异常场测量数据进行预处理；

步骤2、利用径向基函数插值法，将观测平面上的离散地磁异常场测量数据进行网格化，获得观测平面上的网格化地磁异常测量数据；

步骤3、先后对观测平面上的地磁异常场进行行向和列向的一维递归余弦变换，得到地磁异常场的二维余弦变换波谱；

步骤4、由观测平面上的地磁异常场波谱计算延拓平面上的地磁模量梯度场的水平x、水平y和垂直z波谱，具体包括：

按式由观测平面上的地磁异常场的二维余弦变换波谱ΔT^C(u_x,u_y,z₀)计算延拓平面上的地磁模量梯度场的波谱和式中，R_1WC(u_x,u_y,Δz)为余弦变换下的由观测平面地磁异常场计算延拓平面地磁模量梯度场的一步波数域转换因子，其表达式为

步骤5、先后对延拓平面上的地磁模量梯度场的水平x波谱进行行向的移相π/2的一维递归反余弦变换和列向的一维递归反余弦变换，定义为行向正弦、列向余弦反变换，得到延拓平面上的地磁模量梯度场的水平x分量；

步骤6、先后对延拓平面上的地磁模量梯度场的水平y和波谱进行行向的一维递归反余弦变换和列向的移相π/2的一维递归反正弦变换，定义为行向余弦、列向正弦反变换，得到延拓平面上的地磁模量梯度场的水平y分量；

步骤7、先后对延拓平面上的地磁模量梯度场的垂直z波谱进行行向和列向的一维递归反余弦变换，得到延拓平面上的地磁模量梯度场的垂直z分量。