《单片空间后方交会报告》由会员分享，可在线阅读，更多相关《单片空间后方交会报告（12页珍藏版）》请在金锄头文库上搜索。
1、为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划单片空间后方交会报告单片空间后方交会程序设计实习报告班级：测绘一班姓名：刘彦德学号：实习任务：用C或VC+语言实现单片后方交会的计算。并验算所给数据得出所给像片的的外方位元素。实习目的：通过实习编程计算深入理解单片空间后方交会的原理和计算步骤；通过上机调试程序加强动手能力的培养；实习原理：单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应的像片坐标量测值出发，根据共线条件方程，求解该影像在航空摄影时刻的像片外方位元素Xs，Ys，Zs，。共线条件方程如下：x-x0=
2、-f*a1(X-Xs)+b1(Y-Ys)+c1(Z-Zs)/a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)y-y0=-f*a2(X-Xs)+b2(Y-Ys)+c2(Z-Zs)/a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)x,y为像点的像平面坐标；x0,y0,f为影像的外方位元素；Xs,Ys,Zs为摄站点的物方空间坐标；X,Y,Z为物方点的物方空间坐标；程序框图：实习数据：像片内方位元素：f=y0=0；摄影比例尺：1：2500个人验证数据：第二张像片的53个点；试验要求：从中任选1张像片进行空间后方交会。最后解算出外方位元素，并评定各元素精度，和求解单位权中误差。要求：线元素精度m
3、x等高于米；角元素精度高于弧度。同时和提供的外方位元素进行比较，要求线元素误差小于米；角元素误差30秒。输出结果：实习总结：关于此次的上机实习，虽然在实习的初期对单片空间后方交会的理论知识理解的不是很好，但在实习过程中一步步摸索前进，最后对后方交会的计算过程整体的来龙去脉有了一个大致的认识。在此次实习中也碰到过不少的关于部分计算程序实现的难题，但在老师的指导和同学的帮助下还是完成在整个实习，在此表示真心的感谢。现在我更进一步的了解了单片空间后方交会的原理与计算步骤，加深了对课堂上所听到的理论知识的理解与认识，理论联系实践，加深了对其的认识。同时还增强了自己动手编程和程序调试的能力，使自己在程序
4、设计和编辑上更加熟悉。空间后方交会编程实习报告一实习目的用程序设计语言编写一个完整的单片空间后方交会程序，通过对提供的试验数据进行计算，输出像片的外方位元素并评定精度。本实验的目的在于让学生深入理解单片空间后方交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过程。通过上机调试程序加强动手能力的培养，通过对实验结果的分析，增强学生综合运用所学知识解决实际问题的能力。二实习内容利用一定数量的地面控制点，根据共线条件方程求解像片外方位元素。三实习数据已知航摄仪的内方位元素：fk=，x0=y0=，摄影比例尺为1:50000；4个地面控制点的地面坐标及其对应像点的像片坐
5、标：四实习原理如果我们知道每幅影像的6个外方位元素，就能确定被摄物体与航摄影像的关系。因此，如何获取影像的外方位元素，一直是摄影测量工作者所探讨的问题。可采取的方法有：利用雷达、全球定位系统、惯性导航系统以及星相摄影机来获取影像的外方位元素；也可以利用影像覆盖范围内一定数量的控制点的空间坐标与摄影坐标，根据共线条件方程，反求该影像的外方位元素，这种方法称为单幅影像的空间后方交会。单像空间后方交会的基本思想是：以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，解求该影像在航空摄影时刻的外方位元素Xs，Ys，Zs，t，w，k。五实习流程
6、获取已知数据。从摄影资料中查取影像比例尺1/m，平均摄影距离量测控制点的像点坐标并进行必要的影像坐标系统误差改正，得到像点坐标。确定未知数的初始值。单像空间后方交会必须给出待定参数的初始值，在竖直航空摄影且地面控制点大体对称分布的情况下，可按如下方法确定初始值：Zs=H=m*f+Zi/4；Xs=xi/n;Y0s=Yi/n;t=0;式中：m为摄影比例尺分母；(4)计算旋转矩阵R。利用角元素的近似值按下式计算方向余弦值a1,a2,a3,b1,b2,b3,c1,c2,c3，组成R阵。逐点计算像点坐标的近似值。利用未知数的近似值按共线条件方程计算控制点像点坐标的近似值、。按下式逐点计算误差方程式的系数
7、和常数项，组成误差方程。(7)计算法方程的系数矩阵AA与常数项AL，组成法方程；解求外方位元素。根据法方程，解求外方位元素的改正数，并与相应的近似值求和，得到外方位元素新的近似值。检查计算是否收敛。将所求得的外方位元素的改正数与规定的限差比较，通常对t、Xs、Ys、Zs的改正数t，Xs，Ys，Zs给予限差，当改正数小于限差时，迭代结束。否则用新的近似值重复步骤计算，直到满足要求为止。空间后方交会的精度估计：按上述方法所求得的影像外方位元素的精度可以通过法方程式中未知数的系数矩阵的逆阵-1来解求，此时视像点坐标为等精度不相关观测值。因为ATA）-1中第i个主对角线上的元素Qii就是法方程式中第i
8、个未知数的权倒数，若单位权中误差为m0，则第i个未知数的中误差为：mi=Qiim0当参加空间后方交会的控制点有n个时，则单位权中误差可按下式计算：m?vv2n?6TT六主要代码与详解voidR(doublet,doublew,doublek,double*a,double*b,double*c)a0=cos(t)*cos(k)-sin(t)*sin(w)*sin(k);a1=-cos(t)*sin(k)-sin(t)*sin(w)*cos(k);a2=-sin(t)*cos(w);b0=cos(w)*sin(k);b1=cos(w)*cos(k);b2=-sin(w);c0=sin(t)*co
9、s(k)+cos(t)*sin(w)*sin(k);c1=-sin(t)*sin(k)+cos(t)*sin(w)*cos(k);c2=cos(t)*cos(w);/子函数计算旋转矩阵R。利用角元素的近似值按下式计算方向余弦值a1,a2,a3,b1,b2,b3,c1,c2,c3，组成R阵。voidmain()inti,m,num;doublet,w,k,Xs,Ys,Zs,f;/六个外方位元素与焦距fdoublexN=-,-,-,yN=-,-,;/4个像点坐标doubleXN=,YN=,ZN=,;/四个控制点的空间坐标m=50000;f=;/影像比例尺1/m与焦距f/以上主要是已知数据的定义do
10、ubleV6=0,a3,b3,c3,XoN,YoN,ZoN,A6*M,B6*M,lM,C36,D6,E8;for(i=0;iN;i+)xi=xi/;yi=yi/;/像点坐标单位换算/以下为确定未知数的初始值。Xs=(X0+X1+X2+X3)/N;Ys=(Y0+Y1+Y2+Y3)/N;Zs=m*f+(Z0+Z1+Z2+Z3)/N;/Xs，Ys，Zs为摄站点的空间坐标初始值f=f/;t=w=k=;/角元素的近似值for(num=1;num10;num+)/num控制循环迭代次数for(i=0;iN;i+)/下面的循环用来逐点计算误差方程式的系数和常数项，组成误差方程。/以下是用共线条件方程计算控制
11、点像点坐标的近似值、。Xoi=-f*(a0*(Xi-Xs)+b0*(Yi-Ys)+c0*(Zi-Zs)/(a2*(Xi-Xs)+b2*(Yi-Ys)+c2*(Zi-Zs);Yoi=-f*(a1*(Xi-Xs)+b1*(Yi-Ys)+c1*(Zi-Zs)/(a2*(Xi-Xs)+b2*(Yi-Ys)+c2*(Zi-Zs);Zoi=a2*(Xi-Xs)+b2*(Yi-Ys)+c2*(Zi-Zs);/Zo便于后面计算A12*i+0=(a0*f+a2*xi)/Zoi;A12*i+1=(b0*f+b2*xi)/Zoi;A12*i+2=(c0*f+c2*xi)/Zoi;R(t,w,k,a,b,c);/计算
12、旋转矩阵R。利用角元素的近似值按下式计算方向余弦值a1,a2,a3,b1,b2,b3,c1,c2,c3，组成R阵。A12*i+3=yi*sin(w)-(xi*(xi*cos(k)-yi*sin(k)/f+f*cos(k)*cos(w);A12*i+4=-f*sin(k)-xi*(xi*sin(k)+yi*cos(k)/f;A12*i+5=yi;A12*i+6=(a1*f+a2*yi)/Zoi;A12*i+7=(b1*f+b2*yi)/Zoi;A12*i+8=(c1*f+c2*yi)/Zoi;A12*i+9=-xi*sin(w)-(yi*(xi*cos(k)-yi*sin(k)/f-f*sin(
13、k)*cos(w);A12*i+10=-f*cos(k)-yi*(xi*sin(k)+yi*cos(k)/f;A12*i+11=-xi;l2*i=xi-Xoi;l2*i+1=yi-Yoi;/计算lzhuanzhi(A,B,8,6);/转置Axiangchen(B,A,C,6,8,6);/求ATAxiangchen(B,l,D,6,8,1);/求常数项ATLqiuni(C,6);/求ATA的逆xiangchen(C,D,V,6,6,1);/求改正数VXs+=V0;Ys+=V1;Zs+=V2;t+=V3;w+=V4;k+=V5;/结果改正if(fabs(V0)=)&(fabs(V1)=)&(fabs(V2)=)&(fabs(V3)=)&(fabs(V4)=)&(fabs(V5)=)break;/限差控制。检查计算是否收敛。将所求得的外方位元素的改正数与规定的限差比较。当改正数小于限差时，迭代结束。否则用新的近似值重复计算，直到满足要求为止。/以下是计算单位权中误差doubles=0,m0,v8;xiangchen(A,V,E,8,6,1);/计算AXfor(i=0;i8;i+)vi=Ei-li;/计算AX-Ls+=vi*vi;/计算vvm0=sqrt(s/2);/单位权中误差按下式计算：m?R(t,w,k,a,b,c);/输出结果vv2n?6
