线路任意点坐标计算、及任意点对应桩号，左右偏距计算程序
卡西欧4800、4850系列计算器测量计算程序
一、字母含义；
K: 表示拟计算位置的线路桩号 ;
H: 表示计算位置距路线中心线的偏距，左偏为正，右偏为负。
T: 各段线路上作为起算点处的切线
方位角。
S: 拟计算点到起算点的曲线长。
L: 在圆曲线上表示曲线长，在
缓和曲线上表示缓和段长度，在直线上为零。
R: 表示曲线半径，左偏为正，右偏为负
E、F: 起算点的坐标值。
M”X1”N”Y1”: 已知点坐标，求其对应位置桩号及左右偏距。
二、程序
1、 坐标计算(COORD)
{K,H}:KH“L+，R-”:Prog “DATA”:”(X,Y)=”:
X=X+HSinW :Pause 1:Y=Y-HcosW:
2、坐标反算线路桩号(FS ZH)
Fix 4:M”X1”N”Y1”:LbI 1:Prog “DATA”:PoI (M-X,N-Y):Fixm:J<0=>J=J+360:≠>J=J⊿ Abs(Sin(W-J ))=1=>”K=”:K: Pause 1 :
“L+,R- =”: H=ISin(W-J):≠>K=K+Icos(W-J): Goto 1
3、子程序Prog “DATA”
K<（第1段与第2段线路分界处的路线桩号）＝> T=(第1段起算点处的切线方位角值):S＝K－（第1段起算点处的路线桩号）:L＝（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第1段起算点的X坐标值）:F=（第1段起算点的Y 坐标值）
≠ >K<(第2段与第3段线路分界处的路线桩号):=>T=(第2段起算点处的切线方向角值):S=K-（第2段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第2段起算点的X坐标值）F=（第2段起算点的Y坐标值）
≠ >K<(第3段与第4段线路分界处的路线桩号):=>T=(第3段起算点处的切线方向角值):S=K-（第3段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第3段起算点的X坐标值）F=（第3段起算点的Y坐标值）
……
≠ >K<(第n-1段与第n段线路分界处的路线桩号):=>T=(第n-1段起算点处的切线方向角值):S=K-（第n-1段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第n-1段起算点的X坐标值）F=（第n-1段起算点的Y坐标值）
≠>T=(第n段起算点处的切线方向角值):S=K-（第n段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第n段起算点的X坐标值）F=（第n段起算点的Y坐标值）⊿⊿⊿⊿⊿⊿⊿⊿(⊿共n-1个)
L ≠ 0 = > C =90S^2/ЛRL : ≠ > C=0:⊿
L=S => V=2RsinC : Q=T-C : W=T-2C ≠> L=0 =>V=S : Q=T : W=T :≠> Pol(S-S^5/40R^2L^2+S^9/3456R^4L^4,S^3/6RL-S^7/336R^3L^3+S^11/42240R^5L^5):Fixm:V=I:Q=T-J:W=T-CΔΔ
X=E+VcosQ:
Y=F+VsinQ
三、计算示例（浙江
龙丽高速公路某标段线路）
计算坐标:COORD
{K,H}:KH“L+，R+”:Prog “DATA”:”(X,Y)=”:X=X+HSinW Pause 1:
Y=Y-HcosW:
计算桩号、偏距:FS ZH
Fix 4:Deg:M”X1”N”Y1”:LbI 1:Prog “DATA”:PoI (M-X,N-Y):Fixm:J<0=>J=J+360:≠>J=J⊿
Abs(Sin(J- W ))=1=>”K=”:K: Pause 1 :
“L+,R- =”: H=ISin(J-W):≠> K=K+Icos(J-W):Goto 1
数据库:DATA
K<74886.420=>T=151°44´18.1":S=K-73544.990:R=-2500:E=155009.949:F=439722.191: S<0=>L=S:≠>L=0Δ
≠>K<77616.980=>T=131°13´39" : S=K-76676.330:R=5000:E=152435.423: F=441466.078: S<0 =>L=S:≠>L=0Δ
≠>K<80984.800=>T=152°13´22.7" : S=K-79449.180:R=-5000:E=150385.156: F=443302.124:S<0=>L=S:≠>L=0Δ
≠>T=134°00´16.4":S=K-82574.660:R=5000:E=147760.219:F=444968.04: S<0=>L=S:≠>L=0ΔΔΔΔ
L ≠ 0 => C =90S^2/ЛRL : ≠> C=0: Δ
L=S => V=2RsinC : Q=T-C : W=T-2C ≠> L=0 =>V=S : Q=T : W=T :≠> Pol(S-S^5/40R^2L^2+S^9/3456R^4L^4,S^3/6RL-S^7/336R^3L^3+S^11/42240R^5L^5):Fixm:V=I:Q=T-J:W=T-CΔΔ
X=E+VcosQ: Y=F+VsinQ
本回答被网友采纳