erdas使用详细操作教程.doc

erdas专题制图详细操作 实验数据： 1、裁剪144029(行列号)所需矢量文件 2、卫星影像：LE71440292000188SGS01（数据获取来源为usgs网站：/下载所得。） 初步查看分析影像信息： 首先在erdas中用tiff格式打开我们下载的影像: 点击查看我们下载的影像的信息： 可以看到我们下载的影像的分辨率，投影信息等。 图片格式转换： 过程如下： 在import\expot窗口中： 这里的Type我们可以选择tiff 也可以选择geotiff，本人尝试过，选择这两种任意一种类型对输出影像几乎没有区别。（其区别还有待查证） 以此方法，将七个波段的.tiff影像图都转换成.img格式的影像。由于band8 为全色影像，我们的实验中不需要用它，所以这里我不进行转换了。只转换以下七个波段即可。 波段组合：将七个波段组合成.img文件 设置完毕点OK。于是我们就得到了 .img格式的影像文件。 几何精校正: 前面已经知道这张影像的地理信息了，为了验证其实图片投影大致的准确性，我们可以进行几何精校正进行验证。几何精校正校正过程如下： 这里我们使用多项式几何校正。 其中，多项式变换（Polynomial）在卫星图像校正过程中应用较多，在调用多项式模型时，需要确定多项式的次方数（Polynomial Order），通常整景图像选择3次方。次方数与所需要的最少控制点数是相关的，最少控制点数计算公式为(（t+1）*(t+2))/2,式中t为次方数，即1次方最少需要3个控制点，2次方最少需要6个控制点，3次方至少需要10个控制点，依次类推。 设置完毕点close。 会弹出GCP Tool Refrence Setup窗口。 由于我们要使用的是谷歌上的数据，因此我们要选择用键盘输入。 确认信息: 信息正确，点OK。 在谷歌中设置：菜单栏中：工具-选项 3D视图选项卡中，显示纬度\经度项中 使用通用横轴墨卡托投影。 设置完后点应用和确定即可。 在谷歌中以及卫星影像中分别选取同名点， 找到后，点击中的（或者点击窗口中的均可。） 在窗口点击左键。如上图所示。 再将谷歌上该点的坐标复制到erdas的以下窗口中 的对应点的列的下方。 重复此步骤，由于校正时，最佳效果是其最少需要的控制点的2至3倍，而我们使用的三次多项式，则我们至少要选择10，为了效果更佳，因此我一共选择了25个控制点。控制点整体分布如下图：（控制点要分布均匀，应该尽量选择误差较小的点，尽可能的减少误差，以其保证准确性。） 控制点结果如下表 进行控制点检查： 点计算检查点误差，由上图可看到，检查点的总误差为0.6658。因此，我们的校正结果是合格的。 输出配准后图像： 点窗口中的 设置完毕点OK。 输出图像如下图 我们将其与原图叠加观察其与原图像区别：添加原影像，在raster options 窗口中将clear display 前的对勾去掉即可。 两张图同时打开中该窗口中后，再进行如下操作。 弹出以下窗口： 仔细观察上图可以看到图像的上边界有微小变化，为了使效果更清晰，我将几何精校正后的图像进行增强，然后再与原图像进行叠加，来进行粗略的比较，比较结果如图：左图为几何精校正增强后，右图为几何精校正之前，从图中纹理又叠合情况来看，图中纹理变化不大。（如下图） 几何精校正结果分析： 综合整景影像分析校正前后影像变化，本人认为，下载所得原图像已经具有的投影信息已经够我们使用了，所以，做不做几何精校正这一步，对后面我们要进行的操作影像不会很大。 用shp矢量文件来裁剪影像 （2）在ERDAS的View窗口中，选择File—Open—Vector Layer，添加我们的“裁剪矢量”.shp文件，点击OK将其加载到当前窗口中。 （3）在窗口中，点击边界多边形，使其处于选中状态（即变为黄色，ERDAS中默认的未选中状态为淡蓝色），然后在AOI菜单中选则Copy Selection to AOI…，此时，被选中多边形的边界变为虚线。 （4）在窗口中，点击刚才转为AOI的多边形，其周围出现一个矩形框，将整个边界包围起来（即边界的外界矩形），在File菜单中选择Save—AOI Layer as…，在对话框中设置AOI文件的名称及保存路径，点击将AOI文件保存在磁盘上已被后面使用。 （5）关闭当前View窗口，在 ERDAS主菜单上选择DataPrep—Sunset Image…，在弹出的对话框中设定要裁剪的影像位置，裁剪后影像的输出位置，是否在输出统计中忽略零值。点击AOI…按钮，在弹出的对话框中选择AOI File，并指定刚才保存的AOI文件的路径，点击确定回到Subset对话框，点击确定按钮开始裁剪。裁剪操作如下： 裁剪过程如下： 添加完毕点击OK。回到subset窗口中点O