基于参考影像约束和非凸低秩约束的遥感图像重建方法技术

本发明专利技术公开了一种基于参考影像结构约束和非凸低秩约束的遥感图像重建方法，包括如下步骤：首先建立目标影像和参考影像在高阶滤波后的结构相似约束，然后用非凸低秩近似约束代替压缩感知的

【技术实现步骤摘要】

本专利技术涉及多源遥感图像数据的信号重建方法，具体来说，涉及一种基于参考影像结构约束和非凸低秩约束的遥感图像重建方法。
技术介绍
在遥感领域，同一区域通常包含多源、多时相的影像，这些遥感图像具有不同的光谱特性、时间分辨率和空间分辨率。在一些遥感应用中，我们仅有一些观测图像，无法获取某一地区某个时刻的原始图像，如果我们需要这些原始图像的信息，需要通过重建的方法来进行遥感图像重建。空间位置对应的遥感影像在光谱和时间尺度上的冗余信息有助于重建。同一卫星内搭载的不同传感器，彼此间具有光谱重叠或相似的地物反射特性，造成拍摄的图像存在很大的相似性。同时，由于地面地物信息的迁移缓慢，不同历史时刻拍摄的同一位置的影像也具有很大的相似性。它们的这种相似性体现在结构的相似上，即图像中的边缘、轮廓等细节信息。将这些多源遥感影像的结构变化特征作为参考约束信息加入到稀疏系数约束和目标图像的重建过程，可以提高目标图像的重建精度。提取结构信息的典型方法是使用方向滤波器，如梯度算子、Gabor滤波器。但是，用于提取结构信息的滤波器应该具有最大的普适性，即普遍适用于各种图像。满足这一要求的理论是独立成分分析，它可以通过提取图像的结构信息来分离不同的图像内容，同时可训练出用于投影变换的滤波器。基于独立成分分析思想和高阶马尔可夫随机场理论的专家场模型通过从Berkeley分割数据集中提取出2万幅图像补丁训练出了两组滤波器，分别是8个3x3的滤波器和24个5x5的滤波器，表现出了普适性。相比其它滤波器，这些滤波器具有理论支持，能够更好地提取结构信息。单纯将遥感影像的结构信息作为参考约束信息注入到重建图像的主要不足是：注入过程是不加选择的，无法准确刻画光谱图像的各个局部结构、忽略稀疏系数相关性、无法刻画图像的局部稀疏度、引入吉布斯效应等问题。基于组群稀疏的低秩逼近和非凸表达解决了一般稀疏表示和引入为参考约束信息的遥感图像重建算法存在的这些问题，可以避免参考约束信息的过度注入，使注入的边缘结构信息更合理，取得更优的遥感影像纹理细节恢复。针对相关技术中的问题，目前尚未提出有效的解决方案。
技术实现思路
针对相关技术中的上述技术问题，本专利技术提出一种基于参考影像结构约束和非凸低秩约束的遥感图像重建方法。为实现上述技术目的，本专利技术的技术方案是这样实现的：一种基于参考影像结构约束和非凸低秩约束的遥感图像重建方法，包括如下步骤：包括如下步骤：首先建立目标影像和参考影像在高阶滤波后的结构相似约束，然后用非凸低秩近似约束代替压缩感知的l1范数来约束目标影像稀疏系数，建立遥感图像稀疏优化重建模型并求解。进一步的，使用一组大小为5x5或3x3的专家场滤波器组对参考图像做二维滤波，计算与目标影像相匹配的参考影像的稀疏系数，以目标影像与参考影像的专家场滤波系数相似作为约束条件。进一步的，用非凸低秩核范数约束目标影像稀疏系数，将专家场滤波器滤波后的系数相似加入到目标影像的稀疏系数中，构建目标函数。进一步的，通过共轭梯度算法、泰勒一阶近似和奇异值分解迭代求解带目标影像低秩先验信息的非局部图像重建模型。进一步的，使用矩匹配法更新参考图像的均值和标准差，使其与重建图像一致。进一步的，使用高阶滤波器提取结构信息，其所使用的滤波器是从Berkeley分割数据集中提取20000幅图像补丁并由专家场模型加以训练得到。进一步的，获取低秩相似矩阵时，所用到的相似图像块矩阵的位置关系是从参考图像中获得的。进一步的，用非凸低秩核范数约束目标影像稀疏系数，将相似度加入到目标影像的稀疏系数中进行更新，构建了重建模型的目标函数：其中，模型中第一项保证重建结果与观测数据保持匹配约束；第二项为高阶滤波系数的结构约束项，Hk为与第k个滤波器滤波过程等价的矩阵运算表示，λk为用于第k个滤波器的规整化调节系数；第三项为图像进行组群稀疏和相似块低秩约束的正则项，λ表示图像块的稀疏水平，η表示图像块相似程度拟合的权重。本专利技术的有益效果：本专利技术可将将参考影像和重建图像经专家场滤波器滤波得到的稀疏系数相似作为先验约束，将广义非凸低秩核范数作为目标影像稀疏系数约束，利用二者互补优势建立遥感图像重建模型，提高目标影像的重建精度。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是根据本专利技术实施例所述的基于参考影像结构约束和非凸低秩约束的遥感图像重建方法的操作步骤框图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本专利技术保护的范围。如图1所示，根据本专利技术实施例所述的一种基于参考影像结构约束和非凸低秩约束的遥感图像重建方法，包括如下步骤：步骤一，使用一组大小为5x5或3x3的专家场滤波器组对参考影像进行二维滤波处理，计算出与目标影像相匹配的参考影像的稀疏系数，并以目标影像与参考影像的稀疏系数相似作为约束条件；步骤二，构建重建模型的目标函数，其中，利用非凸低秩核范数约束目标影像的稀疏系数，将步骤一中得到的参考影像的稀疏系数加入到目标影像的稀疏系数中以构建重建模型的目标函数；步骤三，通过共轭梯度算法、泰勒一阶近似和奇异值分解迭代求解带目标影像低秩先验信息的非局部图像重建模型；步骤四，通过矩匹配法更新参考影像的均值和标准差，使其与重建图像一致；步骤五，迭代执行步骤一至步骤四。其中，在步骤一中，使用高阶滤波器提取结构信息，其所使用的滤波器是从Berkeley分割数据集中提取20000幅图像补丁并由专家场模型加以训练得到的，所述目标影像的结构特征向量与参考影像的结构特征向量的相似度在该高阶滤波器滤波后以l2范数的形式进行代价评估。在步骤二中，构建目标函数的过程中，获取低秩相似矩阵所用到的相似图像块矩阵的位置关系由参考影像中获得。步骤二中构建的重建模型的目标函数为：其中，模型中第一项保证重建结果与观测数据保持匹配约束；第二项为高阶滤波系数的结构约束项，Hk为与第k个滤波器滤波过程等价的矩阵运算表示，λk为用于第k个滤波器的规整化调节系数；第三项为图像本文档来自技高网...

【技术保护点】
一种基于参考影像结构约束和非凸低秩约束的遥感图像重建方法，其特征在于，包括如下步骤：首先建立目标影像和参考影像在高阶滤波后的结构相似约束，然后用非凸低秩近似约束代替压缩感知的范数来约束目标影像稀疏系数，建立遥感图像稀疏优化重建模型并求解。

【技术特征摘要】
1.一种基于参考影像结构约束和非凸低秩约束的遥感图像重建方法，其特征在于，包
括如下步骤：首先建立目标影像和参考影像在高阶滤波后的结构相似约束，然后用非凸低
秩近似约束代替压缩感知的范数来约束目标影像稀疏系数，建立遥感图像稀疏优化重建
模型并求解。
2.根据权利要求1所述的基于参考影像结构约束和非凸低秩约束的遥感图像重建方
法，其特征在于：使用一组大小为5x5或3x3的专家场滤波器组对参考图像做二维滤波，计算
与目标影像相匹配的参考影像的稀疏系数，以目标影像与参考影像的专家场滤波系数相似
作为约束条件。
3.根据权利要求2所述的基于参考影像结构约束和非凸低秩约束的遥感图像重建方
法，其特征在于：用非凸低秩核范数约束目标影像稀疏系数，将专家场滤波器滤波后的系数
相似加入到目标影像的稀疏系数中，构建目标函数。
4.根据权利要求3所述的基于参考影像结构约束和非凸低秩约束的遥感图像重建方
法，其特征在于：通过共轭梯度算法、泰勒一阶近似和奇异值分解迭代求解带目标影像低秩
先验信息的非局部图像重建模型。
5.根据权利要求4所述的基于参考影像结构约束和非凸低秩约束的遥感图...

【专利技术属性】
技术研发人员：王力哲，卢红阳，魏静波，
申请(专利权)人：中国科学院遥感与数字地球研究所，
类型：发明
国别省市：北京;11