坐标正反算(带高程.docVIP

Lb1 9 100→DimZ:”ZS=1,FS=2”?Q:”K+”?M If Q=1:Then Goto A:Else Goto S:IfEnd Lb1 S “X1=”?S “Y1=”?T Lb1 A If M23285.856（第一缓和曲线起点）:Then 22396.61（起点交点）→Z:3049173.247（起点X坐标）→A:121°1°16.97°（起点方位角）→C: C:236818.413（起点Y坐标）→B:1×1020 （起点半径）→R:0（右转输入0，左转输入1）→F:0（起点缓和曲线长）→L:Goto 0:IfEnd If M23647.847（第二缓和曲线起点或圆曲线终点）:Then 23285.856（第一缓和曲线起点）→Z:3048706.061（起点X坐标）→A:121°1°16.97°（起点方位角）→C:237595.285（起点Y坐标）→B:640（起点半径）→R:0（右转输入0，左转输入1）→F:180（缓和曲线长）→L:Goto 0:IfEnd If M23827.847（第二缓和曲线终点）:Then 3048473.122（起点X坐标）→A:237868.071（起点Y坐标）→B: 145°22°16.81°（线元方位角）→C :(1÷640)（上一段曲率半径） →D: 1÷(1×1020)（下一段曲率半径） →E:23647.847（第二缓和曲线起点或圆曲线终点）→Z:23827.847→H:Goto H:IfEnd Lb1 0 If Q=2:Then 0→U:0→O: Else “U=”?U: “O=”?O:IfEnd M-Z→Z[12] If √(Z[12]2)L:Then Goto 1:IfEnd If √(Z[12]2)≥L:Then Goto 2:IfEnd Lb1 1 Z[12]- Z[12]5÷40÷R2÷L2→Z[18] (-1)F×((Z[12]3÷6÷R÷L)- Z[12] 7÷3366÷R3÷L3) →Z[19] A+ Z[18]cos(C)- Z[19]sin(C) →N B+ Z[18]sin(C)+Z[19]cos(C) →W C+90×(-1)F×Z[12]2◢(LRπ) →Z[14] If Q=2:Then Goto F:IfEnd N+Ucos(Z[14]+O) →X W+Usin(Z[14]+O) →Y “JSGC=1”?→Z[25] If Z[25]=1:Then Prog “GAOCHEN”:IfEnd “X=”:X◢ “Y=”:Y◢ Goto 3 Lb1 2 ((Z[12]-L) ÷R+L÷(2R)) ×180÷π→Z[13] Rsin(Z[13]+L÷2-L3÷240÷R2→Z[18] (-1)F×(R×(1-cos(Z[13]))+ L2÷(24R)- L4÷2688÷R3) →Z[19] A+ Z[18]cos(C)- Z[19]sin(C) →N B+ Z[18]sin(C)+Z[19]cos(C) →W If R=1×1020:Then C→Z[14]:Else C+90×(-1) F×(2Z[12]-L) ÷R÷π→Z[14]:IfEnd If Q=2:Then Goto F:IfEnd N+Ucos((-1)F×Z[13]+C+O) →X W+Usin((-1)F×Z[13]+C+O) →Y “JSGC=1”?→Z[25] If Z[25]=1:Then Prog “GAOCHEN”:IfEnd “X=”:X◢ “Y=”:Y◢ Goto 3 Lb1 H (E-D) ÷Abs(H-Z) →Z[15] Abs(M-Z) →Z[16] Z[15] ×Z[16]→Z[17] (C+90×Z[16](Z[17]+2×D) ÷π) → Z[14] (C+15×Z[16](Z[17] ÷12+2×D) ÷(2π) → Z[1] (C+45×Z[16](3×Z[17] ÷12+2×D) ÷(2π) → Z[2] (C+75×Z[16](5×Z[17] ÷12+2×D) ÷(2π) → Z[3] (C+105×Z[16](7×Z[17] ÷12+2×D) ÷(2π) → Z[4] (C+135×Z[16](9×Z[17] ÷12+2×D) ÷(2π) → Z[5] (C+165×Z[16](11×Z[17] ÷12+2×D) ÷(2π) → Z[6] (C+15×Z[16]( Z[17] ÷6+2×D) ÷π → Z[7] (C+30×Z[16](2 Z[17] ÷6+2×D) ÷π → Z[8] (C+45×Z[16]( 3Z[17] ÷6+2×D) ÷π → Z[9] (C+60×Z[16](4 Z[17] ÷6+2×D) ÷π → Z[10]