CN105488141A - 一种探地雷达数据的多尺度可视方法 - Google Patents
一种探地雷达数据的多尺度可视方法 Download PDFInfo
- 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
Links
- 230000000149 penetrating effect Effects 0.000 title claims abstract description 9
- 238000007794 visualization technique Methods 0.000 title abstract description 5
- 238000009499 grossing Methods 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 10
- 230000000007 visual effect Effects 0.000 claims abstract description 6
- 238000013507 mapping Methods 0.000 claims description 11
- 238000000605 extraction Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 3
- 238000012800 visualization Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/904—Browsing; 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:建立光滑模型,对子剖面进行光滑处理,经过光滑处理后,突出了数据中蕴含的结构特征,光滑模型表示为:
v为光滑近似,T(v)为基于参数v的子剖面,λ为调节参数,为保真项,用来控制模型的拟合度,第一项为正则项;
步骤2:主结构提取:
采用基于FCM的方法获取子剖面的主结构;
步骤3:主题映射;将分布在各个子剖面的主结构汇聚起来,形成多尺度的主题映射,主题映射定义为:
x为主题映射,D(·)表示x与f之间的欧式距离,z为子剖面的主题结构,A为权重矩阵,g为可变尺度,L(·)表示各个子剖面在尺度g下到x的集成变换。
进一步,所述步骤1中的计算:
其中,
这里,R(x,y)和C(x,y)分别为沿X轴和Y轴的变换函数,即
l(x,y)和h(x,y)分别为沿X轴和Y轴的邻域,G(k,σ)为高斯核函数,k和σ为光滑度控制参数,M和N分别为沿X轴和Y轴的邻域范围,i和j为邻域l(x,y)和h(x,y)的元素标识,δ为单位函数。
进一步,所述步骤1中为保真项:
其中,
t(x,y)为方向相关的对比度,mx(x,y)和my(x,y)分别为与X和Y方向相关的均值函数。
本发明的有益效果是将探地雷达数据中不明确的结构特征以多尺度的形式呈现出来,有规律地增强数据的可视深度和广度,帮助数据分析人员更深入地了解结构特征存在的可能性,为综合分析提供技术支持。
附图说明
图1是可变尺度g取值为256至1的可视化结果示意图。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
对原始GPR数据进行分解,形成多个独立的窄带子剖面。
步骤1:建立光滑模型,通过纹理光滑获得显著性区域,这些显著性区域经过融合变换,可以得到主模式。
光滑模型表示为:
这里,v为光滑近似,T(v)为基于参数v的子剖面。λ为调节参数。为保真项,用来控制模型的拟合度。第一项为正则项。由于GPR数据具有异步属性,我们给出以下定义计算
其中,
R(x,y)和C(x,y)分别为沿X轴和Y轴的变换函数,即:
l(x,y)和h(x,y)分别为沿X轴和Y轴的邻域,G(k,σ)为高斯核函数,k和σ为光滑度控制参数。M和N分别为沿X轴和Y轴的邻域范围。i和j为邻域l(x,y)和h(x,y)的元素标识。δ为单位函数。
第二项为保真项,我们给出了相异性的定义:
其中,
这里,t(x,y)为方向相关的对比度,mx(x,y)和my(x,y)分别为与X和Y方向相关的均值函数。
步骤1:对子剖面进行光滑处理,经过光滑处理后,突出了数据中蕴含的结构特征,在此基础上,进行后续的主结构提取。
步骤2:主结构提取;主结构指的是每一个子剖面中具有显著视觉特征的区域表示。采用基于FCM(fuzzyc-means)的方法获取子剖面的主结构。
步骤3:主题映射:
主题映射的目的是将分布在各个子剖面的主结构汇聚起来,形成多尺度的主题映射。主题映射定义为:
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:建立光滑模型,对子剖面进行光滑处理,经过光滑处理后,突出了数据中蕴含的结构特征,光滑模型表示为:
v为光滑近似,T(v)为基于参数v的子剖面,λ为调节参数,为保真项,用来控制模型的拟合度,第一项|▽v|为正则项;
步骤2:主结构提取:
采用基于FCM的方法获取子剖面的主结构;
步骤3:主题映射;将分布在各个子剖面的主结构汇聚起来,形成多尺度的主题映射,主题映射定义为:
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;
其中,
这里,R(x,y)和C(x,y)分别为沿X轴和Y轴的变换函数,即:
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中为保真项:
其中,
t(x,y)为方向相关的对比度,mx(x,y)和my(x,y)分别为与X和Y方向相关的均值函数。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106022339A (zh) * | 2016-05-17 | 2016-10-12 | 华北水利水电大学 | 一种复垦土地浅埋地埋管深度的提取方法 |
Citations (3)
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的探地雷达数据可视化方法 |
-
2015
- 2015-11-25 CN CN201510835612.7A patent/CN105488141B/zh active Active
Patent Citations (3)
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)
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 |