Erdas实验合集.doc

实验一 设置一张三维图像在erdas图标面板工具条中选择“VirtualGIS”“ VirtualGIS Viewer.”，就打开了“VirtualGIS Viewer”窗口:在窗口处选择添加图像，可以“File” “open” “DEM”，也可以用打开按钮打开添加窗口：添加需要的图像，如：DEMmerge-tub.img,然后在此界面选择“Raster Opinions”,选择“DEM”，OK。同样地，添加第二张图像，如XS_truecolor_sub.img，在“Raster Opinions”中选择“Raster Overlay”，Display as选择“true Color”或者根据需要选择，红绿蓝也可以改为3、2、1等，然后OK。然后在“VirtualGIS Viewer”中显示如下图：可以在图中国的红框位置选择不同显示方式，以观看最后的图像：最后效果如下（选择Target并拖动后）：下图是选择“Birds Eye”并调整后的效果：实验二 几何纠正第一步：显示图像文件：打开两个窗口，并列平铺，在其中某一个窗口，点击“view” “Tile Viewers”,两个图像就并列平铺在屏幕上了。然后在Viewer#1窗口中打开需要校正的Landsat TM图像tmAtlanta.img，在Viewer#2窗口中打开作为地理参考的校正过的SPOT图像panAtlanta.img第二步：启动几何校正模块：Viewer#1菜单条：RasterGeometric Correction，打开Set Geometric Model对话框：选择多项式几何校正计算模型：Polynomial OK。系统同时打开Geo Correction Tools 对话框和Polynomial Model Properties窗口，在Polynomial Model Properties窗口中，定义多项式模型参数及投影参数：选择多项式次方（Polynomial Order）：2定义投影参数：ProjectionMap Units：Meters单击Apply，再点击Close，打开GCP Tool Reference Setup对话框：第三步：启动控制点工具在GCP Tools Reference Setup对话框中选择采点模式：选择Existing Viewer，OK。系统打开Viewer Selection Instruction指示器，在显示作为地理参考图象panAtlanta.img的Viewer#2中单击，系统打开Reference Map Information提示框，OK。屏幕自动变化为包含两个主窗口，两个放大窗口，两个关联方框，控制点工具对话框和几何校正工具等。控制点工具被启动，进入控制点采集状态。第四步：采集地面控制点在GCP工具对话框中单击Select GCP图标，进入GCP选择状态：可以GCP数据表中设置输入GCP的颜色（color）：在Viewer#1中移动关联框位置，寻找明显地物特征点，作为GCP；在GCP工具对话框中单击Create GCP图标,并在Viewer#3中单击定点,GCP数据表记录一个GCP,包括编号,标识码,X坐标,Y坐标；在GCP工具对话框中单击Select GCP图标，重新进入GCP选择状态：在Viewer#2中移动关联框位置，寻找与Viewer3中GCP相同的地物特征点，作为GCP；在GCP工具对话框中单击Create GCP图标,并在Viewer#4中单击定点,系统将把参考点的坐标(X Reference,Y Reference)显示在GCP数据表中；在GCP工具对话框中单击Select GCP图标，重新进入GCP选择状态,并将光标移回到Viewer#1,准备采集另一个输入控制点。不断重复以上步骤,采集若干GCP,直到满足所选定的几何校正模型为止。第五步：图像重采样在Geo Correction Tools对话框中单击Image Resample图标,打开Resample (图象重采样)对话框,定义重采样参数:输出图象文件名(Output File)：定义输出象元大小(Output Cell Sizes)：X值30 ，Y值30设置输出统计中忽略零值：选中Ignore Zero in StatsOK（关闭Resample对话框,启动重采样过程）最后结果：实验三 图像镶嵌首先在ERDAS面板工具条中，点击“Dataprep”“Mosaic images” “Mosaic Tool”，打开Mosaic Tool视窗：然后加载图像：“Edit”“Add Images”, 或点击，把需要镶嵌的图像加载到视窗中：选择wasia1_mss.img，然后在“Image Opinions”中选择“Compute Active Area”，OK，把这幅图像添加到了视窗中，用同样的方法，把wasia2_mss.img、wasia3_tm.img都添加到视窗中：利用Mosaic Tool 视窗中的几个工具可以对图像进行调整叠放次序接下来进行图像色彩匹配设置：在Mosaic Tool 视窗工具条中点击 图标，进入设置图像关系模式点击图标，打开Set Overlap Function 对话框：在Set Overlap Function 对话框中，可以对Intersection Type（相交类型）和Select Function（重合区像元灰度函数）进行设置，选择后点击Apply，再点击Close 关闭Set Overlap Function 对话框。另外还可以在Mosaic Tool 视窗工具条中点击图标，打开Color Corrections 对话框：在Color Corrections 对话框中可以选择是否使用直方图匹配（最底下一个），如果选择使用，则选中后点击Set，打开Histogram Matching对话框：在这里可以选择匹配方法：Overlap Areas或者For All Images，选择后点击OK，再在Color Corrections对话框中点击OK。综合利用以上方法就可以对图像进行不同的匹配设置。最后定义输出图像在Mosaic Tool 视窗工具条中点击图标，进入设置输出图像模式；点击图标，打开Output Image Options 对话框：在Define Output Map Areas（定义输出地图区域）处设置：Union of All Inputs，然后OK。以上所有步骤设置完成后就可以输出结果了：在Mosaic Tool 视窗菜单条中点击Process；点击Run Mosaic，打开Run Mosaic 对话框；在Run Mosaic 对话框中，设置以下参数：设置输出文件名；输出图像区域：All；点击OK，打开Mosaic Modeler 进程状态条；点击OK，关闭状态条，完成数据输出。以下几张图为不同设置参数下的图样：（1） 不使用直方图匹配，选择Overlay：（2） 不使用直方图匹配，选择Average：（3） 不使用直方图匹配，选择Minimum：（4） 不使用直方图匹配，选择Maximum：（5） 不使用直方图匹配，选择Feather：（6） 使用直方图匹配，匹配方式选择For All Images，选择Feather：（7） 使用直方图匹配，匹配方式选择Overlap Areas，选择Feather：实验四 图像增强处理傅里叶高通/低通滤波图或效果图 空间增强效果图进行傅里叶高/低通滤波处理之前，首先要进行傅里叶快速变换：点击图标 “Fourier Analysis” “Fourier Transform”,打开“Fourier Transform”对话框，在对话框中选择输入文件，指定输出文件：点击OK进行快速傅里叶变换。然后打开傅里叶编辑器进行傅里叶变换；点击 图标 “Fourier Analysis” “Fourier Transform Editor”，打开“Fourier Editor（傅里叶编辑器）”，在傅里叶编辑器中输入上次装换后的fft文件：然后选择“Mask” “Filters”,打开Low/High Pass Filter对话框，Filter Type选择低通滤波，Window Function选择Ideal，然后OK：效果如下图所示：同样的，在Low/High Pass Filter对话框设置高通滤波后，效果如下图：实验五 光谱增强1、 主成分变换点击 图标 “Spectral Enhancement” “Principal comp”,打开Principal Components对话框，添加输入图像，指定输出位置，完成相关设置OK，原图和处理后的效果如下图：原图处理效果图特征矩阵ortho.mtx:0.630084002043407 -0.7714103194762875 0.088995895267462720.6386525809583911 0.449602837298354 -0.62448392255400120.4417205351850302 0.4503147873418882 0.77595074656649582、主成分逆变换：点击 图标 “Spectral Enhancement” “Inverse Principal comp”,打开Inverse Principal Components对话框，添加输入图像，输入特征矩阵，指定输出位置，完成相关设置，OK，原图和处理后的效果如下图：原图处理后效果图2、自然彩色变换：点击 图标 “Spectral Enhancement” “Natural Color”,打开Natural Color对话框，添加输入图像，指定输出位置，完成相关设置OK，原图和处理后的效果如下图：原图处理效果图3、 缨帽变换点击 图标 “Spectral Enhancement” “Tasseled Cap”,打开Tasseled Cap对话框，添加输入文件，指定输出位置，设置相关参数后，OK，原图和处理后的效果如下图：原图处理效果图4、RGB to HIS变换点击 图标 “Spectral Enhancement” “RGB to HIS”,打开RGB to HIS对话框，添加输入文件，指定输出位置，设置相关参数后，OK，原图和处理后的效果如下图：原图处理后效果图4、 HIS to RGB变换：点击 图标 “Spectral Enhancement” “HIS to RGB”,打开HIS to RGB对话框，添加输入文件，指定输出位置，设置相关参数后，OK，原图和处理后的效果如下图：原图处理后效果图5、 指数计算点击 图标 “Spectral Enhancement” “Indices”,打开Indices对话框，添加输入文件，指定输出位置，设置相关参数后，OK，原图和处理后的效果如下图：原图处理后效果图实验六 非监督分类ERDAS图标面板工具条：点击Classifier图标Unsupervised classificationUnsupervised classification对话框。在Unsupervised classification对话框中，确定输入图像文件，定义输出分类图像文件，选择生成分类模板文件(output Signature file)：feijiandu.sig；选择分类参数决定方法： Use Signature Means；确定初始分类数 (Number of classes)：10，定义最大循环次数 (Maximum Iteration)：30；设置循环收敛阈值 (Convergence Threshold)：0.95；OK，原图和处理后的效果如下图：原图处理后效果图分类后处理：1、 聚类统计点击Interpreter图标GIS AnalysisClumpClump对话框，添加输入图像，指定输出位置，设置相关参数，OK：原图和处理后的效果如下图：原图处理后效果图2、 过滤分析点击Interpreter图标GIS AnalysisSieveSieve对话框，添加输入图像，指定输出位置，设置相关参数，OK：原图和处理后的效果如下图：原图处理后效果图3、 去除分析点击Interpreter图标GIS AnalaysisEliminateEliminate对话框，添加输入图像，指定输出位置，设置相关参数，OK：原图和处理后的效果如下图：原图处理后效果图4、 分类重编码点击Interpreter图标GIS AnalysisRecodeRecode对话框，添加输入图像，指定输出位置，设置相关参数，点击Setup Recode按钮打开Thematic Recode表格：根据需要改变New Value字段的取值(直接输入)在本例中将原来的十类依次两两合并，形成五类。原图和处理后的效果如下图：原图处理后效果图实验七 监督分类1、 打开需要分类的图像2、 打开分类模板编辑器：点击ClassifierSignature Editor：3、打开AIO工具，选择多边形按钮，在图像窗口中选择与类别对应的区域，绘制多边形AIO。4. 建立分类模板，在分类模板编辑器中，点击按钮，将多边形区域加载到分类模板中属性表中。再多绘制一些AIO，分别加载到分类模板属性表中。在分类模板中，可以依次单击AIO，并单击将所选模板合成一个新的模板，定义分类类别名称并设置颜色。5、 重复以上步骤，将所有类型建立分类模板。分类所有水体、纸杯、居民区、道路等：6、 保存分类模板：执行监督分类点击Classifier图标classificationSupervised classification，进入对话框，确定输入原始文件，指定输出位置，确定分类模板文件,完成相关设置：原图和处理后的效果如下图：原图监督分类后图实验八 GIS分析 地形分析GIS分析1、 邻域分析点击Interpreter图标GIS AnalysisNeighborhood Functions，打开对话框，添加输入图像，指定输出位置，设置相关参数，OK：原图和处理后的效果如下图：原图邻域分析后图2、 查找分析：点击Interpreter图标GIS AnalysisSearch，打开对话框，添加输入图像，指定输出位置，设置相关参数，OK：原图依然是“Inlandc.img”,处理后的效果如下图：3、 指标分析：点击Interpreter图标GIS AnalysisIndex，打开对话框，添加输入图像，指定输出位置，设置相关参数，OK：输入的两幅图：处理后的图：4、 叠加分析点击Interpreter图标GIS AnalysisOverlay，打开对话框，添加输入图像，指定输出位置，设置相关参数，OK：输入的两个图像：处理后的图像：地形分析1、 坡度分析点击Interpreter图标Topographic AnalysisSlope，打开对话框，添加输入图像，指定输出位置，设置相关参数，OK：原图和处理后的图效果如下：原图坡度分析后图2、 坡向分析点击Interpreter图标Topographic AnalysisAspec