遥感影像镶嵌在不同软件上的比较——Envi, Erdas, Qmasaic

老板最近给了一个任务，要用到广东珠三角地区2005年，2010年，2015年的土地利用情况，所以笔者开始下载影像进行图像镶嵌，然后分类。由于镶嵌花了笔者大量的时间，所以把镶嵌过程记录下来，以便后来的人不用再花那么多时间。

一、数据准备

1.使用FIleZilla Client客户端下载地理空间数据云上的遥感影像

2.对于2005年，2010年，首先下载了Landsat7ETM+数据，再对数据使用landsat_gapfill去条带工具进行去条带处理。去完条带之后，发现效果不好，条带里损失了很多信息，这样分类精度会大大降低，所以重新下载了2005年，2010年的Landsat TM5影像。注意下载的需要拼接的影像时间间隔一定要在三个月内，最好是同一天的影像，不然在镶嵌过程中会有很大的色差。在选择影像下载时间的时候，要先看那年有没有清晰度合适的可下载影像。2010年珠三角区域应该是没有合适影像的，所以，下载之前一定要先搞清楚。

接着对影像进行辐射定标，大气校正。如果只做分类，那可以不用大气校正，但必须辐射定标，以产生dat数据。

二、影像拼接。

2.1使用Envi进行拼接

第一步.导入数据

第二步.准备工作

如果使用seamless Mosaic进行整体拼接，效果非常差，不能用于分类。

2.2使用Qmosaic进行拼接

第一步.导入数据

第二步.使用羽化功能去黑框

第三步.使用匀色功能进行处理，Qmosaic有几种匀色方式，“色彩均衡”、“直方图匹配”。如果进行整体匹配，效果不好，最好分幅进行匹配，Qmosaic拼接的速度比Envi快很多，也比Erdas要快，但快的不明显。

第四步.生成镶嵌线，按照默认来就可以

第五步.输出结果，可以用AOI作为掩膜。笔者首先使用Qmosaic6.1.2版本，镶嵌需要购买，如果不购买，输出结果会有很多十字，大大影响效果。

所以找到Qmosaic5.7.5进行重新进行镶嵌。

要注意的是，在批处理前要找到submit里的Simultaneous Processes设置并行运算处理的影像数目。

第三步.镶嵌

打开MosaicPro工具，将需要镶嵌的影像导入。

提醒：一开始导入的时候影像不会显示出来。要打开view选项中的show raster才能够使影像显示出来，如果没打开，即使下面vis.选中了，遥感影像也不会显示。

根据需要选择镶嵌的操作，包括重采样、调整影像波长值、色彩校正、叠加和接边线操作。

提醒：set overlap/seam，line function选项原本有四个在设置接边线后会变成只有两个选项。

最后进行镶嵌到磁盘操作。

因为没有事先将海面裁去，所以使用Erdas镶嵌效果仍然不太好，最后还是使用Envi生成的结果。当然镶嵌速度要比Envi快得多。

第四步.裁剪

为了搞清楚在Erdas中如何裁剪也花了一番功夫

首先将shp文件导入，然后需要在vector选项卡中的edit中选择Paste from seleted object，这才能shp文件转化成AOI，接着将AOI文件另存为，然后才可以用Subset&Clip工具进行裁剪。在进行裁剪前，一定要核对shp文件和栅格文件是不是同一坐标系，如果不同，则裁剪会失败。

对于掩膜提取，要保证掩膜和影像图层数相同，如果不同会 错。

三、总结

对Erdas、ENVI、Qmosaic三个软件的镶嵌功能进行比较花了笔者近一个月时间，最后得出的结论是三个软件各有优点，但针对镶嵌而言，ENVI > Qmosaic > Erdas。

希望笔者介绍的经验能帮助到后来的GISer。