1.一种应用SAR影像提升国产光学卫星影像无控几何质量方法,其特征在于,包括以下步骤:
11)构建附加太阳高度角 变量的无控几何定位误差补偿模型;利用光学卫星的初始内外方位参数以及姿轨星历数据构建严密几何定位模型,每隔5000行计算卫星质心与太阳质心连线矢量在卫星本体坐标系Z轴的夹角,根据太阳高度角 构建6阶傅里叶级数形式的无控几何定位误差补偿模型;
12)应用SAR影像提升国产光学卫星影像无控几何质量准备工作;获取待进行误差补偿的长条带国产光学卫星影像,获取长条带影像覆盖区对应的SAR影像和SRTM DEM数据,并完成光学影像和SAR影像的匹配以得到求解误差模型系数的参考数据;
13)应用从SAR影像和SRTM DEM中提取的参考数据对无控几何定位误差补偿模型系数进行求解;对偶特定光学卫星长条带影像,每隔5000逐太阳高度角 构建无控几何定位误差补偿模型并线性化处理,对无控几何定位误差补偿模型中待求补偿系数赋初始值,然后迭代求解,直至计算完太阳高度角范围内的所有误差补偿模型系数。
2.根据权利要求1所述的应用SAR影像提升国产光学卫星影像无控几何质量方法,其特征在于,所述构建附加太阳高度角 变量的无控几何定位误差补偿模型包括以下步骤:
21)构建光学卫星影像的严密几何定位模型;利用影像的内外方位元素以及姿轨星历数据,构建严密几何定位模型如下:上式中,为观测像元所在的影像列数, 为像元 卫星本体坐标系下的指向角,是像元 的对应观测时间, 是卫星本体坐标系到星敏坐标系的旋转矩阵, 是星敏感器坐标系到J2000坐标系的旋转矩阵, 是J2000坐标系到WGS84坐标系的旋转矩阵,且 和 随着不同观测时间 而变换,是比例系数,为相机中心在WGS84坐标系下的坐标矢量,
分别为相机中心在WGS84坐标系下X轴、Y轴、Z轴的坐标, 为观测目标在WGS84坐标系下的坐标矢量,X、Y、Z为观测目标在WGS84坐标系下X轴、Y轴、Z轴的坐标;
22)计算观测时刻太阳与卫星本体质心连线在卫星本体坐标系下与 轴的夹角 ,本发明称该夹角为太阳高度角,具体步骤为:
221)由于太阳高度角在相邻影像行之间变化不大,本发明取每隔5000行计算;根据待计算太阳高度角的行获取成像时间 ,根据卫星轨道模型内插出卫星质心的位置矢量,其中 分别为卫星质心在 时刻位于WGS84坐标系下X轴、Y轴、Z轴的坐标;
222)获取太阳质心在WGS84坐标系下的位置矢量 ,其中 为太阳质心在 时刻位于WGS84坐标系下X轴、Y轴、Z轴的坐标;
223)获取 时刻卫星质心与太阳质心的相对归一化位置矢量 ,其表达式为,其中 符号为取模运算;
224)利用光学卫星影像的严密几何定位模型中的旋转矩阵 将相对归一化位置矢量 由WGS84坐标系转到卫星本体body坐标系下,转换完后的矢量记为 ,其具体表达式如下:;
225)计算 矢量与本体坐标系Z轴的夹角 ,计算公式如下:上式中 为本体坐标系Z轴的基向量 ;
23)构建附加太阳高度角 变量的无控几何定位误差补偿模型;
受星敏视场空间误差、星敏安装矩阵热变形、热弹性误差、翼板热变形等误差源影响,国产光学卫星在传统定标后姿态角仍存在低频误差,从而造成全球无控几何定位慢漂;在定位模型基础上,构建如下补偿模型:上式中, 为补偿外方位元素误差的补偿矩阵,具体形式如下:上式中,为太阳高度角, 、 和 分别表示 变量下的 俯仰姿态角误差补偿值, 翻滚姿态角误差补偿值和 偏航姿态角误差补偿值,具体表达式如下:上式中,K为傅里叶级数阶数,在本发明中K值设定为3, 与 分别表示余弦分量与正弦分量, 为变量 的变化频率, 根据特定光学卫星运行轨道范围内的太阳高度角变量变化最大值和最小值参数进行设定; 、 、 表示太阳高度角变量 下的姿态误差值; , , , , , 表示误差模型系数,为待求未知量。
3.根据权利要求1所述的应用SAR影像提升国产光学卫星影像无控几何质量方法,其特征在于,所述构建附加太阳高度角 变量的无控几何定位误差补偿模型包括以下步骤:
31)国产光学卫星长条带影像获取:由于本发明需对整个卫星运行轨道范围内的太阳高度角构建补偿模型,需获取长条带影像,影像尽可能云量少,且尽可能的覆盖国内外区域;
32)参考影像数据的获取:为了求解无控几何定位误差补偿模型补偿系数,获取长条带影像覆盖区对应的SAR影像和SRTM DEM数据,并对SAR影像进行正射校正和滤波处理;
33)国产光学卫星与SAR卫星影像进行匹配,获取SAR影像平面坐标,并以平面坐标在SRTM DEM影像上内插出高程值,形成三维参考点大地测量坐标 ,其中 为通过匹配从SAR影像上提取的纬度和经度, 为从SRTM DEM对应平面位置提取的高程,并将大地测量坐标 转为地心直角坐标形式 ,其中 分别为参考点在地心直角坐标系下的坐标分量。
4.根据权利要求1所述的应用SAR影像提升国产光学卫星影像无控几何质量方法,其特征在于,所述应用从SAR影像和SRTM DEM中提取的参考数据对无控几何定位误差补偿模型系数进行求解包括以下步骤:
41)转化误差模型形式:将23)构建的附加太阳高度角 变量的无控几何定位误差补偿模型转换为如下形式:上式中,下标 表示为长条带影像I , 为对应长条带影像I严密几何定位模型的X、Y、Z的观测值函数;
42)逐个太阳高度角 对无控几何定位误差补偿模型建立参数方程,以构建误差方程;
分别对无控几何定位误差补偿模型中待求的 , , , , , 误差模型系数,其中 ,求取偏导,则对于误差模型系数在长条带影像I上线性关系式记为:线性化:
,
线性化:
,
线性化:
,
上式中, 为 函数观测值残差向量, 为 函数观测值残差向量, 为 函数观测值残差向量, 、 、 、 、 以及 为 函数关于误差模型系数 , , , , 以及 的偏导数, 、 、 、、 以及 为 函数关于误差模型系数 , , , , 以及的偏导数, 、 、 、 、 以及 为 函数关于误差模型参数 , , , , 以及 的偏导数, 、 、 、 、以及 为误差模型系数 , , , , 以及 的改正数, 、以及 分别表示利用误差模型系数未知数初始值代入41)中公式计算得到的, 、 以及 分别表示观测值函数;
上式记为:
,
上式中,表示SAR影像误差补偿模型系数, 表示误差补偿模型系数相应的线性化后系数矩阵,表示为长条带影像I,为常数项,为观测值残差向量;
43)对待求误差模型系数赋初始值:分别对 , , , , , 误差模型系数,其中 附初值为0;
44)利用第二步提取的参考点坐标求解误差方程,得到误差模型系数的改正数 、、 、 、 、 并将改正值更新补偿误差模型系数中;
45)重复步骤43)至44)步,直到误差模型系数的的改正数 、 、 、 、、 都小于阈值 ,迭代停止;
46)得到的太阳高度角 变量下的误差模型系数,间隔5000行重复42)‑45),计算下一个太阳高度角 的误差模型系数,直至覆盖完光学卫星运行范围内的太阳高度角范围,完成无控几何定位误差补偿模型系数求解;
求解完误差补偿模型系数后,在标准景影像生产过程中将误差补偿模型补偿至几何定位模型中,可有效提升国产光学卫星影像无控几何定位精度。
5.根据权利要求1所述的一种应用SAR影像提升国产光学卫星影像无控几何质量系统,包括以下模块:光学几何定位模型构建模块,用于完成长条带光学卫星影像严密几何定位模型构建;
应用SAR影像提升国产光学卫星影像无控几何质量准备模块:完成SAR影像正射校正,光学影像与SAR影像匹配获取三维参考点功能;无控几何定位误差补偿模型构建模块:用于长条带影像逐太阳高度角 变量的无控几何定位误差补偿模型构建;无控几何定位误差补偿模型系数求解模块:用于完成误差补偿模型系数求解;几何定位误差补偿模型系数补偿模块:用于将几何定位误差补偿模型系数补偿至单景光学影像中,提升无控几何定位精度。