CN105488141A - 一种探地雷达数据的多尺度可视方法 - Google Patents

一种探地雷达数据的多尺度可视方法 Download PDF

Info

Publication number
CN105488141A
CN105488141A CN201510835612.7A CN201510835612A CN105488141A CN 105488141 A CN105488141 A CN 105488141A CN 201510835612 A CN201510835612 A CN 201510835612A CN 105488141 A CN105488141 A CN 105488141A
Authority
CN
China
Prior art keywords
sub
scale
axis
section
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510835612.7A
Other languages
English (en)
Other versions
CN105488141B (zh
Inventor
原达
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Technology and Business University
Original Assignee
Shandong Technology and Business University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong Technology and Business University filed Critical Shandong Technology and Business University
Priority to CN201510835612.7A priority Critical patent/CN105488141B/zh
Publication of CN105488141A publication Critical patent/CN105488141A/zh
Application granted granted Critical
Publication of CN105488141B publication Critical patent/CN105488141B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90Details of database functions independent of the retrieved data types
    • G06F16/904Browsing; Visualisation therefor

Landscapes

  • Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了一种探地雷达数据的多尺度可视方法,建立光滑模型,对子剖面进行光滑处理,经过光滑处理后,突出了数据中蕴含的结构特征,采用基于FCM的方法获取子剖面的主结构;将分布在各个子剖面的主结构汇聚起来,形成多尺度的主题映射。本发明的有益效果是将探地雷达数据中不明确的结构特征以多尺度的形式呈现出来,有规律地增强数据的可视深度和广度,帮助数据分析人员更深入地了解结构特征存在的可能性,为综合分析提供技术支持。

Description

一种探地雷达数据的多尺度可视方法
技术领域
本发明属于雷达探测技术领域,涉及一种探地雷达数据的多尺度可视方法。
背景技术
探地雷达(GroundPenetratingRadar,GPR)是用无线电波来确定地下介质分布规律的探测技术,广泛应用于人类生存和生活服务,包括石油矿产资源勘查、基础工程建设、军事、环境保护和灾害防治等。探地数据分析与解译是将GPR数据转换为有用信息的核心环节。具体来说,就是依据获得的散射场,以及散射数据中的相位、频率和波速等参量,推断出区域内不同介质的空间结构分布及属性特征。
目前GPR数据分析与解译工作仍然停留在强烈依赖解译人员经验的阶段,其结果的产生常常因个人观点和经验的不同而有着很大的差异。原因在于数据处理的不确定性对数据的解译带来了极大的困难。这种困难并不仅仅是由于处理方法的缺陷或技巧造成的,而是观测技术本身也存在着难以逾越的障碍。
GPR系统的动态范围至少为60dB,而计算机屏幕上图像显示的动态范围大约是10~20dB,这意味着仅有少部分可用信息能够以图像的形式直接显示出来,其他附加的信息没有得到很好的显示和利用。本技术对原始数据中内蕴的属性信息进行充分挖掘和关联,使以往难以察觉地的规律显露出来。
目前的探地雷达数据多尺度可视化方法复杂,效率低,精度低。
发明内容
本发明的目的在于提供一种探地雷达数据的可视方法,解决了目前的探地雷达数据多尺度可视化方法复杂的问题。
本发明所采用的技术方案是按照以下步骤进行:
步骤1:建立光滑模型,对子剖面进行光滑处理,经过光滑处理后,突出了数据中蕴含的结构特征,光滑模型表示为:
inf v { T ( v ) = ∫ | ▿ v | + λ F ‾ 2 } ;
v为光滑近似,T(v)为基于参数v的子剖面,λ为调节参数,为保真项,用来控制模型的拟合度,第一项为正则项;
步骤2:主结构提取:
采用基于FCM的方法获取子剖面的主结构;
步骤3:主题映射;将分布在各个子剖面的主结构汇聚起来,形成多尺度的主题映射,主题映射定义为:
min x D ( x , f ) , s . t . L ( x , g ) = Σ i = 1 k A i z i ;
x为主题映射,D(·)表示x与f之间的欧式距离,z为子剖面的主题结构,A为权重矩阵,g为可变尺度,L(·)表示各个子剖面在尺度g下到x的集成变换。
进一步,所述步骤1中的计算:
| ▿ v ( x , y ) | = ( | ▿ v x ( x , y ) | 2 + | ▿ v y ( x , y ) | 2 ) 1 / 2 ;
其中,
| ▿ x ( x , y ) | = | ∂ x R ( x , y ) | ;
| ▿ y ( x , y ) | = | ∂ y C ( x , y ) | ;
这里,R(x,y)和C(x,y)分别为沿X轴和Y轴的变换函数,即
R ( x , y ) = 1 M Σ i , j ∈ l ( x , y ) δ l o g ( v ( i , j ) p ( i , j ) ) ;
C ( x , y ) = 1 N Σ i , j ∈ h ( x , y ) G ( k , σ ) * v ( i , j ) ;
l(x,y)和h(x,y)分别为沿X轴和Y轴的邻域,G(k,σ)为高斯核函数,k和σ为光滑度控制参数,M和N分别为沿X轴和Y轴的邻域范围,i和j为邻域l(x,y)和h(x,y)的元素标识,δ为单位函数。
进一步,所述步骤1中为保真项:
F ‾ ( x , y ) = t ( x , y ) · q ( x , y ) 1 - q ( x , y ) ;
其中, q ( x , y ) = | s ( x , y ) - v ( x , y ) | | s ( x , y ) | ,
t ( x , y ) = | s ( x , y ) - m y ( x , y ) ϵ + s ( x , y ) - m x ( x , y ) | ,
t(x,y)为方向相关的对比度,mx(x,y)和my(x,y)分别为与X和Y方向相关的均值函数。
本发明的有益效果是将探地雷达数据中不明确的结构特征以多尺度的形式呈现出来,有规律地增强数据的可视深度和广度,帮助数据分析人员更深入地了解结构特征存在的可能性,为综合分析提供技术支持。
附图说明
图1是可变尺度g取值为256至1的可视化结果示意图。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
对原始GPR数据进行分解,形成多个独立的窄带子剖面。
步骤1:建立光滑模型,通过纹理光滑获得显著性区域,这些显著性区域经过融合变换,可以得到主模式。
光滑模型表示为:
inf v { T ( v ) = ∫ | ▿ v | + λ F ‾ 2 } ;
这里,v为光滑近似,T(v)为基于参数v的子剖面。λ为调节参数。为保真项,用来控制模型的拟合度。第一项为正则项。由于GPR数据具有异步属性,我们给出以下定义计算
| ▿ v ( x , y ) | = ( | ▿ v x ( x , y ) | 2 + | ▿ v y ( x , y ) | 2 ) 1 / 2 ;
其中,
| ▿ x ( x , y ) | = | ∂ x R ( x , y ) | ;
| ▿ y ( x , y ) | = | ∂ y C ( x , y ) | ;
R(x,y)和C(x,y)分别为沿X轴和Y轴的变换函数,即:
R ( x , y ) = 1 M Σ i , j ∈ l ( x , y ) δ l o g ( v ( i , j ) p ( i , j ) ) ;
C ( x , y ) = 1 N Σ i , j ∈ h ( x , y ) G ( k , σ ) * v ( i , j ) ;
l(x,y)和h(x,y)分别为沿X轴和Y轴的邻域,G(k,σ)为高斯核函数,k和σ为光滑度控制参数。M和N分别为沿X轴和Y轴的邻域范围。i和j为邻域l(x,y)和h(x,y)的元素标识。δ为单位函数。
第二项为保真项,我们给出了相异性的定义:
F ‾ ( x , y ) = t ( x , y ) · q ( x , y ) 1 - q ( x , y ) ;
其中, q ( x , y ) = | s ( x , y ) - v ( x , y ) | | s ( x , y ) | ;
t ( x , y ) = | s ( x , y ) - m y ( x , y ) ϵ + s ( x , y ) - m x ( x , y ) | ;
这里,t(x,y)为方向相关的对比度,mx(x,y)和my(x,y)分别为与X和Y方向相关的均值函数。
步骤1:对子剖面进行光滑处理,经过光滑处理后,突出了数据中蕴含的结构特征,在此基础上,进行后续的主结构提取。
步骤2:主结构提取;主结构指的是每一个子剖面中具有显著视觉特征的区域表示。采用基于FCM(fuzzyc-means)的方法获取子剖面的主结构。
步骤3:主题映射:
主题映射的目的是将分布在各个子剖面的主结构汇聚起来,形成多尺度的主题映射。主题映射定义为:
min x D ( x , f ) , s . t . L ( x , g ) = Σ i = 1 k A i z i ;
x为主题映射,D(·)表示x与f之间的欧式距离,z为子剖面的主题结构,A为权重矩阵。g为可变尺度,L(·)表示各个子剖面在尺度g下到x的集成变换。
图1为可变尺度g取值为256、128、64、32、16、8、4、2、1时的可视化结果示意图。
本发明提出了主模式,主模式将多个子剖面进行合成处理,目的是将数据中隐含的特征呈现出来。主要有以下几个方面:
光滑处理:提出了针对性的光滑模型。
主结构提取:采用FCM方法进行主结构提取。
主题映射:给出了主题映射的概念。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。

Claims (3)

1.一种探地雷达数据的可视方法,其特征在于按照以下步骤进行:
步骤1:建立光滑模型,对子剖面进行光滑处理,经过光滑处理后,突出了数据中蕴含的结构特征,光滑模型表示为:
inf v { T ( v ) = ∫ | ▿ v | + λ F ‾ 2 } ;
v为光滑近似,T(v)为基于参数v的子剖面,λ为调节参数,为保真项,用来控制模型的拟合度,第一项|▽v|为正则项;
步骤2:主结构提取:
采用基于FCM的方法获取子剖面的主结构;
步骤3:主题映射;将分布在各个子剖面的主结构汇聚起来,形成多尺度的主题映射,主题映射定义为:
m i n x D ( x , f ) , s . t . L ( x , g ) = Σ i = 1 k A i z i ;
x为主题映射,D(·)表示x与f之间的欧式距离,z为子剖面的主题结构,A为权重矩阵,g为可变尺度,L(·)表示各个子剖面在尺度g下到x的集成变换。
2.按照权利要求1所述一种探地雷达数据的可视方法,其特征在于:
所述步骤1中|▽v|的计算:
|▽v(x,y)|=(|▽vx(x,y)|2+|▽vy(x,y)|2)1/2
其中,
| ▿ x ( x , y ) | = | ∂ x R ( x , y ) | ;
| ▿ y ( x , y ) | = | ∂ y C ( x , y ) | ;
这里,R(x,y)和C(x,y)分别为沿X轴和Y轴的变换函数,即:
R ( x , y ) = 1 M Σ i , j ∈ l ( x , y ) δ log ( v ( i , j ) p ( i , j ) ) ;
C ( x , y ) = 1 N Σ i , j ∈ h ( x , y ) G ( k , σ ) * v ( i , j ) ;
l(x,y)和h(x,y)分别为沿X轴和Y轴的邻域,G(k,σ)为高斯核函数,k和σ为光滑度控制参数,M和N分别为沿X轴和Y轴的邻域范围,i和j为邻域l(x,y)和h(x,y)的元素标识,δ为单位函数。
3.按照权利要求1所述一种探地雷达数据的可视方法,其特征在于:
所述步骤1中为保真项:
F ‾ ( x , y ) = t ( x , y ) · q ( x , y ) 1 - q ( x , y ) ;
其中, q ( x , y ) = | s ( x , y ) - v ( x , y ) | | s ( x , y ) | ;
t ( x , y ) = | s ( x , y ) - m y ( x , y ) ϵ + s ( x , y ) - m x ( x , y ) | ;
t(x,y)为方向相关的对比度,mx(x,y)和my(x,y)分别为与X和Y方向相关的均值函数。
CN201510835612.7A 2015-11-25 2015-11-25 一种探地雷达数据的可视方法 Active CN105488141B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510835612.7A CN105488141B (zh) 2015-11-25 2015-11-25 一种探地雷达数据的可视方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510835612.7A CN105488141B (zh) 2015-11-25 2015-11-25 一种探地雷达数据的可视方法

Publications (2)

Publication Number Publication Date
CN105488141A true CN105488141A (zh) 2016-04-13
CN105488141B CN105488141B (zh) 2018-02-06

Family

ID=55675116

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510835612.7A Active CN105488141B (zh) 2015-11-25 2015-11-25 一种探地雷达数据的可视方法

Country Status (1)

Country Link
CN (1) CN105488141B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106022339A (zh) * 2016-05-17 2016-10-12 华北水利水电大学 一种复垦土地浅埋地埋管深度的提取方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080079723A1 (en) * 2006-05-16 2008-04-03 David Hanson System and method for visualizing multiple-sensor subsurface imaging data
CN103576127A (zh) * 2012-07-18 2014-02-12 地球物理测勘系统有限公司 用于多天线的地面穿透雷达的合并显示
CN104766090A (zh) * 2015-03-17 2015-07-08 山东工商学院 一种基于bemd和sofm的探地雷达数据可视化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080079723A1 (en) * 2006-05-16 2008-04-03 David Hanson System and method for visualizing multiple-sensor subsurface imaging data
CN103576127A (zh) * 2012-07-18 2014-02-12 地球物理测勘系统有限公司 用于多天线的地面穿透雷达的合并显示
CN104766090A (zh) * 2015-03-17 2015-07-08 山东工商学院 一种基于bemd和sofm的探地雷达数据可视化方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106022339A (zh) * 2016-05-17 2016-10-12 华北水利水电大学 一种复垦土地浅埋地埋管深度的提取方法
CN106022339B (zh) * 2016-05-17 2019-05-24 华北水利水电大学 一种复垦土地浅埋地埋管深度的提取方法

Also Published As

Publication number Publication date
CN105488141B (zh) 2018-02-06

Similar Documents

Publication Publication Date Title
CN102047294B (zh) 对表示网格体的几何不规则网格数据的三维体成像的系统及方法
CN103839292B (zh) 对等值面进行采样并生成高质量三角网格的方法
TWI309027B (en) Computer implemented method and computer system
Ruzinoor et al. A review on 3D terrain visualization of GIS data: techniques and software
CN105279790A (zh) 裂缝网络三维数字岩心建模方法
US20140333616A1 (en) Method for representing terrain, method for creating terrain primitives, and apparatus using the methods
CN104200452A (zh) 基于谱图小波变换的红外与可见光图像融合方法及其装置
CN103323846A (zh) 一种基于极化干涉合成孔径雷达的反演方法及装置
CN104517284A (zh) 基于深度置信网的极化sar图像分割
CN103616494A (zh) 一种模拟岩心孔隙空间流体分布的方法及系统
CN104766090B (zh) 一种基于bemd和sofm的探地雷达数据可视化方法
CN105631932A (zh) 一种轮廓线指导的三维模型重构造方法
CN106599053A (zh) 一种三维模型检索方法
CN113296152A (zh) 断层检测方法及装置
CN109407149A (zh) 基于Hessian矩阵的地震相干数据裂缝检测方法
Zhao et al. Regional W‐phase source inversion for moderate to large earthquakes in China and neighboring areas
US10969507B2 (en) Seismic channel geobody extraction for basin fairway simulation
CN105488141B (zh) 一种探地雷达数据的可视方法
CN110660051B (zh) 一种基于导航金字塔的张量投票处理方法
CN105372718B (zh) 一种2d位场构建3d位场的方法
CN107589449A (zh) 基于曲线Gabor滤波的三维数据断层增强方法
CN107356963A (zh) 一种数据驱动的自适应的地震信号相干体属性分析方法
Waelder An application of the fuzzy theory in surface interpolation and surface deformation analysis
CN106846481A (zh) 一种地质剖面图的生成方法
CN105447314B (zh) 一种探地雷达(gpr)数据分解的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant