[发明专利]一种基于统计分布的卫星多光谱图像波段配准方法

说明书

技术领域

本发明涉及一种卫星多光谱图像配准的处理方法。

背景技术

遥感技术目前已在社会和科技层面转向实际应用，这些应用包括自然灾害的处理、气候变化评估、自然资源管理、环境保护等，所有这些都涉及长期监测地球表面。近些年，图像配准在遥感应用中变的非常重要。图像配准是图像处理中的基本任务，指的是匹配两幅或多幅来自不同时间、不同遥感器、不同视角的同一物体或场景的图像。图像配准用于数字图像处理是为了把两幅或多幅数字图像准确对准以便分析与比较，涉及生理学、计算机视觉、模式识别、影像理解等多个领域知识。精确的配准算法对于支持镶嵌遥感卫星图像、追踪地球表面环境变化、基础科学研究是非常重要的。

图像配准，即通过计算一组变换参数把两幅图像对准，这个问题看似定义简单明了，似乎应该已有清楚的、通用的方法，而实际上远非如此。由于对应于各种不同数据的应用众多，图像配准已经发展成为一项复杂的、具有很强挑战性的、包含许多方法策略的任务。随着遥感、医学、以及其他领域获取图像能力的不断增强，导致在过去20年里对图像配准技术进行了大量研究。但到目前为止,还没有一种配准方法能够解决所有的配准问题,只能根据具体的数据类型和应用来研究相应的算法。图像配准算法常分为基于区域和基于特征的两种方法,但是一般只适用于灰度差别较小的影像，反映数据之间的线性特征，不适用于灰度差别较大的影像之间的配准。

发明内容

本发明的技术解决问题是：克服现有技术的不足，提供一种基于统计分布的卫星多光谱图像波段配准方法。

本发明的技术解决方案是：一种基于统计分布的卫星多光谱图像波段配准方法，步骤如下：

(1)以多光谱图像的任一一个谱段图像为参考图像,其他谱段图像作为待配准图像，分别对参考图像和待配准图像用高斯模板函数进行降采样生成不小于2层的不同分辨率的图像；对每个待配准图像分别按照下述步骤进行处理：

(2)在第1层低分辨率图像层，利用核密度函数计算参考图像和待配准图像的边缘概率分布和联合概率分布，并计算参考图像和待配准图像之间的负互信息值；以负互信息为目标函数进行叠代优化，求取最小的相似性测度值或达到规定的叠代次数时的平移变换参数；

(3)在第二层分辨率图像层，利用上一层得到的平移变换参数对本层中的待配准图像做变换，然后再利用核密度函数计算参考图像和待配准图像的边缘概率分布和联合概率分布，并计算参考图像和待配准图像之间的负互信息值；以负互信息为目标函数进行叠代优化，求取最小的相似性测度值或达到规定的叠代次数时的刚体变换参数；

(4)在剩余高分辨率图像层，从低分辨率开始，分别利用上一层得到的刚体变换参数对本层中的待配准图像做变换，然后再利用核密度函数计算参考图像和待配准图像的边缘概率分布和联合概率分布，并计算参考图像和待配准图像之间的负互信息值；以负互信息为目标函数进行叠代优化，求取最小的相似性测度值或达到规定的叠代次数时的刚体变换参数；在最后一层利用变换参数对待配准图像进行变换、重采样得到配准后的图像。

本发明与现有技术相比的优点在于：

(1)本发明针对卫星多光谱数据的特点，针对灰度差异大的多光谱图像（如蓝波段和近红外波段）配准，从配准的重要方面相似性测度进行突破，采用新的相似性测度准则，即基于图像统计分布信息的方法，有效的解决了传统图像配准方法相似性测度要求图像灰度相似、成线性关系的局限，该方法不需要对影像做分割、特征提取等预处理，适用面比较广，具有较高的精度和较好的鲁棒性。

（2）本发明采用核密度函数估计变量的概率密度。最早的利用直方图近似变量的概率密度的方法具有较大的估计误差，同时存储直方图所需的空间随着样本的特征变量数目成指数增长。

（3）本发明采用分层变换的方法在金字塔不同层采用不同的变换，先从整体上作平移变换迅速缩小搜索范围，然后在原始高分辨率图像层采用仿射变换。

（4）本发明采用多分辨率配准策略通过由粗到精的方式解决配准问题，可以避免互信息局部极值提高配准精度，同时提高了配准算法的速度和鲁棒性。

附图说明

图1为本发明原理示意图；

图2为本发明方法流程图。

具体实施方式

本发明实现了以灰度分布统计信息为基础的互信息作为配准相似性测度的配准方法，它不需要对影像做分割、特征提取等预处理，适用面比较广，具有较高的精度和较好的鲁棒性。同时，为了提高配准精度和速度，本发明将互信息实现方法和图像金字塔算法有机结合起来，使得单一的配准过程变成多分辨率的从粗到细的配准过程。