遥感数据图像处理实验三、遥感图像的几何校正与裁剪

  • 格式:pdf
  • 大小:4.65 MB
  • 文档页数:62

实验三、遥感图像的几何校正与裁剪实验内容:1.图像分幅裁剪(Subset Image)2.图像几何校正(Geometric Correction)3.图像拼接处理(Mosaic Imgaes)4.生成三维地形表面(3D Surfacing)1.图像分幅裁剪在实际工作中,经常需要根据研究工作范围对图像进行分幅裁剪,按照ERDAS IMAGINE 8.4实现图像分幅裁剪的过程,可以将图像分幅裁剪为两类型:规则分幅裁剪,不规则分幅裁剪。

1.1规则分幅裁剪(以c:\Program File\ IMAGINE 8.4\examples\lanier.img为例)规则分幅裁剪是指裁剪图像的范围是一个矩形,通过左上角和右上角两点的坐标可以确定图像的裁剪位置,过程如下:方法一:→ERDAS IMAGINE 8.4 图标面板菜单条:Main→Data Preparation(或单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标)→打开Data Preparation 对话框→单击Subset Image按钮,打开Subset对话框在Subset对话框中需要设置下列参数:→输入文件名(Input File):lanier.img→输出文件名(Output File):lanier_sub.img→坐标类型(Coordinate Type):Map→裁剪范围(Subset Definition):ULX、ULY、LRX、LRY(注:ULX,ULY是指左上角的坐标,LRX,LRY是指右上角的坐标,缺省状态为整个图像范围)→输出数据类型(Output Data Type):Unsigned 8 Bit→输出文件类型(Output Layer Type):Continuous→输出统计忽略零值:Ignore Zero In Output Stats→输出像元波段(Select Layers):2,3,4→OK(关闭Subset对话框,执行图像裁剪)方法二:→ERDAS IMAGINE 8.4图标面板菜单条:Main→Start IMAGINE Viewer(或单击RDAS IMAGINE 8.4图标面板工具条“Viewer”图标)→打开一个二维视窗→单击视窗工具条最左端的“打开文件”图标→打开Select Layer To Add对话框在Select Layer To Add对话框完成以下设置:→Look In:examples→File Name:lanier.img→Files of type:IMAGINE Image→双击OK按钮→在二维视窗中打开lanier.img文件→单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标 →打开Data Preparation对话框→单击Subset Image按钮→打开Subset对话框在Subset对话框中需要设置下列参数:→输入文件名(Input File):lanier.img→输出文件名(Output File):lanier_sub.img→坐标类型(Coordinate Type):Map→输出数据类型(Output Data Type):Unsigned 8 Bit →输出文件类型(Output Layer Type):Continuous→输出统计忽略零值:Ignore Zero In Output Stats→输出像元波段(Select Layers):2,3,4→单击From Inquire Box按钮→打开Invalid Coordinate Type对话框→单击Continue→在显示图像文件lanier.img视窗中单击工具条的“+”按钮,打开Inquire Cursor 对话框,在视窗中移动十字光标,确定裁剪范围左上角和右下角,读取其坐标分别填入Subset Image对话框的ULX,ULY中和LRX,LRY中→单击OK按钮(关闭Subset对话框,执行图像裁剪)方法三:首先在视窗中打开lanier.img文件→AOI→Tools打开AOI工具面板→单击矩形框确定裁剪范围→File→Save→AOI Layer As→打开Save AOI As对话框,输入文件名:2→单击OK(退出Save AOI As对话框)→单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标 →打开Data Preparation对话框→单击Subset Image按钮→打开Subset对话框在Subset对话框中需要设置下列参数:→输入文件名(Input File):lanier.img→输出文件名(Output File):lanier_sub.img→坐标类型(Coordinate Type):Map→输出数据类型(Output Data Type):Unsigned 8 Bit→输出文件类型(Output Layer Type):Continuous →输出统计忽略零值:Ignore Zero In Output Stats →输出像元波段(Select Layers):2,3,4→单击AOI按钮→打开Choose AOI对话框→在Choose AOI对话框作如下设置:→AOI Source:File→AOI File:2→单击OK(退出Choose AOI对话框)→单击OK(退出Subset对话框,执行图像裁剪)→单击OK(退出Modeler对话框,完成图像裁剪)1.2不规则分幅裁剪不规则分幅裁剪是指裁剪图像的边界范围是个任意多边形,无法通过左上角和右下角两点的坐标确定图像的裁剪位置,而必须事先生成一个完整的闭合多边形区域,可以是一个AOI多边形,也可以是ArcInfo的一个Polygon Coverage,针对不同的情况采用不同的裁剪过程。

(一) AOI多边形裁剪(以c:\Program File\ IMAGINE 8.4\examples\lanier.img为例)→在视窗中打开需要裁剪的图像→应用AOI工具绘制多边形可以将多边形AOI文件保存在文件中(*.aoi),也可以暂时不退出视窗,将图像与AOI 多边形保留在视窗中,然后进行以下操作:→ERDAS IMAGINE 8.4 图标面板菜单条:Main→Data Preparation(或单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标),打开Data Preparation 对话框→单击Subset Image按钮,打开Subset对话框在Subset对话框中需要设置下列参数:→输入文件名(Input File):lanier.img→输出文件名(Output File):lanier_sub.img →Coordinate Type:Map→Output:Unsigned 8 bit→Oupput:Continuous→选定Ignore Zero in Output Stats→Select Layers:2,3,4→单击AOI按钮,打开Choose AOI对话框→选定File→选择文件:2.aoi→单击OK(退出Choose AOI对话框)→单击OK(退出Subset对话框,执行图像裁剪)→单击Modeler对话框中的OK按钮完成图像裁剪(二) ArcInfo多边形裁剪(以c:\Program File\ IMAGINE 8.4\examples\gy_resample.img为例)如果是按照行政区划边界或自然区划边界进行图像的分幅裁剪,往往是首先利用ArcInfo或者ERDAS的Vector模块绘制精确的边界多边形(Polygon),然后以ArcInfo的Polygon为边界进行图像裁剪。

对于这种情况,需要调用ERDAS其它模块的功能分两步完成。

第一步:将ArcInfo多边形转换为栅格图像文件打开Vector To Raster对话框有两种方法:方法一:ERDAS IMAGINE 8.4 图标面板菜单条:Main→Image Interpreter (或单击ERDAS IMAGINE 8.4 图标面板工具条“Interpreter”图标)→打开Image Interpreter对话框→单击Utilities→打开Utilities对话框→选择Vector To Raster→打开Vector To Raster对话框方法二:ERDAS IMAGINE 8.4 图标面板菜单条:Main→Vector (或单击ERDAS IMAGINE 8.4 图标面板工具条“Vector”图标)→打开Vector Utilities对话框→选择Vector to Raster→打开Vector to Raster对话框本例中选择前一种方法:在Vector to Raster对话框中设置下列参数:→Input Vector File(输入矢量文件名称):boundgy.shp →Vector Type(确定矢量文件类型):polygon→Use Attribute as Value(使用矢量属性值):GYBOUND_ID→Output File(输出栅格文件名称):raster.img→Data Type(栅格数据类型):Unsigned 8 bit→Layer Type(栅格文件类型):Thematic→Size Difinition(转换范围大小):ULX,ULY,LRX,LRY(缺省条件下为整个图像范围) →Units(坐标单位):Meters→Cell Size(输出像元大小):X:30/Y:30→选择Square Cell(正方形像元)→单击OK(关闭Vector To Raster对话框,执行矢量到栅格的转换)第二步:通过掩模运算(Mask)实现图像的不规则裁剪ERDAS IMAGINE 8.4 图标面板菜单条:Main→Image Interpreter (或单击ERDAS IMAGINE 8.4 图标面板工具条“Interpreter”图标)→Image Interpreter→Utilities→Mask→打开Mask对话框在Mask对话框中设置下列参数:→Input File(输入需要裁剪的图像文件名称):gy_resample.img→Input Mask File(输入掩模文件名称):raster.img→点击Setup Recode按钮设置裁剪区域内新值(New Value)为1,区域外取0值 →选择Intersection(确定掩模区域作交集运算)→Output File(输出图像文件名称):mask.img→Output (输出数据类型):Unsigned 8 bit→单击OK(关闭Mask对话框,执行掩模运算)2.图像几何校正(Geometric Correction)几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程;而将地图坐标系统赋予图像数据的过程,称为地理参考(Georeferencing)。