图像分幅裁剪与拼接

《图像分幅裁剪与拼接》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、图像分幅裁剪（SubsetImage）在实际工作中，经常需要根据研究工作范围对图像进行分幅裁剪，按照ERDAS实现图像分幅裁剪的过程，可以将图像分幅裁剪分为两种类型：规则分幅和不规则分幅。1.规则分幅裁剪（RectangleSubsetImage）规则分幅是指裁剪图像的边界范围是一个矩形，通过左上角和右下角两点的坐标，就可以确定图像的裁剪位置，整个裁剪过程比较简单。ERDAS图标面板菜单条：Main→DataPreparation→DataPreparation菜单（图4.1）→选择SubsetImage→开SubsetImage对话框（图4.8）ERDAS图标面板工

2、具条：点击DataPrep图标→打开DataPreparation菜单（图4.1）→选择SubsetImage→打开SubsetImage对话框（图4.8）图4.8SubsetImage对话框在SubsetImage对话框中需要设置下列参数：→输入文件名称（InputFile）：Lanier.img→输出文件名称（OutputFile）：Lanier_sub.img→坐标类型（CoordinateType）：File→裁剪范围（SubsetDefinition）：输入ULX、ULY、LRX、LRY→输出数据类型（OutputDataType）：Unsigned8bit

3、→输出统计忽略零值：IgnoreZeroinOutputStats→输出像元波段（SelectLayers）：1:7（表示选择从第一波段到第七波段）4,3,2（表示选择4，3，2四个波段）→OK（关闭SubsetImage对话框，执行图像裁剪）说明：①在上述图像裁剪过程中，裁剪范围是通过直接输入左上角坐标和右下角坐标定义的。此外，还可以通过两种方式定义裁剪范围：其一是应用查询框（InquireBox），然后在SubsetImage对话框中选择FromInquireBox功能；其二是应用AOI，然后在SubsetImage对话框中选择AOI功能，打开AOI对话框，并确定

4、AOI区域来自图像视窗即可。②输出波段的选择：一种是用冒号表示从第几波段到第几波段，另外一种是用逗号隔开表示仅仅选择这几个波段。2.不规则分幅裁剪（PolygonSubsetImage）不规则分幅裁剪是指裁剪图像的边界范围是个任意多边形，无法通过左上角和右下角两点的坐标确定图像的裁剪位置，而必须事先生成一个完整闭合多边形区域，可以是一个AOI多边形，也可以是ArcInfo的一个PolygonCoverage。针对不同的情况采用不同裁剪过程。⑴AOI多边形裁剪（PolygonAOISubsetImage）首先在视窗中打开需要裁剪的图像，并应用AOI工具绘制多边形AOI，

5、可以将多边形AOI保存在文件中（*.aoi），也可以暂时不退出视窗，将图像与AOI多边形保留在视窗中，然后：ERDAS图标面板菜单条：Main→DataPreparation→DataPreparation菜单（图4.1）→选择SubsetImage→开SubsetImage对话框（图4.8）ERDAS图标面板工具条：点击DataPrep图标→打开DataPreparation菜单（图4.1）→选择SubsetImage→打开SubsetImage对话框（图4.8）在SubsetImage对话框中需要设置下列参数：→输入文件名称（InputFile）：Lanier.i

6、mg→输出文件名称（OutputFile）：Lanier_sub.img→应用AOI确定裁剪范围：点击AOI按钮→打开选择AOI（ChooseAOI）对话框（图4.9）图4.9ChooseAOI对话框→在ChooseAOI对话框中确定AOI的来源（AOISource）：File或Viewer→输出数据类型（OutputDataType）：Unsigned8bit→输出统计忽略零值：IgnoreZeroinOutputStats→输出像元波段（SelectLayers）：1:7（表示选择从第一波段到第七波段）4,3,2（表示选择4，3，2四个波段）→OK（关闭Subse

7、tImage对话框，执行图像裁剪）⑵ArcInfo多边形裁剪（PolygonCoverageSubsetImage）如果按照行政区划边界或自然区划边界进行图像的分幅裁剪，往往是首先利用ArcInfo或ERDAS的Vector模块绘制精确的边界多边形（Polygon），然后以ArcInfo的Polygon为边界条件进行图像裁剪。对于这种情况，需要调用ERDAS其他模块的操作分两步完成。第一步：将ArcInfo多边形转换成栅格图像文件ERDAS图标面板菜单条：Main→ImageInterpreter→Utilities→VectortoRaster→打