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值。