一种改进型多源遥感数据时空融合方法
【技术领域】
[0001] 一种改进型多源遥感数据时空融合方法，属于数字图像处理技术领域，特别设及 遥感图像融合及遥感图像分类技术。
【背景技术】
[0002] 多源遥感数据融合是利用不同空间分辨率、不同时间分辨率和不同光谱分辨率的 光学遥感数据，甚至是光学遥感数据与雷达数据所各自具有的不同数据特点，融合出能够 综合多源数据特征的融合影像的方法。它对提高遥感影像利用效率和提高遥感应用效果都 具有重要的意义。
[0003] 多源遥感数据时空融合是多源遥感数据融合领域的前沿研究内容之一。目前多源 遥感数据融合主要是全色高空间分辨率影像与多光谱影像间的融合，融合的目的只要为提 高多光谱影像的空间分辨率，融合的方法主要有代数法和基于分量替换的方法等。
[0004] 代数法是通过对待融合影像进行代数运算来实现融合，如利用乘法运算进行SAR 影像和光学影像的融合，利用相关系数加权方法进行全色影像和多光谱影像的融合等。
[0005] 基于分量替换的融合方法是将影像先进行某种变换，然后利用高空间分辨率影像 替换变换空间中的一个分量，在进行逆变换来实现融合。分量替换的融合方法主要包括HIS 变换融合和小波变换融合等。
[0006] HIS变换融合（参见文献；赵英时等.遥感应用分析原理与方法.北京：科学出版 社，2003)是将图像彩色空间RGB空间变换到HIS空间，再用全色波段替换HIS空间中的强 度分量I分量，再进行HIS逆变换实现影像的融合。其目的主要为提高多光谱影像的空间 分辨率。
[0007] 小波变换融合（参见文献；赵英时等.遥感应用分析原理与方法.北京：科学出版 社，2003)是利用小波变换具有变焦性、信息保持性和小波基选择的灵活性等优点，将图像 分解为一些具有不同空间分辨率、频率特性和方向持性的子信号（图像），其中它的分频特 征，相当于高、低双频滤波器，能够将一信号分解为低频信息（图像）和高频细节/纹理信 息（图像），同时又不失原信号所包含的信息，将高空间分辨率影像的小波变换后的高频部 分代替多光谱影像的小波变换高频部分，再进行小波逆变换来实现融合。小波变换融合可 W用于W非线性的对数映射方式融合不同类型的图像数据，使融合后的图像既保留原高分 辨率适感影像的结构信息，又融合多光谱影像丰富的光谱信息，提高影像的解译能力、分类 精度。
[000引基于分量替换的融合方法还包括基于G-S变换的影像融合、基于主成分分析的影 像融合和基于高通滤波影像融合等。
【发明内容】
[0009] 本发明提供一种改进型多源遥感时间的时空融合方法，用W实现利用时序低空间 分辨率影像和高空间分辨率影像，融合生成高时空分辨率影像，解决目前由于遥感卫星数 据获取能力不足和云等天气条件的影响造成的遥感数据缺乏问题。
[0010] 本发明技术方案如下：
[0011] 一种改进型多源遥感数据时空融合方法，其特征在于包括W下步骤：
[0012] 步骤一、丰度矩阵提取。基于高分辨率遥感分类图，W低空间分辨率像元大小提取 每个低空间分辨率像元的丰度矩阵；高分辨率遥感分类图可W是化CD和FROM-GLC分类数 据，也可W基于第1期高分辨率影像采用监督分类方法提取。
[0013] 步骤二、自适应窗口和步长选取。基于线性光谱混合模型，输入时序低空间分辨率 像元反射率、中空间分辨率像元反射率和丰度矩阵，采用线性光谱混合模型解算窗口大小 和步长移动大小自适应选取方法，确定线性光谱混合模型解算的最佳窗口大小和步长移动 大小。
[0014] 步骤=、地物反射率的时间曲线提取。利用时序低空间分辨率影像和丰度矩阵， 采用最小二乘法，根据步骤二确定的线性光谱混合模型解算的最佳窗口大小和步长移动大 小，逐窗口、逐步长解算各个窗口的类别时序平均反射率和第1期的高空间分辨率类别平 均反射率，将类别时序平均反射率作为地物反射率的时间曲线。
[0015] 步骤四、传感器误差校正和高空间分辨率类别平均反射率计算。根据第1期的线 性光谱混合模型解算地物反射率和地物中空间分辨率像元平均反射率，采用线性回归方法 构建线性光谱混合模型解算地物反射率和地物中空间分辨率像元平均反射率的关系模型， W纠正不同传感器差异造成的反射率差异影响，根据线性光谱混合模型解算地物反射率和 地物中空间分辨率像元平均反射率的关系模型，从步骤（3)解算的各类别的地物反射率时 间变化曲线计算出地物中空间分辨率像元平均反射率的时间曲线；
[0016] 步骤五、融合影像生成。利用地物中空间分辨率像元平均反射率的时间曲线，根据 像元反射率和类别平均反射率之间的关系模型，计算出每个像元的时间反射率曲线，融合 出时序高空间影像；
[0017] 步骤六、融合效果精度评价。利用真实中分辨率影像，采用相关系数、方差、平均绝 对误差、MSE和偏差等参数对融合影像，对融合影像进行精度评价。
[001引本发明与现有技术相比所具有的优点；本发明采用自适应窗口和步长选取方法自 动选取每个类别的最佳线性光谱混合模型解算窗口和步长大小，根据线性回归方法纠正了 不同传感器造成的误差；与传统方法相比，本发明不仅可W提高融合影像的空间分辨率，还 可W提高影像的时间分辨率，从而可W解决目前由于卫星传感器数据获取能力不足、W及 光学卫星容易受到云的影响而导致的遥感数据缺失问题。
【附图说明】
[0019] 图1 ;本发明的一种改进型多源遥感时间的时空融合方法流程图；
[0020] 图2 ;融合结果图；
[0021] 图3 ;融合结果精度评价图；
【具体实施方式】
[0022] 为了更好地理解本发明的技术方案，下面结合附图及【具体实施方式】详细介绍本发 明。
[0023] 本发明一种改进型多源遥感时间的时空融合方法，该方法主要包括W下几个步 骤：
[0024] 1、丰度矩阵提取；
[0025] 2、自适应窗口和步长选取；
[0026] 3、地物反射率的时间曲线提取；
[0027] 4、传感器误差校正和高空间分辨率类别平均反射率计算；
[002引 5、融合影像生成；
[0029] 6、融合效果精度评价；
[0030] 本发明的具体实现流程如图1所示，各部分具体实施细节如下：
[0031] 1、丰度矩阵提取
[0032] 根据线性混合模型，每一个低空间分辨率像元的反射率R(i，t)可W看做各组分 平均反射率;：(c，/)的线性组合；
【主权项】
1. 一种改进型多源遥感数据时空融合方法，其特征在于包括w下步骤： (1) 基于高分辨率遥感分类图，W低空间分辨率像元大小提取每个低空间分辨率像元 的丰度矩阵； (2) 基于线性光谱混合模型，输入时序低空间分辨率像元反射率、中空间分辨率像元反 射率和丰度矩阵，采用线性光谱混合模型解算窗口大小和步长移动大小自适应选取方法， 确定线性光谱混合模型解算的最佳窗口大小和步长移动大小； (3) 根据自适应窗口和步长选取方法确定的线性光谱混合模型解算的最佳窗口大小和 步长移动大小，逐个窗口、逐个移动步长移动，解算线性光谱混合模型，获得各类别的地物 反射率时间变化曲线和第1期的地物中空间分辨率像元平均反射率； (4) 根据第1期的线性光谱混合模型解算地物反射率和地物中空间分辨率像元平均反 射率，采用线性回归方法构建线性光谱混合模型解算地物反射率和地物中空间分辨率像元 平均反射率的关系模型，W纠正不同传感器差异造成的反射率差异影响，根据线性光谱混 合模型解算地物反射率和地物中空间分辨率像元平均反射率的关系模型，从步骤（3)解算 的各类别的地物反射率时间变化曲线计算出地物中空间分辨率像元平均反射率的时间曲 线； (5) 利用地物中空间分辨率像元平均反射率的时间曲线，根据像元反射率和类别平均 反射率之间的关系模型，计算出每个像元的时间反射率曲线，融合出时序高空间影像； (6) 利用真实中分辨率影像，采用相关系数、方差、平均绝对误差、MSE和偏差等参数 对融合影像，对融合影像进行精度评价。
2. 根据权利要求1所述步骤似自适应窗口和步长选取，其特征在于包括W下4个步 骤： (1) 初始化；给定线性光谱混合模型解算的窗口最大、最小和窗口改变大小=个参数 的初始值；给定窗口移动步长的最大、最小和步长改变大小=个参数的初始值； (2) 根据初始化的6个参数，根据窗口大小选取出第1期的低空间分辨率像元反射率、 第1期的中空间分辨率像元反射率和丰度矩阵，根据线性光谱混合模型解算该窗口条件下 一个窗口的类别平均反射率和中空间类别平均反射率；根据窗口步长大小，在整个影像上 移动窗口，并计算每个窗口的类别平均反射率和中空间类别平均反射率；计算给定窗口大 小和步长大小条件下的类别平均反射率和中空间类别平均反射率的相关系数； (3) 遍历；根据给定的窗口改变大小和步长改变大小两个参数，改变窗口的大小和步 长大小，并计算给定窗口大小和步长大小条件下的类别平均反射率和中空间类别平均反射 率，W及二者之间的相关系数； (4) 选取类别平均反射率和中空间类别平均反射率相关系数最大的窗口大小和步长大 小最为最佳的窗口大小和步长大小参数。
3. 根据权利要求1所述步骤（4)传感器误差校正和高空间分辨率类别平均反射率时间 曲线计算，其特征在于传统多源遥感数据时空融合方法未考虑传感器差异对融合模型的精 度影响，本分发明采用线性回归方法构建类别平均反射率和中空间类别平均反射率的关系 模型，通过该模型校正传感器差异造成的误差，并计算出高空间分辨率类别平均反射率时 间曲线。
【专利摘要】本发明公开了一种改进型多源遥感数据时空融合方法。该方法利用从分类图上提取的丰度矩阵，采用自适应窗口和移动步长选择方法，逐窗口解算线性光谱混合模型，获得各类别的地物反射率时间变化曲线；采用线性回归方法构建类别平均反射率与地物中空间分辨率像元平均反射率的关系模型，并根据该模型计算时序的地物中空间分辨率平均像元反射率；再结合像元反射率与类别平均反射率间的关系模型，融合出高时空分辨率融合影像；最后利用真实中分辨率影像，采用相关系数、方差、平均绝对误差、RMSE和偏差等参数对融合影像进行精度评价。利用该方法可以从中空间分辨率影像与高时间分辨率影像融合出高时间分辨率影像的方法。
【IPC分类】G06T5-50
【公开号】CN104715467
【申请号】CN201510099032
【发明人】邬明权, 牛铮
【申请人】中国科学院遥感与数字地球研究所
【公开日】2015年6月17日
【申请日】2015年3月6日