中线坐标计算改进版

1、公路中线坐标计算改进版王庆军ZHYHZ主程序1. Deg:Fix3:12DimZ2. “ZHYHZ:Z1-N”?Z93. IfZ9=0:Then “X1”?D:“Y1”?E:“X2”?F:“Y2”?G:“X3”?H:“Y3”?V:“JD2”?W:“Ls1”?L:“Ls2”?T:“R”?R:IfEnd4. Cls:“1:KBXY:2:XYKB”?Z10:Z10=2Goto35. “CZX”?X:“CZY”?Y6. Lbi0:“Kp”?K7. Z90Prog“XLZB”:Prog“AA”8. Lbi1:90Z:“JJ=”?Z:Z=0Goto0:“JB=”?B9.If B0：Then RecAbsB

2、，Q+Z180：Else RecAbsB，Q+Z：IfEnd10.M+IM：N+JN11.PolMX，NY：Prog“BB”：JJ: II12.1C:Prog“ZHYHZ-XS”13.Lbi2:0S:“S=”?S:I-SSIf Abs(s)0.01:Then Goto2:Else Goto1:EfEnd14.Lbi3:“Xb=”?M:“Yb=”?N:“Kp=”?K15.0S:0B:MZ11:NZ1216.Lbi4:Z90Prog“XLZB”:Prog“AA”17.Pol(Z11-M,Z12-N):If J0:Then J+360J:IfEnd:JJ:II18.ICos(J-Q)S:K+SK:I

3、f Abs(S)0.0001:Then Goto4:IfEnd19.ISin(J-Q)B20.2C:Prog“ZHYHZ-XS”:Goto3 ZHYHZ-XS显示子程序1. Cis:If C=1:Then “X=”:Locate4,1,M:“Y=”:Locate4,2,N:Fix4:“FJ=”:Locate4,3,J°:“HD=”:Locate4,4,I2. If C=2:Then “K=”:Locate4,2,K:“B=”:Locate4,3,BIfEndXLZB路由器1. If Z9=1:Then Prog“XLZB-1”:Goto1:IfEnd2. If Z9=2:Then P

4、rog“XLZB-2”:Goto1:IfEnd3. If Z9=3:Then Prog“XLZB-3”:Goto1:IfEnd4. If Z9=4:Then Prog“XLZB-4”:Goto1:IfEnd5. Cls:Locate4,2,“ZHYHZ-ERROR”Sotp6. Lbi1XLZB-1数据库1.IfK×××××××ThenGotoA：ElseIfK×××××××AndK××××××&

5、#215;ThenGotoBElseIfK×××××××AndK×××××××：ThenGotoC：ElseIfK×××××××ThenGotoNEfEndEfEndEfEnd(N个)2.Cls:Locate 6,2,“KPOUT” Stop3LbiA(X1)×××××××D(Y1)×××

6、××××E(X2)×××××××F(Y2)×××××××G(X3)×××××××H(Y3)×××××××V(JD2)×××××××W(LS1)×××××

7、5;×L(LS2)×××××××T(R)×××××××RGoto04LbiB(X1)×××××××D(Y1)×××××××E(X2)×××××××F(Y2)×××××××

8、G(X3)×××××××H(Y3)×××××××V(JD2)×××××××W(LS1)×××××××L(LS2)×××××××T(R)×××××××RGoto05LbiC(X1)&

9、#215;××××××D(Y1)×××××××E(X2)×××××××F)(Y2)×××××××G(X3)×××××××H(Y3)×××××××V(JD2)×××

10、;××××W(LS1)×××××××L(LS2)×××××××T(R)×××××××RGoto0NLbiN(X1)×××××××D(Y1)×××××××E(X2)××××

11、×××F(Y2)×××××××G:(X3)×××××××H(Y3)×××××××V(JD2)×××××××W(LS1)×××××××L(LS2)××××××&#

12、215;T(R)×××××××RGoto0Lbi0AA坐标计算子程序1.PolFD，GE：Prog“BB”：JO：PolHF，VG：Prog“BB”：JU UOZ1：Z1÷AbsZ1P L2÷24RL4÷2688R3Z54.L÷2L3÷240R2Z65.T2÷24RT4÷2688R3Z76.T÷2T3÷240R2Z87.R+Z5tanAbsZ1÷2+Z6Z5Z7÷SimAbsZ1Z38.R+Z7tanAbsZ1

13、47;2+Z8+Z5Z7÷SimAbsZ1Z49.RAbsZ1÷180+L+T÷2Z2：Prog“CC”BBIf J0：Then 360+JJ：IfEndCC If KWZ3：Then Goto1：Else If KWZ3And KWZ3+L：Then Goto2：Else If KWZ3+L And KWZ3+Z2T：Then Goto3：Else If KWZ3+Z2T And KWZ3+Z2：Then Goto4：Else If KWZ3+Z2：Then Goto5：IfEnd：IfEnd：IfEnd：IfEnd：IfEnd Lbi1：WZ3KS：FZ3+S

14、CosOM：GZ3+SSimON：OQ：GotoA Lbi2：KW+Z3S：RLZ1：Prog“KK”：Prog“JJ”：GotoA Lbi3：KW+Z3S：902SL÷RZ1：L÷2L3÷240R2+RSimZ1I：L2÷24R+R1CosZ1J：Prog“JJ”：GotoA Lbi4：WZ3+Z2KS：RTZ1：Prog“KK”：F+Z4ICosU+JCosU+90PM：G+Z4ISimU+JSimU+90PN：UPZ1Q：GotoA Lbi5：KW+Z3Z2S：F+Z4+SCosUM：G +Z4+SSimUN：UQ：GotoA LbiAJJ FZ3

15、ICosO+JCosO+90PM：GZ3ISimO+JSimO+90PNO+PZ1QKK SS5÷40Z12+S9÷3456Z14S13÷599040Z16IS3÷6Z1S7÷336Z13+S11÷4220Z15J 90S2÷(Z1)Z1说明：本程序可计算前HZ后ZH上所有任意中桩，边桩座标：提问：X？Y？为测站点坐标Z？为边桩与路线交角K？为待求点桩号B？为距离中桩距离左“-B”，右“+B”D，E为前一交点桩座标，F，G为本交点桩座标，H，V为后交点桩座标，W为本交点桩号，L为本曲线前缓和曲线长，T为本曲线后缓和曲线长，R

16、为本曲线半径，C为后曲线直缓点桩号。输出：M，N为待求点K所对应的X，Y坐标，J为方位角，I为放线平距。当平曲线无缓和曲线时L，T均输入零。本程序中的数学模型是根据公路工程计算手册中桩大地坐标而建立，为了简化程序，本程序采用了一个主程序ZHYHZ和假设干个子程序组成。匝道全线贯穿程序第三版王庆军TYQXJS主程序1. Deg:Fix3:6Dim2. “TYQXJS:1-N”?Z53. IfZ5=0:Then “X0”?U:“Y0”?V:“S0”?O:“F0”?G:“Ls”?H:“Rs”?P:“Rn”?R:“Q”?Q:If End4. Cls:“1:SZXY,2:XYSZ”?Z6:Z6=2Got

17、o25. “CZX=”?M:“CZY=”?N6. Lbi0:“Kp”K7. Z50Prog“TYQXJS-DATA”:Abs(K-O)W8. Prog“SUB1”9. Lbi1:90Z:“JJ”Z:Z=0Goto0:“JB”B10. If B0:Then Res(AbsB,Z4+Z-180):Else Rec(AbsB,Z4+Z):IfEnd11. X+IX:Y+JY12. Pol(X-M,Y-N):J0J+360J13. 1T:Prog“TYQXJS-XS”14. Lbi3:“S=”?S:I-SS If Abs(S)0.01:Then Goto3:Else Goto1:IfEnd15. Lb

18、i2:“X”?X:“Y”?Y:“Kp”?K16. Z50Prog“TYQXJS-DATA”17. Prog“SUB2”18. 2T:Prog“TYQXJS-XS”SUB1坐标正算子程序11÷PC：P-R÷2HPRD：180÷E2.0.1739274226A:0.3260725774Z1:0.0694318442S:0.3300094782L:1-LZ2:1-SZ33.U+W(ACos(G+QESW(C+SWD)+Z1Cos(G+QELW(C+LWD)+Z1Cos(G+Z2QEW(C+Z2WD)+A Cos(G+Z3QEW(C+Z3WD)X4.V+W(ASin(G+

19、QESW(C+SWD)+Z1Sin(G+QELW(C+LWD)+Z1Sin(G+Z2QEW(C+Z2WD)+A Sin(G+Z3QEW(C+Z3WD)Y5.G+QEW(C+WD)Z4SUB2反算里程和边距子程序1 XI：YJ2 Abs(J-V)Cos(G-90)-(I-U)Sin(G-90)W:0B3 Lbi0:Prog”SUB1”4 (G-90)+QEW(C+WD)L5 (J-Y)CosL-(I-X)SinLB6 If AbsB10(-6):Then Goto1:Else W+BW:Goto0:IfEnd7 Lbi1:0B:Prog“SUB1”8 (J-Y)÷Sin(Z4+90)B

20、9 O+WK:BBTYQXJS-DATA路由器1.If Z5=1:Then Prog“AZD”:Goto1:IfEnd2.If Z5=2:Then Prog“BZD”:Goto1:IfEnd3.If Z5=3:Then Prog“CZD”:Goto1:IfEnd4.If Z5+4:Then Prog“DZD”:Goto1:IfEnd5.Cls:Locate 4,2,“TYQXJS-ERROR” Stop6.Lbi1TYQXJS-XS显示子程序1. Cls:If T=1:Then “X=”:Locate 4,1,X:“Y=”:Locate 4,2,Y:Fix4:“FJ=”:Locate 4,3,

21、J°:“HD”:Locate 4,4,I EfEnd2. If T=2:Then “K=”:Locate 4,1,K:“B=”:Locate 4,2,B IfEndBZD(以B,C匝道为例的数据库)1 If K102.511:Then GotoA:Else If K102.511 And K272.465:Then GotoB:Else If K272.465:Then GotoC:IfEnd:IfEnd:IfEnd2LbiA:2904.645P:58R:102.511H:1Q:0O:3030137.263U:506569.752V:319°1139.7G:Goto03Lb

22、iB:58P:58R:169.954H:1Q:102.511O:3030228.025U:506530.533V:10°5018.1G:Goto04LbiC:58P:1045R:76.371H:1Q:272.465O:3030218.405U:506645.484V:178°4343.5G:Goto05Lbi0CZD1.If K212.645:Then GotoA:Else If K212.645 And K314.337:Thed GotoB:Else If K314.337 And K469.526:Then GotoC:Else IfK469.526:Then Got

23、oD:IfEnd:IfEnd:IfEnd:IfEnd2.LbiA:4171.441P:4171.441R:312.645H:1Q:0O:3029949.692U:506489.917V:36°3306G:Goto03.LbiB:4171.441P:62R:101.692H:1Q:212.645O:3030117.214U:506620.855V:39°2820.6G:Goto04.LbiC:62P:62R:155.189H:1Q:314.337O:3030173.052U:506702.128V:87°931.2G:Goto05.LbiD:62P:2457.229R:185.628H:1Q:469.5260:3030063.238U:506744.577V:230°3420.7G:Goto06.Lbi0本程序为怀通公路竹田BC匝