线路任意点坐标计算、及任意点对应桩号,左右偏距计算程序.tmp

线路任意点坐标计算、及任意点对应桩号,左右偏距计算程
序.tmp
线路任意点坐标计算、及任意点对应桩号，左右偏距计算程序卡西欧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):Fi xm: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<（第１段与第２段线路分界处的路线桩号）＝> T=(第１段起算点处的切线方位角值):Ｓ＝Ｋ－（第１段起算点处的路线桩号）:Ｌ＝（在圆曲线上等于Ｓ；在缓和曲线上等于缓和段长度；在直线上为零。

）: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
本程序在高速公路测量中,由于曲线半径较大,其精度完全能满足设计要求,如果有兴趣的朋友不妨把程序中的Pol(S-S^5/40R^2L^2,S^3/6RL-S^7/336R^3L^3):函数变为
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):可满足更高要求的放样精
度.。