基于平行观测的星载推扫式光学传感器高频误差消除方法技术

本发明专利技术涉及一种基于平行观测的星载推扫式光学传感器高频误差探测方法，属于航空航天领域。该方法通过分析TDICCD器件使当前高分辨率对地光学卫星的主流成像器件，分析得出在利用全球SRTM‑DEM数据获取高程的前提下，平行观测中的同名交会误差可能是由成像几何参数中的高频误差和内方位元素误差造成，并提出内方位元素误差属于较为稳定的系统误差，在一段时间内不会发生变化或者变化不显著，可通过在轨检校消除。通过对平行观测中同名点的提取，计算同名点交会误差，实现高精度的在轨检校，精确恢复相机内方位元素后建立几何定位模型，在SRTM‑DEM辅助下计算同名点交会误差。对交会误差进行频谱分析，即可探测高频误差的频率、振幅等信息。

High Frequency Error Elimination Method for Spaceborne Push-Broom Optical Sensor Based on Parallel Observation

【技术实现步骤摘要】
基于平行观测的星载推扫式光学传感器高频误差消除方法
本专利技术涉及一种基于平行观测的星载推扫式光学传感器高频误差探测方法，属于航空航天领域。
技术介绍
在2010年启动实施的高分专项推动下，天绘系列卫星、资源三号卫星及高分系列卫星陆续发射，我国在航天遥感领域取得了系列突破和可喜成果。但时至今日，我国资源系列、遥感系列卫星的大量数据仍然因不具备定位误差的系统性而得不到较好应用，从而造成堆积浪费。保障国产卫星数据应用效果，是真正解决国内遥感数据自主化的关键。目前，国外研究的对象仅仅针对平台的高频抖动，在抖动的平滑性假设或是正余弦波形描述的条件下从配准误差中恢复高频误差，难以适用于我国现役在轨卫星的高频误差消除。国内现有研究中，没有充分考虑现役光学卫星平台特征，缺乏对高频误差源的完整识别，难以真正应用于现役国产卫星而提升其几何质量。研究高频误差探测及消除，对实现光学卫星几何质量的提升具有重要意义。
技术实现思路
本专利技术的目的是为了解决现有高频误差探测方法精度不高的问题，提供一种基于平行观测的星载推扫式光学传感器高频误差探测方法，该方法基于全球SRTM-DEM数据和几何检校技术剥离出高程误差本文档来自技高网...

【技术保护点】
1.基于平行观测的星载推扫式光学传感器高频误差探测方法，其特征在于：包括如下步骤：步骤1，CCD搭接区特征点配准；采用几何定位模型，计算相邻CCD线阵重叠像素，以重叠像素数作为同名点搜索的最大窗口，以尽量保障为每一行影像获取同名点，以便探测更高频率的误差；匹配过程基于相关系数测度确定像素级配准点位，进一步基于最小二乘匹配获取子像素级点位；给定CCD1线阵上一点(x1,y1)，其同名点，即子像素级点位，获取方法如下：步骤1.1、基于几何定位模型，计算CCD1线阵、CCD2线阵重叠像素数overlap，以该重叠像素数overlap为相关系数配准的最大搜索窗口；步骤1.2、基于几何定位模型，计算(x...

【技术特征摘要】
1.基于平行观测的星载推扫式光学传感器高频误差探测方法，其特征在于：包括如下步骤：步骤1，CCD搭接区特征点配准；采用几何定位模型，计算相邻CCD线阵重叠像素，以重叠像素数作为同名点搜索的最大窗口，以尽量保障为每一行影像获取同名点，以便探测更高频率的误差；匹配过程基于相关系数测度确定像素级配准点位，进一步基于最小二乘匹配获取子像素级点位；给定CCD1线阵上一点(x1,y1)，其同名点，即子像素级点位，获取方法如下：步骤1.1、基于几何定位模型，计算CCD1线阵、CCD2线阵重叠像素数overlap，以该重叠像素数overlap为相关系数配准的最大搜索窗口；步骤1.2、基于几何定位模型，计算(x1,y1)对应的地物坐标(X,Y,Z)，并计算(X,Y,Z)在CCD2线阵上对应的像素坐标步骤1.3、以为中心，以重叠像素数overlap为最大搜索窗口，按下式逐点计算与以(x1,y1)为中心的领域的相关系数；并取最大相关系数点位(x'2,y'2)作为像素级配准点位；式中：i、j分别为像素所在行数、列数；(c，r)为搜索区中心的坐标；w,h分别表示相关系数计算窗口的宽、高；g、g′分别表示CCD1线阵、CCD2线阵影像灰度；分别表示CCD1线阵、CCD2线阵影像灰度的平均值；为(x1+r，y1+c)点的灰度；为(x2+r，y2+c)点的灰度。步骤1.4以点位(x1,y1)、(x'2,y'2)作为初值，按式(2)进行最小二乘匹配，获取子像素级配准点位(x2,y2)，即完成配准；步骤2，同名点粗差探测，去除误匹配点对；由于各线阵CCD上存在辐射差异且部分区域缺乏纹理特征原因，误匹配几乎难以避免，从而对高频误差探测及消除造成较大影响；因此，需要对误匹配点对进行剔除；将同名点对相对位置关系按横坐标作为时间，纵坐标作为像素建立关系坐标系；对同名点原始相对关系进行中值滤波，对比同名点原始相对关系和该中值滤波结果，若同名点原始相对关系与中值滤波结果之差大于预设阈值，则判定为误匹配点，则剔除该点；步骤3，高频误差计算；通过在轨几何检校消除相机的内方位元素误差后，平行观测中的同名点交会误差即为高频误差的具体表现；定义平行观测同名点对分别成像于t0时刻和t1时刻，则同名点交会误差表示为：式中，(ΔXΔYΔZ)T为定位误差；为t0时刻的定位误差；为t1时刻的定位误差；显然，如果t0时刻和t1时刻星上误差不同，则导致定位误差不相等，即从而产生同名点交会误差；通过对该公式的分析能够看出，卫星的高频误差即为高频姿态误差，则基于偏置矩阵对姿态误差的补偿原理，在式(3)中引入姿态补偿矩阵：式中，(ΔXΔYΔZ)T为定位误差；为t时刻对应的姿态补偿矩阵，定义如下：其中，ωt、κt为待求解的姿态补偿角；现引入相对姿态概念，通过求解相对姿态消除高频姿态误差对几何定位模型的影响；令A1代表卫星成像过程中的真实姿态，A0代表星上下传的姿态；如果以前一时刻t0的定位模型为基准，仅通过解求后一时刻t1相对于t0的姿态补偿矩阵来消除同名点交会误差，则实际求解恢复的姿态为A2；通过对比可知，恢复姿态A2的过程能够恢复不同时刻姿态间的相对关系；当利用A0姿态进行几何处理时，姿态误差为随机误差；而当利用A2姿态进行几何处理时，姿态误差为系统误差，则对影像的相对定位精度而言，A2与真实姿态A1具有等效性；因此，通过恢复A2，不仅可以简化求解，也可以求解高频姿态误差；当不考虑时间因素时：根据对同名点基于几何定位模型及全球SRTM-DEM计算对应的地物坐标，得到控制点；当行上控制点数≥2，则可用常量偏置矩阵方法求解t1时刻对应的姿态补偿矩阵Roffset；当行上控制点数<2时，无法求解；当考虑时间因素时：以同名点交会误差的相近程度作为归类测度，对所有同名点进行分组，逐组求解姿态补偿矩阵，降低误匹配的影响；在CCD1、CCD2上获得N对同名点对(xccd1,yccd1,xccd2,yccd2)i,i≤N，并且yccd1＜yccd2，对所有同名点根据相近程度进行分组，并逐组求解姿态补偿矩阵的具体求...

【专利技术属性】
技术研发人员：孙向东，张过，蒋永华，
申请(专利权)人：张过，
类型：发明
国别省市：湖北,42