测量专业基本计算与编程(3)

Text1.text=\：Text2.text=\：Text3.text=\给文本框赋初值 Text1.SetFocus' Text1文本框获得焦点 End sub

然后运行，在Text1文本框与Text2文本框中输入数据，单击计算按钮，在Text3文本框中显示计算结果。

一个较好的应用程序必须要有好的用户界面，在大多应用程序中大多包含着许多通用的东西，如工具栏、状态栏、工具提示、上下文菜单、帮助以及选项卡对话框等。VB 具有把所有这些添加到应用程序中的能力，因此应该很好的学习并加以掌握运用。

二 测量中几种常用的计算 角度化为弧度、弧度化为角度以及坐标方位角的计算在测量数据处理中是经常遇到的计算问题。也是最基本的计算问题。

2.1 角度化为弧度

将度、分、秒形式的角度angle化为弧度，采用函数功能来实现，其函数值为返回的弧度值。其数据格式为: 参数angle的整数为表示度，小数点后两位表示分, 小数点后第三位开始表示秒， 如180度34分54.23秒, 即 angle=180.345423。 Function Radian(ByVal angle As Double) As Double Dim mm As Double, a As Double a = Abs(xx)

a = a + 0.0000001 dd% = Int(a)

ii% = Int((a - dd%) * 100) mm = (a - dd%) * 100 - ii% mm = dd% + ii% / 60 + mm / 36 Radian = pi * mm / 180

Radian = Sgn(angle) * Radian End Function

2.2 弧度化为角度计算

将弧度值radian 化为度、分、秒的角度形式，采用函数功能实现，计算返回角度值。其数据格式为: 函数值qdms的整数为表示度，小数点后两位表示分, 小数点后第三位开始表示秒， 如160度14分23.03秒, 即qdms=160.142303。 Function qdms(ByVal radian As Double) As Double Dim a1, xx , second As Double a1 = Abs(radian) a1 = a1 * 180 / pi degree% = Int(a1)

xx = (a1 - degree%) * 60 minute% = Int(xx)

second = (xx - minute%) * 60

qdms = degree% + minute% / 100 + second / 10000 qdms = Sgn(radian) * qdms End Function

2.3 坐标方位的计算

已知点号为 和 两点的平面坐标 、 、 、 ，求 和 两点的坐标方位角, 先计算 ，在按如下计算公式如下：

Function qiua(ByVal px As Double, ByVal py As Double) As Double Dim PZ As Double If px = 0 Then PZ = pi / 2#

If py < 0 Then PZ = pi * 1.5 Else

PZ = Atn(py / px)

If px < 0 Then PZ = PZ + pi

If py < 0 And px > 0 Then PZ = PZ + pi * 2 End If qiua = PZ End Function

三 空间大地坐标与空间直角坐标之间的换算 3.1 由空间大地坐标计算空间直角坐标

已知椭球 、 ，大地坐标直角坐标X、Y、Z，求大地坐标B、L、H，计算公式如下 式中 。

3.2 空间直角坐标计算空间大地坐标

已知椭球 、 ，大地坐标直角坐标X、Y、Z，求大地坐标B、L、H。 1、 直接解法计算公式（参见熊介编著：《椭球大地测量学》） 计算辅助量 计算大地经度 计算大地纬度

计算大地高

2、 迭代解法计算公式汇编 计算辅助量

， ， 计算大地纬度

计算大地纬度

令 ， ，

迭代循环计算 直到满足 ，以保证 的计算精度至 计算大地高

3.3 计算程序与范例

空间大地坐标与空间直角坐标之间的计算程序界面如下图所示。在窗体XYZ-BLH上设置3个Frame控件，在 Frame1 内放置3个选择按钮 option1(0)、option1(1)、option(2)，控制椭球参数；在Frame2 内放置2个选择按钮 option2(0) 与 option2(1)，控制计算方法；在Frame3 内放置1个MSFlexGrid 控件，4个 CommandButton 命令按钮，其 Caption 属性分别为文件、添加、删除、保存，其功能主要用来编辑数据与保存数据；在窗体 XYZ-BLH上 设置3个 CommandButton 命令按钮，其 Caption 属性分别为计算、显示、退出，其功能为计算、用 Word 显示计算结果、退出。其计算界面直观、方便。 计算范例1：

已知数据 ， ， 。

计算结果 ， ， 。 计算范例2：

已知数据 ， ， 。 计算结果 ， ， 。

图3－1 空间大地坐标与空间直角坐标换算程序界面

上述计算采用直接算法，计算子程序如下： Sub BLH_XYZ(b, l, H, x, y, Z, ByVal K As Integer) EPS = e2 / (1# - e2) BB = a / Sqr(1# + EPS) Select Case K Case 1

p = Sqr(x ^ 2 + y ^ 2) UO = Atn(Z * a / p / BB)

b = Atn((Z + EPS * BB * Sin(UO) ^ 3) / (p - e2 * a * Cos(UO) ^ 3)) p1 = Sqr(1# + (1# - e2) * Tan(b) ^ 2) b = Atn(Z / p + a * e2 * Tan(b) / p / p1) l = Atn(y / x)

If b < 0# Then b = b + 2# * pi If l < 0# Then l = l + 2# * pi U = Atn(BB / a * Tan(b))

H = Sqr((p - a * Cos(U)) ^ 2 + (Z - BB * Sin(U)) ^ 2) If (p - a * Cos(U)) < 0# Then H = -H Case 2

RN = a / Sqr(1# - e2 * Sin(b) * Sin(b))

x = (RN + H) * Cos(b) * Cos(l) y = (RN + H) * Cos(b) * Sin(l) Z = (RN * (1# - e2) + H) * Sin(b) End Select End Sub

四 大地主题问题正反计算

参考椭球面是大地测量计算的基准面。大地坐标是椭球面上的基本坐标系，根据大地测量的观测成果（如距离与方向），从大地原点出发，逐点计算在椭球面上的大地坐标；或根据两点的大地坐标，计算它们之间的大地线长度和大地方位角，这类计算称为大地问题解算，又叫大地主题解算。大地问题解算的用途是多方面的，随着现代空间技术和航空航天、航海等领域的发展，大地问题解算（尤其是大地反算）有着更为重要的作用。鉴于各种用途与要求不同，产生了不同的大地解算方法与公式。

椭球面上大地坐标的解算比平面坐标的解算要复杂得多，正是由于这种复杂性导致了大地问题解算公式的多样化，其解算方法多大几十种。若按解算的距离来分类，一般分为短距离（400公里以内），中距离（400～1000公里）和长距离（1000～2000公里）；若按解法分类可分为直接解法和间接解法；若按解算精度来分类，又可分为精密公式与近似公式等。这里介绍三种方法即高斯平均引数法、白塞尔方法与嵌套系数法的计算过程、步骤与计算程序。

4.1 高斯平均引数大地问题解算

1、高斯平均引数正算计算公式（ km） 1）计算辅助量公式

2) 计算 、 、 的初值 3）计算 、 、 …… 此处隐藏：1958字，全部文档内容请下载后查看。喜欢就下载吧 ……