遥感图像双线性插值小波变换融合的优化方法技术

一种遥感图像双线性插值小波变换融合的优化方法，采用双线性插值算式对高光谱图像进行重采样优化计算，使其空间分辨率等于全色图像，根据小波分解理论对全色图像与重采样后的高光谱图像进行小波分解，获取各自的近似低频和细节高频分量。遥感图像的小波系数融合时，将高光谱图像的小波分解的近似分量作为新图像的近似分量，而细节分量融合，则利用线性计算来获取，根据小波系数融合所得新图像数据块的小波系数，进行小波重建和数据块拼接，生成所需的遥感融合图像。本发明专利技术在保留高光谱图像的光谱信息的同时，能够提高其空间分辨率，同时显著提高了计算机处理速度，获得的融合图像质量更加可靠。（*该技术在2023年保护过期，可自由使用*）

【技术实现步骤摘要】

本专利技术涉及一种，是高光谱遥感图像经双线性重采样后，与高空间分辨率全色图像融合时进行高精度及快速计算的融合优化方法，是多源遥感影像后处理阶段的一项核心技术，在各类军用或民用的遥感信息处理系统、数字城市空间信息系统等领域中均可有广泛的应用。一般地说，遥感图像融合就是图像合成技术，将不同平台(卫星与机载)上的同一或不同传感器获取的不同空间与光谱分辨率图像按特定的算法进行处理，以使所产生的新图像同时具有原来图像的多光谱特性以及高地面分辨率，来实现不同的应用需求。目前，应用较多的遥感图像融合是不同分辨率的多光谱(Multi-Spectral)与PAN(Panchromatic)图像融合或高光谱(Hyper-Spectral)图像与PAN图像的融合，主要方法有IHS变换、主成分分析(PCA)、BROVERY变换以及小波变换合等。对这些融合算法进行归纳，它们的一个共同特征是，融合前均将高光谱图像重采样为与全色图像空间分辨率相等的新图像，以供进一步处理。为此，面临两个方面的问题需要解决1.高光谱图像重采样算法对融合图像的光谱及空间分辨率影响；2.由于高光谱图像的大小一般为几百兆的数量级，若其重采样算法变得复杂时，须考虑计算机的运算速度。目前，进行上述遥感图像融合时，在精度许可的条件下，主要应用近邻插值方法，以满足算法的处理速度要求。这里采用的近邻插值算法的原理是，像素点插值搜索过程中，若高光谱图像中某素点离重采样算法所得的新图像中某像素点的空间位置最近，则将其像素值作为新图像对应点的像素值。该方法通常情况下可以满足图像处理要求。但若要求光谱信息损失较小时，则它难以达到所需的融合精度，究其原因，是由于近邻算法本身的局限性造成的，所以需研究新的算法。为实现这样的目的，本专利技术的技术方案中，利用遥感图像处理与开发平台提供的开发语言，针对高光谱图像的文件大小及计算机的处理速度，采用双线性插值算式对高光谱图像进行重采样优化计算，使其空间分辨率等于全色图像。然后，根据小波分解理论对全色图像与重采样后的高光谱图像进行小波分解，获取各自的近似低频和细节高频分量。遥感图像的小波系数融合时，将高光谱图像的小波分解的近似分量作为新图像的近似分量，而细节分量融合，则利用线性计算来获取，这样有利于保留高光谱图像的光谱信息。最后，根据小波系数融合所得新图像数据块的小波系数，进行小波重建和数据块拼接，生成所需的遥感融合图像。本专利技术的方法包括如下具体步骤1.数据特性的提取根据应用的需求，获得配准后待融合的高光谱和全色遥感图像文件，通过遥感图像处理与开发平台提供的接口函数，取得它们的特性数据，包括文件大小，空间分辨率以及图像高光谱遥感图像的波段数等。2.高光谱图像的双线性重采样以全色图像的空间分辨率为标准，获取计算机的内存与CPU处理速度等指标后，根据高光谱文件对应的数据量的大小将其分割为几个连续的数据块(数目为计算机内存量/高光谱图像数据量)，按双线性插值规则，对各数据块分别逐点进行搜索、插值，生成具有与全色图像空间分辨率相同的新高光谱图像数据块。3.全色图像的数据分块，按照第二步对应的数据分块数，将全色图像分成与经双线性插值后所得的新高光谱图像数据块数目相等的数据块。4.循环调入对应的全色与高光谱图像数据块，进行各数据块的二维小波分解。就是在水平和垂直方向对全色与高光谱图像数据块各自独立地进行一次一维小波变换。相当于先在水平方向对图像分别进行低通和高通滤波，并作下采样，然后在垂直方向上实施同样的操作。最终获得各个尺度下的低频近似分量和三个方向(水平、垂直、斜线)的高频细节分量，它们构成对应数据块的小波系数。5.高光谱与全色图像各数据块的小波系数融合将高光谱图像数据块的低频近似分量作为融合图像数据块的低频分量，而融合图像数据块的高频分量，则使用线性加权算法由全色与高光谱图像数据块的高频细节分量取得，由此构建新图像的小波系数(低频与高频分量)。6.融合图像数据块的小波变换重建及数据块的拼接将获得的融合图像数据块的小波系数按小波变换的重构原理，进行图像重构，获得所需的高空间分辨率和光谱分辨率的数据块，然后将这些数据块按前述高光谱图像的分割顺序进行拼接，组成新图像。本专利技术的特色在于，由于同一地区高光谱与全色遥感图像的空间分辨率及光谱分辨率等特性上相差很大，按双线性规则对高光谱图像进行重采样后，利用小波变换进行融合，在保留高光谱图像的光谱信息的同时，能够提高其空间分辨率，使地物空间结构和纹理特性能更加明确的显示。同时，由于进行了数据分块，将各数据块依次调入计算机内存，能显著地提高计算机处理速度。总之，本专利技术提出的方法可广泛地用于遥感图像数据的融合处理，获得更加可靠的融合图像质量。如附图说明图1所示，首先利用基于Windows的ENVI/IDL遥感图像处理与开发平台提供的接口函数来获取高光谱与全色图像的特性数据。将高光谱图像按一定的比例进行数据分块后对各块再实施双线性重采样，使各数据块的空间分辨率等于全色图像的空间分辨率。此后，按相同的比例对全色图像进行数据分块。循环调入全色与重采样后高光谱图像数据块，对各对应的数据块进行小波分解，获得小波分解的低频近似分量与高频细节分量。循环调入小波分解后的数据块进行小波系数融合，构造新图像数据块的小波系数。最后，对新图像数据块进行小波变换重建和数据块的拼接，得到所需的融合图像。图2为双线性重采样算法所使用的插值模型。图3是全色和高光谱图像对应数据块小波变换融合的计算过程示意图。其中Cap、Chp、Cvp、Cdp分别为全色图像数据块小波分解后的低频分量、高频水平、高频垂直和高频斜线分量；Cam、Chm、Cvm、Cdm分别为高光谱图像重采样后数据块小波分解后的低频分量、高频水平、高频垂直和高频斜线分量。图4为本专利技术采用的方法与其它方法的融合效果对比。其中，图4(a)的Land TM图像为高光谱图像，本实施例取其波段4、3、2，文件大小1007×560，空间分辨率28m。图4(b)的SPOT遥感影像是全色图像，文件大小2820×1569，窨分辨率10m。图4(c)是近邻插值小波融合效果，图4(d)是双线性插值小波融合效果，从中可以看出，双线性插值小波变换融合有利于保留高光谱图像的光谱信息。图1为本专利技术提出的一种用于高光谱与全色图像双线性插值小波变换融合优化方法的总体框图，数据文件的输入是图4(a)的Land TM图像为高光谱图像，取其波段4、3、2，文件大小1007×560，空间分辨率28m。图4(b)的SPOT遥感影像是全色图像，文件大小2820×1569，空间分辨率10m。各部分具体实施细节如下1.数据特性的提取，获得高光谱与全色图像的特性数据利用ENVI/IDL开发平台的ENVI DISPLAY BANDS函数，可以显示出文件的大小、波段数及空间分分辨率等。2.高光谱图像的双线性重采样，利用ENVI/IDLIDL的C语言接口调用WINDOWS2000操作系统提供的API函数，获取计算机的硬件特性，如CPU处理速度、内存大小后，对高光谱图像按图2所示的方法进行双线性重采样以及按内存/高光谱文件数据量的比例进行高光谱数据文件的分块。对于高光谱图像中某像素点(x，y)及其周围的点(0，0)、(0，1)、(本文档来自技高网...

【技术保护点】
一种遥感图像双线性插值小波变换融合的优化方法，其特征在于包括如下具体步骤：１）获得配准后待融合的高光谱和全色遥感图像文件后，通过遥感图像处理与开发平台提供的接口函数，取得它们的特性数据，包括文件大小，空间分辨率以及图像高光谱遥感图像的波段数等；２）以全色图像的空间分辨率为标准，获取计算机的内存与ＣＰＵ处理速度等指标，根据高光谱文件对应的数据量的大小将其按计算机内存量／高光谱图像数据量的比例值分割为几个连续的数据块，对各数据块分别逐点进行双线性重采样，生成具有与全色图像空间分辨率相同的新高光谱图像数据块；３）将全色图像分成与经双线性插值后所得的新高光谱图像数据块数目相等的数据块；４）循环调入对应的全色与高光谱图像数据块，进行各数据块的二维小波分解，即在水平和垂直方向对其各自独立地进行一次一维小波变换，先在水平方向对图像分别进行低通和高通滤波，并作下采样，然后在垂直方向上实施同样的操作，获得各个尺度下的低频分量和水平、垂直、斜线三个方向的高频分量，构成对应数据块的小波系数；５）进行高光谱与全色图像各数据块的小波系数融合，将高光谱图像数据块的低频分量作为融合图像数据块的低频分量，而融合图像数据块的高频分量，则使用线性加权算法由全色与高光谱图像数据块的高频系数取得，由此构建新图像的小波系数；６）将获得的融合图像数据块的小波系数按小波变换的重构原理，进行图像重构，获得所需的空间分辨率和光谱分辨率数据块，然后将这些数据块前述的分割顺序进行拼接，组成新图像。...

【技术特征摘要】

【专利技术属性】
技术研发人员：敬忠良，周前祥，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：31[中国|上海]