13像片纠正与正射影像图汇编.ppt

《 摄 影 测 量 》 像片纠正与正射影像图 数字微分纠正 一、像片纠正的概念 像片平面图或正射影像图是地图的一种，它是用像片上的影像表示地物的形状和平面位置的地图形式。利用中心投影的航摄像片编制像片平面图或正射影像图，是将中心投影转变为正射投影的问题。 当像片水平且地面为水平的情况下，航摄像片就相当于该地区比例尺为1：M(=f/H) 的平面（正射影像）图。 将竖直摄影的航摄像片通过投影变换，消除像片倾斜引起的像点位移，并限制或消除地形起伏引起的像点位移，获得相当于航摄像机物镜主光轴在铅垂位置摄影的水平像片，同时改化规定的成图比例尺，这种作业过程称为像片纠正。 主要是单张像片的像片纠正。 1、数字微分纠正的基本原理 数字纠正: 像素的几何位置和灰度 x = fx (X，Y) ； y = fy (X，Y) X = φx (x，y)； Y = φy (x，y) 反解法 数字微分纠正与光学微分纠正一样，是实现两个二维图像之间的几何变换 正解法 在已知像片内定向参数，外方位元素以及数字高程模型的前提下，按一定的数学模型用控制点解算 2、反解法（间接法）数字微分纠正 在已知像片内定向参数，外方位元素以及数字高程模型的前提下，可以进行数字微分纠正的反解算法，其步骤如下： （1）首先，按照图比例尺的要求将作业范围内的地面模型格网的间隔ΔX和ΔY的值重新确定，根据地面模型格网缩小图比例尺大小，即可得到纠正影像的格网。 （2）、计算地面点坐标 设正射影像上任意一点（象素中心）P的坐标为(X’,Y’)。由正射影像左下角图廓点坐标(X0,Y0)与正射影像比例尺分母M，计算P点所对应的地面坐标(X,Y)。 X ＝ X0 + M· X’ Y ＝ Y0 + M· Y’ ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? DEM (X0,Y0) ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （3）、计算像点坐标(x,y)或像点扫描坐标（I，J） DEM (X0,Y0) 应用共线方程求解原始图像上相应像点坐标(x,y)： 式中，Z是P点的高程，由DEM内插得到。需要注意的是这里算得的坐标不是以行列号表示的扫描坐标，因此还应由内定向将扫描坐标转换为原始图像的像点坐标。 在所有理论讨论中均取用像片坐标系中的像点坐标值x、y。 原始数字影像像元素的位置是以“扫描坐标系”行J（步进方向）、列I（扫描方向）为计量，为此，应利用像片坐标与扫描坐标之关系，求得相应的扫描坐标I、J，也可直接由X、Y、Z求解I、J。由于在象片扫描的数字化过程中，其扫描坐标系一般只是大致平行于象片坐标系，所以同一点的象片坐标x、y与其扫描坐标I、J不相等。由于影像数字化过程中可能产生的仿射变形，所以由下式 得到： 式中的系数L1、L2……L11是内定向参数m1、m2、n1、n2，像主点坐标I0、J0，旋转矩阵元素ai、bi、ci以及摄站坐标XS、YS、ZS的函数。 根据上式即可由纠正后影像坐标X、Y、Z直接获得数字化影像的像元素扫描坐标I、J。 （4）、灰度内插（插值） 由于所求得的原始图像上的像点坐标不一定正好落在其扫描采样的点子上，以此，这个像点的灰度值不能直接读出，必须进行灰度内插，一般采用双线性内插方法求得像点p的灰度值g(x,y)。 优点：由于此法是从输出的某一个有规律的数字地面模型的节点（X，Y）出发，反算其相应的输入影像点的点位（x，y），此时数字地面模型上的节点高程Z是已知的，故可以直接使用。 缺点：由于反算而得的在原始影像上的象元素，一般不会正好落在其扫描采样的点子上，则这点的灰度值不能直接读取，还需进行内插求得，这个过程叫重采样。最后进行赋值。 （5）、灰度赋值 最后将原始图像上像点p的灰度值赋给纠正后的正射影像上相应的像元素P， 即G(X,Y)=g(x,y) 依次对每个纠正像元素进行上述运算，即能获得纠正的数字图像，该方法称之为数字微分纠正的反解算法。 反解法解算流程 X Y 纠正影像 原始影像 x y 灰度内插 反算 3、正解法（直接法）数字微分纠正 原始影像 正射影像 Z? 对于单张未知外方位元素的航摄像片而言， a. 进行单张像片空间后方交会：即根据一张航摄像片上不在一条直线上的三个以上已知点的地面坐标，建立摄影中心、像方点、地物点三点共线的共线方程计算该像片6个外方位元素。 b. 量测像点坐标（x、y） c. 取一高程近似值Z。 d. 将（x、