ZHU Daiyin, ZHANG Ying, YU Xiang, et al. Imaging signal processing technology for miniature synthetic aperture radar[J]. Journal of Radars, 2019, 8(6): 793–803. doi: 10.12000/JR19094
Citation: YUAN Yubing, YE Shengbo, JI Yicai, et al. Fast refocusing algorithm based on three-dimensional wall compensation[J]. Journal of Radars, 2024, 13(4): 822–837. doi: 10.12000/JR24051

Fast Refocusing Algorithm Based on Three-dimensional Wall Compensation

DOI: 10.12000/JR24051
Funds:  The National Key R&D Program of China (2021YFC3002100); Science and Technology on Near-Surface Detection Laboratory (6142414220710)
More Information
  • Ultra-wideband through-wall radar, leveraging its ability to penetrate walls, can be used together with Multiple-Input Multiple-Output (MIMO) technology to image hidden targets behind walls. This approach provides rich information for detecting and locating people within buildings. This paper introduces a closed-loop interferometric calibration method based on a multitransmitter multireceiver ultrawideband wall-penetrating radar system in the Frequency Modulated Continuous Wave (FMCW) regime. This method aims to correct scattering issues caused by internal system errors. The presence of walls causes the target imaging position to deviate from the real position. To address this, this paper derives a three-Dimensional (3D) wall compensation algorithm jointing channels and pixel points. Then, a fast refocusing algorithm is proposed based on the geometric properties of the imaging area. The first step involves removing the influence of walls on delay time and determining the presence of the target. Subsequently, in view of the geometric properties of the region, a spherical coordinate grid division adapted to the region shape is selected. Localized refocusing is then performed in the subregion. This avoids the issue of electromagnetic wave attenuation, causing strong targets to mask weak ones in the imaging results. At the same time, the adoption of spherical coordinates for gridding and localized imaging greatly reduces the overall time consumption by the algorithm. Through simulation analysis and experimental verification, the proposed calibration method can effectively compensate for system errors. The fast refocusing algorithm can be used to realize multitarget 3D localization of the human body behind walls, with the localization accuracy of each dimension surpassing 10 cm and computational speeds improving by five times compared with those of existing algorithms. In terms of target detection probability, the proposed algorithm consistently identifies weak targets that other algorithms may overlook.

     

  • 感兴趣区域(Region of Interest, ROI)的高分辨成像与观测始终是微波遥感技术的关键和研究热点。合成孔径雷达(Synthetic Aperture Radar, SAR)凭借其全天时、全天候和穿透能力强等优势[1,2],能够获取米级至亚米级分辨率的ROI 2维散射图像,广泛应用于超视距探测、定位测量、目标检测与解译等领域。

    随着国内外轻小型飞机、无人飞行器的迅猛发展,搭载于这些小型平台上的微型合成孔径雷达(Miniature Synthetic Aperture Radar, MiniSAR)引起了广泛关注。MiniSAR系统通常采用调频连续波体制,其发射信号的峰值功率因而能够显著降低,并结合固态电子器件设计,以便于实现体积与重量的轻小型化。因此,MiniSAR作为一种机动灵活的成像系统[3],同样具备大带宽和高分辨能力,而且易于实现聚束、条带和圆迹等观测模式。在一定意义上,MiniSAR突破了传统大型机载雷达在时间与空间方面的严格限制,应对自然灾害突发、紧急救援、实时战场侦查等方面突显优势。

    为此,MiniSAR成像系统成为了一种提高ROI响应速度与侦查灵活性的可靠手段。尽管如此,MiniSAR成像信号处理依然面临若干技术难题:(1)低慢速小型平台及其复杂航迹难以实现较长相干积累时间的聚焦成像;(2)ROI非合作动目标成像常常伴随着多普勒频移与距离模糊,可引起不同程度的散焦效应;(3)轻小型平台约束下的GPS/IMU精度无法满足高分辨需求,需要与运动补偿技术相结合;(4)综合考虑算法复杂性、硬件资源和实时性等因素,与之相关联的并行处理架构的设计任务同样重要。文献[4]设计了一种W波段的无人机载SAR系统,采用了基于锁相技术的线性调频源设计,并对毫米波介质集成波导天线、3维集成、运动补偿方法等关键技术进行了研究,系统可实现7.5 cm的成像分辨率。考虑到MiniSAR复杂航迹,文献[5]通过结合快速分解后向投影(Fast Factorized Back-Projection, FFBP)算法与数据驱动运动补偿技术,并依据准极坐标网格的解析频谱推导了相位误差的空间不变性,从而保证了时域成像时自聚焦技术的嵌入与合并。在此基础上,为了进一步克服MiniSAR运动误差的影响,残留相位误差校正一般通过自聚焦技术实现[6],最显著的方法有相位梯度自聚焦(Phase Gradient Autofocus, PGA)和对比度优化自聚焦(Contrast Optimization Autofocus, COA)。文献[7]解析了散射体在距离多普勒域的固有稀疏特性,并构建了一种多任务学习稀疏贝叶斯框架实现稀疏散射系数与相位误差的估计。该方法优势在于通过适当引入估计的非确定性信息来改善聚焦精度,同时增强了在低信噪比和高欠采样比条件下的适用性。此外,MiniSAR图像的动目标散焦效应同样不可忽视,文献[8]以数据驱动运动补偿为切入点,具体涵盖了损失函数梯度距离对准、时间窗优化、COA自聚焦、时频信息重构和横向定标等一系列算法,有效解决了目标非均匀旋转引起的空变模糊。此类基于优化策略的数据驱动形式对于同时改善成像效率与聚焦深度开辟了新途径,也为观测系统处理效率的提高与实时性设计提供了可能。

    综上所述,针对机动平台MiniSAR高分辨成像信号处理技术,需要研究兼具高精度的聚焦算法和高效率的硬件架构设计。据此,本文提出了一套完整的成像信号处理方法及其对应的并行处理硬件架构。所提信号处理方案涵盖了极坐标格式算法(Polar Format Algorithm, PFA)粗聚焦、ROI动目标重聚焦和多子孔径PGA精聚焦算法;同时,该方案可映射至现场可编程门阵列(Field Programmable Gate Array, FPGA)技术,从而实现MiniSAR数据的高速并行处理与实时成像。最后,基于多组X波段MiniSAR实测数据处理结果,验证了所提成像信号处理技术的有效性和可靠性。

    MiniSAR成像算法是信号处理技术的核心。在高分辨条件下,算法需具备:(1)非常高的成像精度;(2)鉴于微小型无人机飞行的复杂工况,算法需具备自动补偿雷达平台非共面运动的能力;(3)回波数据量非常大,算法需简洁高效并易于硬件实现。传统的近似成像算法(如距离-多普勒算法)显然已不能满足成像精度需求。距离徙动算法(Range Migration Algorithm, RMA)和线频调变标算法(Chirp ScAling, CSA)作为常用频域算法,具有良好的计算效率,但这两种算法都要求雷达平台匀速直线飞行,同时RMA算法不适合嵌入运动补偿,CSA算法仍然存在近似误差,因此也不能很好地满足机动平台高分辨率MiniSAR成像需求。基于时域处理的BP算法能够在载机非理想飞行和非均匀采样条件下使用,但其主要缺点在于计算量太大,需要分布式计算处理,由此导致的成像实现硬件代价较大。相比之下,PFA算法同样在时域进行信号调整和补偿,易于校正平台非理想航迹,非常适合在大机动条件下的成像,而且算法同样简洁高效。此外,依据极坐标格式转换的两维尺度变换特性,基于尺度变换原理的PFA算法对成像效率与聚焦精度进行了新的解释[9],避免了直接插值操作,因而也避免了由于插值核长度有限而带来的相位误差。同时,算法映射至并行硬件处理器件,更可显示其高效性。因此,MiniSAR粗聚焦处理可依据PFA成像算法。

    雷达位置由瞬时方位与俯仰角为θφ决定,与之对应的参考方位角为固定值θref=π/2,参考俯仰角φref为孔径中心时刻角度。Rc为航迹到ROI中心点的最短距离。MiniSAR系统采用dechirp形式处理,其回波信号形式为

    s(τ,t)=Arect(τ2r/cTr)exp[j4πkc(τ2Rrefc)RΔ]exp(j4πcfcRΔ)exp(j4πkc2R2Δ)
    (1)

    式中,c表示光速,tτ分别为慢时间和快时间,Tr为脉冲宽度,k为调频斜率,fc为系统载频,r为雷达到目标的瞬时距离。对于条带模式参考距离Rref=Rc,且RΔ=rRc。需要说明的是,针对聚束和条带式统一化处理,均可按照聚束模式进行处理。此时,转为聚束模式处理时参考距离并非Rc,而是天线相位中心(Antenna Phase Center, APC)与ROI中心的距离Ra,且多旋翼无人机录取回波数据时RΔ=rRa

    PFA首先对每一距离向脉冲进行相位补偿

    sref(τ)=exp[j4πkc(τ+fck2Rcc)(RcRa)]exp[j4πkc2(RcRa)2]
    (2)

    并依据尺度变换原理对信号距离向数据乘以二次相位函数

    ϕscl(τ)=exp[jπk(1δr)(τ2Rac)2]
    (3)

    式中,δr(t)=cosφrefcos[θ(t)θref]cosφ(t)为尺度变换因子。所得数据经快速傅里叶变换(Fast Fourier Transform, FFT)后,需乘以匹配滤波函数

    H(fτ)=exp(jπδrkf2τ)exp{j2πfτ[fc(δr1)δrk+2(RaRc)c]}
    (4)

    式中,fτ为距离向采样频率。此时,所得数据经逆向傅里叶变换(Inverse Fast Fourier Transform, IFFT)后,再乘以二次相位函数:

    ϕins(τ)=exp{jπk(δrδ2r)[τ+fc(δr1)δrk2Rcc]2}
    (5)

    此外,方位向处理采用keystone变换,具体形式为

    t=fcfc+fτt
    (6)

    至此,所得2维重采样数据通过FFT操作即可成像,从而实现MiniSAR图像的粗聚焦。

    一般而言,ROI动目标成像常常伴随着多普勒频移与距离模糊,从而导致了MiniSAR图像中不同程度的散焦效应。为此,MiniSAR非合作动目标成像与重聚焦可借鉴逆合成孔径雷达(Inverse Synthetic Aperture Radar, ISAR)相关信号处理技术[10]。ROI动目标散焦数据可生成等效的ISAR回波,即可将MiniSAR非合作动目标重聚焦问题转化为ISAR相关技术。

    MiniSAR由于飞行工况复杂,对地面运动目标成像难度显著增加。借鉴以往机载SAR对地面运动目标成像处理的经验,本文提出对ROI动目标数据采用ISAR距离对准算法实现重聚焦成像的技术方案。该方案依据预设全局准则,迭代并优化性能指标,能够有效抑制突跳和漂移误差。其中,全局最小熵算法表现出较好的鲁棒性,但对准精度和收敛速度通常局限于一个距离单元和信号相关性[11,12]。因而,本文拟采用一种改进型全局最小熵距离对准算法,以提高对准精度和动目标的重聚焦质量。该方法依据平均距离像(Average Range Profile, ARP)熵值,同时ARP可由M个距离像列向量pi(r+Δri)的和包络计算获得。

    pave(r)=Mi=1|pi(r+Δri)|
    (7)

    式中,Δri为脉冲包络偏移量。从而ARP熵值为

    E(Δr)=rpave(r)Slnpave(r)Sdr,S=rpave(r)dr
    (8)

    距离对准需求解式(8)极小值实现,并获得估值

    [pi(r)pave(r)]r|r=Δri=0
    (9)

    ARP熵值可以通过搜索策略获取pi(r)pave(r)相关峰值及其对应的Δri值。然而,当E(Δr)终止迭代时,对准精度只能达到一个距离单元。为了进一步改善距离对准精度,算法考虑了时域亚距离单元移位问题。此处,将式(9)的卷积部分表示为R(Ω),并进行以下傅氏变换和离散化操作。

    R(Ω)FTR(fr)IFTR(Ω)=12πpi(fr)pave(fr)ejfrΩdfr1MTM/2k=M/2pi(2πMTk)pave(2πMTk)ej2πMkΩ
    (10)

    式中,pi(fr)pave(fr)分别对应于pi(r)pave(r)的傅氏变换。依据信号序列性质,式(10)可简化为

    ˜R(ξ)=M/21k=1Re[pi(2πMTk)pave(2πMTk)ej2πMkξ]
    (11)

    因而,式(11)可作为目标函数,并结合黄金分割法估计相关峰位置及其对应的亚距离单元偏移量。

    除此之外,复杂信号的回波相关性通常呈现一定的下降趋势,同样对包络对齐的精度造成影响。此时,式(8)同样可离散化并除去其中无关的常数因子后可得

    E(r)=Nr=1Mi=1pi(r)pave(r)
    (12)

    令内积Λ=pi(r),pave(r),依据施瓦茨不等式,当且仅当两者线性相关时获得最大值。因此,式(7)转化为pave(r)=Mi=1pi(r),其线性组合关系使得

    Λ=pi(r),Mi=1pi(r)=pi(r),p1(r)+pi(r),p2(r)+···+pi(r),pM(r)
    (13)

    显然,pi(r)之间的相关性影响着ARP熵值E(r)的最小值位置,且距离像之间的相关性与算法收敛速度成正相关。因而,该方法以分块形式提高动目标重聚焦精度和效率。

    尽管改进型ISAR距离对准算法的有效性作者已在参考文献[11,12]中进行了检验,然而由于飞行空域的限制,本MiniSAR系统到目前为止尚未对包含动目标的场景进行成像,因此上述方法在MiniSAR系统中的实施效果仍有待进一步检验。本文将在后续试验过程中设法进一步采集包含合作和非合作动目标的数据进行验证。

    高分辨情况下MiniSAR运动测量单元提供的位置信息精度可能无法满足相干性需求,且因传播介质不均匀而导致的回波延迟误差效应同样不可忽视。微型运动平台更容易受到气流扰动的影响而具有更大的轨迹偏离,严重影响聚焦质量。因此,本文针对MiniSAR图像精聚焦技术,采用多子孔径PGA-MD自聚焦,能够显著改善高分辨情况下的聚焦深度[13]

    假设子孔径脉冲数目和间隔分别为Nam0, m为方位坐标,此时对应的子孔径数目为P=[(MNa)/m0]+1,相邻子孔径重叠脉冲数目为Nam0。相位误差采用Q阶多项式模型Qq=0ϑqmq, ϑ为多项式系数,则第p个子孔径相位误差(Subaperture Phase Error, SPE)可描述为

    Φp(m)=rectp(m)Qq=0ϑqmqrectp(m)Qpq=0ϑqp(mpm0)q,0pP1
    (14)

    不同距离单元的多普勒频率通过PGA循环移位操作消除,然而常常伴随着未知线性相位的引入。此时,SPE梯度由式(15)组成

    ˜˙Φp(m)=˙Φp(m)+Δp+˙εp(m)
    (15)

    式中,˙Φp(m), Δp˙εp(m)分别代表实际的SPE梯度、未知常数项和杂噪信息。˜˙Φp(m)中未知线性分量Δp的校正与SPE梯度的相干拼接相关联。

    此时,第ij个子孔径对应于相位历史域的线性相位误差计算为

    Lj(m)Li(m)=(ΔjΔi)m=2πNaδi,jm
    (16)

    式中,方位相对偏移量δi,j通过利用MD算法对P(P1)/2个子图像进行互相关获得。此时,两者关系方程组为

    HΔ=δ
    (17)

    式中,H=[111111111][P(P1)/2]×(P1), Δ=[Δ1Δ2ΔP1], δ=2πNa[δ0,1δ0,P1δ1,2δP2,P1]

    方程组式(17)有唯一的最小均方解

    ˜ΔMD=[˜ΔMD1˜ΔMD2···˜ΔMDP1]T=(HTH)1HTδ
    (18)

    最后,SPE梯度的拼接操作可由式(19)实现。

    ˜˙ΦMDp(m)=˜˙Φp(m)˜ΔMDp,1pP1
    (19)

    此处需要说明的是,针对高分辨MiniSAR图像精聚焦处理,算法选择相邻子图像实现相位误差估计:(1)MiniSAR相邻子孔径信号相关性最优,对应相关峰值估计精度高;(2)子图像的减少能够显著降低H的维数和计算量,以提高相位补偿效率。

    综上所述,本文所提MiniSAR成像信号处理方案可以归结为一幅流程图,涵盖了PFA粗聚焦、ROI动目标重聚焦和PGA-MD精聚焦,如图1所示。

    图  1  所提MiniSAR成像信号处理流程图
    Figure  1.  Proposed flow diagram of MiniSAR imaging signal processing

    微型化、实时成像处理是系统设计和实现的关键。相较于大型机载SAR系统,MiniSAR系统的微型化设计,要在保证系统参数指标的前提下,尽量简化系统的设计结构。调频连续波SAR系统保留了传统脉冲SAR系统高分辨率的优点,又具有调频连续波雷达体积小、功率低、成本低的特点,使其适用于在无人机上进行低空成像。雷达接收机采用dechirp形式,结构相对简单,降低了AD器件的采样频率与设计成本,从而在保证系统分辨率的情况下,能够显著减小处理数据量,有利于后续高分辨率成像算法的实时处理。

    MiniSAR成像信号处理技术复杂度较高、数据吞吐量大,如何实时高效地处理数据具有实际工程价值。一般而言,GPU和FPGA作为常用的并行处理技术,均能可靠地应用于SAR相关信号处理[14]。然而,GPU较为适用于非相关性单指令多数据流的运算,面临大量不规则内存访问和同步扩展时,其运算容易导致低效率和高功耗。而FPGA适用于比特流操作,对系统实现细节具有较高级别的控制。依赖于丰富的查找表和寄存器等资源,并行方式FPGA芯片具备高吞吐量,低功耗,强抗干扰能力等优势。因此,FPGA并行架构具备灵活的处理资源配置和优化能力,可以较好适应MiniSAR成像信号处理技术需求。

    本文设计了一种基于FPGA架构的成像处理系统(暂未考虑动目标重聚焦技术),并如图2所示,具体涉及了3个核心子模块:数据传输子模块,通过以太网实现FPGA开发板与上位机的数据传输;算法处理子模块,对数据进行PFA成像、自聚焦与几何失真校正处理;DDR3 SDRAM读写子模块,实现2维数据的转置与连续地址的高速读写。

    图  2  基于FPGA设计架构的MiniSAR成像处理系统
    Figure  2.  MiniSAR imaging processing system based on FPGA design architecture

    (1) 数据传输模块:系统采用以太网实现主机与板卡之间的通信,并在网络层和传输层的协议控制上选择UDP/IP协议来实现。FPGA实现UDP/IP设计相对简单,资源利用较少。实测数据处理初始化阶段,上位机负责把MiniSAR参数和回波数据通过以太网传输到FPGA模块,数据写入DDR3 SDRAM。同时,也可将最终成像处理结果发送回上位机显示。

    (2) 成像处理算法模块:PFA粗聚焦模块主要由参数计算模块、距离向PCS模块、FFT模块、方位向插值模块组成;自聚焦模块主要由距离降分辨模块、方位相关复用模块、FFT模块、PGA时分复用模块、相位误差计算模块、相位拼接模块、相位补偿模块和子图像拼接模块组成;几何失真校正模块主要由校正网络铺设模块、存储地址计算模块和校正点地址写入模块组成。

    具体实现方式:PFA粗聚焦模块首先将MiniSAR回波数据写入到DDR3内存,同时传输其基本参数至参数计算模块。采用64位高精度浮点计算分别获得距离向、方位向与自聚焦处理所需的具体参数,对回波数据进行距离向PCS处理并以转置形式写入DDR。此外,参考距离补偿同样针对每一个脉冲,且可与距离向重采样一并处理,从而简化操作和节省反复读写DDR3的时间,提高了系统处理速度。最后转置读取每一个距离门数据通过Sinc插值核实现方位向重采样处理。图3为距离与方位实现的具体模块图,距离向处理过程中需要2次FFT(或IFFT)操作,虽然同时例化若干FFT的IP核能简化代码设计,但是极大地占用FPGA片上资源。而本系统中FFT模块和IFFT模块采用时分复用同一个FFT模块的方式,只需要例化一个FFT的IP核,可以保证运算速度的同时大大减少FPGA板内资源消耗。而方位向插值模块主要包含:待插值点判断模块和Sinc卷积模块。插值过程先根据待插点坐标寻找出相邻8个点的坐标位置,将坐标系数与输入的雷达回波数据对应相乘得到插值结果。针对处于样本边界外或不存在相邻8个点的待插值点,本设计中采用标识信号进行判断,使特殊点和正常点处理一体化,从而模块设计更简单、稳定。

    图  3  PFA 2维处理模块
    Figure  3.  2-D processing modules of PFA

    自聚焦模块以PFA极坐标格式转换后的数据开始,在DDR3中按方位向连续存储。第1步进行距离向降分辨率成像,并选取部分距离单元(设为2048)进行方位向到距离向的转置,同时设置FFT模块对距离向数据作2048点的FFT成像,以数据转置形式按方位向连续存放,以便于提高读写速率;第2步利用PGA-MD估计全孔径相位误差。本文设计的PGA模块结构如图4所示,并设置PGA模块的聚焦相关参数,启动PGA模块,直到每个子孔径数据全部处理完毕,并将估计的定点格式的相位梯度值传递给相位拼接模块。然后对聚焦后的相邻子孔径图像对进行方位互相关获得相对偏移量,并传递给相位拼接模块完成子孔径相位梯度拼接及积分得到全孔径相位误差。第3步设置分块PGA的相关参数,再次启动PGA模块,直到每个子图像处理完毕。同时利用方位互相关模块估计PGA聚焦前后的偏移量并完成子图像的拼接。

    图  4  相位梯度自聚焦模块结构图
    Figure  4.  Structure diagram of phase gradient autofocus module

    几何失真校正模块以逐距离向脉冲形式处理,计算第一条距离向校正点对应的PFA图像坐标,然后依据待处理数据在DDR3中的存储方式,将图像坐标转换为DDR3存储地址,再次逐距离向采样点数进行DDR3数据读取操作,完成该条脉冲的几何失真校正并顺序写入DDR3,重复上述操作直至所有脉冲全部处理完毕,即可获得最终成像结果。

    (3) DDR3读写转置模块:用以存储MiniSAR回波数据和算法处理结果,内置的读写控制模块可实现数据处理时所需的连续、转置读写操作,单通道内存为4 G,开启双通道后可实现8 G存储,满足实时处理的海量数据需求。本系统利用一片DDR3内存条,矩阵转置采用的是直接式分段存储方法。具体存储过程如图5所示。同时采用一种读写速率平衡的转置方式,其本质是在SDRAM的同一行中存放相同分量的列数据和行数据,在略微增加写操作换行次数的同时大大减少读操作换行次数,从而提高转置效率。

    图  5  直接式分段存储
    Figure  5.  Directly segmented storage

    所提FPGA成像系统已在Xilinx公司的Virtex7-XC7VX6907开发平台上进行了验证。当系统工作频率为200 MHz,处理8192(方位向)×4096(距离向)像素点的32位单精度浮点数据时,系统处理时间为3.58 s,此时FPGA处理时间低于对应的脉冲积累时间,表明其成像效率可满足实时性处理的需求。此外,表1给出了系统的FPGA资源利用情况。系统利用片内DSP48E1 资源来提高乘加模块,FFT模块的处理速度。使用最多的资源为LUT资源,占用率达到63%,表明了本设计较高的算法控制能力。其中Block RAM与DSP资源均占用30%左右,使用合理。浮点计算能力尚有冗余,完全可以处理较大数据量的数据或者增加其余算法处理模块,有利于系统后续的扩展使用。Block资源主要用于通过存储换取较高的算法处理效率,可以看出本系统的硬件框架通过较少的Block资源得到了较高的算法处理效率并且本系统还有较高的处理效率提升空间,也可采用多脉冲并行处理架构进一步提高处理速度,满足实时性设计需求。

    表  1  资源利用率
    Table  1.  Resource utilization
    资源已用资源可用资源利用率(%)
    Slice Registers35548086640041
    Slice LUTs27354943320063
    Block RAM/FIFO501147034
    DSP48E1s1204360033
    下载: 导出CSV 
    | 显示表格

    为了验证本文所提高分辨MiniSAR成像信号处理方案的有效性,依据图2给出的FPGA系统架构,对X波段MiniSAR实测数据进行了成像处理。MiniSAR系统的关键参数如表2所示,该系统采用调频连续波体制,搭载八旋翼无人机进行试验的场景如图6所示。

    表  2  主要的系统参数
    Table  2.  Main system parameters
    系统参数数值
    带宽(GHz)1.8
    载频(GHz)9.7
    飞行速度(m/s)5
    脉冲宽度(ms)4
    数据采样率(MHz)50
    脉冲重复频率(Hz)250
    下载: 导出CSV 
    | 显示表格
    图  6  MiniSAR试验系统
    Figure  6.  MiniSAR experimental system

    MiniSAR系统距离分辨率为0.1 m,并以聚束和条带模式完成了多组ROI回波数据采集。ROI#1的粗聚焦与精聚焦成像结果如图7所示,所提方法能够校正大机动微型平台引起的相位误差,并获得良好的聚焦效果。图8给出了ROI#2和ROI#3聚束成像与光学影像对比图,场景地物清晰,细节丰富,且对比度和信噪比较高,聚焦质量良好,反映了MiniSAR系统的高分辨成像能力。此外,图9图10分别给出了ROI#4和ROI#5的条带模式成像结果,实现了大测绘场景的成像与观测,从而进一步验证了所研制的高分辨率MiniSAR成像信号处理系统的有效性和可靠性。

    图  7  ROI#1成像处理与聚焦对比
    Figure  7.  Imaging processing and focusing comparison for ROI#1
    图  8  聚束式实测数据ROI#2和ROI#3处理结果
    Figure  8.  Measured data processing results for spotlight ROI#2 and ROI#3
    图  9  条带式实测数据ROI#4处理结果
    Figure  9.  Measured data processing result for stripmap ROI#4
    图  10  条带式实测数据ROI#5处理结果
    Figure  10.  Measured data processing result for stripmap ROI#5

    针对机动平台MiniSAR高分辨成像系统,本文提出了一套完整的信号处理方法及其对应的硬件设计架构。

    (1) 高精度MiniSAR成像信号处理方案涵盖了PFA粗聚焦、动目标重聚焦、PGA-MD精聚焦和几何失真校正算法;

    (2) 基于FPGA硬件设计的高速并行成像处理系统,契合了遥感监视系统的快速响应需求;

    (3) 多组X波段MiniSAR系统实测ROI处理结果验证了本文方法的有效性和可靠性。

  • [1]
    MAITI S and BHATTACHARYA A. Microwave detection of respiration rate of a living human hidden behind an inhomogeneous optically opaque medium[J]. IEEE Sensors Journal, 2021, 21(5): 6133–6144. doi: 10.1109/JSEN.2020.3043846.
    [2]
    YE Shengbo, ZHOU Bin, and FANG Guangyou. Design of a novel ultrawideband digital receiver for pulse ground-penetrating radar[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8(4): 656–660. doi: 10.1109/LGRS.2010.2098844.
    [3]
    FANG Guangyou. The research activities of Ultrawide-band (UWB) radar in China[C]. IEEE International Conference on Ultra-Wideband, Singapore, Singapore, 2007: 43–45. doi: 10.1109/ICUWB.2007.4380912.
    [4]
    QU Xiaodong, GAO Weicheng, MENG Haoyu, et al. Indoor human behavior recognition method based on wavelet scattering network and conditional random field model[J]. IEEE Transactions on Geoscience and Remote Sensing, 2023, 61: 5104815. doi: 10.1109/TGRS.2023.3276023.
    [5]
    金添, 宋勇平, 崔国龙, 等. 低频电磁波建筑物内部结构透视技术研究进展[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.
    [6]
    QIU Lei, JIN Tian, LU Biying, et al. An isophase-based life signal extraction in through-the-wall radar[J]. IEEE Geoscience and Remote Sensing Letters, 2017, 14(2): 193–197. doi: 10.1109/LGRS.2016.2633622.
    [7]
    金添, 宋勇平. 超宽带雷达建筑物结构稀疏成像[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.
    [8]
    YUAN Yubing, JI Yicai, YE Shengbo, et al. A clutter identification and removal method based on long delay lines and cross-correlation in through-wall detection[J]. Applied Sciences, 2024, 14(3): 1299. doi: 10.3390/app14031299.
    [9]
    MOHAMMED I, COLLINGS I B, and HANLY S V. Multiple target localization through-the-wall using non-coherent Bi-static radar[C]. 2019 13th International Conference on Signal Processing and Communication Systems, Gold Coast, Australia, 2019: 1–8. doi: 10.1109/ICSPCS47537.2019.9008415.
    [10]
    刘新, 阎焜, 杨光耀, 等. UWB-MIMO穿墙雷达三维成像与运动补偿算法研究[J]. 电子与信息学报, 2020, 42(9): 2253–2260. doi: 10.11999/JEIT190356.

    LIU Xin, YAN Kun, YANG Guangyao, et al. Study on 3D imaging and motion compensation algorithm for UWB-MIMO through-wall radar[J]. Journal of Electronics & Information Technology, 2020, 42(9): 2253–2260. doi: 10.11999/JEIT190356.
    [11]
    LIN Bo, LI Chao, JI Yicai, et al. Efficient scaling techniques for 2-D sparse MIMO array far-field imaging[J]. IEEE Sensors Journal, 2024, 24(8): 12604–12615. doi: 10.1109/JSEN.2024.3354246.
    [12]
    AFTANAS M, ROVNAKOVA J, DRUTAROVSKY M, et al. Efficient method of TOA estimation for through wall imaging by UWB radar[C]. 2008 IEEE International Conference on Ultra-Wideband, Hannover, Germany, 2008: 101–104. doi: 10.1109/ICUWB.2008.4653361.
    [13]
    AHMAD F, AMIN M G, and KASSAM S A. Synthetic aperture beamformer for imaging through a dielectric wall[J]. IEEE Transactions on Aerospace and Electronic Systems, 2005, 41(1): 271–283. doi: 10.1109/TAES.2005.1413761.
    [14]
    WANG Hanning, LU Biying, ZHOU Zhimin, et al. Through-the-wall imaging and correction based on the estimation of wall parameters[C]. 2011 IEEE CIE International Conference on Radar, Chengdu, China, 2011: 1327–1330. doi: 10.1109/CIE-Radar.2011.6159802.
    [15]
    CUI Guolong, KONG Lingjiang, and YANG Jianyu. A back-projection algorithm to stepped-frequency synthetic aperture through-the-wall radar imaging[C]. 2007 1st Asian and Pacific Conference on Synthetic Aperture Radar, Huangshan, China, 2007: 123–126. doi: 10.1109/APSAR.2007.4418570.
    [16]
    ROVNAKOVA J and KOCUR D. Compensation of wall effect for through wall tracking of moving targets[J]. Radioengineering, 2009, 18(2): 189–195.
    [17]
    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.
    [18]
    LIU Jiangang, KONG Lingjiang, YANG Xiaobo, et al. Refraction angle approximation algorithm for wall compensation in TWRI[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(7): 943–946. doi: 10.1109/LGRS.2016.2555291.
    [19]
    ZHAO Yang, LU Biying, and SUN Xin. Three-dimensional imaging for UWB though-the-wall radar[C]. 2013 IEEE Third International Conference on Information Science and Technology, Yangzhou, China, 2013: 1503–1506. doi: 10.1109/ICIST.2013.6747822.
    [20]
    AFTANAS I M. Through wall imaging with UWB radar system[D]. [Ph.D. dissertation], Technical University of Kosice, 2010.
    [21]
    GU Xiang and ZHANG Yunhua. Autofocus imaging simulation for through-wall radar by using FDTD with unknown wall characteristics[C]. 2010 Asia-Pacific Microwave Conference, Yokohama, Japan, 2010: 1657–1660.
  • Relative Articles

    [1]NIE Jiali, CUI Yuanhao, ZHANG Di, ZHANG Ronghui, MU Junsheng, JING Xiaojun. Vehicle Network Beamforming Method Based on Multimodal Feature Fusion[J]. Journal of Radars. doi: 10.12000/JR24242
    [2]LIU Shengheng, YAN He, LI Xingkang, XU Dazhuan, WANG Dongming, HUANG Yongming. HRT-Net: A Cell-free Cooperative Sensing Method Based on Attention-weighted Hierarchical Network[J]. Journal of Radars. doi: 10.12000/JR24227
    [3]ZHANG Yipeng, LU Dongdong, QIU Xiaolan, LI Fei. Few-shot Ship Classification of SAR Images via Scattering Point Topology and Dual-branch Convolutional Neural Network[J]. Journal of Radars, 2024, 13(2): 411-427. doi: 10.12000/JR23172
    [4]HAN Zhaoyun, CEN Xi, CUI Jiahe, LI Yachao, ZHANG Peng. Self-supervised Learning Method for SAR Interference Suppression Based on Abnormal Texture Perception[J]. Journal of Radars, 2023, 12(1): 154-172. doi: 10.12000/JR22168
    [5]TAN Kaiwen, ZHANG Limin, YAN Wenjun, XU Congan, LING Qing, LIU Hengyan. A Semi-supervised Emitter Identification Method for Imbalanced Category[J]. Journal of Radars, 2022, 11(4): 713-727. doi: 10.12000/JR22043
    [6]LI Wenna, ZHANG Shunsheng, WANG Wenqin. Multitarget-tracking Method for Airborne Radar Based on a Transformer Network[J]. Journal of Radars, 2022, 11(3): 469-478. doi: 10.12000/JR22009
    [7]ZHU Qingtao, YIN Junjun, ZENG Liang, YANG Jian. Polarimetric SAR Image Affine Registration Based on Neighborhood Consensus[J]. Journal of Radars, 2021, 10(1): 49-60. doi: 10.12000/JR20120
    [8]QIN Xianxiang, YU Wangsheng, WANG Peng, CHEN Tianping, ZOU Huanxin. Weakly Supervised Classification of PolSAR Images Based on Sample Refinement with Complex-Valued Convolutional Neural Network[J]. Journal of Radars, 2020, 9(3): 525-538. doi: 10.12000/JR20062
    [9]HU Tao, LI Weihua, QIN Xianxiang, WANG Peng, YU Wangsheng, LI Jun. Terrain Classification of Polarimetric Synthetic Aperture Radar Images Based on Deep Learning and Conditional Random Field Model[J]. Journal of Radars, 2019, 8(4): 471-478. doi: 10.12000/JR18065
    [10]HUA Wenqiang, WANG Shuang, GUO Yanhe, XIE Wen. Semi-supervised PolSAR Image Classification Based on the Neighborhood Minimum Spanning Tree[J]. Journal of Radars, 2019, 8(4): 458-470. doi: 10.12000/JR18104
    [11]ZHANG Xiaoling, ZHANG Tianwen, SHI Jun, WEI Shunjun. High-speed and High-accurate SAR Ship Detection Based on a Depthwise Separable Convolution Neural Network[J]. Journal of Radars, 2019, 8(6): 841-851. doi: 10.12000/JR19111
    [12]LIU Jiaqi, CHEN Bo, JIE Xi. Radar High-resolution Range Profile Target Recognition Based on Attention Mechanism and Bidirectional Gated Recurrent[J]. Journal of Radars, 2019, 8(5): 589-597. doi: 10.12000/JR19014
    [13]ZOU Huanxin, LI Meilin, MA Qian, SUN Jiachi, CAO Xu, QIN Xianxiang. An Unsupervised PolSAR Image Classification Algorithm Based on Tensor Product Graph Diffusion[J]. Journal of Radars, 2019, 8(4): 436-447. doi: 10.12000/JR19057
    [14]CHEN Huiyuan, LIU Zeyu, GUO Weiwei, ZHANG Zenghui, YU Wenxian. Fast Detection of Ship Targets for Large-scale Remote Sensing Image Based on a Cascade Convolutional Neural Network[J]. Journal of Radars, 2019, 8(3): 413-424. doi: 10.12000/JR19041
    [15]Zhao Juanping, Guo Weiwei, Liu Bin, Cui Shiyong, Zhang Zenghui, Yu Wenxian. Convolutional Neural Network-based SAR Image Classification with Noisy Labels[J]. Journal of Radars, 2017, 6(5): 514-523. doi: 10.12000/JR16140
    [16]Zou Huanxin, Luo Tiancheng, Zhang Yue, Zhou Shilin. Combined Conditional Random Fields Model for Supervised PolSAR Images Classification[J]. Journal of Radars, 2017, 6(5): 541-553. doi: 10.12000/JR16109
    [17]Tian Zhuangzhuang, Zhan Ronghui, Hu Jiemin, Zhang Jun. SAR ATR Based on Convolutional Neural Network[J]. Journal of Radars, 2016, 5(3): 320-325. doi: 10.12000/JR16037
    [18]Sun Xun, Huang Pingping, Tu Shangtan, Yang Xiangli. Polarimetric SAR Image Classification Using Multiple-feature Fusion and Ensemble Learning[J]. Journal of Radars, 2016, 5(6): 692-700. doi: 10.12000/JR15132
    [19]Xing Yanxiao, Zhang Yi, Li Ning, Wang Yu, Hu Guixiang. Polarimetric SAR Image Supervised Classification Method Integrating Eigenvalues[J]. Journal of Radars, 2016, 5(2): 217-227. doi: 10.12000/JR16019
    [20]Hua Wen-qiang, Wang Shuang, Hou Biao. Semi-supervised Learning for Classification of Polarimetric SAR Images Based on SVM-Wishart[J]. Journal of Radars, 2015, 4(1): 93-98. doi: 10.12000/JR14138
  • Cited by

    Periodical cited type(3)

    1. 李锵,王旭,关欣. 一种结合三重注意力机制的双路径网络胸片疾病分类方法. 电子与信息学报. 2023(04): 1412-1425 .
    2. 谢雯,滑文强,焦李成,王若男. 采用深度学习的极化SAR地物分类方法综述. 西安电子科技大学学报. 2023(03): 151-170 .
    3. 杨俊,刘灵辉. 加权全极化SAR图像分类下的功能区土地时空变化特征提取. 测绘通报. 2022(07): 138-142+153 .

    Other cited types(4)

  • Created with Highcharts 5.0.7Amount of accessChart context menuAbstract Views, HTML Views, PDF Downloads StatisticsAbstract ViewsHTML ViewsPDF Downloads2024-052024-062024-072024-082024-092024-102024-112024-122025-012025-022025-032025-040102030405060
    Created with Highcharts 5.0.7Chart context menuAccess Class DistributionFULLTEXT: 28.2 %FULLTEXT: 28.2 %META: 65.4 %META: 65.4 %PDF: 6.5 %PDF: 6.5 %FULLTEXTMETAPDF
    Created with Highcharts 5.0.7Chart context menuAccess Area Distribution其他: 2.3 %其他: 2.3 %其他: 0.2 %其他: 0.2 %Central District: 0.1 %Central District: 0.1 %China: 1.7 %China: 1.7 %Finland: 0.0 %Finland: 0.0 %France: 0.1 %France: 0.1 %Germany: 0.0 %Germany: 0.0 %India: 0.0 %India: 0.0 %Iran (ISLAMIC Republic Of): 0.1 %Iran (ISLAMIC Republic Of): 0.1 %Russian Federation: 0.0 %Russian Federation: 0.0 %Seattle: 0.0 %Seattle: 0.0 %Taiwan, China: 0.0 %Taiwan, China: 0.0 %United Kingdom: 0.0 %United Kingdom: 0.0 %United States: 0.3 %United States: 0.3 %[]: 0.4 %[]: 0.4 %三门峡: 0.0 %三门峡: 0.0 %上海: 0.6 %上海: 0.6 %东京: 0.0 %东京: 0.0 %东京都: 0.0 %东京都: 0.0 %东莞: 0.0 %东莞: 0.0 %中卫: 0.1 %中卫: 0.1 %临汾: 0.1 %临汾: 0.1 %丹东: 0.1 %丹东: 0.1 %乌鲁木齐: 0.0 %乌鲁木齐: 0.0 %九江: 0.0 %九江: 0.0 %佛山: 0.2 %佛山: 0.2 %六安: 0.0 %六安: 0.0 %兰州: 0.1 %兰州: 0.1 %兰辛: 0.0 %兰辛: 0.0 %北京: 9.4 %北京: 9.4 %南京: 1.5 %南京: 1.5 %南昌: 0.2 %南昌: 0.2 %卡拉奇: 0.0 %卡拉奇: 0.0 %厦门: 0.2 %厦门: 0.2 %台北: 0.1 %台北: 0.1 %台州: 0.0 %台州: 0.0 %合肥: 0.3 %合肥: 0.3 %吴忠: 0.2 %吴忠: 0.2 %哈尔滨: 0.1 %哈尔滨: 0.1 %哥伦布: 0.1 %哥伦布: 0.1 %大连: 0.0 %大连: 0.0 %天津: 0.4 %天津: 0.4 %太原: 0.1 %太原: 0.1 %威海: 0.0 %威海: 0.0 %宁波: 0.1 %宁波: 0.1 %安庆: 0.0 %安庆: 0.0 %安康: 0.4 %安康: 0.4 %宝鸡: 0.0 %宝鸡: 0.0 %宣城: 0.1 %宣城: 0.1 %宿迁: 0.0 %宿迁: 0.0 %巴音郭楞蒙古自治州: 0.0 %巴音郭楞蒙古自治州: 0.0 %常州: 0.1 %常州: 0.1 %平顶山: 0.0 %平顶山: 0.0 %广州: 0.3 %广州: 0.3 %库比蒂诺: 0.1 %库比蒂诺: 0.1 %张家口: 1.4 %张家口: 1.4 %张家界: 0.1 %张家界: 0.1 %徐州: 0.4 %徐州: 0.4 %德里: 0.2 %德里: 0.2 %惠州: 0.0 %惠州: 0.0 %成都: 1.9 %成都: 1.9 %扬州: 0.2 %扬州: 0.2 %无锡: 0.1 %无锡: 0.1 %昆明: 0.4 %昆明: 0.4 %晋中: 0.1 %晋中: 0.1 %晋城: 0.1 %晋城: 0.1 %朝阳: 0.0 %朝阳: 0.0 %杭州: 0.8 %杭州: 0.8 %查塔努加: 0.0 %查塔努加: 0.0 %柳州: 0.0 %柳州: 0.0 %格兰特县: 0.1 %格兰特县: 0.1 %桂林: 0.0 %桂林: 0.0 %榆林: 0.0 %榆林: 0.0 %武汉: 0.5 %武汉: 0.5 %池州: 0.0 %池州: 0.0 %沈阳: 0.3 %沈阳: 0.3 %沙田: 0.1 %沙田: 0.1 %泰州: 0.1 %泰州: 0.1 %泰米尔纳德: 0.1 %泰米尔纳德: 0.1 %洛阳: 0.0 %洛阳: 0.0 %济南: 0.3 %济南: 0.3 %海口: 0.1 %海口: 0.1 %淮南: 0.0 %淮南: 0.0 %深圳: 1.0 %深圳: 1.0 %温州: 0.1 %温州: 0.1 %渭南: 0.1 %渭南: 0.1 %湖州: 0.4 %湖州: 0.4 %湘潭: 0.0 %湘潭: 0.0 %漯河: 0.3 %漯河: 0.3 %烟台: 0.1 %烟台: 0.1 %牡丹江: 0.0 %牡丹江: 0.0 %石家庄: 0.1 %石家庄: 0.1 %福州: 0.1 %福州: 0.1 %秦皇岛: 0.1 %秦皇岛: 0.1 %纽约: 0.3 %纽约: 0.3 %绍兴: 0.2 %绍兴: 0.2 %维沙卡帕特南: 0.1 %维沙卡帕特南: 0.1 %绵阳: 0.1 %绵阳: 0.1 %美国伊利诺斯芝加哥: 0.0 %美国伊利诺斯芝加哥: 0.0 %自贡: 0.0 %自贡: 0.0 %舟山: 0.0 %舟山: 0.0 %芒廷维尤: 24.9 %芒廷维尤: 24.9 %芝加哥: 0.8 %芝加哥: 0.8 %苏州: 1.0 %苏州: 1.0 %莫斯科: 0.1 %莫斯科: 0.1 %菏泽: 0.1 %菏泽: 0.1 %葫芦岛: 0.2 %葫芦岛: 0.2 %衡阳: 0.0 %衡阳: 0.0 %衢州: 0.2 %衢州: 0.2 %西宁: 35.1 %西宁: 35.1 %西安: 2.4 %西安: 2.4 %诺伊达: 0.0 %诺伊达: 0.0 %贵阳: 0.1 %贵阳: 0.1 %赣州: 0.0 %赣州: 0.0 %达州: 0.1 %达州: 0.1 %运城: 1.1 %运城: 1.1 %通化: 0.1 %通化: 0.1 %邢台: 0.1 %邢台: 0.1 %邯郸: 0.1 %邯郸: 0.1 %郑州: 0.2 %郑州: 0.2 %重庆: 0.1 %重庆: 0.1 %长春: 0.0 %长春: 0.0 %长沙: 0.5 %长沙: 0.5 %阳泉: 0.0 %阳泉: 0.0 %雷恩: 0.1 %雷恩: 0.1 %青岛: 0.3 %青岛: 0.3 %首尔特别: 0.0 %首尔特别: 0.0 %香港: 0.1 %香港: 0.1 %香港特别行政区: 0.2 %香港特别行政区: 0.2 %马鞍山: 0.1 %马鞍山: 0.1 %齐齐哈尔: 0.1 %齐齐哈尔: 0.1 %其他其他Central DistrictChinaFinlandFranceGermanyIndiaIran (ISLAMIC Republic Of)Russian FederationSeattleTaiwan, ChinaUnited KingdomUnited States[]三门峡上海东京东京都东莞中卫临汾丹东乌鲁木齐九江佛山六安兰州兰辛北京南京南昌卡拉奇厦门台北台州合肥吴忠哈尔滨哥伦布大连天津太原威海宁波安庆安康宝鸡宣城宿迁巴音郭楞蒙古自治州常州平顶山广州库比蒂诺张家口张家界徐州德里惠州成都扬州无锡昆明晋中晋城朝阳杭州查塔努加柳州格兰特县桂林榆林武汉池州沈阳沙田泰州泰米尔纳德洛阳济南海口淮南深圳温州渭南湖州湘潭漯河烟台牡丹江石家庄福州秦皇岛纽约绍兴维沙卡帕特南绵阳美国伊利诺斯芝加哥自贡舟山芒廷维尤芝加哥苏州莫斯科菏泽葫芦岛衡阳衢州西宁西安诺伊达贵阳赣州达州运城通化邢台邯郸郑州重庆长春长沙阳泉雷恩青岛首尔特别香港香港特别行政区马鞍山齐齐哈尔

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(24)  / Tables(7)

    Article views(656) PDF downloads(155) Cited by(7)
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    • 表  1  资源利用率
      Table  1.  Resource utilization
      资源已用资源可用资源利用率(%)
      Slice Registers35548086640041
      Slice LUTs27354943320063
      Block RAM/FIFO501147034
      DSP48E1s1204360033
      下载: 导出CSV 
      | 显示表格
    • 表  2  主要的系统参数
      Table  2.  Main system parameters
      系统参数数值
      带宽(GHz)1.8
      载频(GHz)9.7
      飞行速度(m/s)5
      脉冲宽度(ms)4
      数据采样率(MHz)50
      脉冲重复频率(Hz)250
      下载: 导出CSV 
      | 显示表格