单张相片后方交会
目录
●作业任务 (3)
●解算原理 (3)
●具体过程 (4)
●算法描述及程序流程 (4)
●计算结果 (7)
●结果分析 (8)
●心得体会及建议 (8)
●参考文献 (9)
一，作业任务
已知摄影机主距f=153.24mm，四对点的像点坐标与相应地面坐标列入下表：
表1-1
计算近似垂直摄影情况下后方交会解。
二，解算原理
【关键词1】中心投影构像方程
在摄影测量学中，最重要的方程就是中心投影构像方程（图2-1）。这个方程
将地面点在地面摄影测量坐标系中的坐标（物方坐标）和地面点对应像点的像平
面坐标联系起来。在解析摄影测量与数字摄影测量中是极其有用的。在以后将要
学习到的双像摄影测量光束法、解析测图仪原理及数字影像纠正等都要用到该
式。
图2-1
在上述公式中：x和y分别为以像主点为原点的像点坐标，相应地面点坐标
为X,Y,Z，相片主距f以及外方位元素Xs，Ys，Zs，ψ，ω，κ。
而在此次作业中，就是已知四个地面控制点的坐标以及其对应的像点坐标，
通过间接平差原理来求解此张航片的外方位元素。
【关键词2】间接平差
在一个平差问题中，当所选的独立参数X的个数等于必要观测值t时，可将
每个观测值表达成这t个参数的函数，组成观测方程，然后依据最小二乘原理求
解，这种以观测方程为函数模型的平差方法，就是间接平差方法
间接平差的函数模型为：
随机模型为：
平差准则为：VtPV=min
【关键词3】单像空间后方交会
利用至少三个已知地面控制点的坐标A（Xa，Ya，Za）、B（Xb，Yb，Zb）、Z（Xc，
Yc，Zc），与其影像上对应的三个像点的影像坐标a（xa，ya）、b（xb，yb）、c（xc，yc），根据共线方程，反求该像点的外方位元素Xs，Ys，Zs，ψ，ω，κ。这种解算方法是以单张像片为基础，亦称单像空间后方交会。
在此次作业中，就是已知四个控制点在地面摄影测量坐标系中的坐标和对应的像点坐标。由此可以列出8个误差方程，存在两个多余观测数，则n=2。故可利用间接平差里，依据最小二乘法则，进行求解。由于共线条件方程是非线性函数模型，为了便于计算，需把非线性函数表达式用泰勒公式展开成现行形式，即“线性化”。而又因为仅取泰勒级数展开式的一次项，未知数的近似值改正是粗略的，所以必须计算采用逐渐趋近法，解求过程需要反复趋近，直至改正值小于限差为止。
三，具体过程
1，获取已知点数据：从摄影资料中查取像片比例尺1/m，平均航高，内方元素x0，y0，f；从外业测量成果中，获取控制点的地面测量坐标Xt，Yt，Zt，并转换成摄影测量坐标X,Y,Z。
2，量测控制点的像点坐标：将控制点标刺在像片上，利用立体坐标量测仪量测控制点的像框坐标，并经像点坐标改正，得到像点坐标x，y。
3，确定未知数的初始值：在竖直摄影测量情况下，角元素的初始值为0，及ψ=ω=κ=0；
线元素中，Zso =m*f+(Z[0]+Z[1]+Z[2]+Z[3])/4，Xso，Yso的取值可用四个角点上制点坐标的平均值，即：Xso=(X[0]+X[1]+X[2]+X[3])/4;Yso=(Y[0]+Y[1]+Y[2]+Y[3])/4；4，计算旋转矩阵R：利用角元素的近似值计算方向余弦，组成R阵。公式如下：R[0][0]=cos(ψ)*cos(k)-sin(ψ)*sin(w)*sin(k);
R[0][1]=-cos(ψ)*sin(k)-sin(ψ)*sin(w)*cos(k);
R[0][2]=-sin(ψ)*cos(w);
R[1][0]=cos(w)*sin(k);
R[1][1]=cos(w)*cos(k);
R[1][2]=-sin(w);
R[2][0]=sin(ψ)*cos(k)+cos(ψ)*sin(w)*sin(k);
R[2][1]=-sin(ψ)*sin(k)+cos(ψ)*sin(w)*cos(k);
R[2][2]=cos(ψ)*cos(w);
5，逐点计算像点坐标的近似值：利用未知数的近似值按共线方程计算控制点像点坐标的近似值（x）、（y）；
6，组成误差方程式：参照教材（5-8）式、（5-9b）式、（5-4）式逐点计算误差方程的系数阵和常数项。
7，组成法方程：计算法方程的系数矩阵与常数项。
8，解求外方位元素：根据法方程，按间接平差原理解求外方位元素改正值，并与相应的近似值求和，得到外方位元素的新的近似值。
9，检查计算是否收敛：将求得的外方位元素的改正值与规定的限差比较，小于限差则计算终止，否则用新的近似值重复第4至第8步骤计算，直至满足要求为止。
四，算法描述及程序流程。
算法描述（图4-1）：
图4-1
程序流程：
1，读入数据：新建一个txt格式的文本文档（图4-1），在程序中通过读入数据操作将已
知数据导入到程序中。
图4-2
实现这一操作的语句为：
ifstream shuju;
shuju.open(“原始坐标数据.txt”);
if(shuju.is_open())
{
while(!shuju.eof ())
{
for(int i=0;i<4;i++)
{
shuju >>x[i]; shuju;x[i]=(x[i]-x0)/1000;
shuju >>y[i]; shuju;y[i]=(y[i]-y0)/1000;
shuju >>X[i]; shuju;
shuju >>Y[i]; shuju;
shuju >>Z[i]; shuju; //分列读入数据
}
}yuanshi.close();
}
2，给六个外方位元素赋初值，并计算旋转矩阵R。实现此操作的源程序如下：R[0][0]=cos(q)*cos(k)-sin(q)*sin(w)*sin(k);
R[0][1]=-cos(q)*sin(k)-sin(q)*sin(w)*cos(k);
R[0][2]=-sin(q)*cos(w);
R[1][0]=cos(w)*sin(k);
R[1][1]=cos(w)*cos(k);
R[1][2]=-sin(w);
R[2][0]=sin(q)*cos(k)+cos(q)*sin(w)*sin(k);
R[2][1]=-sin(q)*sin(k)+cos(q)*sin(w)*cos(k);
R[2][2]=cos(q)*cos(w);
3，计算系数阵A和常数项L 。实现此操作程序如下：
for(int i=0,J=0,H=0;i<=3;i++)
{
XX=R[0][0]*(X[i]-Xs)+R[1][0]*(Y[i]-Ys)+R[2][0]*(Z[i]-Zs);
YY=R[0][1]*(X[i]-Xs)+R[1][1]*(Y[i]-Ys)+R[2][1]*(Z[i]-Zs);
ZZ=R[0][2]*(X[i]-Xs)+R[1][2]*(Y[i]-Ys)+R[2][2]*(Z[i]-Zs);
xx=-f*XX/ZZ;
yy=-f*YY/ZZ;
l[J][0]=x[i]-xx;
l[J+1][0]=y[i]-yy;
J=J+2;
A[H][0]=(R[0][0]*f+R[0][2]*x[i])/ZZ;
A[H][1]=(R[1][0]*f+R[1][2]*x[i])/ZZ;
A[H][2]=(R[2][0]*f+R[2][2]*x[i])/ZZ;
A[H][3]=y[i]*sin(w)-(x[i]*(x[i]*cos(k)-y[i]*sin(k))/f+f*cos(k))*cos(w);
A[H][4]=-f*sin(k)-x[i]*(x[i]*sin(k)+y[i]*cos(k))/f;
A[H][5]=y[i];
A[H+1][0]=(R[0][1]*f+R[0][2]*y[i])/ZZ;
A[H+1][1]=(R[1][1]*f+R[1][2]*y[i])/ZZ;
A[H+1][2]=(R[2][1]*f+R[2][2]*y[i])/ZZ;
A[H+1][3]=-x[i]*sin(w)-(y[i]*(x[i]*cos(k)-y[i]*sin(k))/f-f*sin(k))*cos(w);
A[H+1][4]=-f*cos(k)-y[i]/f*(x[i]*sin(k)+y[i]*cos(k));
A[H+1][5]=-x[i];
H=H+2;
}
4，根据间接平差原理解算出改正数并与限差（在此程序中，限差是三个角度的改正值的
绝对值小于6秒）进行比较，若小于限差则计算完成，否则计算新的外方位元素的近似值，生成新的R阵，系数阵和L阵，重新平差。为实现这一迭代计算，程序中采用了一个do while循环语句。
While的循环条件是fabs(XF[3][0])*t>=6.0 || fabs(XF[4][0])*t>=6.0 || fabs(XF[5][0])*t>=6.0。（其中XF[3][0],XF[4][0],XF[5][0]分别为外方位元素中三个角元素的改正值，fabs是取绝对值）
2，精度估计及生成解析结果文档（.txt）。
另：矩阵相乘及求逆的子函数见附件。
五，计算结果
2，计算结果：
迭代次数为：3次
此张相片的外方位元素为：
三个线元素：
Xs=39795.5m；Ys=27476.5m；Zs=7572.69m
三个角元素：
ψ=-14.353″；ω=7.61006″；κ=-243.281″
精度评定如下：
Xs: ±1.10739mm；Ys: ±1.24952mm；Zs: ±0.48813mm
ψ:±0.000178625″；ω:±0.000161461″；κ:±7.20381e-005″
2，生成结果的文本文档如下图（5-1）所示：
图5-1
程序调试结果如下图（5-2）：
图5-1
六、结果分析
此次单像后方交会的迭代次数为三次，符合要求。可见此次解算六个外方位元素初始值取得很好。通过对外方位元素的最终结果分析，角度近似垂直，并且线元素改正值也不大，所以此张相片的摄影为近似垂直摄影。
通过精度评估发现，此次观测角度的精度远高于线元素的观测精度。
七，心得体会及建议
1，“尽信书则不如无书”，此次编程中我发现了书上(摄影测量学·武汉大学出版社·2005年9月)的两个错误的公式，一个是P66页上关于Zso确定的公式。书上的是Zso=m*f，而正确的是Zso =m*f+(Z[0]+Z[1]+Z[2]+Z[3])/4。还有一个是P66页关于系数阵的一个公式。书上的公式是a24=-xsinw-【x*(x*cosk-ysink)/f-f*sink】cosw，而正确的公式应该是a24=-xsinw-【y*(x*cosk-ysink)/f-f*sink】cosw。由于这两个公式的错误，给我的编程造成了很大的困难。后来才发现这两个公式有误，程序才能继续进行下去。
2，闭门造车是造不出好车的，要集思广益。此次编程过程中，我与测绘一班的李志鹏同学积极讨论，请教学长。后来通过分析张勇学长的程序，借鉴了张勇学长关于系数阵和L阵的构造。当时自己没有想到通过数组来完成这个操作，只是单纯地输入每个系数，共48个数据，不仅工作量大，而且容易出错。不是明智之举。以后我会把这个思想应用在后来的学习和工作中。当然，还不得不提及关于读入文件和生成结果文件的操作。在程序的开端，我通过文件输入流，读入一个txt格式的文本文档，在程序的末尾，又生成了一个txt的结果文档。其实这个是比较模块化的东西，难度不大。
3，还有不得不提及的一点，这就是在此次编程过程中，关于矩阵相乘和矩阵求逆子函数调用的操作实现。首先，我自己编写了一个关于矩阵相乘的函数，在网上找了一个
矩阵求逆的子函数。可是在实际应用时，老是提示配型不正确。通过讨论，我们发现，将所有矩阵都定义成8*8的，而在调用的时候，通过实参传递给形参实际的矩阵维数，这样就能解决矩阵配型不正确的问题。实在是获益匪浅，这个技巧对于以后的编程很有用。
4，其实，看来学长的程序才知道自己在程序编写方面还不行，现在要好好补补课。虽然在原理方面掌握得比较扎实，但是真正应用起来还是有问题的。当我的程序多次不能调试出正确的结果的时候，自己都开始怀疑原理上是不是有问题了。其实，学长的程序做得很好了，在做出一些什么有新意的东西也不大可能了。只是通过这次编程，我对于函数的调用这块，有很大的收获。并且在构造系数阵的思维方面有了突破。这个算是我的第一个小成就吧，我会再接再厉。在以后的学习中更上一层楼！
八，参考文献
（1）王佩军、徐亚明著·摄影测量学·武汉大学出版社·2005年9月
（2）杨先宏、邵鸿潮等著·摄影测量学·武汉大学出版社·1996年8月
（3）景红编·大学计算机程序设计基础·西南交通大学出版社·2008年3月
（4）谭浩强编·C++程序设计·清华大学出版社·2004年6月
空间后方交会的解算 一． 空间后方交会的目的 摄影测量主要利用摄影的方法获取地面的信息，主要是是点位信息，属性信息，因此要对此进行空间定位和建模，并首先确定模型的参数，这就是空间后方交会的目的，用以求出模型外方位元素。 二． 空间后方交会的原理 空间后方交会的原理是共线方程。 共线方程是依据相似三角形原理给出的，其形式如下 111333222333()()() ()()() ()()()()()()A S A S A S A S A S A S A S A S A S A S A S A S a X X b Y Y c Z Z x f a X X a Y Y a Z Z a X X b Y Y c Z Z y f a X X a Y Y a Z Z -+-+-=--+-+--+-+-=--+-+- 上式成为中心投影的构线方程， 我们可以根据几个已知点，来计算方程的参数，一般需要六个方程，或者要三个点，为提高精度，可存在多余观测，然后利用最小二乘求其最小二乘解。 将公式利用泰勒公式线性化，取至一次项，得到其系数矩阵A ；引入改正数（残差）V ，则可将其写成矩阵形式： V AX L =- 其中 111333222333[,]()()()()()()()()()()()()()()T x y A S A S A S x A S A S A S A S A S A S y A S A S A S L l l a X X b Y Y c Z Z l x x x f a X X a Y Y a Z Z a X X b Y Y c Z Z l y y y f a X X a Y Y a Z Z =-+-+-=-=+-+-+--+-+-=-=+-+-+- 则1()T T X A A A L -= X 为外方位元素的近似改正数， 由于采用泰勒展开取至一次项，为减少误差，要将的出的值作为近似值进行迭代，知道小于规定的误差 三． 空间后方交会解算过程 1. 已知条件 近似垂直摄影
摄影测量原理 单张影像后方交会实习
目录 一实习目的 (3) 二实习原理 (3) 1. 间接平差 (3) 2. 共线方程 (3) 3. 单向空间后方交会 (4) 三计算流程 (4) 1. 求解步骤 (4) 2.计算机框图 (4) 四程序实现 (5) 五结果分析 (6) 1.外方位元素 (6) 2.误差 (6) 3.旋转矩阵R (7) 六实习体会 (7) 1. 平台的选择 (7) 2.问题的解决 (7) 3.心得体会 (8) 七代码展示 (8)
一实习目的 为了增强同学们对后方交会公式的理解，培养同学们对迭代循环编程的熟悉感，本次摄影测量课间实习内容定为用C语言或其他程序编写单片空间后方交会程序，最终输出像点坐标、地面坐标、单位权中误差、外方位元素及其精度。 已知四对点的影像坐标和地面坐标如下。内方位元素fk=153.24mm，x0=y0=0。 本次实习，我使用了matlab2014进行后方交会程序实现。结果与参考答案一致，精度良好。 二实习原理 题干中有四个控制点在地面摄影测量坐标系中的坐标和对应的像点坐标，由此可列出8个误差方程，存在2个多余观测（n=2）。故可利用间接平差的最小二乘法则求解。 由于共线方程是非线性函数模型，为了方便计算，需要将其“线性化”。但如果仅取泰勒级数展开式的一次项，未知数的近似值改正是不精确的。因此必须采用迭代趋近法计算，直到外方位元素的改正值小于限差。 1.间接平差 间接平差为平差计算最常用的方法。在确定多个未知量的最或然值时，选择它们之间不存在任何条件关系的独立量作为未知量组成用未知量表达测量的函数关系、列出误差方程式，按最小二乘法原理求得未知量的最或然值的平差方法。 在一个间接平差问题中，当所选的独立参数X个数与必要观测值t个数相等时，可将每个观测值表达成这t个参数的函数，组成观测方程。 函数模型为：L = BX + d。 2.共线方程 共线方程是中心投影构像的数学基础，也是各种摄影测量处理方法的重要理论基础。 式中： x，y 为像点的像平面坐标； x0，y0，f 为影像的内方位元素； XS，YS，ZS 为摄站点的物方空间坐标； XA，YA，ZA 为物方点的物方空间坐标； ai，bi，ci （i = 1,2,3）为影像的3 个外方位角元素组成的9 个方向余弦。
摄影测量学 第一章 1、摄影测量的定义、任务？ 定义：摄影测量与遥感是从非接触成像和其他传感器系统，通过记录、量测、分析与表达等处理，获取地球及其环境和其他物体可靠信息的工艺，科学与技术。其中摄影测量侧重于提取几何信息，遥感侧重于物理信息。任务：（1）测绘各种比例尺地形图。（2）建立数字地面模型（地形数据库）。 2、摄影测量学：是对研究的对象进行摄影，根据所得的构象信息，从几何方面和物理方面 加以分析研究，从而对所摄对象的本质提供各种资料的一门学科。 3、解决的基本问题：几何定位和影像解译。 4、摄影测量的三个发展阶段及其特点。 1、航摄仪物镜的焦距与其主距有什么不同？ 焦距：自物方主点S1到物方焦点F1的距离称为光学系的物方焦距f1；自像方主点S2到像方焦点F2的距离称为镜头的像方焦距f2。 主距：像主点和摄影机物镜后节点之间的距离称为摄影机主距。 2、量测摄影机与非量测摄影机的区别？ （1）量测摄影机的主距是一个固定的已知值 （2）量测摄影机的承片框上具有框标，即固定不变的承片框上，四个边的中点各安置一个机械标志；框标，其目的是建立像片的直角，框标坐标系。 （3）量测摄影机的内方位元素是已知值。 3、航向重叠：摄影时飞机沿相邻影像之间必须保持一定的重叠度。一般P=50%~65%；P值最小不能小于53%。 旁向重叠：完成一条航线的摄影后，飞机进入另一条航线进行测量摄影，相邻航线影像之间也必须有一定的重叠。一般q=30%~40%，最小不得小于15%。 4、B与近景C之间这一段间隔内的所有景物，在像面上仍可获得清晰的图像，此时，近景 与远景之间的纵深度称为景深。 5、超焦点距离：当物镜向无限远物体对光时，不仅远处的物体构象清晰，而且在离开物镜 不小于某一距离H的所有物体，其构象都很清晰，这个距离H就称为超焦点距离。 第三章 1、航摄像片上特殊的点、线、面。 （1）像主点：摄影中心S在像片平面上的投影点。 （2）像底点：主垂线与像片面P的交点n称为像底点。 （3）等角点：倾角α的平分线与像片面交于点C称C点为等角点。 （4）主纵线：主垂面W与像平面P的交线称为主纵线W。 （5）等比线：过像主点平行于合线的直线称为等比线。 2、摄影测量常用的坐标系统，它们是如何定义的？ （1）像平面坐标系：是以该像片的像主点为坐标原点的坐标系，用来表示像点在像片面上的位置，在实际应用中，常采用框标连线的交点为坐标原点，称为框标平面坐标系。X、y轴的方向按需要而定，常取与航线方向一致的连线为x轴，航线方向为正。 （2）像空间坐标系：以摄影中心S为坐标原点，X轴和Y轴分别与像平面直角坐标系的X轴和Y轴平行，Z轴与主光轴重合，向上为正，像点的像空间坐标系表示为（x、y、-f）。 （3）像空间辅助坐标系：其坐标原点是摄影中心S坐标轴依情况而定，通常有三种方法：a、以每一条航线的第一张像片的像空间坐标系作为像空间辅助坐标系。 b、取u、v、w轴系分别平行于地面摄影测量坐标系D-XYZ，这样同一像点a在像空间坐标系中的坐标为x、y、z=（-f），而在像空间辅助坐标系中的坐标为u、v、w。 c、以每个像片对的左片摄影中
单张相片后方交会
目录 ●作业任务 (3) ●解算原理 (3) ●具体过程 (4) ●算法描述及程序流程 (4) ●计算结果 (7) ●结果分析 (8) ●心得体会及建议 (8) ●参考文献 (9)
一，作业任务 已知摄影机主距f=153.24mm，四对点的像点坐标与相应地面坐标列入下表： 表1-1 计算近似垂直摄影情况下后方交会解。 二，解算原理 【关键词1】中心投影构像方程 在摄影测量学中，最重要的方程就是中心投影构像方程（图2-1）。这个方程 将地面点在地面摄影测量坐标系中的坐标（物方坐标）和地面点对应像点的像平 面坐标联系起来。在解析摄影测量与数字摄影测量中是极其有用的。在以后将要 学习到的双像摄影测量光束法、解析测图仪原理及数字影像纠正等都要用到该 式。 图2-1 在上述公式中：x和y分别为以像主点为原点的像点坐标，相应地面点坐标 为X,Y,Z，相片主距f以及外方位元素Xs，Ys，Zs，ψ，ω，κ。 而在此次作业中，就是已知四个地面控制点的坐标以及其对应的像点坐标， 通过间接平差原理来求解此张航片的外方位元素。 【关键词2】间接平差 在一个平差问题中，当所选的独立参数X的个数等于必要观测值t时，可将 每个观测值表达成这t个参数的函数，组成观测方程，然后依据最小二乘原理求 解，这种以观测方程为函数模型的平差方法，就是间接平差方法 间接平差的函数模型为： 随机模型为： 平差准则为：VtPV=min 【关键词3】单像空间后方交会 利用至少三个已知地面控制点的坐标A（Xa，Ya，Za）、B（Xb，Yb，Zb）、Z（Xc，
Yc，Zc），与其影像上对应的三个像点的影像坐标a（xa，ya）、b（xb，yb）、c（xc，yc），根据共线方程，反求该像点的外方位元素Xs，Ys，Zs，ψ，ω，κ。这种解算方法是以单张像片为基础，亦称单像空间后方交会。 在此次作业中，就是已知四个控制点在地面摄影测量坐标系中的坐标和对应的像点坐标。由此可以列出8个误差方程，存在两个多余观测数，则n=2。故可利用间接平差里，依据最小二乘法则，进行求解。由于共线条件方程是非线性函数模型，为了便于计算，需把非线性函数表达式用泰勒公式展开成现行形式，即“线性化”。而又因为仅取泰勒级数展开式的一次项，未知数的近似值改正是粗略的，所以必须计算采用逐渐趋近法，解求过程需要反复趋近，直至改正值小于限差为止。 三，具体过程 1，获取已知点数据：从摄影资料中查取像片比例尺1/m，平均航高，内方元素x0，y0，f；从外业测量成果中，获取控制点的地面测量坐标Xt，Yt，Zt，并转换成摄影测量坐标X,Y,Z。 2，量测控制点的像点坐标：将控制点标刺在像片上，利用立体坐标量测仪量测控制点的像框坐标，并经像点坐标改正，得到像点坐标x，y。 3，确定未知数的初始值：在竖直摄影测量情况下，角元素的初始值为0，及ψ=ω=κ=0； 线元素中，Zso =m*f+(Z[0]+Z[1]+Z[2]+Z[3])/4，Xso，Yso的取值可用四个角点上制点坐标的平均值，即：Xso=(X[0]+X[1]+X[2]+X[3])/4;Yso=(Y[0]+Y[1]+Y[2]+Y[3])/4；4，计算旋转矩阵R：利用角元素的近似值计算方向余弦，组成R阵。公式如下：R[0][0]=cos(ψ)*cos(k)-sin(ψ)*sin(w)*sin(k); R[0][1]=-cos(ψ)*sin(k)-sin(ψ)*sin(w)*cos(k); R[0][2]=-sin(ψ)*cos(w); R[1][0]=cos(w)*sin(k); R[1][1]=cos(w)*cos(k); R[1][2]=-sin(w); R[2][0]=sin(ψ)*cos(k)+cos(ψ)*sin(w)*sin(k); R[2][1]=-sin(ψ)*sin(k)+cos(ψ)*sin(w)*cos(k); R[2][2]=cos(ψ)*cos(w); 5，逐点计算像点坐标的近似值：利用未知数的近似值按共线方程计算控制点像点坐标的近似值（x）、（y）； 6，组成误差方程式：参照教材（5-8）式、（5-9b）式、（5-4）式逐点计算误差方程的系数阵和常数项。 7，组成法方程：计算法方程的系数矩阵与常数项。 8，解求外方位元素：根据法方程，按间接平差原理解求外方位元素改正值，并与相应的近似值求和，得到外方位元素的新的近似值。 9，检查计算是否收敛：将求得的外方位元素的改正值与规定的限差比较，小于限差则计算终止，否则用新的近似值重复第4至第8步骤计算，直至满足要求为止。 四，算法描述及程序流程。 算法描述（图4-1）：
名词解释 1空中三角测量：利用航摄像片与所摄目标之间的空间几何关系，根据少量像片控制点，计算待求点的平面位置、高程和像片外方位元素的测量方法。 2像点位移：由于在实际航空摄影时，在中心投影的情况下，当航摄的飞行姿态出现较大倾斜即像片有倾斜，地面有起伏时，便会导致地面点在航摄像片上构像相对于在理想情况下的构像，产生了位置的差异，这一差异称为像点位移。 3摄影基线：航线方向相邻两个摄影站点间的空间距离。 4航向重叠：同一条航线上，相邻两张像片应有一定范围的影像重叠，称为航向重叠。 5旁向重叠：相邻航线相邻两像片的重叠度 6同名核线：同一核面与左右影像相交形成的两条核线，其中核面指物方点与摄影基线所确定的平面。 7像片的内方位元素：表示摄影中心与像片之间相互位置的参数，f,x0,y0 8像片的外方位元素：表示摄影中心和像片在地面坐标系中的位置和姿态的参数。 9相对定向：根据立体像对内在的几何关系恢复两张像片之间的相对位置和姿态，使同名光线对对相交，建立与地面相似的立体模型。即确定一个立体像对两像片的相对位置。 10绝对定向元素：描述立体像对在摄影瞬间的绝对位置和姿态的参数。 11单像空间后方交会：利用至少三个已知地面控制点的坐标，与其影像上对应三个像点的影像坐标，根据共线条件方程，反求该像片的外方位元素。 12空间前方交会：由立体像对中两张像片的内、外方位元素和像点坐标来确定相应地面点的地面坐标的方法，称为空间前方交会。 13同名像点：同名光线在左右相片上的构像 填空 1、4D 产品是指 DEM 、DLG 、DRG 、DOM 。 2、摄影测量按用途可分为地形摄影测量、非地形摄影测量。 3、摄影测量学的发展经过了模拟摄影测量、解析摄影测量、数字摄影测量三个阶段。 4、模拟摄影测量是利用光学/机械投影方法实现摄影过程的反转。 5、解析摄影测量以电子计算机为主要手段，通过对摄影像片的量测和解析计算方法的交会方式来研究和确定被摄物体的形状、大小、位置、性质及其相互关系，并提供各种摄影测量产品的一门科学。 6、像点坐标的系统误差改正主要包括底片变形改正，摄影机物镜畸变差改正，大气折光改正和地球曲率改正。 7、共线方程表达的是像点、投影中心与地面点之间关系。 8、立体摄影测量基础是共面条件方程。 9、把一条航线的航摄像片根据地物影像拼接起来，各张像片的主点连线不在一条直线上，而呈现为弯弯曲曲的折线，称航线弯曲。 10、航摄像片为量测像片，有光学框标和机械框标。 11、地图是地面的正射投影，像片是地面的中心投影。 12、在像空间坐标系中，像点的z 坐标值都为-f 。 13、一张像片的外方位元素包括：三个直线元素(Xs 、Ys 、Zs ）：描述摄影中心的空间坐标值；三个角元素(?、ω、κ) ) ：描述像片的空间姿态。 14、相对定向的理论基础、目的、标准是两像片上同名像点的投影光线对对相交。 15、双像解析摄影测量的任务是利用解析计算方法处理立体像对，获取地面点的三维空间信息。 16、在摄影测量中，一个立体像对的同名像点在各自的像平面坐标系的x 、y 坐标之差，分别称为左右视差、上下视差。 17、解析法相对定向的理论基础是同名光线对对相交于核面内。 18、解析绝对定向需要量测 2 个平高和 1 个高程以上的控制点，一般是在模型四个角布设四个控制点。 19、解析空中三角测量按数学模型分为航带法、独立模型法、光束法。 20、像底点上不存在投影差，但存在倾斜误差。倾斜航片上等比线上点的倾斜误差等于零。 21、立体模型空间相对定向时，连续像对的相对定向元素为 ，单独像对的相对定向元素为 。 22、某像点的像平面坐标为(x,y)，摄影仪主距为f ，则该点在像空间坐标系中的坐标为(x ，y ，-f )。 23、摄影测量采用的五种常用坐标系中，地面测量坐标系是左手系。 222 v w b b φωκ、、、、22211ωκ?κ?、、、、
4D 产品是指DEM、DLG、DRG、DOM。 摄影测量学:是利用光学摄影机摄取照片，通过像片来研究和确定被摄物体的形状大小位置和相互关系的一门科学技术摄影测量按远近分为航天摄影测量、航空摄影测量，地面摄影测量，近景摄影测量，显微镜摄影测量。 摄影测量按用途口J分为地形摄影测量、非地形摄影测量。 摄影测量学的发展经过了模拟摄影测量、解析摄影测量、数字摄影测量三个阶段。 2.由于立体像对选取的像空间辅助坐标系的不同分为连续邃戒与里独像对 摄影机的主距：航空摄影物镜中心至底片面的距离是固定值1?摄影比例尺：严格讲，摄影比例尺是指航摄像片上一线段为J与地向上相应线段的水干距L之比。摄影像片的影像比例尺处处均不相等 3?摄影航高：摄影机的物镜中心至该面的距离 2?绝对航高:摄影物镜相对于平均海平而的航高，指摄影物镜在摄影瞬间的真实海拔高度。 3?相对航高：摄影物镜相对于某一基准面或某一点的高度 2 ?制定航摄计划： 1.确定摄区范围； 2.选择航摄仪； 3.确定航摄仪的比例尺；4,确定摄影航高；5,需要像片数，F1期等。 5.摄影基线：航线方向相邻两摄站点间的空间距离称为摄影基线。 2?摄影资料的基本要求：1.影像的色调，2.像片的重叠，3.像片倾角，4.航线弯曲，5,像片旋角 2?像片倾角：空中摄影采用竖直摄影方式，即摄影瞬间摄影机的主光轴近似与地面垂直，它偏离铅垂线的夹角应小于3D,夹角称为像片倾角。 3?航向重叠：同一条航线内相邻像片之间的影像重叠称为航向重叠，一般要求在60%以上。目的：保 证像片立体量测与拼接 4?旁向重叠：相邻航线的重叠称为旁向重叠，重叠度要求在24%以上 5?中心投影：投影光线会聚与一点 7?像主点：摄影机主光轴在框标平面上的垂足 &像底点：主垂线与像片面的交点 2 ?摄影测量常用的坐标系统有哪些? 像平面坐标系；像空间坐标系；像空间辅助坐标系；摄影测量坐标系；地面测量坐标系 3.对于一张航摄像片其内外方位元素为内外方位元素均为常数， 8?内方位元素：内方位元素是表示摄影中心与像片之间相关位置的参数，包括三个参数。即摄影中心 到像片的垂距（主距）f及像主点o在像框标坐标系中的坐标兀。,儿 9?外方位元素：在恢复内方位元素的基础上，确定摄影光束在摄影瞬间的空间位置与姿态的参数称为外方位元素， 外方位角元素：确定像空间坐标系的三轴在地面坐标系中的方向。 14 ?像点在像空间直角坐标系与像空间辅助坐标系的变换关系： U X坷a2 a3X V=R y—久b2伏y W-f°1 C2°3-f 13?同名像点：同名光线在左右相片上的构像 14 ?摄影基线：同一航线内相邻两摄站的连线 15?核线：核面与像片的交线，核线会聚于核点 16?核面：摄影基线与地而点所作平而 17.同名像点：地面上一点在相邻两张像片上的构像
摄影测量后方交会程序（c/c++） 输入数据截图： 结果截图： 程序源代码（其中的矩阵求逆在前面已经有了，链接）：
#include
} printf("解算完毕...\n"); do{ printf("计算结果保存于\"结果.txt\"文件中\n" "请选择操作（输入P打开结果数据，R打开原始数据，其它退出程序）："); fflush(stdin); //刷新输入流 char order=getchar(); if ('P'==order || 'p'==order) { system("结果.txt"); } else if ('R'==order || 'r'==order) { system("data.txt"); } else break; system("cls"); }while(1); system("PAUSE"); return 0; } /********************************************** *函数名：InputData *函数介绍：从文件（data.txt）中读取数据, *文件格式如下： *点数 m(未知写作0) * 内方位元素(f x0 y0) *编号 x y X Y Z *下面是一个实例： 4 0 153.24 0 0 1 -86.15 -68.99 36589.41 25273.3 2 2195.17
摄影测量学习题 一、名词解释： 1、摄影测量学：是对研究的对象进行摄影，根据所获得的构像信息，从几何方面和物理方 面加以分析研究，从而对所摄对象的本质提供各种资料的一门学科。 2、光圈号数 ：相对孔径的倒数 3、景深 ：远景与近景之间的纵深距离称为景深 4、超焦点距离：当物镜向无限远物体对光时，不仅远处的物体构象清晰，而且在离开物镜 不小于某一距离H 的所有物体，其构象都很清晰，这个距离H 就称为超焦点距离或称为无限远起点 5、视场： 将物镜对光于无穷远，在焦面上会看到一个照度不均匀的明亮圆。这个直径为 ab 的明亮圆的范围称为视场 6、视场角 ：物镜的像方主点与视场直径所张的角2α。 7、像场 ：在视场面积内能获得清晰影像的区域 8、像场角； 物镜的像方主点与像场直径所张的角2β。 像主点：摄影机轴在框标平面上的垂足。 11、航向重叠 ：沿飞行方向上相邻像片所摄地面的重叠区。 12、旁向重叠：两相邻航带摄区之间的重叠 主光轴 ：通过诸透镜光轴的轴 主点： 主平面与光轴的交点 13、摄影基线 ：相邻像片摄影站(投影中心)之间的空间连线。 15、内方位元素 确定物镜后节点和像片面相对位置的数据。 16、外方位元素 确定摄影摄影机或像片的空间位置和姿态的参数 焦点 平行光轴的投射光线经物镜后产生折射，该折射线与光轴的交点。 17、像片倾角 航摄仪光轴与通过物镜中心的铅垂线所夹的角称为像片的倾斜角 19、像片旋角 相邻像片的主点连线与像幅沿航线方向两框标连线之间的夹角称为像片的旋 偏角 20、倾斜误差 因像片倾斜引起的像点位移 节点 投射光线与成像光线与光轴的交角u 和u ′相等时，投射光线与成像光线与光轴的交点。 21、投影差 因地形起伏引起的像点位移 22、摄影比例尺 航摄相片上某一线段构成的长度与地面上相应水平距离之比。 23、像片控制点 为联系地面与相片而测定地面坐标的像点。 相对孔径 物镜焦距与有效孔径之比 25、左右视差 同名像点在各自像平面坐标系中的x 坐标之差 26、上下视差 同名像点在各自像平面坐标系中的Y 坐标之差 27、核点 基线延长线与左、右像片的交点k 1、k 2称为核点 28、核线 核面与像片的交线称为核线 29、核面 通过摄影基线S 1S 2与任一地面点A 所作的平面W A 30、投影基线 两摄站的连线 31、像片基线 指相邻两张像片主点的连线 32、解析空中三角测量 即在一条航带几十条像对覆盖的区域或由几条航带几百哥像对构成 的区域内，仅仅由外业实测几个少量的控制点，按一定的数学模型，平差 解算出摄影测量作业过程中所需的全部控制点及每张像片的外方位元素 33、空间后方交会 就是利用地面控制点的已知坐标值反求像片外方位元素 ()()()()(){} 2332233213322232332 1[]Z X Y X Y Y Y X X X Z Y X X Y Z X Y Y X Z X Y X Y =-+-+-+-+--
程序运行环境为Visual Studio2010.运行前请先将坐标数据放在debug 下。 1.单像空间后方交会 C语言程序： #include 单像空间后方交会和双像解析空间后方-前 方交会的算法程序实现 遥感科学与技术 摘要：如果已知每张像片的6个外方位元素，就能确定被摄物体与航摄像片的关系。因此，利用单像空间后方交会的方法，可以迅速的算出每张像片的6个外方位元素。而前方交会的计算，可以算出像片上点对应于地面点的三维坐标。基于这两点，利用计算机强大的运算能力，可以代替人脑快速的完成复杂的计算过程。 关键词：后方交会，前方交会，外方位元素，C++编程 0.引言： 单张像片空间后方交会是摄影测量基本问题之一,是由若干控制点及其相应像点坐标求解摄站参数(X S,Y S,ZS,ψ、ω、κ)。单像空间后方交会主要有三种方法:基于共线条件方程的平差解法、角锥法、基于直接线性变换的解法。而本文将介绍第一种方法，基于共线条件方程反求象片的外方位元素。 而空间前方交会先以单张像片为单位进行空间后方交会，分别求出两张像片的外方位元素，再根据待定点的一对像点坐标，用空间前方交会的方法求解待定点的地面坐标。可以说，这种求解地面点的坐标的方法是以单张像片空间后方交会为基础的，因此，单张像片空间后方交会成为解决这两个问题以及算法程序实现的关键。 1.单像空间后方交会的算法程序实现： （1）空间后方交会的基本原理：对于遥感影像,如何获取像片的外方位元素,一直是摄影测量工作者探讨的问题,其方法有:利用雷达(Radar)、全球定位系统(GPS)、惯性导航系统(I N S)以及星像摄影机来获取像片的外方位元素;也可以利用一定数量的地面控制点,根据共线方程,反求像片的外方位元素,这种方法称为单像空间后方交会(如图1所示)。 图中,地面坐标X i、Yi、Zi和对应的像点坐标x i、yi是已知的,外方位元素XS、Y S、ZS(摄站点坐标),ψ、ω、κ(像片姿态角)是待求的。 （2）空间后方交会数学模型：空间后方交会的数学模型是共线方程, 即中心投影的构像方程： 式中X、Y、Z是地面某点在地面摄影测量坐标系中的坐标,x,y是该地面点在像片上的构像点的像片坐标,对 于空间后方交会而言它们是已知的,还有主距f是已知的。而9个方向余弦a 1,a 2,a3;b1,b 2,b 3;c 1,c2,c 3是未知的，具体表达式可以取 摄影测量学习题 一、名词解释: １、摄影测量学：是对研究的对象进行摄影，根据所获得的构像信息，从几何方面和物理方 面加以分析研究，从而对所摄对象的本质提供各种资料的一门学科。 ２、光圈号数?:相对孔径的倒数 ３、景深 :远景与近景之间的纵深距离称为景深 ? ４、超焦点距离:当物镜向无限远物体对光时,不仅远处的物体构象清晰，而且在离开物镜不 小于某一距离H 的所有物体,其构象都很清晰,这个距离Ｈ就称为超焦点距离或称为无限远起点 5、视场:?将物镜对光于无穷远,在焦面上会看到一个照度不均匀的明亮圆。这个直径为ａb 的明亮圆的范围称为视场 6、视场角 :物镜的像方主点与视场直径所张的角2α。 ? 7、像场?：在视场面积内能获得清晰影像的区域 ? 8、像场角; 物镜的像方主点与像场直径所张的角2β。 像主点:摄影机轴在框标平面上的垂足。 11、航向重叠?:沿飞行方向上相邻像片所摄地面的重叠区。 ? 12、旁向重叠:两相邻航带摄区之间的重叠? 主光轴?:通过诸透镜光轴的轴? 主点: 主平面与光轴的交点 13、摄影基线?：相邻像片摄影站（投影中心)之间的空间连线。 15、内方位元素 确定物镜后节点和像片面相对位置的数据。? 16、外方位元素 确定摄影摄影机或像片的空间位置和姿态的参数 焦点 平行光轴的投射光线经物镜后产生折射,该折射线与光轴的交点。 １7、像片倾角 航摄仪光轴与通过物镜中心的铅垂线所夹的角称为像片的倾斜角 19、像片旋角 相邻像片的主点连线与像幅沿航线方向两框标连线之间的夹角称为像片的旋 偏角? 2０、倾斜误差 因像片倾斜引起的像点位移 节点 投射光线与成像光线与光轴的交角u 和u ′相等时,投射光线与成像光线与光轴的交点。 2１、投影差 ?因地形起伏引起的像点位移 2２、摄影比例尺?航摄相片上某一线段构成的长度与地面上相应水平距离之比。 ２3、像片控制点? 为联系地面与相片而测定地面坐标的像点。 相对孔径 ?物镜焦距与有效孔径之比 25、左右视差 同名像点在各自像平面坐标系中的x 坐标之差 2６、上下视差?同名像点在各自像平面坐标系中的Y 坐标之差 27、核点 基线延长线与左、右像片的交点ｋ１、k ２称为核点 ? 28、核线 核面与像片的交线称为核线 29、核面?通过摄影基线S 1S 2与任一地面点Ａ所作的平面W Ａ?? 30、投影基线?两摄站的连线? 3１、像片基线?指相邻两张像片主点的连线 ３2、解析空中三角测量 即在一条航带几十条像对覆盖的区域或由几条航带几百哥像对构 成的区域内,仅仅由外业实测几个少量的控制点，按一定的数学模型，平 差解算出摄影测量作业过程中所需的全部控制点及每张像片的外方位元 素 ()()()()(){}2332233213322232332 1[]Z X Y X Y Y Y X X X Z Y X X Y Z X Y Y X Z X Y X Y =-+-+-+-+-- 摄影测量学空间后方交会实验报告测绘101徐斌摄影测量学实验报告 实验一、单像空间后方交会 学院: 建测学院 班级: 测绘101 姓名: 徐斌 学号: 26 一( 实验目的 1.深入了解单像空间后方交会的计算过程; 2.加强空间后方交会基本公式和误差方程式，法线方程式的记忆; 3.通过上机调试程序加强动手能力的培养。 二(实验原理 以单幅影像为基础，从该影像所覆盖地面范围内若干控制点和相应点的像坐标量测值出发，根据共线条件方程，求解该影像在航空摄影时刻的相片外方位元素。 三(实验内容 1.程序图框图 2.实验数据 (1)已知航摄仪内方位元素f,153.24mm，Xo,Yo,0。限差0.1秒 (2)已知4对点的影像坐标和地面坐标: 影像坐标地面坐标 x(mm) y(mm) X(m) Y(m) Z(m) 1 -86.15 -68.99 36589.41 25273.32 2195.17 2 -53.40 82.21 37631.08 31324.51 728.69 3 -14.78 -76.63 39100.97 24934.98 2386.50 4 10.46 64.43 40426.54 30319.81 757.31 3.实验程序 Form1.cs 程序 using System; using System.Collections.Generic; using https://www.doczj.com/doc/0d10765222.html,ponentModel; using System.Data; using System.Drawing; using System.Linq; using System.Text; using System.Windows.Forms; using System.IO; namespace 后方交会1 { public partial class Form1 : Form { public Form1() { InitializeComponent(); } public double f,m, Xs, Ys, Zs, a1, a2, a3, b1, b2, b3, c1, c2, c3, q, w, k; public static int N,s; public double[] x = new double[4]; 摄影比例尺：摄影比例尺越大，像片地面的分辨率越高，有利于影像的解译与提高成图精度 摄影航高：相对航高：绝对航高： 摄影测量生产对摄影资料的基本要求：影像的色调、像片倾角（摄影机主光轴与铅垂线的夹 角，α＝ 0 时为最理想的情形）像片重叠：航向重叠：同一航线内相邻像片应有一定的影 像重叠；旁向重叠：相邻航线也应有一定的重叠；航线弯曲：一条航线内各张像片的像主点 连线不在一条直线上；像片旋角：相邻两像片的主点的连线与像片沿航线方向的两框标连线 之间的夹角；像片旋角过大会减小立体相对的有效观察范围 中心投影：所有投射线或其延长线都通过一个固定点的投影 阴位：投影中心位于物和像之间。（距摄影中心f ） 阳位：投影中心位于物和像同侧。（距摄影中心f ） 像方坐标系：像平面坐标系（像主点o 为原点） 像空间坐标系(x 、y 、-f) 像空间辅助坐标系S-uvw 物方坐标系：地面测量坐标系T-XYZ （高斯平面坐标+高程）左手系 地面摄影测量坐标系D-XYZ 内方位元素: x 0，y 0，f 作用： 1、像点的框标坐标系向像空间坐标系的改化； 2、确定摄影光束的形状； 外方位元素：确定摄影光束在摄影瞬间的空间位置和姿态的参数 线元素（X S ，Y S ，Z S ） 角元素（航向倾角?、 旁向倾角ω、 像片旋角κ） 共线条件方程（摄影中心、像点、地面点） 像点位移：因像片倾斜引起的像点位移 同摄站同主距的倾斜像片和水平像片沿等比线重 合时,地面点在倾斜像片上的像点与相应水平像片上像点之间的直线移位 像点位于等比线上，无像片倾斜引起的像点位移 等比线上部的像点的像片倾斜误差方向向着等角点 等比线下部的像点的像片倾斜误差方向背向等角点 （1） 当 时， ，即等比线上的点不会因像片倾斜产生像点位移 （2）当 ，像点位移朝向等角点(一、二像限) （3）当 ，像点位移背向等角点(三、四像限) （4）当 时，主纵线上点的位移最大 像片纠正：因像片倾斜产生的影像变形改正 因地面起伏引起的像点位移（投影差）：当地面有起伏时，高于或低于所选定的基准面 的地面点的像点，与该地面点在基准面上的垂直投影点的像点之间的直线移位 地形起伏像点位移的符号与该点的高差符号相同，像片上任何一点都存在像点位移 物镜畸变、大气折光、地球曲率及底片变形等一些因素均会导致像点位移 航摄像片：中心投影，平均比例尺，影像有变形，方位发生变化 地形图：正射投影，比例尺固定，图形形状与实地完全相似，方位保持不变 在表示方法上：地形图是按成图比例尺，用各种规定的符号、注记和等高线表示地物地 貌；航片则是通过影像的大小、形状和色调表示。 在表示内容上：在地形图上用相应的符号、文字、数字注记表示，在像片上这些是不存 ??? ????-+-+--+-+--=-+-+--+-+--=)Z Z (c )Y Y (b )X X (a )Z Z (c )Y Y (b )X X (a f y )Z Z (c )Y Y (b )X X (a )Z Z (c )Y Y (b )X X (a f x S S S S S S S S S S S S 333222 333111 { System; System.Collections.Generic; System.Linq; System.Text; namespace 单像空间后方交会 { class Program { static void Main( string [] args) for (j = 0; j < 5; j++) if (j < 3) "请输入第 {0} 个点的第 {1} 个地面坐标： ", i + 1, j + 1); double .Parse( Console .ReadLine()); "请输入第 {0} 个点的第 {1} 个像点 坐标： ", i + 1, j - 2); double .Parse( Console .ReadLine()); Console .WriteLine(); // 归算像点坐标 (i = 0; i < 4; i++) for (j = 3; j < 5; j++) if (j == 3) zuobiao[i, j] = zuobiao[i, j] - x0; else zuobiao[i, j] = zuobiao[i, j] - y0; // 计算和确定初值 double zs0 = m * f, xs0 = 0, ys0 = 0; for (i = 0; i < 4; i++) else using using using using x0 = y0 = int x0, y0, i, j; double f, m; Console .Write( " 请输入像片比例尺： "); double .Parse( Console .ReadLine()); Console .Write( " 请输入像片的内方位元素 x0:" ); // 均以毫米为单 位 int .Parse( Console .ReadLine()); Console .Write( " 请输入像片的内方位元素 y0:" ); int .Parse( Console .ReadLine()); Console .Write( " 请输入摄影机主距 f:" ); double .Parse( Console .ReadLine()); Console .WriteLine(); // 输入坐标数据 double [,] zuobiao = new double [4, 5]; (i = 0; i < 4; i++) for Console .Write( zuobiao[i, j] = Console .Write( zuobiao[i, j] = for 一． 实验目的： 掌握摄影测量空间后方交会的原理，利用计算机编程语言实现空间后方交会外方位元素的解算。 二． 仪器用具及已知数据文件： 计算机windows xp 系统，编程软件（VISUAL C++6.0），地面控制点在摄影测量坐标系中的坐标及其像点坐标文件shuju.txt 。 三． 实验内容： 单张影像的空间后方交会：利用已知地面控制点数据及相应像点坐标根据共线方程反求影像的外方位元素。 数学模型：共线条件方程式： )(3)(3)(3)(1)(1)(1Zs Z c Ys Y b Xs X a Zs Z c Ys Y b Xs X a f x -+-+--+-+--= )(3)(3)(3)(2)(2)(2Zs Z c Ys Y b Xs X a Zs Z c Ys Y b Xs X a f y -+-+--+-+--= 求解过程： （1）获取已知数据。从航摄资料中查取平均航高与摄影机主距；获取控制点的地面测量坐标并转换为地面摄影测量坐标。 （2）量测控制点的像点坐标并做系统改正。 （3）确定未知数的初始值。在竖直摄影且地面控制点大致分布均匀的情况下，按如下方法确定初始值，即： n X X S ∑=0,n Y Y S ∑=0,n Z mf Z S ∑=0 φ =ω=κ=0 式中；m 为摄影比例尺分母；n 为控制点个数。 （4）用三个角元素的初始值，计算个方向余弦，组成旋转矩阵R 。 （5）逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地面 坐标代入共线方程式，逐点计算像点坐标的近似值（x ）、（y ）。 （6）逐点计算误差方程式的系数和常数项，组成误差方程式。 （7）计算法方程的系数矩阵A A T 和常数项l A T ，组成法方程式。 （8）解法方程，求得外方位元素的改正数dXs ,S dY ，s dZ ，d φ,d ω,d κ。 （9）用前次迭代取得的近似值，加本次迭代的改正数，计算外方位元素 的新值。 摄影测量学复习资料 (全) 一、名词解释 1、解析相对定向：根据同名光线对对相交这一立体相对内在的几何关系，通过量测的像点坐标，用解析计算方法解求相对定向元素，建立与地面相似的立体模型，确定模型点的三维坐标。 2、GPS辅助空中三角测量：将基于载波相位观测量的动态 GPS 定位技术获取的摄影中心曝光时刻的三维坐标作为带权观测值，引入光束法区域网平差中，整体求解影像外方位元素和加密点的地面坐标，并对其质量进行评定的理论和方法。 3、主合点：地面上一组平行于摄影方向线的光束在像片上的构像 4、核线：立体像对中，同名光线与摄影基线所组成核面与左右像片的交线。 5、航向重叠：同一条航线上相邻两张像片的重叠度。 6、旁向重叠：两相邻航带摄区之间的重叠。 7、影像匹配：利用互相关函数，评价两块影像的相似性以确定同名点 8、影像的内方元素：是描述摄影中心与像片之间相关位置的参数。 9、影像的外方元素：描述像片在物方坐标的位置和姿态的参数。 10、景深：远景与近景之间的纵深距离称为景深 11、空间前方交会：由立体像对中两张像片的内、外方位元素和像点坐标来确定相应地面点的地面坐标的方法，称为空间前方交会。 12、空间后方交会：利用一定数量的地面控制点，根据共线条件方程或反求像片的外方位元素这种方法称为单张像片的空间后方交会。 13、摄影基线：相邻两摄站点之间的连线。 14、像主点：像片主光轴与像平面的交点。 15、立体像对：相邻摄站获取的具有一定重叠度的两张影像。 16、数字影像重采样：当欲知不位于采样点上的像素值时，需进行灰度重采样。 17、核面：过摄影基线与物方任意一点组成的平面。 18、中心投影：所有投影光线均经过同一个投影中心。 19、单模型绝对定向：相对定向所构建的立体模型经平移、缩放、旋转后纳入到地面坐标系中的过程相对定向：根据立体像对内在的几何关系恢复两张像片之间的相对位置和姿态，使同名光线对对相交，建立与地面相似的立体模型。即确定一个立体像对两像片的相对位置。 20、数字影像内定向：同一像点的像平面坐标与其扫描坐标不相等，需要加以换算，这种换算称为数字影像内定向。 21、像主点：摄影机主光轴在框标平面上的垂足 22、内部可靠性：一定假设条件下，平差系统所能发现的模型误差的下界值 22、外部可靠性：一定显著性水平和检验功效下，平差系统不能发现的模型误差对平差结果的影响。 23、摄影学：利用光学摄影机摄取相片，通过相片来研究和确定被摄物体的形状，大小，位置和相互关系的一门学科技术。 24、影像信息学：是一门记录、储存、传输、量测、处理、解译、分析和显示由非接触传感器影响获得的目标及其环境信息的科学技术和经济实体。单像空间后方交会和双像解析空间后方-前方交会的算法程序实现
摄影测量学考试知识点汇总
摄影测量学空间后方交会实验报告测绘101徐斌
摄影测量学__考前知识点整理
摄影测量学单像空间后方交会程序设计作业
空间后方交会程序
摄影测量学复习资料(全)复习过程
相关主题
文本预览