标准景资源三号图像融合
分享
上一篇文章中我们介绍了在有很多背景0值的情况下,可以根据像元值范围生成掩膜文件,让背景不参与融合可以得到更高的融合结果。
标准景的资源三号全色和多光谱数据有部分区域未重叠(如下图),我们可以获取多光谱和全色的重叠区,以重叠区作为掩膜区域参与融合,最终得到重叠区的融合结果。
图1:不规则的图像
重叠区计算工具下载:http://pan.baidu.com/s/1eQlCU8y
将get_overlay_shapefile.sav文件拷贝到…\Exelis\ENVI51\classic\save_add 文件夹中。 (1)在开始菜单中启动 ENVI 5.1 classic 。资源三号数据的全色和多光谱图像都经过正射校正和自动几何配准(6度夹角校正:
http://blog.sina.com.cn/s/blog_764b1e9d0101btnb.html )。 (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)在 Selec t 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反变换得到融合图像。
文章来源:http://blog.sina.com.cn/s/blog_764b1e9d0102uz82.html
标准景的资源三号全色和多光谱数据有部分区域未重叠(如下图),我们可以获取多光谱和全色的重叠区,以重叠区作为掩膜区域参与融合,最终得到重叠区的融合结果。
图1:不规则的图像
重叠区计算工具下载:http://pan.baidu.com/s/1eQlCU8y
将get_overlay_shapefile.sav文件拷贝到…\Exelis\ENVI51\classic\save_add 文件夹中。 (1)在开始菜单中启动 ENVI 5.1 classic 。资源三号数据的全色和多光谱图像都经过正射校正和自动几何配准(6度夹角校正:
http://blog.sina.com.cn/s/blog_764b1e9d0101btnb.html )。 (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)在 Selec t 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反变换得到融合图像。
文章来源:http://blog.sina.com.cn/s/blog_764b1e9d0102uz82.html
0 个评论
相关问题
- ArcGIS如何求栅格图像的“质心”
- [一些技术]ArcGIS中对图像的配准和矢量化
- 怎样用arcgis对已经矢量化的不同时期同一地区的土地利用图像进行时空变化分析?
- landsat8 融合后无法定标?
- arcmap中如何将线图层和面图层融合或叠加在一个图层
- envi 两幅图像图像镶嵌后有明显的色差及接缝线,该怎么解决?
- 我配置arcgis data store ,请问发布服务的时候怎么还是说是资源将被复制到服务器呢?求教大佬解答一下
- c# IFeatureClass.Search超出系统资源
- 缓冲区重合部分,如何融合为一个,并计算两个融合部分相同字段值的和,或者说他们两部分重合部分的叠加影响?
- ENVI中图像分辨率 (Resolution)与比例尺 (Sacle) 之间的如何换算?
- 利用ArcGIS将tif发布切片服务,怎么将tif图像周边的“白色区域”发布为透明切片