vb或MATLAB编写摄影测量前方交会-后方交会，可打支付宝

先后方交会求外方位元素

Private Sub Command1_Click()

   Dim x(4) As Double, y(4) As Double

   Dim a(8, 6) As Double

   Dim Xd(4) As Double, Yd(4) As Double, Zd(4) As Double

   Dim n As Integer, i As Integer

   Dim Xy(4) As Double, Yy(4) As Double, x_x(4) As Double, y_y(4) As Double, z_z(4) As Double

   Dim Zs0 As Double, Xs0 As Double, Ys0 As Double, ψ As Double, ω As Double, κ As Double, f As Double, dx(6, 1) As Double

   Dim L(8, 1) As Double

   Dim a1 As Double, a2 As Double, a3 As Double, b1 As Double, b2 As Double, b3 As Double, c1 As Double, c2 As Double, c3 As Double

       Dim v As Variant

   ψ = 0

   ω = 0

   κ = 0

   n = 0

   f = 153.24

   Xy(1) = Text5: Xy(2) = Text7: Xy(3) = Text9: Xy(4) = Text11

   Yy(1) = Text6: Yy(2) = Text8: Yy(3) = Text10: Yy(4) = Text12

   Xd(1) = Text13: Xd(2) = Text16: Xd(3) = Text19: Xd(4) = Text22

   Yd(1) = Text14: Yd(2) = Text17: Yd(3) = Text20: Yd(4) = Text23

   Zd(1) = Text15: Zd(2) = Text18: Zd(3) = Text21: Zd(4) = Text24

   Zs0 = 2000

   Xs0 = (Xd(1) + Xd(2) + Xd(3) + Xd(4)) / 4

   Ys0 = (Yd(1) + Yd(2) + Yd(3) + Yd(4)) / 4

   Do

   a1 = Cos(ψ) * Cos(κ) - Sin(ψ) * Sin(ω) * Sin(κ)

   a2 = -Cos(ψ) * Sin(κ) - Sin(ψ) * Sin(ω) * Cos(κ)

   a3 = -Sin(ψ) * Cos(ω)

   b1 = Cos(ω) * Sin(κ)

   b2 = Cos(ω) * Cos(κ)

   b3 = -Sin(ω)

   c1 = Sin(ψ) * Cos(κ) + Cos(ψ) * Sin(ω) * Sin(κ)

   c2 = -Sin(ψ) * Sin(κ) + Cos(ψ) * Sin(ω) * Cos(κ)

   c3 = Cos(ω) * Cos(ψ)

   For i = 1 To 4

     x_x(i) = a1 * (Xd(i) - Xs0) + b1 * (Yd(i) - Ys0) + c1 * (Zd(i) - Zs0)

     y_y(i) = a2 * (Xd(i) - Xs0) + b2 * (Yd(i) - Ys0) + c2 * (Zd(i) - Zs0)

     z_z(i) = a3 * (Xd(i) - Xs0) + b3 * (Yd(i) - Ys0) + c3 * (Zd(i) - Zs0)

     x(i) = -f * x_x(i) / z_z(i)

     y(i) = -f * y_y(i) / z_z(i)

     L(2 * i - 1, 1) = Xy(i) - x(i)

     L(2 * i, 1) = Yy(i) - y(i)

     a(2 * i - 1, 1) = (a1 * f + a3 * Xy(i)) / z_z(i)

     a(2 * i - 1, 2) = (b1 * f + b3 * Xy(i)) / z_z(i)

     a(2 * i - 1, 3) = (c1 * f + c3 * Xy(i)) / z_z(i)

     a(2 * i - 1, 4) = Yy(i) * Sin(ω) - ((Xy(i) / f) * (Xy(i) * Cos(κ) - Yy(i) * Sin(κ)) 衫消+ f * Cos(κ)) * Cos(ω)

     a(2 * i - 1, 5) = -f * Sin(κ) - (Xy(i) / f) * (Xy(i) * Sin(κ) + Yy(i) * Cos(κ))

     a(2 * 扰历i - 1, 6) = Yy(i)

     a(2 * i, 1) = (a2 * f + a3 * Yy(i)) / z_z(i)

     a(2 * i, 2) = (b2 * f + b3 * Yy(i)) / z_z(i)

     a(2 * i, 3) = (c2 * f + c3 * Yy(i)) / z_z(i)

     a(2 * i, 4) = -Xy(i) * Sin(ω) - ((Yy(i) / f) * (Xy(i) * Cos(κ) - Yy(i) * Sin(κ)) - f * Sin(κ)) * Cos(ω)

     a(2 * i, 5) = -f * 或李知Cos(κ) - Yy(i) / f * (Xy(i) * Sin(κ) + Yy(i) * Cos(κ))

     a(2 * i, 6) = -Xy(i)

   Next i

Call EquationResult(a, L, dx)

  ψ = ψ + dx(4, 1): ω = ω + dx(5, 1): κ = κ + dx(6, 1)

  Xs0 = Xs0 + dx(1, 1): Ys0 = Ys0 + dx(2, 1): Zs0 = Zs0 + dx(3, 1)

  Text25 = Xs0: Text26 = Ys0: Text27 = Zs0

  n = n + 1

Loop While Abs(dx(4, 1)) > 0.00003 Or Abs(dx(5, 1)) > 0.00003 Or Abs(dx(6, 1)) > 0.00003 Or Abs(dx(1, 1)) > 0.001 Or Abs(dx(2, 1)) > 0.001 Or Abs(dx(3, 1)) > 0.001

MsgBox "迭代次数为" & n & "次"

End Sub

在前方交会求待定点的坐标

2017摄影测量实习报告【一】

一、实习目的与要求

本次实习是在摄影测量的教学基础上，理论实际相联系的动手 *** 作实习，是我们在学习测量专业的一个重要的实习环节。一方面是培养我们的实践 *** 作能力和运用软件解算数据的能力，另一方面培养我们在今后遇到问题应该如何去解决的能力，通过实习发现自己在实践动手方面的不足并想办法解决，为以后的工作实践打下扎实的基础。使我们熟练地掌握摄影测量及遥感的原理，信息获取的途径，数字处理系统和应用处理方法。并进一步巩固和深化理论知识，使理论与实践相结合。切实加强我们大家的实践动手能力，提高大家对这门新技术的认识和把握，全面培养我们的应用能了、创新能力和探索精神。

二、实习地点

桂林市雁山区大埠乡

桂林理工大学博文管理学院机房

三、实习用具

小比例尺航片两张、画图板一个、透明纸两张、铅笔、橡皮电子计算机、ENVI遥感图像处理系统、编程软件(MATLAB、Visual Basic)

四、实习任务与要求

掌握航片调绘的方法步骤

掌握使用编程软件设计解算移动曲面法数字高程模型内插子程序

掌握使用编程软件设计解算空间后方交会

掌指余握使用ENVI遥感图像处理系统处理遥感影像

五、实习步骤

航片调绘

1

本次实习的遥感图像调丛逗旁绘主要判读航片测区地物属性，在透明纸上勾出边界，必要时进行清绘。在进行野外调绘之前，将调绘航片平放在画图板上，然后再将比调绘图稍大一些的透明纸盖于调绘航片上，用胶带粘好，连同调绘航片用夹子固定于画图板。

方程的线性化形式，计算近似垂直摄影情况下像片的外方位元素。

解算步骤：

获取已知数据：m,x0,y0,f,Xtp,Ytp,Ztp

量测控制点像点坐标：x，y

确定未知数初值

组成误差方程式：若P=I，X=(ATA)?TL

解求外方位元素改正数、外方位元素的近似值

检查迭代是否收敛，是否需要重复计算。

使用ENVI系统处理遥感影像

主要要求学会使用ENVI系统对遥感影像进行监督分类和非监督分类

监督分类

制作分类模版：打开一幅遥感影像，在影像窗口打开Overlay-Region of interest，在Zoom窗口依次绘制可识别地物类别的区域

监督分类：Classification-Supervised-Minimum Distance

渗橡监督分类后处理：Classification-post Classification——Clump classes。

非监督分类

打开一幅遥感影像，单击主菜单Classification-Unsupervised-Isodata，得到非监督分类的结果

点击Classification-post Classification-Combine classes合并相同或相似类别进行监督分类后处理

六、实习心得

摄影测量是一门专业的测绘学科，也是一门应用很广的学科，随着遥感技术的不断发展，这门学科正从几何学向信息科学发展。它的发展及运用对我们测绘来说是很有帮助的。而摄影测量实习则可以提高我们对摄影测量知识的理解，加强我们的实际运用能力。因此学校安排了三周的摄影测量实习，这对提升我们的摄影测实际 *** 作能力是很有帮助的。

此次实习分四个板块，分别是全数字立体测图数字摄影测量的编程遥感影像自动分类像片的判读与调绘等。全数字立体测图是利用计算机代替解析测图仪、用数字影像代替模拟像片、用数字光标代替光学光标，直接在计算机上进行数字化测图的作业方法。这个实习要求我们学会使用ENVY软件构建立体模型，制作测区的DEM、DOM和等高线图，同时熟练使用交互式数字影像测图系统在立体影像上量测不同类地物， 并时行地物数据采集及我，生成数字测图文件，按标准的制图符号将之输出为矢量地形图。数字摄影测量的编程则要求我们学习使用Matlab进行摄影测量编程，掌握移动曲面法数字高程模型内插子程序的设计方法和空间后方交会程序的设计方法。遥感影像自动分类则是让我们了解并掌握督与非督分类的过程和方法，并利用监督分类结果制作一幅影像地图。像片的判读与调绘则是让我们利用学过的几类常用遥感影像的判读技术与方法，完成航空像片或彩红外片的判读和外业调绘工作，掌握全野外调绘的基本技能。

2017摄影测量实习报告【二】

一、 实习任务

利用自己所熟悉的一种编程语言，实现单像空间后方交会，解求此张像片的6个外方位元素 ， ， ， ，ω，κ ，范文之实习报告:摄影测量实习报告。

二、 实习目的

1、 深刻理解单张像片空间后方交会的原理与意义

2、 在存在多余观测值时，利用最小二乘平差方法，经过迭代，求的外方位元素的最佳值

3、 熟悉VC编程方法，利用编程实现计算。

三、 实习原理

以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，求解该影象在航空摄影时刻的像片外方位元素 ， ， ， ，ω，κ共线条件方程如下：

x-x0=-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为影像的外方位元素

， ，为摄站点的物方空间坐标X,Y,Z为物方点的物方空间坐标

旋转矩阵R为

由于此共线条件方程是非线性方程，先对其进行线性化，利用泰勒展开得：

=(x)-x++++++++

=(y)-y++++++++

像点观测值一般视为等权，即P=I

矩阵形式：V=AX-L，P=I

通过间接平差，为提高精度，增加多余观测方程，根据最小二乘平差原理，可计算出外方位元素的改正数。经过迭代计算，每次迭代用未知数的近似值与上次迭代计算的改正数之和作为新的近似值，重复计算，求出新的改正数，这样反复趋近，直到改正数小于某个限值为止。

四、 程序框图

输入原始数据

归算像点坐标x，y

计算并确定初值 ， ， ， ，

组成旋转矩阵R

计算(x)(y)和

逐点组成误差方程式并法化

所有点完否?

解法方程，求未知数改正数

计算改正后的外方位元素

未知数改正数<限差否?

整理并输出计算结果

正常结束

非正常结束

输出中间结果和出错信息

迭代次数是否小于限差否?

否

否

否

是

五、计算结果

1、像点坐标，地面坐标

点数

像点编号 x y X Y Z

2像片内方位元素：f = 153.840 x0=y0=0

摄影比例尺：1：2500

运算结果:

六、 数据分析

选取第六张像片进行计算，迭代次数为2次。经过比较发现，计算出的6个外方位元素与所给参考值相比，相差很小，计算结果符合要求：线元素误差小于0.5米角元素误差30秒。

计算其精度，可以通过法方程式中未知数的系数矩阵的逆阵(A)-1来求解，此时，视像点坐标为等精度不相关观测值。因为逆阵中第i个主对角线上元素Qii就是法方程式中第i个未知数的权倒数，若单位权中误差为m0，则第i个未知数的中误差为：

mi=

摄影测量实习报告4篇汇集摄影测量实习报告4篇汇集当参加空间后方交会的控制点有n个时，则单位权中误差可按下式计算：

m0=

要求：线元素精度mx等，高于0.05米角元素精度高于0.00003弧度。计算结果都达到标准。

在此次计算中，我运用了所给的全部控制点，而空间后方交会所运用的控制点，应该避免位于一个圆柱面上，否则会出现解不唯一的情况。选点时，还需要避免选择的点过于聚集在一起，或位于一条直线上，所选控制点最好分布在像片的四角和中央。并且数量充足，这样有利于提高解算精度。

迭代时，所选择控制条件不同，迭代次数略有不同，所以最后结果也会略有不同。一般设置为线元素改正数小于0.01m，角元素改正数小于0.1’。

所提供X Y Z为地面测量坐标，带入共线方程时，需要转换为地面摄影测量坐标，最简单的方法为互换XY的数值，即可达到转换坐标目的。并且其单位为米，而像点坐标的单位为厘米，需要统一坐标单位。

这次实习持续时间很长，经历了几次数据的更改，所以程序也几经修改，由最初的直接输入数据到后来可以以自行读入数据，并且可以选择计算的像片，功能有所完善，我也在实践的过程中，对空间后方交会有了更深的理解。深刻理解了共线条件方程的运用，各个量的意义，受益匪浅。

2017摄影测量实习报告【三】

一.实习目的：

1.联系水准仪的安置、整平、瞄准与读书和怎样测定地面两点间的高程

2.掌握经纬仪对中，整平，瞄准与读书等基本 *** 作要领

3.掌握钢尺量距的一般方法

4.练习用经纬仪配合小平板仪测绘地形图

5.培养学生综合应用测量理论知识分析解决土建施工放样中一般问题的能#from 本文来自高考资源网http://www.gkstk.com end#力。

二.实习要求：

1.掌握水准仪、经纬仪、平板仪等一些主要一起的性能和如何 *** 作使用

2.掌握数据的计算和处理方法

3.掌握地形图测绘的基本方法，具有初步测绘小区域大比例尺地形图的工作能力

4.掌握测设的基本方法。

三.实习任务

1.控制点高程测量

2.导线长度测量

3.水平角度测量

摄影测量实习报告4篇汇集文章摄影测量实习报告4篇汇集出自http://www.gkstk.com/article/1428157571170.html，转载请保留此链接!4.闭合导线内业计算

5.测绘地形图

6.建筑物放样，道路中线测量和管线测量。

四.实习项目与程序：

1.外业测量：

(1)测量控制点高程

(2)测量控制点间距离

(3)测量闭合导线内角。

2.内业计算：

(1)计算控制点间高差，推算各点间高程

(2)计算个控制点间距离及相对误差

(3)计算个内角闭合差及内角

(4)根据以上计算数据推算个点坐标。

3.测绘地形图：

(1)将坐标范围内的控制点标定到图纸上

(2)根据控制点周围的地物地貌测量某些点的高程，再标在图纸上。

4.施工测量：圆曲线，缓和曲线测量。

五.实习方法：水准仪、经纬仪、平板仪的使用。

1.测量的方法：

(1)水准仪架在两个控制点的中间，距离两点大致相等。在前后两点各立水准尺一把。

(2)望远镜对准水准尺并推动，再将水准仪调平，调节三个脚螺旋，使得圆水准器旗袍居中，然后微调倾螺旋，从左边的窗口看到水准管的气泡闭合。

(3)调水平微动螺旋，使得十子丝在水准尺上测得后视读数和前视读数并记录下来。

(4)三脚架架腿抬高或降低，重新测量后视读数和前视读数并记录下来，测得高差不得超过5mm，否则重测。

2.角度测量的方法：

(1)经纬仪架在控制点上，用脚螺旋进行对中，再伸缩架腿调节圆水准气泡居中，然后调节脚螺旋使得水准管气泡也居中。通过对中器观察是否对中，否则反复调平。

(2)望远镜调成盘左，对准左面的目标并制动，调节微倾和微动螺旋，使得十字丝瞄准目标，把配置度盘的按钮拔出，记下读数。顺时针转动照准部，对准右面的目标并制动，读出右面的读数，记录读数。

(3)望远镜调成盘右，对准右面的目标并制动，调节调节微倾和微动螺旋，是的十字丝瞄准目标，把配置度盘的按钮拔出，记下读数。逆时针转动照准部，对准左面的目标并制动，读出左边的读数，记录读数。

(4)两次测量角之差不能超过40秒，否则重测。

3.距离测量的方法：

(1)用前面的方法将经纬仪对中整平，再进行定线。

(2)然后用钢尺沿着路线测出导线长度。

(3)往返各测一次，两次距离的相对误差不能超过三千分之一，否则重测。

4.平板仪的使用：

(1)立上三脚架，将平板固定，把图纸也固定在平板上。

(2)将平板仪的一边靠在两个控制点上，瞄准地面上的点，然后进行对中整平。

(3)整平后进行测绘。量出控制点到某地物的距离并且紧靠建筑物立标杆，通过平板仪瞄准标杆则确定了这个方向。根据比例尺换算成图上距离，将地物地貌画在图上。

(4)将所有坐标范围内的地物地貌全都画在图上，并用规定符号表示。

六.测量精度：

1.距离往返测量相对误差不超过1/3000

2.水准仪高差测量中高差闭合差在容许值±12vnmm或±40Lmm范围内

3.测内角时一测回中上、下半测回角值之差不得超过±40``。

七.计算成果和示意图见测设数据计算表。

八.体会：

通过本次实习，巩固了以前所学知识，掌握了水准仪、经纬仪的基本 *** 作，还有怎样施工放样，学会了地形图的绘制方法。从而积累了许多经验。(1)、立标尺时，标尺除立直外还要选在重要的地方.因此，选点就非常重要，点一定要选在有代表性的地方.同时要注意并点非越多越好.相反选取的无用点过多不但会增加测量，计算和绘图的劳动量和多费时间，而且会因点多而杂乱产生较大的误差。(2)、要先将道路和主要建筑物确定下来，然后在添加其余次要方面，这样不但条理清楚，利于作图而且更有利于作图的准确和随时进行实物和图形的对比从而检验测量数据的准确与否.但这两周实习也给了我们不少教训：由于某个数据的读错、记错及算错都给我们带来了不少麻烦，从而让我们知道了做任何事都要认真。还有一个组的团结也是至关重要的，他关系到整个组的进度。先前我们组由于配合不够默契，分工也不够合理，整体进度受到极大的影响，后来通过组内的交流，彻底解决了以上问题。实习进度有了很大的改观，进度和效果自然就提上来了。这告诉我们团结就是力量，我们以后工作的时候也是一样，只有团结才能把事办好。

---