前方交会实验报告(含VB程序代码)

需要源程序设计的的可以联系我~

立体像对前方交会实验报告

一、实验目的

在掌握前方交会原理的基础上，自己编写前方交会程序，在计算机上调试，输出计算结果并对计算结果进行检验。通过上机调试程序加强动手能力的培养，通过对实验过程的掌握以及对实验结果的分析，增强学生综合运用所学知识解决实际问题的能力。

二、实验仪器

计算机，VB6.0

三、实验数据

1. 模拟像片一对：左片号23，右片号24； 2. 航摄机主距：f=150mm； 3. 左片23号片外方位元素：

= 0°25′00″ = 1°00′00″ = 0°10′00″ Xs=103007.006117 Ys=139998.994849 Zs=4801.9989994 （m） 右片24号片外方为元素： =1°39′59″ = 0°10′00″ =0°40′00″ Xs=106002.023762 Ys=140005.002780 Zs=4797.009648 （m）

待求像点坐标如下表：

四、实验内容

利用所给立体像对两张像片的内、外方位元素，编写空间前方交会程序，根据所给像对中若干同名像点在左右像片上的坐标，解求其对应的地面点的物方坐标，实现空间前方交会的过程。 五、实验成果

程序流程图：

需要源程序设计的的可以联系我~

程序设计界面：

需要源程序设计的的可以联系我~

程序运行界面：

需要源程序设计的的可以联系我~

运行结果：

（注：表上显示地面点坐标依次是：7,9,4,6,5,8）

附：excel进行角度转换：

六、程序如下：

Dim f#, x1#, y1#, x2#, y2#, i%, j%, u1#, u2#, v1#, v2#, w1#, w2#, fai1#, kab1#, omg1#, fai2#, kab2#, omg2# Dim a12#, a13#, b11#, b12#, b13#, c11#, c12#, c13#, a21#, a22#, a23#, b21#, b22#, b23#, c21#, c22#, c23# Dim n1#, n2#, bu#, bv#, bw#

Dim xs1#, xs2#, ys1#, ys2#, zs1#, zs2#

需要源程序设计的的可以联系我~

Private Sub Command1_Click()

fai1 = Val(Text1.Text): kab1 = Val(Text3.Text): omg1 = Val(Text2.Text) fai2 = Val(Text11.Text): kab2 = Val(Text10.Text): omg2 = Val(Text7.Text) f = Val(Text13.Text)

xs1 = Val(Text4.Text): xs2 = Val(Text9.Text): ys1 = Val(Text5.Text): ys2 = Val(Text8.Text): zs1 = Val(Text6.Text): zs2 = Val(Text12.Text) End Sub

Private Sub Command3_Click()

x1 = Val(InputBox("输入23片坐标x的值")) y1 = Val(InputBox("输入23片坐标y的值")) x2 = Val(InputBox("输入24片坐标x的值")) y2 = Val(InputBox("输入24片坐标y的值")) End Sub

Private Sub Command2_Click() Text1.Visible = False Text2.Visible = False Text3.Visible = False Text4.Visible = False Text5.Visible = False Text6.Visible = False Text7.Visible = False Text8.Visible = False Text9.Visible = False Text10.Visible = False Text11.Visible = False Text12.Visible = False Text13.Visible = False Label1.Visible = False Label2.Visible = False Label3.Visible = False Label4.Visible = False Label5.Visible = False Label6.Visible = False Label7.Visible = False Label8.Visible = False Label9.Visible = False Label10.Visible = False Label11.Visible = False Label12.Visible = False Label13.Visible = False Label14.Visible = False

a11 = Cos(fai1) * Cos(kab1) - Sin(fai1) * Sin(omg1) * Sin(kab1) a12 = -Cos(fai1) * Sin(kab1) - Sin(fai1) * Sin(omg1) * Cos(kab1) a13 = -Sin(fai1) * Cos(omg1) b11 = Cos(omg1) * Sin(kab1)

需要源程序设计的的可以联系我~

b12 = Cos(omg1) * Cos(kab1) b13 = -Sin(omg1)

c11 = Sin(fai1) * Cos(kab1) + Cos(fai1) * Sin(omg1) * Sin(kab1) c12 = -Sin(fai1) * Sin(kab1) + Cos(fai1) * Sin(omg1) * Cos(kab1) c13 = Cos(fai1) * Cos(omg1)

a21 = Cos(fai2) * Cos(kab2) - Sin(fai2) * Sin(omg2) * Sin(kab2) a22 = -Cos(fai2) * Sin(kab2) - Sin(fai2) * Sin(omg2) * Cos(kab2) a23 = -Sin(fai2) * Cos(omg2) b21 = Cos(omg2) * Sin(kab2) b22 = Cos(omg2) * Cos(kab2) b23 = -Sin(omg2)

c21 = Sin(fai2) * Cos(kab2) + Cos(fai2) * Sin(omg2) * Sin(kab2) c22 = -Sin(fai2) * Sin(kab2) + Cos(fai2) * Sin(omg2) * Cos(kab2) c23 = Cos(fai2) * Cos(omg2)

u1 = a11 * x1 + a12 * y1 - a13 * f v1 = b11 * x1 + b12 * y1 - b13 * f w1 = c11 * x1 + c12 * y1 - c13 * f

u2 = a21 * x2 + a22 * y2 - a23 * f v2 = b21 * x2 + b22 * y2 - b23 * f w2 = c21 * x2 + c22 * y2 - c23 * f

bu = xs2 - xs1 bv = ys2 - ys1 bw = zs2 - zs1

n1 = (bu * w2 - bw * u2) / (u1 * w2 - u2 * w1) n2 = (bu * w1 - bw * u1) / (u1 * w2 - u2 * w1)

Dim xx#, yy#, yy1#, yy2#, zz#

xx = Fix((xs1 + u1 * n1) * 1000 + 0.5) / 1000 yy1 = Val(Text5.Text) + v1 * n1 yy2 = Val(Text8.Text) + v2 * n2

yy = Fix(((yy1 + yy2) / 2) * 1000 + 0.5) / 1000

zz = Fix((Val(Text6.Text) + w1 * n1) * 1000 + 0.5) / 1000 Form1.Print "地面坐标为:"

Form1.Print "xx="; xx, "yy="; yy; "zz="; zz End Sub

七、实验心得

通过本次前方交会实验，熟悉掌握了前方交会原理，并利用计算机程序对其进行了解决，收益颇多！