一种CSAR地面动目标重聚焦成像方法与流程

本发明属于合成孔径雷达(syntheticapertureradar，sar)地面动目标指示(groundmovingtargetindication，gmti)领域，涉及一种适用于曲线合成孔径雷达(curvesar，csar)的地面动目标重聚焦(groundmovingtargetrefocusing，gmtr)成像方法。
背景技术：
：合成孔径雷达(syntheticapertureradar，sar)是一种可对观测场景进行高分辨微波成像的雷达技术，csar是一种新的sar成像体制，即雷达平台围绕观测场景做大曲线或宽角度圆弧运动，并且波束始终指向目标场景进行观测成像的雷达系统。csar可以实现地面静止目标亚波长级分辨率，并且能够多角度全方位观测目标，获取更加详细的目标方位散射信息。为了提高战场感知能力，人们希望csar在完成静止目标成像侦察的同时，能够实现对地面动目标的侦察探测，即具有gmti功能。csar-gmti能够同时完成对静止/运动目标成像侦察，极大的拓展了sar技术的使用范围。然而，常规csar成像方法只能对静止目标进行精确成像，地面动目标自身的运动会引入额外的距离走动以及距离弯曲，导致动目标在csar图像中发生散焦和位置偏移现象。地面动目标重聚焦成像能够提供动目标的准确情报信息，因此地面动目标重聚焦成像是gmti重要研究内容之一。地面动目标一般为非合作目标，很难准确获取其运动参数信息，增加了重聚焦成像的难度。此外，csar的特殊运动几何在带来全方位多孔径角观测优势的同时，也引入了复杂的方位和距离耦合，这使得传统直线sar下的动目标重聚焦方法失效。如何在csar工作模式下，实现非合作地面动目标精确重聚焦成像是一个亟待解决的技术问题。技术实现要素：本发明的目的在于提供一种适用于csar的地面动目标重聚焦方法，以提高csar-gmti性能及其实用价值。为了实现上述技术目的，本发明的技术方案是：一种csar地面动目标重聚焦成像方法，包括以下步骤：s1、接收动目标的雷达回波信号，并进行预处理；s2、将预处理后的雷达回波信号沿方位向划分成若干个子孔径回波信号；s3、对于某个子孔径回波信号，将子孔径回波信号变换到距离频域，逐频点提取信号沿方位向的相位；s4、根据提取的相位计算得到距离估计值，对距离估计值用二次曲线拟合，根据拟合曲线构造相位补偿因子；s5、利用相位补偿因子进行运动补偿；s6、利用map-drift自聚焦方法补偿相位估计误差；s7、对步骤s6的结果分别沿方位向做傅里叶变换和沿距离向做逆傅里叶变换，获得该子孔径下动目标重聚焦结果；s8、重复步骤s3-s7，逐一对所有子孔径回波信号进行处理，获得多子孔径下动目标重聚焦图像序列。优选地，所述预处理包括距离向脉冲压缩和地杂波抑制。优选地，所述步骤s3的具体过程为：设n表示动目标散射点总个数，n表示散射点序号，n取值范围为0,1,2,…,n-1；n、n均为整数，选取n＝0散射点为参考点，an表示动目标第n个散射点回波信号幅度，rn(ta)表示雷达与动目标第n个散射点之间的瞬时斜距，δrn(ta)＝rn(ta)-r0(ta)，f为信号频率，c表示光速，π表示圆周率，rect(·)表示矩形窗函数，ta表示慢时间，tsub表示第k个子孔径的持续时间，ta,k表示第k个子孔径的中心时刻，则第k个子孔径回波距离频域信号ssk(f,ta)为：将δrn(ta)利用小角近似处理，得到：δrn(ta)＝rn(ta)-r0(ta)≈yn-xnω(ta-ta,k)，其中ta∈[ta,k-tsub/2,ta,k+tsub/2]，ω表示相对旋转角速度，则进一步地，ssk(f,ta)的相位估计值为：其中∠(·)表示计算复信号幅角，unwrap(·)为解缠绕操作。优选地，所述步骤s4的具体过程为：根据提取的相位计算距离估计值为：其中为的离散化形式，nr表示距离频率采样点数，表示距离频率采样点，即距离频率f的离散采样对应的频率，nr取值为0,1,2,…,nr-1，利用二次曲线对进行拟合，得到r0(ta)的估计结果分别表示常数项、一次项、二次项拟合系数，利用上式构造相位补偿因子h(f,ta)：优选地，所述步骤s4的具体过程为：将相位补偿因子与第k个子孔径回波信号相乘，得到运动补偿之后的第k个子孔径回波信号。为了更好的理解本发明技术方案，下面对相关原理和推导过程作详细说明。如图1所示，以笛卡尔坐标系原点为探测场景中心，包含多个散射点的动目标在场景中以速度运动，vx,vy分别表示x轴、y轴方向的速度分量；动目标第n个散射点的位置为xn,yn表示空间坐标值，其中ta表示慢时间，t表示向量转置符号。雷达平台以速度围绕探测场景做曲线运动，慢时间ta时刻的雷达平台的位置坐标为h表示雷达平台与探测场景的垂直距离。雷达与动目标第n个散射点之间的瞬时斜距可以表示为其中||·||2表示2范数。不失一般性，选n＝0散射点为参考点，则雷达发射信号为线性调频信号(linearfrequencymodulation，lfm)，中心频率为fc，带宽为b。接收到的csar动目标回波信号经正交解调、脉冲压缩以及抑制地杂波之后，信号表示为s(r,ta)：其中r表示斜距，sinc(·)表示辛格函数，c表示光速，n表示动目标散射点总个数，an表示动目标第n个散射点回波信号幅度，j为虚数单位，exp表示以自然对数e为底指数函数，相应的距离频域信号g(f,ta)表示为其中f∈[fc-b/2,fc+b/2]，表示傅里叶变换。为便于动目标成像处理，雷达全孔径回波被均匀划分为k个子孔径回波。第k个子孔径回波为：其中tsub和ta,k分别表示第k个子孔径的持续时间和中心时刻，k为整数且0≤k≤k-1。rect(·)表示矩形窗函数。如图2所示，在子孔径下，以参考点为圆心，建立新的坐标系o'x'y'，在不引起混淆的情况下，仍然用(xn,yn)表示第n个散射点在新坐标系o'x'y'下的位置。在o'x'y'坐标系下，利用小角近似，δrn(ta)表示为：δrn(ta)＝rn(ta)-r0(ta)≈yn-xnω(ta-ta,k)(6)其中ta∈[ta,k-tsub/2,ta,k+tsub/2]，ω表示雷达与目标之间相对旋转角速度。将(6)代入(5)中可得为便于推导，将(7)重写为其中式(11)中s(f,ta)和分别表示式(7)中求和项的幅度和相位。相位由多个散射点之间的差分距离共同决定。相应的，定义由于雷达与目标之间距离r0(ta)远大于目标散射点之间的差分距离，因此远小于将r0(ta)泰勒展开，可以得到相应的其中分别表示一阶导数、二阶导数。α1与α2由雷达与动目标之间的相对运动决定。当子孔径持续时间tsub较小时，上式中高阶项可以忽略。由式(13)可以看出，是关于ta的一个二次函数，再结合上文中远小于的结论，可以利用二次曲线拟合的方法，从ssk(f,ta)的相位中估计得到ssk(f,ta)的相位可以由下式得到其中表示估计值，∠(·)表示计算复信号幅角，unwrap(·)为解缠绕操作。现代数字信号处理中，信号一般会离散化。用表示距离频率f的离散采样点，式(14)离散化形式可以表示为相应的距离估计为：其中nr表示距离频率采样点数。利用二次曲线对进行拟合，即可得到r0(ta)的估计结果。r0(ta)的估计为分别表示常数项、一次项、二次项拟合系数。利用上式构造相位补偿因子将式(17)与式(8)相乘，得到运动补偿之后的回波信号为其中上式中，只有二次项相位会影响成像质量，利用map-drift算法，补偿该残余二次相位误差后，雷达回波为随后沿f(距离频率)傅里叶变换，沿ta(方位时间)做傅里叶变换，得到第k个子孔径下动目标重聚焦结果其中fd与慢时间ta互为傅里叶变换对，表示多普勒频率。按上述方法，逐个处理所有子孔径回波，即可获得多个子孔径观测下动目标重聚焦图像序列。与现有技术相比，本发明具有以下优点：1)本发明采用csar模式，能够获得多子孔径(多角度)动目标重聚焦成像结果；2)通过直接在距离频域提取相位，本发明能够在目标运动参数未知情况下，校正距离徙动，实现运动补偿。附图说明图1为本发明的应用场景示意图；图2是某一子孔径下观测几何示意图；图3是本发明方法流程图；图4是实施例中某一实测动目标回波在不同子孔径下重聚焦成像结果图。具体实施方式为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。图1为本发明的应用场景示意图，图中雷达平台沿圆周飞行，速度为波束始终指向探测区域，对探测区域内动目标(汽车)进行长时间观测。图2为子孔径下观测几何示意图，选取动目标某一散射点为参考点进行运动补偿。图3为本发明的流程图。一种适用于csar的地面动目标重聚焦方法，包括以下步骤：s1、接收动目标的雷达回波信号，进行预处理。预处理包括距离向脉冲压缩以及地杂波抑制。距离向脉冲压缩可以通过匹配滤波实现。根据动目标多普勒频偏是否大于静止杂波多普勒展宽可以将动目标分为快速运动目标和慢速运动目标。针对快速运动目标，地杂波抑制可以通过简单的多普勒滤波实现。针对慢速运动目标，利用多通道技术如dpca、stap等可以实现高性能杂波抑制。以快速运动目标为例，将回波变换到多普勒域，通过多普勒滤波，根据地杂波多普勒展宽确定截止频率，滤除低频带杂波，实现地杂波抑制，获取动目标回波；s2、子孔径回波划分：将接收到的雷达回波沿方位划分成若干个子孔径用于后续处理；s3.逐频点提取相位：将子孔径动目标回波变换到距离频域，逐频点提取信号沿方位向的相位；s4.二次曲线拟合：对提取出的相位进行二次曲线拟合，估计动目标参考点距离徙动曲线；s5.运动补偿：利用拟合获得到二次曲线，构造相位补偿因子，实现运动补偿；s6.md自聚焦：利用map-drift(md)自聚焦方法补偿相位估计误差，提高动目标重聚焦成像质量。s7.二维傅里叶变换：分别沿方位向做傅里叶变换(fft)和沿距离向做逆傅里叶变换(ifft)，获得子孔径下动目标重聚焦结果；s8.对所有子孔径重复s2-s7步骤，获得多子孔径下动目标重聚焦图像序列。本发明方法通过实测雷达回波数据进行了验证，实验结果证明了本发明的有效性。在实验中，本发明中的系统参数如下表1所示。表1发射信号中心频率(fc)ka波段发射信号带宽(b)900mhz雷达平台速度(v)80m/s圆周轨迹半径6000m雷达平台飞行高度3000m图4是某一非合作地面动目标(大型卡车)在不同子孔径下重聚焦结果示意图，子孔径持续时间为0.6s。从上到下，图(a)、(c)、(e)、(g)，分别对应子孔径1、子孔径20、子孔径40及子孔径60。左侧图(a)、(c)、(e)、(g)为按照静止目标成像方法对动目标进行成像的结果，右侧图(b)、(d)、(f)、(h)为使用本发明方法进行重聚焦成像结果。可以看出，使用本发明方法处理后，动目标聚焦结果有很大提升，动目标轮廓明显，聚焦质量高。以上所述仅是本发明的优选实施方式，本发明的保护范围并不局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。当前第1页12