《近景摄影测量光束法平差报告.doc》由会员分享，可在线阅读，更多相关《近景摄影测量光束法平差报告.doc（14页珍藏版）》请在金锄头文库上搜索。
1、近景摄影测量光束法平差报告2011 年 6 月 4 日1 作业目的 -32 外业控制点的观测与解算 -33 近景影像获取 -44 LPS刺点点位 -45 光束法平差与精度评定 -56 总结 -111 作业目的以近景摄影测量大实习为基础，对所摄取近景相片解析处理，以外业控制点的解算成果以及内业LPS平差结果为依据，编写光束法平差程序，由22个控制点的像素坐标及5个“已知控制点”的三维坐标求解其余17个控制点的三维坐标，并评定精度。2 作业条件及数据点号像素坐标列（J）像素坐标行（I）XYZ左片：2650.9892114.93497.4532353.7473299.895382792.491225
2、9.531508.8008342.3524298.6832102791.483740.514508.8138342.3548307.0717163928.5592120.49520.2969353.7531300.1146214890.5842130.45527.9857353.5821300.10371648.6242765.5820003660.4521441.4110004728.563816.58500051965.8952557.99600061910.1051210.0700072767.4553044.53100092774.0591493.061000123319.011266
3、5.417000133312.2861986.582000143298.4681284.901000154055.0522705.029000173808.9851539.018000183715.006962.032000193836.444706.426000204883.392691.6510002247541681000234825.4091018.545000右片：2670.9482129.967497.4532353.7473299.895382346.4432264.542508.8008342.3524298.6832102361.448691.079508.8138342.3
4、548307.0717164088.4192115.427520.2969353.7531300.1146205203.4412736.112527.9857353.5821300.10371666.1032764.8820003685.4031472.5740004754.414860.65600051652.4312568.50300061600.0581207.01400072312.0273077.96400092334.4721470.473000123083.1932691.257000133083.0661976.987000143072.4281240.419000154230
5、.5272741.493000173956.4451498.102000183852.033888.02000194075.943613.315000215214.4922119.571000225144.4631507.538000235139.982903.989000外方位元素初始值：(左)Xs1 = 497.9149,Ys1 = 301.2754,Zs1 = 297.2430,Q1=14.9560,W1=4.7765,K1=-0.0308,(右)Xs2 = 509.6501,Ys2 = 301.3448,Zs2 = 297.4727,Q2=2.1777 ,W2=4.5969,K2=0.
6、0791,相片主距：50mm3 平差思想31 共线条件方程本次作业中，光束法平差基于如下共线条件方程：该共线方程是描述摄影中心S、像点a以及物点A位于一直线上的关系式。像点a在成像过程中存在某种系统误差，其改正数（x，y）添加在上式左边，并有：将共线条件方程线性化，其中：32 像点坐标误差方程式由于本作业采用光束法平差，两张照片共44个点（每张22个），每个点可列2个像点坐标误差方程（vx，vy），故总共88个误差方程。而每个误差方程中：有12个外方位元素改正（每张相片6个），6个内方位元素改正（每张相片3个），51个加密点坐标改正（5个控制点坐标改正为0，故只有有17个加密点的三维坐标改正：
7、dX，dY,dZ），4个畸变参数改正（每张相片2个，k1，k2）。33 平差由88个像点坐标误差方程式组成方程组：,其中：V阵为像点坐标误差改正，88行*1列；A阵为系数阵，88行*73列；X阵伟未知数阵，73行*1列（未知数包括12个外方位元素、6个内方位元素、51个加密点三维坐标、4个畸变参数）L阵为常数项阵，88行*1列由式X=（AtA）-1 ATL解求未知数即可4 程序#include #include #include #include#includeconst int N=90;/求转置矩阵templatevoid Transpose(T1*mat1,T2*mat2,int a,i
8、nt b)int i,j;for(i=0;ia;i+)for(j=0;jb;j+)mat2ji=mat1ij;return;/求矩阵的乘积templatevoid Array_mul(T1*mat1,T2 * mat2,T2 * result,int a,int b,int c) int i,j,k;for(i=0;ia;i+)for(j=0;jc;j+)resultij=0;for(k=0;kb;k+)resultij+=mat1ik*mat2kj;return;/求逆矩阵inverse(double ANN,int m)int i=0,j=0,k=0;double C2020,B2020;
9、for(i=0;i2*m;i+) for(j=0;j2*m;j+) if(i=j) Cij=1.0; else Cij=0.0;for(i=0;im;i+) for(j=0;jm;j+)Bij=Aij;for(i=0;im;i+) for(j=m;j2*m;j+)Bij=Cij-m;cout.precision(5);for(k=0;km;k+)for(i=k;i=0;j-)if(Bik!=0)Bij=Bij/Bik; for(i=m-1;ik;i-)if(Bik=0)continue; for(j=0;j2*m;j+) Bij=Bij-Bkj; for(k=1;km;k+) for(i=0;
10、i=i;j-) Bij=Bij-Bik*Bkj;for(i=0;im;i+) for(j=0;jm;j+)Aji=Bjm+i; return 1; void main()/定义两张相片共个控制点和加密点的像素坐标（J1,I1,J2,I2）和像平面直角坐标（x1,y1,x2,y2），及地面坐标（X,Y,Z）double J1N=0.0,I1N=0.0,J2N=0.0,I2N=0.0;double x1N=0.0,x2N=0.0,z1N=0.0,z2N=0.0;double XN=0.0,YN=0.0,ZN=0.0; /导入控制点坐标数据ifstream infile;infile.open(控制点坐标数据.txt); if(infile.is_open() while(!infile.eof ()for(int i=0;iJ1i;infile.ignore(1); infileI1i;infile.ignore(1); infileXi;infile.ignore(1); infileYi;infile.ignore(1);
