一种基于变分方法的高光谱图像融合方法技术

本发明专利技术提供一种基于变分方法的高光谱图像融合方法，它是一种基于变分方法的高光谱图像和高分辨率可见光图像的融合方法，有三大步骤：步骤一：计算机在ＭＡＴＬＡＢ７．１环境下读取高光谱图像数据；步骤二：建立基于变分方法的泛函表达式，选择合适的权重系数；步骤三：利用分裂Ｂｒｅｇｍａｎ迭代算法对上述泛函求解极值函数，得到融合结果。本发明专利技术克服了现有技术不宜直接移植到高光谱图像融合、光谱保持性不好的缺点，取得了较好的融合结果，在高光谱图像融合技术领域里具有广阔的应用前景。

【技术实现步骤摘要】

本专利技术涉及，属于高光谱遥感图像融合

技术介绍
高光谱遥感是高光谱分辨率遥感(Hyperspectral Remote Sensing)的简称，它 是在电磁波谱的可见光和红外波段范围内，获取许多光谱连续的影像数据的技术。利用成 像光谱仪，可以采集到高光谱数据，所得高光谱数据包含了丰富的空间和光谱信息，即所谓 “图谱合一”，其突出特点在于其极高的光谱分辨率(一般在IOnm左右)，精细的光谱分辨率 反映了地物光谱的细微特征，使得在光谱域内进行遥感定量分析成为可能。与传统遥感技 术相比，高光谱遥感技术可以以较窄的波段区间、较多的波段数量提供遥感信息，在军事、 地质、环境研究等领域中得到了广泛的应用。但是，由于高光谱数据一般是来自于空基或 者天基的遥感平台，其空间分辨率一般在几米、几十米甚至上百米不等，空间分辨率较低。 另一方面，利用其它的遥感手段，我们可以得到一些空间分辨率很高的图像，比如可见光图 像。因此利用图像融合技术，在保持上述两种源图像优点的基础上，将两种图像结合，得到 光谱分辨率和空间分辨率都比较高的结果图像，具有较高的价值。目前存在的图像融合方法，大多为在像素级别针对多光谱图像的融合方法。其中 经典的方法有直接加权、主成分分析(Principal ComponentAnalysis,PCA)、小波变换等方 法。它们在多光谱图像融合领域被研究的比较充分，但是在高光谱图像融合领域，相关的研 究还比较少，多是将现有的方法直接移植，由于移植过程中一般不对数据的光谱做针对性 的研究，结果不是很理想。本专利技术针对上述情况，提出一种改进的基于相似性度量的变分方 法以实现高光谱图像融合，融合效果较好。
技术实现思路
1、目的本专利技术的目的是提供，该方法 利用光谱约束和图像细节约束来进行高光谱图像融合。2、技术方案本专利技术是通过以下技术方案实现的本专利技术一种基于变分的高光谱图像融合方法，它包括如下步骤步骤一计算机读取数据。首先计算机在MATLAB 7. 1下读取成像光谱仪所得到的 高光谱图像数据(本专利技术高光谱数据来源于中国HJ-IA卫星所搭载成像光谱仪采集到的青 海省托素湖遥感图像)，再读取已经手工配准过的高分辨率图像数据(本专利技术数据来源于 Google-Map提供的高分辨率影像)。本专利技术中的高光谱数据为100谱段*256像素*200像素大小的三维数据立方体， 其中100表示谱段数，256*200表示每一谱段图像大小，以下用Hn表示第η个谱段的高光谱 源图像，用Un表示第η个谱段融合的高光谱图像，用u表示一幅普通的图像；高分辨率数据 大小为256像素*200像素，以下用M表示。4步骤二 建立泛函表达式。基于变分方法处理图像的基本框架，首先建立高光谱图 像融合的基本泛函表达式，它包括三项第一项为图像的细节保持项Ea，该项用以实现对高光谱数据几何形状的保持。形 式如下toon] 式中，δ，η为预先设定的常系数，Un表示所求第η个谱段的融合图像，θ表示图像中每个像素点处的单位法向量。其定义如下 Vm表示对u求梯度，图像的水平和垂直方向分别用χ和y表示，定义如下 div(6 )表示对θ求散度，具体求解时采用下式进行，div( θ ) = θ χ+ θ y以上两式中，ux，uy, θ χ，θ y分别表示u对χ，y的偏导数，θ对x，y的偏导数。对 于图像而言，偏导数的计算可以用差分代替。第二项为基于与源图相似性度量的保真项Eb，用以实现对高光谱数据灰度的保 持，形式如下 式中，υ，ρ为常系数，P用以平衡源数据中高光谱数据与高分辨率图像数据在 融合结果中所占的不同比重。M表示高分辨率图像，Hn表示第η个谱段源高光谱数据。第三项为光谱保持项Ε。，该式用以实现对高光谱数据的光谱约束。形式如下 式中，μ为预先设定常数，Ui, Uj表示融合后第i，j谱段高光谱数据的图像，Hi, Hj 表示原始高光谱数据第i，j两个谱段的图像。将以上三项直接加和，建立的泛函表达式为E = Ea+Eb+Ec步骤三对此能量泛函求极值函数，则得到的极值函数即为所求融合图像，本专利技术 利用分裂Bregman迭代算法对此泛函求解极值函数。具体步骤如下(Dnz1，开始循环(2)参量初始化 (3)计算第k+Ι次迭代得到的u (即融合图像结果) 式中，λ为选定常数，参量d和b为迭代中产生的变量。<这样取值，若j <n，则取心=Uk/1，若 j > n，则取< =Ukj。(4)进行参量d和b的更新 (5)迭代停止条件。当融合结果量u在相邻两次迭代的过程中变化不大时停止迭 代，本专利技术中停止条件为 其中uk+1表示当前迭代结果，Uk表示前一次迭代结果。本专利技术中，取tol = 0. 4。 得到第η个谱段的融合结果，存入新的数据立方体中。(6)η = η+1，继续第(2)步，直到η > N，得到融合的高光谱数据立方体，至此完成 对高光谱图像与高分辨率可见光图像的融合工作。3、优点及功效本专利技术的优点是针对常见图像融合方法不宜直接移植到高光谱 图像融合，而已有基于变分方法的融合手段建模复杂，求解繁琐，对图像配准要求高的特 点，通过引入基于与源图相似性度量量来进行改进，通过对高光谱数据的光谱保持从光谱 角度进行约束，以相似性度量对图像质量进行约束，最终得到了图像和光谱质量都比较好 的融合结果。通过随机选择的实际两像素点融合前后的光谱曲线可看出，使用本方法进行高光谱图像融合，光谱保持性非常好，事实上，几乎所有的像素点都保持了融合前光谱曲线的基本特征，利用光谱角度进行光谱约束，得到结果图像与源高光谱图像各像素点中的光谱角 度余弦平均值为0. 9941，同时融合图像的视觉质量也比较好，因而利用专利技术的变分方法处 理高光谱图像融合问题效果比较好。附图说明图1本专利技术所述方法融合流程图2 (a)本专利技术实验高光谱数据(129，100)点处融合前后的光谱曲线，细虚线为融 合前光谱曲线，粗实线为融合后光谱曲线图2(b)本专利技术实验高光谱数据(50，128)点处融合前后的光谱曲线，细虚线为融 合前光谱曲线，粗实线为融合后光谱曲线具体实施例方式为了更好地理解本专利技术的技术方案，以下结合附图对本专利技术的实施方式作进一步 描述本专利技术在MATLAB7. 1语言环境下实现。计算机读取高光谱遥感图像数据和配准好 的高分辨率图像数据，首先基于与源图像的相似性度量建立正确的能量泛函表达式，选择 合适的权重系数，然后对泛函表达式求解极值函数，极值函数即为所求的高光谱融合图像。 在求解过程中，本专利技术延续已有方法中的分裂Bregman迭代算法。本专利技术，该方法的流程见图1所示，计 算机配置采用=Intel dual-core处理器，主频1. 6GHz，内存2GB。该融合方法包括以下步 骤步骤一、用计算机读取数据。计算机在MATLAB7. 1环境下读取高光谱图像数据和 高分辨率图像数据。步骤二、建立基于变分方法的能量泛函。建立基于与源图像相似性度量的能量泛 函，选择合适的权重系数。本专利技术实验中选择的参数如下Y = 1，η = 0. 5，ν = 2，P = 4, μ = 2, λ = 1。步骤三、对此泛函求解极值函数。利用分裂Bregman迭代求解。(l)n=l，开始迭代(2)参量初始化 (3)本文档来自技高网...

【技术保护点】
一种基于变分方法的高光谱图像融合方法，其特征在于：将高光谱图像融合的问题转化为对泛函求解极值函数的优化问题，同时考虑了光谱约束和图像的细节约束，该方法具体步骤如下：步骤一：计算机读取数据；首先计算机在ＭＡＴＬＡＢ７．１下读取成像光谱仪所得到的高光谱图像数据，再读取已经手工配准过的高分辨率图像数据；其中的高光谱数据为１００谱段＊２５６像素＊２００像素大小的三维数据立方体，１００表示谱段数，２５６＊２００表示每一谱段图像大小，以下用Ｈ↓［ｎ］表示第ｎ个谱段的高光谱源图像，用ｕ↓［开始循环（２）参量初始化：ｕ↑［０］＝Ｈ↓［ｉ］，（ｄ↓［ｎ］↑［０］）↓［ｘ］＝（ｄ↓［ｎ］↑［０］）↓［ｙ］＝（ｂ↓［ｎ］↑［０］）↓［ｘ］＝（ｂ↓［ｎ］↑［０］）↓［ｙ］＝０（３）计算第ｋ＋１次迭代得到的ｕ（即融合图像结果）：（２υ＋２υ＊ρ＋２μ＊（Ｈ↓［ｊ］）↑［２］－λΔ）ｕ↓［ｎ］↑［ｋ＋１］＝２υ＊Ｈ↓［ｎ］＋２υ＊ρ＊Ｍ－ηｄｉｖ（θ）＋２μ＊Ｈ↓［ｎ］（＊ｕ↓［ｊ］↑［ｋ↑［＊］］＊Ｈ↓［ｎ］）＋λｄｉｖ（ｄ↓［ｎ］↑［ｋ］－ｂ↓［ｎ］↑［ｋ］）式中，λ为选定常数，参量ｄ和ｂ为迭代中产生的变量。ｕ↓［ｊ］↑［ｋ↑［＊］］这样取值，若ｊ＜ｎ，则取ｕ↓［ｊ］↑［ｋ↑［＊］］＝ｕ↓［ｊ］↑［ｋ＋１］，若ｊ＞ｎ，则取ｕ↓［ｊ］↑［ｋ↑［＊］］＝ｕ↓［ｊ］↑［ｋ］。（４）进行参量ｄ和ｂ的更新：（ｄ↓［ｎ］↑［ｋ＋１］）↓［ｘ］＝ｍａｘ（ｓ↑［ｋ］－１／λ，０）＊（▽↓［ｘ］ｕ↓［ｎ］↑［ｋ＋１］＋（ｂ↓［ｎ］↑［ｋ］）↓［ｘ］）／ｓ↑［ｋ］（ｄ↓［ｎ］↑［ｋ＋１］）↓［ｙ］＝ｍａｘ（ｓ↑［ｋ］－１／λ，０）＊（▽↓［ｙ］ｕ↓［ｎ］↑［ｋ＋１］＋（ｂ↓［ｎ］↑［ｋ］）↓［ｙ］）／ｓ↑［ｋ］（ｂ↓［ｎ］↑［ｋ＋１］）↓［ｘ］＝（ｂ↓［ｎ］↑［ｋ］）↓［ｘ］＋（▽↓［ｘ］ｕ↓［ｎ］↑［ｋ＋１］－（ｄ↓［ｎ］↑［ｋ＋１］）↓［ｘ］）（ｂ↓［ｎ］↑［ｋ＋１］）↓［ｙ］＝（ｂ↓［ｎ］↑［ｋ］）↓［ｙ］＋（▽↓［ｙ］ｕ↓［ｎ］↑［ｋ＋１］－（ｄ↓［ｎ］↑［ｋ＋１］）↓［ｙ］）其中ｓ↑［ｋ］＝＊＊＊（５）迭代停止条件。当融合结果量ｕ在相邻两次迭代的过程中变化不大时停止迭代，本专利技术中停止条件为：｜｜ｕ↑［ｋ＋１］－ｕ↑［ｋ］‖＞ｔｏｌ其中ｕ↑［ｋ＋１］表示当前迭代结果，ｕ↑［ｋ］表示前一次迭代结果。本专利技术中，取ｔｏｌ＝０．４。得到第ｎ个谱段的融合结...

【技术特征摘要】

【专利技术属性】
技术研发人员：史振威，安振宇，姜志国，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：11[中国|北京]