摄影测量实验报告（空间后方交会—前方交会）.pdf

摄摄影影测测量量实实验验报报告告（（空空间间后后⽅⽅交交会会—前前⽅⽅交交会会））

空间后⽅交会-空间前⽅交会程序编程实验⼀．实验⽬的要求

掌握运⽤空间后⽅交会-空间前⽅交会求解⾯点的空间位置。学会运⽤空间后⽅交会的原理，根据所给控制点的⾯摄影测量

坐标系坐标以及相应的像平⾯坐标系中的坐标，利⽤计算机编程语⾔实现空间后⽅交会的过程，完成所给像对中两张像⽚各⾃

的外⽅位元素的求解。然后根据空间后⽅交会所得的两张像⽚的内外⽅位元素，利⽤同名像点在左右像⽚上的坐标，求解其对

应的⾯点在摄影测量坐标系中的坐标，并完成精度评定过程，利⽤计算机编程语⾔实现此过程。

⼆．仪器⽤具

计算机、编程软件（MATLAB）

三．实验数据

实验数据包含四个⾯控制点（GCP）的⾯摄影测量坐标及在左右像⽚中的像平⾯坐标。此四对坐标运⽤最⼩⼆乘法求解左

右像⽚的外⽅位元素，即完成了空间后⽅的过程。另外还给出了5对⾯点在左右像⽚中的像平⾯坐标和左右像⽚的内⽅位元

素。实验数据如下：

内⽅位元素：f=152.000mm，x0=0，y0=0

四．实验框图

此过程完成空间后⽅交会求解像⽚的外⽅位元素，其中改正数⼩于限差（0.00003，相当于0.1’的⾓度值）为⽌。在这个过程中

采⽤迭代的⽅法，是外⽅位元素逐渐收敛于理论值，每次迭代所得的改正数都应加到上⼀次的初始值之中。

在空间后⽅交会中运⽤的数学模型为共线⽅程

确定Xs，Ys，Zs的初始值时，对于左⽚可取⾯左边两个GCP的坐标的平均值作为左⽚Xs和Ys的初始值，取右边两个GCP

的坐标平均值作为右⽚Xs和Ys的初始值。Zs可取⾯所有GCP的Z坐标的平均值再加上航⾼。

空间前⽅交会的数学模型为：

五．实验源代码

functinMain_KJQHFJH()

glbalRg1g2mGacb1b2;

m=10000;a=5;c=4;

feval(@shuru);%调⽤shuru（）shurujcp（）函数完成像点及feval(@shurujcp);%CCP有关数据的输⼊

XYZ=feval(@MQZqianfangjh);%调⽤MQZqianfangjh（）函数完成空间前⽅、%%%%%%单位权中误差%%%%%后⽅交会

计算解得外⽅位元素

glbalV1V2;%由于以上三个函数定义在外部⽂件中故需VV=[];%⽤feval（）完成调⽤过程

fri=1:2*c

VV(i)=V1(i);VV(2*i+1)=V2(i);

end

m0=sqrt(VV*(VV)/(2*c-6));

disp(单位权中误差m0为正负:);disp(m0);%计算单位权中误差并将其输出显⽰

输⼊GCP像点坐标及⾯摄影测量坐标系坐标的函数和输⼊所求点像点坐标函数：

functinshurujcp()

glbalcm;

m=input(摄影⽐例尺：);%输⼊GCP像点坐标数据函数并分别将其c=input(GCP的总数=);%存⼊到不同的矩阵之中

disp(GCP左⽚像框标坐标：);

glbalg1;g1=zers(c,2);

i=1;

whilei=c

m=input(x=);

n=input(y=);

g1(i,1)=m;g1(i,2)=n;

i=i+1;

end

disp(GCP右⽚像框标坐标：);

glbalg2;g2=zers(c,2);

i=1;

whilei=c

m=input(x=);

n=input(y=);

g2(i,1)=m;g2(i,2)=n;

i=i+1;

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functinshuru()

glbala;

a=input(计算总像对点数=);%完成想计算所需的像平⾯坐标glbalb1;%坐标输⼊，存⼊不同的矩阵中b1=zers(a,2);

disp(左⽚像点坐标：)

i=1;

whilei=a

m=input(x=);

n=input(y=);

b1(i,1)=m;b1(i,2)=n;

i=i+1;

end

%%

glbalb2;

b2=zers(a,2);

disp(右⽚像点坐标：)

i=1;

whilei=a

m=input(x=)