遥感实验五_影像镶嵌、裁剪、融合

在Mosaic Tool视窗菜单条中，点击Edit/setOverlap Function—打开setOverlap Function对话框，如图1.2.6示，设置以下参数：
.设置相交关系(Intersection Method):No Cutline Exists。
.设置重叠图像元灰度计算(select Function):Average。
图2.1.5
点击DataPrep,在弹出的下拉菜单中单击Subset Images，在Input File中输入裁切的底图xianqiang.img，在Output File中设置输出文件路径和文件名，这里保存名为jianqie3.img。
单击From Inquire Box,然后点击AOI，在弹出的Choose AOI中点击Viewer,点击OK。，最后在subset点击OK,步骤如图2.1.6示。
图2.1.6
图2.1.7
在新视图窗口中打开裁切结果，如图2.1.8示。
图2.1.8
同理对全色影像进行剪切。
操作步骤如图2.2.1—2.2.3示。
图2.1.1
图2.2.2
图2.2.3
全色影像裁切效果如图2.2.4示。
图2.2.4
2.3.按已有图像范围裁切(掩膜)
按已有图像的范围从一幅较大图像中裁切一部分图像时，按下图所示方法操作：其中4处为较大图像文件（即待裁切图像），5处为限定范围的图像文件（即裁切范围），6处为结果文件（即裁切后图像），如图2.3.1示。
.Apply—close。
图像拼接线设置，在Mosaic Tool视窗菜单条中选择Set Mode For Intersection按钮 ，两幅图像之间将出现叠加线，单击两幅图像的相交区域，重叠区域将被高亮显示。根据实际需要，选择拼接线模式：
.
：自动生成重叠区域的拼接线
. ：手动绘制拼接线，可以是矢量文件，或AOI文件
图3.1.2
图3.1.3
计算机自动执行merge命令。在新的视图窗口中打开融合图像resolutionj.img，如下图3.1.4示。
图3.1.4
注：上图是4、3、2波段多光谱影像与全色影像融合后影像。
4、对比分析
分别打开仙游县全色、多光谱（4、3、2波段）和融合后影像，通过对局部地物信息对比，可以看出融合后影像具备全色高空间分辨和多光谱信息量丰富，可视性强的特点，如图4.1.1—4.1.2示。
图1.3.5
点击OK后计算机自动执行Mosaic命令，如图1.3.6示。
图1.3.6
打开一新视图窗口，显示镶嵌后全色影像xiangqian2,如图1.3.7示。
图1.3.7
注意：进行镶嵌的图像必须具有相同的投影系统，否则见无法实现镶嵌拼接。
2、图像裁切。
2.1.规则裁切
点击DataPrep,在弹出的下拉菜单中单击Subset Images，在Input File中选择经过镶嵌处理过的数据xianqiang.img。在Output File中设置保存路径和文件名，文件名位jianqie2.img，如图2.1.1。
1.1.投影定义和转换
在ERDAS中，点击DataPrep,在下拉选项卡中点击Rejection Images,在Input File中输入需要进行投影转换的影像数据——福建某地区2000年5月30米分辨率的的多光谱影像（本例以TM4、3、2波段为例）。在Output File设置保存路径和输出文件名。在Categories中点击右侧的小地球标志进行投影定义。投影参数设置如下图1.1示，点击OK，完成投影转换。本图及以下各图均将WGS-84投影转换成Gauss Kruger投影。
图1.2.2
加载拼接图像的同时，切换到Image Area Options选项卡，设置拼接影像范围选择方法为：Use Entire Image。
图1.2.3
图1.2.4
图像匹配设置，在Mosaic Tool工具条选择Display Color Correction按钮 ，或在Mosaic Tool视窗菜单条中点击Edit /Color Corrections—打击Color Corrections对话框，根据需要设置颜色匹配方法，如图1.2.5示。
投影变换结果如下图所示。
图1.8.1
图1.8.2
图1.8.3
图1.8.4
说明：在ERDAS中进行投影转换时，由于待定义或转换投影区域与所选的投影类别不在同一半球（南北半球）时，会产生“倒像”坐标值会显示为“负值”。因此在定义投影时要根据遥感影像的实际区域确定投影种类（categories）,比如说本次实验中若选择阿根廷、澳大利亚等南半球国家投影种类，则将显示一幅“倒像”。
遥感实验五数字图像处理
------------图像镶嵌、裁切及融合
一、实验目的
学会图像镶嵌、图像裁切及图像融合等技术，通过实际影像的操作，制作可用于实际工作的某区域遥感图像，为下一次实验准备数据。
二、实验数据
某区域的遥感图像：11942E20000504.rar、11943E20010304.rar；
图1.2.5
注意：如果输入的拼接影像自身存在较大的亮度差异（例如，中间暗周围亮或者一边亮一边暗），需要首先利用匀光处理（Use Image Dodging）或色彩平衡（Use Color Balancing）去除单幅影像自身的亮度差异；如果输入的影像间直方图差异比较大，可通过直方图匹配（Use Histogram Matching）执行影像的色彩调整。
1.2.影像镶嵌
打开ERDAS软件，点击DataPrep，在下拉选项中点击Mosaic Images，点击Mosaic Tool，在弹出的窗口中点击 ，加载需要进行镶嵌的影像数据。本例中加载经过投影转换的多光谱影像数据1.img和2.img。
图1.2.1
分别加载1.img、2.img影像，如图1.2.2—1.2.4所示。
图1.5
经过arcgis投影转换过的全色影像如图1.6示，图像属性数据显示投影坐标系统为Gauss Kruger,Datum为D_Krasovsky_1940。
图1.6
同上述在erdas中投影转换的方法对全色图像p119r043_7×20010304p_sub.Img进行投影转换，如图1.7。
图1.7
图1.3
打开arcgis,在工具箱中找到define projection工具，输入需投影转换的影像数据tm_11942_8,设置投影类型，选择modify选项，在图左侧点击selece,在弹出的窗口中选择Krasovskv_1940.prj，点击add,完成投影转换，步骤详见图1.4、1.5。
图1.4
图1.1
同理，对裁切的多光谱小图进行投影转换，原理及步骤亦同上，图1.2示。
图1.2
对全色波段影像数据tm11942_8进行投影变换，原理同多光谱影像投影变换，但在erdas中进行投影转化时由于在选择categories时，选择了南半球国家投影类别发生错误，结果显示为一“倒像”，故tm11942_8影像采用ArcGIS软件进行投影转换，转换目的主要是讲投影信息中的Datum转成Krasovskv。
15米全色影像
融合后ቤተ መጻሕፍቲ ባይዱ像
30米假彩色影像（4、3、2波段）
图4.1.1
以下以多光谱5、4、3波段与全色影像融合为例进行融合操作。操作步骤如下，这里我们只需在mask.img视图窗口中点击Raster,在下拉选项中点击band combinations,进行显示波段调整，这里设为5、4、3波段。
图4.1.2
图2.1.2
图2.1.3
计算机自动执行裁切命令，在新视图窗口中打开裁切图像，如图2.1.4示。
图2.1.4
2.2.利用AOI裁切
新建窗口，打开仙游县的矢量图层xianyou.ship.利用xianyou.ship建立AOI感兴趣区域。
在镶嵌后的图像上点击AOI，选择AOI Tool,弹出AOI工具条，在仙游县矢量图层上单击鼠标，使图像以黄色显示。然后点击 ，生成AOI感兴趣区域，以虚线显示在xiangqian.img图像上，过程如图2.1.5示。
. ：重叠区域内拼接线两侧区域选择切换
. ：无拼接线
注意：该步骤中如果选择了有拼接线的拼接模式，则上一步中重叠区域像元灰度计算方式自动失效，拼接后影像像元灰度值直接采用原拼接图像的像元灰度值。
图1.2.6
点击process,点击Run mosaic,执行镶嵌操作，设置保存文件名位mosaic1，点击OK。
图4.1.5
图4.1.6
说明：如上图所示，融合后影像兼有全色影像和多光谱影像的优点。如，在表现河流下游江心洲及河汊信息时，融合后影像能够分辨河道细节信息诸如河汊和江心洲，比较多光谱
（5、4、3）影像，明显显示出高分辨率特点。
图2.3.1
图2.3.2
在新视图窗口中打开mask后的影像，效果图如图2.3.3示。
图2.3.3
同理对全色影像进行mask处理，操作步骤如图2.3.4示。
图2.3.4
在新视图窗口中打开mask2影像，如图2.3.5示。
图2.3.5
3、图像融合
图像融合是在裁切的基础上，对仙游县全色和多光谱影像进行融合复合，以达到图像增强的目的，方便目视解译。
某区域的范围：xianyou.shp
三、实验内容及主要步骤
1、图像镶嵌
注：要镶嵌的两幅或多幅影像要求具有相同的投影信息，如果不同，则需要首先统一。ERDAS IMAGINE中提供了投影转换的工具，点击 、选择Reproject Images；或者，也可以在ArcGIS的ArcToolbox中选择Projections and Transformations/Raster/Project Raster进行转换。以下以ERDAS IMAGINE软件为例进行投影转换。