本发明涉及一种几何定标方法，特别是一种星载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
其中，
利用步骤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
其中，
利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，根据式(7)，采用最小二乘法，解算距离向斜距改正参数δtr和方位向时间改正参数δta。
本发明由于采取以上技术方案，其具有以下优点：1、本发明是通过基准影像和待标定影像的配准来获取待标定影像几何定标参数解算所需的控制点，不需要依赖地面定标场的控制数据。2、采用本发明可以对星上成像参数进行定期监测，提高影像几何定位精度。3、采用本发明可以实现sar卫星的常态化、短周期几何定标，不需要在地面提前布设靶标，大大节约了人力和财力成本。
附图说明
图1为本发明流程图。
图2为同名点两次成像示意图。
具体实施方式
如图1所示，本发明星载sar影像几何交叉定标方法，具体包括以下步骤：
1)建立待标定影像大气延迟改正模型；电磁波信号在大气中传播遵循fermat定律，即信号在空间任意两点之间传播时间最小化。已知雷达观测信号在大气中传播速度v、折射率n与真空中的光速c的关系有n＝c/v，由此可推得在大气中雷达观测信号发射到地面z点的传播时延为：
大气延迟改正模型多写成大气天顶延迟和高度角相关映射函数的乘积：
其中，n(z)为天顶方向大气折射率，m(ε)是与入射角ε相关的映射函数。
采用以下步骤计算大气天顶延迟：
①利用美国国家环境预报中心(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是脉冲重复频率。
3)获取同名点对
通过影像匹配的方式或者人工刺点的方式从待标定影像和基准影像上选取3对以上同名点(x,y)和(x',y')，同名点两次成像的示意图如图2所示，θ1，θ2分别为待标定影像和基准影像的入射角，δr1,δr2分别为高程误差δh在待标定影像和基准影像上的投影，δs为待标定影像和基准影像的高程投影差之差。
4)大气延迟改正值计算
根据待标定影像成像时间，利用nationalcentersforenvironmentalprediction(ncep)提供的全球大气数据和欧洲定轨中心(code)提供的电离层电子含量分布数据，根据步骤1)建立的大气延迟改正模型计算待标定影像的大气延迟改正值tdelay。
5)定标参数解算
利用基准影像距离多普勒模型及srtm-dem数据通过迭代计算(x',y')对应的地面坐标(x,y,z)，则得到待标定卫星影像控制点(x,y,x,y,z)，几何定标模型可以写成如下形式：
对上式建立误差方程如下：
v＝bx-l
其中，
利用获取的若干对同名点和步骤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
其中，
利用步骤3)获取的若干对同名点和步骤4)得到的大气延迟改正值代入步骤2)的建立的几何定标模型，根据式(7)，采用最小二乘法，解算距离向斜距改正参数δtr和方位向时间改正参数δta。
各模块和各步骤的实现方式相对应，本发明不予撰述。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。