遥感图像几何校正,介绍了几何误差源,并进行了校正

1、12l6.1 遥感图像坐标系统l6.2 遥感图像几何畸变l6.3 遥感图像几何纠正方法l6.4 软件实习31)遥感器坐标系统遥感器坐标系统S-UVWU轴:遥感器飞行方向V轴:垂直于U轴W轴:垂直于UV平面2) 地面坐标系统地面坐标系统O-XYZZ轴:原点处天顶方向XY平面垂至于Z轴3)图像坐标系统图像坐标系统o-xyfx y f 分别平行于UVW轴6.1 遥感图像坐标系统4WVUAZYXZYXPo地面点 P 在地面坐标中的坐标遥感器坐标系的原点 o 在地面坐标中的坐标P 点在遥感器坐标系中的坐标,不同的遥感器不同321321321cccbbbaaaA6.1 遥感图像坐标系统(续1)5sinsi

2、nsincoscos1acossinsinsincos2acossin3asincos1bcoscos2 bsin3bsinsincoscossin1ccossincossinsin2ccoscos3 c6.1 遥感图像坐标系统(续2)6l6.2.1遥感器本身引起的畸变 l6.2.2外部因素引起的畸变 l6.2.3处理过程中引起的畸变 7l遥感器本身引起的几何畸变与遥感器的结构、特性和工作方式不同而异。这些因素主要包括：l1） 透镜的辐射方向畸变像差；l2） 透镜的切线方向畸变像差；l3） 透镜的焦距误差；l4） 透镜的光轴与投影面不正交；l5） 图像的投影面非平面；l6） 探测元件排列不整齐

3、；l7） 采样速率的变化；l8） 采样时刻的偏差;l9) 扫描镜的扫描速度变化 。8例如扫描形式成像的MSS，产生的几何畸变主要是由于扫描镜的非线性振动和其它一些偶然因素引起的。在地面上影响可达395米。 全景畸变：9l影响图像变形的外部因素包括：l1） 地球的曲率l2） 大气密度差引起的折光l3） 地形起伏l4） 地球自传l5） 遥感器轨道位置和姿态等 101112中心投影13多中心投影例如MSSTM等14l遥感图像再处理过程中产生的误差，主要是由于处理设备产生的噪声引起的。l传输、复制、l 光学l 数字 15l遥感图像的几何粗处理和精处理。l遥感图像的几何纠正按照处理方式分为光学纠正和数字

4、纠正。l光学纠正主要用于早期的遥感图像的处理中，现在的应用已经不多。除了对框幅式的航空照片（中心投影）可以进行比较严密的纠正以外，对于大多数动态获得的遥感影像只能进行近似的纠正。l主要介绍数字图像的几何纠正。16 数字图像几何纠正：通过计算机对离散结构的数字图像中的每一个像元逐个进行纠正处理的方法。 这种方法能够精确地改正动态扫描图像的误差。 基本原理：利用图像坐标和地面坐标（另一图像坐标、地图坐标等）之间的数学关系，即输入图像和输出图像间的坐标转换关系实现。 17),(ppXpyxFX ),(ppYpyxFY ),(PPxpYXfx ),(PPypYXfy 其中，（xp ,yp）（XP,YP

5、）分别是任意一个像元在原始图像和纠正后图像中的坐标。直接间接18直接纠正方法直接纠正方法:从原始图像，依次对每个像元根据变换函数 F（），求得它在新图像中的位置。并将灰度值付给新图像的对应位置上。间接纠正法间接纠正法: :6-3和 6-4是反解变换公式。从新图像中依次每个像元，根据变换函数 f () 找到它在原始图像中的位置，并将图像的灰度值赋予新图像的像元。19l纠正后图像和原始图像的形状、大小、方向都不一样。所以在纠正过程的实施之前，必须首先确定新图像的大小范围。l根据公式61，62求出原始图像四个角点（a, b, c, d）在纠正后图像中的对应点（a, b, c, d）的坐标（Xa,Ya

6、）(Xb,Yb) (Xc,Yc) (Xd,Yd)，l然后求出最大值和最小值。 20lX1 = min (Xa, Xb, Xc, Xd)lX2 = max (Xa, Xb, Xc, Xd)lY1 = min (Ya, Yb, Yc, YXd)lY2 = max (Ya, Yb,Yc, Yd)21l目的是确定新图像宽度和高度；l根据精度要求，在新图像的范围内，划分网格，每个网格点就是一个像元。l新图像的行数 M(Y2-Y1)/Y+1;l新图像的列数 N(X2-X1)/X+1;l新图像的任意一个像元的坐标由它的行列号唯一确定。22l纠正后的新图像的每一个像元，根据变换函数，可以得到它在原始图像上的位

7、置。如果求得的位置为整数，则该位置处的像元灰度就是新图像的灰度值。l如果位置不为整数，则有几种方法：l1) 最近邻法l2) 双线性内插法l3）三次卷积法 231) 最近邻法：距离实际位置最近的像元的灰度值作为输出图像像元的灰度值；242）双线性法：以实际位置临近的4个像元值，确定输出像元的灰度值。公式为： 4141432144332211,iiiiinmpgpppppgpgpgpgpg）（ G2 G3 G4 G1 253）三次卷积法以实际位置临近的16个像元值，确定输出像元的灰度值。公式为： 161161,iiiiinmpgpg）（ PSF 三次样条函数 sinc函数26几种采样方法的优缺点：

8、1）最近邻法：计算简单，一般情况下可接受，灰度级细微变化，保留原始数据。2) 双线性插值：计算，效果，细节丧失；斜率27总上所述： 纠正的函数可有多种选择：多项式方法、共线方程方法、随机场内插方法等等。其中多项式方法的应用最为普遍。 准准备备 工工作作 输输入入原原始始图图象象 建建立立纠纠正正函函数数 确确定定输输出出图图象象的的范范围围 逐逐个个像像元元进进行行几几何何变变化化 灰灰度度的的重重采采样样 输输出出纠纠正正后后的的图图象象 效效果果 评评价价 28l多项式纠正的基本思想：图像的变性规律可以看作是平移、缩放、旋转、仿射、偏扭、弯曲等形变的合成。一般的公式为： .54321022

9、),(iiiiiiiixiYcXcYXcYcXccYXfx.54321022),(iiiiiiiiyiYdXdYXdYdXddYXfy.54321022),(iiiiiiiixiycxcyxcycxccyxFX.54321022),(iiiiiiiiyiydxdyxdydxddyxFY2) 直接法 1) 间接法 利用有限的控制点的已知坐标，解求多项式的系数，确定变换函数。然后将各个像元带入多项式进行计算，得到纠正后的坐标。291) 表征空间位置的可靠性，道路交叉点，标志物，水域的边界，山顶，小岛中心，机场等。同名控制点要在图像上均匀分布；清楚辨认；数量应当超过多项式系数的个数（(n+1)*(n

10、+2)/2）。当控制点的个数超过多项式的系数个数时，采用最小2乘法进行系数的确定，使得到的系数最佳。 （续1）控制点的选择原则:ENVI: (degree + 1)230l最小二乘法最早称为回归分析法。由著名的英国生物学家、统计学家道尔顿（F.Gallton）所创。 ( 道尔顿研究英国男子中父亲们的身高与儿子们的身高之间的关系时，创立了回归分析法。 )l探索变量之间关系最重要的方法，用以找出变量之间关系的具体表现形式。l后来，回归分析法从其方法的数学原理误差平方和最小（二乘是平方的意思）出发，改称为最小二乘法。 31l1. 为了精确地描述Y与X之间的关系，必须使用这两个变量的所有观察值，才不至

11、于以“点”概面（作到同步与全面）。l2Y与X之间是否是直线关系（用协方差或相关系数衡量）？若是，将用一条直线描述它们之间的关系。l3在Y与X的散点图上画出直线的方法很多。l任务？找出一条能够最好地描述Y与X（代表所有点）之间关系的直线。l4什么是最好？找出判断“最好”的原则。l最好指的是找这么一条直线，使得所有点到该直线的纵向距离的和（平方和）最小。l5三种距离32点到直线的距离点到直线的垂直线的长度。横向距离点沿（平行）X轴方向到直线的距离。纵向距离点沿（平行）Y轴方向到直线的距离。也就是实际观点的Y坐标减去根据直线方程计算出来的Y的拟合值。33l纵向距离是Y的实际值与拟合值之差，差异大拟合

12、不好，差异小拟合好。l将所有纵向距离平方后相加，即得误差平方和。所以，“最好”直线就是使误差平方和最小的直线。l运用求极值的原理，将求最好拟合直线问题转换为求误差平方和最小的问题。3435最小二乘法的性质:之一：即拟合直线过y和x的平均数点。之二： 残差与自变量x的乘积和等于0，即残差与自变量x两者不相关。36l6.3.1 图像到图像的几何纠正（配准）l6.3.2 图像的拚接(图像的镶嵌)37预备知识：打开两个图像文件，其中： BASE IMAGE： 输出图像的参照图像。 WARP IMAGE：原始图像，要被纠正的图像 GCP: Ground Control Points, 文件类型为 *.p

13、ts38l图像拚接(镶嵌MOSAIK)是将多个具有重叠部分的图像制作成一个没有重叠的新图像。l有基于像元的拚接和基于地理坐标的拚接。l本节介绍基于像元的图像拚接。 39l准备工作：l1） 要参加拚接的图像必须具有统一的坐标系，即首先进行图像的几何纠正。l2） 图像灰度的调整；l对于彩色图像，需要从红绿蓝三个波段分别进行灰度的调整；l对于多波段的图像文件，进行一一对应的多个波段的灰度调整。l灰度调整的方法：进行交互式的图像拉伸，进行图像直方图的规定化(匹配)，或者进行更加复杂的类似变化（如线性回归）。40l确定输出图像的大小；l确定参加图像拚接的方法；有无羽化（Feathering）。 参加拚接

14、的图像如没有羽化，则用上面的图像灰度值直接替代重复部分的灰度值；l如果有羽化，则重复部分的灰度值由两个图像的灰度值确定。l羽化方法：一是边缘羽化，二是用户自定义羽化。 图像拚接的操作：411） 边缘羽化边缘羽化 用户需要键入羽化的距离。羽化部分的灰度值采取两个图像的线性加权值。边缘处Image #2 0%, Image #1 100%, Averaging Distance 处， Image#2占 100%, Image1 占0； 中间部位线性组合。422) 自定义羽化自定义羽化首先必须在图像上使用ANN定义Cutline.灰度值的计算方法：Cutline处，Image #2 占100%, Image #1 占0%, Averaging Distance 处， Image#2占0%, Image1 占100；中间地带，线性插值。43 CUT LINE 是由一根从图像一个边缘到另一个边缘的POLYLINE，以及要切掉部分上的一个SYMBOL定义的。 CUT LINE的定义步骤为：1)打开IMAGE窗口上的菜单命令 FuctionsOverlaysAnnotation，则打开Annotation窗口；2) 选择要在那个图像窗口上画ANN, IMAGE? SCROLL? ZOOM? 3) 在窗口上，选择菜单命令ObjectPolyline；4) 使用鼠标画线，按右键