遥感数据图像处理实验三遥感图像的几何校正与裁剪.docx
《遥感数据图像处理实验三遥感图像的几何校正与裁剪.docx》由会员分享,可在线阅读,更多相关《遥感数据图像处理实验三遥感图像的几何校正与裁剪.docx(49页珍藏版)》请在冰豆网上搜索。
遥感数据图像处理实验三遥感图像的几何校正与裁剪
实验三、遥感图像的几何校正与裁剪
实验内容:
1.图像分幅裁剪(SubsetImage
2.图像几何校正(GeometricCorrection
3.图像拼接处理(MosaicImgaes
4.生成三维地形表面(3DSurfacing
1.图像分幅裁剪
在实际工作中,经常需要根据研究工作范围对图像进行分幅裁剪,按照ERDASIMAGINE8.4实现图像分幅裁剪的过程,可以将图像分幅裁剪为两类型:
规则分幅裁剪,不规则分幅裁剪。
1.1规则分幅裁剪
(以c:
\ProgramFile\IMAGINE8.4\examples\lanier.img为例
规则分幅裁剪是指裁剪图像的范围是一个矩形,通过左上角和右上角两点的坐标可以确定图像的裁剪位置,过程如下:
方法一:
→ERDASIMAGINE8.4图标面板菜单条:
Main→DataPreparation(或单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标
→打开DataPreparation对话框
→单击SubsetImage按钮,打开Subset对话框
在Subset对话框中需要设置下列参数:
→输入文件名(InputFile:
lanier.img
→输出文件名(OutputFile:
lanier_sub.img
→坐标类型(CoordinateType:
Map
→裁剪范围(SubsetDefinition:
ULX、ULY、LRX、LRY
(注:
ULX,ULY是指左上角的坐标,LRX,LRY是指右上角的坐标,缺省状态为整个图像范围
→输出数据类型(OutputDataType:
Unsigned8Bit
→输出文件类型(OutputLayerType:
Continuous
→输出统计忽略零值:
IgnoreZeroInOutputStats
→输出像元波段(SelectLayers:
2,3,4
→OK(关闭Subset对话框,执行图像裁剪
方法二:
→ERDASIMAGINE8.4图标面板菜单条:
Main→StartIMAGINEViewer(或单击RDASIMAGINE8.4图标面板工具条“Viewer”图标
→打开一个二维视窗
→单击视窗工具条最左端的“打开文件”图标
→打开SelectLayerToAdd对话框
在SelectLayerToAdd对话框完成以下设置:
→LookIn:
examples
→FileName:
lanier.img
→Filesoftype:
IMAGINEImage
→双击OK按钮
→在二维视窗中打开lanier.img文件
→单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标→打开DataPreparation对话框
→单击SubsetImage按钮
→打开Subset对话框
在Subset对话框中需要设置下列参数:
→输入文件名(InputFile:
lanier.img
→输出文件名(OutputFile:
lanier_sub.img
→坐标类型(CoordinateType:
Map
→输出数据类型(OutputDataType:
Unsigned8Bit→输出文件类型(OutputLayerType:
Continuous
→输出统计忽略零值:
IgnoreZeroInOutputStats
→输出像元波段(SelectLayers:
2,3,4
→单击FromInquireBox按钮
→打开InvalidCoordinateType对话框
→单击Continue
→在显示图像文件lanier.img视窗中单击工具条的“+”按钮,打开InquireCursor对话框,在视窗中移动十字光标,确定裁剪范围左上角和右下角,读取其坐标分别填入SubsetImage对话框的ULX,ULY中和LRX,LRY中
→单击OK按钮(关闭Subset对话框,执行图像裁剪
方法三:
首先在视窗中打开lanier.img文件
→AOI→Tools打开AOI工具面板
→单击矩形框确定裁剪范围
→File→Save→AOILayerAs
→打开SaveAOIAs对话框,输入文件名:
2
→单击OK(退出SaveAOIAs对话框
→单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标→打开DataPreparation对话框
→单击SubsetImage按钮
→打开Subset对话框
在Subset对话框中需要设置下列参数:
→输入文件名(InputFile:
lanier.img
→输出文件名(OutputFile:
lanier_sub.img
→坐标类型(CoordinateType:
Map
→输出数据类型(OutputDataType:
Unsigned8Bit
→输出文件类型(OutputLayerType:
Continuous→输出统计忽略零值:
IgnoreZeroInOutputStats→输出像元波段(SelectLayers:
2,3,4
→单击AOI按钮
→打开ChooseAOI对话框
→在ChooseAOI对话框作如下设置:
→AOISource:
File
→AOIFile:
2
→单击OK(退出ChooseAOI对话框
→单击OK(退出Subset对话框,执行图像裁剪
→单击OK(退出Modeler对话框,完成图像裁剪
1.2不规则分幅裁剪
不规则分幅裁剪是指裁剪图像的边界范围是个任意多边形,无法通过左上角和右下角两点的坐标确定图像的裁剪位置,而必须事先生成一个完整的闭合多边形区域,可以是一个AOI多边形,也可以是ArcInfo的一个PolygonCoverage,针对不同的情况采用不同的裁剪过程。
(一AOI多边形裁剪
(以c:
\ProgramFile\IMAGINE8.4\examples\lanier.img为例
→在视窗中打开需要裁剪的图像
→应用AOI工具绘制多边形
可以将多边形AOI文件保存在文件中(*.aoi,也可以暂时不退出视窗,将图像与AOI多边形保留在视窗中,然后进行以下操作:
→ERDASIMAGINE8.4图标面板菜单条:
Main→DataPreparation(或单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标,打开DataPreparation对话框
→单击SubsetImage按钮,打开Subset对话框
在Subset对话框中需要设置下列参数:
→输入文件名(InputFile:
lanier.img
→输出文件名(OutputFile:
lanier_sub.img→CoordinateType:
Map
→Output:
Unsigned8bit
→Oupput:
Continuous
→选定IgnoreZeroinOutputStats
→SelectLayers:
2,3,4
→单击AOI按钮,打开ChooseAOI对话框
→选定File
→选择文件:
2.aoi
→单击OK(退出ChooseAOI对话框
→单击OK(退出Subset对话框,执行图像裁剪
→单击Modeler对话框中的OK按钮完成图像裁剪
(二ArcInfo多边形裁剪
(以c:
\ProgramFile\IMAGINE8.4\examples\gy_resample.img为例
如果是按照行政区划边界或自然区划边界进行图像的分幅裁剪,往往是首先利用ArcInfo或者ERDAS的Vector模块绘制精确的边界多边形(Polygon,然后以ArcInfo的Polygon为边界进行图像裁剪。
对于这种情况,需要调用ERDAS其它模块的功能分两步完成。
第一步:
将ArcInfo多边形转换为栅格图像文件
打开VectorToRaster对话框有两种方法:
方法一:
ERDASIMAGINE8.4图标面板菜单条:
Main→ImageInterpreter(或单击ERDASIMAGINE8.4图标面板工具条“Interpreter”图标
→打开ImageInterpreter对话框
→单击Utilities
→打开Utilities对话框
→选择VectorToRaster
→打开VectorToRaster对话框
方法二:
ERDASIMAGINE8.4图标面板菜单条:
Main→Vector(或单击ERDASIMAGINE8.4图标面板工具条“Vector”图标
→打开VectorUtilities对话框
→选择VectortoRaster
→打开VectortoRaster对话框
本例中选择前一种方法:
在VectortoRaster对话框中设置下列参数:
→InputVectorFile(输入矢量文件名称:
boundgy.shp→VectorType(确定矢量文件类型:
polygon
→UseAttributeasValue(使用矢量属性值:
GYBOUND_ID
→OutputFile(输出栅格文件名称:
raster.img
→DataType(栅格数据类型:
Unsigned8bit
→LayerType(栅格文件类型:
Thematic
→SizeDifinition(转换范围大小:
ULX,ULY,LRX,LRY(缺省条件下为整个图像范围→Units(坐标单位:
Meters
→CellSize(输出像元大小:
X:
30/Y:
30
→选择SquareCell(正方形像元
→单击OK(关闭VectorToRaster对话框,执行矢量到栅格的转换
第二步:
通过掩模运算(Mask实现图像的不规则裁剪
ERDASIMAGINE8.4图标面板菜单条:
Main→ImageInterpreter(或单击ERDASIMAGINE8.4图标面板工具条“Interpreter”图标
→ImageInterpreter
→Utilities
→Mask
→打开Mask对话框
在Mask对话框中设置下列参数:
→InputFile(输入需要裁剪的图像文件名称:
gy_resample.img
→InputMaskFile(输入掩模文件名称:
raster.img
→点击SetupRecode按钮设置裁剪区域内新值(NewValue为1,区域外取0值→选择Intersection(确定掩模区域作交集运算
→OutputFile(输出图像文件名称:
mask.img
→Output(输出数据类型:
Unsigned8bit
→单击OK(关闭Mask对话框,执行掩模运算
2.图像几何校正(GeometricCorrection
几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程;而将地图坐标系统赋予图像数据的过程,称为地理参考(Georeferencing。
由于所有地图投影系统都遵从于一定的地图坐标系统,所以几何校正过程包含了地理参考过程。
在ERDASIMAGINE8.4中进行图像几何校正,通常有两种途径启动几何校正模块:
(1数据预处理途径
(以c:
\ProgramFi