实习二遥感图像的几何校正.doc

遥感数字图像处理试验报告二 （2011-2012第一学期） 题目：遥感图像的几何校正班级： 姓名： 学号： 实习二 遥感图像的几何校正1、 实验目的：通过实习操作，掌握遥感图像几何校正的基本方法和步骤，深刻理解遥感图像几何校正的意义。2、 实验数据美国Georgia州Atalanta的一景SPOT全色卫片；三幅哈萨克斯坦的LANDSAT（MSS和TM）影像wasial_mss.img、wasia2_mss.img、wasia3_mss.img；美国西雅图的影像seattl.img；美国加利福利亚州圣迭戈的一景Landsat TM影像dmtm.img3、 实验内容资源卫星数据的校正；图像的仿真变换；航片的正射校正；图像的镶嵌；图像投影变换等。4、 实验步骤2.1 资源卫星数据Landsat的校正以美国Georgia州Atalanta的一景SPOT全色卫片来对同一地区的Landsat影像做图像对图像的几何校正，SPOT卫片已经具有State Plane地图投影，基本步骤为：（1） 显示图像文件 1打开两个视窗Viewer#1和Viewer#2，在ERDAS图标面板菜单条中点击Session/Tile Viewers,将Viewer#1和Viewer#2平铺放置2在视窗Viewer#1中打开需要校正的Landsat TM影像tmAtlanta.img；Viewer#2中打开作为地理参考的校正过的SPOT影像：panAtlanta.img（2） 启动几何校正模块1在Viewer#1中点击Raster/Geometric Correction，选择多项式几何校正模型：Polynomial 点击OK，同时打开如下两个对话框：如图 2-1所示 图 2-12点击Appy，在点击Close，打开如下对话框（ 图2-2）： 图2-2 图2-3注意：该实例采样视窗采点模式，作为地理参考模式的SPOT图像已经含有投影信息，这里无需定义投影信息。如果不是采用视窗采点模式，或者参考图像没有包含投影信息，则必须在这里定义投影类型及其对应的投影参数。（3） 启动控制点工具如上对话框，点击OK，打开Viewer Selection Istructions指示器，在Viewer#2中点击左键，打开如图2-3所示的对话框，点击OK，整个屏幕自动变为两个主视窗、两个放大窗口（Viewer#2和Viewer#3）、两个关联窗口（分别位于两个视窗中，指示放大视窗与主视窗的关系）、控制点工具对话框、几何校正工具等，如图2-4，表面控制点工具被启动，进入控制点采集状态。 图2-4（4） 采集地面控制点在GCP工具启动时，系统缺省设置为自动编辑状态， 图标被激活。1 点击，进入GCP选择状态，在GCP CellArray中右键点击Color栏设置颜色。2 在Viewer#1中移动关联方框位置，寻找明显的地物特征点，如道路交叉处等作为输入GCP，然后点击，并在Viewer#3点击左键定点，GCP将记录刚输入的一个GCP，包括编号、标识码、X坐标、Y坐标。3 在Viewer#2中移动方框位置，寻找相应的地物特征点，作为参考GCP。然后点击，并在Viewer#4点击左键定点，系统将自动将参考点的坐标（X Reference,Y Reference）显示在GCP CellArray中。4 不断重复以上步骤，采用若干GCP，直到满足所选定的几何校正模型为止，而后每采集一个Input GCP,系统将自动生成一个Ref.GCP可以逐步优化校正模型。控制点的总数一般为方程 最高次方的46倍，要删去选中的GCP时，只需在GCP CellArray中右键单击Point#栏，选择Delete Selection即可。所有的输入GCP和参考GCP都可以直接保存在图像图像文件或参考图像文件中（Save Input / Save Reference），也可以保存在控制点文件中（Save Input As/ Save Reference As）便于以后调用，如果是保存在图像文件中，只要打开GCP工具，GCP点即会出现在视窗中；如果保存在GCP文件中，可以通过加载调用（Load Input/Load Reference）。（5） 采集地面检查点地面控制点用于控制点计算，建立转换模多项式方程，而地面检查点Check用于检验所建立的转换方程的精度和实用性。1 在GCP Tool菜单条单击Edit，选择Set Point Type/Check；GCP Tool/Edit/Point Mathing，在Point Mathing对话框中定义以下参数：确定匹配参数：最大搜索半径：3 搜索窗口大小：X：5 ，Y：5 确定约束参数：相关阈值：0.8 删除不匹配的点：选择Discard Unmatced Point复选框2 确定地面检查点：在GCP Tool工具条选择图标，并点击图标，锁住Create GCP功能，如同创建控制点那样，分别在Viewer#1和Viewer#2中定义5个检查点，定义完毕后点击unlock图标，解除Create GCP功能。3 计算检查点误差：在GCP Tool工具条中点击，检查点的误差就会显示在GCP Tool的上方，只有所有检查点的总误差小于一个像元（Pixel），才能继续进行合理的重采样。一般而言，如果控制点定位选择比较准确的话，检查点匹配会较好。ERDAS采用均方根RMS Error 来描述GCP点的单点误差和总误差。（6） 计算转换模型：点击，对话框中Parameters：改变多项式次数；Transformation:提供转换矩阵系数，并记录转换模型；Projection:查看和改变投影参数。（7） 图像重采样：像元属性值的重采样是依据未校正的图像的灰度值，采用某种方法估算校正后的图像像元灰度值的过程。ERDAS IMAGING提供三种重采样算法：最近邻插值（将最邻近像元值直接赋予输出像元）、双线性插值（用双线性方程和2*2窗口计算输出像元）、三次卷积插值（用三次方程和4*4窗口计算输出像元值）。 点击，打开Resample对话框，Resample Method：Bilinear Interpolation ；Output Cell Size：X：30、Y：30，设置输出统计中忽略零值：Ignore Zero in Stats.（8） 保存几何校正模块：定义模式文见（*.gms）（9） 检验校正结果：1右键Viewer#1，点击Geo.Link/Unlink,左键单击Viewer#2,建立与Viewer#1的连接。2右键Viewer#1，点击Inquire Cursor，打开光标查询对话框，移动光标，观测在两屏幕这噢噢那个的位置和匹配程度，并注意查询光标在对话框中的数据变化。2.2 图像的仿真变换图像的仿射变换就是对图像进行旋转、平移、翻转、拉伸等一次线性变换，以便图像的北方向真正朝上。1打开tmAtalanta.img,点击Raster/Geometric Correction/Affine，设置参数。2 点击，进行图像重采样操作。2.3图像镶嵌：图像镶嵌是将具有地理参考的若干相邻图像拼接成一幅图像或一组图像，需要拼接的输入图像必须含有地图投影信息，或者说输入图像必须经过几何校正处理或进行过校正标定。所以输入的图像可以具有不同的投影类型、不同像元大小，但必须具有系统的波段数。以三幅哈萨克斯坦的LANDSAT（MSS和TM）影像wasial_mss.img、wasia2_mss.img、wasia3_mss.img为例。（1） 显示图像在同一视窗中分别加载wasial_mss.img、wasia2_mss.img、wasia3_mss.img，并不选Clear Display,选择Background Transparent 和Fit to Fame.（2） 启动图像拼接工具：点击，选择Mosaic Images（3） 加载Mosaic图像： 点击，分别添加wasial_mss.img、wasia2_mss.img、wasia3_mss.img，并设置图像镶嵌区域：Compute Active Area（4） 图像重叠组合：点击选择与该模式对应的调整图像叠置次序的图标进行上下调整。（5） 图像色彩匹配设置：点击，将Matching Method：Overlap Area（匹配重叠区域），点击 /，设置参数：Intersection Type(相交类型)：Cutline Exists(已有裁切线)；Select Function(重叠区像元灰度值函数)：Cut/Feather(羽化)（6） 定义输出图像：点击/，打开Output Image Options对话框：Define Output Map Area(定义输出地图区域)：Union of All Inputs（7） 运行Mosaic工具：点击Process/Run Masaic，弹出的对话框设置输出图像区域：All2.4 图像重投影变换图像的重投影变换是将图像文件从一中地图投影类型转换到另外一种投影类型，一图像几何校正过程中的投影变换相比，这种直接的投影变换可以避免多项式近似值的拟合，对大范围图像的地理参考是非常有意义的。以美国西雅图的影像为例，进行图像重投影变换练习。（1） 启动图像重投影变换：点击，选择Reproject Images.（2） 图像重投影变换操作：设置参数如下图2-5： 图2-5 图2-6 2.5图像子集裁剪：（1） 矩形子集裁剪：矩形子集的裁剪是指裁剪图像的边界范围是一个矩形，通过左上角和右下角两点的坐标确定图像的裁剪位置。以美国加利福利亚州圣迭戈的一景Landsat TM影像dmtm.img为例介绍其操作方法。点击，选择Subset Image，定义参数，如图2-6： （2） 多边形子集裁剪：1AOI多边形裁剪2多边形Coverage子集裁剪按照行政区划或者自然区划边界进行图像的裁剪，往往是先利用ArcGIS或ERDAS的矢量模块绘制精确的边界多边形Polygon，然后以此Polygon为边界进行图像裁剪，过程：l 将多边形转换成栅格图像文件在ERDAS IMAGING图标面板菜单条中点击Main/Image Interpreter/Utilities/Vector to Raster,打开如图2-7所示对话框，并设置参数： 图2-7 图2-8l 然后通过掩膜运算实现图像不规则裁剪在ERDAS IMAGING图标面板菜单条中点击Main/Image Interpreter/Utilities/Mask,设置参数如图2-8，点击OK。校正结果对比： 校正前图像 校正后的图像 参考图像 5、 实习心得 由于遥感平台位置和运动状态的变化、地形起伏、地球表面曲率、大气折射、地球自转等因素的影响，遥感图像在几何位置上会发生变化，产生诸如行列不均匀，像元大小与地面大小对应不准确，地物形状不规则变化等畸变，称为遥感图像的几何畸变。产生的畸变给定量分析及位置配准造成困难，因此在遥感数据接收后需要对图像进行几何校准以使其接近真实的地理状况。 遥感影像相对于地图投影坐标系统进行配准校正，需要找到遥感影像与地图投影系统之间的数学函数关系通过这种函数关系可计算出原遥感影像中每个像元在地图投影坐标系统上位置从而得到校正后的图像。Erads软件中提供了7种几何校正模型，具体如下： 模型 功能 Affine图像彷射变换（不做投影变换） Polynomial多项式变换（同时做投影变换） Reproject投影变换（转换调用多项式变换） Rubber Sheeting非线性变换、非均匀变换 Camera航空影像正射校正在本次试验中采用Polynomial(多项式变换)的模型，通过在遥感影像和参考影像上分别选取相应的控制点，求出二元二次多项式函数得到变换后的图像坐标（x,y）与参考图像坐标的关系，从而对图像进行几何校正。 几何校正后，我们进一步判断所选控制点的准确性，一般要求RMS要小于0.5，若没有达到要求应重选点或者对所选点进行调整。在控制点采集过程中，随着控制点采集的完成，转换模型就自动生成，单击Geometric Correction Tools对话框中的图标，对画框中Transformatin可查阅多项式参数。如下图所示： 几何校正后，通过三幅图像的对比，可以发现校正后的图像相对于校正前的影像房屋、道路等影像信息发生了变化，像元的坐标发生了变化，而其变化后的地物特征更加接近参考图像的地物特征。通过本次上机实验及实验报告的撰写，对遥感图像几何校正的原理、方法和具体步骤有了一个清晰的认识，几何校正能在很大程度上对遥感影像在几何位置上发生的变化进行调整，减少地形起伏、地球曲率等因素对遥感影像的影响，从而为各种科学研究提供反应真是情况的遥感影像。在几何校正的过程中应注意操作的规范性和态度的严谨，尤其在控制点的采点过程中