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

《erdas影像校正剪切及不同时期影像地物动态变化方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、遥感工程技术应用实习报告作者：石宁卓单位：西安科技大学实习内容（一）数据来源：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.利用遥感影像对地形图进行更新（四）采用软件

2、：在本次实习中采用的软件是遥感影像处理软件ERDASIMAGINE9.2。一、对于扫描地形图的影像纠正、裁切、拼接由于扫描地形图是利用已有的纸质地形图扫描而成的，在扫描过程中会存在着各种变形，所以在利用扫描地形图进行图像处理及数据采集之前需要对其进行纠正。纠正原理：通过图廓坐标对扫描地形图进行纠正。采用ERDAS中图像预处理模块中的几何纠正功能，在扫描地形图上选择一定数量的控制点，然后通过图上坐标判读，在控制点的参考坐标中输入读取的坐标值，并进行重采样，从而对扫描地形图进行纠正。（一）步骤：1.格式转换IMPORT模块将tif的地形

3、图转换成为img格式。tif可以在ERDAS中直接打开，但转换格式之后可以使其操作方便。2.分别对九幅图做几何校正以为九幅地形图是纸质打印扫描出来的，纸质扫描和图纸时间久了出现图像误差故要对九幅图纸一一做几何坐标匹配。步骤如下设置：session——preferences----viewer-----cleardisplay的钩去掉，使多幅图像可在一个窗口中打开。点击Viewers模块——打开地形图；点Raster----GeometricCorrection-----在SetGeometricModel中，选Polynomial（

4、多项式变换），点OK.在PolynomialModelProperties中：PolynomialOrder（多项式次数）设为2点Projection选项：MapUnits：Meters点Add/ChangeProjection，点Custom，进行设置ProjectionType（投影系统）：GaussKrugerSpheroidName（参考椭球）：IAG-75DatumName（基准面）:xian80Longitudeofcentralmeridian：108:00:00.000000E（1：10000图用的3度带，3*36=

5、108，若为36度带）Latitudeoforiginofprojection:0:00:00.000000NFalseeasting：500000metresFalsenorthing:0.000000metres设置完后,点保存，同时也可以将自己设置的坐标投影信息另存为一个文件，方便每次调用。如果另存为文件可以在Standard中选择刚才保存的文件，确定。点SetProjectionformGCPTool，在GCPToolTeferenceSetup中选择KeyboardOnly,确定。在PolynomialModelPrope

6、rties中点Apply,然后关闭在图上选择7个控制点，并根据经纬线输入其校正后x、y值，观察其精度情况.校正：OutputFiles:输入文件名ResampleMethod：BilinearInterpolation(双线性插值法)OutputCellSizes：X:0.85Y:0.85OK下图为具体步骤：步骤一步骤二步骤三步骤四3.裁切、拼接九幅图在主窗口中选择DataPrep——SubsetImages，裁切地形图；在主窗口中选择DataPrep——MosaicImages,Edit——AddImages,装入九幅图Proce

7、ss——RunMosaic,输拼接后文件名，开始拼接（结果如下图）4.用拼接后的地形图对SPOT影像进行几何校正步骤一步骤二拼接后的地形图二、spot遥感影像几何校正1、设置：session——preferences----viewer-----cleardisplay的钩去掉2、点Viewers,在Viewers1中打开遥感影像，再点击Viewers，在Viewer2中打开地形图；点session——tileviewers，遥感影像和地形图并排显示；1、在遥感影像的视窗中，点Raster----GeometricCorrectio

8、n-----在SetGeometricModel中，选Polynomial（多项式变换）；在PolynomialModelProperties中：PolynomialOrder：3控制点最低个数(t+1)(t+2)/2，本作业选择了1