遥感影像室内判读和监督分类.docVIP

实习内容一 遥感影像室内判读和监督分类 （可用Erdas或Envi） 主要数据： (1)分辨率30m的TM遥感影像（待校正）：taian100427tm.img, taian100427tm.tif (2)研究区地形图（已校正）：taianunion2009.img，taianunion2009.tif （3）分辨率2.5m的SPOT全色数据（待校正）：taian090921spot.img和taian090921spot.tif，可作为分类参考。 第一部分、遥感影像的几何精校正——Erdas操作 校正方法：用已校正的地形图（taianunion2009.img）去校正遥感影像（taian100427tm.img），步骤如下： 注：由于地形图的坐标为xian80，而erdas中没有IAG75椭球体和xian80坐标，先要进行以下操作： 在ERDAS安装目录下找到spheroid.tab文件(ERDAS8x版本一般在C:\Program Files \imagine\etc\; ERDAS9x版本一般在C:\Program Files\Leica Geosystems\Geospatial Imaging 9.2\etc\)，用记事本打开，在spheroid.tab文件末尾加入如下语句即可，假设spheroid.tab文件中最后一个椭球体序号为74（可以在文件最后一个椭球体中读出序号），则加入： IAG75 { 75 6378140 6356755.2882 xian 80 0 0 0 0 0 0 0 } 经过以上的操作，重启ERDAS，IAG75 椭球就会出现在ERDAS的椭球选择列表中。 1、打开两个viewer窗口，在一个窗口中打开待校正的遥感影像，另一个窗口中打开地形图； 2、在遥感影像的viewer窗口中选择raster→→geometric correction 1）选择多项式几何校正计算模型：polynomial（多项式变换（同时做投影变换））→→ok Polynomial model properties对话框中： polynomial order:3 注：最少控制点数计算公式为（（t+1）*(t+2)）/2, 即3次方至少需要10个控制点 2）Projection: mapunits选择meters; add/change projection→→custom Projection type 选择Gauss kruger Spheroid name（椭球体）选择IAG75，Datum name 选择xian80，longitude of central meridian 改为117：00：00.000000E，False easting改000000meters，其他默认，ok→→close 注：只要参考地形图的坐标设置好了，这一步可以省略，在第三步完成后，projection自动产生。 3）Gcp tool reference setup对话框：existing viewer→→ok, 出现viewer selection instructions后，用鼠标点一下打开地形图的viewer窗口范围内，ok 4）开始（检查点）控制点采集： 用Geo correction tool中的create GCP图标（圆形）进行采集，找准遥感影像与地形图的对应点后，用create GCP图标分别在两个viewer窗口中对应的位置采点，在GCP Tool表中会看到：采集遥感影像点时出现xinput 和yinput 的坐标值，采集地形图点时出现xref 和yref的坐标值。 5） 采集到十个点以后，会自动在坐标值后出现转换模型的误差计算值RMS Error，此时注意：Gcp tool对话框首先设置edit→point prediction(按照转换模型计算下一点位置)→取消predict next；edit→point matching(借助像元的灰度值匹配控制点)→取消match next。 采集到的gcp点，一定要保存：file→→save input as和save reference as,保存在两个文件中，用的时候可以用load input 和load reference。此处采集至少二十个控制点。控制点要选在明显的地物点，如道路与道路的交叉点、河流桥梁与道路的交叉点、道路的拐角、堤坝的拐角处等。控制点的空间分布力求均匀，而且要尽量覆盖全区。 所有点采集完后，点击unlock图标（表工具栏中的锁，变为锁住状） 6）计算检查控制点的误差， 确保所有点的误差小于一个像元，如果大于，则将其纪录选中后，edit →→set point type→→改为check，这个点就不会参与转