《用实现测量坐标转换系统》由会员分享，可在线阅读，更多相关《用实现测量坐标转换系统（13页珍藏版）》请在金锄头文库上搜索。
1、坐标转换系统（VB ）东华理工大学Theory北京54高斯坐标转换西安80高斯坐标转换系统1.0版-匸Fj-Si-Ifi冋需点赴理打幵处共点戲爵njFwmunm厂UF-一MhinBOJS点处理* J吒址标系I i *a0D坐棣系高聊坐标k高斷坐标F高厮生际y酬h1 1桃h多点处理结果二甲央子午卿经度R尹幵弐址空矗|町-54坐様点导Xi-yTH0413762520 59409497513 333MG403TB2M9 14123媳7760 964館34-换3TT822 127IB49fcT30 43803343TT9033 3791996706 3191434G37377C236.35930495
2、$S3 0451733.鬲3778505 76197496250 39T0933d3117197 8295987 93TC34537T77S7 35212136603 96T393G333T7B420.6632149M5S 7232334.323776392 1945697141, J24SB34的1377513 fC&t帕7蛙7羽嗨I34303T7B09O 110004 歯KE $312134377356. &1256498753.50420T 1专t曽qgi rffvnncrrw? 8仲BJ-54 加 i inSO 草点 it理监20换系核心源代码:高斯坐标转换成大地坐标过程Public
3、Sub GausReverse(a As Double, f As Double, x() As Double, y() As Double, RB() AsDouble, RL() As Double, k As In teger)Dim i As In teger, fxb As Double, fxbl As Double, fybl As Doublee = Sqr(2 * f - f A 2)C = a / Sqr(1 - e a 2)e2 = e / Sqr(1 - e a 2)beitaO = 1 - (3 / 4) * e2 a 2 + (45 / 64) * e2 a 4 -
4、 (175 / 256) * e2 a 6 + (11025 / 16384) * e2 a 8beita2 = beitaO - 1beita4 = (15 / 32) * e2 a 4 - (175 / 384) * e2 a 6 + (3675 / 8192) * e2 a 8beita6 = (-35 / 96) * e2 a 6 + (735 / 2048) * e2 a 8beita8 = (315 / 1024) * e2 a 8Fori = 1 To kBO = x(i) / (C * beitaO)a1 = a * Cos(BO) / Sqr(1 - e a 2 * (Si
5、n(B0) a 2)lO = y(i) / a1Dofxb = 0fxbl = 0fybl = 0t = Tan (BO)yita = e2 * Cos(BO)n = a / Sqr(1 - e a 2 * (Si n(B0) a 2)a2 = (1 / 2) * n * Sin(B0) * Cos(BO)a3 = (1 / 6) * n * (Cos(BO) a 3 * (1 - t a 2 + yita a 2)a4 = (1 / 24) * n * Sin (BO) * (Cos(BO) a 3 * (5 - t a 2 + 9 * yita a 2 +4*yita a 4)a5 = (
6、1 / 120) * n * (Cos(BO) a 5 * (5 - 18 * t a 2 + t a 4 + 14 * yita a2-58 *yita a 2 * t a 2)a6 = (1 / 720) * n * Sin (BO) * (Cos(BO) a 5 * (61 - 58 * t a 2 + t a 4)fxb = fxb + (C * beita6 + C * beita8 * (Cos(BO) a 2) * (Cos(BO) a 2fxb = (fxb + C * beita4) * (Cos(BO) A 2fxb = (fxb + C * beita2) * Sin (
7、BO) * Cos(B0)fxbl = a2 * 10 a 2 + a4 * 10 a 4 + a6 * 10 a 6fybl = a3 * 10 a 3 + a5 * 10 a 5RB(i) = (x(i) - fxb - fxbl) / (C * beita0)al = (a * Cos(RB(i) / Sqr(1 - e a 2 * (Si n( RB(i) a 2)RL(i) = (y(i) - fyb1) / alIf Abs(RB(i) - B0) = 0.0000000001 And Abs(RL(i) - 10) 0 And SX(i) 0 The nRL(i) = pi -
8、RL(i)End IfNO = aHO = Sqr(SX(i) a 2 + SY(i) a 2 + SZ(i) a 2) - Sqr(a * sb)BO = Atn(SZ(i) / (Sqr(SX(i) a 2 + SY(i) a 2) * (1 - e a 2 * NO / (NO + HO)DoNi = a / Sqr(1 - e a 2 * (Si n( BO) a 2)RH(i) = Sqr(SX(i) a 2 + SY(i) a 2) / Cos(BO) - NiRB(i) = Atn( SZ(i) / (Sqr(SX(i) a 2 + SY(i) a 2) * (1 - e a 2 * Ni / (Ni + RH(i)If Abs(RB(i) - BO)
