滤波方法及装置、计算机设备及计算机可读存储介质制造方法及图纸

本发明专利技术公开了一种滤波方法及装置、计算机设备及计算机可读存储介质，该方法包括：对原始物探散点数据进行滤波处理，确定第二物探散点数据；根据物探散点数据中对应位置数据点的相对误差，利用相对误差不小于预设相对误差阈值的第二物探散点数据中的数据点替换原始物探散点数据中对应位置的数据点，得到第三物探散点数据；对第三物探散点数据进行滤波处理，确定第五物探散点数据。本发明专利技术根据原始物探散点数据与第二物探散点数据中每个对应位置的数据点的相对误差，利用相对误差不小于预设相对误差阈值的第二物探散点数据中的数据点替换原始物探散点数据中对应位置的数据点，消除原始物探散点数据中的畸变点，从而提高滤波效果。果。果。

【技术实现步骤摘要】
滤波方法及装置、计算机设备及计算机可读存储介质


[0001]本专利技术涉及地球物理勘探
，尤其涉及滤波方法及装置、计算机设备及计算机可读存储介质。

技术介绍

[0002]本部分旨在为权利要求书中陈述的本专利技术实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
[0003]三维物探散点数据用于地质解释或者三维成图时，一般先采用滤波方法对物探散点数据进行滤波处理，将三维物探散点数据变得平滑。直接对三维物探散点数据进行滤波处理的话，数据容易产生偏差，滤波效果难以控制。
[0004]在三维物探散点数据存在畸变点的情况下，虽然滤波后结果变得平滑，但不可避免的会受到畸变点的影响；特别是当畸变点的数值较大的情况下，这种影响会进一步加剧，最终影响地质解释或三维成图的效果。另外，中值滤波方法较少受到畸变点的影响，但中值滤波方法会使得极值和极值点的位置发生变化，从而影响滤波效果。目前，现有的平滑滤波方法或多或少都会受到畸变点的影响，即使滤波后畸变点数据明显变好，但滤波后畸变点周围的数据也会受到畸变点的污染和影响，导致滤本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种滤波方法，其特征在于，包括：对原始物探散点数据进行滤波处理，确定第二物探散点数据；第二物探散点数据为原始物探散点数据在原始位置的滤波结果；根据原始物探散点数据与第二物探散点数据中每个对应位置的数据点的相对误差，利用相对误差不小于预设相对误差阈值的第二物探散点数据中的数据点替换原始物探散点数据中对应位置的数据点，得到第三物探散点数据；数据点的相对误差反映数据点的畸变程度；对第三物探散点数据进行滤波处理，确定第五物探散点数据；第五物探散点数据为第三物探散点数据在原始位置的滤波结果。2.如权利要求1所述的方法，其特征在于，对原始物探散点数据进行滤波处理，确定第二物探散点数据，包括：利用规则插值网格对原始物探散点数据进行克里金插值，得到第一插值物探散点数据；采用第一余弦滤波窗口对第一插值物探散点数据进行第一余弦滤波，得到第一物探散点数据；第一余弦滤波窗口的滤波系数为第一预设滤波系数区间中任一滤波系数，第一余弦滤波窗口的中心位置数据不参与滤波；对第一物探散点数据进行反插值，得到原始物探散点数据在原始位置的滤波结果第二物探散点数据。3.如权利要求2所述的方法，其特征在于，规则插值网格的网格间隔不大于原始物探散点数据的数据点间隔，规则插值网格的节点密度为原始物探散点数据的节点密度的预设倍数。4.如权利要求1所述的方法，其特征在于，根据原始物探散点数据与第二物探散点数据中每个对应位置的数据点的相对误差，利用相对误差不小于预设相对误差阈值的第二物探散点数据中的数据点替换原始物探散点数据中对应位置的数据点，得到第三物探散点数据，包括：确定原始物探散点数据与第二物探散点数据中每个对应位置的数据点的相对误差；在每个对应位置的数据点的相对误差小于预设相对误差阈值时，将原始物探散点数据中对应位置的数据点作为第三物探散点数据对应位置的数据点；在每个对应位置的数据点的相对误差不小于预设相对误差阈值时，将第二物探散点数据中对应位置的数据点作为第三物探散点数据对应位置的数据点。5.如权利要求2所述的方法，其特征在于，对第三物探散点数据进行滤波处理，确定第五物探散点数据，包括：利用规则插值网格对第三物探散点数据进行克里金插值，得到第二插值物探散点数据；对第二插值物探散点数据进行高斯滤波，得到第四物探散点数据；高斯滤波的滤波系数为第二预设滤波系数区间中任一滤波系数；对第四物探散点数据进行反插值处理，得到第三物探散点数据在原始位置的滤波结果第五物探散点数据。6.如权利要求5所述的方法，其特征在于，第一余弦滤波窗口在各个方向上的节点数相
等且均为奇数，第一余弦滤波窗口在各个方向上的节点数不超过预设节点数阈值；和/或，高斯滤波在各个方向上的节点数相等且均为奇数，高斯滤波在各个方向上的节点数不超过预设节点数阈值。7.如权利要求5或6所述的方法，其特征在于，通过如下公式确定高斯滤波中每个数据的权系数：W＝exp(-factor
×
(fx
i
×
fx
i
+fy
j
×
fy
j
+fz
k
×
fz
k
))；fx
i
＝i/Wx，i＝-Wx/2～Wx/2；fy
j
＝j/Wy，j＝-Wy/2～Wy/2；fz
k
＝k/Wz，k＝-Wz/2～Wz/2；其中，W表示高斯滤波中数据(i,j,k)的权系数，i表示数据(i,j,k)的在x方向上的数据节点编号，j表示数据(i,j,k)在y方向上的数据节点编号，k表示数据(i,j,k)在z方向上的数据节点编号，Wx表示高斯滤波在x方向上的节点数，Wy表示高斯滤波在y方向上的节点数，Wz表示高斯滤波在z方向上的节点数，fx
i
表示在x方向上的数据节点i的权系数，fy
j
表示在y方向上的数据节点j的权系数，fz
k
表示在z方向上的数据节点k的权系数，factor表示高斯滤波的滤波系数。8.一种滤波装置，其特征在于，包括：第一滤波模块，用于对原始物探散点数据进行滤波处理，确定第二物探散点数据；第二物探散点数据为原始物探散点数据在原始位...

【专利技术属性】
技术研发人员：陶德强，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：