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

1、遥感数据图像处理实验三遥感图像的几何校正与裁剪实验三、遥感图像的几何校正与裁剪实验内容:1.图像分幅裁剪(Subset Image2.图像几何校正(Geometric Correction3.图像拼接处理(Mosaic Imgaes4.生成三维地形表面(3D Surfacing1.图像分幅裁剪在实际工作中,经常需要根据研究工作范围对图像进行分幅裁剪,按照ERDAS IMAGINE 8.4实现图像分幅裁剪的过程,可以将图像分幅裁剪为两类型:规则分幅裁剪,不规则分幅裁剪。1.1规则分幅裁剪(以c:Program File IMAGINE 8.4exampleslanier.img为例规则分幅裁剪是

2、指裁剪图像的范围是一个矩形,通过左上角和右上角两点的坐标可以确定图像的裁剪位置,过程如下:方法一:ERDAS IMAGINE 8.4 图标面板菜单条:MainData Preparation(或单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标打开Data Preparation 对话框 单击Subset Image按钮,打开Subset对话框 在Subset对话框中需要设置下列参数:输入文件名(Input File:lanier.img输出文件名(Output File:lanier_sub.img坐标类型(Coordinate Type:Map裁剪范围(Subs

3、et 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图标面板菜单条:MainStart IMAGINE Viewer(或单击RDAS IMAGIN

4、E 8.4图标面板工具条“Viewer”图标打开一个二维视窗单击视窗工具条最左端的“打开文件”图标打开Select Layer To Add对话框在Select Layer To Add对话框完成以下设置:Look In:examplesFile Name:lanier.imgFiles of type:IMAGINE Image双击OK按钮在二维视窗中打开lanier.img文件 单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标 打开Data Preparation对话框 单击Subset Image按钮打开Subset对话框 在Subset对话框中需要设置下列

5、参数:输入文件名(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视

6、窗中单击工具条的“+”按钮,打开Inquire Cursor 对话框,在视窗中移动十字光标,确定裁剪范围左上角和右下角,读取其坐标分别填入Subset Image对话框的ULX,ULY中和LRX,LRY中单击OK按钮(关闭Subset对话框,执行图像裁剪方法三:首先在视窗中打开lanier.img文件 AOITools打开AOI工具面板 单击矩形框确定裁剪范围 FileSaveAOI Layer As打开Save AOI As对话框,输入文件名:2 单击OK(退出Save AOI As对话框单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标 打开Data Prepa

7、ration对话框 单击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对话

8、框 在Choose AOI对话框作如下设置:AOI Source:FileAOI File:2 单击OK(退出Choose AOI对话框单击OK(退出Subset对话框,执行图像裁剪单击OK(退出Modeler对话框,完成图像裁剪1.2不规则分幅裁剪不规则分幅裁剪是指裁剪图像的边界范围是个任意多边形,无法通过左上角和右下角两点的坐标确定图像的裁剪位置,而必须事先生成一个完整的闭合多边形区域,可以是一个AOI多边形,也可以是ArcInfo的一个Polygon Coverage,针对不同的情况采用不同的裁剪过程。(一 AOI多边形裁剪(以c:Program File IMAGINE 8.4exam

9、pleslanier.img为例在视窗中打开需要裁剪的图像 应用AOI工具绘制多边形 可以将多边形AOI文件保存在文件中(*.aoi,也可以暂时不退出视窗,将图像与AOI 多边形保留在视窗中,然后进行以下操作:ERDAS IMAGINE 8.4 图标面板菜单条:MainData Preparation(或单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标,打开Data Preparation 对话框 单击Subset Image按钮,打开Subset对话框 在Subset对话框中需要设置下列参数:输入文件名(Input File:lanier.img输出文件名(Ou

10、tput File:lanier_sub.img Coordinate Type:MapOutput:Unsigned 8 bitOupput:Continuous选定Ignore Zero in Output StatsSelect Layers:2,3,4 单击AOI按钮,打开Choose AOI对话框 选定File选择文件:2.aoi 单击OK(退出Choose AOI对话框单击OK(退出Subset对话框,执行图像裁剪单击Modeler对话框中的OK按钮完成图像裁剪(二 ArcInfo多边形裁剪(以c:Program File IMAGINE 8.4examplesgy_resampl

11、e.img为例如果是按照行政区划边界或自然区划边界进行图像的分幅裁剪,往往是首先利用ArcInfo或者ERDAS的Vector模块绘制精确的边界多边形(Polygon,然后以ArcInfo的Polygon为边界进行图像裁剪。对于这种情况,需要调用ERDAS其它模块的功能分两步完成。第一步:将ArcInfo多边形转换为栅格图像文件打开Vector To Raster对话框有两种方法:方法一:ERDAS IMAGINE 8.4 图标面板菜单条:MainImage Interpreter (或单击ERDAS IMAGINE 8.4 图标面板工具条“Interpreter”图标打开Image Inte

12、rpreter对话框 单击Utilities打开Utilities对话框 选择Vector To Raster打开Vector To Raster对话框 方法二:ERDAS IMAGINE 8.4 图标面板菜单条:MainVector (或单击ERDAS IMAGINE 8.4 图标面板工具条“Vector”图标打开Vector Utilities对话框 选择Vector to Raster打开Vector to Raster对话框 本例中选择前一种方法:在Vector to Raster对话框中设置下列参数:Input Vector File(输入矢量文件名称:boundgy.shp Vec

13、tor Type(确定矢量文件类型:polygonUse Attribute as Value(使用矢量属性值:GYBOUND_IDOutput File(输出栅格文件名称:raster.imgData Type(栅格数据类型:Unsigned 8 bitLayer Type(栅格文件类型:ThematicSize Difinition(转换范围大小:ULX,ULY,LRX,LRY(缺省条件下为整个图像范围 Units(坐标单位:MetersCell Size(输出像元大小:X:30/Y:30选择Square Cell(正方形像元 单击OK(关闭Vector To Raster对话框,执行矢量

14、到栅格的转换 第二步:通过掩模运算(Mask实现图像的不规则裁剪ERDAS IMAGINE 8.4 图标面板菜单条:MainImage Interpreter (或单击ERDAS IMAGINE 8.4 图标面板工具条“Interpreter”图标Image InterpreterUtilitiesMask打开Mask对话框 在Mask对话框中设置下列参数:Input File(输入需要裁剪的图像文件名称:gy_resample.imgInput Mask File(输入掩模文件名称:raster.img点击Setup Recode按钮设置裁剪区域内新值(New Value为1,区域外取0值

15、选择Intersection(确定掩模区域作交集运算Output File(输出图像文件名称:mask.imgOutput (输出数据类型:Unsigned 8 bit 单击OK(关闭Mask对话框,执行掩模运算 2.图像几何校正(Geometric Correction几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程;而将地图坐标系统赋予图像数据的过程,称为地理参考(Georeferencing。由于所有地图投影系统都遵从于一定的地图坐标系统,所以几何校正过程包含了地理参考过程。在ERDAS IMAGINE 8.4中进行图像几何校正,通常有两种途径启动几何校正模块:(1 数据预处理途径(以c:Program Fi