2018年erdas使用详细操作.docVIP

erdas使用详细操作 　　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