平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法制造技术

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

One step upward continuation of plane modulus gradient field of plane geomagnetic anomaly field by recursive cosine transform

The present invention provides a recursive cosine transform method for one step upward continuation of plane modulus gradient field of plane geomagnetic anomaly field. The first plane of geomagnetic field measurement data preprocessing; using radial basis function interpolation on the observation plane irregular distribution of geomagnetic anomaly field discrete data grid, obtain geomagnetic grid abnormal measurement data. Then the use of recursive cosine transform geomagnetic grid abnormal measurement data is converted to the cosine transform spectrum; using step wavenumber domain conversion factor, calculate the modulus of each component of the geomagnetic field gradient spectrum by cosine transform spectrum of geomagnetic anomaly. Finally, the horizontal inverse x component, the horizontal Y component and the vertical Z component are carried out, and the phase shifted phase shifted back cosine transform of the pion /2 is obtained, and the component of the geomagnetic modulus gradient field in the continuation plane is obtained. The calculation amount of the one-step method is about half of the calculation amount of the two step method; the cosine transform method reduces the Gibbs boundary effect and improves the algorithm accuracy; it can be applied to arbitrary length data sequences.

【技术实现步骤摘要】
平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法
本专利技术涉及的是一种地磁位场数据处理方法，具体涉及到一种平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法。
技术介绍
地磁延拓是一种重要的地磁位场数据处理方法，在地磁匹配定位、地磁观测及其资料解释等方面具有重要应用。地磁基准图是实现地磁匹配定位的数据基础。为进行地磁匹配定位，一般需要获得航行深度面上的地磁基准图。不同航行器的航行高度各不相同，通过直接测量方式得到各航行高度上的地磁数据是不现实的。大面积的地磁测量值多数为地磁异常场数据，地磁模量梯度测量虽已开始，但匹配定位需要大范围导航区域的基准图，现有的地磁模量梯度测量数据难以满足这一要求。根据已有的航测地磁数据或者海测地磁数据，通过延拓和导数的同时计算一步获得地磁模量梯度场，是获得航行面地磁模量梯度场的有效途径之一。向上延拓还可融合不同高度的航磁数据，使之归化到同一高度面。地磁观测及其资料解释方面，利用向上延拓可以进行区域场与局部场、深源场与浅源场的分离。地磁模量梯度场具有以下优点：第一，地磁模量梯度与不同磁传感器间的姿态无关，降低了多个传感器测量轴配置精度的要求；本文档来自技高网...

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

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

【专利技术属性】
技术研发人员：黄玉，武立华，吕振川，齐瑞云，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江,23