45°旋转扫描方式空间相机的图像几何校正方法
【技术领域】
[0001] 本发明涉及光学遥感应用领域，具体涉及45°旋转扫描方式空间相机的图像几何 校正方法。
【背景技术】
[0002] 扫描成像是空间相机最常见的一种成像模式，而45°旋转扫描方式由于其扫描镜 尺寸小、幅宽大、稳定性好、能观测冷空间等诸多优点，被广泛用于空间相机中，典型的如中 分辨率光谱成像仪、扫描辐射计等。通过45°旋转扫描反射镜结合多元并扫的方式成像，而 45°扫描镜带来的像旋转，可以由K镜消旋系统或软件消旋方法予以消除。然而，在45° 旋转扫描反射镜推扫过程中，其在垂直于星下点轨迹方向的大幅宽会带来图像中的几何畸 变，表现为地面同一目标可能在图像中对应不同的像素点，对图像如实反映地物分布造成 了严重影响，故而需要对畸变进行校正。
[0003] 几何畸变校正的目的是要建立正确的畸变图像像素位置与校正图像像素位置之 间的对应关系。根据所需要先验信息的情况，目前的校正方法可以分为两类：严格几何校 正和近似几何校正。严格几何校正需要卫星轨道参数、拍摄时刻的卫星姿态参数、传感器参 数、地表参数、环境参数等。当这些参数无法获得的时候，只能采用近似几何校正方法，回 避真实的成像过程，而直接对产生的几何畸变进行模拟，常用的有多项式模型、有理函数模 型、局部区域模型等。严格几何校正中的很多参数实际很难获得，而近似几何校正的精度和 稳定性受模型对畸变的描述能力限制。针对45°旋转扫描方式空间相机的图像几何畸变， 目前尚未有现成的校正方法可以对其处理。
【发明内容】
[0004] 针对现有技术中的缺陷，本发明的目的是提供一种45°旋转扫描方式空间相机的 图像几何校正方法。
[0005] 本发明通过对45°旋转扫描方式的空间相机原理及成像过程进行分析，在遥感图 像中的位置与其对应的地面位置之间，建立了合理的位置对应模型，进一步构建了原始图 像与校正图像的坐标转换关系。
[0006] 根据本发明提供的一种45°旋转扫描方式空间相机的图像几何校正方法，包括如 下步骤：
[0007] 步骤1 :建立原始图像与校正图像的坐标转换关系：
[0009] 其中：（i，j)、（i'，j'）表示同一目标点分别在原始图像、校正图像中的行列号坐 标；R表示地球半径；H表示轨道高度；L。表示星下点像元分辨率；N表示原始图像列数;A表 示观测角度范围；原始图像为遥感图像；
[0010] 步骤2 :通过所述坐标转换关系，对原始图像进行空间坐标投影；
[0011] 步骤3 :对投影后得到的离散点，利用Delaunay三角剖分算法进行重采样；首先构 造 Delaunay三角网格，然后对于采样网格上的每一点，利用该点所在的三角形顶点的数据 进行自然邻域插值；从而获得原始图像经畸变校正后的数字校正图像。
[0012] 优选地，多个波段之间共用所述坐标转换关系、空间坐标投影关系以及Delaunay 三角网格。
[0013] 对于共用主光学系统的对地遥感设备，具有相同地面像元间隔的不同波段，位置 对应模型是完全相同的；地面像元间隔不同的波段之间，位置对应模型也是相同的，仅存在 参数上的差异。所以，多个波段之间可以共用坐标投影关系以及Delaunay三角网格，从而 可以简便地对所有波段的数据进行校正。
[0014] 与现有技术相比，本发明具有如下的有益效果：
[0015] 1、本发明图像几何校正方法，通过对45°旋转扫描方式的空间相机原理及成像过 程进行分析，建立了合理的坐标转换关系。
[0016] 2、本发明通过空间坐标投影以及基于Delaunay三角剖分算法的重采样，得到了 校正图像。
[0017] 3、本发明通过坐标转换关系的共用，可以便捷地对不同波段实施校正处理。因此， 本发明方法合理、实施简易，能够对图像几何畸变的有明显改善效果，应用前景广泛。
【附图说明】
[0018] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、 目的和优点将会变得更明显：
[0019] 图1为45°旋转扫描方式的空间相机成像过程示意图；
[0020] 图2为45°旋转扫描方式的空间相机图像几何畸变示意图。
[0021] 图3A-图3F分别为45°旋转扫描反射镜转过三个不同角度下的原始图像与校正 图像。
[0022] 图中：
[0023] 1-45°旋转扫描反射镜
[0024] 2-望远镜系统
[0025] 3-焦平面
[0026] 4-地面像元
【具体实施方式】
[0027] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术 人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术 人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明 的保护范围。
[0028] 本发明涉及遥感图像处理领域。本发明的目的在于提供一种适用于45°旋转扫 描方式下的空间相机图像几何校正方法，解决图像失配问题。本发明主要的改进点在于：通 过对45°旋转扫描方式的空间相机原理及成像过程进行分析，在遥感图像中的位置与其对 应的地面位置之间，建立了合理的位置对应模型，并构建了原始图像与校正图像的坐标转 换关系。依据坐标转换关系对原始遥感图像进行空间坐标投影。对投影后的离散点，利用 Delaunay三角剖分算法进行重采样。对于共用主光学系统的对地遥感设备，多个波段之间 可以共用坐标投影关系以及Delaunay三角网格，从而可以简便地对所有波段的数据进行 校正。本发明方法合理、实施简易，对图像畸变的改善效果明显。
[0029] 对本发明的理论分析基础简介如下：
[0030] 如附图1所示，通过45°旋转扫描反射镜1的旋转推扫，将地面较宽范围内，结合 卫星的相对运动获得二维图像。在45°旋转扫描反射镜旋转一圈的时间内，卫星星下点沿 轨运动多元探测器对应的地面距离，保证相邻两转内捕获的数据相接。
[0031] 将地球作为半径为R的球体，当在45°旋转扫描反射镜距离天底旋过α弧度时， 采样位置偏离星下点的地面距离LxS
[0033] 其中H为轨道高度，R为地球半径。
[0034] 若多元并扫探测器的间隔为d，光学系统焦距为f，那么物距L远处，在沿轨方向地 面像元间隔Ly可以由如下公式确定：
[0036] 依据成像几何关系，在45°旋转扫描反射镜距离天底旋过α弧度时的物距L满足
[0038] 公式（3)中，当α - 0时，L - Η，即星下点的地面像元间隔为
[0040] 若扫描观测角为土Α，当0〈 α〈Α时，Ly>L。，这是造成中分辨率光谱成像仪的遥感图 像存在畸变的最主要原因。附图2所示即为几何畸变示意图，遥感图像中的各像素坐标实 际在地球分布并不均匀。
[0041] 在分析了成像中的几何关系之后，利用其对畸变图像进行校正，该校正方法包 括：
[0042] 以垂轨方向为X方向，沿轨方向为y方向，对每转推扫获得的图像带进行分析，以 图像中心点为坐标原点，畸变校正的目的就是要建立畸变图像像素位置（i，j)与校正图像 像素位置之间的对应关系。对于多元探测器，记探测器元数为P，对于尺寸为PXN 的中分辨率光谱成像仪1级图像，若扫描镜电机匀速转动，那么位置（i，j)旋过的角度a i 为
[0046] 每转推扫中对星下点轨迹的位置即为（0，j),根据公式 (2)~公式（5)可以得出，
[0048] 依据对应关系公式（6)与公式（7)，可以将原始位置坐标投影到校正图像中，其中 涉及的固有参数包含R、H、L。、A、N，其中地球半径R~637lkm，N可以从图像产品知悉，而参 数H、L。、A均为中分辨率光谱成像仪的常见指标参数，容易获取，故而该方法在遥感图像应 用领域是可行的。
[0049] 根据位置对应关系进行投影之后，位置坐标点并非均匀采样网格位置，故而需要 进行重采样。由于投影之后的坐标点非均匀，故而一般插值方法不能适用，需要进行散点插 值方法。常用的散点插值方法为Delaunay三角剖分算法，首先构造 Delaunay三角网格，然 后对于均匀采样网格上的每一点，利用其所在的三角形顶点的数据进行自然邻域插值。
[0050] 对于相同分辨率的不同波段，其位置对应关系是相同的，由散点坐标构造出的 Delaunay三角网格也是相同的，只是三角顶点处的数据不相同而已。所以，具有相同分辨率 的不同波段可以共用三角网格简化计算。
[0051] 附图3为本发明对某型卫星中分辨率光谱成像仪遥感图像实施校正前后的对比 图。其中图3A、图3C和图3E分为为α在0.52、0.70、0.96(30°、40°、55°)时的原始图 像，图3B、图3D和图3F分别为对应的校正图像。由图3A-图3F可以看出，当45°旋转扫 描反射镜旋过的角度a由0变大时，原始图像中存在的畸变也越来越严重，而本发明方法 很好地解决了图像中由于畸变导致的失配。
[0052] 以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述 特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影 响本发明的实质内容。
【主权项】
1. 一种45°旋转扫描方式空间相机的图像几何校正方法，其特征在于，包括如下步 骤： 步骤1 :建立原始图像与校正图像的坐标转换关系：其中a'，j'）表示同一目标点分别在原始图像、校正图像中的行列号坐标；R表示地球半径；H表示轨道高度；L。表示星下点像元分辨率；N表示原始图像列数;A表示观 测角度范围；原始图像为遥感图像； 步骤2 :通过所述坐标转换关系，对原始图像进行空间坐标投影； 步骤3 :对投影后得到的离散点，利用Delaunay三角剖分算法进行重采样；首先构造Delaunay三角网格，然后对于采样网格上的每一点，利用该点所在的三角形顶点的数据进 行自然邻域插值；从而获得原始图像经畸变校正后的数字校正图像。2. 根据权利要求1所述的45°旋转扫描方式空间相机的图像几何校正方法，其特征在 于，多个波段之间共用所述坐标转换关系、空间坐标投影关系以及Delaunay三角网格。
【专利摘要】本发明提供了一种45°旋转扫描方式空间相机的图像几何校正方法，其通过对45°旋转扫描方式的空间相机原理及成像过程进行分析，在遥感图像中的位置与其对应的地面位置之间，建立了合理的位置对应模型，并构建了原始图像与校正图像的坐标转换关系。依据坐标转换关系对原始遥感图像进行空间坐标投影。对投影后的离散点，利用Delaunay三角剖分算法进行重采样。对于共用主光学系统的对地遥感设备，多个波段之间可以共用坐标投影关系以及Delaunay三角网格，从而可以简便地对所有波段的数据进行校正。本发明方法合理、实施简易，对图像畸变的改善效果明显。
【IPC分类】G06T5/00
【公开号】CN105046667
【申请号】CN201510465971
【发明人】郭玲玲, 赵其昌, 马文佳, 杨勇, 吴泽鹏, 薛伶玲, 李叶飞
【申请人】上海卫星工程研究所
【公开日】2015年11月11日
【申请日】2015年7月31日