Processing math: 100%

低频电磁波建筑物内部结构透视技术研究进展

金添 宋勇平 崔国龙 郭世盛

金添, 宋勇平, 崔国龙, 等. 低频电磁波建筑物内部结构透视技术研究进展[J]. 雷达学报, 2021, 10(3): 342–359. doi: 10.12000/JR20119
引用本文: 金添, 宋勇平, 崔国龙, 等. 低频电磁波建筑物内部结构透视技术研究进展[J]. 雷达学报, 2021, 10(3): 342–359. doi: 10.12000/JR20119
JIN Tian, SONG Yongping, CUI Guolong, et al. Advances on penetrating imaging of building layout technique using low frequency radio waves[J]. Journal of Radars, 2021, 10(3): 342–359. doi: 10.12000/JR20119
Citation: JIN Tian, SONG Yongping, CUI Guolong, et al. Advances on penetrating imaging of building layout technique using low frequency radio waves[J]. Journal of Radars, 2021, 10(3): 342–359. doi: 10.12000/JR20119

低频电磁波建筑物内部结构透视技术研究进展

DOI: 10.12000/JR20119
基金项目: 国家自然科学基金(61971430)
详细信息
    作者简介:

    金 添(1980–),男,国防科学技术大学教授,博士生导师,研究方向为超宽带雷达成像、隐蔽目标检测识别等。全国百篇优秀博士论文获得者,入选教育部“新世纪优秀人才支持计划”,获国际无线电科学联盟青年科学家奖,评为中国电子学会优秀科技工作者。“信号处理与系统”国家精品课程和资源共享课主讲教师,信号处理系列课程国家级教学团队主要成员。已发表论文100余篇,获授权国家发明专利7项,出版专著3部、译著1部、教材1部。《雷达学报》、《信号处理》、《雷达科学与技术》等期刊编委

    宋勇平(1989–),男,空军预警学院雷达士官学校讲师,博士,研究方向为穿墙探测、MIMO雷达成像、微弱目标检测

    崔国龙(1982–),男,电子科技大学教授,青年长江学者,博士生导师,《雷达学报》编委。研究方向为最优化理论和算法、雷达目标检测理论、波形多样性以及城市环境目标探测等

    郭世盛(1991–),男,电子科技大学特聘副研究员,硕士生导师,研究方向为城市环境遮蔽目标探测、基于雷达的人体行为识别

    通讯作者:

    金添 tianjin@nudt.edu.cn

  • 责任主编:孔令讲 Corresponding Editor: KONG Lingjiang
  • 中图分类号: TN957

Advances on Penetrating Imaging of Building Layout Technique Using Low Frequency Radio Waves

Funds: The National Natural Science Foundation of China (61971430)
More Information
  • 摘要: 在进入陌生建筑物内部之前获取其内部结构信息,能够为反恐作战、灾害救援等多种应用提供服务,具有重要的现实意义和研究价值。低频电磁波能够穿透常见建筑物材料传播,进而安全、稳定、隐蔽地获取墙后目标信息。利用低频电磁波获取墙后信息因此成为建筑物内部结构透视领域的研究重点。为获知该领域的发展脉络,并预测未来可能的发展趋势,该文对21世纪初以来该领域国内外公开文献进行了归纳总结。相关文献的梳理结果表明,利用低频电磁波进行建筑物内部结构穿透探测的技术目前主要包括3类:基于反射波测量的穿墙雷达成像技术、基于透射波测量的射频层析成像技术、基于多径信号的墙体位置估计技术。这3类技术均已取得一定具有实际意义的研究成果。该文围绕这3类技术所涵盖主要内容的发展历程进行了梳理,主要包括穿墙雷达墙后静止目标成像原理、基于穿墙雷达的建筑物内部结构观测模式、基于穿墙雷达成像的建筑物内部结构重建技术、基于射频层析成像的建筑物内部结构反演技术、基于多径信号的墙体位置估计技术,并以此对该领域的发展趋势进行了探讨。总结近20年以来低频电磁波建筑物内部结构透视技术的发展历程,可以发现建筑物内部结构穿透探测平台已由传统的机载、车载平台转向微型机器人、无人机等新型平台,而对应的建筑物内部结构信息重建方法,则由传统的雷达成像技术,发展成包含图像增强、稀疏重构等在内的多种新型方法。这些结果表明,建筑物内部结构透视技术正朝着系统化、多样化、智能化的方向发展。

     

  • 超宽带雷达具有穿透能力强和距离向分辨率高等优点[1],它广泛应用于地质勘测、生命救援、穿墙目标跟踪等多种军用和民用领域[24]。传统的超宽带雷达是窄脉冲体制,它实现简单,但抗干扰能力差,在跟踪距离较远、干扰强等复杂的环境中,这种体制的雷达难以满足实际的要求[5,6]。新体制的超宽带雷达有线性调频体制、频率步进体制和伪随机编码体制,它们的工作方式类似,发射端发射持续时间较长,带宽足够,幅度相对较小的大时间带宽积信号;接收端接收到回波后先进行脉冲压缩或者是傅里叶变换处理,然后才得到雷达的脉冲响应函数。线性调频体制和频率步进体制雷达的信源设计复杂,系统线性度要求高,目前难以在地质勘测、生命救援、穿墙目标跟踪等应用中普及。伪随机编码信号可以由数字芯片直接合成,回波信号经过脉冲压缩而实现高信噪比和强抗干扰特性[7,8]。上世纪90年代中期,欧洲的火星探测计划促进了伪随机编码体制探地雷达的研制[9];德国伊尔梅瑙科技大学的学者研制了发射无载频基带信号的伪随机编码雷达;英国Utsi Electronics公司在2007年推出了中心频率为15 MHz的伪随机编码雷达Groundvue-6[10]。2012年开始,中科院电子所研制了超深编码探地雷达系统[11]和伪随机编码体制的穿墙雷达、生命探测雷达[12]

    近年来,超宽带雷达日益发展成熟,在工程应用中得到普及,但是单发单收雷达日渐显示出自身的不足之处。比如,在地质勘测方面,单发单收雷达观测到的目标信息通常认为是位于测线的正下方;同样地,穿墙目标跟踪和生命探测方面,单发单收雷达探测到的信息通常认为是位于雷达的正前方。虽然合成孔径雷达Synthetic Aperture Radar (SAR)方法能获得目标的距离向和方位向等多维信息[13],但是,平台运动的不确定性导致数据处理复杂,另一方面,平台的机械运动使观测信息难以实时融合处理,这使SAR方法难以推广到生命救援雷达和穿墙目标跟踪雷达等超宽带雷达的应用中去。2004年,美国学者Eran Fishler等人讨论总结了多发多收雷达的发展前景,在论文中指出,多发多收系统可以很大程度地提高雷达的探测性能[14]。多发多收超宽带雷达系统的多个收发通道协同工作,它能实时探测到目标的距离向和方位向信息[15]。美国的田纳西州大学研制了一款1发8收的多发多收雷达系统,该系统接收通道通过射频开关切换而实现8个通道分时输入[16,17]。但由于系统基于脉冲体制,脉冲峰值功率和射频开关耐压的矛盾而制约着雷达系统的性能。GSSI公司的多通道雷达系统SIR30可以实现8个通道同时工作,但是该系统结构复杂、成本高、功耗大,在8个通道同时工作时会造成数据量过大而难以实时处理。早在2002年,美国Time Domain公司就成功研制了一款多发多收穿墙雷达样机[15],这是第1款伪随机编码体制的多发多收超宽带雷达系统。对于伪随机编码体制的多发多收雷达系统,目前这方面的文献相对较少。

    本文提出一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计方法。为了降低雷达系统的成本和复杂性,同时为了尽可能地提高扫描率,采用射频开关切换的方法来实现多发多收雷达系统,该系统的最大通道能实现4发16收。选择码元长度为1023的二相调制m序列作为雷达的发射信号,信号中心频率为1.06 GHz,由于射频开关的耐压值有限,发射峰值功率为18 dBm。通过混合采样方法实现超宽带信号的采样接收,由200 MSPS的实时采样率通过64轮混合采样实现12.8 GSPS的混合采样率,并且调用FPGA内部的DSP硬核来进行实时脉冲压缩,以缓解后续数据处理运算量大的问题。实验结果表明,该雷达系统在距离向上有较远的探测距离,并且能实时探测到目标的方位向和距离向信息。

    伪随机编码多发多收雷达系统的工作过程如图1所示,图中为mn收雷达工作过程示意图。图1中Tm 表示第m号发射天线,Rn 表示第n号接收天线;Snm(t)表示第n号发射天线发射信号、第m号接收天线接收到的回波信号,回波信号已经过脉冲压缩处理。

    图  1  雷达工作过程示意图
    Figure  1.  Schematic diagram of radar operating

    伪随机编码雷达回波经过脉冲压缩处理后得到窄脉冲信号,如图1所示,因此伪随机编码雷达的探测原理与脉冲体制雷达相同。对于mn收的超宽带雷达系统,每个发射天线都有n组接收天线,所以mn收雷达系统总共有m×n组不同的接收数据。由于多个收发天线的空间相对位置不同,对这m×n组数据进行聚集处理便可以得到观测目标的距离向和方位向等信息。

    多发多收雷达可以通过单发单收雷达扩展得到,事实上,把多个独立的单发单收雷达组合起来工作就可以形成多发多收雷达系统,但这时会涉及到雷达之间的时间同步和数据融合问题;另一方面,多个独立的雷达会受到环境因素等的影响而难以保证各个雷达保持一致性。所以,通常会把多发多收雷达做成一个独立的雷达系统,这样不但解决了同步问题,也使数据融合处理显得简单可行。

    在超宽带雷达的应用中,观测目标通常是缓慢变化的,因此可以考虑以开关切换的方法实现多发多收超宽带雷达。对于mn收雷达,为了得到一帧完整的数据,发射开关要依次切换m轮,同时,发射开关每切换一轮,接收开关要依次切换n轮。即经过m×n轮切换才能得到一帧完整的m×n数据,雷达的成像速率比扫率慢m×n倍,因此要尽可能提高雷达系统的扫描率。在开关切换过程中,只要通道的切换速度足够快,就可以近似认为各个通道同时独立工作,并不丢失目标的实时信息。

    目前常用的信号切换开关有机械开关和射频开关两种,机械开关通过机械式接触实现不同通道的切换,它的耐压值高、隔离度好,但是切换速度慢,限制了多发多收雷达系统的切换速度。射频开关由半导体器件和微波电路组成,由电压控制来实现不同通道的切换,它的耐压值有限、隔离度适中,但是切换速度快。对于脉冲超宽带雷达,由于脉冲峰值过高而难以由射频开关切换实现多发多收雷达。伪随机编码体制雷达由于峰值功率小而可以很好地应用于基于射频开关切换的多发多收超宽带雷达。

    为了便于伪随机编码信号的发射,通常以二相调制的伪随机编码信号作为雷达的发射信号。为了得到二相调制信号,一种方法可以借鉴通信中的调制解调思想,采用混频方法把基本信号调制频带上,但是,由于超宽带信号带宽与中心频率相当,器件的宽带调制实现困难。另一种方法是基于FPGA和高速Digital to Analog Convertor (DAC)器件,通过数字方法直接合成伪随机编码雷达发射信号,研究和实验表明,这种方法能方便、可控地生成所需要的信号。

    由于FPGA芯片内部主频一般只能工作在500 MHz以下,为了产生高速的伪随机编码信号,结合FPGA内部的并串转换模块、双倍数据速率Double Data Rate (DDR)接口和多路数据接口的DAC芯片,通过FPGA内部低速时钟产生外部高速伪随机编码信号。图2是高速伪随机编码信号产生原理图,选取四路数据接口的高速DAC,其中fS是FPGA内部DAC采样数据存储器ROM1~ROM16的工作时钟,fD是DDR接口的工作时钟,fDAC是外部高速DAC的工作时钟。根据并串转换、DDR接口和四端口DAC芯片的时钟关系,有fDAC=8fDfD=2fS,即FPGA内部只需要低速时钟便能产生高速的伪随机编码信号。DAC产生的信号通过信号频带内的低通或带通滤波器后便能产生平滑的编码信号,编码信号经放大器放大后便是雷达天线的发射信号。

    图  2  高速伪随机编码信号产生原理图
    Figure  2.  Schematic diagram of high-speed pseudo random coded signal

    对于超宽带信号的采集,目前市场上能满足超宽带信号采样率要求的器件很少见,所以通常会采用混合采样对超宽带信号进行采集[18]。混合采样的思想是牺牲时间来换取更高的采样率,在多个雷达脉冲重复周期Pulse Repetition Period (PRP)内通过多次发送相同的信号来实现一道完整雷达回波的采集。图3是混合采样过程的时序图,采样过程中Analog to Digital Converter (ADC)工作时钟fADC依次比上一轮的混合采样时钟延迟 Δt。则混合采样的采样率为:

    fH=1Δt (1)

    采集一道完整回波的混合采样次数为:

    n=1fADCΔt (2)
    图  3  混合采样时序图
    Figure  3.  Timing diagram of hybrid sampling

    假设每一轮混合采样的采样点数为NP,则n轮混合采样得到的数据是一个n×NP矩阵,如式(3)所示。矩阵中每一轮的采样数据按行存储,而按列读出即是相邻时间间隔为 Δt的采样数据。因此,矩阵中的元素按依次从上到下、从左到右读取便可以得到一道顺序排列好的完整回波。

    S=[S1(0)S1(1)S1(NP1)S2(0)S2(1)S2(NP1)Sn(0)Sn(1)Sn(NP1)] (3)

    在伪随机编码体制的超宽带雷达中,原始回波信号的实时脉冲压缩是信号处理的首要和关键步骤。由于超宽带雷达的采样率高、数据量大,而微处理器处理速度和DSP芯片运算能力有限。本文基于时间域的互相关算法,通过调用FPGA中的DSP硬核进行并行计算,并结合流水线模式,实现快速实时脉冲压缩。伪随机编码信号的脉冲压缩公式为:

    R[l]=Nr1i=0s[i+l]r[i], l=0,1,2,···,NR1 (4)

    式中,r是参考码信号,它通常为雷达的发射信号,长度为Nr个采样点;s是雷达回波信号,长度为Ns个采样点。根据式(4),调用FPGA中的Nr个DSP硬核进行乘法运算,并通过流水线方式进行累加操作,构造图4所示的脉冲压缩并行运算框架。图中D表示1个时钟周期的延时,DSP硬核乘法器的一个输入端分别为各自固定的参考码,另一个输入端同时输入回波信号,压缩结果由最后一级输出端Rsr[i]输出。

    图  4  脉冲压缩并行运算框架
    Figure  4.  Structure of pulse compression computed parallelly

    由于FPGA中的DSP硬核数量有限,如果参考码长度Nr大于乘加器数量,脉冲压缩并行运算将不能实现。为了解决上面的问题,把长度为Nr的参考码分解成M段,每段长度为Nr0=Nr/M;分段条件满足Nr0小于FPGA中DSP硬核数量,同时通过对参考码补0或者适当选择M值使Nr0为整数。由分段方法,可以得到分段后的第j段参考码与原参考码对应关系为:

    rj[i]=r[jNr0+i] (5)

    式中,j=0, 1, ···,M–1,i=0, 1, ···,Nr0–1。分段后,式(4)可以转化成参考码分段并行运算脉冲压缩,如式(6):

    R[l]=M1j=0Rj[l] (6)

    式中,l=0, 1, ···,NR–1,Rj 表示第j段参考码的压缩结果,为:

    Rj[l]=Nr01i=0s[jNr0+l+i]rj[i] (7)

    即先把各段参考码采用图4的并行运算框架进行脉冲压缩,再把压缩的结果累加,容易证明式(4)和式(6)是等价的。

    多发多收雷达系统的原理图如图5所示,主要由4部分组成,从上到下依次为:时钟同步单元、发射机单元、接收机单元和开关切换单元,图中分别用不同的填充色标注。图6是多发多收雷达系统样机,该系统的主要参数如表1所示。

    图  5  多发多收雷达系统原理图
    Figure  5.  Schematic diagram of MIMO radar system
    图  6  多发多收雷达系统样机
    Figure  6.  MIMO radar system
    表  1  多发多收雷达系统主要参数
    Table  1.  Main parameter of MIMO radar
    参数 数值
    通道数量 mn收(最大4发16收)
    码元形式 m序列
    码元长度 1023
    单个码元宽度(ns) 0.9375
    发射信号幅值(Vpp) 5
    信号中心频率(GHz) 1.06
    –10 dB信号带宽(GHz) 0.2~2.1
    实时采样率(MSPS) 200
    混合采样率(GSPS) 12.8
    平均次数(times) 32
    接收机灵敏度(dBm) –97
    脉冲重复频率(KHz) 333.3
    采样时窗(ns) 300
    扫描率(track/s) 164
    成像速率(frame/s) 164mn
    总功耗(W) 20
    下载: 导出CSV 
    | 显示表格

    雷达系统需要考虑收发同步问题,特别在混合采样时,对同步精度要求是ns甚至ps级别的,所以雷达系统通常工作在同一时钟基准上。如图5第1部分所示,参考时钟是一个高稳的50 MHz晶振,50 MHz参考时钟经过频率综合器倍频后产生3.2 GHz时钟。3.2 GHz时钟扇出两路,一路作为DAC的数据采样时钟,另一路16分频后生成200 MHz的时钟。200 MHz时钟扇出两路,一路经过可编程延时线,作为ADC的工作时钟,另一路作为FPGA的输入时钟。FPGA的输入时钟由内部锁相环Phase Lock Loop (PLL)生成100 MHz的时钟,100 MHz作为FPGA系统的工作时钟。从上面分析得到,雷达系统的所有时钟都是以50 MHz高稳晶振作为基准,即雷达系统各部分始终是同步的。

    为了提高扫描率同时考虑信号的主旁瓣比,选择码元长度为1023的m序列作为发射信号,单个码元的宽度为0.9375 ns。结合DAC器件的特性同时考虑天线尺寸,选择二相调制信号作为雷达发射信号,信号的载波频率为1.06 GHz。DAC的工作频率为3.2 GHz,对每个码元采样3个有用点,理论上,此时DAC将会产生三角波,但受器件带宽和系统带宽的限制,实际波形会是平滑的双极性脉冲波形;经过低通滤波器后,波形将会变成类似于正弦波的二相调制信号。伪随机编码体制雷达通常发射小功率的信号,因此只需要用低噪声放大器将小信号放大,放大后的信号如图7所示。

    图  7  雷达发射信号
    Figure  7.  Radar transmitting signal

    接收天线输出的回波信号先通过低通滤波器,把不需要的高频成份滤去;由于回波信号通常很微弱,为了达到ADC芯片最低检测电平,低通滤波器后级联低噪声放大器把回波信号放大18 dB。ADC的工作时钟为200 MHz,由64轮混合采样得到采样率为12.8 GSPS,此时通过FPGA控制可编程延时线使ADC工作时钟每一轮依次比上一轮延迟78 ps。为了进一步提高回波信号质量,对每个通道的回波信号各自进行32次累加平均,即平均次数32次。

    发射信号的时间为960 ns,采样时窗为300 ns,实际采样时间为1260 ns。对于12.8 GSPS的采样率,参考码采样点数为12288点,回波采样点数为16128点,把参考码分成154段调用80个DSP硬核来完成实时脉冲压缩。图8是雷达系统闭环测试得到的脉冲压缩结果图,由于发射单元和接收单元都级联低噪放,为了防止在测试时接收机饱和,测试时在发射和接收机之间级联了40 dB的衰减器。脉冲压缩的旁瓣始终在–30 dB附近,它与m序列的特性一致,所以接收机单元满足雷达系统设计的需求。

    图  8  脉冲压缩结果
    Figure  8.  Results of pulse compression

    选择Mini-Circuits公司的一切四射频开关模块,射频开关的详细参数如表2所示。为了实现最大阵列为4发16收的的多发多收雷达系统,发射通道采用1个上述的一切四射频开关,接收通道采用5个一切四的开关由树型结构级联成一切十六射频开关,图9是收发开关的级联结构图。对于以上6个开关的控制,采用独立的开关控制板,开关控制板与雷达系统通信,由开关控制板实现任意收发组合的阵列形式。

    表  2  射频开关详细参数
    Table  2.  Main parameter of MIMO radar
    参数 数值
    带宽 0~3 GHz
    插损 1.5 dB@0.5~2 GHz
    1 dB压缩点 25 dBm@0.5~2 GHz
    隔离度 37 dB@0.5~2 GHz
    切换速度 25 ns/Turn on; 45 ns/Turn off
    下载: 导出CSV 
    | 显示表格
    图  9  收发开关级联结构图
    Figure  9.  Structure of radio frequency switches

    为了测试4发16收各个通道的一致性,对雷达系统4发16收开关的任意收发组合进行闭环测试,最终闭环测试得到的64组数据结果如图10所示。从图10可以知道,64组数据的脉冲压缩波形近似相同、幅度近似相等,除了在回波延时存在一定的差别外,其他特性基本相一致。而对于回波的延时差异,通过雷达系统校准可以在数据后处理时补偿回来。

    图  10  各个通道一致性测试结果
    Figure  10.  Results of consistency of all channels

    对系统的接收灵敏度进行测试,任意选择收发组合通道,发射信号和接收信号都经过低噪声放大器,中间级联可调衰减器后闭环。通过调节衰减器的大小,测试最小能检测到的信号电平。不同衰减值时的测试结果如图11所示,当衰减120 dB时虽然看到峰值点,事实上难以分辨。所以可以认为衰减115 dB时得到系统接收机的最小检测电平,由发射处输出电平为5 Vpp,通过计算得到接收机的灵敏度为–97 dBm。

    图  11  不同衰减值的闭环测试结果
    Figure  11.  Results of different attenuation

    为了验证雷达系统性能,选择单发单收模式下进行穿墙动目标测试,测试墙体为一般楼房的砖墙,厚度为25 cm,并同实验室研制的同频段脉冲雷达进行对比实验。同频段脉冲雷达的中心频率是1.2 GHz, –10 dB带宽为0.4~2.0 GHz,脉冲幅值为36 Vpp,采样率为13.6 GHz,线性平均次数为64次。对比实验测试结果如图12所示,从图中可以知道,伪随机编码雷达在15 m内跟踪结果非常清晰,而同频段的脉冲雷达的跟踪范围在10 m以内。脉冲雷达随着探测目标距离增大,信噪比迅速下降;而伪随机编码雷达可以通过脉冲压缩而进一步提高系统的信噪比。

    图  12  不同体制雷达的穿墙测试结果
    Figure  12.  Results of different radar systems though the wall

    为了方便验证多发多收雷达的性能,选用2发4收系统作测试,此时天线口径为1.8 m的线阵,其中2个发射天线在两侧,4个接收天线在中间;两侧发射天线和接收天线的间距是22.5 cm,中间相邻接收天线的间距为45 cm。首先在空气中对2发4收系统进行测试,主要是对空气中目标的距离向和方位向信息进行识别。图13是同时对空气中两个直径为7 cm的可比拟金属杯的测试场景,图14是两个目标中心点在不同间距时的测试结果图。

    图  13  2发4收对空气中两个目标测试场景
    Figure  13.  Experimental scene using 4-input 2-output radar to detect two targets in the air
    图  14  两个目标在不同间距时的测试结果图
    Figure  14.  Results of different distance of two targets

    图14是由BP算法处理得到的结果,从图中可以得到,两个目标的距离向固定在2 m处。当两个目标相距35 cm以上时,可以很容易地把它们区分开来,如图14(a)图14(b);当目标相距25 cm时,测试结果中两个目标开始重合,但仍然能区分开来,如图14(c);当目标相距20 cm时,两个目标的测试结果已经重合在一起,这时并不能把它们区分开来。

    为了进一步验证多发多收雷达的性能,仍然用以上空气中的2发4收阵列做穿墙实验,探测墙体后目标的距离向和方位向信息。图15是穿墙测试场景,此时直接对墙体后的人进行识别,图16是对墙体后目标的测试结果。测试实验增加了对背景数据(即无目标时)的记录,并将其与有目标时数据相减得到相关结果。

    图  15  2发4收穿墙测试场景
    Figure  15.  Experimental scene using 4-input 2-output radar to detect human through the wall
    图  16  墙体后目标的测试结果图
    Figure  16.  Results of human behind the wall

    测试墙体为一般楼房的砖墙,厚度为25 cm,图16是同样由BP算法处理得到的结果。图16(a)是20 cm×20 cm的金属板在墙后1 m处的测试结果,主要用来对系统进行标定;图16(b)是对墙后正前方人体测试的结果,目标在墙后4.2 m处;图16(c)是对墙后右前方的人体测试的结果,目标在墙后距离向4.2 m,方位向–1.7 m处;图16(d)是对墙后右前方的人体测试的结果,目标在墙后距离向3.8 m,方位向1.3 m处。从图中可以看出,2发4收雷达能对墙体后的目标进行方位向和距离向的区分。

    本文提出一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计方法。采用射频开关切换的方法来实现多发多收雷达系统,基于FPGA和DAC产生伪随机编码信号,由混合采样方法实现高速超宽带信号的接收采样,并直接在FPGA中进行实时脉冲压缩。实验结果表明,该雷达系统能实时探测到目标的方位向和距离向信息,在地质勘测、生命救援、穿墙目标跟踪等领域具有重要的应用价值。

  • 图  1  穿墙探测的一般信号传播模型

    Figure  1.  General signal propagation model for through-the-wall detection

    图  2  SIRE在多视角探测模式下的实测数据成像结果

    Figure  2.  Imaging results of SIRE’s measured data in multi-view detection mode

    图  3  单边双点观测模式及其实验结果[18]

    Figure  3.  Single-side two-location mode and its experiment result[18]

    图  4  SAPPHIRE雷达系统及其实验结果[19,20]

    Figure  4.  SAPPHIRE radar system and its experiment result[19,20]

    图  5  小型移动平台上的穿墙雷达系统及其实验结果[23]

    Figure  5.  Through-the-wall radar system mounted on a small mobile platform and its experiment result[23]

    图  6  全方向穿墙探测[28]

    Figure  6.  All-dricetions through-wall detection[28]

    图  7  多系统结合的穿墙探测示意图[31,32]

    Figure  7.  Schematic of multi-system combined through-wall detection[31,32]

    图  8  利用Hough变换提取墙体的直线特征[35]

    Figure  8.  Using Hough transformation to extract the straight-line features of the wall[35]

    图  9  利用M-N-K检测器提取建筑物结构[38]

    Figure  9.  Using M-N-K detector to extract building structure[38]

    图  10  基于多方位多尺度的建筑布局成像融合方法[40]

    Figure  10.  Multi-azimuth and multi-scale building layout imaging fusion method[40]

    图  11  原始模糊的建筑物结构图像与生成重建的建筑物结构图像[41]

    Figure  11.  Original blurred building structure images and generated building structure images[41]

    图  12  基于MST的建筑物内部结构估计[43]

    Figure  12.  Estimation of building internal structure based on MST[43]

    图  13  直线结构辅助下的墙体稀疏重构[46]

    Figure  13.  Sparse reconstruction of wall with the aid of linear structure[46]

    图  14  CF加权下的建筑物结构稀疏成像[47]

    Figure  14.  Sparse imaging of building structure under CF weighting[47]

    图  15  基于TV正则化的建筑物结构稀疏重构[48]

    Figure  15.  Sparse reconstruction of building structure based on TV regularization[48]

    图  16  基于WiFi信号的机器人平台穿墙探测及场景反演结果[57]

    Figure  16.  WiFi signal-based robot platform through-the-wall detection and scene inversion results[57]

    图  17  无人机载平台的WiFi探测实验与场景反演结果[59]

    Figure  17.  Experiments on Unmanned Aerial Vehicles (UAVs) and scene inversion results[59]

    图  18  基于RETINA的RTI算法[60]

    Figure  18.  RTI algorithm based on RETINA[60]

    图  19  基于椭圆协方差约束与叠加像素衰减非负约束的RTI结果对比[62]

    Figure  19.  Comparison of RTI results based on elliptic shape and superimposed pixel attenuation non-negative constraints[62]

    图  20  基于TV约束与PIC-ART的RTI结果对比[64]

    Figure  20.  Comparison of RTI results by TV constraint and by PIC-ART[64]

    图  21  基于传播时延的建筑结构层析成像[66]

    Figure  21.  Building structure tomography based on propagation time delay[66]

    图  22  基于内部辐射源的建筑物结构全息成像[67]

    Figure  22.  Holographic imaging of building structure based on internal radiation source[67]

    图  23  基于多径信号的目标定位与建筑布局联合估计[9]

    Figure  23.  Joint target location and building layout estimation based on multipath signals[9]

    图  24  利用目标运动轨迹反推墙体位置[68]

    Figure  24.  Use trajectory to infer wall position[68]

  • [1] NKWARI P K M, SINHA S, and FERREIRA H C. Through-the-wall radar imaging: A review[J]. IETE Technical Review, 2018, 35(6): 631–639. doi: 10.1080/02564602.2017.1364146
    [2] NGUYEN L, RESSLER M, and SICHINA J. Sensing through the wall imaging using the Army Research Lab ultra-wideband synchronous impulse reconstruction (UWB SIRE) radar[C]. The SPIE 6947, Radar Sensor Technology XII, Orlando, USA, 2008: 69470B. doi: 10.1117/12.776869.
    [3] AHMAD F and AMIN M G. Multi-viewpoint image fusion for urban sensing applications[C]. The SPIE 6968, Signal Processing, Sensor Fusion, and Target Recognition XVII, Orlando, USA, 2008: 69680L. doi: 10.1117/12.782302.
    [4] SENG C H, BOUZERDOUM A, PHUNG S L, et al. A two-stage image fusion method for enhanced through-the-wall radar target detection[C]. 2011 IEEE RadarCon (RADAR), Kansas City, USA, 2011: 643–647. doi: 10.1109/RADAR.2011.5960616.
    [5] SENG C H, BOUZERDOUM A, TIVIVE F H C, et al. Fuzzy logic-based image fusion for multi-view through-the-wall radar[C]. 2010 International Conference on Digital Image Computing: Techniques and Applications, Sydney, Australia, 2010: 423–428. doi: 10.1109/DICTA.2010.78.
    [6] SÉVIGNY P and DIFILIPPO D J. A multi-look fusion approach to through-wall radar imaging[C]. 2013 IEEE Radar Conference (RadarCon13), Ottawa, Canada, 2013: 1–6. doi: 10.1109/RADAR.2013.6586154.
    [7] SENG C H, BOUZERDOUM A, AMIN M G, et al. Probabilistic fuzzy image fusion approach for radar through wall sensing[J]. IEEE Transactions on Image Processing, 2013, 22(12): 4938–4951. doi: 10.1109/TIP.2013.2279953
    [8] PATWARI N and AGRAW P. Effects of correlated shadowing: Connectivity, localization, and RF tomography[C]. 2008 International Conference on Information Processing in Sensor Networks (ipsn 2008), St. Louis, USA, 2008: 82–93. doi: 10.1109/IPSN.2008.7.
    [9] LI Li and KROLIK J L. Simultaneous target and multipath positioning[J]. IEEE Journal of Selected Topics in Signal Processing, 2014, 8(1): 153–165. doi: 10.1109/JSTSP.2013.2289949
    [10] 张驰, 李悦丽, 周智敏. 基于独立分量分析法的稀疏阵列穿墙成像雷达直达波干扰抑制[J]. 雷达学报, 2014, 3(5): 524–532. doi: 10.3724/SP.J.1300.2014.14066

    ZHANG Chi, LI Yueli, and ZHOU Zhimin. Wall clutter mitigation in through-the-wall imaging radar with sparse array antenna based on independent component analysis[J]. Journal of Radars, 2014, 3(5): 524–532. doi: 10.3724/SP.J.1300.2014.14066
    [11] TAN Yunhua, WANG Libo, and LI Lianlin. A novel probability model for suppressing multipath ghosts in GPR and TWI Imaging: A numerical study[J]. Journal of Radars, 2015, 4(5): 509–517. doi: 10.12000/JR15056
    [12] JIN Tian, CHEN Bo, and ZHOU Zhimin. Image-domain estimation of wall parameters for autofocusing of through-the-wall SAR imagery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(3): 1836–1843. doi: 10.1109/TGRS.2012.2206395
    [13] FALCONER D G and UEBERSCHAER R M. Floor-plan radar[C]. The SPIE 4033, Radar Sensor Technology V, Orlando, USA, 2000, doi: 10.1117/12.391842.
    [14] DOGARU T and LE C. SAR images of rooms and buildings based on FDTD computer models[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(5): 1388–1401. doi: 10.1109/TGRS.2009.2013841
    [15] DOGARU T and LE C. Through-the-wall radar simulations for complex room imaging[R]. ARL-TR-5205, 2010.
    [16] LE C, DOGARU T, NGUYEN L, et al. Ultrawideband (UWB) radar imaging of building interior: Measurements and predictions[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(5): 1409–1420. doi: 10.1109/TGRS.2009.2016653
    [17] LE C, NGUYEN L, and DOGARU T. Radar imaging of a large building based on near-field Xpatch model[C]. 2010 IEEE Antennas and Propagation Society International Symposium, Toronto, Canada, 2010: 1–4. doi: 10.1109/APS.2010.5562119.
    [18] JIA Yong, ZHONG Xiaoling, LIU Jiangang, et al. Single-side two-location spotlight imaging for building based on MIMO through-wall-radar[J]. Sensors, 2016, 16(9): 1441. doi: 10.3390/s16091441
    [19] DE WIT J J M, VAN ROSSUM W L, and SMITS F M A. SAPPHIRE: A novel building mapping radar[C]. 2009 European Radar Conference (EuRAD), Rome, Italy, 2009: 625–628.
    [20] SMITS F M A, DE WIT J J M, VAN ROSSUM W L, et al. 3D mapping of buildings with SAPPHIRE[C]. The 6th Electro-Magnetic Remote Sensing - Defence Technology Centre (EMRS DTC) Technical Conference, Edinburgh, UK, 2009.
    [21] BROWNE K E, BURKHOLDER R J, and VOLAKIS J L. A novel low-profile portable radar system for high resolution through-wall radar imaging[C]. 2010 IEEE Radar Conference, Washington, USA, 2010: 333–338. doi: 10.1109/RADAR.2010.5494602.
    [22] SONG Yongping, ZHU Jiahua, HU Jun, et al. Grating lobes suppression for ultra-wideband MIMO radar imaging[C]. The 2017 3rd IEEE International Conference on Computer and Communications (ICCC), Chengdu, China, 2017: 957–961. doi: 10.1109/CompComm.2017.8322685.
    [23] BRAGA A J and GENTILE C. An ultra-wideband radar system for through-the-wall imaging using a mobile robot[C]. 2009 IEEE International Conference on Communications, Dresden, Germany, 2009: 1–6. doi: 10.1109/ICC.2009.5198740.
    [24] FALCONER D G, FICKLIN R W, and KONOLIGE K G. Detection, location, and identification of building occupants using a robot-mounted through-wall radar[C]. The SPIE 4037, Digitization of the Battlespace V and Battlefield Biomedical Technologies II, Orlando, USA, 2000, doi: 10.1117/12.395028.
    [25] SONG Yongping, ZHU Jiahua, HU Jun, et al. Enhanced imaging of building interior for portable MIMO through-the-wall radar[J]. Journal of Physics, 2018, 960: 012023.
    [26] SAKAMOTO T and SATO T. A method of estimating a room shape using a single antenna in a multipath environment[C]. The 4th European Conference on Antennas and Propagation, Barcelona, Spain, 2010: 1–5.
    [27] YEKTAKHAH B and SARABANDI K. All-directions through-the-wall radar imaging using a small number of moving transceivers[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(11): 6415–6428. doi: 10.1109/TGRS.2016.2585112
    [28] YEKTAKHAH B and SARABANDI K. All-directions through-the-wall imaging using a small number of moving omnidirectional bi-static FMCW transceivers[J]. IEEE Transactions on Geoscience and Remote Sensing, 2019, 57(5): 2618–2627. doi: 10.1109/TGRS.2018.2875695
    [29] YEKTAKHAH B and SARABANDI K. A method for detection of walls and large flat surfaces in through-the-wall SAR imaging[C]. 2019 IEEE International Symposium on Antennas and Propagation and USNC-URSI Radio Science Meeting, Atlanta, USA, 2019: 815–816. doi: 10.1109/APUSNCURSINRSM.2019.8889357.
    [30] ABEDI H and ZAKERI B. Through-the-multilayered wall imaging using passive synthetic aperture radar[J]. IEEE Transactions on Geoscience and Remote Sensing, 2019, 57(7): 4181–4191. doi: 10.1109/TGRS.2018.2890027
    [31] BARANOSKI E J. VisiBuilding: Sensing through walls[C]. The 4th IEEE Workshop on Sensor Array and Multichannel Processing, Waltham, USA, 2006: 1–22. doi: 10.1109/SAM.2006.1706221.
    [32] BARANOSKI E J. Through-wall imaging: Historical perspective and future directions[J]. Journal of the Franklin Institute, 2008, 345(6): 556–569. doi: 10.1016/j.jfranklin.2008.01.005
    [33] SUBOTIC N, KEYDEL E, BURNS J, et al. Parametric reconstruction of internal building structures via canonical scattering mechanisms[C]. 2008 IEEE International Conference on Acoustics, Speech and Signal Processing, Las Vegas, USA, 2008: 5189–5192. doi: 10.1109/ICASSP.2008.4518828.
    [34] BUONANNO A, D’URSO M, PRISCO G, et al. An improved radar system for inside-building awareness[C]. 2012 Tyrrhenian Workshop on advances in radar and remote sensing (TyWRRS), Naples, Italy, 2012: 85–89. doi: 10.1109/TyWRRS.2012.6381109.
    [35] AFTANAS I M. Through wall imaging with UWB Radar system[D]. [Ph. D. dissertation], Technical University of Košice, 2009: 68–75.
    [36] AFTANAS M and DRUTAROVSKÝ M. Imaging of the building contours with through the wall UWB radar system[J]. Radioengineering, 2009, 18(3): 258–264.
    [37] 贾勇, 孔令讲, 马静, 等. 穿墙雷达多视角建筑布局成像[J]. 电子与信息学报, 2013, 35(5): 1114–1119. doi: 10.3724/SP.J.1146.2012.01056

    JIA Yong, KONG Lingjiang, MA Jing, et al. Through-wall-radar multi-view imaging for architectural layout of building[J]. Journal of Electronics &Information Technology, 2013, 35(5): 1114–1119. doi: 10.3724/SP.J.1146.2012.01056
    [38] JIA Yong, CUI Guolong, KONG Lingjiang, et al. Multichannel and multiview imaging approach to building layout determination of through-wall radar[J]. IEEE Geoscience and Remote Sensing Letters, 2014, 11(5): 970–974. doi: 10.1109/LGRS.2013.2283778
    [39] 姚雪. 透视雷达建筑布局多视角成像算法研究[D]. [硕士论文], 电子科技大学, 2015: 41–44.

    YAO Xue. Perspective dectection research about the multi-view imaging of building layout[D]. [Master dissertation], University of Electronic Science and Technology of China, 2015: 41–44.
    [40] 唐茜, 晋良念. 多方位多尺度建筑物布局成像融合方法[J]. 雷达科学与技术, 2019, 17(1): 83–88, 93. doi: 10.3969/j.issn.1672-2337.2019.01.015

    TANG Qian and JIN Liangnian. Multi-azimuth and multi-scale building layout imaging fusion method[J]. Radar Science and Technology, 2019, 17(1): 83–88, 93. doi: 10.3969/j.issn.1672-2337.2019.01.015
    [41] JIA Yong, CHEN Shengyi, SONG Ruiyuan, et al. Building layout determination based on generative adversarial nets in through-wall radar imaging[C]. 2019 IEEE International Conference on Signal, Information and Data Processing, Chongqing, China, 2019: 1–6.
    [42] 陈波. 基于图理论的穿墙雷达建筑物布局重构技术研究[D]. [博士论文], 国防科学技术大学, 2014: 77–82.

    CHEN Bo. Study on building layout reconstruction techniques based on graph theory from through-the-wall radar[D]. [Ph. D. dissertation], National University of Defense Technology, 2014: 77–82.
    [43] CHEN Bo, JIN Tian, LU Biying, et al. Building interior layout reconstruction from through-the-wall radar image using MST-based method[J]. EURASIP Journal on Advances in Signal Processing, 2014, 2014(1): 1–9. doi: 10.1186/1687-6180-2014-1
    [44] 陈波, 金添, 陆必应, 等. 穿墙雷达中基于最小生成树的建筑物内部结构重构技术[J]. 电子学报, 2015, 43(9): 1682–1688. doi: 10.3969/j.issn.0372-2112.2015.09.002

    CHEN Bo, JIN Tian, LU Biying, et al. Building interior structure reconstruction from through-the-wall radar image using MST-based method[J]. Acta Electronica Sinica, 2015, 43(9): 1682–1688. doi: 10.3969/j.issn.0372-2112.2015.09.002
    [45] LAGUNAS E, AMIN M G, AHMAD F, et al. Compressive sensing for through wall radar imaging of stationary scenes using arbitrary data measurements[C]. The 2012 11th International Conference on Information Science, Signal Processing and their Applications (ISSPA), Montreal, Canada, 2012: 1347–1352. doi: 10.1109/ISSPA.2012.6310503.
    [46] LAGUNAS E, AMIN M G, AHMAD F, et al. Sparsity-based radar imaging of building structures[C]. The 20th European Signal Processing Conference (EUSIPCO), Bucharest, Romania, 2012: 864–868.
    [47] 金添, 宋勇平. 超宽带雷达建筑物结构稀疏成像[J]. 雷达学报, 2018, 7(3): 275–284. doi: 10.12000/JR18031

    JIN Tian and SONG Yongping. Sparse imaging of building layouts in ultra-wideband radar[J]. Journal of Radars, 2018, 7(3): 275–284. doi: 10.12000/JR18031
    [48] ZHAO Jifang, JIN Liangnian, and LIU Qinghua. Through-the-wall radar sparse imaging for building walls[J]. The Journal of Engineering, 2019, 2019(21): 7403–7405. doi: 10.1049/joe.2019.0541
    [49] 杨俊刚. 利用稀疏信息的正则化雷达成像理论与方法研究[D]. [博士论文], 国防科学技术大学, 2013: 1–2.

    YANG Jungang. Research on sparsity-driven regularization radar imaging theory and method[D]. [Ph. D. dissertation], National University of Defense Technology, 2013: 1–2.
    [50] SUNDARAM R, MARTIN R, and ANDERSON C. Regularization in radio tomographic imaging[C]. The SPIE 8753, Wireless Sensing, Localization, and Processing VIII, Baltimore, USA, 2013: 87530O. doi: 10.1117/12.2012167.
    [51] WILSON J and PATWARI N. Radio tomographic imaging with wireless networks[J]. IEEE Transactions on Mobile Computing, 2010, 9(5): 621–632. doi: 10.1109/TMC.2009.174
    [52] LIU Heng, WANG Zhenghuan, BU Xiangyuan, et al. Image reconstruction algorithms for radio tomographic imaging[C]. 2012 IEEE International Conference on Cyber Technology in Automation, Control, and Intelligent Systems (CYBER), Bangkok, Thailand, 2012: 48–53. doi: 10.1109/CYBER.2012.6392525.
    [53] MOSTOFI Y and SEN P. Compressed mapping of communication signal strength[C]. MILCOM 2008 - 2008 IEEE Military Communications Conference, San Diego, USA, 2008: 1–7. doi: 10.1109/MILCOM.2008.4753442.
    [54] MOSTOFI Y. Compressive cooperative sensing and mapping in mobile networks[J]. IEEE Transactions on Mobile Computing, 2011, 10(12): 1769–1784. doi: 10.1109/TMC.2011.31
    [55] MOSTOFI Y. Cooperative wireless-based obstacle/object mapping and see-through capabilities in robotic networks[J]. IEEE Transactions on Mobile Computing, 2013, 12(5): 817–829. doi: 10.1109/TMC.2012.32
    [56] GONZALEZ-RUIZ A, GHAFFARKHAH A, and MOSTOFI Y. An integrated framework for obstacle mapping with see-through capabilities using laser and wireless channel measurements[J]. IEEE Sensors Journal, 2014, 14(1): 25–38. doi: 10.1109/JSEN.2013.2278394
    [57] DEPATLA S, BUCKLAND L, and MOSTOFI Y. X-ray vision with only WiFi power measurements using Rytov wave models[J]. IEEE Transactions on Vehicular Technology, 2015, 64(4): 1376–1387. doi: 10.1109/TVT.2015.2397446
    [58] DEPATLA S, KARANAM C R, and MOSTOFI Y. Robotic through-wall imaging: Radio-frequency imaging possibilities with unmanned vehicles[J]. IEEE Antennas and Propagation Magazine, 2017, 59(5): 47–60. doi: 10.1109/MAP.2017.2731302
    [59] KARANAM C R and MOSTOFI Y. 3D through-wall imaging with unmanned aerial vehicles using WiFi[C]. The 2017 16th ACM/IEEE International Conference on Information Processing in Sensor Networks (IPSN), Pittsburgh, USA, 2017: 131–142.
    [60] HAMILTON B R, MA Xiaoli, BAXLEY R J, et al. Radio frequency tomography in mobile networks[C]. 2012 IEEE Statistical Signal Processing Workshop (SSP), Ann Arbor, USA, 2012: 508–511. doi: 10.1109/SSP.2012.6319745.
    [61] BECK B, BAXLEY R and MA Xiaoli. Regularization techniques for floor plan estimation in radio tomographic imaging[C]. 2013 IEEE Global Conference on Signal and Information Processing, Austin, USA, 2013: 177–180. doi: 10.1109/GlobalSIP.2013.6736844.
    [62] XU Shengxin, LIU Heng, GAO Fei, et al. Experimental verification: Enabling obstacle mapping based on radio tomographic imaging[C]. The 2018 IEEE 88th Vehicular Technology Conference (VTC-Fall), Chicago, USA, 2018: 1–5. doi: 10.1109/VTCFall.2018.8690867.
    [63] GUO Qichang, LI Yanlei, LIANG Xingdong, et al. Through-the-wall image reconstruction via reweighted total variation and prior information in radio tomographic imaging[J]. IEEE Access, 2020, 8: 40057–40066. doi: 10.1109/ACCESS.2020.2976769
    [64] ZHANG Yang, CUI Guolong, RAN Qingxin, et al. Knowledge-assisted building layout reconstruction for through-the-wall radar imaging[C]. 2020 IEEE Radar Conference, Florence, Italy, 2020: 1–5.
    [65] CHEN Jiahui, CUI Guolong, LIU Xinyu, et al. Strong scatterer reconstruction based on chirp-pulse microwave computed tomography[C]. 2019 IEEE International Conference on Signal, Information and Data Processing, Chongqing, China, 2019: 1–4.
    [66] GUO Qichang, LI Yanlei, LIANG Xingdong, et al. A novel CT-mode through-the-wall imaging method based on time delay estimation[J]. IEEE Geoscience and Remote Sensing Letters, 2020, in press. doi: 10.1109/LGRS.2020.3000423
    [67] HOLL P M and REINHARD F. Holography of Wi-Fi radiation[J]. Physical Review Letters, 2017, 118(18): 183901. doi: 10.1103/PhysRevLett.118.183901
    [68] LIU Jiangang, CUI Guolong, JIA Yong, et al. Sidewall detection using multipath in through-wall radar moving target tracking[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(6): 1372–1376. doi: 10.1109/LGRS.2015.2403133
    [69] 刘剑刚. 穿墙雷达隐蔽目标成像跟踪方法研究[D]. [博士论文], 电子科技大学, 2017: 95–106.

    LIU Jiangang. Imaging-tracking technology for hidden targets of a through-the-wall radar[D]. [Ph. D. dissertation], University of Electronic Science and Technology of China, 2017: 95–106.
    [70] GUO Shishen, CUI Guolong, KONG Lingjiang, et al. Multipath analysis and exploitation for MIMO through-the-wall imaging radar[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018, 11(10): 3721–3731. doi: 10.1109/JSTARS.2018.2865746
    [71] YANG Xiaqing, FAN Shihao, GUO Shisheng, et al. NLOS target localization behind an L-shaped corner with an L-band UWB radar[J]. IEEE Access, 2020, 8: 31270–31286. doi: 10.1109/ACCESS.2020.2973046
    [72] LI Songlin, CUI Guolong, GUO Shisheng, et al. On the electromagnetic diffraction propagation model and applications[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2020, 13: 884–895. doi: 10.1109/JSTARS.2020.2974529
  • 加载中
图(24)
计量
  • 文章访问数: 
  • HTML全文浏览量: 
  • PDF下载量: 
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-27
  • 修回日期:  2020-10-19
  • 网络出版日期:  2020-11-02
  • 刊出日期:  2021-06-28

目录

/

返回文章
返回