用(VB)实现测量坐标转换系统.docVIP

下载本文档

11
0
约8.09千字
约 13页
2019-01-16 发布于安徽
举报
版权申诉

用(VB)实现测量坐标转换系统.doc

1、本文档共13页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

实用标准文案精彩文档坐标转换系统（VB） ????????????????????????????????????????????????????????????????????????东华理工大学 Theory 北京54高斯坐标转换西安80高斯坐标转换系统1.0版多点处理结果 ?核心源代码：高斯坐标转换成大地坐标过程Public Sub GausReverse(a As Double, f As Double, x() As Double, y() As Double, RB() As Double, RL() As Double, k As Integer)Dim i As Integer, fxb As Double, fxbl As Double, fybl As Doublee = Sqr(2 * f - f ^ 2)C = a / Sqr(1 - e ^ 2)e2 = e / Sqr(1 - e ^ 2)beita0 = 1 - (3 / 4) * e2 ^ 2 + (45 / 64) * e2 ^ 4 - (175 / 256) * e2 ^ 6 + (11025 / 16384) * e2 ^ 8beita2 = beita0 - 1beita4 = (15 / 32) * e2 ^ 4 - (175 / 384) * e2 ^ 6 + (3675 / 8192) * e2 ^ 8beita6 = (-35 / 96) * e2 ^ 6 + (735 / 2048) * e2 ^ 8beita8 = (315 / 1024) * e2 ^ 8For i = 1 To k??? B0 = x(i) / (C * beita0)??? a1 = a * Cos(B0) / Sqr(1 - e ^ 2 * (Sin(B0)) ^ 2)??? l0 = y(i) / a1 Do??? fxb = 0??? fxbl = 0??? fybl = 0??? t = Tan(B0)??? yita = e2 * Cos(B0)??? n = a / Sqr(1 - e ^ 2 * (Sin(B0)) ^ 2)??? a2 = (1 / 2) * n * Sin(B0) * Cos(B0)??? a3 = (1 / 6) * n * (Cos(B0)) ^ 3 * (1 - t ^ 2 + yita ^ 2)??? a4 = (1 / 24) * n * Sin(B0) * (Cos(B0)) ^ 3 * (5 - t ^ 2 + 9 * yita ^ 2 + 4 * yita ^ 4)??? a5 = (1 / 120) * n * (Cos(B0)) ^ 5 * (5 - 18 * t ^ 2 + t ^ 4 + 14 * yita ^ 2 - 58 * yita ^ 2 * t ^ 2)??? a6 = (1 / 720) * n * Sin(B0) * (Cos(B0)) ^ 5 * (61 - 58 * t ^ 2 + t ^ 4)??? fxb = fxb + (C * beita6 + C * beita8 * (Cos(B0)) ^ 2) * (Cos(B0)) ^ 2??? fxb = (fxb + C * beita4) * (Cos(B0)) ^ 2??? fxb = (fxb + C * beita2) * Sin(B0) * Cos(B0)??? fxbl = a2 * l0 ^ 2 + a4 * l0 ^ 4 + a6 * l0 ^ 6??? fybl = a3 * l0 ^ 3 + a5 * l0 ^ 5??? RB(i) = (x(i) - fxb - fxbl) / (C * beita0)??? a1 = (a * Cos(RB(i))) / Sqr(1 - e ^ 2 * (Sin(RB(i)) ^ 2))??? RL(i) = (y(i) - fybl) / a1??? If Abs(RB(i) - B0) = 0.0000000001 And Abs(RL(i) - l0) = 0.0000000001 Then??????? RL(i) = zrl + l0??????? Exit Do??? Else??????? B0 = RB(i)??????? l0 = RL(i)??? End IfLoopNext iEnd Sub 大地坐标B，L转换为高斯坐标x,y的过程Public Sub BLHGaus(RB() As Double, RL() As Double, GX