整理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,利用遥感影像对地形图进行更新（四）采用软件：在本次实习中采用的软件是遥感影像处理软件ERDASI

2、MAGINE9.2。一、对于扫描地形图的影像纠正、裁切、拼接由于扫描地形图是利用已有的纸质地形图扫描而成的，在扫描过程中会存在着各种变形，所以在利用扫描地形图进行图像处理及数据采集之前需要对其进行纠正。纠正原理：通过图廓坐标对扫描地形图进行纠正。采用ERDAS中图像预处理模块中的几何纠正功能，在扫描地形图上选择一定数量的控制点，然后通过图上坐标判读，在控制点的参考坐标中输入读取的坐标值，并进行重采样，从而对扫描地形图进行纠正。（一）步骤：1,格式转换IMPORT模块将tif的地形图转换成为img格式。tif可以在ERDAS中直接打开，但转换格式之后可以使其操作方便。精品文档精品文档infoSH

3、Hiaoaoispotimg*141014.imgS42014.img®42015.img*142016.img44014.img44015.inng44016.imgaoispot.imgaoispot2.iiTigcail014imgcai20l4.img自1caiGOI5.imgjP|cai3016.img§cai4014.imgK山口xiao$F|Kiao43014img43015.img43018Jmgcai2015.img扈cai2016.imgBHcai3014.imgcai4015.img占1cai4U16.img占1mosaicl.img&pini

4、mgxiao41014.imgxiao42014.imgxiao42015.img%xiaom山口KiaoKiaoizuiz2,分别对九幅图做几何校正以为九幅地形图是纸质打印扫描出来的，纸质扫描和图纸时间久了出现图像误差故要对九幅图纸一一做几何坐标匹配。步骤如下设置：sessionpreferences-viewercleardisplay的钩去掉，使多幅图像可在一个窗口中打开。点击Viewers模块打开地形图；点Raster-GeometricCorrection在SetGeometricModel中，选Polynomial（多项式变换），点OK.在PolynomialModelProper

5、ties中：PolynomialOrder（多项式次数）设为2点Projection选项：MapUnits:Meters点Add/ChangeProjection,点Custom,进行设置ProjectionType（投影系统）：GaussKrugerSpheroidName（参考椭球）：IAG-75精品文档精品文档DatumName（基准面）：xian80Longitudeofcentralmeridian：108:00:00.000000E（1:10000图用的3度带，3*36=108,若为36度带）Latitudeoforiginofprojection:0:00:00.000000NF

6、alseeasting500000metresFalsenorthing:0.000000metres设置完后，点保存，同时也可以将自己设置的坐标投影信息另存为一个文件，方便每次调用。如果另存为文件可以在Standard中选择刚才保存的文件，确定。点SetProjectionformGCPTool,在GCPToolTeferenceSetup中选择KeyboardOnly,确定。在PolynomialModelProperties中点Apply,然后关闭在图上选择7个控制点，并根据经纬线输入其校正后x、y值，观察其精度情况.校正：OutputFiles:输入文件名ResampleMethodB

7、ilineaEnterpolation（双线性插值法）OutputCellSizes：X:0.85Y:0.85OK下图为具体步骤：步骤一精品文档精品文档步骤二步骤三aia.如1Ftrfcl：FwinAC*>-匐3TRYTot词-TRMitiHJR麻E减CM13疝h.111.翎口叼34W可卬耳21J04=EM?1qJ卬H21H1蟠NK畤13B3叼建1.54?CM.步骤四3,裁切、拼接九幅图在主窗口中选择DataPrepSubsetImages裁切地形图；在主窗口中选择DataPrepMosaicImages,EditAddImages装入九幅图Process-RunMosaic,输拼接后文

8、件名，开始拼接（结果如下图）4,用拼接后的地形图对SPOT影像进行几何校正精品文档精品文档:1:口：口：。:*:。:"：。：口：。：门二,:，：口:口LOX：!&：&?X'.：Q,二,工口：,：11-1"""i-i"""I:-:-:;:":-:-:-:”.:I:*:、，：；*:。:，:。.:*:“:*:*;*:*?:：*：'：'：j；：x：；：；,:1：11：:.；.,1:二：二：；:口：1二1：口:"：fJ："：,：口：1二1：口：口：J；：；：；.：；

9、：.：；：.：.：.：；：；：三：；：.：；：；：1Z3A567B91UOdeiFITImaqeName育3ReMHiDleRMSOnlneExcludeAreasIrnaaeDodqtdCdbfBakmcde：版卷币画果崛实习独据国安地图如口ErdireNh0.0000X5R将君1币授果桂酶习数疆碣声地圉痴EMiieNKI0.DOODXe小东君J币授深前虻白J朝也南安地图而切43bEdi七NHD0000X后饰都被课殛感实习敷捐酒安城图加MEntireNhlQODODXaf礴电利费集崛实习数据的安时图痴aM4口EhlireNNftOOOOX总瞬老师邮卷!幽句劭瑁牺友地囹布加打力ErfliieN

10、H.DODDXiolu.d东老师圈序M第曲rJ林捉/由宣地图描皿43NN0.0000X£iHIfa.I步骤二拼接后的地形图二、spot遥感影像几何校正1、设置：sessionpreferences-viewercleardisplay的钩去掉2、点Viewers,在Viewersl中打开遥感影像，再点击Viewers,在Viewer2中打开地形图；点sessiontileviewers,遥感影像和地形图并排显示；精品文档精品文档3、在遥感影像的视窗中，点Raster-GeometricCorrection在SetGeometricModel中，选Polynomial(多项式变换)；在

11、PolynomialModelProperties中:PolynomialOrder：3控制点最低个数(t+1)(t+2)/2,本作业选择了11个控制点其中1个作为多余检核。在PolynomialModelProperties中：点SetProjectionformGCPTool,为ExistingViewers,OK,然后在右侧地形图视窗点一下在PolynomialModelProperties中点Supply,再Close4、在GCPTool视窗中，点Edit,先SetPointTypecontrol可以改变控制点和参考点的颜色便于观察选区情况，点击color进行颜色选择。点的(Creat

12、eGCP),分别在遥感影像和地形图上选择同名点，点击一次，选取一对同名点，重复点击进行控制点选取；因为选择的是三次多项式，至少需要11个控制点。选女？11个点后，在GeoCorrectionTools中，点进行几何校正纠正精度：对于遥感影像，纠正的误差应该小于1像元1.5像元(图像分辨率为10米)，但为了有好的纠正精度，尽量控制TotalRMS在一个像元范围内，如果在选择完控制点后发现RMS过大，可以将单点中误差的点删去，来保证精度。5、Resample视窗中，OutputFiles:输入文件名ResampleMethodBilinearInterpolation(双线性才f值法)Output

13、CellSizes：X:10Y:10点击OK以下为实际操作步骤：步骤一精品文档精品文档步骤二步骤三6、打开校正好后的遥感影像。利用viewer中的Utility下的swipe命令可以把匹配好的spot影像与地形图嵌套来检查匹配的图上精度。精品文档精品文档注：左侧为校正前影像，右侧为校正后影像注：匹配好的spot影像与地形图嵌套。注：匹配好的spot影像与地形图嵌套来检查匹配的图上精度。精品文档精品文档三、地形图修测1、在Viewer#1中导入纠正好的Spot影像，导入纠正好的地形图。2、点击巾lefnew-vector.新建矢量图层。通过vector1001s来对新的地物进行描绘，从而达到更新

14、地形图的目地。3、以下绘制出了五处基本都在南部位置一条环城高速和四处大型建筑物原地形图没有的地物见下图标注。地形图上补测地物精品文档精品文档Spot上的现有地物四、精度分析作优:'向：|枷。闻疝叫：Rsfef/i魏上山长吗气跖小手FartlPWID>岫*1皿决Tft刈心'r'fettis那讪WbNsth1EL-S嗡闻£交翔麻西罐瞋圆Csnbi噂礴QE35Q相2EFl!碗战Mai姆吃您ConteiM曲OQ知Q256G0P*31架期4期安礴期潮即反非(jDtl碗郦Q期1珈4GOPM印®«3ES物即仃砌1孤0.3514)刑as4Q朝5GtP

15、tS阿1疆4H28K蝠55提CaKri败WQ/QE3?6ECPft愧感岫.夏礴道隔3JWSECMi0.5414)155闽。由G0P17咻那碎和讥邮烟CMriq库幅僻1循aGOPffi嘴捌-WB1I黝痂附附眼CWiQ3i涧Q到与9GtPft级用物覆3SJSI9&1M吸Csnbi幡:用。,中喻10GOO皿1期如.耦翻1利加触邮刚凶Q本uktmiOM5却一明3mlm地H策w蝴-D7F0ES31.H7ODO12HI21GridC+SPOT-2卫星是法国SPOT卫星的第二颗卫星，1990年1月发射，至2006年还在运行。星上载有两台完全相同的高分辨率可见光遥感器（HRV）,SPOT-2卫星的HR

16、V传感器有4个波段，在可见光和近红外波段成像，多光谱XS20米分辨率，B1,B2,B3。全色P,10米分辨率。因为此遥感影像的分辨率太低，其几何校正精度较低。所以在对遥感影像的纠正中，纠正的误差保持小于1像元1.5像元（图像分辨率为10米），但为了有好的纠正精度，尽量控制TotalRMS在一个像元范围内，其中，RMS误差精品文档精品文档（均方根）是指GCP的输入（原）位置和逆转换的位置之间的距离，它是在用转换矩阵对一个GCP作转换时所期望输出的坐标与实际输出的坐标之间的偏差。如果在选择完控制点后发现RMS过大，可以将单点中误差的点删去，来保证精度。下图为遥感影像纠正中的RMS值，都在精度范围之

17、内。从本次试验的数据得出的中误差都小于1,并且控制点对应X,Y分量的误差小于都1符合要求。此外，影响几何精校正的因素，主要表现在GCP的数量、分布和定位精度。以及校正方法（包括数学模型、确定亮度值的方法）不同，影像的纠正精度也不同。GCP的数量与分布。GCP数量的增加可以提高校正精度，但也使得寻找GCP的难度加大。GCP的数量也不宜过多，因为过多不会显著提高校正精度，却使计算量大大增加。几何精校正要求GCP均匀分布。若GCP分布不均匀，则会出现较大误差。本实验影像分辨率较低，选点时不太容易时控制点分散和均匀分布，致使图像精度一定程度受到影响。五、结论此次实习用的是1999年航空摄影测量方法制作的1:1万地形图的几何精度，并且应用2002年获取的10米分期率SPOT-2影像对1999年完成的1:10000地形图进行更新。运用遥感手段能够实时动态的更新地形图，对原有地形图修补来说遥感手段减少了人力、物力和财力，方便又不失准确性。从本次实验课以西安为例，从1999年到2002年西安作为一座国际性大城市三年的时间城市建设发生了很大变化。南部地区发展迅猛，城市建筑面积逐渐扩大。六、实验心得通过本次实验我学习到了erdas关于图像的转换格式，裁剪、拼接、几何校正等操作方法，明白了遥感在地