当前位置:
文档之家› 遥感实验五_影像镶嵌、裁剪、融合
遥感实验五_影像镶嵌、裁剪、融合
在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软件为例进行投影转换。