《摄影测量学》
课程设计报告书
1
目录
一、作业任务 ............................................. 3
二、计算原理 ............................................. 4
2.1 点投影系数法计算原理 .............................. 4
2.2 共线方程严密法计算原理 ............................ 6
三、编程实现 ............................................. 7
四、计算结果 ............................................. 7
4.1 点投影系数法计算结果 .............................. 7
4.2 共线方程严密法计算结果 ............................ 8
五、学习体会 ............................................. 8
六、附录 ................................................. 9
6.1 点投影系数法源代码及运行截图 ...................... 9
6.2 共线方程严密法源代码及运行截图 ................... 11
2
一、作业任务
已知一个立体像对:
左 像 片 的 外 方 位 元 素 为 : (14922.9778, 11532.4077, 3230.3383)m, (-0.0566,
-0.1215, -0.6456)°,
右像片的外方位元素为: (16296.9601, 11524.2256, 3239.6070)m, (0.0429, -0.5148,
-0.1538)°,
相机的内方位元素(x0, y0, f)为(-0.004,-0.008, 152.72)mm。
现量测了 6 个像点的框标坐标,如下(单位 mm):
左片: 1(85.7255, 69.6561), 2(0.3206, 59.1841), 3(83.7954, -1.9617),
4(-2.2603, -15.1519), 5(88.2371,-73.7021), 6(2.4792,-67.5320)
右片: 1 ( 0.6313, 70.1653), 2 ( -84.4506, 60.4553), 3( -3.3174, -1.3205)
4 ( -97.0211, -13.7425), 5(2.7516, -72.5590), 6( -82.9161, -65.7818)
采用点投影系数法和共线方程严密法分别计算上述 6 个点的地面坐标。
3
二、计算原理
2.1 点投影系数法计算原理
1、如下图所示,航摄机在两相邻摄站 和 分别拍摄了一张像片,构成立体像对。地
面上任意一点 A 在左、右像片上的构像分别为 和 。为了确定像点和地面点的数学关
系,建立地面摄影测量坐标系
, 轴与航向基本一致。过左测站 建立
与地面摄影测量坐标系平行的像空间辅助坐标系
, 再过右摄站 也建立与
平行的像空间辅助坐标系
。
2、设地面点 在地面摄影测量坐标系
中的坐标为
,对应像点
、 在各自的像空间坐标系中的坐标为
和
,在像空间辅助坐
标系中的坐标分别为
和
。现已知两张像片的外方位元素,就可以
4
1S2S1a2atPtPtPDXYZtPZ1S1111SXYZ2StPtPtPDXYZ2222SXYZAtPtPtPDXYZ,,AAAXYZ1a2a11,,xyf22,,xyf111,,XYZ222,,XYZ
由像点的像空间坐标计算出该点的像空间辅助坐标,即
上式中, 、 为由已知的外方位角元素算得的左、右像片的旋转矩阵。
3、右摄站 在
中的坐标,即摄影基线 的三个分量 、 、 可由外
方位直线元素算得
4、因为左、右像空间辅助坐标系与地面摄影测量坐标系相互平行,且摄站点、像点、
地面点、三点共线。则由上图可得
式中: 、 为左、右像点的投影系数。一般情况下,不同的点有不同的投影系数
值。
5、根据上式可以得到前方交会法计算地面点坐标的公式,即
式
中, 、 仍然未知,为此,结合式子
,有
6、由式子
中的一、三两式联立求解投影系数得
5
112211121212,XxXxYRyYRyZfZf(211)1R2R2S1111SXYZBXBYBZB212121XSSYSSZSSBXXBYYBZZ(212)11122211111112222222ASASASASASASXXYYZZSANSaXYZXXYYZZSANSaXYZ(213)1N2N121212112211221122SSAASSASSXXXNXNXYYNYYNYZNZNZZZ(214)(214)1N2N(212)112211221122XYZBNXNXBNYNYBNZNZ(215)(215)
2.2 共线方程严密法计算原理
1、一个地面点 分别对应左右两张像片两个点 、 ,每个像点对应两个共线方程,
即一个地面未知点对应四个共线方程 ,即
上式中, 、 、 分别表示旋转矩阵中的 中的 、 、 元素。
2、将四个共线方程分别对 X、Y、Z 进行线性化得到关于 X、Y、Z 三个未
知数的四个误差线性方程,即
上式中各系数和常数项分别为
6
22112211121221XZXZBZBXNXZXZBZBXNXZXZ(216)A1a2a111111111111222111111101313131212121013131321212120()()()()()()()()()()()()()()()ASASASASASASASASASASASASASASASaXXbYYcZZxxfaXXbYYcZZaXXbYYcZZyyfaXXbYYcZZaXXbYYcZZxxfa22222222223232322222220232323()()()()()()()()()ASASASASASASASASASXXbYYcZZaXXbYYcZZyyfaXXbYYcZZ(221)ijaijbijciRjajbjc11221112131415162122232425260000xyxylXlYlZllXlYlZllXlYlZllXlYlZl(222)
其中, =1,2。
3、用矩阵形式表示误差方程式为:
其中
4、根据最小二乘原理,可以解的
三、编程实现
采用 C++语言进行了点投影系数法和共线方程严密法求地面点坐标的程序编写,源
代码和运行结果截图详见附录。
四、计算结果
4.1 点投影系数法计算结果
7
11032103310342035203620311103030322()()()()()()()()()iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiixiSiSiSiiSiiSiiSyiSilfaxxalfbxxblfcxxclfayyalfbyyblfcyyclfaXfbYfcZxxaXxxbYxxcZlfaXfbY2030303()()()iiiiiSiSiiSiiSiiSfcZyyaXyybYyycZiVAXL(223)1122111213141516212223242526,TxyxyllllllALllllllllll1()()TTXAAAL(224)
X
16312.0739
14936.5546
16246.0630
14885.6112
16340.7515
14948.3586
Y
12632.1585
12483.0367
11481.5677
11308.3244
10313.7418
10436.3904
点 号
点 1
点 2
点 3
点 4
点 5
点 6
4.2 共线方程严密法计算结果
X
16312.0739
14936.5543
16246.0622
14885.6096
16340.7505
14948.3600
Y
12632.1585
12483.0373
11481.5677
11308.3238
10313.7432
10436.3929
点 号
点 1
点 2
点 3
点 4
点 5
点 6
Z
770.0758
761.4998
812.7402
1016.9372
750.0083
766.0779
Z
770.0758
761.4975
812.7429
1016.9302
750.0117
766.0846
五、学习体会
通过本次实验,我对双像片的后方交会—前方交会方法求地面点坐标的方法有了更
深刻的认识,在方法和原理上都有了进一步的理解,同时也锻炼了编程能力。本次编程
主要延续了第三次作业的方法,虽然完成了作业任务,但还有一些问题,比如没能实现
对文件的读取以及窗口的设计,算法的设计也不够简便,而且在算法设计过程也出现了
一些错误,花了不少时间才改正过来,感觉的调试程序的能力还不够。所以我感觉编程
能力还有待提高,这是以后要努力的方向。
8