Erdas基础教程数据预处理,校正、配准、镶嵌、裁剪

Erdas基础教程: 数据预处理

来源：陕西师范大学旅游与环境学院

在ERDAS中，数据预处理模块为Data preparation。图标面板工具条中，点击图标——Data Preparation菜单

1、图象几何校正

第一步：显示图象文件

在视窗中打开需要校正的Landsat TM图象：lanzhoucity.img.

第二步：启动几何校正模块

在Viewer#1的菜单条中，选择Raster|Geometric Correction

.打开Set Geometric Model对话框

.选择多项式几何校正模型Polynomial——OK

.程序自动打开Geo Correction Tools对话框和Polynomial Model Properties对话框

.先选择Close关闭Polynomial Model Properties对话框.程序自动打开GCP Tool Reference Setup对话框

.选择Keyboard Only

.OK

.程序自动打开Reference Map Information提示框。

.选择Map Units: Meters

.添加地图投影参数，如下图：

.选择OK 确定地图投影参数，并关闭上图。

.选择OK，确定Reference Map Information，并关闭提示框。

.并自动打开采集控制点对话框。

GCP的具体采集过程:

在图像几何校正过程中，采集控制点是一项非常重要和相当重要的工作，具体过程如下：.在GCP工具对话框中点select GCP 图标，进行GCP选择状态。

.在view#1中移动关联方框位置，寻找明显地物特征点，作为输入GCP。

.在GCP工具对话框中点击Great GCP图标，并在view#2中点击左键定点，GCP数据表将记录一个输入的GCP，包括编号、标识码、X、Y坐标。

.在GCP工具对话框中输入地图参数坐标X、Y。

.不断重复上述步骤，采集若干GCP，直到满足所选是的几何校正模型为止。

采集地面检查点

以上所采集的GCP为控制点，用于建立转换方模型及多项式方程，地面检查点，则用于检验所建立的转换方程的精度和实用性，具体过程如下：

.在GCP TOOL菜单条中选择GCP类型：Edit/Set Point Type —check。

.在GCP TOOL 菜单条中确定GCP匹配参数：Edit/Point matching —打开GCP Matching 对话框，并确定参数。

.确定地面检查点，其操作与选择控制点完全一样。

.计算检查点误差：在GCP TOOL工具条中，点击Compute Error 图标，检查点的误差就会显示在GCP TOOL 的上方(如下图)，只有所有检查的误差小于一个像元时，才能进行以下的步骤。

.在Geo Correction Tools 对话框中选择Model Properties 图标打开

选择或检查参数，然后选择close关闭。

图像重采样

.在Geo Correction Tools 对话框中选择Image Resample 图标—打开Image Resample 对话框，并定义重采样参数。

.输出图像文件名(output file):rectify.img

.选择重采样方法(Resample Method):Nearest Neighbor

.定义输出图像在图与像元大小。

.设查输出统计中忽略零值。

.选择OK启动重采样进程，并关闭Image Resample 对话框。

2.图象拼接处理

本练习将同一区域机邻的三幅TM图象进行拼接。

其过程如下：

.启动图象拼接工具,在ERDAS图标面板工具条中，点击Dataprep/Data preparation/Mosaicc lmages—打开Mosaic Tool 视窗。

.加载Mosaic图像，在Mosaic Tool视窗菜单条中，Edit/Add images—打开Add Images for Mosaic 对话框。依次加载窗拼接的图像。

.在Mosaic Tool 视窗工具条中，点击set Input Mode 图标，进入设置图象模式的状态，利用所提供的编辑工具，进行图象叠置组合调查。

.图象匹配设置，点击Edit /Image Matching —打击Matching options 对话框，设置匹配方法：Overlap Areas。

.在Mosaic Tool视窗菜单条中，点击Edit/set Overlap Function—打开set Overlap Function对话框

设置以下参数：

.设置相交关系(Intersection Method):No Cutline Exists。

.设置重叠图像元灰度计算(select Function):Average。

.Apply —close完成。

.运行Mosaic 工具在Mosaic Tool视窗菜单条中,点击Process/Run Mosaic —打开Run Mosaic 对话框。

设置下列参数：确定输出文件名：mosaic.img

确定输出图像区域：ALL

OK进行图像拼接。

3、图象分幅裁剪

在实际工作中，经常根据研究区的工作范围进行图像分幅裁剪，利用ERDAS 可实现两种图像分幅裁剪：规则分幅裁剪，不规则分幅裁剪。

1、规则分幅裁剪

即裁剪的边界范围为一矩形，其具体方法如下：

在ERDAS图标面板工具条中，点击DataPrep/Data preparation/subset Image—打开subset Image 对话框，并设置参数如下：

说明：裁剪范围输入：

①通过直接输入右上角、右下角的坐标值;

②先在图像视窗中放置查询框，然后在对话框中选择From Inquire Box

③先在图像视窗中绘制AOL区域，然后在对话框中选择AIO功能，利用此方法也可实现不规则裁剪。

2、不规则分幅裁剪

①用AOI区域裁剪，与上述的②的方法相同。

②用Arclnfo的多边形裁剪。

.将Arclnfo多边形转换成网格图象。

.ERDAS图标面板工具条中，点击Vector/vector to Raster—打开vector to Raster对话框，并设置参数，并实现转换。

.通过掩膜运算实现图像不规则裁剪。

.ERDAS 图标面板工具条中，点击Interpreter/Utilities /Mask--打开Mask对话框，并设置参数如下：

setup Recode设置裁剪区域内新值为1，区域外取0值。