【技术实现步骤摘要】
一种Gram-Schmdit融合快速处理算法
本专利技术涉及一种算法,尤其涉及一种Gram-Schmdit融合快速处理算法。
技术介绍
对于光学系统的遥感影像,其空间分辨率和光谱分辨率一直存在着不可避免的矛盾。通过将较低空间分辨率的多光谱影像和较高空间分辨率的全色影像的融合,可以产生多光谱和高空间分辨率的影像。因此,各种基于不同算法的融合方法得到了迅速发展和广泛应用。随着遥感技术的发展,由于对图像解译和反演目标参数的需要,一些简单的融合方法在很大程度上已无法满足对于光谱信息保持、空间纹理信息增加的迫切需求。Gram-Schmdit(施密特正交化)融合算法借助ENVI遥感图像处理软件,并依靠组件更替战略,可使得融合影像在光谱信息及空间信息上均具有较好的保真度,因此得到了广泛推广和应用。但现有的Gram-Schmdit融合算法流程的所有中间结果均需要保存成为临时文件,以降低对普通机器内存的占用,并作为后续流程的输入再次读取使用,特别是所处理的影像为高分辨率、大数据量的全色或者较多波段多光谱数据时,大量的IO读写过程,将使得融合处理效率严重降低。
技术实现思路
为了解决上述 ...
【技术保护点】
1.一种Gram‑Schmdit融合快速处理算法,其特征在于:所述算法的整体步骤为:步骤一、简化Gram‑Schmdit变换公式:GS变换公式如下:
【技术特征摘要】
1.一种Gram-Schmdit融合快速处理算法,其特征在于:所述算法的整体步骤为:步骤一、简化Gram-Schmdit变换公式:GS变换公式如下:其中,GST是第i行第j列像素GS变换后产生的第T个分量,BT是原始多光谱影像第i行第j列像素在第T个波段的灰度值,uT是原始多光谱影像在第T个波段灰度值的均值;是第T个波段的原始多光谱影像与第l个GS分量的协方差值再与第l个GS分量方差值的比值,l的取值为1~T-1;GS逆变换公式如下:其中,BT是融合后的多光谱影像第i行第j列像素在第T个波段的灰度值,GST是第i行第j列像素GS变换后产生的第T个分量,uT是原始多光谱影像在第T个波段灰度值的均值;是第T个波段的原始多光谱影像与第l个GS分量的协方差值再与第l个GS分量方差值的比值,l的取值为1~T-1;将GS变换和GS逆变换公式进行合并,并将除了均值、第一个GS分量和替换的GS分量相关参数以外的其他参数约掉,形成新的GS融合公式,即公式3:其中,为融合后的第i行第j列的像素值,BT是第T波段的多光谱影像第i行第j列的像素值,GS1是GS第一分量,是修正后的GS第一分量,a1、a2均是融合公式系数;公式3表明,仅通过求取融合公式系数、GS第一分量、修正后的GS第一分量即可得到GS融合结果;a1、a2的具体表示如下:其中,为第T个波段的原始多光谱影像与GS第一分量的协方差再与GS第一分量方差的比值;为第T个波段的原始多光谱影像与修正后GS第一分量的协方差再与修正后GS第一分量方差的比值;步骤二、对融合公式系数进行约等化:对GS融合公式的系数a1、a2进行约等化处理的推理过程如下:在数学上,求解协方差及方差的公式分别如公式6、公式7所示:其中,X(i,j...
【专利技术属性】
技术研发人员:吴方才,周馨,白晓辉,李艳艳,
申请(专利权)人:航天星图科技北京有限公司,
类型:发明
国别省市:北京,11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。