从国际科学数据服务平台上下载了一幅2011的影像数据,结果条带很严重。经谷歌,是因为Landsat-7 ETM+机载扫描行校正器故障导致2003年5月31日之后的图像都出现数据条带丢失。
修复方法有:
A.差值修复:利用同一景影像完好的数据部分,地数据缝隙进行插值;
B.回归修复:a.利用故障前的正常数据对数据缝隙进行填充;b.或者是利用多景(2-4)不同时相的异常数据生成一景缝隙填充的数据产品。
本文将要尝试的方法:
1.国科平台在线两种方法
2.Envi 去条带补丁
3.Erdas 去条带模块
1.而国科平台上提供了两种修复方法,分别是:多影像局部自适应回归分析模型(2-b)和多影像固定窗口回归分析模型(2-b)。有在线试过第一种回归分析模型,结果显示如下:可以看到条带缺失得到很大的改善,但条带插值的部分和周围的像元有明显的区别,插值效果不太理想。
2.在ESRI社区找到一种条带修复方法:http://bbs.esrichina-bj.cn/ESRI/viewthread.php?tid=81687
修复补丁在ENVI上适用。将算法(.sav)放在..\ITT\IDL71\products\envi47\save_add(这是针对4.7版本的安装路径,4.8同样试用),重新启动ENVI后在Basic Tools → Preprocessing → Data-Specific Utilities→Landsat TM 产生一个菜单Landsat Gapfill,点击该菜单,按需要选择相应项,如选择第一项Single file gap fill,即可进行条带修复。
在这里,解释一下上图三个选项:第一个单个文件缝隙填充(三角剖分)用的是三角插值,后两个为双波段缝隙填充,区别在于匹配的时候是用全局直方图,还是用局部直方图。并且提示如果选择局部插值会花费更长的时间。
处理了第一个波段,效果明显,显示如下:
同样,可以处理六个波段,如果影像有元数据的话,按以下方式打开即可。(File——open external file——lansat——geotiff with metadata)
另外有人写的补丁,源代码在第四页,看评论效果都说还不错,哪天试试。http://bbs.esrichina-bj.cn/ESRI/viewthread.php?tid=52263&extra=pageD4&page=1
3.另外,在Erdas的辐射增强那一部分,有一个模块就是“去条带处理”,英文是Destripe TM Data。
尝试了一下,可以看一下结果:
几乎没什么变化,太让我失望了。处理过程还慢到不行,花了半个多小时,这期间还什么都不能做,很占内存和磁盘空间。