打开Erdas2022软件,依次点击【Raster】→【Subset & Chip】面板,启动【创建子集图像】工具,erdas提供的提取工具如下:
下面主要结合实际案例,讲述常用的【Create Subset Image】和【Mask】工具的使用方法。
一、矩形裁剪
加载被裁剪的多波段卫星影像(数据位于实验数据包erdas005.rar中,订阅专栏后从私信查收完整实验数据包及软件)如下:
矩形裁剪工具的原理是基于两点或者四点坐标,构建一个矩形区域,进而用区域进行裁剪。
1. 指定坐标
主要详细参数设置如下:
Input File(输入文件名):太和镇卫星影像.tif;
Output File(输出文件名):sub5.img,可以指定img格式或者tif等多种栅格格式;
Coodinate Type(坐标类型):File或者Map;
Subset Definition(裁剪范围):有 Two Corners(两角)和 Four Corners (四角)两种方式,这里选择前者,输入 ULX(左上角 X 坐标):563059.3199、ULY(左上角 Y 坐标):3324787.5699 和 LRX(右下角 X 坐标):563268.2528、LRY(右下角 y 坐标):3324646.5085;
Output Data Type(输出数据类型):Unsigned 8 bit,可以与原数据保持一致,通常可以提高位数;
Continuous;(输出波段):1:7(表示从 1 到 7 波段)
参数设置如下图所示:
点击OK,完成矩形区域裁剪,结果如下:
2. 查询框
矩形区域的确定,除了有指定的坐标,也可以从查询框获取。点击Inquire Box,弹出坐标窗口,同时可以直接拖动范围大小,如下。
在确定了查询框大小后,可以在Subset窗口直接点击【From Inquire Box】,点击之后,坐标值会自动读取。点击OK即可完成裁剪。
查询框裁剪结果:
注意:使用坐标时,不能再点击使用查询框,否则坐标值会替换为查询框的坐标值。
二、多边形裁剪
1. 多边形shp
通常情况下,使用已有矢量(可以是shp,dxf,kml等格式)范围对影像进行裁剪。
(1)shp转raster
加载提供的shp格式的矢量数据。
默认显示样式。
你可以在Style面板下对矢量数据符号化:
在Erdas中,不能直接使用矢量范围对影像进行裁剪,需要将矢量转为栅格,再用栅格进行掩膜提取。点击【Manage Data】面板下的【Vector to Raster】工具,设置如下参数。
转换后的栅格:
可以看到,在erdas中,非规则的矢量转换为栅格后,栅格为规则的矩形,显然不能满足非规则形状提取的要求。此时,可以使用ArcGIS软件进行矢栅转换。
其中,一个很重要的参数是【值字段】,此时选择值为1的Id字段。如果shp数据id的值为0,需要编辑为1。
ArcGIS转换后的栅格数据,其栅格值为1。
erdas加载栅格掩膜效果:
可以使用Inquire工具查询栅格值,从下面的属性表中可以看出值为1。
接下来就可以使用erdas中的mask工具进行栅格掩膜裁剪了。
(2)mask
启动【Mask】工具。
Input Mask File中选择转换得到的栅格掩膜文件bound.tif。
掩膜裁剪结果:
2. AOI裁剪
Erdas中提供了方便灵活的AOI(Area Of Intrest,感兴趣区)工具,关于AOI的创建与保存,请点击专栏相应的文章进行专项学习,本案例中的AOI如下图:选择掩膜数据。
启动【Create Subset Image】工具,重点是点击左下角的AOI工具,在打开的【Choose AOI】窗口中,选择对应的Viewer中的AOI或者本地AOI文件即可。
点击OK,执行裁剪,AOI裁剪结果:
《Erdas遥感图像处理实验教程》专栏介绍:
以最新版erdas2022为主,讲述遥感图像显示与数据输入/输出、遥感图像的几何校正、遥感图像的增强、遥感图像分类、图像分析、矢量处理、专题制图、空间建模、处理国产影像等。订阅专栏,提供配套实验数据及erdas软件等资料,持续同步更新!