一种小面阵星载tdi ccd相机的快速几何精校正方法
【技术领域】
[0001]本发明属于遥感对地观测技术领域，涉及一种小面阵星载TDI (XD相机的快速几 何精校正方法。
【背景技术】
[0002] 由于星上星敏感器测姿不准，新一代卫星通过一个对地观测的CCD相机，通过与地 面控制点影像匹配获取地面的精确指向数据。
[0003] 由于卫星姿态控制与测量技术的限制，CCD相机成像时的姿态测量结果存在一个 较大偏差，直接使用卫星自带系统测得的原始姿态数据对图像进行几何校正，在水平方向 上结果会存在较大的系统性偏差。而由原始测量姿态的偏差，在解算地理坐标时通过DEM获 得的像点高程也不完全精确，各像点的高程估算会出现错误，则校正后的图像会出现不规 则的变形，用图像匹配的方式得到的控制点，很可能达到几个像元的误差。
[0004] 提供准确的足印相机几何精校正结果对于立体测图卫星数据处理具有非常重要 的意义。一是在激光分解时，足印相机几何校正的不规则变形可能会给全波形的激光回波 组分带来干扰，影响波形分解的结果，甚至会导致分解后的波形与地物的关联错误。第二个 影响是激光指向依靠足印相机与参照的匹配，若使用不精确的几何校正结果进行匹配，很 可能会有几个像元的偏差甚至更大的误匹配结果，造成激光指向的计算结果不准确，最终 也影响的广义控制点的生成。
[0005] 立体测图卫星数据处理系统需要在较短的时间内处理百万级的激光数据，而每一 个激光光斑都需要一次几何精校正的结果，这就需要几何精校正的算法具有极高的效率才 能满足整个数据处理系统的速度需求。
【发明内容】
[0006] 本发明所要解决的问题是：提供一种小面阵星载TDI CCD相机的快速几何精校正 方法。
[0007] 本发明提供的技术方案采用如下步骤：
[0008] (1)相机像主点及临近像点坐标计算;具体方法步骤如下：
[0009] (1-1)获取(XD图像成像时卫星平台的瞬时姿态数据、(XD图像的像元畸变参数，获 得测区的DEM数据，对于像主点为中心的128X128小窗口图像，使用直接定位法得到小图像 每个像元的地面平面坐标点；
[0010] (1-2)将坐标点转为WGS84坐标系下；
[0011] (1-3)将计算得到的坐标值都加一个统一的平面地理偏移，第一次迭代时初值设 为〈〇,〇>;
[0012] (1-4)根据添加地理偏移后的平面坐标从DEM上获取高程，重新计算像点坐标； [0013] ⑵根据得到的坐标，使用最近邻法对像主点及临近像点进行重采样;使用双线性 插值法，根据每一个像点的WGS84坐标，对图像进行重采样，重采样后的分辨率与参照图像 的分辨率相同；
[0014] (3)将重采样图像与参考双线阵图像进行地理约束下的匹配;方法步骤如下：
[0015] (3-1)用校正重采样图像与参照图像的子图进行模板匹配，得到其在参照图像上 的位置；
[0016] (3-2)根据步骤（1)中（1-1)得到的像主点位置确定其在重采样图像上的像点位 置；
[0017] (3-3)根据步骤(2)得到的像主点在重采样图像上位置、根据(3-2)得到的重采样 图像在参照图像上的位置，计算得到像主点在参照图像上的位置；
[0018] (3-4)根据参照自身的几何定位结果，计算像主点在参照上的位置的地理坐标；
[0019] (4)计算重采样图像的平面误差偏移；用计算得到的像主点真实地理坐标减去像 主点匹配得到的地理坐标，得到像主点的平面误差偏移；
[0020] (5)算法迭代控制，判断新的平面偏移误差与上一次迭代得到的平面偏移是否相 同，若两次结果相同，或者差值小于1个像元，则认为算法趋于稳定，进入下一步处理；否则 进行步骤(2)处理;根据限差判断算法进行下一步或者跳回步骤(1)重新进行计算；
[0021] (6)使用得到的平面误差偏移对图像进行几何精校正，步骤如下：
[0022] (6-1)获取CCD图像成像时卫星平台的瞬时姿态数据、CCD图像的像元畸变参数，获 得测区的DEM数据，对于像主点为中心的128X128小窗口图像，使用直接定位法得到小图像 每个像元的地面平面坐标点；
[0023] (6-2)将坐标点转为WGS84坐标系下；
[0024] (6-3)将计算得到的坐标值都加一个统一的平面地理偏移，第一次迭代时初值设 为〈〇,〇>;
[0025] (6-4)根据添加地理偏移后的平面坐标从DEM上获取高程，重新计算像点坐标；
[0026] (6-5)使用双线性插值法，根据每一个像点的WGS84坐标，对图像进行重采样，得到 最终的几何精校正结果。
[0027]本发明具有良好的可扩展性，可满足所有类似的星载高分辨率TDI CCD相机精校 正的业务需求，具有快速、高精度的特点。
【附图说明】
[0028]图1是卫星TDI C⑶相机的快速几何精校正方法流程图。
[0029] 图2是像元畸变示意图，图（a)为姿态误差引起的CCD图像几何变形图(b)为像点位 置偏差示意图。
[0030] 图3是计算像点平面地理坐标的流程图。
[0031]图4是使用偏移距离对相机图像精确校正流程图。
【具体实施方式】
[0032] 为更好地阐述本发明的技术方案和优点，下面将结合附图对本发明的实施过程做 进一步的描述。
[0033] 本发明的理论基础是小视场角图像各像元水平偏移量近似相同。图2是像元水平 误差变化情况，图2(a)是由于星上姿态测量设备的误差导致的几何校正后图像的地理偏 移。图中显示了姿态不准会导致所有像元往同一个方向偏移。
[0034]图2(b)是CCD相机的任意两个像点的变形情况。其中0是相机中心，G是星下点，A1、 B1是地面上的两个点，由于姿态的误差，其校正后的位置在A2、B2处。
[0035]其中A1与A2之间的几何畸变情况可表示为：
[0036] =//^[Uin( ^(+Δ6* -tan 6^,]
[0037]由于相机视场角较小，即Θα ? 0,，ΛΘ ? 〇则：
[0038] La * H*tanA9
[0039] 则可认为点A与点B的几何变形长度相同：
[0040] La * H*tanA9 ~ Lb
[0041] 根据上述公式的推导，由于相机的视场角较小，各个像元经过几何校正后其地理 偏差都可认为近似相同。若能够找到像主点的平面偏移量，则可对图像所有的像点的坐标 进行改正，得到精确的校正结果。该方法不需要对相机的指向进行精确的检校，能够快速实 现相机的精校正。
[0042] 本发明的第一步涉及CCD图像的粗校正，图3为本发明使用共线方程，根据内方位 元素、卫星姿态、相机安置角、初始高程计算像点平面坐标的原理图。包含以下步骤：
[0043]获取CCD图像成像时卫星平台的瞬时姿态数据、CCD图像的像元畸变参数，获得测 区的DEM数据。
[0044]使用直接定位法得到像主点周围128X128的窗口内每个像元的地面平面坐标点。 首先根据焦距、像元尺寸、像元畸变检校值计算该像点的内方位元素；
[0045] 使用星下点的DEM的高程值作为直接法定位的初始值，使用共线方程、根据内方位 元素、卫星姿态、相机安置角、初始高程计算像元的平面坐标；
[004