1、5852.099 527.925 GCP2 88.56 81.134 -5.252 78.184 5780.02 5906.365 571.549 GCP3 13.362 -79.37 -79.122 -78.879 5210.879 4258.446 461.81 GCP4 82.24 -80.027 -9.887 -80.089 5909.264 4314.283 455.484 1 51.758 81.555 -39.953 78.463 2 14.618 -0.231 -76.016 0.036 3 49.88 -0.792 -42.201 -1.022 4 86.243 -1.346
2、 -7.706 -2.112 5 48.135 -79.962 -44.438 -79.736 内方位元素:f=152.000mm,x0=0,y0=0四 实验框图输入GCP的像点坐标 xy确定初始值=0,Xs,Ys,Zs计算旋转矩阵R计算像点在像空间坐标系中的近似值(x),(y),并组成误差方程的常数项l计算误差方程的系数项组成系数矩阵A组成法方程式,计算系数AA常数项AL求解外方位元素计算、Xs、Ys、Zs改正后的值改正数是否小于限差? 大于小于计算完毕此过程完成空间后方交会求解像片的外方位元素,其中改正数小于限差(0.00003,相当于0.1的角度值)为止。在这个过程中采用迭代的方法,是外
3、方位元素逐渐收敛于理论值,每次迭代所得的改正数都应加到上一次的初始值之中。输入所需计算点的像平面坐标x1,y1;x2,y2根据后方交会所得的旋转矩阵R1,R2计算像点在左右像空间辅助坐标系中的坐标X1 Y1 Z1,X2 Y2 Z2 计算摄影基线的三个坐标分量Bx By Bz 计算个点在左右像片中的的投影系数N1 N2计算地面所求点在地面摄影测量坐标系中的坐标XA YA ZA计算完毕,精度评定在空间后方交会中运用的数学模型为共线方程确定Xs,Ys,Zs的初始值时,对于左片可取地面左边两个GCP的坐标的平均值作为左片Xs 和Ys的初始值,取右边两个GCP的坐标平均值作为右片Xs 和Ys的初始值。Z
4、s可取地面所有GCP的Z坐标的平均值再加上航高。空间前方交会的数学模型为: 五 实验源代码function Main_KJQHFJH()global R g1 g2 m G a c b1 b2;m=10000;a=5;c=4;feval(shuru); %调用shuru()shurujcp()函数完成像点及feval(shurujcp); %CCP有关数据的输入 XYZ=feval(MQZqianfangjh); %调用MQZqianfangjh()函数完成空间前方、% 单位权中误差 % %后方交会计算解得外方位元素global V1 V2; %由于以上三个函数定义在外部文件中故需VV=; %
5、用feval()完成调用过程for i=1:2*c VV(i)=V1(i);VV(2*i+1)=V2(i);endm0=sqrt(VV*(VV)/(2*c-6);disp(单位权中误差m0为正负:);disp(m0); %计算单位权中误差并将其输出显示输入GCP像点坐标及地面摄影测量坐标系坐标的函数和输入所求点像点坐标函数:function shurujcp()global c m;m=input(摄影比例尺: %输入GCP像点坐标数据函数并分别将其c=input(GCP的总数= % 存入到不同的矩阵之中GCP左片像框标坐标:global g1;g1=zeros(c,2);i=1;while
6、i=c m=input(x= n=input(y= g1(i,1)=m;g1(i,2)=n; i=i+1;end GCP右片像框标坐标:global g2;g2=zeros(c,2); g2(i,1)=m;g2(i,2)=n; %function shuru() global a;a=input(计算总像对点数= %完成想计算所需的像平面坐标global b1; %坐标输入,存入不同的矩阵中b1=zeros(a,2);左片像点坐标:)=a b1(i,1)=m;b1(i,2)=n;end %global b2;b2=zeros(a,2);右片像点坐标: b2(i,1)=m;b2(i,2)=n;g
7、lobal c;GCP摄影测量系坐标:global G;G=zeros(3,c);X=Y= v=input(Z= G(i,1)=m;G(i,2)=n;G(i,3)=v;%空间前方交会和后方交会函数:function XYZ=MQZqianfangjh()global R1 R2 a f b1 b2 Ra Rb;global X1 X2;R1=Ra;R2=Rb;R1=zeros(3,3);R2=zeros(3,3);global g1 g2 G V1 V2 V WF c QXX QXX1 QXX2;xs0=(G(1,1)+G(3,1)/2;ys0=(G(1,2)+G(3,2)/2;Xs1,Ys1
8、,Zs1,q1,w1,k1 R=houfangjh(g1,xs0,ys0); %对左片调用后方交会函数R1=R;V1=V;WF1=WF;QXX1=QXX;save 左片外方位元素为.txt WF -ascii %将计算所得的外方位元素存入到.txt % 文件中c g1(i,1)=g1(i,1)+V1(2*i-1);g1(i,2)=g1(i,2)+V1(2*i);左片像点坐标.txt g1 -ascii xs0=(G(2,1)+G(4,1)/2;ys0=(G(2,2)+G(4,2)/2;Xs2,Ys2,Zs2,q2,w2,k2 R=houfangjh(g2,xs0,ys0); %对右片调用后方交
9、会函数 R2=R;V2=V;WF2=WF;QXX2=QXX;右片外方位元素为.txt WF ascii %将计算所得的外方位元素存入到.txt % 文件中c g2(i,1)=g2(i,1)+V2(2*i-1);g2(i,2)=g2(i,2)+V2(2*i);右片像点坐标.txt g2 -ascii X1=zeros(a,3);X2=zeros(a,3);xx=zeros(3,1);xxx=zeros(3,1);a ss=b1(i,1);b1(i,2);-f;dd=b2(i,1);b2(i,2); xx=R1*ss; X1(i,:)=xx;xxx=R2*dd;X2(i,:)=xxxglobal
10、Xs1 Xs2 Ys1 Ys2 Zs1 Zs2; BX=Xs2-Xs1;BY=Ys2-Ys1;BZ=Zs2-Zs1;global N1 N2;N1=zeros(1,a);N2=zeros(1,a); N1(1,i)=(BX*X2(i,3)-BZ*X2(i,1)/(X1(i,1)*X2(i,3)-X2(i,1)*X1(i,3); N2(1,i)=(BX*X1(i,3)-BZ*X1(i,1)/(X1(i,1)*X2(i,3)-X2(i,1)*X1(i,3);end %计算投影系数,并计算五点的三维坐标global XYZ;XYZ=zeros(a,3); XYZ(i,1)=Xs1+N1(1,i)*X1(i,1); XYZ(i,3
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1