一种利用小子域识别算法确定磁性目标位置的方法技术

技术编号:13428364 阅读:57 留言:0更新日期:2016-07-29 18:30
本发明专利技术公开了一种利用小子域识别算法确定磁性目标位置的方法。设置初始阈值K=0.1,计算多个子域中对应平均值Δgi的平均值ΔSg和对应的均方差Sσ,比较Sσ和阈值K的大小,当Sσ<K时,窗口的中心位置的识别输出值设置为0;当Sσ>K时,窗口的中心位置的识别输出值设置为1。将窗口滑动到下一点,重复上述计算,直至完成整个区域的第一次计算。将整个区域中最大的Sσ赋予阈值K。用上述方法对整改区域重新进行循环计算,直至满足循环截止条件。本发明专利技术根据磁异常信息,通过设计特定的算法,判断磁性目标的空间位置、大小、形状等特征参数,可实现目标监控和识别。

【技术实现步骤摘要】

本专利技术属于磁性目标探测领域,尤其涉及一种在地磁测量及地磁反演应用中利用小子域识别算法确定磁性目标位置的方法
技术介绍
地磁场一般可分为地磁正常场和地磁异常场。在实际测量过程中,正常场与异常场的划分往往是相对的。实际上测得的地球磁场强度和理论磁场强度是有区别的,这种区别称为地磁异常。它主要是由磁性物体受地磁场磁化而产生的附加磁场。所以,通过对磁异常的研究可以帮助人们有效的寻找磁性物体。根据磁异常信息,通过一定的数学解释手段,判断磁性目标的空间位置、大小、形状等特征参数,可实现目标监控和识别。Euler反褶积方法是以欧拉齐次方程为基础,推断磁性目标中心的空间位置。然而反演场源水平位置和深度时,Euler齐次方程是一个灵敏度极高的方程式,式中异常导数的计算精度等稍有偏差,便会在反演中产生较强的干扰。小子域识别算法是基于滑动平均法原理进行设计的算法,不仅保留了压缩异常区域的功能,能以较高的分辨率保留异常中的区域特征,对具有梯级带特征的断裂构造位置的确定具有较高的分辨率,同时可以有效的识别窗口是否处于异常区域内,这样就可以在压缩异常区域过程中,将异常区域的中心位置显示出来。举例说明算法的实现过程:当磁场梯度带从左向右增大时,小子域的滤波窗口从左向右滑动,各小子域中以平缓系数最小域内的平均值作为窗口中心位置的值。当窗口右侧的子域处于磁异常的梯度带上,这些子域的均方差较大;而左测的子域处于磁异常的平缓区时,这些子域的均方差较小。因此窗口中心位置的值由较大值变为较小值,这样相当于窗口左边低值的平缓区界限右移。同时,计算多个子域中对应平均值Δgi的平均值ΔSg和对应的均方差Sσ,当Sσ较大时,说明子域中磁场变化剧烈,该窗口处于异常区域中;当Sσ较小时,说明子域中磁场变化平缓,该窗口处于异常区域外。同理,当窗口左侧的子域处于磁异常的梯度带上,而右侧的子域处于磁异常的平缓区时,通过计算,实现窗口右边低值的平缓区界限左移。通过比较Sσ和阈值K的大小,最终达到了压缩梯级带来确定异常区域的中心位置。
技术实现思路
本专利技术的目的是提供一种准确的,具有高可靠识别性的,利用小子域识别算法确定磁性目标位置的方法。一种利用小子域识别算法确定磁性目标位置的方法,包括以下步骤,步骤一:利用磁力仪,获取区域的局部地磁图数据;步骤二:构建5X5数据点的窗口,将窗口划分成9个小子域,设定初始阈值K=0.1;步骤三:分别计算9个子域里的磁异常方差σi,将9个子域中的方差最小对应的小子域的磁异常均值Δgmin作为窗口中心的值;Δgmin=min{Δg1,Δg2,…,Δg9本文档来自技高网
...
一种利用小子域识别算法确定磁性目标位置的方法

【技术保护点】
一种利用小子域识别算法确定磁性目标位置的方法,其特征在于:包括以下步骤,步骤一:利用磁力仪,获取区域的局部地磁图数据;步骤二:构建5X5数据点的窗口,将窗口划分成9个小子域,设定初始阈值K=0.1;步骤三:分别计算9个子域里的磁异常方差σi,将9个子域中的方差最小对应的小子域的磁异常均值Δgmin作为窗口中心的值;其中,Ti(p,q)为第i个子域内第(p,q)个点上的异常值,σi为第i个子域的方差值;步骤四:计算整个窗口的方差Sσ:Sσ=1p×qΣp=1,q=15(T(p,q)-Δgmin)2;]]>步骤五:比较整个窗口的方差Sσ和阈值K的大小,当Sσ<K时,窗口中心位置的识别输出值设置为0;当Sσ>K时,窗口中心位置的识别输出值设置为1;步骤六:窗口滑动到下一点,重复步骤二~步骤五,直至完成全部的区域计算;将整个区域中最大的整个窗口的方差Sσ赋予阈值K;步骤七:重复步骤二~步骤六,判断整个窗口的方差Sσ是否和阈值K相等,如果相等,则进行步骤八;如果不相等,重复步骤二~步骤五,直至整个窗口的方差Sσ和阈值K相等;步骤八:得到地磁图中识别标识为1的位置为磁性目标所处位置。...

【技术特征摘要】
1.一种利用小子域识别算法确定磁性目标位置的方法,其特征在于:包括以下步骤,
步骤一:利用磁力仪,获取区域的局部地磁图数据;
步骤二:构建5X5数据点的窗口,将窗口划分成9个小子域,设定初始阈值K=0.1;
步骤三:分别计算9个子域里的磁异常方差σi,将9个子域中的方差最小对应的小子域的
磁异常均值Δgmin作为窗口中心的值;
其中,Ti(p,q)为第i个子域内第(p,q)个点上的异常值,σi为
第i个子域的方差值;
步骤四:计算整个窗口的方差Sσ:
Sσ=1p×qΣp=1,q=15(T(...

【专利技术属性】
技术研发人员:康崇樊黎明张晓峻王明郑权周健
申请(专利权)人:哈尔滨工程大学
类型:发明
国别省市:黑龙江;23

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1