[发明专利]一种重磁数据处理方法及系统

权利要求书

1.一种重磁数据处理方法，其特征在于，所述方法包括：

获取重磁异常数据；

将所述重磁异常数据中有数据缺失的位置补零，获得补零后重磁位场数据；

计算所述补零后重磁位场数据的归一化径向平均功率谱；所述计算所述补零后重磁位场数据的归一化径向平均功率谱，具体包括：

采用快速傅里叶变换方法求解所述补零后重磁位场数据的重磁位场数据功率谱；

以所述重磁位场数据功率谱的中心为圆心，以基频的整数倍为半径作圆，生成多个环带；每个所述环带对应一个不同的频率；

求解各个所述环带内所有点的功率谱值的平均值，获得多个平均功率谱值；

以径向波数为横坐标，以各个所述环带内的所述平均功率谱值的对数为纵坐标，生成径向平均功率谱；

对所述径向平均功率谱作归一化处理，生成所述重磁位场数据的归一化径向平均功率谱；

根据所述归一化径向平均功率谱确定截止波数；所述根据所述归一化径向平均功率谱确定截止波数，具体包括：

根据所述归一化径向平均功率谱的曲线特征，选取所述归一化径向平均功率谱下降预设比例a所对应的波数为截止波数c_K；

根据所述截止波数对所述补零后重磁位场数据进行去噪、插值和扩边处理，获得处理后的重磁数据；所述根据所述截止波数对所述补零后重磁位场数据进行去噪、插值和扩边处理，获得处理后的重磁数据，具体包括：

对所述补零后重磁位场数据作二维傅里叶变换，获得所述补零后重磁位场数据的频谱G(u,v)；

根据所述截止波数确定低通滤波函数其中c_k表示第k次迭代的截止波数，1≤k≤K且满足c₁≤c₂≤…≤c_K和1＜c_k≤min(M,N)，M和N是插值扩边后重磁位场数据的尺寸；D(u,v)＝[(u-M/2)²+(v-N/2)²]^1/2；

采用所述低通滤波函数T_k(u,v)对所述频谱G(u,v)进行低通滤波处理，获得低通滤波后留下的谱分量；

对所述谱分量作二维傅里叶反变换，获得反变换后的重磁位场数据；

利用所述反变换后的重磁位场数据重置所述补零后重磁位场数据，获得重置后的重磁位场数据；

判断迭代次数k是否不大于预设迭代次数K，获得第一判断结果；

若所述第一判断结果为所述迭代次数k不大于所述预设迭代次数K，令k＝k+1，返回所述根据所述截止波数确定低通滤波函数的步骤；

若所述第一判断结果为所述迭代次数k大于所述预设迭代次数K，输出所述重置后的重磁位场数据作为处理后的重磁数据。

2.根据权利要求1所述的重磁数据处理方法，其特征在于，所述利用所述反变换后的重磁位场数据重置所述补零后重磁位场数据，获得重置后的重磁位场数据，具体包括：

将所述补零后重磁位场数据中没有数据的位置用所述反变换后的重磁位场数据替换，有数据的位置保留，获得所述重置后的重磁位场数据。