文档介绍：摄影测量实验报告(空间后方交会-前方交会)空间后方交会--空间前方交会求解地面点的空间位置。学会运用空间后方交会的原理,根据所给控制点的地面摄影测量坐标系坐标以及相应的像平面坐标系中的坐标,利用计算机编程语言实现空间后方交会的过程,完成所给像对中两张像片各自的外方位元素的求解。然后根据空间后方交会所得的两张像片的内外方位元素,利用同名像点在左右像片上的坐标,求解其对应的地面点在摄影测量坐标系中的坐标,并完成精度评定过程,利用计算机编程语言实现此过程。、编程软件(MATLAB)(GCP)的地面摄影测量坐标及在左右像片中的像平面坐标。此四对坐标运用最小二乘法求解左右像片的外方位元素,即完成了空间后方的过程。另外还给出了5对地面点在左右像片中的像平面坐标和左右像片的内方位元素。实验数据如下:内方位元素:f=,x0=0,y0=,其中改正数小于限差(,’的角度值)为止。在这个过程中采用迭代的方法,是外方位元素逐渐收敛于理论值,每次迭代所得的改正数都应加到上一次的初始值之中。在空间后方交会中运用的数学模型为共线方程确定Xs,Ys,Zs的初始值时,对于左片可取地面左边两个GCP的坐标的平均值作为左片Xs和Ys的初始值,取右边两个GCP的坐标平均值作为右片Xs和Ys的初始值。Zs可取地面所有GCP的Z坐标的平均值再加上航高。空间前方交会的数学模型为:()globalRg1g2mGacb1b2;m=10000;a=5;c=4;feval(***@shuru);%调用shuru()shurujcp()函数完成像点及feval(***@shurujcp);%CCP有关数据的输入XYZ=feval(***@MQZqianfangjh);%调用MQZqianfangjh()函数完成空间前方、%%%%%%单位权中误差%%%%%后方交会计算解得外方位元素globalV1V2;%由于以上三个函数定义在外部文件中故需VV=[];%用feval()完成调用过程fori=1:2*cVV(i)=V1(i);VV(2*i+1)=V2(i);endm0=sqrt(VV*(VV’)/(2*c-6));disp(‘单位权中误差m0为正负:’);disp(m0);%计算单位权中误差并将其输出显示输入GCP像点坐标及地面摄影测量坐标系坐标的函数和输入所求点像点坐标函数:functionshurujcp()globalcm;m=input(‘摄影比例尺:’);%输入GCP像点坐标数据函数并分别将其c=input(‘GCP的总数=‘);%存入到不同的矩阵之中disp(‘GCP左片像框标坐标:’);globalg1;g1=zeros(c,2);i=1;whileim=input(‘x=‘);n=input(‘y=‘);g1(i,1)=m;g1(i,2)=n;i=i+1;enddisp(‘GCP右片像框标坐标:’);globalg2;g2=zeros(c,2);i=1;whileim=input(‘x=‘);n=input(‘y=‘);g2(i,1)=m;g2(i,2)=n;i=i+1;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functionshuru()globala;a=input(‘计算总像对点数=‘);%完成想计算所需的像平面坐标globalb1;%坐标输入,存入不同的矩阵中b1=zeros(a,2);disp(‘左片像点坐标:’)i=1;whileim=input(‘x=‘);n=input(‘y=‘);b1(i,1)=m;b1(i,2)=n;i=i+1;end%%globalb2;b2=zeros(a,2);disp(‘右片像点坐标:’)i=1;whileim=input(‘x=‘);n=input(‘y=‘);b2(i,1)=m;b2(i,2)=n;i=i+1;end%%globalc;c=input(‘GCP的总数=‘);disp(‘GCP摄影测量系坐标:’)globalG;G=zeros(3,c);i=1;whileim=input(‘X=‘);n=input(‘Y=‘);v=input(‘Z=‘);G(i,1)=m;G(i,2)=n;G(i,3)=v;i=i+1;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%空间前方交会和后方交会函数:functionXYZ=MQZqianfangjh()globalR1R2afb1b2RaRb;glob