图像分幅裁剪和拼接.doc

图像分幅裁剪（Subset Image）在实际工作中，经常需要根据研究工作范围对图像进行分幅裁剪，按照ERDAS实现图像分幅裁剪的过程，可以将图像分幅裁剪分为两种类型：规则分幅和不规则分幅。1.规则分幅裁剪（Rectangle Subset Image）规则分幅是指裁剪图像的边界范围是一个矩形，通过左上角和右下角两点的坐标，就可以确定图像的裁剪位置，整个裁剪过程比较简单。ERDAS 图标面板菜单条：Main Data Preparation Data Preparation菜单（图4.1）选择Subset Image开Subset Image对话框（图4.8）ERDAS 图标面板工具条：点击Data Prep图标打开 Data Preparation菜单（图4.1）选择Subset Image打开Subset Image对话框（图4.8）图4.8 Subset Image对话框在Subset Image对话框中需要设置下列参数：输入文件名称（Input File）：Lanier.img输出文件名称（Output File）：Lanier_sub.img坐标类型（Coordinate Type）：File裁剪范围（Subset Definition）：输入ULX、ULY、LRX、LRY输出数据类型（Output Data Type）：Unsigned 8 bit输出统计忽略零值：Ignore Zero in Output Stats输出像元波段（Select Layers）：1:7（表示选择从第一波段到第七波段）4,3,2（表示选择4，3，2四个波段）OK（关闭Subset Image对话框，执行图像裁剪）说明：在上述图像裁剪过程中，裁剪范围是通过直接输入左上角坐标和右下角坐标定义的。此外，还可以通过两种方式定义裁剪范围：其一是应用查询框（Inquire Box），然后在Subset Image对话框中选择From Inquire Box功能；其二是应用AOI，然后在Subset Image对话框中选择AOI功能，打开AOI对话框，并确定AOI区域来自图像视窗即可。输出波段的选择：一种是用冒号表示从第几波段到第几波段，另外一种是用逗号隔开表示仅仅选择这几个波段。2.不规则分幅裁剪（Polygon Subset Image）不规则分幅裁剪是指 裁剪图像的边界范围是个任意多边形，无法通过左上角和右下角两点的坐标确定图像的裁剪位置，而必须事先生成一个完整闭合多边形区域，可以是一个AOI多边形，也可以是ArcInfo的一个Polygon Coverage。针对不同的情况采用不同裁剪过程。AOI多边形裁剪（Polygon AOI Subset Image）首先在视窗中打开需要裁剪的图像，并应用AOI工具绘制多边形AOI，可以将多边形AOI保存在文件中（*.aoi），也可以暂时不退出视窗，将图像与AOI多边形保留在视窗中，然后：ERDAS 图标面板菜单条：Main Data Preparation Data Preparation菜单（图4.1）选择Subset Image开Subset Image对话框（图4.8）ERDAS 图标面板工具条：点击Data Prep图标打开 Data Preparation菜单（图4.1）选择Subset Image打开Subset Image对话框（图4.8）在Subset Image对话框中需要设置下列参数：输入文件名称（Input File）：Lanier.img输出文件名称（Output File）：Lanier_sub.img应用AOI确定裁剪范围：点击AOI按钮打开选择AOI（Choose AOI）对话框（图4.9）图4.9 Choose AOI对话框在Choose AOI对话框中确定AOI的来源（AOI Source）：File或Viewer输出数据类型（Output Data Type）：Unsigned 8 bit输出统计忽略零值：Ignore Zero in Output Stats输出像元波段（Select Layers）：1:7（表示选择从第一波段到第七波段）4,3,2（表示选择4，3，2四个波段）OK（关闭Subset Image对话框，执行图像裁剪）ArcInfo多边形裁剪（Polygon Coverage Subset Image）如果按照行政区划边界或自然区划边界进行图像的分幅裁剪，往往是首先利用ArcInfo或ERDAS的Vector模块绘制精确的边界多边形（Polygon），然后以ArcInfo的Polygon为边界条件进行图像裁剪。对于这种情况，需要调用ERDAS其他模块的操作分两步完成。第一步：将ArcInfo多边形转换成栅格图像文件ERDAS 图标面板菜单条：Main Image Interpreter Utilities Vector to Raster 打开Vector to Raster对话框（图4.10）ERDAS 图标面板工具条：Main Vector Vector to Raster 打开Vector to Raster对话框（图4.10）ERDAS 图标面板工具条：点击Interpreter图标 Utilities Vector to Raster 打开Vector to Raster对话框（图4.10）ERDAS 图标面板工具条：点击IVectot图标 Vector to Raster 打开Vector to Raster对话框（图4.10）图4.10 Vector to Raster对话框（a）图4.10 Vector to Raster对话框（b）在Vector to Raster对话框中需要设置下列参数：确定矢量文件名称（Input Vector File）：zone88确定输出文件名称（Output Raster File）：raster.img确定矢量文件类型（Vector Type）：Polygon转换范围大小（Subset Definition）：ULX、ULY、LRX、LRY输出像元大小（Pixel Size）：30栅格数据类型（Output Pixel Data Type）：Unsigned 8 bit使用矢量属性值（An Item as Pixel Value）：ZONE88-ID栅格文件类型（Output Image Type）：Thematic是否创建图层金字塔（Create Pyramid Layers）输出统计忽略零值：Ignore Zero in Output StatsOK（关闭Vector to Raster对话框，执行失栅转换）第二步：通过掩模运算（Mask）实现图像不规则裁剪ERDAS 图标面板菜单条：Main Image Interpreter Utilities Mask 打开Mask对话框（图4.11）ERDAS 图标面板工具条：点击Interpreter图标 Utilities Maskr 打开Mask对话框（图4.11）图4.11 Mask对话框在Mask对话框中需要设置下列参数：确定输入文件名称（Input Vector File）：Lanier.img确定掩模文件名称（Input Mask File）：raster.img点击Set Recode设置裁剪区域内新值（New Value）为1，区域外取0值确定掩模区域做交集运算：Intersection输出图像文件名称（Output File）：mask.img输出数据类型（Output Data Type）：Unsigned 8 bitOk（关闭Mask对话框，执行掩模运算）实现图像不规则裁剪。图像拼接处理（Mosaic Image）在遥感图像的应用中，当研究区处于几幅图像的交界处或研究区较大需要多幅图像才能图像拼接处理是要将就有地理参考的若干相邻图像合并成一幅图像或一组图像，需要拼接的输入图像必须含有地图投影信息，或者说输入图像必须经过几何校正处理或进行过校正标定。虽然所有的输入图像可以具有不同的投影类型、不同的像元大小，但必须具有相同的波段数。在进行图像拼接时，需要确定一幅参考图像，参考图像将作为输出拼接图像的基准，决定拼接图像的对比度匹配以及输出图像的地图投影、像元大小和数据类型。1.图像拼接功能（Image Mosaic Function）图像拼接处理工具可以通过下列两种途径启动：ERDAS图标面板菜单条：Main Data Preparation Data Preparation菜单选择Mosaic Image打开Mosaic Tool视窗（图4.30）ERDAS图标面板工具条：点击Data Prep图标 打开Data Preparation菜单选择Mosaic Image打开Mosaic Tool视窗（图4.30）图4.30 Mosaic Tool视窗Mosaic Tool视窗由菜单条、工具条、图形窗口和状态条四部分组成，其中菜单条中的菜单命令及其功能如表4.6所列，工具条中的图标及其功能如表4.7所列。表4.6 Mosaic Tool视窗菜单条命令及其功能命令功能FileNewOpenSaveSave AsAnnotationCloseClose All文件操作打开新的图像拼接工具打开图像拼接工程文件（*.mos）保存图像拼接过程文件（*.mos）重新保存图像拼接工程文件将拼接图像轮廓保存为注记文件关闭当前图像拼接工具关闭所有图像拼接工具Edit：Image ListAdd ImageDelete ImageImage MatchingIntersection ListSet Overlay FunctionOutput ListOutput OptionsDelete Outputs编辑操作：显示输入图像列表向图像拼接视窗加载图像删除图像拼接工程中的图像设置拼接图像的匹配参数设置拼接图像的拼接关系设置拼接图像的叠加关系显示输出图像列表设置输出图像参数输出输出图像Process：Run Mosaic处理操作：执行图像拼接处理表4.7 Mosaic Tool视窗工具图标及其功能图标命令功能Add Image向图像拼接视窗加载图像Set Input Mode：设置输入图像模式Image Resample打开图像重采样对话框Image Matching打开图像匹配对话框Send Image to Top将选择图像置于最上层Send Image Up One将选择图像上移一层Send Image to Buttom将选择图像置于最下层Send Image Down One 将选择图像下移一层Reverse Image Order将选择图像次序颠倒Set Intersection Mode：设置图像交接关系Next Intersection选择下一种相交方式Previouse Intersection选择前一种相交方式Overlap Function打开叠加功能对话框Default Cutlines设置缺省相交截切线AOI Cutlines设置AOI区域截切线Toggle开关截切线的应用模式Delete Cutlines删除相交区域截切线Cutlines Selection Viewer：打开截切线选择视窗Auto Cutlines Mode设置截切线自动摸索Set Output Mode：设置输出图像模式：Output Options打开输出图像设置对话框Reset Canvas改变图面尺寸以适应拼接图像Scale Canvas改变图像比例以适应选择对象Select Point选择一个点进行查询Select Area选择一个区域进行查询Zoom Image In By 2两倍放大图形窗口Zoom Image Out By 2两倍缩小图形窗口Select Area for Zoom选择一个区域进行放大Roam Canvas图形窗口漫游Image List拼接图像列表2.卫星图像拼接（Satellite Image Mosaic）本节通过三幅陆地资源卫星图像（wasia1_mss.img、wasia2_mss.img、wasia3_mss.img）的拼接处理，介绍卫星图像的拼接过程。启动图像拼接工具（Start Mosaic Tools）图像拼接工具可以通过下列两种途径启动：ERDAS图标面板菜单条：Main Data Preparation Data Preparation菜单选择Mosaic Image打开Mosaic Tool视窗（图4.30）ERDAS图标面板工具条：点击Data Prep图标 打开Data Preparation菜单选择Mosaic Image打开Mosaic Tool视窗（图4.30）加载Mosaic图像（Add Images for Mosaic）Mosaic Tool视窗菜单条：Edit Add Images打开Add Images for Mosaic对话框（4.31）Mosaic Tool视窗工具条：点击Add Image图标打开Add Images for Mosaic对话框（4.31）图4.31 Add Images for Mosaic对话框在Add Images for Mosaic对话框中，需要设置以下参数：拼接图像文件（Image File Name）：wasia1_mss.img图像拼接区域（Image Area Option）：Compute Active Area（Edge）Add（图像wasia1_mss.img被加载到Mosaic视窗中）同样的过程加载wasia2_mss.img和wasia3_mss.imgClose（关闭Add Images for Mosaic）3.图像叠置组合（Images Stacking）Mosaic Tool视窗工具条：点击Set Input Mode图像，进入设置输入图像模式的状态，Mosaic Tool视窗工具条会出现与该模式对应的调整图像叠置次序的编辑图标（如表4.7所列），充分利用系统所提供的编辑工具，提供需要进行上下层调整，这些调整工具包括：Send Image to Top：将选择图像置于最上层Send Image Up One：将选择图像上移一层Send Image to Buttom：将选择图像置于最下层Send Image Down One：将选择图像下移一层Reverse Image Order：将选择图像次序颠倒调整完成后，在Mosaic Tool视窗图形窗口点击一下，退出图像叠置组合状态。4.图像匹配设置（Image Matching）Mosaic Tool视窗菜单条：EditImage Matching打开Matching Options对话框（图4.32）Mosaic Tool视窗工具条：点击Set Input Mode图标，进入设置输入图像模式（图4.32）点击Image Matching图标打开Matching Options对话框（图4.32）图4.32 Matching Options对话框图4.33Set Overlap Function对话框在Matching Options对话框中做如下设置：设置匹配方法（Matching Method）Overlap Area（重