Research on ultra wideband radar life detection technology
-
摘要: 超宽带（UWB）雷达生命探测技术具有功耗低、穿透性好、保密性高等优点，有利于提高灾后受困人员的生存率。系统总结了UWB雷达生命探测技术的国内外研究进展及现状。根据发射信号形式不同，将UWB雷达生命探测技术分为连续波雷达生命探测技术和脉冲波雷达生命探测技术，分别介绍了2种探测技术的原理与应用优势。基于连续波雷达生命探测技术和脉冲波雷达生命探测技术的各自特点，从探测信号发射、回波信号预处理、生命信号提取与分析3个角度，分析了UWB雷达生命探测关键技术，总结了3种关键技术的研究现状。对UWB雷达生命探测技术的研究进行展望：突破生命探测仪收发机硬件性能，提升发射信号带宽，优化射频功率放大技术，以增大穿墙探测距离；综合利用多种特征提取方法和智能模式分类方法，以及人工智能、大数据、云计算等新一代信息技术，提高目标识别的精确度；研制基于多输入多输出雷达的人体目标辨识与定位装备和高精度分布式组网全极化UWB雷达生命探测仪，提升探测结果维度。Abstract: Ultra wideband (UWB) radar life detection technology has the advantages of low power consumption, good penetration, and high confidentiality. It is beneficial for improving the survival rate of trapped personnel after disasters. This paper systematically summarizes the research progress and current status of UWB radar life detection technology both domestically and internationally. According to the different forms of transmitted signals, UWB radar life detection technology is divided into continuous wave radar life detection technology and pulse wave radar life detection technology. The principles and application advantages of the two detection technologies are introduced respectively. Based on the respective features of continuous wave radar life detection technology and pulse wave radar life detection technology, this paper analyzes the key technologies of UWB radar life detection from three perspectives: detection signal transmission, echo signal preprocessing, and life signal extraction and analysis, and summarizes the research status of the three key technologies. The paper proposes prospects for the research on UWB radar life detection technology. The technology breaks through the hardware performance of life detector transceivers, improves the transmission signal bandwidth, and optimizes RF power amplification technology to increase the detection distance through walls. The technology comprehensively utilizes multiple feature extraction methods and intelligent pattern classification methods, as well as new generation information technologies such as artificial intelligence, big data, and cloud computing, to improve the precision of target recognition. The technology develops a human target recognition and positioning equipment based on multi input multi output radar and a high-precision distributed networked fully polarized UWB radar life detector to enhance the dimension of detection results.
-
0. 引言
煤炭产业是我国能源体系的基石[1]。近年来，随着人工智能、图像处理等技术的蓬勃发展，计算机视觉技术成为一种感知煤矿井下环境的有效手段[2]。尤其是井下发生安全事故时，在救援人员无法进入的情况下，机器人可利用同步定位与地图构建（Simultaneous Localization and Mapping，SLAM）技术，在不依赖先验信息的条件下，估计出自身位姿，并构建出坍塌矿井的地图模型[3]。图像匹配是视觉SLAM技术中极为重要的一环，用于根据图像之间的变换关系确定相机位姿。由于矿井环境特殊，拍摄的图像往往亮度及对比度低[4]。同时，避难硐室、巷道密集管线等图像局部存在较多相似纹理，导致特征的区分度较低。在图像采集时，不可避免存在快速平移与旋转，使得待匹配图像存在较大视差。特殊的光照环境、局部相似纹理的干扰及视差图像本身的缺点共同导致井下图像匹配困难，影响视觉SLAM定位精度与建图结果，因此对井下图像的高质量匹配研究具有重要意义。
目前主流的图像匹配方法包括基于点特征的匹配算法和基于线特征的匹配算法[5]。基于点特征的匹配算法虽然可在井下应用，但是井下光照条件差，巷道图像纹理弱，可能造成特征点的提取数量少、分布不均匀、匹配精度低等[6]。而基于线特征的匹配算法具有较强的鲁棒性和抗噪能力[7]，更加适用于井下图像匹配。
基于线特征的匹配方法主要利用线特征的几何信息[8]及邻域内灰度、梯度的相似性[9]，结合几何约束[10]进行匹配，侧重于搜索范围约束[11]、描述符构建[12]及匹配核验[13]等方面。Li Gang等[14]提出一种基于语义不变量的点线特征匹配方法，为待匹配线特征添加语义约束，降低了线特征的失配率，但该方法依赖语义信息提取网络，难以在井下直接应用。Zheng Xianwei等[15]提出一种基于全局投影变换的线匹配方法，匹配精度较高，但是局限在城市街道等室外场景，在视差变换大时会失效。Wang Qiang等[16]利用位置和方向系统构建出图像间的变换关系来获取匹配线对，消除了跨视角图像间的几何形变，匹配精度较高，但对于宽基线图像匹配仍有提升空间。Shen Liang等[17]提出一种线段误匹配去除方法，解决对线段长度的敏感性和断裂问题，提高了匹配召回率，但是该方法不适用于存在仿射等非刚性变换的图像中。刘肃艳等[18]提出了一种结合线对几何条件约束及单线描述符的方法，该方法性能良好，但涉及的参数较多，依赖手动调参。张珊等[19]提出了一种结合网状描述符和单应约束的近景影像直线匹配算法，该算法能够获得较高的匹配精度，且鲁棒性良好，但是依赖同名特征点约束，在同名特征点稀少的情况下影响匹配结果。随着深度学习的发展，学者们将卷积神经网络（Convolutional Neural Network，CNN）强大的特征提取能力应用在线特征匹配上。A. Vakhitov等[20]提出一种全卷积神经网络构建的可学习型线描述符，用于线特征匹配。M. Lange等[21]提出一种基于机器学习和小波增强的线特征描述符，预处理阶段利用小波变换从图像中提取特征作为描述符。R. Pautrat等[22]首次提出联合线描述与检测的深度学习网络，即自我监督遮挡感知线描述和检测（Self-supervised Occlusion-aware Line Description and Detection，SOLD2），用于在线特征检测与匹配。
基于深度学习的线描述符对线段遮挡等场景具有较高鲁棒性，性能优于传统描述符，但CNN架构的描述符将可变长度线段抽象为固定维进行描述，仍然不利于线段长度及视差变化较大图像的匹配。针对该问题，本文提出一种基于直线段检测法（Line Segment Detector，LSD）和LT（Line Transformers）描述符的矿井图像线特征匹配算法。采用改进单参数同态滤波算法和对比度受限的自适应直方图均衡化（Contrast Limited Adaptive Histogram Equalization，CLAHE）算法进行图像增强，在此基础上进行LSD线特征提取，用LT描述符构建LSD线特征向量[23]，最后利用欧氏距离及最近邻准则进行匹配线对筛选，完成匹配。
1. 算法流程
基于LSD和LT描述符的矿井图像线特征匹配算法流程如图1所示。
![]()
图 1 基于LSD和LT描述符的矿井图像线特征匹配算法流程Figure 1. Process of mine image line feature matching algorithm based on line segment detector(LSD) and line transformers(LT) descriptor在图像的RGB空间，利用改进单参数同态滤波算法进行滤波，在频域降低图像的照射分量，并提高反射分量，增强后图像的亮度与对比度均有提高。在图像的YUV空间，利用CLAHE算法对图像的亮度分量进行均衡[24]，使亮度更加均匀，更好地保留图像的细节信息。再逆变换至RGB空间提取LSD线特征，增强图像后，提取的LSD线特征数量更多。为了克服LSD线特征匹配常用的线二进制描述符（Line Binary Descriptor，LBD）等鲁棒性差、匹配精度低的缺点，采用基于Transformer架构的LT描述符构建LSD线特征向量，并进行L2范数归一化，计算出目标直线与候选直线的特征向量的欧氏距离，结合最小欧氏距离和最近邻准则完成匹配，得到可靠的同名直线。
2. 算法原理
2.1 改进同态滤波算法
同态滤波是一种在频域应用的图像处理算法[25]。根据照射−反射模型，图像密度函数$ f(x,y) $（x，y为像素空间坐标）可表示为自身的照射分量$ i(x,y) $和反射分量$ r(x,y) $的乘积。
$$ f(x,y) = i(x,y)r(x,y) $$ (1) 对式（1）进行对数变换，得
$$ \ln f(x,y) = \ln i(x,y) + \ln r(x,y) $$ (2) 对式（2）进行傅里叶变换，然后应用$ H(u,v) $（u，v为x，y对应的离散频率变量）同态滤波器，得
$$ H(u,v)F(u,v) = H(u,v)I(u,v) + H(u,v)R(u,v) $$ (3) 式中$ F(u,v) $，$ I(u,v) $，$ R(u,v) $分别为$ f(x,y) ， i(x,y) ，r(x,y) $的傅里叶变换。
$ H(u,v) $通常选用高斯型同态滤波传递函数：
$$ H(u,v) = ({P_{\mathrm{H}}} - {P_{\mathrm{L}}})\left[ {1 - {\exp{ \left(- c{{\left(\dfrac{{D(u,v)}}{{2{D_0}}}\right)}^{2n}}\right)}}} \right] + {P_{\mathrm{L}}} $$ (4) $$ D(u,v) = \sqrt {{{(u - {u_0})}^2} + {{(v - {v_0})}^2}} $$ (5) 式中：PH，PL分别为高频增益和低频增益；$ c $为锐化系数；$ D(u,v) $为某频率$ (u,v) $到频率中心$ ({u_0},{v_0}) $的欧氏距离；$ {D_0} $为截止频率；$ n $为滤波器的阶数。
高斯型同态滤波传递函数需要手动调节参数$ {P_{\mathrm{H}}} $，$ {P_{\mathrm{L}}} $，$ c $，$ n $，不利于算法的泛化。因此引入一种单参数同态滤波传递函数[26]：
$$ H(u,v)=\frac{1}{1+D(u\mathit{\mathrm{,}{v}})^{-k}} $$ (6) 式中k为滤波器的调节参数。
式（6）中只存在一个参数k，可降低调参的复杂性。高斯型同态滤波传递函数与单参数同态滤波传递函数的三维结构如图2、图3所示。可看出单参数同态滤波器由中心频率到高频的过渡相较传统高斯型同态滤波器更平缓，斜率更小，因此滤波更加均匀。
使用改进后的单参数同态滤波器对$ F(u,v) $进行滤波后，再用傅里叶逆变换将图像变换至空间域，得
$$ {h_{\mathrm{f}}}(x,y) = {h_{\mathrm{i}}}(x,y) + {h_{\mathrm{r}}}(x,y) $$ (7) 式中$ {h_{\mathrm{f}}}(x,y) $，$ {h_{\mathrm{i}}}(x,y) $，$ {h_{\mathrm{r}}}(x,y) $分别为逆变换后的图像密度函数、照射分量和反射分量。
对式（7）进行指数变换，得到输出图像密度函数：
$$ g(x,y) = {\exp\left({{h_{\mathrm{i}}}(x,y)}\right)}{\exp\left({{h_{\mathrm{r}}}(x,y)}\right)} $$ (8) 采用单参数同态滤波算法在频域对井下图像进行增强，可初步提高井下图像的对比度与亮度。
2.2 CLAHE算法
单参数同态滤波虽然能够提升亮度及对比度，但是仍然存在亮度分布不均匀、局部细节损失的现象，可能会导致线特征断裂、漏检。因此在图像经过同态滤波增强后，采用CLAHE算法对YUV空间的亮度分量进行均衡。均衡后的图像亮度分布更加均匀，细节更明显，可进一步提升线段提取的质量。
CLAHE算法在直方图均衡的基础上，通过抑制部分灰度级被过多合并所引起的噪声放大和局部对比度增强现象，使图像更加清晰[27]。CLAHE算法将输入图像分为若干个大小相等且不重叠的子域，并计算每个子域的直方图。由于算法引入了对比度限制，在计算累积分布函数前，使用预先设定的限制阈值对每个子域直方图进行剪切，从而限制放大幅值，将超出阈值的部分重新均匀分布到其他部分，如图4所示。
![]()
图 4 CLAHE原理Figure 4. Principle of contrast limited adaptive histogram equalization(CLAHE)2.3 LSD线特征向量构建方法
LT描述符是基于Transformer架构的新型描述符，摒弃了传统的CNN和循环神经网络（Recurrent Neural Network，RNN），整个网络融合了自注意力机制。LT描述符仿照自然语言处理（Natural Language Processing，NLP）中的思想，将点视为“单词”，将线段视为“句子”。将NLP的结果作为整条线段的描述符，通过自注意力机制关注线段上关键点来理解线段的上下文，自适应地将各种长度的线抽象为固定大小的描述符。同时，采用了组描述符的思想，将线的几何属性共享到邻域，形成线签名网络，通过在邻域内共享线的相对几何结构，使描述符学习到邻域的几何属性。LT描述符适合应用在具有较大视差变化及较多相似纹理干扰的图像匹配中。
基于LSD和LT描述符的矿井图像线特征匹配算法模型如图5所示。进行LSD线检测，同时使用CNN得到整张图像的密集描述符地图（Descriptor map），包含图像中每个像素的特征点置信度得分与该点描述符。对LSD线均匀取关键点，作为点标记，表示为$ \boldsymbol{p}=[x\; \ y\; \ C] $，其中$ C $为关键点置信度，在密集描述符地图中查找每个关键点对应的描述符向量，实现点标记对应的点嵌入（Point embedding）提取，点嵌入$ {\boldsymbol{E}} \in {{\bf{R}}^{1 \times w}} $，其中$ w $为点描述符的维度。特殊嵌入[LINE]是线描述符的初始状态，表示为$ {{\boldsymbol{E}}_{{\mathrm{line}}}} \in {{\bf{R}}^{(n + 1) \times w}} $，其权值在训练过程中学习。在多层感知器（Multilayer Perceptron，MLP）中利用每个关键点的位置得到位置嵌入（Positional embedding），表示为$ {{\boldsymbol{E}}_{{\mathrm{pos}}}} \in {{\bf{R}}^{(n + 1) \times w}} $。结合构建好的点嵌入与位置嵌入，使用Transformers编码器对线描述符进行建模。
![]()
图 5 基于LSD和LT描述符的矿井图像线特征匹配算法模型Figure 5. Model of mine image line feature matching algorithm based on LSD and LT descriptorTransformer编码器由多头自注意（Multi-head Self-Attention，MSA）层和MLP层组成，具有残差连接与层归一化（Layer Normalization，LN）功能，将Transformer编码器堆叠$ L $次，得
$$ \left\{\begin{split} & {{\boldsymbol{Z}}_0} =\left[ {{{\boldsymbol{E}}_{{\mathrm{line}}}}\;{{\boldsymbol{E}}_1}\;{{\boldsymbol{E}}_2} \;\cdots \;{{\boldsymbol{E}}_{N}}} \right] + {{\boldsymbol{E}}_{{\mathrm{pos}}}} \\ & {\boldsymbol{Z}}_{I - 1}^\prime = {O_{\mathrm{LN}}}(O_{\mathrm{MSA}}({{\boldsymbol{Z}}_{I - 1}},{\boldsymbol{m}}) + {{\boldsymbol{Z}}_{I - 1}}) \\ & {{\boldsymbol{Z}}_I} = {O_{\mathrm{LN}}}({O_{\mathrm{ALP}}}({\boldsymbol{Z}}_{I - 1}^\prime ) + {\boldsymbol{Z}}_{I - 1}^\prime ) \\ & {\boldsymbol{d}} = {\boldsymbol{Z}}_L^0 \end{split}\right. $$ (9) 式中：$ {{\boldsymbol{Z}}_0} $为Transformer编码器的初始输入；$ {\boldsymbol{Z}}_{I - 1}^\prime $为计算第I$(I = 1,2, \cdots, L) $次堆叠时Transformer编码器的输入$ {{\boldsymbol{Z}}_I} $的中间变量；OLN，OMSA，OMLP分别为归一化及通过MSA层、MLP层的操作；$ {{\boldsymbol{m}}} $为掩码向量；d为描述符；$ {\boldsymbol{Z}}_L^0 $为第L次堆叠后Transformer编码器的输出。
当第I−1次堆叠时的Transformer编码器输入$ {{\boldsymbol{Z}}_{I - 1}} $进入MSA层时，为了解决直线特征长度不同的问题，在网络中加入掩码向量m，去除相关性比较低的特征点。
将直线a的中点位置$ ({x_a},{y_a}) $、直线与主方向的角度信息（$ \cos \;{\theta _a},\sin \;{\theta _a} $）、线长$ {l_a} $输入MLP层，得到线属性嵌入，并添加至描述符$ {{\boldsymbol{d}}_a} $，得到含有线属性的描述符$ {\boldsymbol{d}}_a^\prime $。对图像中所有M个线描述符进行构建，得到初始含有线属性的描述符向量集合$ {{\boldsymbol{s}}_0} $。设sI为签名网络结果，将sI−1输入MSA层后得到的计算结果与sI−1进行级联，作为MLP层的输入，堆叠M次后得到sI。
$$ \left\{\begin{split} {{\boldsymbol{d}}}_{a}^{\prime }=&{{\boldsymbol{d}}}_{a}+O_{\mathrm{MLP}}({x}_{a},{y}_{a},{l}_{a},\mathrm{cos}\;{\theta }_{a},\mathrm{sin}\;{\theta }_{a})\\ {{\boldsymbol{s}}}_{0}=&\left[{{\boldsymbol{d}}}_{1}^{\prime }\;{{\boldsymbol{d}}}_{2}^{\prime }\;\cdots\; {{\boldsymbol{d}}}_{M}^{\prime }\right]\\ {{\boldsymbol{s}}}_{I}=&{{\boldsymbol{s}}}_{I-1}+O_{\mathrm{MLP}}({{\boldsymbol{s}}}_{I-1}||O_{\mathrm{MSA}}({{\boldsymbol{s}}}_{I-1}))\quad \end{split}\right. $$ (10) 3. 实验与结果分析
实验采用的编程语言为Python3.7，使用计算机视觉库OpenCV 3.4.18.65，集成开发环境为PyCharm 2021.3.1，GPU为NVIDIA 2080Ti 16 GiB，深度学习框架Pytorch1.10。数据集源于在模拟矿井实验室进行SLAM过程中采集的RGB帧，分辨率为640×480，共1 802张图像。依次进行图像增强实验、线特征提取实验、线特征匹配实验。LT描述符的训练参数见表1。
表 1 LT描述符的训练参数Table 1. Training parameters of the LT descriptor参数 值 学习率 0.001 训练轮次 1 000 图像大小 640×480 线长度最小阈值 16 最大Token数 21 Token间距 8 描述符维度 256 注意力头数量 4 编码器特征维数 [32 64 128 256] 线段描述层数量 12 前馈层内部维度 1 024 签名网络层数 7 Transformer编码器层数 12 3.1 图像增强实验
分别采用改进同态滤波算法、EnlightenGAN算法[28]和本文算法进行实验。图像增强结果和对应的灰度分布如图6与图7所示。
可看出，3种算法对井下图像均有增强作用。使用EnlightenGAN算法后，图像对比度显著增强，但是亮度明显存在过度增强，且灰度级分布不均匀，大多分布在180～250。使用改进同态滤波算法增强后，图像的亮度与对比度有一定提高，像素灰度级分布较均匀。本文算法结合了同态滤波和CLAHE算法的优点，增强图像的亮度适中，对比度良好，灰度分布更加均匀。
图像增强的结果统计见表2。标准差是图像对比度的评价指标，标准差越大，代表图像对比度越大；均值是图像的亮度评价指标，均值越大，代表图像亮度越高；信息熵是反映图像信息量的评价指标，信息熵越大，代表图像信息量越大；峰值信噪比（Peak Signal-to-Noise Ratio，PSNR）表征图像的保真性，图像失真越小，该值越大。可看出EnlightenGAN算法增强图像的均值提高最大，从视觉主观来看是因为亮度被过度增强，且EnlightenGAN算法增强图像的PSNR仅为7.598，说明该算法增强图像失真较为严重，且信息熵也不及本文算法与改进同态滤波算法，说明图像部分信息丢失。经本文算法处理后图像的对比度、信息量、保真性均为最优，且亮度提升适中。从主观比较与客观数据综合来看，本文算法对井下低照度图像的增强效果最好。
表 2 图像增强结果统计Table 2. Statistics of image enhancement results算法 标准差 均值 信息熵 PSNR 改进同态滤波算法 63.583 123.820 7.712 10.014 EnlightenGAN算法 63.312 153.775 7.683 7.598 本文算法 64.108 129.174 7.798 10.178 3.2 线特征提取实验
采用增强前后的图像进行LSD线特征提取实验，分别在巷道、水房、避难硐室、工作面处选取视频帧共1 802帧。实验结果如图8所示，其中绿色表示增强前后共有的线特征，蓝色表示图像增强后新增的线特征。可看出图像增强后提取出的线特征更多。将部分断裂的线特征修复合并，可在一定程度上减少误匹配现象[29]。
LSD线段提取数量见表3，由于短线提取效率低，实验时滤掉了长度小于20像素的线段。本文算法提取的线段数量平均增长32.92%，说明在线段提取前进行图像增强处理这一关键步骤有效，可使LSD线特征质量更高，提取的线段更多。
表 3 LSD线段提取数量Table 3. LSD line segment extraction quantity采集点 帧数 LSD线平均数/条 增长率/% 原图 本文算法增强图像 巷道 480 94.8 125.6 32.49 水房 218 120.5 166.4 38.09 工作面 324 138.8 181.5 30.85 避难硐室 780 184.4 244.8 32.75 平均值 144.6 192.2 32.92 3.3 线特征匹配实验
3.3.1 图像选择及相关设置
采用视觉SLAM技术进行定位与建图时，最主要的运动包括平移与旋转，因此，挑选同时具有视角变化与旋转变化的井下低照度图像，分别选取变化程度小和变化程度大的2组代表性图像进行实验，以验证本文算法在不同程度视差干扰时的匹配性能。视觉SLAM技术经常应用于井下发生坍塌事故时的救援任务，场景局部密集的线缆、堆积的矿石及受困人员所处避难硐室的墙壁等均包含较多相似纹理，相似纹理越多，线特征的特征向量区分度越低，给匹配带来干扰。因此，选取相似纹理较少和相似纹理较多的2组代表性图像进行实验，以验证本文算法在相似纹理占比不同时的性能。如果在具有较大视差及较多相似纹理干扰的情况下，匹配性能仍然较高，说明算法的鲁棒性良好。
实验图像如图9所示，图像属性见表4。分别采用LBD、LBD_NNDR、LT、本文算法进行实验。LBD_NNDR对OpenCV视觉库中LSD线匹配的最近邻策略进行了优化：在目标与候选直线的配对中，首先计算出目标直线与每条候选直线特征向量的汉明距离，汉明距离越小，说明2个特征向量越相似；筛选出目标直线到候选直线的最小汉明距离${X_{\mathrm{m}}}$与次小汉明距离${X_{{\mathrm{sm}}}}$，如果${X_{\mathrm{m}}}$小于阈值$ {T_{\mathrm{d}}} $，则在此基础上根据最近邻距离比准则（Nearest Neighbor Distance Ratio，NNDR）进行二次筛选；若${X_{\mathrm{m}}}$与${X_{{\mathrm{sm}}}}$的比值小于阈值$ {T_{\mathrm{n}}} $，则认为该线段为最佳匹配，$ {T_{\mathrm{d}}} $与$ {T_{\mathrm{n}}} $选取工程上常用的30与0.8。LT描述符使用文献[24]推荐的权重。本文算法的匹配策略是计算出目标直线与每条候选直线L2范数归一化描述符向量之间的欧氏距离，利用欧氏距离表征2个线对的相似度，选择相似度最大的2个线对作为候选同名直线。根据工程经验，相似度阈值选取0.8能够同时兼顾精度与同名直线数量。
表 4 图像属性Table 4. Image attributes图像序号 采集位置 旋转与平移程度 相似纹理占比 1 巷道 较小 较小 2 避难硐室 较小 较大 3 巷道 较大 较小 4 避难硐室 较大 较大 3.3.2 实验结果
图像1匹配结果如图10所示。未匹配直线标为蓝色，匹配错误线对标为红色，正确匹配线对标为绿色，匹配正误由人工目视判别。由于短线匹配效率低，实验时滤掉了长度小于20像素的线特征。图像1旋转平移程度小、相似纹理较少，4种算法的性能均良好，精度分别为85.48%，89.65%，90.19%，92.06%，正确匹配数量分别为53对、52对、46对、58对。图像2、图像3匹配结果如图11、图12所示。
图像2相对于图像1存在相似纹理较多的干扰。图像3相对于图像1存在较大的视差干扰。从匹配结果可看出，本文算法性能最优，图像2的正确匹配数量高达100对，远超过LBD、LBD_NNDR、LT对应的28对、25对、52对。图像2的匹配精度为85.47%。图像3的正确匹配数量为39对，精度为92.86%。精度与正确匹配数量均维持在较高水平，说明本文算法在较大视差干扰下或者较大平移旋转视差下的鲁棒性良好。
图像4相对于图像1，同时具备了较大的视差及相似纹理较多的双重干扰，匹配结果如图13所示。LBD对图像4的匹配效果极差，几乎无法正确匹配，精度仅为15.00%，正确匹配数量仅为9对。LBD−NNDR能够剔除部分误匹配对，但是精度与正确匹配数量依然很低。LT的正确直线匹配数量与精度有一定提升。本文算法性能最优，线特征正确匹配数量达65对，精度达76.92%，精度与正确匹配数量维持在较高水平。这是因为本文算法在图像增强后，得到了更多可靠的LSD线特征，而且采用的LT描述符自适应地将变化较大的直线抽象成了固定维度，更加适合线段变化较大时的匹配，并且线的签名网络获取到邻域的几何属性，LT描述符构建出的特征向量包含更多位置及几何信息，能够在一定程度上克服相似纹理干扰。
3.3.3 统计分析
线特征匹配实验数据统计见表5。可看出对于任一类图像，本文算法得到的线特征数量、同名直线数量、正确匹配数量、精度均优于LBD，LBD_NNDR，LT。
表 5 线特征匹配实验数据统计Table 5. Statistics of experimental data of line feature matching图像
序号变化
程度相似纹
理占比算法 线特征数量/条 同名直
线数/对正确匹
配数/对匹配
精度/%左 右 1 较小 较小 LBD 125 105 62 53 85.48 LBD_NNDR 125 105 58 52 89.65 LT 125 105 51 46 90.19 本文算法 149 130 63 58 92.06 2 较小 较大 LBD 180 182 50 28 56.00 LBD_NNDR 180 182 33 25 75.76 LT 180 182 62 52 83.87 本文算法 327 335 123 100 85.47 3 较大 较小 LBD 100 181 20 12 60.00 LBD_NNDR 100 181 15 11 73.33 LT 100 181 37 34 91.89 本文算法 137 258 42 39 92.86 4 较大 较大 LBD 287 252 60 9 15.00 LBD_NNDR 287 252 18 8 44.44 LT 287 252 55 38 69.09 本文算法 350 385 65 50 76.92 4种算法的平均精度与平均正确匹配数分别如图14、图15所示。可看出本文算法平均精度为86.83%，较LBD，LBD_NNDR，LT分别提升32.71%，16.03%，3.07%。本文算法平均正确匹配数为61.75对，是LBD的2.422倍、LBD_NNDR的2.572倍、LT的1.453倍。统计分析结果说明本文算法性能优良，能够满足井下图像的稳健匹配需求。
4. 结论
1） 针对井下图像线特征匹配精度低、正确匹配量少、鲁棒性差的问题，提出了一种基于LSD和LT描述符的矿井图像线特征匹配算法，采用改进单参数同态滤波算法和CLAHE算法进行图像增强，在图像增强的基础上使用LT描述符对提取的LSD线进行描述和匹配。
2） 实验结果表明：本文算法的平均精度为86.83%，较LBD，LBD_NNDR，LT分别提升32.71%，16.03%，3.07%；平均正确匹配数为61.75对，是LBD的2.422倍、LBD_NNDR的2.572倍、LT的1.453倍；本文算法提取的LSD线数量更多，质量更好，在不同干扰下鲁棒性良好。
3） 下一步将继续优化本文算法的匹配策略，研究如何剔除更多的误匹配线对。
-
![]()
图 2 YSR25矿用防爆型雷达生命探测仪模拟现场应用
Figure 2. Simulated field application of YSR25 mine-used explosion-proof radar life detector
表 1 部分雷达生命探测设备
Table 1 Part of radar life detectors
型号 产地 性能指标 Lifelocator 3+ 美国 频率范围100～700 MHz，探测角度±60°，可探测6 m内的呼吸活动及10 m内的移动目标 SRI UWB−SAR 美国 中心频率450 MHz，距离分辨率1 m，带宽250 MHz，最大探测距离25 m Xaver−400 以色列 频率范围3～10 GHz，检测张角80°，穿墙范围4，8，20 m，分辨率≤1 m，最大探测距离100 m Po−400 GPR− Detector 俄罗斯 中心频率400 MHz，穿墙后可检测7 m以外的移动人体，最大探测距离20 m GPR−Detector 英国 中心频率400 MHz，探测深度5 m，穿墙后可测20 m内移动目标，最大探测深度5 m LEADERSCAN 法国 中心频率400 MHz，穿墙后最大探测距离10 m，最大探测距离30 m 
下载: 导出CSV
-
[1] 李斌,傅唯威,王勇. 基于周期抽样的超宽带生命探测雷达信号处理方法[J]. 仪器仪表学报,2010,31(9):1979-1985. LI Bin,FU Weiwei,WANG Yong. Signal processing for ultra-wideband life detection radar based on periodic sampling[J]. Chinese Journal of Scientific Instrument,2010,31(9):1979-1985.
[2] 费元春. 超宽带雷达与理论技术[M]. 北京:国防工业出版社,2010. FEI Yuanchun. Ultra wideband radar and theoretical technology[M]. Beijing:National Defense Industry Press,2010.
[3] 何昊. 雷达目标模拟器快速测频模块设计[D]. 成都:电子科技大学,2020. HE Hao. Design of fast frequency measurement module in radar target simulator[D]. Chengdu:University of Electronic Science and Technology of China,2020.
[4] SEKI T,SAITO K,TAKAHASHI M,et al. Proposition of UWB antenna wearing arm for WBAN[J]. ITE Technical Report,2009,33(10):1-4.
[5] HAFNER N,MOSTAFANEZHAD I,LUBECKE V M,et al. Non-contact cardiopulmonary sensing with a baby monitor[C]. The 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society,Lyon,2007:2300-2302.
[6] SHIKHSARMAST F M,LYU Tinging,LIANG Xiaolin,et al. Random-noise denoising and clutter elimination of human respiration movements based on an improved time window selection algorithm using wavelet transform[J]. Sensors,2018. DOI: 10.3390/s19010095.
[7] SOLDOVIERI F,CATAPANO I,CROCCO L,et al. A feasibility study for life signs monitoring via a continuous-wave radar[J]. International Journal of Antennas & Propagation,2014,12(6):1072-1075.
[8] 王健琪,董秀珍,王海滨,等. 基于毫米波的呼吸、心率非接触检测实验[J]. 第四军医大学学报,2001(2):180-182. WANG Jianqi,DONG Xiuzhen,WANG Haibin,et al. Experimental study on non-contact detection of breathing and heartbeat based on millimeter wave[J]. Journal of the Fourth Military Medical University,2001(2):180-182.
[9] 王保生. 超宽带信号穿透塌方体衰减特性及生命特征检测方法研究[D]. 哈尔滨:黑龙江科技大学,2017. WANG Baosheng. Research on attenuation characteristics and vital feature detection method of ultra wide band signal penetrating landslides [D]. Harbin:Heilongjiang University of Science and Technology,2017.
[10] 陈绍黔,王湘新,韩明华. 警用超宽带雷达式生命探测仪的研制[J]. 警察技术,2013(2):74-77. CHEN Shaoqian,WANG Xiangxin,HAN Minghua. Development of a police ultra wideband radar type life detector[J]. Police Technology,2013(2):74-77.
[11] 李秀贵,倪原. 一种UWB雷达脉冲信号发生器的设计[J]. 电子科技,2011,24(11):6-8. LI Xiugui,NI Yuan. Design of an ultra wide band radar pulse generator[J]. Electronic Science and Technology,2011,24(11):6-8.
[12] 胡巍. 基于多普勒雷达的非接触式生命体征检测技术研究[D]. 合肥:中国科学技术大学,2014. HU Wei. Research on non-contact vital signs detection technology based on doppler radar [D]. Hefei:University of Science and Technology of China,2014.
[13] 陈瑞鼎,鹿琪,单子涵,等. 基于卡尔曼滤波的超宽带穿墙雷达移动目标探测[J]. 地球物理学进展,2017,32(4):1758-1763. CHEN Ruiding,LU Qi,SHAN Zihan,et al. Moving target detection with the UWB through-wall radar based on Kalman filter[J]. Progress in Geophysics,2017,32(4):1758-1763.
[14] 周小龙. 无载波超宽带雷达人体动作识别理论方法研究[D]. 桂林:桂林电子科技大学,2019. ZHOU Xiaolong. Research on the theoretical method of human motion recognition with carrier-free UWB radar[D]. Guilin:Guilin University of Electronic Science and Technology,2019.
[15] VITEBSKIY S,CARIN L,RESSLER M A,et al. Ultra-wideband,short-pulse ground-penetrating radar:simulation and measurement[J]. IEEE Transaction on Geosci & Remote Sensing,1997,35(3):762-772.
[16] 饶俊峰,章薇,李孜,等. 雪崩三极管串联的纳秒脉冲发生器[J]. 强激光与粒子束,2018,30(9):77-83. RAO Junfeng,ZHANG Wei,LI Zi,et al. Nanosecond pulse generator with avalanche transistors in series[J]. High Power Laser and Particle Beams,2018,30(9):77-83.
[17] 何兴坤,钱鹏程. 基于雪崩三极管的脉冲发生器设计[J]. 电子测量技术,2019,42(21):65-69. HE Xingkun,QIAN Pengcheng. Pulse generator design based on avalanche transistor[J]. Electronic Measurement Technology,2019,42(21):65-69.
[18] RONGEN M,SCHAUFEL M. Design and evaluation of a versatile picosecond light pulser[J]. Journal of Instrumentation,2018,13(6). DOI: 10.1088/1748-0221/13/06/P06002.
[19] 张岩松,张亚东,梁步阁,等. 基于场效应管与阶跃恢复二极管的皮秒级脉冲源设计[J]. 强激光与粒子束,2017,29(4):81-85. ZHANG Yansong,ZHANG Yadong,LIANG Buge,et al. Design of picosecond pulse source based on MOSFET and step recovery diode[J]. High Power Laser and Particle Beams,2017,29(4):81-85.
[20] 王蕾,冯谞浩,孟庆鹤. 采用双阶跃恢复二极管的高幅度双极性窄脉冲探地雷达脉冲源的设计[J]. 电子技术与软件工程,2016(22):124-125,162. WANG Lei,FENG Xuhao,MENG Qinghe. Design of high amplitude bi-polar narrow pulse ground penetrating radar pulse source using double step recovery diode[J]. Electronic Technology & Software Engineering,2016(22):124-125,162.
[21] 王亚杰,何鹏军,荆晓鹏,等. 基于漂移阶跃恢复二极管开关的脉冲源仿真计算[J]. 强激光与粒子束,2018,30(9):92-97. WANG Yajie,HE Pengjun,JING Xiaopeng,et al. Simulation and calculation of pulsed power source based on drift step recovery diode switching[J]. High Power Laser and Particle Beams,2018,30(9):92-97.
[22] 周铀. 射频毫米波集成电路中功率放大技术研究[D]. 成都:电子科技大学,2018. ZHOU You. Research on power amplifying technology of RF/milimeter wave integrated circuits[D]. Chengdu:University of Electronic Science and Technology of China,2018.
[23] 孔德钰. 基于人工神经网络的低功耗射频收发系统研究[D]. 成都:电子科技大学,2019. KONG Deyu. Low power radio frequency transceiver system research based on artificial neural network[D]. Chengdu:University of Electronic Science and Technology of China,2019.
[24] 张彦婷,林文斌,唐晋生. 用于探地雷达的超宽带单极子天线设计[J]. 电子测量技术,2019,42(15):65-68. ZHANG Yanting,LIN Wenbin,TANG Jinsheng. Ultra-wideband monopole antenna design for ground penetrating radar[J]. Electronic Measurement Technology,2019,42(15):65-68.
[25] 郑艳. 基于CPW馈电的超宽带平面单极子天线研究[D]. 成都:西南交通大学,2019. ZHENG Yan. Research on ultra wideband planar monopole antenna based on CPW feeding[D]. Chengdu:Southwest Jiaotong University,2019.
[26] SHOME P P, KHAN T, LASKAR R H. CSRR-loaded UWB monopole antenna with electronically tunable triple band-notch characteristics for cognitive radio applications[J]. Microwave and Optical Technology Letters, 2020, 62（9）. DOI: 10.1002/mop.32394. SHOME P P,KHAN T,LASKAR R H. CSRR-loaded UWB monopole antenna with electronically tunable triple band-notch characteristics for cognitive radio applications[J]. Microwave and Optical Technology Letters,2020,62(9). DOI: 10.1002/mop.32394.
[27] 高冬冬,徐玉. 基于24 GHz多普勒雷达的非接触呼吸监测系统设计[J]. 自动化应用,2019(2):53-56. GAO Dongdong,XU Yu. Design of non-contact respiratory monitoring system based on 24 GHz doppler radar[J]. Automation Application,2019(2):53-56.
[28] 唐贵基,王晓龙. 参数优化变分模态分解方法在滚动轴承早期故障诊断中的应用[J]. 西安交通大学学报,2015,49(5):73-81. TANG Guiji,WANG Xiaolong. Parameter optimized variational mode decomposition method with application to incipient fault diagnosis of rolling bearing[J]. Journal of Xi'an Jiaotong University,2015,49(5):73-81.
[29] 刘长良,武英杰,甄成刚. 基于变分模态分解和模糊C均值聚类的滚动轴承故障诊断[J]. 中国电机工程学报,2015,35(13):3358-3365. LIU Changliang,WU Yingjie,ZHEN Chenggang. Rolling bearing fault diagnosis based on variational mode decomposition and fuzzy C means clustering[J]. Proceedings of the CSEE,2015,35(13):3358-3365.
[30] KUMAR A,ZHOU Yuqing,XIANG Jiawei. Optimization of VMD using kernel-based mutual information for the extraction of weak features to detect bearing defects[J]. Measurement,2021,168. DOI: 10.1016/j.measurement.2020.108402.
[31] 黄沁元,谢罗峰,殷国富,等. 基于变分模态分解和天牛须搜索的磁瓦内部缺陷声振检测[J]. 振动与冲击,2020,39(17):124-133. HUANG Qinyuan,XIE Luofeng,YIN Guofu,et al. Acoustic-vibration detection for internal defects of magnetic tile based on VMD and BAS[J]. Journal of Vibration and Shock,2020,39(17):124-133.
[32] WANG Ran,XU Lei,LIU Fengkai. Bearing fault diagnosis based on improved VMD and DCNN[J]. Journal of Vibroengineering,2020,22(5):1055-1068. DOI: 10.21595/jve.2020.21187
[33] 康守强,胡明武,王玉静,等. 基于特征迁移学习的变工况下滚动轴承故障诊断方法[J]. 中国电机工程学报,2019,39(3):764-772,955. KANG Shouqiang,HU Mingwu,WANG Yujing,et al. Fault diagnosis method of a rolling bearing under variable working conditions based on feature transfer learning[J]. Proceedings of the CSEE,2019,39(3):764-772,955.
[34] 马增强,李亚超,刘政,等. 基于变分模态分解和Teager能量算子的滚动轴承故障特征提取[J]. 振动与冲击,2016,35(13):134-139. MA Zengqiang,LI Yachao,LIU Zheng,et al. Rolling bearings' fault feature extraction based on variational mode decomposition and Teager energy operator[J]. Journal of Vibration and Shock,2016,35(13):134-139.
[35] 乔美英,闫书豪,兰建义,等. 基于VMD−TEO窗和DBiLSTM的早期轴承故障诊断[J]. 中山大学学报(自然科学版),2020,59(5):66-77. QIAO Meiying,YAN Shuhao,LAN Jianyi,et al. Early bearing fault diagnosis based on VMD-TEO window function and DBiLSTM[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni,2020,59(5):66-77.
[36] 吴小涛,江敏,孙洪军,等. 蝗虫优化相关向量机模型在径流预测中的应用[J]. 水电能源科学,2020,38(9):24-27,35. WU Xiaotao,JIANG Min,SUN Hongjun,et al. Application of relevant vector machine based grasshopper optimization algorithm in runoff prediction[J]. Water Resources and Power,2020,38(9):24-27,35.
[37] LIU Chongdang,ZHANG Linxuan,NIU Jiahe,et al. Intelligent prognostics of machining tools based on adaptive variational mode decomposition and deep learning method with attention mechanism[J]. Neurocomputing,2020,417:239-254. DOI: 10.1016/j.neucom.2020.06.116
[38] LI Jing , ZENG Zhaofa, SUN Jiguang, et al. Through-wall detection of human being's movement by UWB radar[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9（6）: 1-5. LI Jing ,ZENG Zhaofa,SUN Jiguang,et al. Through-wall detection of human being's movement by UWB radar[J]. IEEE Geoscience and Remote Sensing Letters,2012,9(6):1-5.
[39] 戴舜,朱方,徐艳云,等. 基于PCA与EMD的超宽带雷达生命信号检测算法[J]. 电子学报,2012,40(2):344-349. DAI Shun,ZHU Fang,XU Yanyun,et al. Vital signal detection method based on principal component analysis and empirical mode decomposition for ultra wideband radar[J]. Acta Electronica Sinica,2012,40(2):344-349.
[40] LI Jing, LIU Lanbo, ZENG Zhaofa, et al. Advanced signal processing for vital sign extraction with applications in UWB radar detection of trapped victims in complex environments[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 7: 783-791. LI Jing,LIU Lanbo,ZENG Zhaofa,et al. Advanced signal processing for vital sign extraction with applications in UWB radar detection of trapped victims in complex environments[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2014,7:783-791.
[41] 王亮. 基于超宽带雷达生命探测算法研究[D]. 杭州:浙江大学,2017. WANG Liang. Research of life detection for UWB radar[D]. Hangzhou:Zhejiang University,2017.
[42] 唐良勇,赵恒,张亚菊. 基于WA−EMD算法的脉冲式超宽带雷达多目标生命体征检测[J]. 南京理工大学学报,2017,41(2):198-206. TANG Liangyong,ZHAO Heng,ZHANG Yaju. Multiple-subject vital sign detection for impulse-radio ultra-wideband radar based on WA-EMD[J]. Journal of Nanjing University of Science and Technology,2017,41(2):198-206.
[43] 崔丽辉. 基于超宽带雷达的生命体征信息检测研究[D]. 济南:山东建筑大学,2017. CUI Lihui. Research on the detection of vital sign based on UWB radar[D]. Jinan:Shandong Jianzhu University,2017.
[44] 郭继坤,王保生,郝维来,等. 基于超宽带信号的矿井塌方体下生命特征的检测方法[J]. 黑龙江科技大学学报,2017,27(1):73-76,96. GUO Jikun,WANG Baosheng,HAO Weilai,et al. Research on detection method of life under mine based on ultra-wideband signal[J]. Journal of Heilongjiang University of Science and Technology,2017,27(1):73-76,96.
[45] 陈德肯. 基于雷达信号的非接触体征检测技术研究[D]. 北京:北京邮电大学,2019. CHEN Deken. Research on non-contact vital signs detection technology using radar signals[D]. Beijing:Beijing University of Posts and Telecommunications,2019.
[46] 郑学召,杨卓瑞,郭军,等. 灾后救援生命探测仪的现状和发展趋势[J]. 工矿自动化,2023,49(6):104-111. ZHENG Xuezhao,YANG Zhuorui,GUO Jun,et al. The current status and development trend of post-disaster rescue life detectors[J]. Journal of Mine Automation,2023,49(6):104-111.
[47] 文虎,周博,郑学召,等. UWB雷达在矿山钻孔救援中的应用研究[J]. 工矿自动化,2023,49(6):88-94. WEN Hu,ZHOU Bo,ZHENG Xuezhao,et al. Research on the application of UWB radar in mine drilling rescue[J]. Journal of Mine Automation,2023,49(6):88-94.
-
期刊类型引用(1)
1. 张科学，吕鑫淼，郑庆学，王晓玲，李小磊，刘昇，刘伟，李鑫磊，闫星辰，许雯，尹宇航. 基于eNSP的智慧矿山网络防火墙技术研究. 中国煤炭. 2024(08): 94-103 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 1602
- HTML全文浏览量: 115
- PDF下载量: 78
- 被引次数: 1
