[发明专利]一种多时次遥感影像优化合成方法及系统

说明书

本发明公开了一种多时次遥感影像优化合成方法及系统，涉及遥感影像处理领域。该方法包括：获取光谱波段数据、太阳天顶角和观测天顶角；对应生成矩阵以及三个赋值为0的权重矩阵γ₁、γ₂以及γ₃；根据预设规则对三个矩阵进行赋值；将三个矩阵叠加计算，生成最终权重矩阵γ，对所述最终权重矩阵γ中每个像元按波段层取最大值，并根据选取结果生成指数矩阵λ；依次读取n个光学遥感影像文件中同一波段的光谱波段数据，生成数值矩阵η，根据所述指数矩阵λ，提取所述数值矩阵η中的数据，将所述数据作为该波段的输出数据；重复上述过程，直至所有波段均完成上述步骤，生成完整输出文件。本发明能够达到提取出最优质影像的效果。

技术领域

本发明涉及遥感影像处理领域，尤其涉及一种多时次遥感影像优化合成方法及系统。

背景技术

极轨卫星依倾斜轨道对地球进行光学成像，获得的影像在几何纠正以后，在低纬度地区，前后两个轨道影像间的重叠区域比较小，有时甚至还有缝隙，而在地球极区之外的高纬度地区，影像间的重叠区域很大，一日内有时可获得多达5幅影像，只有把这些影像合成为一幅高质量的影像，才能保持与低纬度地区的影像一致，生成一幅完整的全球影像。

具有地理坐标的同一地区多幅影像可以通过图像拼接的方式进行合成，最终的影像到底是取自哪一幅，是否取自质量最好的影像，没有办法控制。为了去除云的影响，通常采用植被指数最大值法将10天的影像合成为一幅影像，但是高纬度地区遥感影像的光照条件变化比较大、观测几何带来的畸变也比较严重，植被指数最大值法也不是把质量最好的影像取出来。且现有技术方案是适当约束卫星观测角，把卫星观测角较小的影像数据取出来。但并没有考虑太阳的光照条件。

发明内容

本发明所要解决的技术问题是针对现有技术的不足，提供一种多时次遥感影像优化合成方法及系统。

本发明解决上述技术问题的技术方案如下：一种多时次遥感影像优化合成方法，包括：

步骤1，获取待处理区域的待合成的n个光学遥感影像文件，所述光学遥感影像文件包括：光谱波段数据、太阳天顶角和观测天顶角；

步骤2，根据每个光学遥感影像文件中的光谱波段数据的同一波段数据生成光谱波段矩阵Ref，根据每个光学遥感影像文件中的太阳天顶角生成太阳天顶角矩阵SZA，根据每个光学遥感影像文件中的观测天顶角生成观测天顶角矩阵VZA，同时生成三个赋值为0的权重矩阵γ₁、γ₂以及γ₃；

步骤3，计算正午时刻的太阳天顶角SZA₁₂，将所述太阳天顶角矩阵SZA中各像元的数值与所述正午时刻的太阳天顶角SZA₁₂进行比较，根据第一组阈值VZA₁对第一权重矩阵γ₁进行赋值；

步骤4，将所述观测天顶角矩阵VZA中各像元的数值与第二组阈值VZA₂进行比较，根据比较结果，对第二权重矩阵γ₂进行赋值；

步骤5，对所述光谱波段矩阵Ref中每一波段层进行数值有效性计算，根据计算结果对第三权重矩阵γ₃进行赋值；

步骤6，将所述第一权重矩阵γ₁、所述第二权重矩阵γ₂以及所述第三权重矩阵γ₃进行叠加计算，生成最终权重矩阵γ，对所述最终权重矩阵γ中每个像元按波段层取最大值，并根据选取结果生成指数矩阵λ；

步骤7，依次读取n个光学遥感影像文件中同一波段的光谱波段数据，生成数值矩阵η，根据所述指数矩阵λ，提取所述数值矩阵η中的数据，将所述数据作为该波段的输出数据；

步骤8，重复步骤2至步骤7，直至所有波段均完成上述步骤，生成完整输出文件。