一种基于中频信号调制和补偿的高保真遥感图像融合方法技术

本发明专利技术提供一种基于中频信号调制和补偿的遥感图像融合方法，主要解决遥感图像融合中光谱保持性能与空间保持性能之间存在矛盾、光谱保持性能较好的算法严重受制于配准精度和调制传递函数（MTF）计算精度的问题。本发明专利技术的核心机制是利用高分辨率全色数据和低分辨率多光谱数据的中频成分进行空间细节调制与注入，并以二者的线性组合进行光谱信息/空间结构信息的补偿。融合结果具有很好的光谱保持能力，对配准偏差较不敏感，可以有效抑制边缘的虚假结构。融合方法的健壮性高，光谱保持性能和空间保持性能在较高水平上达到平衡，并可根据融合任务需求通过单个参数进行性能调整，可适用于遥感及医学影像融合。

【技术实现步骤摘要】
一种基于中频信号调制和补偿的高保真遥感图像融合方法
本专利技术涉及一种基于中频信号调制和补偿的遥感图像融合技术,属于遥感影像处理数据融合

技术介绍
卫星遥感技术可以对地球表面进行大规模全覆盖的观测，在测绘、气象、海洋、农业、自然资源调查、灾害监测、国防安全等诸多方面均发挥巨大作用，世界各国也都十分重视。由于传感器信噪比的限制，光学遥感卫星所获取图像的光谱分辨率和空间分辨率之间存在矛盾，因此目前光学遥感卫星一般配置一个高分辨率全色波段和多个低分辨率多光谱波段。遥感图像融合技术是指将高空间分辨率低光谱分辨率的全色图像与低空间分辨率高光谱分辨率的多光谱图像集成，生成兼具高空间分辨率和高光谱分辨率的遥感图像的图像处理技术，又称为像素级遥感数据融合技术。以IHS、PCA、BT变换为代表的第一代遥感图像融合主要侧重于对多光谱图像进行视觉效果增强，所采用的理论假设较强，光谱扭曲比较严重；以小波变换、HPF等为代表的第二代遥感图像融合面向计算机图像处理和定量遥感应用，采用从全色图像提取空间细节，经调制后注入多光谱图像的策略，光谱保持性能远胜于第一代方法，但容易引起频谱混叠，表现为图像发虚以及空间分辨率增强不足等，全色与多光谱图像之间的配准偏差对融合结果亦会造成较大影响。目前学术界致力于构建第三代遥感图像融合框架，利用稀疏表达、压缩传感等新理论，以及机器学习等新技术，将遥感图像融合问题视为优化问题予以解决。目前第三代融合方法距离达到实用水平还有相当大的差距。第一代融合方法多属于分量替换方法(ComponentSubstitution,CS)，即将多光谱数据投影变换到某一特征空间，利用全色数据或其某种变换形式替换其中某一相关性较高的特征，并进行逆变换获得融合图像。第二代融合方法则多采用多分辨率分析技术(Multi-ResolutionAnalysis,MRA)，对全色和多光谱数据分别进行空间频谱的分解，然后依据一定的波段间空间结构模型，将全色数据中包含的空间细节成分注入到低分辨率多光谱数据中。CS方法可以较好的保持全色数据的空间细节信息，但会严重扭曲多光谱数据的光谱信息，因此融合图像的空间性能好，光谱性能差；MRA方法则相反，光谱性能好，空间性能差。有相当数量的研究试图将CS和MRA方法结合起来，构造混合模型以吸取二者的优势，但通过理论分析可知混合模型通常都与某种单一方法等价。因此也有学者认为，光谱性能和空间性能是图像融合中的一对矛盾。虽然这一猜想源于空间分辨率和光谱分辨率的矛盾，理论上确实存在对二者共同提升的上限，但现有方法呈现出这种矛盾平衡距离理论上限还有相当大的提升空间。提升MRA方法的空间性能，或改善CS方法的光谱保持性能，是目前比较可行的解决方案。目前对MRA方法的改进，主要着眼于基于MTF的滤波对MRA方法几何性能的提高，但实际的MTF值受到诸多因素影响，并且基于MTF滤波的方法并不能解决配准误差的问题。如何构建稳健的、对配准偏差较不敏感的MRA方法是一个需要解决的问题。
技术实现思路
本专利技术为解决现有基于MRA的遥感图像融合方法空间性能较差的问题，在采用适配MTF的滤波器基础上，通过添加空间/光谱补偿项，降低MTF的不确定性对融合结果的影响，减弱图像配准偏差造成的虚轮廓效应；在保持一定的光谱保持能力的基础上，将融合图像的空间分辨率提高到与CS方法相近的程度。本专利技术的技术方案为一种基于遥感图像中频信息进行空间细节调制和空间补偿的遥感影像融合方法，包含以下步骤：步骤1：对原始多光谱图像MS通过立方卷积法上采样至与全色图像P相同的尺寸，得到插值图像MSl；步骤2：根据多光谱传感器相对全色传感器的相对MTF构造高斯低通滤波器G；步骤3：图像成分提取，包括以下步骤，步骤3.1：利用G对P进行高通滤波，获取高频成分PH；步骤3.2：对PH进行低通滤波，获取中频成分PM；步骤3.3：对MSl进行高通滤波，获取高频成分步骤4：根据和PM的计算波段间结构参数W；步骤5；根据W将PH注入MSl，得到未经空间补偿的融合图像步骤6：根据和PM构造空间补偿项步骤7：根据和自定义性能参数k对进行补偿，并进行异常处理，得到融合图像。步骤1中可采用其他图像插值核，但必须保证采用的是图像插值方法，而不是信号插值方法，否则会造成半个像元的偏移。步骤2中使用相对MTF值构造高斯滤波器。令全色与多光谱数据的MTF值分别为MP和MM,则相对MTF值M＝MP/MM.由于大部分数据厂商提供的全色图像已经过MTF补偿，所以M可直接采用MM.如果数据厂商未提供数据的MTF值，可通过其他方法对图像MTF值进行测量，或采用典型值，一般为0.15到0.3.对于给定的M值，滤波器G定义为其中x,y为离散值，K为归一化系数,σ为尺度参数。σ由下式获得其中M为多光谱传感器相对全色传感器的相对MTF，r为多光谱和全色图像像元尺寸的比值，r＞1.给定滤波器抽头数p,则归一化系数为步骤3.1采用步骤2构造的滤波器进行高通滤波，通过下式实现步骤3.2和3.3分别对进行高通滤波和相应的低通滤波，给定低通滤波核H,步骤3.2的滤波过程如下：步骤3.3的滤波过程如下：H推荐采用G,也可采用其他自定义滤波核。步骤4根据和PM的计算波段间结构参数W如下其中σM为的标准差，σP为PM的标准差。步骤5根据遥感图像像素级融合统一理论框架进行细节注入，得到未经空间补偿的融合图像，细节注入方法如下式所示：步骤6根据和PB构造空间补偿项构造方法如下：步骤7根据对进行空间补偿，得到融合图像。空间补偿过程采用下式其中Fused是融合图像，k为自定义的补偿强度系数，推荐值为0.5，取值范围为(-1,1).配准偏差较大时可采用较大的k值，光谱性能要求很高时可采用较小的k值。本专利技术与现有技术相比，改进主要体现在：首先，本专利技术充分保留了MRA方法的优势，具有良好的光谱保持性能；其次，采用中频空间信息构建波段间关系，提高了空间细节调制的准确度；再次，利用中频空间信息进行信息补偿，采用缺省参数可以有效降低MTF的不确定性对融合结果的影响，提高了融合方法的稳健度；采用小于0的补偿系数可进行光谱信息增强；最后，本专利技术总体效果可以减弱图像配准偏差造成的虚轮廓效应，在保持良好的光谱保持能力的基础上，将融合图像的空间分辨率提高到与CS方法相近的程度，从而使得本专利技术的方法在空间性能和光谱性能上达到了较高水平的平衡。具体实施方式下面实施例对本专利技术进一步说明，本专利技术包括但不仅限于下述实施例。采用真实国产遥感卫星高分二号(GF-2)星载遥感多光谱和全色光图像，多光谱图像包含四个波段(MSS1～MSS4)，全色图像(PAN)为单波段。多光谱图像的空间分辨率为4.0m，大小为2500行×2500列。全色图像空间分辨率为1.0m本文档来自技高网...

【技术保护点】
1.一种基于中频信号调制和补偿的遥感图像融合方法，其特征在于包括下述步骤：/n步骤1：对原始多光谱图像MS通过立方卷积法上采样至与全色图像P相同的尺寸，得到插值图像MS

【技术特征摘要】
1.一种基于中频信号调制和补偿的遥感图像融合方法，其特征在于包括下述步骤：
步骤1：对原始多光谱图像MS通过立方卷积法上采样至与全色图像P相同的尺寸，得到插值图像MSl；
步骤2：根据多光谱传感器相对全色传感器的相对MTF构造高斯低通滤波器G；
步骤3：图像成分提取，包括以下步骤，
步骤3.1：利用G对P进行高通滤波，获取高频成分PH；
步骤3.2：对PH进行低通滤波，获取中频成分PM；
步骤3.3：对MSl进行高通滤波，获取高频成分
步骤4：根据和PM的计算波段间结构参数W；
步骤5；根据W将PH注入MSl，得到未经空间补偿的融合图像
步骤6：根据和PM构造空间补偿项
步骤7：根据和自定义性能参数k对进行补偿，并进行异常处理，得到融合图像。


2.根据权利要求1所述的一种基于中频信号调制和补偿的遥感图像融合方法，其特征在于：步骤2中采用多光谱传感器相对全色传感器的相对MTF，即多光谱和全色传感器MTF值的比值，构造高斯低通滤波器G。


3.根据权利要求1或2所述的一种基于中频信号调制和补偿的遥感图像融合方法，滤波器G定义为K为归一化系数,σ为尺度参数。


4.根据权利要求3所述的一种基于中频信号调制和补偿的遥感图像融合方法，其特征在于：σ由下式获得



其中M为多光谱传感器相对全色传感器的相对MTF，r为多光谱和全色图像像元尺寸的比值，r＞1。
...

【专利技术属性】
技术研发人员：窦闻，戚浩平，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏;32