基于GNSS的无源双基地雷达地面动目标成像方法研究.pdf

DOI：10.12265/j.gnss.2023042基于 GNSS 的无源双基地雷达地面动目标成像方法研究唐涛，王鹏波，陈杰，周新凯，曾虹程(北京航空航天大学电子信息工程学院,北京 100191)摘要：基于全球卫星导航系统(GNSS)的无源双基地雷达是一种极具潜力的天基雷达系统，具有功耗低、隐蔽性好等优点.但由于信号功率密度低，需要通过长时间积累提升目标信噪比(SNR)，而目标运动会造成距离徙动和多普勒扩散，传统的脉冲多普勒雷达处理方法不再适用.针对此问题，本文提出一种基于拉东傅里叶变换的动目标成像方法.首先通过在传统拉东傅里叶变换中增加去多普勒调频率补偿因子，联合搜索并估计动目标的三维参数，然后通过参数补偿实现目标在距离多普勒域的聚焦成像.该方法既考虑了多普勒相位调制的影响，又在距离频域快速部署实现，可以达到比传统方法更好的聚焦效果.最后基于 GPS 卫星照射源的实测实验对所提方法的性能进行了验证.关键词：全球卫星导航系统(GNSS)；无源双基地雷达；动目标成像；拉东傅里叶变换；长时间相干积累中图分类号：P228.4；TP721文献标志码：A文章编号：1008-9268（2023）03-0024-090引言无源雷达是指使用机会辐射源作为雷达发射机，来进行目标探测、定位和跟踪的雷达系统，其具有设备简单、成本低、功耗低、抗干扰能力强、隐蔽性好等优点1.无源雷达适用于目标监控、国界监控等，在军事和民用领域都具有广阔的应用前景.近年来，全球导航卫星系统(GNSS)越来越多被应用于遥感领域2.在各类机会辐射源中，全球导航卫星信号具有独特的优势：一方面，GNSS 卫星提供全球覆盖，可以在地球表面的任何位置同时使用多颗GNSS卫星(例如GPS、GLONASS、Galileo 和北斗卫星导航系统(BDS)；另一方面，借助GNSS 授时服务可以轻松实现雷达时相同步.基于 GNSS 的遥感应用主要有两类：一类是GNSS-Reflectometry(GNSS-R)技术，主要应用包括反演海洋风速3、土壤水分4等；另一方面是基于GNSS 的双基地雷达成像5-6及动目标探测7-9.卫星导航信号的全球覆盖，为基于 GNSS 的无源双基地雷达系统的部署提供了绝佳的条件.但是，基于 GNSS 的雷达信号处理最大的难点在于卫星发射机距离地面很远，导致信号落地功率密度很低，动目标回波信噪比(SNR)太低，使得目标检测性能下降.有效的处理方案是通过长时间相干积累实现动目标成像，但目标的移动会引起信号积累时间内目标回波信号的距离偏移和多普勒频移现象，从而限制了处理时间.一系列基于检测前聚焦10-12的方法被应用于微弱目标的长时间相干积累数据处理中.文献 13 提出基于拉东傅里叶变换(RFT)的信号积累方法，利用目标速度和距离的联合搜索来实现目标回波的长时间积累，补偿了动目标的距离徙动和多普勒一次相位.文献 14 提出了一种广义的 RFT(GRFT)方法，可以实现运动目标回波的长时间相干积累，该方法在处理中涉及高维参数搜索，计算复杂度较高.文献 15提出了一种分段混合积累方法，它将较长的积累时间分成几个子段时间分别处理后再联合分析，可以避免参数遍历，但是信号的 SNR 增益会随着处理时间的延长降低.针对上述提到的回波信号 SNR 过低导致目标检收稿日期：2023-03-13资助项目：北京市自然科学基金(4222006)；国家自然科学基金(62101014)通信作者：王鹏波E-mail：第 48卷第3期全球定位系统Vol.48，No.32023年6月GNSS World of ChinaJune，2023测困难的问题，结合现有的信号积累方案，本文提出一种基于 RFT 的地面动目标成像方法.借鉴 RFT 的思路，在其中增加多普勒调频率参数搜索的步骤，通过方位去调频滤波器，增加信号相干积累的时长.就单目标而言，在参数空间域进行峰值提取，可以有效提取出动目标的多普勒频率和多普勒调频率参数，最终完成动目标精确聚焦和成像.该方案可以实现低SNR 下运动目标的检测前聚焦.通过动目标参数估计和补偿，该方法在距离多普勒域完成了动目标的聚焦和成像，为将基于 GNSS信号的无源双基地雷达系统应用到动目标探测领域提供了一定的理论和方法探索.本文的其余部分组织如下：第一节介绍了基于GNSS 的无源双基地雷达系统，建立了地面动目标成像的信号模型；第二节提出了基于 RFT 的动目标成像方法，通过联合搜索动目标的多普勒参数，完成参数估计与补偿，实现动目标聚焦成像；第三节详细介绍和分析了基于 GPSL5 信号的汽车目标实测实验，通过数据处理结果评估所提方案的有效性和性能；最后，对全文工作进行了总结和展望.1信号建模 1.1 基于 GNSS 的双基地雷达地面动目标探测系统在主动雷达、光电和红外摄像机等各类传感器系统中，无源双基地雷达被认为是一种有潜力的传感器系统，有助于在特定区域实现对动目标的持续监测.考虑到 GNSS 系统的特性，借助其作为照射源的无源双基地雷达系统可以直接安装在城市或者郊区的重点区域，实现隐蔽条件的环境监测.如图 1 所示，在设定的场景中，基于 GNSS 信号的地面动目标探测系统由三个主要部分组成，以提供对地面特定区域内动目标的连续监测.在该系统中，导航卫星作为机会辐射源，接收端可以采用地基或空基，以及信号处理模块完成最终的动目标成像任务.为了简洁，本文采用地基固定平台作为接收机，进行信号模型推导和后续的实测实验.1.2 动目标信号模型建立建立基于 GNSS 的无源双基地雷达回波信号模型是进行信号处理与分析的基础.以图 1 的几何构型为基础，以 GPS 卫星为照射源建立信号模型.需要说明的是，虽然每个时刻可以接收到多颗 GPS 卫星的信号，但是本文暂不考虑多卫星融合等，只挑选合适的单颗卫星数据进行动目标成像处理.单颗 GPS 卫星发射的连续信号可以建模为p(t)=n=g(tnTprt)=n=(D(tnTprt)C(tnTprt)sin(2fL5(tnTprt)+L5).(1)D()C()fL5L5Tprt式中：和分别为数据码和测距码；和分别为 GPSL5 信号的载波和初始相位；为 C/A 码的周期.GPS 信号的 C/A 码即导航信号测距码，重复周期为 1ms.目标反射后的导航卫星信号可以表示为S(t)=n=Wart(t)g(tnTprt(t).(2)War式中：为目标后向散射系数；为接收天线方向图.假设信号传播过程满足“停走”模型，即假设雷达平台在发射和接收信号时是静止的，则信号传播延迟可以表示为(t)=RS(t)+RR(t)c=R(t)c.(3)RS()RR()c式中：为卫星发射机到目标的距离；为目标到接收机的距离；为信号传播速度.最后，接收机经过混频、中频放大、单边带滤波、正交相干检波等，输出的复数信号为S(t)=n=WartR(t)cD(t)CA(tnTprtR(t)c)exp(j2R(t).(4)式中，为信号波长.1.3 双通道信号预处理基于 GNSS 的无源双基地雷达通过双天线接收，完成信号的获取.其中直射通道天线接收导航卫星的地面动目标地基接收机GNSS(雷达发射机)RSZYXRrefRR图 1 系统场景示意图第3期唐涛，等：基于 GNSS 的无源双基地雷达地面动目标成像方法研究25直达信号，反射通道天线用于接收运动目标的反射回波，完成运动目标的参数估计及聚焦成像.两个通道相互独立，并互相辅助来完成整个系统的运动目标探测任务.首先要通过双通道信号预处理，完成信号同步和卫星运动补偿等操作，再利用本文提出的动目标成像方法进行后续处理.双通道预处理的主要流程如图 2 所示，可分为直射通道信号处理和反射通道信号处理两步.直射通道信号处理反射通道信号处理中频解调距离维脉冲压缩数据码解调卫星信号捕获导航卫星跟踪与定位中频解调距离维脉冲压缩数据码解调二维数据生成卫星运动补偿后续的长时间相干积累和动目标聚焦成像图 2 双通道信号预处理流程首先对式(4)进行脉冲压缩和数据编码解调，可以得到：S(t)=n=WartR(t)cP(tnTprtR(t)c)exp(j2R(t).(5)P()WarRWar(t)War(nTprt)R(t)R(nTprt)式中，为测距码的自相关函数.可以认为和关于时间 t 是缓慢变化的，即，.通过变量替换和数据重排，可以将上述一维信号转化为二维形式，即：S(,)=War()P(R()c)exp(j2R().(6)式中：为慢时间；为快时间.可以知道回波能量的距离徙动和多普勒扩散主要来自两个方面，一是卫星相对于接收机的运动，二是目标的运动.因为目标的运动速度远小于卫星运动速度，所以要先补偿卫星运动.此时，信号变为S(,)=War()P(R()Rref()c)exp(j2R()Rref().(7)Rref式中，为卫星和接收器之间的距离.此时得到的信号，就可以通过后文提出的基于 RFT 的长时间积累方法进行处理成像.2基于 RFT 的动目标成像基于 GNSS 的地面动目标探测系统，由于信号落地功率低，要想实现动目标的有效探测，脉冲之间常常需要长时间的相干积累.此时目标回波往往具有明显的距离徙动和多普勒扩散.针对这一问题，本节提出了一种基于 RFT 的动目标成像方法，在传统 RFT的基础上，增加了对目标多普勒调频率的搜索，可以在长时间相干积累中达到更好的聚焦效果.在参数空间通过门限判决解算动目标参数，完成对动目标参数的补偿，最终实现在距离多普勒域的动目标的聚焦成像.下面从动目标参数估计和参数补偿及成像两方面对提出的动目标成像方法进行介绍.2.1 动目标参数估计2.1.1动目标回波信号解析R()Rref()将公式(7)中的在信号积累中心时刻进行泰勒展开，保留到二次项，可以得到：R()Rref()R0fd0.5fr2.(8)R0,fd,fr式中：分别为积累中心时刻的总距离.由于目标运动引入的信号多普勒频率和多普勒调频率.由于卫星导航信号的带宽较窄，距离分辨率一般只有数十米到数百米.式(8)中距离徙动的二次项往往小于距离分辨单元，可以忽略.同时可以忽略多普勒相位的常数项，并将其代回公式(7)中，则回波信号表达式变为S(,)=War()P(R0fdc)exp(j2fd0.5fr2).(9)为了快速部署，运用大量快速傅里叶变换(FFT)操作加快算法实现效率.将上述回波信号变到距离频域如下：S(,f)War()Pf(f)exp(j2f+fcc(fd0.5fr2).(10)Pf()P()fc式中：为的傅里叶变换；为信号载频.26全 球 定 位 系 统第48卷2.1.2基于 RFT 的参数搜索传统的利用 RFT 进行信号积累的方法，主要考虑了目标能量在距离维的线性距离徙动，并且考虑目标回波距离走动和多普勒频率之间的耦合关系.通过对多普勒频率参数的搜索实现目标能量的相干积累，在常规的信号积累中，因为积累时间较短，可以认为目标相对接收机做匀速运动，不会存在问题.但是在本文探讨的问题中，从 2.1.1 节中可以看到，由于导航信号的微弱和超长的积累时长，目标回波的多普勒徙动不可忽略.因此本文基于 RFT 的思路，又在其中引入了对目标多普勒调频率参数的搜索.为了提高算法效率，同时考虑距离徙动校正和多普勒相位的补偿，在回波的距离频域用一个补偿滤波器进行操作，该方法的处理表达式如下：Spara=Ta/2rTa/2IFFT(S(,f)Hc(f,fdi,fri)Hc(f,fdi,fri)=exp(j2(fc+f)c(fdita+0.5frita2).(11)式中：Ta为总积累时长；ta为方位时间；IFFT 为距离维的逆快速傅里叶变换.fdifriNRNfdNfrSpara上述操作的物理含义是：沿给定的轨迹曲线提取信号能量，然后对提取得到的信号进行多普勒相位解调制，最后对相位解调后的信号进行积分得到目标能量的积累结果.由于目标距离徙动轨迹和多普勒相位的耦合，对于每个多普勒频率搜索值和多普勒调频率搜索值，均需要更新能量提取路径.记三个参数的搜索次数分别为，和，将对应一个三维参数空间矩阵也就是上式(11)中的.2.1.3动目标参数解算当搜寻参数与目标的真实距离、多普勒频率、多普勒调频率等参数匹配时，提出的方法将输出一个极大值，该极大值就是目标回波能量的相干积累结果，其对应的搜寻参数即是目标真实运动参数的一个最优估计值.由三个搜索参数构成的参数空间的示意图如图 3 所示.多普勒调频率距离多普勒调频率fdSpara(R,fd,fr)maxRfr图 3 三维参数空间动目标参数解算为了方便，以单目标为例，这时可以直接选取参数空间域峰值对应的三维坐标作为目标的最优参数估计，即：Spara(R,fd,fr)max=Spara(Rest,fd_est,fr_est).(12)Rest,fd_est,fr_est当搜索参数取到的与动目标参数匹配时，参数空间达到峰值.2.2 动目标参数补偿及成像估算出动目标的多普勒频率和多普勒调频率后，可以通过频域补偿因子完成目标的运动补偿，从而完成距离徙动校正和多普勒扩散消除.首先通过在回波信号的距离频域乘以补偿因子 Hr完成距离徙动校正：Hr=exp(j2fd_esttafc),(13)然后通过距离维度 IFFT，将回波信号变换到距离时域，在方位向乘以去多普勒调频因子 Ha补偿掉二次相位调制；最后通过方位多普勒滤波器组处理，实现最终的距离多普勒域成像.Ha表达式如下所示：Ha=exp(jfr_estt2a).(14)整个基于 RFT 的动目标成像算法流程如图 4所示.第3期唐涛，等：基于 GNSS 的无源双基地雷达地面动目标成像方法研究27距离向脉冲压缩卫星运动补偿三维参数搜索:R,fd,fr变换到距离频域参数空间峰值坐标解算雷达回波数据距离徙动补偿多普勒徙动补偿方位多普勒滤波器组动目标距离多普勒域聚焦成像结果动目标参数补偿动目标参数估计图 4 基于 RFT 的动目标成像流程 2.3 算法性能分析fdfr本文提出的处理方法是在传统 RFT 方法的基础上，增加了对目标多普勒调频率的搜索.理论上，相干积累的增益只与相干积累的脉冲数(积累时长)有关，在动目标的多普勒频率和多普勒调频率均不可忽略时，可以通过最大相干积累时长来比较算法的处理增益.下面从计算复杂度和最大相干积累时长两方面对提出的方法和传统 RFT 方法分析对比.2.3.1计算复杂度对比NpNr,Nfd,Nfr用表示积累脉冲数，分别表示距离、多普勒频率和多普勒调频率搜索数.则传统RFT 方法的计算复杂度为IRFT=NrNfdNpIm+NrNfd(Np1)Ia+NfdIFFT(Nr).(15)ImIaIFFT(N)式中：和分别为复数乘法和复数加法的计算复杂度；表示 N 点 FFT 的计算复杂度.本文提出方法的计算复杂度为Iproposed=NrNfdNfrNpIm+NrNfdNfr(Np1)Ia+NfdNfrIFFT(Nr).(16)2.3.2最大相干积累时长对比Tint_RFTfr传统 RFT 方法的最大相干处理时长主要受到由多普勒调频率引起的多普勒徙动的约束，推导可得：Tint_RFT=1fr.(17)Tint_proposedfr本文提出方法的最大相干处理时长主要受到由多普勒调频率引起的距离徙动的约束，推导可得：Tint_proposed=8f0B1fr=8f0BTint_RFT.(18)式中：B 为信号带宽，f0为系统载频.通过以上分析可以看出，比起传统 RFT 方法，本文提出的方法由于增加了对多普勒调频率的搜索，增加了一定的计算复杂度，但是显著提升了回波信号的相干积累时长.对于非合作微弱目标的长时间积累检测问题，具有更好的处理效果.3实验验证 3.1 实测实验设置为了验证提出的基于 RFT 的动目标成像方法的有效性，本文开展了基于 GNSS 信号的运动目标长时间相干积累验证实验.该实验在中国北京昌平区北沙河中路上的某十字路口开展，为了使目标在接收天线波束内驻留时间足够长，本次实验的运动目标选择为行驶中的汽车，沿着公路自北向南朝接收机做减速运动，选择 GPSL5 波段信号为照射源，实验场景示意图如图 5 所示.行驶中的汽车NEWS接收机GPS PRN1图 5 基于 GNSS 的无源双基地雷达汽车探测示意图实验现场配置如图 6 所示，其中接收机位于目标南侧，迎着目标驶来方向.反射天线增益 10dB，指向目标方向，直射天线为全向 GPS 天线，指向天空以捕获尽可能多的导航卫星.通过直射通道完成 GPS 信号的捕获、跟踪和接收机定位工作，反射通道完成目标回波的长时间相干积累和动目标参数估计、补偿28全 球 定 位 系 统第48卷及聚焦成像.表 1 列出了实测实验中系统参数的详细信息.实验时刻的 GPS 星空图如图 7 所示，圆点是接收机所在位置，GPS 卫星 PRN1 被选为发射机.此时，雷达接收机和发射机位于待测目标的同一侧(正南侧)，构成一个相对标准的基于后向散射的双基地几何结构，有利于后续的信号处理与分析.天线和接收机待测汽车目标图 6 实测实验现场及设备照片表 1 实测实验参数设置参数取值参数取值信号源GPSL5卫星PRN1载频1176.45MHz等效PRF1000Hz系统采样率62MHz信号带宽10.23MHz天线增益10dBi积累时长4s选取卫星PRN136027090GPSGLONASSGalileoBeiDouQZSSIRBSSCutoff180(0)图 7 实验时刻 GPS 星空图(PRN1 被选择)3.2 数据处理及分析利用提出的动目标成像方法进行数据处理，可以对汽车目标在距离多普勒域中实现聚焦和成像.首先，对直达波信号进行处理，如图 8 所示，可以看到接收机直射通道成功捕获到了所选卫星.通过双通道处理，反射通道信号经过卫星运动补偿和距离压缩后的结果如图 9 所示，可以看到 4s 内汽车由远及近靠近接收机的运动轨迹.值得说明的是，为了更清楚地观测距离压缩的效果，此图进行了低通滤波，减小了观测的带宽范围.1789 1116 18232730卫星编号012345678捕获信号强度值图 8 直射通道卫星捕获列表012345004003002001000100双站距离/m方位时间/s200 300 400 500图 9 经过卫星运动补偿后的距离脉冲压缩结果直接对 4s 的回波数据用传统的方位多普勒滤波器组进行处理，可以得到目标的距离多普勒域粗成像结果，如图 10 所示.可以看到，在零多普勒区域存在很强的地面杂波和直射信号干扰，在目标检测前必须对其进行处理.由于地面杂波和直射信号的多普勒带宽仅由卫星的运动引起，而 GPS 卫星的轨道很高，因此地面杂波和直射信号的多普勒带宽很低，直接采用滤波器容易滤除.此外可以看到目标存在加速度，在做减速运动，目标回波能量在距离和多普勒维度都有明显的徙动，即回波能量在距离多普勒域存在二维扩散.这种情况下传统的脉冲多普勒雷达的处理算法和基本的 RFT 方法都无法得到动目标的高 SNR 聚焦结果.采用提出的基于 RFT 的方法，对动目标的距离、多普勒频率、多普勒调频率三维参数进行联合搜索.第3期唐涛，等：基于 GNSS 的无源双基地雷达地面动目标成像方法研究29Rfd=63.5 Hzfr=10.4 Hz/s图 11 显示了参数空间的 4s 相干积累结果，取出峰值处的二维切片(固定距离)，此时峰值对应的三维参数对应于动目标的距离和多普勒参数估计值.针对于本次实验，可以得到汽车对应的参数分别是双站距离=67.7m，多普勒频率，多普勒调频率.5004003002001000100双站距离/m200 300 400 500100多普勒频率/Hz80604020020406080100图 10 传统处理算法得到的动目标距离多普勒域成像结果251051208040040801202520151050多普勒调频率参数搜索/Hz多普勒频率参数搜索/Hz参数空间幅度/dB图 11 经过 4 s 相干积累处理的峰值处二维切面将估算出的动目标参数带入公式(13)(14)中处理，即可得到动目标距离多普勒域的聚焦成像结果，如图 12 所示.可以看到，因为经过了动目标参数补偿，目标信号回波在距离多普勒域聚焦结果良好，能量的二维扩散现象被消除.25105幅度/dB20015010050050100150200500300100100300500双站距离/m多普勒频率/Hz图 12 经过参数补偿后的动目标距离多普勒域成像结果本文对提出的动目标成像方法处理前后的距离多普勒结果进行了对比评估，结果如图 13 所示，可以看到聚焦后的 SNR 有了明显提升，提升了约 14dB，说明本文提出的动目标成像方案的有效性.3025201510505101 0008006004002000200双站距离/m幅度/dB400 600 800 1 000传统方法本文提出方法(a)距离维多普勒频率/Hz302520151050510幅度/dB20020406080100(b)多普勒维传统方法本文提出方法图 13 本文提出的方法处理前后动目标二维剖面对比4结束语基于 GNSS 信号的无源双基地雷达动目标检测及成像被认为是基于导航信号的遥感领域中有潜力的应用.本文针对长时间积累时动目标能量的距离徙动和多普勒扩散问题，提出一种基于 RFT 的动目标参数估计方法；然后通过参数补偿，实现动目标在距离多普勒域的聚焦成像.该方法关键步骤是联合搜索目标的距离、多普勒频率和多普勒调频率参数，所有操作在距离频域部署，估计出动目标多普勒参数，最终提升动目标成像效果.为了验证所提出的动目标成像方法的可行性，进行了以 GPSL5 信号作为机会照射源的汽车目标实测试验，可以看到经过参数估计和补偿后，成像结果良好.特别是通过对 RFT 方法的补充，积累后的动目标 SNR 显著增加，验证了所提方案的有效性.30全 球 定 位 系 统第48卷然而，为了实现复杂背景条件下运动目标检测成像及后续处理，仍有很多技术难题待突破，如信号长时间积累时效率与效果的兼顾、运动目标的多帧检测跟踪等，这些将是后续研究的重点.目前本文所做的基于 GNSS 信号的实测实验是有限的，在未来的研究工作中，将开展更多的实验，探索基于 GNSS 的无源双站雷达系统的应用.参考文献KUSCHEL H,CRISTALLINI D,OLSEN K E.Tutorial:passive radar tutorialJ.IEEE aerospace and electronicsystems magazine,2019,34(2):2-19.DOI:10.1109/MAES.2018.1601461杨东凯,李晓辉,王峰.GNSS反射信号海洋遥感应用现状分析J.无线电工程,2019,49(10):843-848.2刘原华,何孟然,牛新亮.GNSS-R海面风速反演技术研究J.全球定位系统,2020,45(2):55-59.3严颂华,龚健雅,张训械,等.GNSS-R测量地表土壤湿度的地基实验J.地球物理学报,2011,54(11):2735-2744.4范学祯,刘飞峰,张天,等.基于北斗导航卫星的BiSAR大场景成像J.信号处理,2017,33(4):546-550.5田卫明,曾涛,胡程.基于导航信号的BiSAR成像技术J.雷达学报,2013,2(1):39-45.6CHERNIAKOVM,SAINIR,ZUOR,etal.Space-surfacebistatic synthetic aperture radar with global navigationsatellite system transmitter of opportunity-experimentalresultsJ.IET radar,sonar&navigation,2007,6(1):447-458.DOI:10.1049/iet-rsn:200601727李中余,黄川,武俊杰,等.基于GNSS的无源雷达海面目标检测技术综述J.雷达科学与技术,2020,18(4):404-416.8王博,朱云龙,高超群,等.GNSS-R双基SAR运动目标成像技术探讨J.导航定位学报,2019,7(1):59-64.9MAH,ANTONIOUM,PASTINAD,etal.MaritimemovingtargetindicationusingpassiveGNSS-basedbistaticradarJ.10IEEEtransactionsonaerospaceandelectronicsystems,2018,54(1):115-130.PIERALICEF,SANTIF,PASTINAD,etal.GNSS-basedpassiveradarformaritimesurveillance:longintegrationtimeMTI techniqueC/IEEE Radar Conference(RadarConf),2017.11何振宇,杨扬,陈武,等.基于GPS-R的双基地SAR移动目标成像方法研究J.遥感学报,2022,26(12):2555-2567.12XUJ,YUJ,PENGYN,XIAXG,etal.Radon-Fouriertransform for radar target detection,I:generalized dopplerfilterbankJ.IEEEtransactionsonaerospaceandelectronicsystems,2011,47(2):1186-1202.DOI:10.1109/TAES.2011.575125113XUJ,XIAXG,PENGSB,etal.Radarmaneuveringtargetmotion estimation based on generalized Radon-FouriertransformJ.IEEEtransactionsonsignalprocessing,2012,60(12):6190-6201.DOI:10.1109/TSP.2012.221713714XU J,ZHOU X,QIAN L C,et al.Hybrid integration forhighly maneuvering radar target detection based ongeneralizedradon-fouriertransformJ.IEEEtransactionsonaerospace and electronic systems,2016,52(5):2554-2561.DOI:10.1109/TAES.2016.15007615作者简介唐涛(1996)，男，博士，研究方向为无源雷达信号处理技术等.王鹏波(1979)，男，博士，副教授，研究方向为星载 SAR 仿真与成像技术、无源雷达信号处理技术等.陈杰(1973)，男，博士，教授，研究方向为星载SAR 仿真与成像技术、无源雷达信号处理技术等.周新凯(1992)，男，博士，研究方向为无源雷达信号处理技术等.曾虹程(1989)，男，博士，副教授，研究方向为星载 SAR 仿真与成像技术、无源雷达信号处理技术等.Ground moving target imaging method study using GNSS-basedpassive bistatic radarTANGTao，WANGPengbo，CHENJie，ZHOUXinkai，ZENGHongcheng(School of Electronics and Information Engineering,Beihang University,Beijing 100191,China)Abstract：PassivebistaticradarbasedonGlobalNavigationSatelliteSystem(GNSS)isapotentialspace-basedradarsystem,whichhastheadvantagesoflowpowerconsumptionandgoodconcealment.Duetothelowsignalpowerdensity,itisnecessarytoincreasethetargetsignal-to-noiseratiothroughlong-timeintegration.However,targetmotionwillcauserangemigrationandDopplerspread,andthetraditionalpulse第3期唐涛，等：基于 GNSS 的无源双基地雷达地面动目标成像方法研究31dopplerradarprocessingmethodisnolongerapplicable.Aimingatthisproblem,thispaperproposesamovingtargetimagingmethodbasedonRadonFouriertransform(RFT).Firstly,byaddingtheDopplerratecompensationfactortothetraditionalRFTmethod,thethree-dimensionalparametersofthemovingtargetarejointlysearchedandestimated.AndthenthroughparameterscompensationtocompletethetargetimagingintherangeDopplerdomain.ThismethodnotonlyconsiderstheinfluenceofDopplerphasemodulation,butalsorealizesrapiddeploymentintherangefrequencydomain,whichcanachieveabetterimagingeffectthantraditionalmethods.TheperformanceoftheproposedmethodisverifiedbasedontheactualexperimentusingGPSsatelliteastheilluminatorofopportunity.Keywords：GlobalSatelliteNavigationSystem(GNSS);passivebistaticradar;movingtargetimaging;RadonFouriertransform;long-timecoherentintegration(上接第 11 页)WARNOCK A,RUF C.Response to variations in riverflowratebyaspaceborneGNSS-RriverwidthestimatorJ.Remote sensing,2019,11(20):2450.DOI:10.3390/rs1120245071ZEIGERP,FRAPPARTF,DARROZESJ,etal.SNR-Basedwaterheightretrievalinrivers:applicationtohighamplitudeasymmetric tides in the garonne riverJ.Remote sensing,2021,13(9):1856.DOI:10.3390/rs1309185672王峰,杨鹏瑜,杨东凯.GNSS-I/MR河流边界及水位测量:73理论与仿真J/OL.(2023-01-10)2023-01-20.北京航空航天大学学报,2023:1-15.作者简介杨东凯(1972)，男，博士，教授，研究方向为GNSS 反射信号与遥感应用.李杰(1993)，男，博士研究生，研究方向为GNSS 反射信号与陆表参数遥感应用.A review of the application of GNSS reflected signal on landsurface remote sensingYANGDongkai，LIJie(Beihang University,Beijing 100191,China)Abstract：GlobalNavigationSatelliteSystem(GNSS)hasbeendevelopedformorethanhalfacentury,anditcannotonlyprovideuserswithnavigation,positioningandtimingservices,butalsobeusedforremotesensingoftheEarth,anditsapplicationsarebeyondimagination.BasedontheapplicationofGNSS-R,thispapersystematicallyintroducesthebasicconceptsandsignalcharacteristicsofGNSS-R,focusesonthelandsurface remote sensing applications,and analyzes the possible future development directions,which canprovideanimportantreferencefortheapplicationofGNSS-R.Keywords：GlobalNavigationSatelliteSystem;reflectedsignal;landremotesensing;soilmoisture;ve-getationwatercontent;snowdepth;waterbodyidentification;riverboundary32全 球 定 位 系 统第48卷