非局部均值滤波在重力异常数据处理中的应用
|
摘要
根据小子域滤波原理,小子域滤波产生异常边界扭曲主要表现为高频干扰,提出小子域滤波与非局部均值滤波联合的方法以克服小子域滤波存在的异常边界扭曲问题。利用小子域滤波良好的边界增强能力和非局部均值滤波高频细节保护能力,在压制小子域滤波缺陷的同时尽可能保护小子域滤波的边界增强效果。通过模拟计算和莱州湾海域实测重力数据处理表明,本文方法在增强异常边界的同时压制小子域滤波产生的边界扭曲效应,获得了更好的重力异常处理效果。与单一小子域滤波相比,本文方法不需要考虑异常边界形态,具有更强的稳定性和普适性。
According to the small subdomain filtering principle,the small subdomain filtering produces abnormal boundary distortion mainly as high frequency interference. A combination of a small subdomain filtering and non-local mean filtering method is proposed to overcome the abnormal boundary distortion problem in small subdomain filtering. The good boundary enhancement ability of small sub-domain filtering and the high-frequency detail protection ability of non-local mean filtering are used to suppress the small sub-domain filtering defect while protecting the boundary enhancement effect of the small sub-domain filtering as much as possible. The simulation calculation and the measured gravity data processing in the Laizhou Bay area show that the proposed method suppresses the boundary distortion effect generated by the small sub-domain filtering while enhancing the anomaly boundary,and obtains better gravity anomaly processing result. Compared with single small domain filtering,the proposed method does not need to consider the abnormal boundary shape,and is more stable and universal.
According to the small subdomain filtering principle,the small subdomain filtering produces abnormal boundary distortion mainly as high frequency interference. A combination of a small subdomain filtering and non-local mean filtering method is proposed to overcome the abnormal boundary distortion problem in small subdomain filtering. The good boundary enhancement ability of small sub-domain filtering and the high-frequency detail protection ability of non-local mean filtering are used to suppress the small sub-domain filtering defect while protecting the boundary enhancement effect of the small sub-domain filtering as much as possible. The simulation calculation and the measured gravity data processing in the Laizhou Bay area show that the proposed method suppresses the boundary distortion effect generated by the small sub-domain filtering while enhancing the anomaly boundary,and obtains better gravity anomaly processing result. Compared with single small domain filtering,the proposed method does not need to consider the abnormal boundary shape,and is more stable and universal.
引文
[1]穆石敏,申宁华,孙运生.区域地球物理数据处理方法及其应用[M].长春:科学技术出版社,1990:1-292.MU Shi-min,SHEN Ning-hua,SUN Yun-sheng. Regional geophysical data processing methods and their applications[M]. Changchun:Science and Technology Press,1990:1-292.
[2]李红雨,杨长保,吴燕冈,等.小波变换在位场资料去噪和位场分离中的应用[J].世界地质,2014,33(1):200-208.LI Hong-yu,YANG Chang-bao,WU Yan-gang,et al.Application of wavelet transform for de-noising and potential field separation in gravity and magnetic data processing[J]. Global Geology,2014,33(1):200-208.
[3]杨高印.位场数据处理的一项新技术─小子域滤波法[J].石油地球物理勘探,1995,30(2):240-244.YANG Gao-yin. A new technique for potential-field data processing:small subdomain filtering[J]. Oil Geophysical Prospecting,1995,30(2):240-244.
[4]马涛,王铁成,王雨.一种改进的网格数据保持梯度滤波方法[J].石油地球物理勘探,2007,42(2):198-201.MA Tao,WANG Tie-cheng,WANG Yu. An improved gradient-preserved grid data filtering approach[J]. Oil Geophysical Prospecting,2007,42(2):198-201.
[5]张凤旭,张凤琴,刘财,等.断裂构造精细解释技术:三方向小子域滤波[J].地球物理学报,2007,50(5):1543-1550.ZHANG Feng-xu,ZHANG Feng-qin,LIU Cai,et al. A technique for elaborate explanation of faulted structures:three-directional small subdomain filtering[J]. Chinese Journal of Geophysics,2007,50(5):1543-1550.
[6]肖锋,吴燕冈,孟令顺.位场数据处理中小子域滤波的改进[J].石油地球物理勘探,2010,45(1):136-139.XIAO Feng,WU Yan-gang,MENG Ling-shun. An improvement on small subdomain filtering in potential field data processing[J]. Oil Geophysical Prospecting,2010,45(1):136-139.
[7] Jiang F Y,Gao L K. Edge enhancement of gravity gradient tensor data using small sub-domain filtering method[J]. Applied Geophysics,2012,9(2):119-130.
[8]马国庆,杜晓娟,李丽丽.优化小子域滤波方法研究及其应用[J].石油地球物理勘探,2013,48(4):658-662.MA Guo-qing,DU Xiao-juan,LI Li-li. An optimized small sub-domain filtering method[J]. Oil Geophysical Prospecting,2013,48(4):658-662.
[9]马国庆,孟令顺,杜晓娟.多环带中值滤波法实现重力区域场与局部场的划分[J].世界地质,2010,29(1):113-117.MA Guo-qing,MENG Ling-shun,DU Xiao-juan. Division of regional gravity field and local gravity field by multi-ring median filtering method[J]. Global Geology,2010,29(1):113-117.
[10]祝严刚,张桂梅.一种改进的非局部均值图像去噪算法[J].计算机工程与应用,2017,53(18):192-198.ZHU Yan-gang,ZHANG Gui-mei. Improved non-local means denoising algorithm[J]. Computer Engineering and Applications,2017,53(18):192-198.
[11] Buades A,Coll B,Morel J M. A non-local algorithm for image denoising[C]//2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition(CVPR'05),2005,2:60-65.
[12] Hemalata B,Subhasis C. Novel speed-up strategies for nonlocal means denoising with patch and edge patch based dictionaries[J]. IEEE Transactions on Image Processing,2014,23(1):356-365.
[13]张权,罗立民,桂志国,等.一种基于优化参数的非局部均值滤波算法[J].计算机应用与软件,2012,29(3):78-81.ZHANG Quan,LUO Li-min,GUI Zhi-guo,et al. A non-local mean filtering algorithm based on optimum parameter[J]. Computer Applications and Software,2012,29(3):78-81.
[14]郭贝贝,易三莉,贺建峰,等.改进的非局部均值滤波算法[J].计算机工程,2016,42(7):227-231.GUO Bei-bei,YI San-li,HE Jian-feng,et al. Improved non-local means filtering algorithm[J]. Computer Engineering,2016,42(7):227-231.
[15]苏超,杜晓娟,马国庆,等.几何平均重力归一化总梯度在山东招远金矿采空区的应用[J].世界地质,2014,33(4):889-894.SU Chao,DU Xiao-juan,MA Guo-qing,et al. Appication of normalized full gradient method with geometric mean in mined-out area of Zhaoyuan gold deposit in Shandong[J]. Global Geology,2014,33(4):889-894.
[2]李红雨,杨长保,吴燕冈,等.小波变换在位场资料去噪和位场分离中的应用[J].世界地质,2014,33(1):200-208.LI Hong-yu,YANG Chang-bao,WU Yan-gang,et al.Application of wavelet transform for de-noising and potential field separation in gravity and magnetic data processing[J]. Global Geology,2014,33(1):200-208.
[3]杨高印.位场数据处理的一项新技术─小子域滤波法[J].石油地球物理勘探,1995,30(2):240-244.YANG Gao-yin. A new technique for potential-field data processing:small subdomain filtering[J]. Oil Geophysical Prospecting,1995,30(2):240-244.
[4]马涛,王铁成,王雨.一种改进的网格数据保持梯度滤波方法[J].石油地球物理勘探,2007,42(2):198-201.MA Tao,WANG Tie-cheng,WANG Yu. An improved gradient-preserved grid data filtering approach[J]. Oil Geophysical Prospecting,2007,42(2):198-201.
[5]张凤旭,张凤琴,刘财,等.断裂构造精细解释技术:三方向小子域滤波[J].地球物理学报,2007,50(5):1543-1550.ZHANG Feng-xu,ZHANG Feng-qin,LIU Cai,et al. A technique for elaborate explanation of faulted structures:three-directional small subdomain filtering[J]. Chinese Journal of Geophysics,2007,50(5):1543-1550.
[6]肖锋,吴燕冈,孟令顺.位场数据处理中小子域滤波的改进[J].石油地球物理勘探,2010,45(1):136-139.XIAO Feng,WU Yan-gang,MENG Ling-shun. An improvement on small subdomain filtering in potential field data processing[J]. Oil Geophysical Prospecting,2010,45(1):136-139.
[7] Jiang F Y,Gao L K. Edge enhancement of gravity gradient tensor data using small sub-domain filtering method[J]. Applied Geophysics,2012,9(2):119-130.
[8]马国庆,杜晓娟,李丽丽.优化小子域滤波方法研究及其应用[J].石油地球物理勘探,2013,48(4):658-662.MA Guo-qing,DU Xiao-juan,LI Li-li. An optimized small sub-domain filtering method[J]. Oil Geophysical Prospecting,2013,48(4):658-662.
[9]马国庆,孟令顺,杜晓娟.多环带中值滤波法实现重力区域场与局部场的划分[J].世界地质,2010,29(1):113-117.MA Guo-qing,MENG Ling-shun,DU Xiao-juan. Division of regional gravity field and local gravity field by multi-ring median filtering method[J]. Global Geology,2010,29(1):113-117.
[10]祝严刚,张桂梅.一种改进的非局部均值图像去噪算法[J].计算机工程与应用,2017,53(18):192-198.ZHU Yan-gang,ZHANG Gui-mei. Improved non-local means denoising algorithm[J]. Computer Engineering and Applications,2017,53(18):192-198.
[11] Buades A,Coll B,Morel J M. A non-local algorithm for image denoising[C]//2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition(CVPR'05),2005,2:60-65.
[12] Hemalata B,Subhasis C. Novel speed-up strategies for nonlocal means denoising with patch and edge patch based dictionaries[J]. IEEE Transactions on Image Processing,2014,23(1):356-365.
[13]张权,罗立民,桂志国,等.一种基于优化参数的非局部均值滤波算法[J].计算机应用与软件,2012,29(3):78-81.ZHANG Quan,LUO Li-min,GUI Zhi-guo,et al. A non-local mean filtering algorithm based on optimum parameter[J]. Computer Applications and Software,2012,29(3):78-81.
[14]郭贝贝,易三莉,贺建峰,等.改进的非局部均值滤波算法[J].计算机工程,2016,42(7):227-231.GUO Bei-bei,YI San-li,HE Jian-feng,et al. Improved non-local means filtering algorithm[J]. Computer Engineering,2016,42(7):227-231.
[15]苏超,杜晓娟,马国庆,等.几何平均重力归一化总梯度在山东招远金矿采空区的应用[J].世界地质,2014,33(4):889-894.SU Chao,DU Xiao-juan,MA Guo-qing,et al. Appication of normalized full gradient method with geometric mean in mined-out area of Zhaoyuan gold deposit in Shandong[J]. Global Geology,2014,33(4):889-894.