第五章摄影测量解析基础.pptx

《第五章摄影测量解析基础.pptx》由会员分享，可在线阅读，更多相关《第五章摄影测量解析基础.pptx（106页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、5.1 5.1 像点坐标量测像点坐标量测5.2 5.2 单像空间后方交会单像空间后方交会5.3 5.3 立体像对的前方交会立体像对的前方交会5.4 5.4 立体像对的解析法相对定向立体像对的解析法相对定向5.5 5.5 立体模型的解析法绝对定向立体模型的解析法绝对定向5.6 5.6 双像解析的光束法严密解双像解析的光束法严密解第1页/共106页5.1 像点坐标量测第2页/共106页用于量测像点坐标的仪器称为用于量测像点坐标的仪器称为立体坐标量测仪立体坐标量测仪。分为。分为普通普通立体坐标量测仪立体坐标量测仪和和精密精密立体坐标量测仪立体坐标量测仪。量测精度达到量测精度达到3-53-5umum的

2、立体坐标量测仪称为的立体坐标量测仪称为精密精密立体坐标量测仪。立体坐标量测仪。量测的成果：量测的成果：u左片上的像点坐标（左片上的像点坐标（x,yx,y）、右片上同名）、右片上同名点的点的左右视差左右视差p p及及上下视差上下视差q qu 左右片上同名像点的坐标（左右片上同名像点的坐标（x x1 1,y,y1 1）、）、（x x2 2,y,y2 2）第3页/共106页左右视差：同名像点在各自的像平面坐标系的x坐标之差上下视差：同名像点在各自的像平面坐标系的y坐标之差常用的仪器：常用的仪器：steko 1818steko 1818、HCT-1HCT-1（20 20 ）PSK2PSK2（1 1 ）

3、等）等第4页/共106页Steko 1818 型立体坐标量测仪X手轮x读数鼓Y手轮上下视差环左右视差手轮上下视差（q）读数鼓左右视差（p）读数鼓第5页/共106页u 准备工作准备工作：标出像片坐标系，标出要量测的：标出像片坐标系，标出要量测的像点位置（定向点、加密点、外业控制点）并像点位置（定向点、加密点、外业控制点）并标明点号标明点号u 像片归心像片归心：使像片坐标系的：使像片坐标系的原点原点位于仪器坐位于仪器坐标系的已知位置（标系的已知位置（框标连线框标连线交点交点与旋转与旋转中心中心重合重合）u 像片定向像片定向：像片平面坐标系与仪器坐标系的：像片平面坐标系与仪器坐标系的轴系轴系平行平行

4、,记下仪器的记下仪器的X X0 0、Y Y0 0、P P0 0、Q Q0 0。立体坐标量测立体坐标量测步骤步骤第6页/共106页u 像点坐标量测：像点坐标量测：利用利用X X、Y Y、P P、Q Q手轮的转动，使测手轮的转动，使测标立体切准待量测的点，并记下读数。标立体切准待量测的点，并记下读数。左片上像点坐标：左片上像点坐标：右片上像点坐标：右片上像点坐标：第7页/共106页PSK-2精密立体坐标量测仪第8页/共106页BC2解析测图仪第9页/共106页获取获取六个外方位元素六个外方位元素的方法的方法?雷达、雷达、GPSGPS、INSINS、星象相机、星象相机地面控制点反算（地面控制点反算（

5、单像空间后方交会单像空间后方交会）POS系统5.2 单像空间后方交会第10页/共106页单像空间后方交会：单像空间后方交会：利用航片上利用航片上三个三个以上的以上的像点坐标及对应的地面控制点坐标，计算像像点坐标及对应的地面控制点坐标，计算像片片外方位元素外方位元素的工作。的工作。XYZaxyzs(Xs,Ys,Zs)ACBbc第11页/共106页后方交会摄影测量的后方交会像片的外方位元素第12页/共106页第13页/共106页第14页/共106页第15页/共106页空间后方交会的基础方程：共线方程第16页/共106页第17页/共106页共线条件方程式共线条件方程式二、基本关系式二、基本关系式u已

6、知值已知值 ：x x0 0,y,y0 0,f,m,f,m,X,Y,ZX,Y,Z（控制点）（控制点）u观测值：观测值：u待求：待求：Xs,Ys,Zs,Xs,Ys,Zs,(由外方位角元素由外方位角元素 ,确定）确定）像片的外方位元像片的外方位元素素非线性函数模型非线性函数模型，线性化，线性化第18页/共106页线性化：按泰勒公式展开，取小值一次项线性化：按泰勒公式展开，取小值一次项偏导数，系数偏导数，系数外方位元素初始值的改正数外方位元素初始值的改正数像点坐标近似值，将外方位元素的初始值代入共线条件方程的计算值像点坐标近似值，将外方位元素的初始值代入共线条件方程的计算值 一个控制点可以列两个方程，

7、至少要一个控制点可以列两个方程，至少要三个控制点三个控制点解六个外方位解六个外方位元素，有多余观测用平差的方法计算元素，有多余观测用平差的方法计算第19页/共106页三、误差方程式和法方程式三、误差方程式和法方程式(通常在像片的四个角上选取通常在像片的四个角上选取四个地四个地面控制点面控制点)观测值：观测值：像点坐标像点坐标第20页/共106页第21页/共106页第22页/共106页第23页/共106页 在竖直摄影的情况下，角元素都很小在竖直摄影的情况下，角元素都很小（3度）可以近似地用度）可以近似地用 及及 ，各系数可简化为：，各系数可简化为：第24页/共106页N个点的误差方程式的矩阵形式

8、：个点的误差方程式的矩阵形式：第25页/共106页法方程式（最小二乘原理）法方程式（最小二乘原理）逐步趋近计算逐步趋近计算怎样进怎样进行？行？第26页/共106页第27页/共106页第28页/共106页1 1、获取已知数据获取已知数据:包括平均航高，内方位元素，从外业包括平均航高，内方位元素，从外业测量成果中，获取控制点的测量成果中，获取控制点的地面测量坐标地面测量坐标，并转化成地，并转化成地面面摄影测量摄影测量坐标。坐标。2 2、量测控制点的像点坐标：量测控制点的像点坐标：将控制点刺在像片上，利用将控制点刺在像片上，利用立体坐标量测仪量测控制点在框标坐标系中的坐标。并立体坐标量测仪量测控制点

9、在框标坐标系中的坐标。并转化成以像主点为坐标原点的坐标。转化成以像主点为坐标原点的坐标。3 3、确定未知数的初始值确定未知数的初始值 ：在竖直摄影情况下，三个角在竖直摄影情况下，三个角元素的初始值为元素的初始值为0 0，即，即 线元素为线元素为空间后方交会的步骤空间后方交会的步骤第29页/共106页4 4 4 4、计算旋转矩阵计算旋转矩阵计算旋转矩阵计算旋转矩阵R R R R：利用角元素的近似值按计算方利用角元素的近似值按计算方利用角元素的近似值按计算方利用角元素的近似值按计算方向余弦，组成旋转矩阵。向余弦，组成旋转矩阵。向余弦，组成旋转矩阵。向余弦，组成旋转矩阵。5 5 5 5、逐点计算像点

10、坐标近似值：逐点计算像点坐标近似值：逐点计算像点坐标近似值：逐点计算像点坐标近似值：利用未知数的近似利用未知数的近似利用未知数的近似利用未知数的近似值代入共线方程，计算控制点像点坐标的近似值值代入共线方程，计算控制点像点坐标的近似值值代入共线方程，计算控制点像点坐标的近似值值代入共线方程，计算控制点像点坐标的近似值（x x）（）（）（）（y y）6 6 6 6、组成误差方程式：组成误差方程式：组成误差方程式：组成误差方程式：按公式组成误差方程式，然按公式组成误差方程式，然按公式组成误差方程式，然按公式组成误差方程式，然后组成法方程式，解算未知数的改正数；后组成法方程式，解算未知数的改正数；后组

11、成法方程式，解算未知数的改正数；后组成法方程式，解算未知数的改正数；7 7 7 7、改正数小于指定值，则完成；否则将解算的未、改正数小于指定值，则完成；否则将解算的未、改正数小于指定值，则完成；否则将解算的未、改正数小于指定值，则完成；否则将解算的未知数加上初始值，作为新的初始值，重复知数加上初始值，作为新的初始值，重复知数加上初始值，作为新的初始值，重复知数加上初始值，作为新的初始值，重复4-64-64-64-6步。步。步。步。第30页/共106页四、空间后方交会的解算过程获取已知数据(像片比例尺、航高、内方位元素、控制点的物方坐标)量测控制点的像平面坐标系坐标确定未知数的初始值 计算旋转矩

12、阵R R 计算像点坐标的近似值(x),(y)组成误差方程式和法方程式:解求外方位元素结束是否改正数小于限差否？第31页/共106页五、空间后方交会的精度未知数的协因数阵：未知数的中误差：六、空间后方交会的不定性：控制点不能位于同一个圆柱面上，否则解不唯一后方交会像片的外方位元素解求相应地面点的坐标一张像片？第32页/共106页5.3 5.3 立体像对的前方交会立体像对的前方交会一、立体像对前方交会的概念一、立体像对前方交会的概念对单张像片而言，知道了该像片的外方位元素后，对单张像片而言，知道了该像片的外方位元素后，不能由像片上的像点坐标求相应的物点坐标。不能由像片上的像点坐标求相应的物点坐标。

13、但如果知道了一对立体像对的外方位元素，则由立但如果知道了一对立体像对的外方位元素，则由立体像对上的体像对上的同名像点同名像点的坐标可以求出该对同名像点的坐标可以求出该对同名像点相对应的物点坐标。相对应的物点坐标。第33页/共106页第34页/共106页由立体像对由立体像对中两张像片中两张像片的内、外方的内、外方位元素和像位元素和像点坐标来确点坐标来确定相应地面定相应地面点在物方空点在物方空间坐标系中间坐标系中坐标的方法坐标的方法XYZa1(x1,y1)x1y1z1S1A(X,Y,Z)a2(x2,y2)z2y2x2S2立体像对前方交会的定义立体像对前方交会的定义第35页/共106页第36页/共1

14、06页式中式中 （X,Y,Z）物点的物方坐标（通常为地面物点的物方坐标（通常为地面摄测坐标）摄测坐标）投影中心的物方坐标投影中心的物方坐标 像点的像空间辅助坐标像点的像空间辅助坐标 N 投影系数投影系数投影系数的计算式投影系数的计算式:二、空间前方交会基本关系式二、空间前方交会基本关系式第37页/共106页地面点地面点A A在左右像片上的构像分别为在左右像片上的构像分别为a a1 1和和a a2 2A A在在D-XYZD-XYZ中的坐标为：中的坐标为：（X XA A，Y YA A，Z ZA A）S S1 1在在D-XYZD-XYZ中的坐标为：中的坐标为：（X XS1S1，Y YS1S1，Z Z

15、S1S1）S S2 2在在D-XYZD-XYZ中的坐标为：中的坐标为：（X XS2S2，Y YS2S2，Z ZS2S2）U1U2V2W2W1V1第38页/共106页摄影基线摄影基线B B u=Xs2 Xs1B v=Ys2 Ys1B w=Zs2 Zs1s1s2第39页/共106页a a1 1的像空间坐标为（的像空间坐标为（x x1 1，y y1 1，-f-f），像空间辅助坐标即在），像空间辅助坐标即在S S1 1-U-U1 1V V1 1W W1 1中的坐标为（中的坐标为（u u1 1，v v1 1，w w1 1）a a2 2的像空间坐标为（的像空间坐标为（x x2 2，y y2 2，-f-f）

16、，像空间辅助坐标即在），像空间辅助坐标即在S S2 2-U-U2 2V V2 2W W2 2中的坐标（中的坐标（u u2 2，v v2 2，w w2 2）R R1 1，R R2 2为由已知的外方位角元素计算的左、右像片的为由已知的外方位角元素计算的左、右像片的旋转矩阵。旋转矩阵。第40页/共106页根据三点共线有：根据三点共线有：点投影系点投影系数数点投影系点投影系数数Aa1s1v1u1w1YA-Ys1XA-Xs1ZA-Zs1a2s2A其中其中N N1 1称为左投影系数，称为左投影系数，N N2 2为右投影系数。为右投影系数。第41页/共106页变形为：变形为：根据基线分量：由于外方位元素已知

17、，因此是已知的根据基线分量：由于外方位元素已知，因此是已知的Y Y坐标应取平均值：坐标应取平均值：第42页/共106页点投影法前方交会点投影法前方交会(1)、(3)式联立求解第43页/共106页前方交会的基本步骤及计算公式：前方交会的基本步骤及计算公式：1 1、由已知外方位角元素及像点坐、由已知外方位角元素及像点坐标计算像点在像空间辅助坐标系标计算像点在像空间辅助坐标系中的坐标。中的坐标。2 2、由外方位线元素计算投影基线分量。、由外方位线元素计算投影基线分量。3 3、计算点投影系数。、计算点投影系数。4 4、计算地面点的坐标。、计算地面点的坐标。第44页/共106页三、双像解析的空间后方与前

18、交方法三、双像解析的空间后方与前交方法一、野外像片控制测量一、野外像片控制测量二、用立体坐标量测仪量测像二、用立体坐标量测仪量测像点坐标（左、右像片上）点坐标（左、右像片上）三、空间后方交会计算两张像片外方位元素三、空间后方交会计算两张像片外方位元素四、空间前方交会计算待定点地面坐标四、空间前方交会计算待定点地面坐标第45页/共106页像片控制测量获得像控点物方坐标像片控制测量获得像控点物方坐标像点坐标量测，得（像点坐标量测，得（x x1,y,y1 1,x,x2 2,y,y2 2）i;(xi;(x1 1,y,y1 1,x,x2 2,y,y2 2)j;)j;单像空间后方交会，分别得左右像片的外方

19、位元素单像空间后方交会，分别得左右像片的外方位元素 双像空间前方交会得待定点物方坐标双像空间前方交会得待定点物方坐标（X,Y,ZX,Y,Z）j j单像空间后交+双像空间前交作业流程第46页/共106页5.4 5.4 立体像对的解析法相对定向立体像对的解析法相对定向一、解析法相对定向的概念一、解析法相对定向的概念 用于描述两张像片相对位置和姿态的参数称为相对定向元素。用解析用于描述两张像片相对位置和姿态的参数称为相对定向元素。用解析计算的方法解求相对定向元素的过程称为解析相对定向从而建立地面立体模计算的方法解求相对定向元素的过程称为解析相对定向从而建立地面立体模型。型。第47页/共106页 相对

20、定向元素是描述立体像对中两张像片的相对定向元素是描述立体像对中两张像片的相对位置和姿态关系的元素，因此，可以把两相对位置和姿态关系的元素，因此，可以把两张像片各自相对于选定的同一个像空间辅助坐张像片各自相对于选定的同一个像空间辅助坐标系来讨论相对定向元素。为便于讨论，仿照标系来讨论相对定向元素。为便于讨论，仿照外方位元素的定义，引入外方位元素的定义，引入“相对方位元素相对方位元素概概念将像片在选定的像空间辅助坐标系中的位置念将像片在选定的像空间辅助坐标系中的位置（摄影中心（摄影中心S S的坐标，用的坐标，用 表示表示 ，和姿态（像片的姿态角，用，和姿态（像片的姿态角，用 表示）定义为像片的相对

21、方位元素坐标系的表示）定义为像片的相对方位元素坐标系的选择通常有两种形式：连续像对相对定向坐标选择通常有两种形式：连续像对相对定向坐标系和单独像对相对定坐标系，相应的相对定向系和单独像对相对定坐标系，相应的相对定向元素分为连续像对相对定向元素和单独像对相元素分为连续像对相对定向元素和单独像对相对定向元素。对定向元素。第48页/共106页第49页/共106页第50页/共106页第51页/共106页 1 1 连续像对相对定向元素连续像对相对定向元素 连续像对相对定向是以左方像片为基准，求出连续像对相对定向是以左方像片为基准，求出右方像片相对于左方像片的相对方位元素、选定右方像片相对于左方像片的相对

22、方位元素、选定像空间辅助坐标系像空间辅助坐标系S S1 1U U1 1V V1 1W W1 1使得左像片在使得左像片在S S1 1U U1 1V V1 1W W1 1中的相对方位元素均为已知值。为简便讨论，中的相对方位元素均为已知值。为简便讨论，以左像片的像空间坐标系作为像空间辅助坐标系以左像片的像空间坐标系作为像空间辅助坐标系如下图。此时，左、右像片的相对方位元素为：如下图。此时，左、右像片的相对方位元素为：左像片：左像片：右像片：右像片：由于由于 ，因此，相对定向需要解求的元素只，因此，相对定向需要解求的元素只有有5 5个，即个，即 称为连续像对相对定向称为连续像对相对定向元素。元素。第5

23、2页/共106页第53页/共106页2 2 单独像对相对定向元素单独像对相对定向元素 单独像对相对定向是以摄影基线作为像空间辅单独像对相对定向是以摄影基线作为像空间辅助坐标系的助坐标系的X X轴，以左摄影中心轴，以左摄影中心S S1 1为原点，左像片为原点，左像片主光轴与摄影基线主光轴与摄影基线B B组成的主核面为组成的主核面为UZUZ平面，构成平面，构成右手直角坐标系右手直角坐标系S S1 1U U1 1V V1 1W W1 1。此时，左、右像片的。此时，左、右像片的相对方位元素为：相对方位元素为：左像片：左像片：右像片：右像片：由于由于 只影响相对定向后建立的模型大小，只影响相对定向后建立

24、的模型大小，而不影响模型的建立，因此，相对定向需要解求而不影响模型的建立，因此，相对定向需要解求的元素只有的元素只有5 5个，即个，即 称称为单独像对相对定向元素。为单独像对相对定向元素。第54页/共106页第55页/共106页二、解析法相对定向原理二、解析法相对定向原理 从两个摄站对同一地面摄取一个立体像对时，同名射从两个摄站对同一地面摄取一个立体像对时，同名射线对对相交于地面点，见下图，此时，若保持两张像片线对对相交于地面点，见下图，此时，若保持两张像片之间相对位置和姿态关系不变将两张橡片整体移动时，之间相对位置和姿态关系不变将两张橡片整体移动时，同名射线对对相交的特性也不发生变化。反过来

25、，若完同名射线对对相交的特性也不发生变化。反过来，若完成了相对定向，恢复两张像片的相对定向元素，就能实成了相对定向，恢复两张像片的相对定向元素，就能实现同名射线对对相交，建立相对立体模型。因此，同名现同名射线对对相交，建立相对立体模型。因此，同名射线对对相交是相对定向的理论基础。射线对对相交是相对定向的理论基础。第56页/共106页二、解析法相对定向的共面条件二、解析法相对定向的共面条件 如图所示，如图所示，和和 为一对同名射线。其矢量用为一对同名射线。其矢量用 和和 表示，摄影基线矢量用表示，摄影基线矢量用 表示。同名射表示。同名射线对对相交，表明射线线对对相交，表明射线 位位于同一平面内，

26、亦即三矢量共面。根据矢量代数，三于同一平面内，亦即三矢量共面。根据矢量代数，三矢量共面，它们的混合积等于零，即矢量共面，它们的混合积等于零，即 上式即为共面条件方程，其值为零的条件是完成上式即为共面条件方程，其值为零的条件是完成相对定向的标准，用于解求相对定向元素。相对定向的标准，用于解求相对定向元素。第57页/共106页解析法相对定向的共面条件解析法相对定向的共面条件1 1、理论基础：同名射线对对相交、理论基础：同名射线对对相交2 2、共面条件、共面条件改成坐标形式为：改成坐标形式为：其中：其中：第58页/共106页三、解求相对定向元素的关系式三、解求相对定向元素的关系式 1.1.连续像对相

27、对定向元素解求式连续像对相对定向元素解求式 连续像对相对定向是以左像片为基础，求出右像片相对于左像片的五个相对连续像对相对定向是以左像片为基础，求出右像片相对于左像片的五个相对定向元素。定向元素。2.2.单独像对相对定向元素解求式单独像对相对定向元素解求式第59页/共106页1.1.连续像对相对定向连续像对相对定向 连续像对相对定向是以连续像对相对定向是以左像片左像片为基准，求出右像片为基准，求出右像片相对于左像片的相对于左像片的5 5个定向元素，即个定向元素，即 以左像片的像空间坐标系作为像空间辅助坐标系，以左像片的像空间坐标系作为像空间辅助坐标系，记为记为 S S1 1U U1 1V V1

28、 1W W1 1，过右摄影中心作另一像空间辅助坐，过右摄影中心作另一像空间辅助坐标系标系S S2 2U U2 2V V2 2W W2 2 ，两者的相应坐标轴相互，两者的相应坐标轴相互平行平行。此时，。此时，a a1 1,a,a2 2在各自的像空间辅助坐标系中的坐标分别为在各自的像空间辅助坐标系中的坐标分别为：（u1,v1,w1）和（）和（u2,v2,w2）,S S2 2在在 S S1 1U U1 1V V1 1W W1 1 中的坐中的坐标为标为 ,则上述共面条件方程式可以用坐标则上述共面条件方程式可以用坐标表示为；表示为；第60页/共106页 其中，其中，R R是是右像右像片片相对于相对于像空

29、间辅像空间辅助坐标系助坐标系的三个角的三个角元素的元素的 函数由于函数由于 只只涉及模型比例尺，涉及模型比例尺，相对定向中可给予相对定向中可给予定值，因此，只有定值，因此，只有5 5个相对定向元素个相对定向元素 为未知值。为未知值。W2W1V2U2V1U1第61页/共106页Bs1s2bwbvbu为了统一，将为了统一，将b bv v,b,bw w也化为用角度表示，由图可知也化为用角度表示，由图可知第62页/共106页进而，共面条件方程变为：进而，共面条件方程变为：上式是非线性函数，按照泰勒级数展开，取到一次项，上式是非线性函数，按照泰勒级数展开，取到一次项，使之线性化。使之线性化。F F0 0

30、为用未知数为用未知数(相对定向元素）的相对定向元素）的近似值近似值及给定的及给定的b bu u代代入计算出的。入计算出的。第63页/共106页要求出上式中的偏导数，必须先求得偏导数要求出上式中的偏导数，必须先求得偏导数因推导过程只考虑一次项，所以坐标变换矩阵可以用旋转矩阵的因推导过程只考虑一次项，所以坐标变换矩阵可以用旋转矩阵的小值一次项来表示：小值一次项来表示：对上式分别对对上式分别对 求导数，得：求导数，得：第64页/共106页由上式则可以求得各个偏导数：由上式则可以求得各个偏导数：第65页/共106页将求得的偏导数代入上式，得将求得的偏导数代入上式，得将上式展开后，除以将上式展开后，除以

31、b bu u，并略去二次项，得，并略去二次项，得第66页/共106页由前方交会公式的推导可知由前方交会公式的推导可知如果再近似取如果再近似取 ，并将上个幻灯，并将上个幻灯片中的片中的(a)(a)式左右两侧同乘以式左右两侧同乘以 则有则有第67页/共106页其中，其中，N1N1为左投影点为左投影点a1a1的投影系数的投影系数第68页/共106页连续法连续法相对定向中相对定向中常数项的几何意义常数项的几何意义X1Y1Z1X2Y2Z2Aa1s1Y1X1Z1a2s2N1v1N2v2Q为定向点上模型上下视差当一个立体像对完成相对定向，Q0当一个立体像对未完成相对定向，即同名光线不相交，Q0第69页/共1

32、06页 上上两式两式便是连续像对相对定向的作业公式。在立体便是连续像对相对定向的作业公式。在立体像对中每量测一对同名像点的像点坐标像对中每量测一对同名像点的像点坐标(x(x1 1,y,y1 1)和和(x(x2 2,y,y2 2)，就可以列出一个，就可以列出一个Q Q方程式。由于式有方程式。由于式有5 5个未知个未知数数 ,因此，至少需要量测因此，至少需要量测5 5对同名像点、对同名像点、当有多余观测值时，将当有多余观测值时，将Q Q视为观测值，由此得到误差视为观测值，由此得到误差方程式；方程式；利用误差方程式，接最小二乖原理组成法方程，解利用误差方程式，接最小二乖原理组成法方程，解求出未知数，

33、即求出未知数，即5 5个相对定向元素的改正数。个相对定向元素的改正数。第70页/共106页由于线性化是取用泰勒展开式的一次项，因此要趋近运算，直到改正数达到所要求的精度为止。最后得到各未知数的值为：第71页/共106页获取已知数据 x0,y0,f ,x1,y1,x2,y2 假定摄影基线 bu=x1-x2设定相对定向元素的初值 2 2 20 由相对定向元素计算像空间辅助坐标 u1,v1,w1,u2,v2,w2逐点计算误差方程式的系数和常数项并法化解法方程，求相对定向元素改正数求相对定向元素的新值判断迭代是否收敛（限差0.01=3100.01=3104 4）相对定向元素计算相对定向元素计算第72页

34、/共106页第73页/共106页2.2.单独像对相对定向元素解求式单独像对相对定向元素解求式共面条件的坐标表达式：U1U2W1V1W2V2第74页/共106页线性化后的误差方程式为：线性化后的误差方程式为：第75页/共106页四、相对定向元素解算过程四、相对定向元素解算过程 摄影测量中，相对定向常用摄影测量中，相对定向常用6 6个标准个标准点位来解求，点位分布如右图，并按点位来解求，点位分布如右图，并按图中位置命名为图中位置命名为1 1，2 2，3 3，4 4，5 5，6 6点、点、其中，其中，l l，2 2点位于像主点点位于像主点O O1 1，O O2 2邻近邻近的明显点，各点距边界的距离应

35、大于的明显点，各点距边界的距离应大于1.5cm1.5cm，而且，而且，1 1，3 3，5 5三点和三点和2 2，4 4，6 6三点尽量位于与三点尽量位于与O O1 1O O2 2连线垂直的直线连线垂直的直线上。上。利用利用6 6对相对定向点的像点坐标，可以组成对相对定向点的像点坐标，可以组成误差方程式，并解之。误差方程式，并解之。第76页/共106页若模型点在像空间辅助坐标系若模型点在像空间辅助坐标系S S1 1-U-U1 1V V1 1W W1 1中的中的坐标为（坐标为（U U，V V，W W），其计算过程为：），其计算过程为：（1 1）根据相对定向元素计算像点的像空间辅）根据相对定向元素计

36、算像点的像空间辅助坐标（助坐标（u u1 1，v v1 1，w w1 1）及（）及（u u2 2，v v2 2，w w2 2）；）；（2 2）计算左右像点的投影系数）计算左右像点的投影系数N N1 1、N N2 2 （3 3）求模型点在像空间辅助坐标系的坐标）求模型点在像空间辅助坐标系的坐标实际计算中，将获得的模型点在像空间辅助实际计算中，将获得的模型点在像空间辅助坐标系的坐标再乘以坐标系的坐标再乘以摄影比例尺的分母摄影比例尺的分母，其，其模型放大成约为实地后，再进行绝对定向。模型放大成约为实地后，再进行绝对定向。五、模型点坐标的计算五、模型点坐标的计算第77页/共106页 相相对对定定向向建

37、建立立的的立立体体模模型型，是是一一个个以以相相对对定定向向中中选选定定的的像像空空间间辅辅助助坐坐标标系系为为基基准准的的模模型型，比比例例尺尺也也是是未未知知的的。要要确确定定立立体体模模型型在在地地面面测测量量坐坐标标系系中中的的正正确确位位置置，则则需需要要把把模模型型点点的的摄摄影影测测量量坐坐标标转转化化为为地地面面测测量量坐坐标标。这这一一工工作作需需要要借借助助于于地地面面测测量量坐坐标标为为已已知知值值的的地地面面控控制制点点来来进进行行，称称为为立立体体模模型型的的绝绝对对定定向向。所所以以，解解析析法法绝绝对对定定向向的的目目的的就就是是将将相相对定向后求出的摄影测量坐标

38、变换为地面测量坐标。对定向后求出的摄影测量坐标变换为地面测量坐标。由由于于地地面面测测量量坐坐标标系系为为左左手手直直角角坐坐标标系系，而而摄摄影影测测量量坐坐标标系系为为右右手手直直角角坐坐标标系系，因因此此，应应首首先先将将地面测量坐标系转换为地面摄影测量坐标系。地面测量坐标系转换为地面摄影测量坐标系。5.5 立体模型的解析法绝对定向立体模型的解析法绝对定向第78页/共106页 我我们们知知道道，一一个个像像对对的的两两张张像像片片有有十十二二个个外外方方位位元元素素，相相对对定定向向求求得得五五个个元元素素后后，要要恢恢复复像像对对的的绝绝对对位位置置，还还要要解解求求七七个个绝绝对对定

39、定向向元元素素：包包括括模模型型的的旋旋转转、平平移移和和缩缩放放。它它需需要要地地面面控控制制点点来来解解求求，这这种种坐坐标标变变换换，在在数数学学上为一个不同原点的三维空问相似变换，其公式为：上为一个不同原点的三维空问相似变换，其公式为：绝对定向元素：绝对定向元素：模型缩放比例因子模型缩放比例因子模型旋转因子模型旋转因子坐标原点平移量坐标原点平移量解析法绝对定向，就是利用已知的地面控制点，从绝对定向解析法绝对定向，就是利用已知的地面控制点，从绝对定向的关系式出发，解求上述七个绝对定向元素。的关系式出发，解求上述七个绝对定向元素。第79页/共106页第80页/共106页第81页/共106页

40、 上上式式中中地地面面摄摄影影测测量量坐坐标标为为已已知知值值，模模型型的的像像空空间间辅辅助助坐坐标标（U、V、W）为为经经相相对对定定向向后后计计算算出出的的值值，由由于于上上式式是是一一个个非非线线性性函函数数，为为便便于于计计算算，首首先先将将其其线线性性化化。为为此此 ，引引入入七七个个绝绝对对定定向向元元素素的的初初始始值及改正数：值及改正数：二、绝对定向公式的线性化及绝对定向元素的解算二、绝对定向公式的线性化及绝对定向元素的解算第82页/共106页代入基本公式，并按泰勒级数展开，取一次小项得代入基本公式，并按泰勒级数展开，取一次小项得F F0 0为用绝对定向元素近似值代入得到的近

41、似值为用绝对定向元素近似值代入得到的近似值，考虑到考虑到小角度的情况，绝对定向基本公式的近似式为小角度的情况，绝对定向基本公式的近似式为第83页/共106页对上式求微分，代入取小值一次项，经整理得：对上式求微分，代入取小值一次项，经整理得：线性化后的公式：线性化后的公式：U UV VW WX XY YZ ZU UV VW WX XSOSOY YSOSOZ ZSOSO上式中有上式中有7 7个未知数，至少需列个未知数，至少需列7 7个方程，则至个方程，则至少需要少需要两个两个平高控制点和平高控制点和一个一个高程控制点，而高程控制点，而且三个控制点不能在同一条直线上。且三个控制点不能在同一条直线上。

42、第84页/共106页将摄影测量坐标视为观测值，列出误差方程式：将摄影测量坐标视为观测值，列出误差方程式：式中式中U UV VW W-W-W0 0U U0 0-W-WV V-V-VU U0 0第85页/共106页法方程的建立与求解法方程的建立与求解量测 2 个平高和 1 个高程以上的控制点可以按最小二乘平差原理求绝对定向元素第86页/共106页 实际绝对定向计算中，为了简便计算，常选模型的重实际绝对定向计算中，为了简便计算，常选模型的重心为坐标系的原点，用心为坐标系的原点，用g g表示，以重心为原点的坐标表示，以重心为原点的坐标称为重心化坐标。称为重心化坐标。重心坐标重心坐标重心化坐标重心化坐标

43、目的目的减减少少模模型型点点坐坐标标在在计计算算过过程程中中的的有有效效位位数数，以以保保证证计计算算的的精度精度使使法法方方程程的的系系数数简简化化，个个别别项项数数值值变变为为零零，以以提提高高计计算算速速度度第87页/共106页重心化坐标代入绝对定向的基本公式重心化坐标代入绝对定向的基本公式第88页/共106页得到用重心化坐标表示的误差方程式：得到用重心化坐标表示的误差方程式：其中，其中，第89页/共106页 求出绝对定向元素以后，将未知点的重求出绝对定向元素以后，将未知点的重心化摄影测量坐标，求出相应的重心化的心化摄影测量坐标，求出相应的重心化的地面摄影测量坐标，然后反求出地面摄影地面

44、摄影测量坐标，然后反求出地面摄影测量坐标，最后将地面摄影测量坐标转回测量坐标，最后将地面摄影测量坐标转回到地面测量坐标，提交成果。到地面测量坐标，提交成果。第90页/共106页获取控制点的两套坐标 Xp,Yp,Zp,Xtp,Ytp,Ztp 给绝对定向参数的初值 1,0,X0,Y0,Z0 绝对定向绝对定向参数的计算参数的计算u计算重心化坐标u计算误差方程式的系数和常数项u解法方程，求绝对定向参数改正数u计算绝对定向参数的新值u判断迭代是否收敛第91页/共106页第92页/共106页四、双像解析的相对定向四、双像解析的相对定向绝对定向法绝对定向法 立体像对相对定向立体像对相对定向-绝对定向法解求模

45、型点的绝对定向法解求模型点的地面坐标，其过程为：地面坐标，其过程为：用连续像对或单独像对的相对定向元素的误用连续像对或单独像对的相对定向元素的误差方程式解求像对的差方程式解求像对的相对定向元相对定向元素；素；由相对定向元素组成左右像片的旋转矩阵，由相对定向元素组成左右像片的旋转矩阵，并利用前方交会式求出模型点在像空间辅助并利用前方交会式求出模型点在像空间辅助坐标系中的坐标；坐标系中的坐标；根据已知地面控制点的坐标，按绝对定向元根据已知地面控制点的坐标，按绝对定向元素的误差方程式解求该立体模型的素的误差方程式解求该立体模型的绝对定向绝对定向元素；元素；按绝对定向公式，将所有待定点的坐标纳入按绝对

46、定向公式，将所有待定点的坐标纳入地面摄影测量坐标中。地面摄影测量坐标中。第93页/共106页 像片控制测量，获得像控点物方坐标像片控制测量，获得像控点物方坐标 像点坐标量测，像点坐标量测，解析相对定向，得（连续法）解析相对定向，得（连续法）解析相对定向+解析绝对定向作业流程计算流程计算流程第94页/共106页计算模型点坐标计算模型点坐标 确定未知数初始值确定未知数初始值第95页/共106页 计计算算 R R计计算算 组组成成误误差差方方程程第96页/共106页 组组成成法法方方程程并并解解第97页/共106页 小于限差计算计算第98页/共106页5.6 5.6 双像解析的光束法严密解双像解析的

47、光束法严密解 光束法：用已知的少数控制点以及待求的地面点，在立体像对内同时解求两光束法：用已知的少数控制点以及待求的地面点，在立体像对内同时解求两像片的外方位元素和待定点的坐标，俗称一步定向法。像片的外方位元素和待定点的坐标，俗称一步定向法。第99页/共106页一、基本公式一、基本公式除了有六个外方位元素外，待定点的地面摄影测量坐标除了有六个外方位元素外，待定点的地面摄影测量坐标X X、Y Y、Z Z也是未知数。也是未知数。将共线方程按泰勒公式展开，取一次小项得：将共线方程按泰勒公式展开，取一次小项得：第100页/共106页对于右片上的像点，可列出右片的误差方程式：对于右片上的像点，可列出右片

48、的误差方程式：对于控制点而言：对于控制点而言：对于待求点而言，式中除了对于待求点而言，式中除了1212个外方位元素个外方位元素外，还有外，还有三个待定三个待定点坐标改正数点坐标改正数。若一对像对中有四个控制点及。若一对像对中有四个控制点及n n个待定点，则可个待定点，则可列出的误差方程个数为：列出的误差方程个数为：16+4n12+3n未知数的个数为：未知数的个数为：(1)(2)将（将（1 1）（）（2 2）式写成矩阵形式为：）式写成矩阵形式为：第101页/共106页第102页/共106页用矩阵形式表示的总的误差方程式为：用矩阵形式表示的总的误差方程式为：法方程式为：法方程式为：用新的符号表示为

49、：用新的符号表示为：第103页/共106页改化法方程式：消去其中的一类未知数所得到的法方程式。改化法方程式：消去其中的一类未知数所得到的法方程式。消去待定点的坐标改正数消去待定点的坐标改正数t t，得改化法方程式为：，得改化法方程式为：若消去外方位元素，则得到另一组改化法方程式：若消去外方位元素，则得到另一组改化法方程式：第104页/共106页三种方法的比较三种方法的比较方法方法严密性严密性控制点要求控制点要求使用场合使用场合后交前交法后交前交法点位精度决定于外方点位精度决定于外方位元素精度位元素精度没有利用多余条件没有利用多余条件3个平高点个平高点已知像片外方位已知像片外方位元素元素相对定向绝对相对定向绝对定向法定向法点位精度决定于相对点位精度决定于相对定向绝对定向精度定向绝对定向精度2个平高点个平高点 1个高程点个高程点航带法加密航带法加密一步定向法一步定向法理论最严密理论最严密精度最高精度最高2个平高点个平高点 1个高程点个高程点光束法加密光束法加密第105页/共106页感谢您的观看！第106页/共106页