一种探地雷达地下目标位置检测方法
【技术领域】
[0001] 本发明涉及一种探地雷达地下目标位置检测方法，属于探地雷达探测领域。
【背景技术】
[0002] 探地雷达是近几十年迅速发展起来的一种有效的浅层地下目标探测技术，它是一 种非破坏性探测手段，具有探测速度快、分辨率高、操作方便灵活、探测成本低等诸多优点， 已被广泛应用于地下目标，如空洞、管道、地雷等的探测及定位。
[0003] 探地雷达探测的二维回波数据称为B-Scan数据，它是后续雷达信号处理、目标识 别及解译的数据基础，探地雷达目标定位技术也要基于B-Scan数据。对实现目标准确定位 影响最大的是探地雷达B-Scan数据中的"杂波"。探地雷达杂波可看作是除了目标回波以 外的各种回波，通常包括天线直达波、地表回波、地下非均匀介质产生的回波、以及伪目标 所产生的回波等等。探地雷达杂波使得对地下目标的准确探测变得困难，尤其对于浅层埋 地目标，目标回波与地表回波相比是较弱的成分，并且目标回波与地表回波间的时延很小， 目标回波易被地表强回波这类杂波所淹没。因此有效的抑制杂波方法是实现探地雷达目标 准确定位的首要任务。
[0004] 目前国内外常见的定位方法主要是基于B-Scan图像的双曲线提取，根据提取到 的双曲线进行速度估计再计算目标深度。目前常用的地下目标定位方法主要有以下五种： 1、基于神经网络对双曲线的提取，需要较多的数据进行训练，不易实现在线检测；2、釆用模 糊聚类的模式识别方法，对于金属管线和非金属管线都可能存在的浅层探测而言，容易产 生虚警，且容易漏掉非金属管线目标；3、基于图像分割和霍夫变换的方法，应用在浅层探测 管线的时候，不能有效区分较强的杂波和目标回波；4、基于图像分割和模板匹配的方法应 用在浅层探测管线时候，由于管径的大小可能多变，从而对应的模版也较多，导致算法运算 时间较长；5、基于形态学的曲线检测，是根据图像的灰度值进行检测判断，能够判断目标的 区域但是得到是多根曲线，进行下一步计算还需对曲线进行处理。
[0005] 另外，单演小波（Monogenic)信号是一维解析信号二维延拓的结果，该方法是将 实系数的小波基和它复值Riesz变换结合起来构成多尺度单演信号分析，从而使得每个小 波系数都有振幅，相位和方向信息。
【发明内容】
[0006] 本发明的目的是提供一种探地雷达地下目标位置检测方法，用以解决传统的地下 目标定位方法存在着诸多弊端的问题。
[0007] 为实现上述目的，本发明的方案包括一种探地雷达地下目标位置检测方法，所述 位置检测方法包括以下步骤：
[0008] (1)、将探地雷达B-Scan数据进行多层分解提取，提取多尺度单演信号调制特征 分量，包括振幅、相位和方向特征分量；
[0009] (2)、提取第一个尺度，即从多层分解中选择第一层的振幅、相位和方向特征分量， 利用HSB通道表示方法结合成彩色图像；
[0010](3)、从振幅特征分量中找出目标感兴趣区域；
[0011] (4)、对目标感兴趣区域内的彩色图像进行灰度化，并对所述目标感兴趣区域的振 幅图像进行边缘提取；
[0012](5)、进行目标双曲线的定位；
[0013] (6)、根据所述目标双曲线实现目标的位置确定。
[0014] 所述步骤⑴具体为：
[0015] 1)、构造小波域单演信号分析框架：假设f是原始信号，解析小波杓〇〇,得到的 小波系数ω^]:
[0017]
:，其中，rjk]是Riesz变换 的实部，r2,Jk]是Riesz变换的虚部；
[0018] 2)、计算局部结构矩阵：
[0019]
其中，n，m= 1，2rjk]和r2[k]分别是 Riesz变换的实部和虚部；
[0020] 基于张量的方法得到方向向量为：
[0022] 根据选定的方向，相应的方向Hilbert变换是：
[0023]qi[k]= r!,; [k]cos Θ +r2,; [k]sin Θ或者
[0024]3)、计算振幅特征分量和相位特征分量：
[0026] 所述步骤（2)中，所述利用HSB通道表示方法结合成彩色图像具体为：所述彩色图 像中的色度分量用方向特征分量表不，饱和度分量用相位特征分量表不，亮度分量用相位 特征分量表示。
[0027] 所述步骤⑶具体为：
[0028] 1)、针对每一道A-Scan数据，利用能量和方差两个统计量，选取能量和方差都为 峰值的区域，并结合第一设定阈值，确定出目标所在的深度范围，用时间窗口表示为：ta~ tb；
[0029] 2)、让每道A-scan数据中除ta~、时间窗口之外的数据全为0,然后根据A-scan 数据的能量变化情况曲线，选择曲线中最大值处的位置，并结合第二设定阈值，确定目标的 目标感兴趣区域。
[0030] 所述步骤⑷具体为：
[0031] 1)、对目标感兴趣区域内的彩色图像进行灰度化；
[0032]2)、采用高斯滤波进行去噪和剔除虚假目标；
[0033]3)、利用Canny算子将所述目标感兴趣区域的正负图像进行边缘提取转换为二值 图像。
[0034] 所述步骤（5)包括以下两个步骤：
[0035]1)、利用最小二乘法估计目标双曲线的顶点坐标；
[0036] 2)、根据估计的目标双曲线的顶点坐标，采用霍夫变换来求取波速，进而得到目标 双曲线的定位。
[0037] 所述步骤（5)中的步骤1)具体为：设定所述二值图像中值为1的点为待拟合的 点，优化条件是：
[0038]
[0039] 其中，（X。，t。）是估计系数的向量，\和ti是二值图像中η个值为1的点 数据行向量的横纵坐标值，i= 1，…，n，f(X。，t。，X)是
的变形，
[0040] 然后采用对多元函数求极值的方法，对X。求偏导，让偏导等于0;对t。求 偏导，让偏导等于0;进而得到关于参数（X(]，t。）的方程，通过解所述关于参数（X(]，t。）的方 程得到拟合的曲线；
[0041] 最后，取该拟合好的曲线的纵坐标的极大值点为（X。，t。）中的坐标X。，同时得到坐 标t〇;
[0042]所述步骤（5)中的步骤2)具体为：根据估计出的顶点坐标（X。，t。)，采用一维的霍 夫变换求取波速v:根据探地区域设定传播速度v在一定范围内变换，在此范围内变换参数 v的值，此时参数空间将有聚焦点，聚焦点的集中的值即为波速V。
[0043] 所述目标的定位包括：
[0044]目标深度定位：H=t。·V，其中，t。为双曲线顶点的反射回波时延，v为电磁波在 介质中的传播速度，Η为目标的深度位置；
[0045]目标水平位置定位：S=X。·ΔX，X。为双曲线顶点的横坐标，ΔX为天线的移动步 长，S为目标的水平位置。
[0046] 首先，本发明提供的探地雷达地下目标位置检测方法，不需要进行数据的训练，所 以容易进行在线检测；能够较完整保留目标信息，从而提升目标定位的精度，对于金属管线 和非金属管线都可能存在的浅层探测而言，不易产生虚警，所以不易漏掉非金属管线目 标；能够提升杂波抑制效果，有效区分较强的杂波和目标回波；基于的算法较为简单，算法 运算时间短，保证了快速进行目标定位。
【附图说明】
[0047] 图1是本发明提供的探地雷达地下目标位置检测方法的流程图；
[0048] 图2是探地雷达实测的B-Scan回波图；
[0049] 图3是彩色图像二值化和R0I算法限定区域后得到的图；
[0050] 图4是采用雷达天线与目标B-Scan回波的几何位置关系图。
【具体实施方式】
[0051] 下面结合附图对本发明做进一步详细的说明。
[0052] 单演小波（Monogenic)信号是一维解析信号二维延拓的结果，该方法是将实系数 的小波基和它复值Riesz变换结合起来构成多尺度单演信号分析，从而使得每个小波系数 都有振幅，相位和方向信息，因此，采用单演小波进行目标回波振幅特征的提取，能够在较 完整保留目标?目息的同时提升杂波抑制效果，从而提升目标定位的精度。
[0053] 本发明是一种基于多尺度单演（Monogenic)信号分析的方法来进行探地雷达地 下目标的位置检测和确定，如图1所示，具体如下：
[0054] 1、首先将图2所示的探地雷达回波B-Scan图像进行多层分解提取，提取多尺度单 演信号调制特征分量，包括振幅、相位和方向特征分量。具体包含以下步骤：
[0055] 1. 1、针对探地雷达的B-Scan探测回波数据f(X)，其中X= (X，y)，求其Riesz 变换：
[0056]
[0057] 其中货即表示f(x)的Riesz变换，= 为Riesz变换的实部， = //?/'(7))为Riesz变换的虚部。
[0058] 1.2、构造小波域单演信号分析框架。假设f是原始信号，解析小波Mx)，得到的 小波系数ω^]:
[0061] 其中，是cojk]的Riesz变换的实部，r2,i[k]是cojk]的Riesz变换的虚 部。
[0062] 1.3、计算局部结构矩阵：
[0064] 其中，n，m= 1，2rjk]和r2[k]分别是Riesz变换的实部和虚部；
[0065] 基于张量的方法得到方向向量为：
[0067] u = (cos θ，sin θ ) (6)
[0068] 根据选定的方向，相应的方向Hilbert变换可以用下式表示：
[0069]q；[k]=rh；[k] cos Θ +r2>；[k] sin Θ (7)
[0070] 或者：
[0072] 1. 4、计算得到的振幅特征分量和相位特征分量：
[0075] 从而实现了将探地雷达B-Scan数据进行多层分解提取，在小波子频带中由小波 振幅提取出多尺度的振幅、相位和方向特征分量。
[0076] 2、从上述多尺度的信息中，即多层分解中，提取第一个尺度的振幅、相位和方向特 征分量，即从多层分解中选择第一层的振幅、相位和方向特征分量，利用HSB通道表示方法 结合成彩色图像，该彩色图像中，色度分量用方向特征分量表示，饱和度分量用相位特征分 量表示，亮度分量用相位特征分量表示。
[0077] 3、采用基于能量统计的目标感兴趣区域提取方法从上述振幅分量图像中提取目 标区域。具体包括以下两个步骤：
[0078] 3.1获取