后方交会求外方位元素.docx

《后方交会求外方位元素.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、西南交通大学摄影测量学课程设计报告单像空间后方交会求解外方位元素目录：作业任务…………………………………………………………3计算原理…………………………………………………………3算法流程…………………………………………………………6源程序……………………………………………………………7计算结果…………………………………………………………12结果分析…………………………………………………………13心得体会…………………………………………………………13作业任务：已知条件摄影机主距f=153.24mm，x0=0，y0=0,像片比例尺为1:40000，有四对点的像点坐标与相应的地面坐标如

2、下表。点号像点坐标地面坐标　x(mm)y(mm)X(m)Y(m)Z(m)1-86.15-68.9936589.4125273.322195.172-53.4082.2137631.0831324.51728.693-14.78-76.6339100.9724934.982386.50410.4664.4340426.5430319.81757.31以单像空间后方交会方法，求解该像片的外方位元素。计算原理：如图所示，物点A和摄影中心S在地面摄影测量坐标系中的坐标依次是（X,Y,Z）、（XS,YS,ZS）；像点a在像空间坐标系中的坐标是（x,y,-f）。那么由共线条件方程知，其中ai，

3、bi，ci是只含三个独立参数ψ，ω，κ的九个方向余弦。在方程中共有六个未知参数XS,YS,ZS,ψ,ω,κ，所以有三个不在一条直线上的已知地面点坐标就可以求出像片的这六个外方位元素。由于共线条件方程是非线性方程，为了便于迭代计算，需要把方程用泰勒级数展开，取一次项得到线型表达式，有用新的符号表示各偏导数后为其中（x）、（y）是函数近似值，,是外方位元素近似值的改正数，它们的系数为函数的偏导数。为了便于推导，令===那么有对于系数，其严密算法（以,为例）如下：对于竖直摄影而言，像片的角方位元素都是小值，因而各系数的近似值为为了提高精度和可靠性，通常需要测量四个或更多的地面控制点和对应

4、的像点坐标，采用最小二乘平差方法解算。此时像点坐标（x,y）作为观测值，加入相应的偶然误差改正数,，可列出每个点的误差方程式用矩阵表示为那么由算法流程：Ø获取已知数据m,x0,y0,f,Xtp,Ytp,ZtpØ量测控制点像点坐标x,yØ确定未知数初值Xs0,Ys0,Zs0,j0,w0,k0Ø组成误差方程式并法化Ø解求外方位元素改正数Ø检查迭代是否收敛具体的流程程序框图如下：输入原始数据像点坐标计算，系统误差改正确定外方位元素初始值组成旋转矩阵R结束并显示错误信息是迭代次数小于n逐点组成误差方程式并法化所以像点完否否否解法方程，求外方位元素改正数是计算改正后的外方位元素外方位元素改正

5、数是否小于限差否是输出计算成果，计算结束源程序：程序代码：#include#includeusingnamespacestd;constintN=4;constintn=6;/*---------------矩阵相乘---------------------*/voidmult(double*m1,double*m2,double*result,inti_1,intj_12,intj_2){inti,j,k;for(i=0;i

6、;k++)result[i*j_2+j]+=m1[i*j_12+k]*m2[j+k*j_2];}return;}/*---------------矩阵求逆---------------------*/voidinverse(doublec[n][n]){inti,j,h,k;doublep;doubleq[n][12];for(i=0;i

7、0;k0;k--,h--)//消去对角线以上的数据for(i=k-1;i>=0;i--){if(q[i][h]==0)continue;p=q[k][h]/q[i][h];for