一种星载SAR影像无场几何定标方法技术

本发明专利技术提供一种星载SAR影像无场几何定标方法，其特征在于：根据影像同名点定位一致性建立无场几何定标模型，所述无场几何定标模型由距离方程和多普勒方程组成，在距离方程上补偿斜距改正值和大气延迟改正值；获取同名点对，计算无场几何定标模型中参数，计算待标定影像的大气延迟改正值；通过约束同名点交会残差最小来探测得到几何定标参数。本发明专利技术能够解决SAR卫星的快速精确标定及短周期标定难题，不需要依赖地面定标场的控制数据和高精度参考数据，可以提高影像几何定位精度，不需要在地面提前布设靶标。

【技术实现步骤摘要】
一种星载SAR影像无场几何定标方法
本专利技术涉及一种几何定标方法，特别是一种星载SAR(SyntheticApertureRadar)影像无场几何定标方法。
技术介绍
几何定标是指利用地面高精度控制数据精确标定星上成像几何参数，对卫星影像的几何精度提升具有重要意义。传统的几何定标方法依赖于几何定标场的高精度控制数据，因此，卫星发射后，需要收集定标场区域的影像完成几何定标，只有当卫星成功获取到定标场区域影像才能进行几何定标，这无法满足卫星常态化、快速几何定标需求。另外，固定可用的几何定标场数量过少，卫星对定标场的拍摄频次过低，无法利用定标技术监测星上几何成像参数的变化，导致卫星产品精度降低。另一方面，新近出现的交叉几何定标方法虽然为常态化定标提供了一种解决思路，但是交叉几何定标的完成仍然需要高精度的参考数据，定标效率受到限制。
技术实现思路
因此，针对上述问题，本专利技术的目的在于提供一种不依赖地面控制数据和高精度参考数据的无场几何定标方法，根据影像同名点定位一致性建立定标模型，同时解算同名点地面坐标与几何定标参数，通过约束同名点交会残差最小来探测得到几何定标参数，解决SAR卫星的快速精确标定及短周期标定难题。为实现上述目的，本专利技术采取的技术方案为一种星载SAR影像无场几何定标方法，根据影像同名点定位一致性建立无场几何定标模型，所述无场几何定标模型由距离方程和多普勒方程组成，在距离方程上补偿斜距改正值Rs和大气延迟改正值Ratmo；获取同名点对，计算无场几何定标模型中参数，计算待标定影像的大气延迟改正值；通过约束同名点交会残差最小来探测得到几何定标参数。而且，所述无场几何定标模型如下式，其中，R为天线相位中心与未知地面点之间的距离，(Xs,Ys,Zs)为卫星位置矢量，(X,Y,Z)为未知地面点坐标，Rs为待探测的斜距系统误差，Ratmo为大气延迟误差，fD为目标点对应的多普勒中心频率，(Xv,Yv,Zv)为未知地面点成像时刻卫星的速度矢量，λ为雷达波长。而且，为解算几何定标参数，定标影像最少为3景，获取同名点对不少于3个。而且，计算无场几何定标模型中参数，包括根据根据定标影像的近距Rnear和像素间隔D计算天线相位中心与未知地面点之间的距离R，根据同名点对应的地面目标点成像时刻的时间ta，计算卫星的位置矢量Xs,Ys,Zs和速度矢量Xv,Yv,Zv。而且，所述计算待标定影像的大气延迟改正值，根据NCEP提供的对流层大气数据和CODE提供的电离层数据作为外部输入实现。而且，探测几何定标参数时，对无场几何定标模型的解算涉及双重迭代，第一层迭代以斜距改正值阈值作为迭代收敛的条件，第二层迭代以方位向时间改正值阈值作为迭代收敛的条件。本专利技术解决了SAR卫星的快速精确标定及短周期标定难题，其具有以下优点：1、本专利技术是通过约束同名点交会残差最小来探测得到几何定标参数，不需要依赖地面定标场的控制数据和高精度参考数据。2、采用本专利技术可以对星上成像参数进行定期监测，提高影像几何定位精度。3、采用本专利技术可以实现SAR卫星的常态化、短周期几何定标，不需要在地面提前布设靶标，大大节约了人力和财力成本。附图说明图1为本专利技术实施例的流程图。图2为本专利技术实施例的无场几何定标原理示意图具体实施方式以下结合附图和实施例详细说明本专利技术的技术方案。本专利技术提供一种星载SAR影像无场几何定标方法，根据影像同名点定位一致性建立几何定标模型，同时解算同名点地面坐标与几何定标参数，通过约束同名点交会残差最小来探测得到几何定标参数。如图1所示，本专利技术实施例提供的星载SAR影像无场几何定标方法，包括以下步骤：1)建立无场几何定标模型建立待标定影像的无场几何定标模型为：其中，R为天线相位中心与未知地面点之间的距离，(Xs,Ys,Zs)为卫星位置矢量，(X,Y,Z)为未知地面点坐标，Rs为待探测的斜距系统误差，Ratmo为大气延迟误差，fD为目标点对应的多普勒中心频率，(Xv,Yv,Zv)为未知地面点成像时刻卫星的速度矢量，λ为雷达波长。可见，上述几何定标模型由距离方程和多普勒方程组成，在距离方程上补偿斜距改正值Rs和大气延迟改正值Ratmo。2)获取同名点对如图2所示，S1、S2和S3是三张SLC影像对应SAR卫星天线相位中心，不规则曲线为真实的地表，点A为地表上的地物点。当不存在任何误差时候，S1与物方点A之间的斜距为R1，S2与物方点A之间的斜距为R2，S3与物方点A之间的斜距为R3，当不存在任何误差的时候，三张影像的交会点为A。实际上，由于存在影像系统误差(如斜距误差△R)，S2和S1的前方交会点为B，S3和S1的前方交会点为C，和S2的前方交会点为D。如果斜距测量的误差是由物方点坐标误差引起的，则前方交会点B、C和D应该是同一个点，三张SLC影像之间同名点的交会残差可以作为几何定标参数求解的一个准则。几何定标参数可以通过约束交会残差最小来探测得到。无场定标需要至少三幅定标影像，通过影像匹配的方式或者人工刺点的方式从三幅定标影像上选取3对以上同名点。具体实施时，同名点不能少于3个，定标影像最少为3景。3)计算无场几何定标模型中各参数①从影像提供的辅助文件中读取该定标影像的近距Rnear、像素间隔D、方位向起始时间的测量值ta0、脉冲重复频率prf、多普勒中心频率参数fD和雷达波长λ。②根据同名点在影像上的行列号坐标(i,j)，计算同名点对应的地面目标点成像时刻的时间ta和天线相位中心与未知地面点之间的距离R。计算公式如下式(2)，该计算实现为现有技术，本专利技术不予赘述。③根据同名点对应的地面目标点成像时刻的时间ta，计算卫星的位置矢量Xs,Ys,Zs和速度矢量Xv,Yv,Zv。该计算实现为现有技术，本专利技术不予赘述。4)大气延迟改正值计算建立大气延迟改正模型；电磁波信号在大气中传播遵循Fermat定律，即信号在空间任意两点之间传播时间最小化。已知雷达观测信号在大气中传播速度c1、折射率n与真空中的光速c的关系有n＝c/c1，由此可推得在大气中雷达观测信号发射到地面点的传播时延ttrop为：大气延迟改正模型为大气天顶延迟和入射角相关映射函数的乘积：其中，ΔL为待计算的大气延迟改正值，z为传输路径，n(z)为天顶方向大气折射率，m(ε)是与入射角ε相关的映射函数。实施例中，大气延迟改正值计算采用美国气象环境预报中心(NCEP)提供的对流层大气数据和欧洲定轨中心(CODE)提供的电离层数据作为外部输入。采用以下步骤计算大气天顶延迟：①利用美国国家环境预报中心(NCEP)提供的全球大气数据，从中获取影像区域的干大气压强Pd，地面温度T，经验常数k1确定干大气延迟分量Δdry：②根据NCEP大气数据，从中获取影像区域的湿大气压强Pw和地面温度T，经验常数k2，经验常数k3确定湿大气延迟分量Δwet：③根据影像辅助文件中提供的雷达信号频率f、欧洲定轨中心(CODE)提供的天顶方向电子含量tec和经验常数k4计算电离层天顶方向延迟分量Δiono：④确定映射函数m(ε)形式为：⑤将步骤①确定的干大气延迟分量Δdry、步骤②确定的湿大气延迟分量Δwet、步骤③确定的电离层天顶方向延迟分量Δiono、步骤④确定的映射函数m(ε)代入大气延迟改正模型，得到大气延迟改正模型为本文档来自技高网...

【技术保护点】
1.一种星载SAR影像无场几何定标方法，其特征在于：根据影像同名点定位一致性建立无场几何定标模型，所述无场几何定标模型由距离方程和多普勒方程组成，在距离方程上补偿斜距改正值Rs和大气延迟改正值Ratmo；获取同名点对，计算无场几何定标模型中参数，计算待标定影像的大气延迟改正值；通过约束同名点交会残差最小来探测得到几何定标参数。

【技术特征摘要】
1.一种星载SAR影像无场几何定标方法，其特征在于：根据影像同名点定位一致性建立无场几何定标模型，所述无场几何定标模型由距离方程和多普勒方程组成，在距离方程上补偿斜距改正值Rs和大气延迟改正值Ratmo；获取同名点对，计算无场几何定标模型中参数，计算待标定影像的大气延迟改正值；通过约束同名点交会残差最小来探测得到几何定标参数。2.如权利要求书1所述星载SAR影像无场几何定标方法，其特征在于：所述无场几何定标模型如下式，其中，R为天线相位中心与未知地面点之间的距离，(Xs,Ys,Zs)为卫星位置矢量，(X,Y,Z)为未知地面点坐标，Rs为待探测的斜距系统误差，Ratmo为大气延迟误差，fD为目标点对应的多普勒中心频率，(Xv,Yv,Zv)为未知地面点成像时刻卫星的速度矢量，λ为雷达波长。3.如权利要求书2所述星载SAR影像无场几何定标方法，其特征在于：为解算几何定...

【专利技术属性】
技术研发人员：张过，邓明军，赵瑞山，郭风成，徐凯，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北,42