基于星载SAR方位向NLCS成像的目标定位方法与流程

本发明涉及sar成像技术，具体涉及一种基于星载sar方位向nlcs成像的目标定位方法。

背景技术：

合成孔径雷达(syntheticapertureradar，简称sar)是继光学成像雷达之后又一重要的遥感观测系统，借助脉冲压缩技术和合成孔径原理同时实现距离和方位的二维高分辨率成像。相对于光学遥感技术而言，微波遥感技术具有全天候、全天时、多波段、多极化、可穿透等优点。因此，sar被应用于军事侦察、环境保护、灾害控制、资源勘探、海洋观测、地质测绘、农业勘查等诸多方面，成为国际上应用最为广泛的航空航天对地观测途径之一。

cs算法是继rd算法之后又一个被广泛应用于星载sar图像目标定位的经典算法，其利用线性scaling函数对距离信号进行时域扰动，使距离线性调频信号的零频位置产生相对于参考距离呈线性关系的平移，抵消了距离徙动曲线沿距离向的线性空变性。然而，cs算法使用相同的多普勒调频率进行方位向匹配压缩，已经无法适应目前卫星分辨率、幅宽、合成孔径时间对图像处理精度的要求。方位向nlcs成像处理以传统为cs基础，在距离处理后的距离多普勒域增加一维方位时域扰动，配合时域线性距离徙动校正，对同一距离门内不同斜距目标的方位调频率进行一致化处理，消除了方位信号的方位空变特性，保证了频域方位压缩的高效性。但是，一维方位时域扰动在实现方位调频率一致化的同时，也会引起附加相位，使得多普勒中心频率平移，导致方位压缩位置发生变化。若仍然使用常规的方法进行几何校正，目标定位将出现偏差。

技术实现要素：

本发明的目的在于提出一种基于星载sar方位向nlcs成像的目标定位方法，以提高目标几何定位的精度。

实现本发明目的的技术解决方案为：一种基于星载sar方位向nlcs成像的目标定位方法，包括如下步骤：

步骤1：计算方位向nlcs成像后由扰动因子引起的方位向时间偏差；

步骤2：根据方位向时间偏差计算准确的方位时间轴；

步骤3：根据距离-多普勒法，结合准确的方位时间，进行目标几何定位。

本发明与现有技术相比，其显著优点为：分析了方位向nlcs处理后的方位向时间偏差，计算了准确的方位时间轴，并在sar严密几何定位模型中引入重新生成的方位时间，提高了目标几何定位的精度。

附图说明

图1为本发明基于星载sar方位向nlcs成像的目标定位方法的流程图。

具体实施方式

下面结合附图和具体实施例，进一步描述本发明方案。

基于星载sar方位向nlcs成像的目标定位方法，包括如下步骤：

步骤1：计算方位向nlcs成像后由扰动因子引起的方位向时间偏差；

1)采用局部拟合的方法求取调频率扰动函数；

对于同一距离门内的点，选择方位中心位置为参考方位位置，计算在该距离门内不同方位位置的点与参考点之间调频率的差：

δfr(ηi)＝fr(ηi)-fr(ηref)(1)

其中，ηi:第i的点的方位时间；

ηref:中心方位时间；

fr(ηi):第i的点的方位时间对应的调频率；

fr(ηref):中心方位时间对应的调频率；

得到δfr(ηi)和方位时间ηi的关系曲线后，采用多项式拟合的方法拟合曲线δfr(ηi)，得到其多项式表达式阶数n取大于等于3的整数。

将上述拟合多项式对方位时间进行两次积分，得到方位调频率扰动函数的相位为：

其中，η:方位时间；

n:拟合次数；

pn:n次项的系数；

位调频率扰动函数的相位；

π：圆周率；

则扰动函数表达式为由于扰动函数式φ(η)是η的n次多项式(n≥3)，所以，在其矫正调频率变化的同时会引入一些给成像质量带来的负面影响的相位误差，并且误差大小随η变化。

对于不同方位位置的目标，在方位时刻为ηd时，在其合成孔径时间ta内，由φ(η)引入的相位为：

由于扰动函数引入线性相位该相位会引起频谱偏移，偏移量为从而导致聚焦后的目标在方位向上产生如下移动：

其中，δy(ηd)：目标在方位向的位移；

v：雷达运动速度；

δf(ηd)：方位频谱偏移量；

扰动函数的线性项；

所以，方位向扰动处理后方位向时间偏差为：

步骤2：根据方位向时间偏差计算准确的方位时间轴；

其中，η：方位时间；

δη：方位时间偏移；

ηnlcs：考虑方位时间偏移的时间轴；

n为拟合多项式的阶数；

fr(ηref):中心方位时间对应的调频率。

步骤3：替换sar严密几何定位模型中方位时间轴，实现精确几何定位。

合成孔径雷达卫星定位原理是利用等距离线、等多普勒线在地球等高面上的交点确定影像像元位置。距离-多普勒法完全是从sar成像几何的角度来探讨像点与物点之间的对应关系，就是sar的严密几何模型。为计算简单，认为运行在高空中的卫星平台，运行比较稳定，其位置的变化用加速度、速度、位置状态变量就可以比较准确的描述。图像中某一时刻下卫星的位置和速度为：

其中，xs、ys、zs：卫星在x、y、z方向的位置；

x0、y0、z0：卫星在x、y、z方向的初始位置；

卫星在x、y、z方向的初始速度；

ax、ay、az：卫星在x、y、z方向的加速度；

vx、vy、vz：卫星在x、y、z方向的速度；

η：方位时间；

在常规的几何校正中，时间η是和方位向扫描行密切相关的，即η＝y/prf，y是图像的行数，prf为sar的脉冲重复频率，是一个固定的已知量。时间序列是一个线性变化的量。在方位扰动的nlcs成像处理后，方位时间轴被改变，准确的时间轴为：

将公式7中的η替换为ηnlcs，得到一个距离向的扫描行内，斜距图像上任意一点应该满足下列的斜距和多普勒条件：

式中：r0，近地点斜距；

λ，雷达波长；

mx，斜距的分辨率；

fd，多普勒中心频率；

(x,y):影像坐标；

(xp,yp,zp):与影像点对应的物点的地面三维坐标；

上述式子组成了完备的系统，r0、λ、mx、fd在成像处理后均是已知的参数。因此已知影像坐标(x,y)的情况下，求解与影像点对应的物点的地面三维坐标(xp,yp,zp)，也就完成了sar图像的精确几何校正。