一种星载SAR影像几何交叉定标方法和系统技术方案

本发明专利技术涉及一种星载SAR影像几何交叉定标方法和系统，包括步骤1)，建立待标定影像的大气延迟改正模型；步骤2)，建立待标定影像的几何定标模型；步骤3)，从待标定影像和基准影像上选取若干对同名点；步骤4，根据待标定影像成像时间，利用NCEP提供的全球大气数据和欧洲定轨中心(CODE)提供的电离层电子含量分布数据，通过步骤1)建立的大气延迟改正模型计算待标定影像的大气延迟改正值；步骤5)，利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，完成几何定标参数解算。采用本发明专利技术可以实现SAR卫星的常态化、短周期几何定标，不需要在地面提前布设靶标，大大节约了人力和财力成本。

A geometric cross calibration method and system for Spaceborne SAR images

The invention relates to a spaceborne SAR image geometric cross-calibration method and system, including step 1, establishing an atmospheric delay correction model for the image to be calibrated, step 2, establishing a geometric calibration model for the image to be calibrated, step 3, selecting several pairs of identical points from the image to be calibrated and the reference image, step 4, according to the image to be calibrated; For imaging time, the atmospheric delay correction model is established to calculate the atmospheric delay correction value of the image to be calibrated by using the global atmospheric data provided by NCEP and the ionospheric electron content distribution data provided by the European Orbital Determination Center (CODE) through step 1; step 5; and several pairs of identical points and step 4 obtained by step 3). The geometric calibration model established by gas delay correction value in step 2) completes the calculation of geometric calibration parameters. The invention can realize the normalization of the SAR satellite and the short period geometric calibration, and does not need to lay the target in advance on the ground, thus greatly saving the human and financial costs.

【技术实现步骤摘要】
一种星载SAR影像几何交叉定标方法和系统
本专利技术涉及一种几何定标方法，特别是一种星载SAR(SyntheticApertureRadar)影像几何交叉定标方法和系统。
技术介绍
几何定标是指利用地面高精度控制数据精确标定星上成像几何参数，对卫星影像的几何精度提升具有重要意义。传统的几何定标方法依赖于几何定标场的高精度控制数据，因此，卫星发射后，需要收集定标场区域的影像完成几何定标，只有当卫星成功获取到定标场区域影像才能进行几何定标，这无法满足卫星常态化、快速几何定标需求。另外，固定可用的几何定标场数量过少，卫星对定标场的拍摄频次过低，无法利用定标技术监测星上几何成像参数的变化，导致卫星产品精度降低。
技术实现思路
因此，针对上述问题，本专利技术的目的在于提供一种不依赖地面控制数据的几何交叉定标方法，在相近入射角情况下，基于同名点定位一致性约束来实现几何交叉定标，解决SAR卫星的快速精确标定及短周期标定难题。为实现上述目的，本专利技术采取以下技术方案：一种星载SAR影像几何交叉定标方法，其具体包括以下步骤：步骤1)，建立待标定影像的大气延迟改正模型；步骤2)，建立待标定影像的几何定标模型；步骤3)，从待标定影像和基准影像上选取若干对同名点；步骤4，根据待标定影像成像时间，利用NCEP提供的全球大气数据和欧洲定轨中心(CODE)提供的电离层电子含量分布数据，通过步骤1)建立的大气延迟改正模型计算待标定影像的大气延迟改正值；步骤5)，利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，完成几何定标参数解算。进一步的，所述步骤1)中，建立大气延迟改正模型的具体实现方式如下；①利用美国国家环境预报中心(NCEP)提供的全球大气数据，从中获取影像区域的干大气压强Pd，地面温度T，经验常数k1，确定干大气延迟分量：②根据NCEP大气数据，从中获取影像区域的湿大气压强Pw和地面温度T，经验常数k2，经验常数k3，确定湿大气延迟分量：③根据影像辅助文件中提供的雷达信号频率f、欧洲定轨中心(CODE)提供的天顶方向电子含量tec和经验常数k4计算电离层天顶方向延迟分量：④确定映射函数形式为：其中，m(ε)是与入射角ε相关的映射函数；⑤将步骤①确定的干大气延迟分量Δdry、步骤②确定的湿大气延迟分量Δwet、步骤③确定的电离层延迟改正、步骤④确定的映射函数代入大气延迟改正模型，得到大气延迟改正模型为：进一步的，所述步骤2)中，建立待标定影像的几何定标模型的具体实现方式如下；①根据距离-多普勒定位模型的反算方法确定目标点成像时的距离向快时间tr和方位向慢时间ta；②从影像提供的辅助文件中确定该定标景影像的距离向起始时间的测量值tr0和方位向起始时间的测量值ta0；③通过影像匹配的方式或者人工刺点的方式确定目标点的像平面列坐标x和像平面行坐标y，建立待标定影像的几何定标模型为：式中，tdelay为大气延迟改正值，Δtr是距离向斜距改正参数，Δta是方位向时间改正参数，fs是雷达采样频率，prf是脉冲重复频率。进一步的，所述步骤4)中，通过影像匹配的方式或者人工刺点的方式从待标定影像和基准影像上选取同名点(x,y)和(x',y')，同名点的对数大于3。进一步的，所述步骤5)中解算几何定标参数的具体实现方式如下，利用基准影像距离多普勒模型及SRTM-DEM数据通过迭代计算(x',y')对应的地面坐标(X,Y,Z)，则得到待标定卫星影像目标点(x,y,X,Y,Z)，其中(x,y)和(x',y')分别为待标定影像和基准影像上的同名点，几何定标模型可以写成如下形式：对上式建立误差方程如下：V＝BX-L其中，X＝[dΔtr,dΔta]T,利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，根据式(7)，采用最小二乘法，解算距离向斜距改正参数Δtr和方位向时间改正参数Δta。本专利技术还提供一种星载SAR影像几何交叉定标系统，包括如下模块：大气延迟改正模型建立模块，用于建立待标定影像的大气延迟改正模型；几何定标模型建立模块，用于建立待标定影像的几何定标模型；同名点选取模块，用于从待标定影像和基准影像上选取若干对同名点；大气延迟改正值计算模块，用于根据待标定影像成像时间，利用NCEP提供的全球大气数据和欧洲定轨中心(CODE)提供的电离层电子含量分布数据，通过步骤1)建立的大气延迟改正模型计算待标定影像的大气延迟改正值；几何定标参数解算模块，用于利用同名点选取模块中获取的若干对同名点和大气延迟改正值计算模块中得到的大气延迟改正值代入几何定标模型建立模块中建立的几何定标模型，完成几何定标参数解算。进一步的，所述大气延迟改正模型建立模块中，建立大气延迟改正模型的具体实现方式如下；①利用美国国家环境预报中心(NCEP)提供的全球大气数据，从中获取影像区域的干大气压强Pd，地面温度T，经验常数k1确定干大气延迟分量：②根据NCEP大气数据，从中获取影像区域的湿大气压强Pw和地面温度T，经验常数k2，经验常数k3确定湿大气延迟分量：③根据影像辅助文件中提供的雷达信号频率f、欧洲定轨中心(CODE)提供的天顶方向电子含量tec和经验常数k4计算电离层天顶方向延迟分量：④确定映射函数形式为：其中，m(ε)是与入射角ε相关的映射函数；⑤将步骤①确定的干大气延迟分量Δdry、步骤②确定的湿大气延迟分量Δwet、步骤③确定的电离层延迟改正、步骤④确定的映射函数代入大气延迟改正模型，得到大气延迟改正模型为：进一步的，所述几何定标模型建立模块中，建立待标定影像的几何定标模型的具体实现方式如下；①根据距离-多普勒定位模型的反算方法确定目标点成像时的距离向快时间tr和方位向慢时间ta；②从影像提供的辅助文件中确定该定标景影像的距离向起始时间的测量值tr0和方位向起始时间的测量值ta0；③影像匹配的方式或者人工刺点的方式确定目标点的像平面列坐标x和像平面行坐标y，建立待标定影像的几何定标模型为：式中，tdelay为大气延迟改正值，Δtr是距离向斜距改正参数，Δta是方位向时间改正参数，fs是雷达采样频率，prf是脉冲重复频率。进一步的，所述同名点选取模块中，通过影像匹配的方式或者人工刺点的方式从待标定影像和基准影像上选取同名点(x,y)和(x',y')，同名点的对数大于3。进一步的，所述几何定标参数解算模块中解算几何定标参数的具体实现方式如下，利用基准影像距离多普勒模型及SRTM-DEM数据通过迭代计算(x',y')对应的地面坐标(X,Y,Z)，则得到待标定卫星影像目标点(x,y,X,Y,Z)，其中(x,y)和(x',y')分别为待标定影像和基准影像上的同名点，几何定标模型可以写成如下形式：对上式建立误差方程如下：V＝BX-L其中，X＝[dΔtr,dΔta]T,利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，根据式(7)，采用最小二乘法，解算距离向斜距改正参数Δtr和方位向时间改正参数Δta。本专利技术由于采取以上技术方案，其具有以下优点：1、本专利技术是通过基准影像和待标定影像的配准来获取待标定影像几何定标参数解算所需的控制点，不需要依赖地面定标场的控本文档来自技高网...

【技术保护点】
1.一种星载SAR影像几何交叉定标方法，其特征在于，包括如下步骤：步骤1)，建立待标定影像的大气延迟改正模型；步骤2)，建立待标定影像的几何定标模型；步骤3)，从待标定影像和基准影像上选取若干对同名点；步骤4，根据待标定影像成像时间，利用NCEP提供的全球大气数据和欧洲定轨中心(CODE)提供的电离层电子含量分布数据，通过步骤1)建立的大气延迟改正模型计算待标定影像的大气延迟改正值；步骤5)，利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，完成几何定标参数解算。

【技术特征摘要】
1.一种星载SAR影像几何交叉定标方法，其特征在于，包括如下步骤：步骤1)，建立待标定影像的大气延迟改正模型；步骤2)，建立待标定影像的几何定标模型；步骤3)，从待标定影像和基准影像上选取若干对同名点；步骤4，根据待标定影像成像时间，利用NCEP提供的全球大气数据和欧洲定轨中心(CODE)提供的电离层电子含量分布数据，通过步骤1)建立的大气延迟改正模型计算待标定影像的大气延迟改正值；步骤5)，利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，完成几何定标参数解算。2.如权利要求1所述的一种星载SAR影像几何交叉定标方法，其特征在于：所述步骤1)中，建立大气延迟改正模型的具体实现方式如下；①利用美国国家环境预报中心(NCEP)提供的全球大气数据，从中获取影像区域的干大气压强Pd，地面温度T，经验常数k1，确定干大气延迟分量：②根据NCEP大气数据，从中获取影像区域的湿大气压强Pw和地面温度T，经验常数k2，经验常数k3，确定湿大气延迟分量：③根据影像辅助文件中提供的雷达信号频率f、欧洲定轨中心(CODE)提供的天顶方向电子含量tec和经验常数k4计算电离层天顶方向延迟分量：④确定映射函数形式为：其中，m(ε)是与入射角ε相关的映射函数；⑤将步骤①确定的干大气延迟分量Δdry、步骤②确定的湿大气延迟分量Δwet、步骤③确定的电离层延迟改正、步骤④确定的映射函数代入大气延迟改正模型，得到大气延迟改正模型为：3.如权利要求书1所述的一种星载SAR影像几何交叉定标方法，其特征在于：所述步骤2)中，建立待标定影像的几何定标模型的具体实现方式如下；①根据距离-多普勒定位模型的反算方法确定目标点成像时的距离向快时间tr和方位向慢时间ta；②从影像提供的辅助文件中确定该定标景影像的距离向起始时间的测量值tr0和方位向起始时间的测量值ta0；③通过影像匹配的方式或者人工刺点的方式确定目标点的像平面列坐标x和像平面行坐标y，建立待标定影像的几何定标模型为：式中，tdelay为大气延迟改正值，Δtr是距离向斜距改正参数，Δta是方位向时间改正参数，fs是雷达采样频率，prf是脉冲重复频率。4.如权利要求书1所述的一种星载SAR影像几何交叉定标方法，其特征在于：所述步骤4)中，通过影像匹配的方式或者人工刺点的方式从待标定影像和基准影像上选取同名点(x,y)和(x',y')，同名点的对数大于3。5.如权利要求书3所述的一种星载SAR影像几何交叉定标方法，其特征在于：所述步骤5)中解算几何定标参数的具体实现方式如下，利用基准影像距离多普勒模型及SRTM-DEM数据通过迭代计算(x',y')对应的地面坐标(X,Y,Z)，则得到待标定卫星影像目标点(x,y,X,Y,Z)，其中(x,y)和(x',y')分别为待标定影像和基准影像上的同名点，几何定标模型可以写成如下形式：对上式建立误差方程如下：V＝BX-L其中，X＝[dΔtr,dΔta]T,利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，根据式(7)，采用最小二乘法，解算距离向斜距改正参数Δtr和方位向时间改正参数Δta。6.一种星载SAR影像几何交叉定标系统，其特征在于，包括如下模块：大气延迟改正...

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