欢迎来到天天文库
浏览记录
ID:21075141
大小:140.59 KB
页数:22页
时间:2018-10-19
《近景摄影测量光束法平差报告》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、近景摄影测量光束法平差报告近景摄影测量光朿法平差报告2011年6月4日1作业冃的2夕卜射空制点的观测与解算3近景影像获取4LPS刺点点位5光束法平差与精度评定6总结33445111作业目的以近景摄影测量大实习为基础,对所摄取近景相片解析处理,以外业控制点的解算成果以及内业LPS平差结果为依据,编写光束法平差程序,由22个控制点的像素坐标及5个“己知控制点”的三维坐标求解其余17个控制点的三维坐标,并评定精度。2作业条件及数据点号像素坐标列0)像素坐标行(I)XYZ左片:2650.98982792.491102791.483163928.559214890.5841
2、648.6243660.4524728.56351965.89561910.10572767.45592774.059193836.444204883.39224754234825.409右片:2670.94882346.443102361.448164088.419205203.4411666.1033685.4034754.41451652.4312114.932259.531740.5142120.492130.452765.5821441.411816.5852557.9961210.073044.5311493.0612665.4171986.58212
3、84.9012705.0291539.018962.032706.4262691.65116811018.5452129.9672264.542691.0792115.4272736.1122764.8821472.574860.6562568.503497.4532508.8008508.8138520.2969527.985700000000000000000497.4532508.8008508.8138520.2969527.98570000353.7473342.3524342.3548353.7531353.5821000000000000000003
4、53.7473342.3524342.3548353.7531353.58210000299.8953298.6832307.0717300.1146300.103700000000000000000299.8953298.6832307.0717300.1146300.10370000679121314151718192122231600.0582312.0272334.4723083.1933083.0663072.4284230.5273956.4453852.0334075.9435214.4925144.4635139.9821207.0143077.9
5、641470.4732691.2571976.9871240.4192741.4931498.102888.02外方位元素初始值:(左)Xsl=497.9149,Ysl=301.2754,Zsl=297.2430,Ql=14.9560,Wl=4.7765,Kl=-0.0308,(右)Xs2=509.6501,Ys2=301.3448,Zs2=297.4727^2=2.1777,W2=4.5969,K2=0.0791,相片主距:50mm3平差思想3.1共线条件方程木次作业中,光束法平差基于如下共线条件方程:该共线方程是描述摄影中心S、像点a以及物点A位于一直线上的关
6、系式。像点a在成像过程中存在某种系统误差,其改正数(Ax,Ay)添加在上式左边,并有:将共线条件方程线性化,其中:3.2像点坐标误差方程式由于本作业采用光朿法平差,两张照片共44个点(每张22个),每个点可列2个像点坐标误差方程(vx,vy),故总共88个误差方程。而每个误差方程中:宥12个外方位元素改正(每张相片6个),6个内方位元素改正(每张相片3个),51个加密点坐标改正(5个控制点坐标改正为0,故只有有17个加密点的三维坐标改正:dX,dY,dZ),4个畸变参数改正(每张相片2个,kl,k2)。3.3平差由88个像点坐标误差方程式组成方程组:V?AX?L,
7、其屮:A阵为系数阵,88行*73列;X阵伟未知数阵,73行*1列(未知数包括12个外方位元素、6个内方位元素、51个加密点三维坐标、4个畸变参数)L阵为常数项阵,88行*1列由式X=(AtA)-1ATL解求未知数即可4程序#include<iomanip.h>include<stdlib.h>include<math.h>#include<fstream.h>#include<iostream.h>constintN=90;//求转置矩阵template<typenameTl,typenameT2&
8、gt;vo
此文档下载收益归作者所有