CN113406636B - 一种曲线运动双基sar的nlcs成像方法 - Google Patents

一种曲线运动双基sar的nlcs成像方法 Download PDF

Info

Publication number
CN113406636B
CN113406636B CN202110677996.XA CN202110677996A CN113406636B CN 113406636 B CN113406636 B CN 113406636B CN 202110677996 A CN202110677996 A CN 202110677996A CN 113406636 B CN113406636 B CN 113406636B
Authority
CN
China
Prior art keywords
azimuth
function
distance
point
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110677996.XA
Other languages
English (en)
Other versions
CN113406636A (zh
Inventor
谭鸽伟
李梦慧
杨晶晶
徐熙毅
李彪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaqiao University
Original Assignee
Huaqiao University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huaqiao University filed Critical Huaqiao University
Priority to CN202110677996.XA priority Critical patent/CN113406636B/zh
Publication of CN113406636A publication Critical patent/CN113406636A/zh
Application granted granted Critical
Publication of CN113406636B publication Critical patent/CN113406636B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004—SAR image acquisition techniques
    • G01S13/9019—Auto-focussing of the SAR signals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004—SAR image acquisition techniques
    • G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004—SAR image acquisition techniques
    • G01S13/9017—SAR image acquisition techniques with time domain processing of the SAR signals in azimuth
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9058—Bistatic or multistatic SAR

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明针对需要多角度观测目标等特殊场景需求，提出了一种曲线运动双基SAR的NLCS成像方法。首先利用Chebyshev多项式分解三维速度和加速度下的斜距函数，目的是提高斜距函数的近似精度；再通过LRCMC补偿线性距离走动项，使大斜视等效为小斜视或正侧视成像，简化后期算法的处理过程；针对三维速度和加速度带来的方位时间空变问题，提出方位向非线性变标算法，利用扰动函数均衡不同点目标的多普勒调频率，使同一距离单元的方位滤波可在相同参数下进行，从而解决了多普勒调频率随方位向的空变问题，并有效提高了点目标的聚焦深度和成像分辨率。

Description

一种曲线运动双基SAR的NLCS成像方法
技术领域
本发明属于雷达信号处理领域，特别涉及合成孔径雷达成像技术中的一种曲线运动双基SAR的NLCS成像方法。
背景技术
双基合成孔径雷达(Bistatic Synthetic Aperture Radar，BiSAR)指发射平台和接收平台分置工作的雷达成像系统。与传统的单基SAR相比，在运用上有更大的灵活性，作用距离更远，能够获取更丰富的信息。但同时收发平台分置的结构也带来了一些技术挑战，双基SAR斜距函数是两个双曲线之和，这使频谱函数的求解更加困难，增加了成像算法的复杂性。
针对实际应用场景中需要多角度观测目标等特殊场景需求时，双基SAR的载体平台往往是存在三维速度和加速度而做曲线轨迹运动，曲线运动轨迹双基SAR的特殊构型具有更大的灵活性、自主性和机动性。三维速度和加速度不仅降低了双基SAR的斜距的精度，且增加了对成像算法的要求。国内张爱丽等人利用Taylor级数近似斜距函数，这种方法相较于正交多项式的方法斜距误差较大；余海等人是针对双基SAR的NLCS算法，但只考虑了平台的三维速度，忽略了因大气湍流等引起的三维加速度带来的误差影响；吴俊杰等人将KT(Keystone Transform)变换应用于NLCS算法中，修正了双基SAR静止场景中运动目标的LRCM，但忽略了剩余RCM，使得聚焦效果在方位和距离向有比较大的衰减；电子科技大学的李忠宇等人提出了一种自适应NLCS技术，校正了高阶RCM并平衡了多普勒参数的空变性，但是该方法只适用于前视SAR的运动目标模式，适用范围受限。李雨等人是将斜距采用Taylor级数展开，方位向通过NLCS算法均衡了不同点目标的调频率，但对合成孔径时间有约束，大的合成孔径时间下会发生频谱混叠现象。
基于上述存在的问题，克服该现有方法所存在的缺陷是本技术领域亟待解决的问题。
发明内容
本发明的目的在于克服现有技术的不足，针对实际应用场景中需要多角度观测目标等特殊场景需求，提出了一种曲线运动双基SAR的NLCS成像方法，有效解决了三维速度和加速度带来的多普勒调频率的方位向空变性问题，并提高点目标的聚焦性能和成像质量。
本发明的技术方案为：
一种曲线运动双基SAR的NLCS成像方法，其特征在于，针对需要多角度观测目标等特殊场景需求，构造曲线运动双基SAR的斜距函数表达式，建立曲线运动双基SAR回波模型；对构造的斜距函数表达式，利用Chebyshev正交多项式近似为关于方位时间的四阶级数形式；其它步骤包括如下：
1)对点目标在距离频域方位时域完成距离压缩，并对LRCM和多普勒线性相位进行有效校正和补偿；
2)经过距离向的预处理后，进行距离向Fourier逆变换，变换到二维时域；
3)提出非线性变标函数，与距离向预处理后的二维时域回波信号相乘，均衡不同点目标的多普勒调频率；
4)作方位向Fourier变换，在距离-多普勒域进行方位向匹配滤波，最后进行方位Fourier逆变换得到双基SAR的点目标聚焦图像。
优选地，针对需要多角度观测目标等特殊场景需求，构造曲线运动双基SAR的斜距函数表达式，建立曲线运动双基SAR回波模型，所述斜距函数表达式表示如下：
RA(t)＝RT(t)+RR(t)
Figure BDA0003121588630000021
Figure BDA0003121588630000022
式中场景区域中心点为A(x0,y0,0)，T和R分别表示发射平台和接收平台，它们分别沿不同路径，以不同三维速度(vTx,vTy,vTz)和(vRx,vRy,vRz)飞行，且具有不同的三维加速度(aTx,aTy,aTz)和(aRx,aRy,aRz)。经过t时刻，发射平台的位置坐标为(xT,yT,zT)，其中xT＝vTxt+aTxt2/2，yT＝vTyt+aTyt2/2，zT＝vTzt+aTzt2/2+HT接收平台的位置坐标为(xR,yR,zR)，其中xR＝vRxt+aRxt2/2，yR＝vRyt+aRyt2/2，zR＝vRzt+aRzt2/2+HR。ai0＝x0+y0+Hi，ai1＝-x0vix-y0viy+Hiviz，
Figure BDA0003121588630000023
ai3＝vixaix+viyaiy+vizaiz，ai4＝(aix+aiy+aiz)2/4。
曲线运动双基SAR回波信号表示为
Figure BDA0003121588630000024
式中A0是点目标的后向散射系数，τ和t∈[-Ts/2,Ts/2]分别为距离向快时间和方位向慢时间，Ts是合成孔径时间；ωr(·)和ωaz(·)分别为距离窗函数和方位窗函数；c为电磁波传播速度；fc为载波频率；Kr为发射信号的调频率；RA(t)为中心点目标A的瞬时斜距总历程。。
优选地，对构造的斜距函数表达式，利用Chebyshev正交多项式近似为关于方位时间的四阶级数形式步骤中，具体包括：
对方位向时间t作归一化处理，即令
Figure BDA0003121588630000031
x∈[-1,1]，斜距函数表示为：
Figure BDA0003121588630000032
/>
分别对RT(x)和RR(x)进行Chebyshev正交分解，并按x的幂级数整理为：
RT(x)＝αT0+αT1x+αT2x2+αT3x3+αT4x4
RR(x)＝αR0+αR1x+αR2x2+αR3x3+αR4x4
其中，分解系数
Figure BDA0003121588630000033
αi1＝Ci1-3Ci3，αi2＝2Ci2-8Ci4，αi3＝4Ci3,αi4＝8Ci4(式中Cij是Chebyshev系数，j＝0,1,2,3,4，i指T或R)，且有
Figure BDA0003121588630000034
其中，Tj(x)是Chebyshev多项式递推式，且有T0(x)＝1，T1(x)＝x，...，Tj(x)＝2xTj-1(x)-Tj-2(x)，
Figure BDA0003121588630000035
为变量节点，n＝4是展开阶数，
Figure BDA0003121588630000036
将
Figure BDA0003121588630000037
代入归一化后的斜距函数表达式，并按t的四阶幂级数整理，得到恒加速度作用下双基SAR斜距的Chebyshev分解式为
R(t)＝λ0+λ1t+λ2t2+λ3t3+λ4t4+…
其中，λ0＝αT0+αR0，
Figure BDA0003121588630000038
Figure BDA0003121588630000039
式中
λ0＝[(CT0/2-CT2+CT4)+(CR0/2-CR2+CR4)]
Figure BDA0003121588630000041
Figure BDA0003121588630000042
Figure BDA0003121588630000043
其中λ0是常数项；λ1是距离徙动项系数，用于线性相位去除；λ2是距离弯曲项系数，用于控制方位频率调制(FM)速率，同时也可以用于扰动函数的推导；λ3和λ4是高阶项系数，用于使用级数反演法推导二位频谱并定义方位向匹配滤波器。
优选地，在距离频域方位时域完成距离压缩，并对LRCM和多普勒线性相位进行有效校正和补偿步骤中，具体包括：
双向的距离偏移量在假设t＝0处计算：
Figure BDA0003121588630000044
/>
LRCMC和距离压缩对应补偿函数：
Figure BDA0003121588630000045
此时中心点A的斜距函数为：
RA1(t)＝λ0+λ2t2+λ3t3+λ4t4+…
选出参考点C进行后续的分析：
RC1(t)＝λ0+λ2(t-tC)2+λ3(t-tC)3+λ4(t-tC)4+…
线性相位去除对应的补偿函数：
Figure BDA0003121588630000046
优选地，经过距离向的预处理后，进行距离向Fourier逆变换，变换到二维时域步骤中，具体包括：
经过LRCMC和线性相位去除后的信号为：
Figure BDA0003121588630000047
对A点和C点预处理后的信号做距离向IFFT，得到时域回波信号为：
Figure BDA0003121588630000051
Figure BDA0003121588630000052
优选地，提出非线性变标函数，与距离向预处理后的二维时域回波信号相乘，均衡不同点目标的多普勒调频率步骤中，具体包括：
提出一个扰动函数exp{jπαt3}，它是关于方位时间t的三次指数函数，将扰动函数与回波信号相乘从而实现均衡不同点目标调频率的目的。
将A点的回波信号与扰动函数相乘，并忽略斜距的常数项和高于二次项的系数可得：
Figure BDA0003121588630000053
同样，将参考点C的回波信号与扰动函数相乘，并忽略斜距的常数项和高于二次项的系数可得：
Figure BDA0003121588630000054
求解出扰动函数系数为
Figure BDA0003121588630000055
优选地，作方位向Fourier变换，在距离-多普勒域进行方位向匹配滤波，最后进行方位Fourier逆变换得到双基SAR聚焦图像步骤中，具体包括：
目标A的高阶项可通过执行方位FFT(忽略常数项)进行处理：
Figure BDA0003121588630000056
利用驻留相位原理求得驻留相位点，并将驻留相位点代入方位相位中，可得到方位向补偿函数为：
Figure BDA0003121588630000057
Figure BDA0003121588630000058
经过方位相位匹配滤波后，将回波信号变换到二维时域，即得到双基SAR的点目标聚焦图像。
采用上述方案，本发明的有益效果：
在需要多角度观测目标等特殊应用场景中，BiSAR载体平台往往存在三维速度和加速度而作曲线轨迹运动，三维速度和加速度不仅降低了BiSAR的斜距的精度，且增加了对成像算法的要求。本发明提出了一种曲线运动双基SAR的NLCS成像方法，利用Chebyshev多项式分解三维速度和加速度下的斜距函数，相比于Taylor近似斜距的方法，大大提高了斜距函数的近似精度；通过LRCMC补偿线性距离走动项和线性相位，使大斜视等效为小斜视或正侧视成像，简化后期算法的处理过程，降低了距离-方位耦合；针对三维速度和加速度带来的方位时间空变问题，提出方位向非线性变标算法，利用扰动函数解决了多普勒调频率的方位向空变问题，使同一距离单元点目标可以采用相同的方位补偿函数进行匹配滤波。
NLCS算法是时频域混合算法，算法只包含FFT和复数乘法运算，不涉及插值处理，因此易于工程实时应用。扰动函数是关于方位时间的三次函数，在时域与回波信号相乘，需要逐点进行处理，运算效率低于频域算法。通过仿真实验证明了NLCS算法不仅对中心点能够良好成像，而且对边缘点也可以实现高分辨率的成像效果，可见该算法可以很好的改善边缘点目标的成像质量；且基于Chebyshev分解的BiSAR的NLCS算法不仅适用于曲线运动BiSAR，在直线运动模式下，仍然能够实现高分辨率成像。综上所述，基于Chebyshev分解的NLCS算法应用场景广泛，成像质量高，是一种值得广泛推广的成像算法。
附图说明
图1是本发明具体实施方案采用的双基SAR几何构型图；
图2是本发明双基SAR的NLCS成像方法流程图；
图3是LRCMC和线性相位去除示意图；
图4是扰动函数均衡调频率的原理图；
图5是本发明方法的多点成像图；
图6是本发明方法中心点和边缘点等高线图；
图7是本发明方法中心点和边缘点的距离向脉冲响应图；
图8是本发明方法中心点和边缘点的方位向脉冲响应图。
具体实施方式
本发明主要采用仿真实验的方法进行验证，所有步骤、结论都在Matlab2016a上验证正确。下面结合发明附图和具体实施方式对本发明作进一步的详细描述。
如图2所示，本发明是一种基于Chebyshev正交分解的曲线运动双基SAR的NLCS(Non-Linear Chirp Scaling)成像方法，针对需要多角度观测目标等特殊场景需求，构造曲线运动双基SAR的斜距函数表达式，建立曲线运动双基SAR回波模型。
具体的，本发明具体实施方案采用的双基SAR几何构型如图1所示，其中：
场景区域中心点为A(x0,y0,0)，T和R分别表示发射平台和接收平台，它们分别沿不同路径，以不同三维速度(vTx,vTy,vTz)和(vRx,vRy,vRz)飞行，且具有不同的三维加速度(aTx,aTy,aTz)和(aRx,aRy,aRz)。经过t时刻，发射平台的位置坐标为(xT,yT,zT)，其中xT＝vTxt+aTxt2/2，yT＝vTyt+aTyt2/2，zT＝vTzt+aTzt2/2+HT接收平台的位置坐标为(xR,yR,zR)，其中xR＝vRxt+aRxt2/2，yR＝vRyt+aRyt2/2，zR＝vRzt+aRzt2/2+HR。
曲线运动轨迹双基SAR斜距函数表达式如下：
RA(t)＝RT(t)+RR(t)
Figure BDA0003121588630000071
Figure BDA0003121588630000072
式中ai0＝x0+y0+Hi，ai1＝-x0vix-y0viy+Hiviz，ai2＝vi2x+vi2y+vi2z-x0aix-y0aiy+Hiaiz，ai3＝vixaix+viyaiy+vizaiz，ai4＝(aix+aiy+aiz)2/4。
曲线运动双基SAR回波信号表示为
Figure BDA0003121588630000073
其中，A0是点目标的后向散射系数，τ和t∈[-Ts/2,Ts/2]分别为距离向快时间和方位向慢时间，Ts是合成孔径时间；ωr(·)和ωaz(·)分别为距离窗函数和方位窗函数；λ为载波波长；c表示光速；Kr为发射信号的调频率；RA(t)为中心点目标A的瞬时斜距总历程。
进一步的，对构造的斜距函数表达式，利用Chebyshev正交多项式近似为关于方位时间的四阶级数形式。
具体的，将Chebyshev多项式应用在曲线运动双基SAR斜距函数上，首先对方位向时间t作归一化处理，即令
Figure BDA0003121588630000074
x∈[-1,1]，斜距函数表示为：/>
Figure BDA0003121588630000081
然后分别对RT(x)和RR(x)进行Chebyshev正交分解，并按x的幂级数整理为：
RT(x)＝αT0+αT1x+αT2x2+αT3x3+αT4x4
RR(x)＝αR0+αR1x+αR2x2+αR3x3+αR4x4
其中，分解系数
Figure BDA0003121588630000082
αi1＝Ci1-3Ci3，αi2＝2Ci2-8Ci4，αi3＝4Ci3,αi4＝8Ci4(式中Cij是Chebyshev系数，j＝0,1,2,3,4，i指T或R)，且有
Figure BDA0003121588630000083
其中，Tj(x)是Chebyshev多项式递推式，且有T0(x)＝1，T1(x)＝x，...，Tj(x)＝2xTj-1(x)-Tj-2(x)，
Figure BDA0003121588630000084
为变量节点，n＝4是展开阶数，
Figure BDA0003121588630000085
将
Figure BDA0003121588630000086
代入归一化后的斜距函数表达式，并按t的四阶幂级数整理，得到三维加速度作用下的曲线运动双基SAR斜距的Chebyshev分解式为
R(t)＝λ0+λ1t+λ2t2+λ3t3+λ4t4+…
其中，λ0＝αT0+αR0，
Figure BDA0003121588630000087
Figure BDA0003121588630000088
本发明的其它步骤包括包括如下：
1)对点目标在距离频域方位时域完成距离压缩，并对LRCM和多普勒线性相位进行有效校正和补偿。
具体的，由于点目标的斜距轨迹中一般包含线性距离徙动分量(LRCM)和非线性距离徙动分量，当复合天线波束斜视时，会出现线性分量，而大部分的距离-多普勒耦合来自线性分量，因此需要通过LRCMC操作去除LRCM分量，使得目标轨迹与方位向坐标轴对齐，以便于方位压缩。
对原始回波信号作距离向FFT，得到距离频域-方位时域信号为：
Figure BDA0003121588630000091
LRCM的斜率随着距离的变化而变化，因为斜视角随着目标距离的变化而变化，在随方位变化的双基地情况下，斜视角也随方位变化。为了解决这个问题，需要在不变性区域内进行LRCMC，以保持斜视角度的变化较小。双向的距离偏移量在假设t＝0处计算：
Figure BDA0003121588630000092
LRCMC和距离压缩对应补偿函数为，其中fτ表示距离向频率：
Figure BDA0003121588630000093
此时中心点A的斜距函数为：
RA1(t)＝λ0+λ2t2+λ3t3+λ4t4+…
选出参考点C进行后续的分析：
RC1(t)＝λ0+λ2(t-tC)2+λ3(t-tC)3+λ4(t-tC)4+…
其中，其中tC表示参考点C的方位时间。
LRCMC只是去除了斜距函数中的线性平移分量，多普勒频移引起的线性相位项仍然存在，因此应该将此线性相位项去除，以便于MSR的应用，从而消除距离和方位向的耦合。
线性相位去除对应的补偿函数：
Figure BDA0003121588630000094
经过LRCMC和线性相位去除后的信号为，其中fc表示载波频率：
Figure BDA0003121588630000095
2)经过距离向的预处理后，进行距离向Fourier逆变换，变换到二维时域。
具体的，对A点和C点预处理后的信号做距离向IFFT，得到时域回波信号为：
Figure BDA0003121588630000096
Figure BDA0003121588630000101
其中ta表示A点的方位时间。
3)提出非线性变标函数，与距离向预处理后的二维时域回波信号相乘，均衡不同点目标的多普勒调频率。
具体的，经过LRCMC和距离压缩后，不同方位向调频率的点目标位于相同的距离单元内。为了能够在多普勒频域中实现方位向压缩，在时域中应用NLCS算法，以均衡不同点目标的调频率，提出一个扰动函数exp{jπαt3}，它是关于方位时间t的三次指数函数，将扰动函数与信号相乘从而实现均衡不同点目标调频率的目的。
将A点的回波信号与扰动函数相乘，并忽略斜距的常数项和高于二次项的系数可得：
Figure BDA0003121588630000102
由于上式中的近似仅包括二阶相位项，因此NLCS算法只对目标的线性调频分量进行均衡。
同样，将参考点C的回波信号与扰动函数相乘，并忽略斜距的常数项和高于二次项的系数可得：
Figure BDA0003121588630000103
令tn＝t-tC，并带入spertC(τ,ta)中得：
Figure BDA0003121588630000104
为了求出扰动系数α，令tC为0，同时令
Figure BDA0003121588630000105
和/>
Figure BDA0003121588630000106
相加为零，即
Figure BDA0003121588630000107
可以求解出扰动系数/>
Figure BDA0003121588630000108
式中
RT(0)+RR(0)＝λ1tC＝2[(CT1-3CT3)+(CR1-3CR3)]tC/Ts
则
Figure BDA0003121588630000111
spertC(τ,tn)中第一个相位项是由微扰动过程引起的，第二项是小的多普勒频移，第三项是常数项，取决于目标的位置，这对聚焦过程没有影响，第四项是线性调频调制项，对所有目标都一样。所以，最终方位相位调制包含在第一个和第四个指数项中，并且这种调制对于同一距离单元中的所有目标都是相同的。
5)作方位向Fourier变换，在距离-多普勒域进行方位向匹配滤波，最后进行方位Fourier逆变换得到双基SAR的点目标聚焦图像。
具体的，目标A的高阶项可通过执行方位FFT(忽略常数项)进行处理：
Figure BDA0003121588630000112
利用驻留相位原理和级数反演法求得驻留相位点：
Figure BDA0003121588630000113
令
Figure BDA0003121588630000114
得到fa与方位时间t的关系式为：
Figure BDA0003121588630000115
利用级数反演法求得驻留相位点为：
Figure BDA0003121588630000116
其中，
Figure BDA0003121588630000117
将驻留相位点代入方位相位中，可得到方位向补偿函数为：
Figure BDA0003121588630000118
Figure BDA0003121588630000119
经过方位相位的匹配滤波后作方位向IFFT得到点目标聚焦图像。
下面将通过实验仿真对本发明方法加以说明。
图2是本发明NLCS算法流程图。首先对点目标进行预处理，在距离频域作LRCMC和线性相位去除，为了均衡不同点目标的调频率，提出扰动函数，将预处理后的点目标回波信号与扰动函数相乘，实现均衡不同点目标调频率的目的，最后在方位频域进行方位相位的补偿实现点目标的高分辨率成像。
图3是LRCMC和线性相位去除示意图。如图3(a)所示，C与E有相同的最近斜距，B与D有相同的最近斜距，B和E有相同的波束中心。经过LRCMC和线性相位去除后，如图3(b)所示，A，C，D位于相同的距离单元内，且具有相同的斜率，LRCMC操作的目的是将具有不同方位向调频率的点目标移动到同一个距离单元中，以便于利用扰动函数均衡不同点目标得调频率，从而补偿多普勒调频率的方位向空变形。
图4是利用扰动函数均衡调频率的原理图。为了简化示意图，图中没有显示由斜视引起的多普勒偏移量。图4中，(a)为三个点目标的方位向调频率信号的实部，图(b)为各方位向信号对应的相位。从图中可知，三个曲线具有不同的方位向调频率，这里调频率指信号实部曲线的二阶导数，该参数随着目标的原始距离位置的变化而变化。为了均衡目标的调频率，提出了关于方位时间三次方的扰动函数，它的曲线图如图(c)所示。将此立方相位添加到图(b)相位中，得到图(d)的三个目标相位，根据指数函数的性质，点目标的相位被改变，并具有相同的调频率。图(e)是经过扰动之后的点目标信号的实部。
为了验证NLCS算法在BiSAR系统的有效性和可行性，采用表1中的数据进行实验仿真。发射平台和接收平台的三维速度和加速度在表2中列出。
表1仿真参数
Figure BDA0003121588630000121
表2速度和加速度
Figure BDA0003121588630000122
图5是基于Chebyshev近似的NLCS算法25点目标等高线图，从图中可以看出，不管是场景中心点，还是场景边缘点，都可以实现良好的成像效果，可见该算法可以改善边缘点的聚焦效果。同时，还可以发现一个规律，位于同一个距离单元的点目标的成像效果相同，这是因为首先进行了线性距离徙动校正，将不同距离单元的点目标校正在了相同的距离单元，最后在距离单元中应用相同的方位向匹配滤波器，从而实现了NLCS算法的高分辨率成像。
图6是中心点P13、近端点P23和边缘点P25的成像等高线图。从图中可以看出，即使在三维速度和加速度的作用下，中心点和边缘点都有着良好的成像质量和聚焦效果。
图7是中心点P13、近端点P23和边缘点P25的距离向脉冲响应图。从图中可以看出，中心点和边缘点的距离向脉冲响应都接近理想的冲激响应；图8是中心点P13、近端点P23和边缘点P25的方位向脉冲响应图。从图中可以看出，中心点和边缘点的方位向脉冲响应都较为理想，且中心点P13的主瓣宽度较窄，说明方位分辨率较高。综上所述，本发明基于Chebyshev正交分解的NLCS算法对目标的距离向和方位向都可以实现良好的聚焦效果。
为了评估本发明的基于Chebyshev多项式的NLCS算法的成像性能，分别计算了中心点P13、近端点P23和边缘点P25的距离和方位向的积分旁瓣比(Integral Sidelobe Ratio,ISLR)和峰值旁瓣比(Peak Sidelobe Ratio,PSLR)，以及分辨率。如表3所示，通过数值分析和对比，本发明算法的PSLR和ISLR测量值接近理论值，中心点和边缘点的距离向的分辨率近似，中心点的方位向分辨率最高。
表3曲线运动模式下两种方法成像质量评估
Figure BDA0003121588630000131
综上所述，本发明的方法，利用Chebyshev多项式分解三维速度和加速度下的斜距函数，相比于Taylor近似斜距的方法，大大提高了斜距函数的近似精度；通过LRCMC补偿线性距离走动项和线性相位，使大斜视等效为小斜视或正侧视成像，简化后期算法的处理过程，降低了距离-方位耦合；针对三维速度和加速度带来的方位时间空变问题，提出方位向非线性变标算法，利用扰动函数解决了多普勒调频率的方位向空变问题，使同一距离单元点目标可以采用相同的方位补偿函数进行匹配滤波。通过仿真实验证明了，NLCS算法不仅对中心点能够良好成像，而且对边缘点也可以实现高分辨率的成像效果，可见该算法可以很好的改善边缘点目标的成像质量，且应用范围广泛。
上述仅为本发明的具体实施方式，但本发明的设计构思并不局限于此，凡利用此构思对本发明进行非实质性的改动，均应属于侵犯本发明保护范围的行为。

Claims (5)

1.一种曲线运动双基SAR的NLCS成像方法，其特征在于，针对需要多角度观测目标场景需求，构造曲线运动双基SAR的斜距函数表达式，建立曲线运动双基SAR回波模型；对构造的斜距函数表达式，利用Chebyshev正交多项式近似为关于方位时间的四阶级数形式；其余步骤如下：
1)对点目标在距离频域方位时域完成距离压缩，并对LRCM和多普勒线性相位进行有效校正和补偿；
2)经过距离向的预处理后，进行距离向Fourier逆变换，变换到二维时域，具体包括如下：
经过LRCMC和线性相位去除后的信号为，其中fc为载波频率：
Figure FDA0004171404970000011
对A点和C点预处理后的信号做距离向IFFT，得到时域回波信号为，其中ta表示A点的方位时间：
Figure FDA0004171404970000012
Figure FDA0004171404970000013
3)提出非线性变标函数，与距离向预处理后的二维时域回波信号相乘，均衡不同点目标的多普勒调频率，具体包括如下：提出一个扰动函数exp{jπαt3}，它是关于方位时间t的三次指数函数，将扰动函数与时域回波信号相乘从而实现均衡不同点目标调频率的目的；将A点的时域回波信号与扰动函数相乘，并忽略斜距的常数项和高于二次项的系数可得：
Figure FDA0004171404970000014
同样，将参考点C的回波信号与扰动函数相乘，并忽略斜距的常数项和高于二次项的系数可得：
Figure FDA0004171404970000015
求解出扰动函数系数为
Figure FDA0004171404970000021
4)作方位向Fourier变换，在距离-多普勒域进行方位向匹配滤波，最后进行方位Fourier逆变换得到双基SAR的点目标聚焦图像，具体包括如下：
目标A的高阶项可通过执行方位FFT进行处理：
Figure FDA0004171404970000022
利用驻留相位原理求得驻留相位点，并将驻留相位点代入方位相位中，可得到方位向补偿函数为：
Figure FDA0004171404970000023
Figure FDA0004171404970000024
其中，
Figure FDA0004171404970000025
表示方位补偿函数的相位，t(fa)表示方位向的驻留相位点，fa表示方位向频率，经过方位相位匹配滤波后，将回波信号变换到二维时域，即得到双基SAR的点目标聚焦图像。
2.根据权利要求1所述的一种曲线运动双基SAR的NLCS成像方法，其特征在于，所述斜距函数表达式表示如下：
RA(t)＝RT(t)+RR(t)
Figure FDA0004171404970000026
Figure FDA0004171404970000027
式中场景区域中心点为A(x0,y0,0)，T和R分别表示发射平台和接收平台，它们分别沿不同路径，以不同三维速度(vTx,vTy,vTz)和(vRx,vRy,vRz)飞行，且具有不同的三维加速度(aTx,aTy,aTz)和(aRx,aRy,aRz)；经过t时刻，发射平台的位置坐标为(xT,yT,zT)，其中xT＝vTxt+aTxt2/2，yT＝vTyt+aTyt2/2，zT＝vTzt+aTzt2/2+HT接收平台的位置坐标为(xR,yR,zR)，其中xR＝vRxt+aRxt2/2，yR＝vRyt+aRyt2/2，zR＝vRzt+aRzt2/2+HR，ai0＝x0+y0+Hi，ai1＝-x0vix-y0viy+Hiviz，
Figure FDA0004171404970000028
ai3＝vixaix+viyaiy+vizaiz，ai4＝(aix+aiy+aiz)2/4,i＝T,R。
3.如权利要求1所述的一种曲线运动双基SAR的NLCS成像方法，其特征在于，所述曲线运动双基SAR回波模型表示为
Figure FDA0004171404970000031
其中，A0是点目标的后向散射系数，τ和t∈[-Ts/2,Ts/2]分别为距离向快时间和方位向慢时间，Ts是合成孔径时间；ωr(·)和ωaz(·)分别为距离窗函数和方位窗函数；c是光速；λ为载波波长；Kr为发射信号的调频率；RA(t)为中心点目标A的瞬时斜距总历程。
4.根据权利要求2所述的一种曲线运动双基SAR的NLCS成像方法，其特征在于，对构造的斜距函数表达式，利用Chebyshev正交多项式近似为关于方位时间的四阶级数形式步骤中，具体包括：
对方位向时间t作归一化处理，即令
Figure FDA0004171404970000032
Ts表示合成孔径时间，斜距函数表示为：/>
Figure FDA0004171404970000033
分别对RT(x)和RR(x)进行Chebyshev正交分解，并按x的幂级数整理为：
RT(x)＝αT0+αT1x+αT2x2+αT3x3+αT4x4
RR(x)＝αR0+αR1x+αR2x2+αR3x3+αR4x4
其中，分解系数
Figure FDA0004171404970000034
αi1＝Ci1-3Ci3，αi2＝2Ci2-8Ci4，αi3＝4Ci3,αi4＝8Ci4，式中Cij是Chebyshev系数，j＝0,1,2,3,4，i指T或R，且有
Figure FDA0004171404970000035
其中，Tj(x)是Chebyshev多项式递推式，且有T0(x)＝1，T1(x)＝x，...，Tj(x)＝2xTj-1(x)-Tj-2(x)，
Figure FDA0004171404970000036
为变量节点，n＝4是展开阶数，
Figure FDA0004171404970000037
将
Figure FDA0004171404970000041
代入归一化后的斜距函数表达式，并按t的四阶幂级数整理，得到恒加速度作用下双基SAR斜距的Chebyshev分解式为
R(t)＝λ0+λ1t+λ2t2+λ3t3+λ4t4+…
其中，λ0＝αT0+αR0，
Figure FDA0004171404970000042
Figure FDA0004171404970000043
5.根据权利要求4所述一种曲线运动双基SAR的NLCS成像方法，其特征在于，步骤1)具体包括如下：
双向的距离偏移量在假设t＝0处计算：
Figure FDA0004171404970000044
LRCMC和距离压缩对应补偿函数为，其中fτ表示距离向频率：
Figure FDA0004171404970000045
此时中心点A的斜距函数为：
RA1(t)＝λ0+λ2t2+λ3t3+λ4t4+…
选出参考点C进行后续的分析，其中tC表示参考点C的方位时间：
RC1(t)＝λ0+λ2(t-tC)2+λ3(t-tC)3+λ4(t-tC)4+…
线性相位去除对应的补偿函数：
Figure FDA0004171404970000046
/>
CN202110677996.XA 2021-06-18 2021-06-18 一种曲线运动双基sar的nlcs成像方法 Active CN113406636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110677996.XA CN113406636B (zh) 2021-06-18 2021-06-18 一种曲线运动双基sar的nlcs成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110677996.XA CN113406636B (zh) 2021-06-18 2021-06-18 一种曲线运动双基sar的nlcs成像方法

Publications (2)

Publication Number Publication Date
CN113406636A CN113406636A (zh) 2021-09-17
CN113406636B true CN113406636B (zh) 2023-05-26

Family

ID=77681356

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110677996.XA Active CN113406636B (zh) 2021-06-18 2021-06-18 一种曲线运动双基sar的nlcs成像方法

Country Status (1)

Country Link
CN (1) CN113406636B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7482965B1 (en) * 2004-12-17 2009-01-27 Mayo Foundation For Medical Education And Research Digital chirp waveform generator and method
CN102176016A (zh) * 2011-01-25 2011-09-07 北京航空航天大学 一种大斜视滑动聚束sar成像处理方法
EP2725382A1 (de) * 2012-10-26 2014-04-30 Astrium GmbH Synthetisches Apertur Radar zur gleichzeitigen Bildaufnahme und Bewegtzielerkennung
CN104931965A (zh) * 2015-05-04 2015-09-23 电子科技大学 一种固定发射站双基正侧视sar成像方法
CN105005032A (zh) * 2015-07-21 2015-10-28 电子科技大学 一种基于级数反演成像算法的sar移频干扰方法
CN108710111A (zh) * 2018-02-05 2018-10-26 西安电子科技大学 一种机载双基前视sar方位相位的二维空变校正方法
CN111208514A (zh) * 2020-01-06 2020-05-29 华侨大学 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7482965B1 (en) * 2004-12-17 2009-01-27 Mayo Foundation For Medical Education And Research Digital chirp waveform generator and method
CN102176016A (zh) * 2011-01-25 2011-09-07 北京航空航天大学 一种大斜视滑动聚束sar成像处理方法
EP2725382A1 (de) * 2012-10-26 2014-04-30 Astrium GmbH Synthetisches Apertur Radar zur gleichzeitigen Bildaufnahme und Bewegtzielerkennung
CN104931965A (zh) * 2015-05-04 2015-09-23 电子科技大学 一种固定发射站双基正侧视sar成像方法
CN105005032A (zh) * 2015-07-21 2015-10-28 电子科技大学 一种基于级数反演成像算法的sar移频干扰方法
CN108710111A (zh) * 2018-02-05 2018-10-26 西安电子科技大学 一种机载双基前视sar方位相位的二维空变校正方法
CN111208514A (zh) * 2020-01-06 2020-05-29 华侨大学 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
High-resolution squinted sar imaging algorithm applied in the near-field environment;Wang S等;《 Complexity》;第 1-13页 *
基于贝叶斯统计学习的雷达超分辨成像技术研究;岳云泽 等;《万方学位论文在线期刊库》;第1-97页 *
基于运动补偿和正交解耦合的双基SAR成像算法;李梦慧 等;《信号处理》;第37卷(第01期);第75-85页 *
机载非平飞双基SAR高分辨成像算法研究;陈国瑾;《中国优秀硕士学位论文全文数据库 信息科技辑》(第02期);第I136-1386页 *

Also Published As

Publication number Publication date
CN113406636A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
CN106249237B (zh) 一种曲线轨迹下大斜视sar频域成像方法
CN102393518B (zh) 一种适用于大斜视角的机载sar成像方法
CN108459321B (zh) 基于距离-方位圆模型的大斜视高分辨率sar成像方法
CN106970386A (zh) 一种雷达多普勒波束锐化的优化方法
CN102749620B (zh) 一种弹载/机载雷达单脉冲前视成像处理方法
CN109270528B (zh) 基于全解析距离模型的一站固定式双站sar成像方法
CN113589285B (zh) 一种飞行器sar实时成像方法
CN104865573A (zh) 一种双基前视高机动平台sar扩展场景成像优化方法
CN109507667B (zh) 基于导航卫星信号的双站sar接收机运动误差补偿方法
CN109613532B (zh) 一种机载雷达实时多普勒波束锐化超分辨成像方法
CN106291489B (zh) 适用于多种发射信号波形的合成孔径雷达回波仿真方法
CN108061890B (zh) 一种sar成像方法
CN103018739A (zh) 一种校正多通道幅相误差的微波三维成像方法
CN108226925A (zh) 一种适用于弹载大前斜视时变参数sar的重叠子孔径算法
CN106950565A (zh) 星载sar成像抖动补偿方法、成像方法
CN114545411A (zh) 一种基于工程实现的极坐标格式多模高分辨sar成像方法
CN103064084A (zh) 基于距离频域的解模糊方法
CN111722225B (zh) 基于先验相位结构信息的双基sar两维自聚焦方法
CN110879391B (zh) 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法
CN111487614A (zh) 基于子孔径的曲线航迹弹载sar波前重建成像方法及系统
CN111208514A (zh) 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法
CN107271996B (zh) 一种机载cssar地面运动目标成像方法
CN110441771B (zh) 基于方位时间重采样的高速机动飞行sar成像方法
CN112285658A (zh) 基于递推计算的星载sar欺骗干扰信号快速生成方法
Jun et al. Principle and methods on bistatic SAR signal processing via time correlation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant