标准景资源三号图像融合

在有0值背景情况下的高分一号/资源三号图像融合一文中我们介绍了在有很多背景0值的情况下，可以根据像元值范围生成掩膜文件，让背景不参与融合可以得到更高的融合结果。

    标准景的资源三号全色和多光谱数据有部分区域未重叠（如下图），我们可以获取多光谱和全色的重叠区，以重叠区作为掩膜区域参与融合，最终得到重叠区的融合结果。

图1：不规则的图像

     重叠区计算工具下载：http://pan.baidu.com/s/1eQlCU8y

    将get_overlay_shapefile.sav文件拷贝到…\Exelis\ENVI51\classic\save_add 文件夹中。

      （1）在开始菜单中启动ENVI 5.1 classic。资源三号数据的全色和多光谱图像都经过正射校正和自动几何配准（6度夹角校正：资源三号卫星影像融合处理）。

        （2）选择Basic tools->Get Overlay Shapefile,如下图所示，输入全色和多光谱图像，选择输出矢量格式和路径。 

       注：如果选择shapefile格式输出，选择选择Fille->open vector，安装默认参数打开。

图4：重叠区矢量

      （3）选择主菜单->Transform->Image Sharpening->Gram-Schmidt Spectral Sharpening。
      （4）在Select Low Spatial Resolution Multi Band Input File对话框中选择低分辨率多光谱图像。同时在Select Mask Band选项中选择Mask Options->Build Mask。

 

图5：选择文件和构建掩膜文件

         （5）在Mask Definition面板中，选择Options->Import EVFs，在文件选择框中选择前面打开的重叠区矢量以及多光谱图像文件。
         （6）选择掩膜文件和路径输出。

注：如果Mask Definition面板中已有一个默认掩膜，单击Delete All Items按钮删除。

图6：通过像元值范围构建掩膜文件

           （7）在Select High Spatial Resolution Pan Input Band对话框中选择全色图像。
           （8）在弹出的Gram-Schmidt Spectral Sharpening面板中，需要选择降低高分辨率全色波段的方法，选择Average of Low Resolution Multispectral File方法。

图7：Gram-Schmidt Spectral Sharpening面板

   （9）选择重采样方法和输入路径及文件名，单击OK输出。

     如下图所示，使用掩膜的结果更清晰，色彩更真实。

图8：Gram-Schmidt融合结果

注意：Gram-Schmidt方法的大致流程：第一步，从低分辨率的波段中复制出一个全色波段。第二步，对复制出的全色波段和多波段进行Gram-Schmidt变换，其中全色波段被作为第一个波段。第三步，用高空间分辨率的全色波段替换Gram-Schmidt 变换后的第一个波段。最后，应用Gram-Schmidt反变换得到融合图像。