摄影测量实验报告(空间后方交会—前方交会).docx

空间后方交会-空间前方交会程序编程实验实验目的要求掌握运用空间后方交会-空间前方交会求解地面点的空间位置。学会运用空间后方交会的原理根据所给控制点的地面摄影测量坐标系坐标以及相应的像平面坐标系中的坐标利用计算机编程语言实现空间后方交会的过程完成所给像对中两张像片各自的外方位元素的求解。然后根据空间后方交会所得的两张像片的内外方位元素利用同名像点在左右像片上的坐标求解其对应的地面点在摄影测量坐标系中的坐标并完成精度评定过程利用计算机编程语言实现此过程。仪器用具计算机、编程软件（MATLAB）实验数据实验数据包含四个地面控制点（GCP）的地面摄影测量坐标及在左右像片中的像平面坐标。此四对坐标运用最小二乘法求解左右像片的外方位元素即完成了空间后方的过程。另外还给出了5对地面点在左右像片中的像平面坐标和左右像片的内方位元素。实验数据如下：点号左片右片地面摄影测量坐标xyxyXYZGCP116.01279.963-73.9378.7065083.2055852.099527.925GCP288.5681.134-5.25278.1845780.025906.365571.549GCP313.362-79.37-79.122-78.8795210.8794258.446461.81GCP482.24-80.027-9.887-80.0895909.2644314.283455.484151.75881.555-39.95378.463214.618-0.231-76.0160.036349.88-0.792-42.201-1.022486.243-1.346-7.706-2.112548.135-79.962-44.438-79.736内方位元素：f=152.000mmx0=0y0=0实验框图输入GCP的像点坐标xy确定初始值ψ=ω=κ=0XsYsZs计算旋转矩阵R计算像点在像空间坐标系中的近似值（x）（y）并组成误差方程的常数项l计算误差方程的系数项组成系数矩阵A组成法方程式计算系数A’A常数项A’L求解外方位元素计算ψ、ω、κ、Xs、Ys、Zs改正后的值改正数是否小于限差？大于小于计算完毕此过程完成空间后方交会求解像片的外方位元素其中改正数小于限差（0.00003相当于0.1’的角度值）为止。在这个过程中采用迭代的方法是外方位元素逐渐收敛于理论值每次迭代所得的改正数都应加到上一次的初始值之中。输入所需计算点的像平面坐标x1y1；x2y2根据后方交会所得的旋转矩阵R1R2计算像点在左右像空间辅助坐标系中的坐标X1Y1Z1X2Y2Z2计算摄影基线的三个坐标分量BxByBz计算个点在左右像片中的的投影系数N1N2计算地面所求点在地面摄影测量坐标系中的坐标XAYAZA计算完毕精度评定在空间后方交会中运用的数学模型为共线方程确定XsYsZs的初始值时对于左片可取地面左边两个GCP的坐标的平均值作为左片Xs和Ys的初始值取右边两个GCP的坐标平均值作为右片Xs和Ys的初始值。Zs可取地面所有GCP的Z坐标的平均值再加上航高。空间前方交会的数学模型为：实验源代码functionMain_KJQHFJH()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(