任意里程中边桩坐标进行正反算(改版1)

《任意里程中边桩坐标进行正反算(改版1)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、任意里程中边桩坐标进行正反算一、程序功能    本程序可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、止点里程、起点切线方位角、起点曲率、止点曲率、道路中线的前进方向）及里程边距夹角或坐标，对该曲线段范围内(如不在该范围则自动转入下一线元)任意里程中边桩坐标进行正反算。二、源程序1.主程序(FQXZB)453字节Lbl0:N"1.LC=>XY,XY=>LC":S:“LC”:prog“QXYS”:U"X0":V"Y0":O"LC1":H"LC2":G"Fi":P"R1":R"

2、R2":Q"ZX[Z:0,L:-1,R:+1]":Z[1]=(1/R-1/P)/(2Abs(H-O)):Z[2]=180/π:N=1=>Goto1:≠>Goto 3⊿←┘Lbl 1:{ABK}:S"LC":A"B-L":B“B-R””K"KJ[Q(+),H(-)":W"JA":D"PC[L(-1),R(+1)]":S≥H=>S=H⊿Z[3]=Abs(S-O):Prog"SUB1":“X=”:X◢“Y=”:Y◢“i=”:F=G+QZ[2]Z[3](1/P+Z[1]Z[3]):F≥0=>F→DMS◢Z:≠>F=F+36

3、0:F→DMS◢Prog"DXD":Goto2←┘≠>F=F+360:F→DMS⊿Prog"DXD":Goto 2Lbl 2："X[L]=":I=X-Acos(F+W)+KcosF+Dcos(F+90)◢"Y[L]=":J=Y-Asin(F+W)+KsinF+Dsin(F+90)◢Pol(I-E,J-L):J＜0=>J=J+360⊿"FWJ[L]=":J→DMS◢"S[L]=":I◢"X[Z]=":I=X+KcosF+Dcos(F+90)◢"Y[Z]=":I=X+KsinF+Dcos(F+90)◢Pol(I-E,J

4、-L):J＜0=>J=J+360⊿"FWJ[Z]=":J→DMS◢"S[Z]=":I◢"X[R]=":I=X+Bcos(F+W)+KcosF+Dcos(F+90)◢"Y[R]=":J=Y+Bsin(F+W)+KsinF+Dsin(F+90)◢Pol(I-E,J-L):J＜0=>J=J+360⊿"FWJ[R]=":J→DMS◢"S[R]=":I◢S=H≠>Goto4←┘Lbl 3:{XY}:XY:I=X:J=Y:Prog"SUB2":"LC=":S=O+Z[3]◢"B[L(-1),R(+1)]=":Z◢Goto 3L

5、bl 4:U=X:V=Y:O=H:G=G+180Q(Z[1]Z[3]+1/P)Z[3]/π5:P=RGoto 02. 子程序(SUB1)414字节Z[4]=0.1739274226:Z[5]=0.3260725774:Z[6]=0.0694318442:Z[7]=0.3300094782:F=1-Z[7]:M=1-Z[6]:Z[8]=G+QZ[2]MZ[3](1/P+MZ[1]Z[3]):X=U+Z[3](Z[4]cos(G+QZ[2]Z[3]Z[6](1/P+Z[1]Z[3]Z[6]))+Z[5]cos(G+QZ

6、[2]Z[3]Z[7](1/P+Z[1]Z[3]Z[7]))+Z[5]cos(G+QZ[2]FZ[3](1/P+FZ[1]Z[3]))+Z[4]cosZ[8]):Y=V+Z[3](Z[4]sin(G+QZ[2]Z[3]Z[6](1/P+Z[1]Z[3]Z[6]))+Z[5]sin(G+QZ[2]Z[3]Z[7](1/P+Z[1]Z[3]Z[7]))+Z[5]sin(G+QZ[2]FZ[3](1/P+FZ[1]Z[3]))+Z[4]sinZ[8])3. 子程序(SUB2)180字节T=G-90:Z[3]=Abs((Y

7、-V)cosT-(X-U)sinT):Z=0:Lbl 0:Prog"SUB3":Z[7]=T+QZ[2]Z[3](1/P+Z[1]Z[3]):Z=(J-Y)cosZ[7]-(I-X)sinZ[7]:AbsZ<1E-6=>Goto1:≠>Z[3]=Z[3]+Z:Goto 0⊿←┘Lbl 1:Z=0:Prog"SUB3":Z=(J-Y)÷sinF4. 子程序(SUB3)55字节Prog"SUB1":F=G+QZ[2]Z[3](1/P+Z[1]Z[3])+90:X=X+ZcosF:Y=Y+ZsinF5.子程序(QXYS)

8、S≥0=>GOTO0←┘Lbl0:S≤313.137=>U=555357.153:V=4797782.079:O=0:H=313.137:G=106°6′35.8″:P=10^45:R=10^45:Q=0:Goto 1：⊿S＜443.137=>U=555270.263:V=480082.919:O=313.137:H=443.137:G=106°5′35.8