CN105005042B - 一种探地雷达地下目标定位方法 - Google Patents

Abstract

本发明涉及一种探地雷达地下目标定位方法，(1)、将探地雷达B‑Scan数据进行多尺度单演信号分析提取信号的振幅分量；(2)、从振幅分量图像中找出目标感兴趣区域；(3)、对目标感兴趣区域的振幅图像进行边缘提取；(4)、利用霍夫变换进行目标双曲线的定位；(5)、根据目标双曲线实现目标的定位。该方法不需要进行数据的训练，容易进行在线检测；能够较完整保留目标信息，从而提升目标定位的精度，对于金属管线和非金属管线都可能存在的浅层探测而言,不易产生虚警,所以不易漏掉非金属管线目标；能够提升杂波抑制效果，有效区分较强的杂波和目标回波；基于的算法较为简单，算法运算时间短，保证了快速进行目标定位。

Description

一种探地雷达地下目标定位方法

技术领域

本发明涉及一种探地雷达地下目标定位方法，属于探地雷达探测领域。

背景技术

探地雷达是近几十年迅速发展起来的一种有效的浅层地下目标探测技术，它是一种非破坏性探测手段，具有探测速度快、分辨率高、操作方便灵活、探测成本低等诸多优点，已被广泛应用于地下目标，如空洞、管道、地雷等的探测及定位。

探地雷达探测的二维回波数据称为B-Scan数据，它是后续雷达信号处理、目标识别及解译的数据基础，探地雷达目标定位技术也要基于B-Scan数据。对实现目标准确定位影响最大的是探地雷达B-Scan数据中的“杂波”。探地雷达杂波可看作是除了目标回波以外的各种回波,通常包括天线直达波、地表回波、地下非均匀介质产生的回波、以及伪目标所产生的回波等等。探地雷达杂波使得对地下目标的准确探测变得困难，尤其对于浅层埋地目标，目标回波与地表回波相比是较弱的成分，并且目标回波与地表回波间的时延很小，目标回波易被地表强回波这类杂波所淹没。因此有效的抑制杂波方法是实现探地雷达目标准确定位的首要任务。

目前国内外常见的定位方法主要是基于B-Scan图像的双曲线提取，根据提取到的双曲线进行速度估计再计算目标深度。目前常用的地下目标定位方法主要有以下五种：1、基于神经网络对双曲线的提取，需要较多的数据进行训练，不易实现在线检测；2、釆用模糊聚类的模式识别方法，对于金属管线和非金属管线都可能存在的浅层探测而言，容易产生虚警，且容易漏掉非金属管线目标；3、基于图像分割和霍夫变换的方法，应用在浅层探测管线的时候，不能有效区分较强的杂波和目标回波；4、基于图像分割和模板匹配的方法应用在浅层探测管线时候，由于管径的大小可能多变，从而对应的模版也较多，导致算法运算时间较长；5、基于形态学的曲线检测，是根据图像的灰度值进行检测判断，能够判断目标的区域但是得到是多根曲线，进行下一步计算还需对曲线进行处理。

另外，单演小波(Monogenic)信号是一维解析信号二维延拓的结果，该方法是将实系数的小波基和它复值Riesz变换结合起来构成多尺度单演信号分析，从而使得每个小波系数都有振幅，相位和方向信息。

发明内容

本发明的目的是提供一种探地雷达地下目标定位方法，用以解决传统的地下目标定位方法存在着诸多弊端的问题。

为实现上述目的，本发明的方案包括一种探地雷达地下目标定位方法，包括以下步骤：

(1)、将探地雷达B-Scan数据进行多尺度单演信号分析提取信号的振幅分量；

(2)、从所述振幅分量图像中找出目标感兴趣区域；

(3)、对所述目标感兴趣区域的振幅图像进行边缘提取；

(4)、利用霍夫变换进行目标双曲线的定位；

(5)、根据所述目标双曲线实现目标的定位。

所述步骤(1)具体为：

1)、构造小波域单演信号分析框架：假设f是原始信号，解析小波ψ_i(x)，得到的小波系数ω_i[k]：

其中，r_1，i[k]是Riesz变换的实部，r_2，i[k]是Riesz变换的虚部；

2)、基于张量的方法得到方向向量为：u＝(cosθ，sinθ)，其解析变换形式为：

q_i[k]＝r_1，i[k]cosθ+r_2，i[k]sinθ或者

3)、计算幅值和相位：

所述步骤2)具体为：

1)、针对每一道A-Scan数据，利用能量和方差两个统计量，选取能量和方差都为峰值的区域，并结合第一设定阈值，确定出目标所在的深度范围，用时间窗口表示为：t_a～t_b；

2)、让每道A-scan数据中除t_a～t_b时间窗口之外的数据全为0，然后根据A-scan数据的能量变化情况曲线，选择曲线中最大值处的位置，并结合第二设定阈值，确定目标的ROI位置。

所述步骤(3)具体为：

1)、采用高斯滤波进行去噪和剔除虚假目标；

2)、利用Canny算子将所述目标感兴趣区域的正负图像进行边缘提取转换为二值图像。

所述步骤(4)具体为：

1)、针对目标感兴趣区域的边缘提取的结果，将每个边缘像素点(x,t)都定义为目标的中心点(x₀,t₀)；

2)、将每个所述边缘像素点代入公式：中，结合不同的速度值v变换到霍夫平面；

3)、霍夫平面内的出现频数最大的三个参数：x₀,t₀及v对应的值即为该双曲线的顶点坐标(x₀,t₀)及波速v。

所述目标的定位包括：

目标深度定位：H＝t₀·v，其中，t₀为双曲线顶点的反射回波时延，v为电磁波在介质中的传播速度，H为目标的深度位置；

目标水平位置定位：S＝x₀·Δx，x₀为双曲线顶点的横坐标，Δx为天线的移动步长，S为目标的水平位置。

首先，本发明提供的探地雷达地下目标定位方法，不需要进行数据的训练，所以容易进行在线检测；能够较完整保留目标信息，从而提升目标定位的精度，对于金属管线和非金属管线都可能存在的浅层探测而言,不易产生虚警,所以不易漏掉非金属管线目标；能够提升杂波抑制效果，有效区分较强的杂波和目标回波；基于的算法较为简单，算法运算时间短，保证了快速进行目标定位。

附图说明

图1是本发明提供的探地雷达地下目标定位方法的流程图；

图2是探地雷达实测的B-Scan回波图；

图3是利用多尺度单演信号提取振幅后的图像；

图4是采用霍夫变换进行目标双曲线的定位的原理图；

图5是提取到的单演振幅图像依次经边缘检测和霍夫变换的处理后产生的图像；

图6是霍夫平面上的各个曲线图。

具体实施方式

下面结合附图对本发明做进一步详细的说明。

单演小波(Monogenic)信号是一维解析信号二维延拓的结果，该方法是将实系数的小波基和它复值Riesz变换结合起来构成多尺度单演信号分析，从而使得每个小波系数都有振幅，相位和方向信息，因此，采用单演小波进行目标回波振幅特征的提取，能够在较完整保留目标信息的同时提升杂波抑制效果，从而提升目标定位的精度。

本发明是一种基于多尺度单演(Monogenic)信号分析的方法来进行探地雷达地下目标的定位，如图1所示，具体如下：

1、首先将图2所示的探地雷达回波B-Scan图像进行多尺度单演信号分析提取信号的振幅分量。具体包含以下步骤：

1.1、针对探地雷达的B-Scan探测回波数据f(χ)，其中χ＝(x，y)，求其Riesz变换：

其中即表示f(χ)的Riesz变换，为Riesz变换的实部，为Riesz变换的虚部。

1.2、构造小波域单演信号分析框架。假设f是原始信号，解析小波ψ_i(x)，得到的小波系数ω_i[k]：

其中，r_1，i[k]是ω_i[k]的Riesz变换的实部，r_2，i[k]是ω_i[k]的Riesz变换的虚部。

1.3、基于张量的方法得到方向向量为：

u＝(cosθ，sinθ) (4)

其相应的解析变换形式可以用下式表示：

q_i[k]＝r_1，i[k]cosθ+r_2，i[k]sinθ (5)

或者：

1.5、计算得到的幅度和相位：

从而实现了在小波子频带中由小波振幅提取出了探地雷达B-Scan探测数据中的幅度分量，如图3所示。

2、采用基于能量统计的目标感兴趣区域提取方法从上述幅度分量图像中提取目标区域。具体包括以下两个步骤：

2.1获取时间窗口：

探地雷达二维B-Scan回波图像是由多道一维A-Scan回波数据组成的。在步骤1中提取出的B-Scan图像振幅分量的基础上，针对每一道A-Scan数据，根据其在时间t方向的特点可知，在非目标区域A-Scan在t方向的方差较小，能量较小；在目标区域，A-scan的方差较大，能量较大。利用能量和方差两个统计量，选取能量和方差都为峰值的区域，并结合合适的阈值，确定出目标所在的深度范围，用时间表示t_a～t_b，即时间窗口。(其中，本段中的阀值可由人工预设或迭代计算方式获得，这是常规技术，这里不做赘述)

2.2加窗能量法提取ROI：

选择时间窗口后，对A-scan数据进行归一化处理，让每道A-scan中除t_a～t_b之外的数据全为0。加窗后，在测线方向上每道A-scan数据都减少了许多干扰信号，剩下的信号成分中如果无目标，则每个采样点的幅值接近于零；如有目标，则有较大幅值。由此，可以根据加窗后每个A-scan数据的能量变化情况，得到对应的能量变化情况曲线，选择能量曲线中最大值处的位置，并结合阈值的选取确定目标的ROI位置。(其中，本段中的阀值也可由人工预设或迭代计算方式获得，这是常规技术，这里不做赘述)

上述为本发明提供的一种ROI提取方法，现有技术中的常用的ROI提取方法有时频分析、神经网络、模板匹配法。时频分析方法对信噪比不高的区域，特别是有相关频率的干扰的区域容易产生误判；神经网络法需要有大的已知数据做训练，其计算速度也不能很好地适用于现场探测；模板匹配法需要对测量数据有很多的先验知识，因此限制了它的应用。

3、利用Canny算子针对上述ROI区域的振幅图像进行边缘提取。具体包含以下两个步骤：

3.1、采用高斯滤波对ROI区域的振幅图像进行去噪和剔除虚假目标；

3.2、利用Canny算子将已通过高斯滤波处理后的ROI区域的振幅图像进行边缘提取转换为二值图像。

4、采用霍夫(Hough)变换进行目标双曲线的定位。

由探地雷达原理知：

其中，x表示天线位置，x₀表示目标的水平位置，v表示电磁波在介质中的传播速度，t₀表示目标反射回波时延。因此，求出顶点坐标(x₀,t₀)以及波速v即可准确定位目标，如图4所示。

采用霍夫(Hough)变换进行目标双曲线的定位的基本原理是将图像空间中的曲线变换到参数空间中,通过检测参数空间中的极值点,确定出该曲线的描述参数,从而提取图像中的规则曲线。具体步骤如下：

4.1、首先针对ROI区域的Canny边缘检测结果，将每个边缘像素点(x,t)都定义为目标的中心点(x₀,t₀)；

4.2、根据实际探测介质的电特性，给出速度v的大致取值范围；

4.3、将每个边缘像素点(x,t)代入公式(9)，在保证t₀为大于零的实数的情况下，结合不同的速度值v变换到霍夫平面；

4.4、用4.3中的方法遍历各条Canny边缘图像，统计各组Canny边缘图像中的三个参数：x₀,t₀及v出现的频数，每组参数中出现频数最大的值即为该双曲线的顶点坐标(x₀,t₀)及波速v，如图5和6所示。

5、实现目标的定位。具体包含以下两个步骤：

5.1、目标深度定位：根据探地雷达原理知，双曲线顶点的纵坐标表示最短回波时延，即在这一测点探地雷达距离目标最近，因此双曲线的纵坐标就代表目标对应的深度。将提取到的双曲线顶点的反射回波时延乘以电磁波在介质中的传播速度v，就得到目标的深度位置H，即：

H＝t₀·v (10)

5.2、目标水平位置定位：同理，双曲线顶点的横坐标即代表目标对应的水平位置，因此，将提取到的双曲线顶点的横坐标乘以天线的移动步长Δx，就得到目标的水平位置S，即：

S＝x₀·Δx (11)

采用本发明提供的定位方法的效果如图6所示。

说明书中未详细说明的部分属于本领域技术人员公知的现有技术。

以上给出了具体的实施方式，但本发明不局限于所描述的实施方式。本发明的基本思路在于上述基本方案，对本领域普通技术人员而言，根据本发明的教导，设计出各种变形的模型、公式、参数并不需要花费创造性劳动。在不脱离本发明的原理和精神的情况下对实施方式进行的变化、修改、替换和变型仍落入本发明的保护范围内。

Claims (1)

1.一种探地雷达地下目标定位方法，其特征在于，所述定位方法包括以下步骤：

(1)、将探地雷达B-Scan数据进行多尺度单演信号分析提取信号的振幅分量；所述步骤(1)具体为：1)、构造小波域单演信号分析框架：假设f是原始信号，解析小波ψ_i(x)，得到的小波系数ω_i[k]：其中，r_1，i[k]是Riesz变换的实部，r_2，i[k]是Riesz变换的虚部；2)、基于张量的方法得到方向向量为：u＝(cosθ，sinθ)，其解析变换形式为：q_i[k]＝r_1，i[k]cosθ+r_2，i[k]sinθ或者3)、计算幅值和相位：

(2)、从所述振幅分量图像中找出目标感兴趣区域；所述步骤2)具体为：1)、针对每一道A-Scan数据，利用能量和方差两个统计量，选取能量和方差都为峰值的区域，并结合第一设定阈值，确定出目标所在的深度范围，用时间窗口表示为：t_a～t_b；2)、让每道A-scan数据中除t_a～t_b时间窗口之外的数据全为0，然后根据A-scan数据的能量变化情况曲线，选择曲线中最大值处的位置，并结合第二设定阈值，确定目标的ROI位置；

(3)、对所述目标感兴趣区域的振幅图像进行边缘提取；所述步骤(3)具体为：1)、采用高斯滤波进行去噪和剔除虚假目标；2)、利用Canny算子将所述目标感兴趣区域的正负图像进行边缘提取转换为二值图像；

(4)、利用霍夫变换进行目标双曲线的定位；所述步骤(4)具体为：1)、针对目标感兴趣区域的边缘提取的结果，将每个边缘像素点(x,t)都定义为目标的中心点(x₀,t₀)；2)、将每个所述边缘像素点代入公式：中，结合不同的速度值v变换到霍夫平面；3)、霍夫平面内的出现频数最大的三个参数：x₀,t₀及v对应的值即为该双曲线的顶点坐标(x₀,t₀)及波速v；

(5)、根据所述目标双曲线实现目标的定位；所述目标的定位包括：目标深度定位：H＝t₀·v，其中，t₀为双曲线顶点的反射回波时延，v为电磁波在介质中的传播速度，H为目标的深度位置；目标水平位置定位：S＝x₀·Δx，x₀为双曲线顶点的横坐标，Δx为天线的移动步长，S为目标的水平位置。