公路中桩坐标计算(按线路前进方向执行)

公路中桩坐标计算（按线路前进方向执行）
主程序
MAIN----------CD
A”X[0]“：B“Y[0]”：F“F[0]”：E“K[0]”⍓
P=1=>Prog“ZX”:≠>P=2=> Prog“YQX”: ≠>P=3=> Prog“HHQX1”: ≠>P=4=> Prog“HHQX2”
X[0]---起点坐标X
Y[0]---起点坐标Y
F[0]---起点方位角
K[0]---起点桩号
P-----线形判断(P=1直线段、P=2圆曲线段、P=3第一段缓和曲线、P=4第二段缓和曲线)
子程序1
ZX-------------CD
LbI 0：{K}:X=A+Rec((K-E),F) ⍓
Y=B+J⍓
Goto 0
K---计算点桩号
E---起算点桩号
子程序2
YQX-------------CD
Lb1 1：{K}：Q=180（K-E）÷（ЛR）：W=0=>X=A+Rec((2Rsin(Q÷2))，（F-Q÷2）)⍓
Y=B+J⍓
H=F-(K-E)×180÷（ЛR）⍓
Goto1⍓
≠>W=1=>X=A+Rec((2Rsin(Q÷2)),(F+Q÷2))⍓
Y=B+J⍓
H=F+(K-E)×180÷（ЛR）⍓
Goto1
W---判断曲线偏向（W=0偏右、W=1偏左）
H---计算点该点的切线方位角
子程序3
HHQX1------------CD
Lbl2{K}:S=K-E-(K-E)^5/(40R^2L^2)+(K-E)^9/(3456(R^4L^4):T=(K-E) ^3/(6RL)-(K-E)^7/(336(R^3L^3)+(K-E)^11/(42240R^5L^5)
G=180(K-E)^2/(6ЛRL)
W=0=>X=A+ScosF+TsinF⍓
Y=B+SsinF-TcosF⍓
H=F-3G⍓
Goto2
≠>W=1=>X=A+ScosF-TsinF⍓
Y=B+SsinF+TcosF⍓
H=F+3G⍓
Goto2
子程序4
HHQX2------------CD
Lbl
3:{K}:S=(K-E)-(K-E)^3/(6R^2)+(K-E)^4/(8R^2L)+(L^2-3R^2)(K-E)^5/( 120R^4L^2)-(K-E)^6/(72R^4L)+(K-E)^7/(112R^4L^2)-(K-E)^8/(384R^ 4L^3)+(K-E)^9/(3456R^4L^4)
T=(K-E)^2/(2R)-(K-E)^3/(6RL)-(K-E)^4/(24R^3)+(K-E)^5/(20R^3L)+( L^2-15R^2)(K-E)^6/(720R^5L^2)+(R^2-L^2)(K-E)^7/(336R^5L^3)+(K-E)^8/(384R^5L^2)-(K-E)^9/(864R^5L^3)+(K-E)^10/(3840R^5L^4)-(K-E)^11/(42240R^5L^5)
G=180(K-E)^2/(6ЛRL)
W=0=>X=A+ScosF+TsinF⍓
Y=B+SsinF+TcosF⍓
H=F-3G⍓
Goto3
≠>W=1=>X=A+ScosF-TsinF⍓
Y=B+SsinF+TcosF⍓
H=F+3G⍓
Goto3
符号同上本程序适应自由边桩，边桩程序由各人习惯而定，在子程序
后加上边桩程序即可。

本程序特点在于不论在线路逐桩坐标计算时，线形改变而计算的路线不会改变。