《任意里程中边桩坐标进行正反算(改版1)》由会员分享，可在线阅读，更多相关《任意里程中边桩坐标进行正反算(改版1)（5页珍藏版）》请在金锄头文库上搜索。
1、1任意里程中边桩坐标进行正反算任意里程中边桩坐标进行正反算一、程序功能一、程序功能本程序可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型） 的线元要素（起点坐标、起点里程、止点里程、起点切线方位角、起点曲率、 止点曲率、道路中线的前进方向）及里程边距夹角或坐标，对该曲线段范围内 (如不在该范围则自动转入下一线元)任意里程中边桩坐标进行正反算。 二、源程序二、源程序 1.1. 主程序主程序(FQXZB)(FQXZB) 453453 字节字节Lbl0:N“1.LC=XY,XY=LC“:S:“LC”:prog“QXYS”:U“X0“:V“Y0“:O“LC1“:H“LC 2“:G“Fi“:P“R1
2、“:R“R2“:Q“ZXZ:0,L:-1,R:+1“:Z1=(1/R-1/P)/(2Abs(H- O):Z2=180/:N=1=Goto 1:Goto 3 Lbl 1:ABK:S“LC“:A“B-L“:B“B-R”K“KJQ(+),H(-)“:W“JA“:D“PCL(-1),R (+1)“:SH=S=H Z3=Abs(S-O):Prog“SUB1“:“X=”:X “Y=”:Y “i=”:F=G+Q Z2Z3(1/P+Z1Z 3):F0=FDM S Z: F=F+360:FDMS Prog“DXD“:Goto 2 F=F+360: FDMS Prog“DXD“:Goto 2 Lbl 2：“XL=
3、“:I=X-Acos(F+W)+KcosF+Dcos(F+90) “YL=“:J=Y-Asin(F+W)+KsinF+Dsin(F+90) Pol(I-E,J- L):J0=J=J+360 “FWJL=“:JDMS “SL=“:I “XZ=“:I=X+KcosF+Dcos(F+90) “YZ=“:I=X+KsinF+Dcos(F+90) Pol(I- E,J-L):J0=J=J+360 “FWJZ=“:JDMS “SZ=“:I “XR=“:I=X+Bcos(F+W)+KcosF+Dcos(F+90) “YR=“:J=Y+Bsin(F+W)+KsinF+Dsin(F+90) Pol(I-E,J-
4、 L):J0=J=J+360 “FWJR=“:JDMS “SR=“:I S=H Goto 4 Lbl 3:XY:XY:I=X:J=Y:Prog“SUB2“:“LC=“:S=O+Z3 “BL(-1),R(+1)= “:Z2Goto 3 Lbl 4:U=X:V=Y:O=H:G=G+180Q(Z1Z3 +1/P)Z3/:P=R Goto 0 2.2.子程序子程序(SUB1)(SUB1) 414414 字节字节Z4=0.1739274226:Z5=0.3260725774:Z6=0.0694318442:Z7=0.33000947 82:F=1-Z7:M=1- Z6:Z8=G+QZ2MZ3(1/P+M
5、Z1Z3):X=U+Z3(Z4cos(G+QZ2Z3Z6( 1/P+Z1Z3Z6)+Z5cos(G+QZ2Z3Z7(1/P+Z1Z3Z7)+Z5cos (G+QZ2FZ3(1/P+FZ1Z3)+Z4cosZ8):Y=V+Z3(Z4sin(G+QZ2Z 3Z6(1/P+Z1Z3Z6)+Z5sin(G+QZ2Z3Z7(1/P+Z1Z3Z7)+Z5sin(G+QZ2FZ3(1/P+FZ1Z3)+Z4sin Z8) 3.3. 子程序子程序(SUB2)(SUB2) 180180 字节字节T=G-90:Z3=Abs(Y-V)cosT-(X- U)sinT):Z=0:Lbl 0:Prog“SUB3“:Z7=
6、T+QZ2Z3(1/P+Z1Z3):Z=(J- Y)cosZ7-(I-X)sinZ7:AbsZGoto1: Z3=Z3+Z:Goto 0 Lbl 1:Z=0:Prog“SUB3“:Z=(J-Y)sinF 4.4.子程序子程序(SUB3)(SUB3) 5555 字节字节Prog “SUB1“:F=G+QZ2Z3(1/P+Z1Z3)+90:X=X+ZcosF:Y=Y+ZsinF 5.5. 子程序子程序(QXYS)(QXYS) S0=GOTO 0 Lbl 0:S313.137=U=555357.153:V=4797782.079:O=0:H=313.137:G =106635.8:P= 1045:R
7、=1045:Q=0:Got o 1： S443.137=U=555270.263:V=480082.919:O=313.137:H=443.137:G=106 535.8: P=1045:R=500:Q=-1: Goto 1： Goto 1 Lbl 1 6.6. 子程序子程序(DXD)(DXD) C“DXD“:C=1=“XZD-5=“:E=491893.218 “YZD-5=“:L=410802.542 “HZD- 5=“:11.239 C=2=“XZD-9=“:E=491883.398 “YZD-9=“:L=410730.632 “HZD- 5=“:10.896 C=3=“XZX-5=“:E
8、=492673.735 “YZX- 5=“:L=410471.185 “HZX-5=“:9.574 C=4=“XZX-12=“:E=492523.316 “YZX-12=“:L=410520.292 “HZX- 5=“:10.045 C=0=“XZX-316=“:“YZX-16=“ 三、使用说明三、使用说明 1 1、规定：、规定：(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元 往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z 取负值；当位于中 线中线右侧时，Z 取正值。(3) 当线元为直线时，其起点、
9、止点的曲率为 0。(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率均等 于圆弧的半径。(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率为 10E45；与圆 曲线相接时，曲率等于圆曲线的半径的倒数。止点与直线相接时，曲率为 10E45；与圆曲线相接时，曲率等于圆曲线的半径的倒数。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率等于设计规定 的值；与圆曲线相接时，曲率等于圆曲线的半径的倒数。止点与直线相接时， 曲率等于设计规定的值；与圆曲线相接时，曲率等于圆曲线的半径的倒数。 2 2、输入与显示说明、输入与显示说明输入部分：输入部分：N ? 选择计算方式，输入 1 表示进行由里程、边距计算坐标 ；输入非 1 值表示由坐标反算里程和边距。X0 ？线元起点的 X 坐标 Y0 ？线元起点的 Y 坐标 LC1 ？线元起点里程Fi ？线元起点切线方位角LC2 ？线元止点里程R1 ？线元起点曲率R2 ？线元止点曲率Q ？线元左右偏标志(左偏 Q=-1，右偏 Q=1，直线段 Q=0
