erdas影像校正剪切及不同时期影像地物动态变化方法..docVIP

遥感工程技术应用实习报告 作者：石宁卓 单位：西安科技大学 实习内容 （一）数据来源： 1. 应用1999年航空摄影测量方法制作的1：1万地形图的几何精度，对较新SPOT卫星影像进行几何校正。 2. 应用2002年获取的10分辨率SPOT-2影像对 1：10000地形图进行更新。 （二）现有数据： 1. 1999年西安地区1：10000栅格地形图（9幅） 2. 2002年西安地区10米分辨率SPOT-2全色影像。 （三）主要工作过程： 1. 对于扫描地形图的影像纠正、裁切、拼接 2. 利用拼接好的地形图对遥感影像进行几何纠正 3. 利用遥感影像对地形图进行更新 （四）采用软件： 在本次实习中采用的软件是遥感影像处理软件ERDAS IMAGINE9.2。 一 、对于扫描地形图的影像纠正、裁切、拼接 由于扫描地形图是利用已有的纸质地形图扫描而成的，在扫描过程中会存在着各种变形，所以在利用扫描地形图进行图像处理及数据采集之前需要对其进行纠正。 纠正原理：通过图廓坐标对扫描地形图进行纠正。采用ERDAS中图像预处理模块中的几何纠正功能，在扫描地形图上选择一定数量的控制点，然后通过图上坐标判读，在控制点的参考坐标中输入读取的坐标值，并进行重采样，从而对扫描地形图进行纠正。 （一）步骤： 1 .格式转换 IMPORT模块 将tif的地形图转换成为img格式。tif可以在ERDAS中直接打开，但转换格式之后可以使其操作方便。 2. 分别对九幅图做几何校正 以为九幅地形图是纸质打印扫描出来的，纸质扫描和图纸时间久了出现图像误差故要对九幅图纸一一做几何坐标匹配。步骤如下 设置：session——preferences----viewer-----clear display的钩去掉，使多幅图像可在一个窗口中打开。 点击Viewers模块——打开地形图； 点Raster----Geometric Correction-----在Set Geometric Model中，选Polynomial（多项式变换），点OK. 在Polynomial Model Properties中： Polynomial Order（多项式次数）设为2 点Projection选项：Map Units：Meters 点Add/Change Projection， 点Custom，进行设置 Projection Type（投影系统）：Gauss Kruger Spheroid Name（参考椭球）：IAG-75 Datum Name（基准面）: xian 80 Longitude of central meridian：108:00:00.000000E（1：10000图用的3度带，3*36=108，若为36度带） Latitude of origin of projection:0:00:00.000000N False easting：500000metres False northing:0.000000metres 设置完后,点保存，同时也可以将自己设置的坐标投影信息另存为一个文件，方便每次调用。如果另存为文件可以在Standard中选择刚才保存的文件，确定。 点Set Projection form GCP Tool，在GCPTool Teference Setup中选择Keyboard Only,确定。 在Polynomial Model Properties中点Apply,然后关闭 在图上选择7个控制点，并根据经纬线输入其校正后x、y值，观察其精度情况. 校正： Output Files: 输入文件名 Resample Method：Bilinear Interpolation(双线性插值法) Output Cell Sizes：X:0.85 Y:0.85 OK 下图为具体步骤： 步骤一 步骤二 步骤三 步骤四 3. 裁切、拼接九幅图 在主窗口中选择DataPrep——Subset Images，裁切地形图； 在主窗口中选择DataPrep——Mosaic Images, Edit——Add Images,装入九幅图 Process——Run Mosaic,输拼接后文件名，开始拼接（结果如下图） 4 . 用拼接后的地形图对SPOT影像进行几何校正 步骤一 步骤二 拼接后的地形图 二、spot 遥感影像几何校正 设置：session——preferences----viewer-----clear display的钩去掉 点Viewers,在Viewers1中打开遥感影像，再点击Viewers，在Viewer2中