卡西欧5800坐标转换

CASIO5800坐标正反算计算程序CASIO 9860/5800坐标正反算计算程序1. 坐标正算主程序(命名为ZBZS)第1行：Lbl 0:”K=”?K:”BIAN=”? Z:”α=”?B第2行：Prog “A”第3行：”X=”:N+Zcos(F+B)◢第4行：”Y=”:E+Zsin(F+B)◢第5行：”F=”:F◢第6行：Goto 0K——计算点的里程BIAN——计算点到中桩的距离（左负右正）α——取前右夹角为正2. 坐标反算桩号和偏距主程序(命名为ZBFS)第1行：”X1=”? C:”Y1=”?D: ”K1=”?K第2行：Lbl 0:Prog “A”第3行：Pol(C-N,D-E)第4行：List Ans[1]→I第5行：List Ans[2]→J第6行：Icos(F-J)→S:K+S→K第7行：Abs(S)>0.0001=>Goto 0第8行：”K1=”:K◢第9行：”BIAN=”:Isin(J-F)→Z◢X1——取样点的X坐标Y1——取样点的Y坐标K1——输入时为计算起始点(在线路内即可)，输出时为反算点的桩号Z——偏距（左负右正）3. 计算坐标子程序(命名为XYF)为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1行：K-A→S：(Q-P)÷L→I第2行：N+∫(cos(F+X(2P+XI)×90÷π),0,S)→N第3行：E+∫(sin(F+X(2P+XI)×90÷π),0,S)→E第4行：F+S(2P+SI)×90÷π→F4. 数据库（命名为A）第1行：K≤175.191=>Stop第2行：175.191→A:428513.730→N:557954.037→E:92°26′40″→F:0→P:1/ 240→Q:70.417→L:K≤A+L =>GoTo 1（第一缓和曲线，圆半径为240）第3行：245.607→A: 428507.298→N:558024.092→E: 100°50′59.4″→F: 1/240→P:1/240→Q:72.915→L: K≤A+L =>Goto 1（第圆曲线，半径为240）第4行：318.522→A: 428482.988→N:558092.538→E: 118°15′25.2″→F: 1/240→P: 0→Q: 55.104→L: K≤A+L =>Goto 1（第二缓和曲线，圆半径为240）第5行：373.627→A:428453.283→N:558138.912→E:124°50′4.5″→F:0→P:-1/180→Q:67.222→L:K≤A+L=>Goto 1:Stop（第一缓和曲线，圆半径为180）第6行：Lbl 1:Prog “XYF”A——曲线段起点的里程N——曲线段起点的x坐标E——曲线段起点的y坐标F——曲线段起点的坐标方位角P——曲线段起点的曲率(左负右正)Q——曲线段终点的曲率(左负右正)L——曲线段长度（尽量使用长度，为计算断链方便）说明：（1）在9860中，程序中所有公式和部分函数结果均存储在ListAns列表数组中，要想多次调用最好随公式取出结果，并赋给变量。

一、坐标正算基本公式 (02)二、坐标反算原理 (04)三、高程数据库录入变换 (05)四、计算器程序 (07)01、ZBZS （坐标正算） (07)02、ZBFS （坐标反算） (08)03、GCJF （高程积分） (09)04、PJFY （坡脚放样） (10)05、JFCX （积分程序） (11)06、ZBFY （坐标放样） (11)07、DT （递推） (12)08、HP （横坡） (13)09、LK （路宽） (14)10、SJKl （平面数据库） (14)IK SJK2 （纵面数据库） (14)12、SJK3 （左路宽度数据库） (15)13、SJK4 （右路宽度数据库） (15)14、SJK5 （横坡数据库） (16)15、SJK6 （下边坡数据库） (16)16、SJK7 （左上边坡数据库） (17)17、SJK8 （右上边坡数据库） (18)五、后记 (19)CASIO 5800计算器公路工程测量程序4 一、正算所涉及的计算公式图表1 在图1中，A点为回旋曲线起点，B点为回旋曲线止点，I点为所求坐标点。

设：A点的X坐标为XA , Y坐标为Y A, A点的切线方位角为α, A点的曲率为PA点的里程为L A, B点的曲率为p 8, B点的里程为L B, I点的曲率为p∣, I点的里A,程为L I O I点的切线角为B o由于回旋线上各点曲率半径Ri和该点至曲线起点的距离L成反比。

故此任意点的曲率为；Q =丄=土（C为常数）（1 ）R i C由式（1）可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为L点的曲率为；iL-I,、Pi = PA + （PB - PJ × I, _ /（2）当曲线右偏时P 8、PA取正值，反之取负值。

设：则有:P i = PA + ML在I 点处取一微段，则有：d β=- = P l d l （单位为弧度）Ri因已知回旋曲线起点A 的切线方位角α ,则里程为Li 点的切线方位角为:a, =a + βi将式（7）代入式（8）得：a i = a + + （单位为弧度）（9） 2对于式（9）,当P A =0, M 二0时，则a i =a ,式（9）变成计算直线段上任意 点切线方位角的计算公式；当P A=C （C 为常数）,M=O 时，则CIi 二a + P A L,式（9） 变成计算圆曲线上任意点切线方位角的计算公式。

大地坐标与施工坐标相互转换（CASIOfx-5800P）刚刚接触5800计算器，看到大地坐标与施工坐标相互转换是两个分开的。

现在为了学习一下5800程序输入，把它组合成一个程序文件，请大家试用指点。

程序清单:Fix 4:”ALFA=”?→J: ”AXO=”? →C: ”BYO=”? →D↙Lbl 3:”1→AG,2→GA”?Q↙Q=1＝＞Goto 1: Q=2＝＞Goto 2↙Lbl 1:”AA=”?→A:”BA=”?→B: ↙”XG=”:C+ACOS(J)-BSIN(J)▲”YG=”:D+ASIN (J)+BCOS (J)▲Goto 3↙Lbl 2:”XG=”?→X:”YG=”?→Y: ↙”AA=”( X-C) COS(J)+(Y-D)SIN(J)▲”BA=”: ( C-X) SIN (J)+(Y-D)COS (J)▲Goto 3↙说明：ALFA＝大地坐标与施工坐标夹角。

AXO，BYO＝施工坐标系的原点的大地坐标AA，BA＝施工坐标X,Y＝大地坐标例：运行程序显示输入ALFA＝? 大地坐标与施工坐标夹角13'0'0'AXO＝？施工坐标系的原点的大地坐标X2870340.477BYO＝？施工坐标系的原点的大地坐标Y610780.2221→AG,2→GA?1为施工转换大地，2为大地转换施工1AA=?要转换成大地坐标点的施工坐标A54.9BA=?要转换成大地坐标点的施工坐标B82.5计算结果XG=2870375.411YG=610872.9573回车1→AG,2→GA?1为施工转换大地，2为大地转换施工2XG=? 要转换成施工坐标点的大地坐标X 2870375.411YG=? 要转换成施工坐标点的大地坐标Y 610872.9573计算结果A=54.8995B=82.5001（欢迎参与右下角评论）公办中专（中央广播电视中等专业学校）常年滚动招生，滚动注册：专业齐全（中专），建筑类专业（可考二级建造师、安全员、造价员必备学历、征兵入伍可用）药剂类专业（可考执业药师）（注：费用低、出证快，教育部直属中专，可做学历认证，国家学信网永久查询。

卡西欧5800公路坐标正反算程序卡西欧5800公路坐标正反算程序是一种用于测量和计算公路坐标的工具。

它基于卡西欧5800全站仪的技术和功能，能够准确地测量和计算公路的坐标信息。

下面我将详细介绍卡西欧5800公路坐标正反算程序的工作原理和使用方法。

一、工作原理：卡西欧5800全站仪是一种高精度的测量仪器，它能够通过测量角度和距离来确定地点的坐标。

在公路测量中，我们通常会在公路上设置一系列的控制点，然后使用全站仪测量这些控制点的坐标。

通过这些测量数据，我们可以计算出公路上其他点的坐标。

卡西欧5800公路坐标正反算程序基于这些测量数据，通过一系列的计算和算法，能够准确地计算出公路上其他点的坐标。

具体的计算方法包括正算和反算。

正算是指已知控制点的坐标，通过测量数据和算法计算出其他点的坐标。

反算是指已知某个点的坐标，通过测量数据和算法计算出该点在公路上的位置。

二、使用方法：1. 设置控制点：首先，在公路上选择一些具有代表性的控制点，并使用全站仪进行测量。

在测量过程中，需要确保全站仪的水平和垂直仪器准确，以获得准确的测量数据。

2. 导入测量数据：将测量数据导入卡西欧5800公路坐标正反算程序。

该程序通常会提供一个界面，用于导入和管理测量数据。

3. 正算：如果需要计算公路上其他点的坐标，可以选择正算功能。

在正算功能中，需要输入已知控制点的坐标和测量数据。

程序将根据这些数据进行计算，并给出其他点的坐标。

4. 反算：如果需要计算某个点在公路上的位置，可以选择反算功能。

在反算功能中，需要输入已知点的坐标和测量数据。

程序将根据这些数据进行计算，并给出该点在公路上的位置。

5. 结果输出：计算完成后，卡西欧5800公路坐标正反算程序将给出计算结果。

结果通常以表格或图形的形式呈现，以便用户查看和使用。

三、注意事项：1. 测量准确性：为了获得准确的测量结果，使用卡西欧5800全站仪时需要注意测量准确性。

在测量过程中，需要保持仪器的稳定和准确，并避免测量误差。

1.3坐标转换(ZBZH)＂K=1－NO，≠YES＂？KIf K=1：Then Goto1：Else Goto 2：If EndLb1 1：＂XA=＂？A：＂YA=＂？B：＂XB=＂？C：＂YB=？D：＂XI=＂？E：＂YI=？F：＂XJ=＂？G：＂YJ=＂？HPol (C－A , D－B) ：J→O：Pol (G－E , H－F) ：J→U：O－U→TIf T﹤0：Then T+360→T：If End＂X0=＂：A-E×cos(T)+F×sin(T) →M▲“Y0=”：B－E×SIN T－F×COS T→N▲＂TA=＂：T DMS▲Goto 3Lbl 2：＂X0=＂？M：＂Y0=＂？N：＂TA=＂？TLbl 3：＂XX=＂？R：＂YY＂？V：＂P=1－SG,≠＂？PIf P=1：Then Goto 4：Else Goto 5：If EndLbl 4：＂X=＂：(R－M)×cos(T)+(V－N)×sin(T) →X▲＂Y=＂：－(R－M)×sin(T)+(V－N)×cos(T) →Y▲Goto 3Lbl 5：＂X=＂：R×cos(T) －V×sin(T)+M→X▲＂Y=＂：R×sin(T)+V×cos(T)+N→Y▲Goto 3说明：1、输入判断值K：a: K＝1时，输入A、B两点在两坐标系中的各自坐标A（XA，YA）、B（XB，YB）、A（XI，YI）、B （XJ，YJ），计算出施工坐标系的原点坐标在大地坐标系中的坐标值和X轴的夹角。

b…..K≠1时，直接输入施工坐标系原点在大地坐标系中的坐标（X0，Y0）和X轴的夹角TA。

2、输入计算点坐标P（XX，YY）。

3、输入判断值P：a: P=1时，计算出P点在施工坐标系中的坐标。

b: P≠1时，计算出P点在大地坐标系中的坐标。

卡西欧fx-5800 坐标正算程序“SO”算方位Lbl 0“D”：Pol(X-A，Y-B) →D◢If J≥0:Then J◢Else “J”:J+360◢IfEnd“ZZ”正算直线坐标A+Rec(Abs(P-K),Q) →XB+J→Y“X”：X+Rec(D,Q+N×90) →X◢“Y”:Y+J→Y◢P?所求桩号 A，B：起点x,y K:起点桩号 Q：方位角D：边桩距离 N？左偏“-1”右偏：“+1”X= Y= 所求点坐标“YY”正算圆曲线坐标Q+C×（P-K）÷S×180÷∏→Z“X”：A+C×S×（Sin(Z)-Sin(Q)）+Rec(D,Z+N×90) →X◢“Y”:B-C×S×（Cos(Z)-Cos(Q)）+J→Y◢P:所求点桩号 A，B：起点x,y K：起点桩号Q：方位角 S：曲线半径 D？边桩距离N？左偏“-1”，右偏“+1”X= y= 所算点坐标 C：线路左转“-1”右转“+1”“H1”正算缓和曲线坐标Abs(P-K) →LT2÷L→HL2×180÷∏÷2÷T2→VQ+C×V→ZL-L∧（3）÷40÷H2+L∧（5）÷3456÷H∧（4）-L∧（7）÷H∧（6）÷599040→EL2÷6÷H- L∧（4）÷336÷H∧（3）+ L∧（6）÷42240÷H∧（5）- L∧（8）÷H∧（7）÷（6×1612800）→FIf C<0：Then-F→F：IfEndA+E×Cos(Q)-F×Sin(Q) →X“X”：X+Rec(D,Z+M×N×90) →X◢B+E×Sin(Q)+F×Cos(Q) →Y“Y”：Y+J◢P？所求点桩号 A，B：起点x,y K:起算点里程桩号Q：前切线方位角 T（是圆曲线要素的A1）=√（根号）缓曲总长×半径C=切线线路方向左转“-1”，右转“+1“M=第一缓曲“+1”，第二缓曲”-1“N：算边桩左偏“-1”，右偏“+1”D：边桩距离坐标正算数据库例子“肇兴隧道右线数据库”Lbl 0“P”?P:“N”？N：“D”？DIf P≥16836.5 And P≤17009.624:Then 66881.803→A：14635.782→B：16836.5→K:280032030.50→Q:902.044→T:1→C: 1→M:Prog“H1”：IfEnd (第一缓曲线)If P≥17009.624 And P≤17529.681:Then 66914.52→A：14465.781→B：17009.624→K:281035049.370→Q:4700→S：Prog“YY”:IfEnd （圆曲线）If P≥17539.681 And 17702.805≥P:Then 67102.343→A：13799.117→B：17702.805→K:108059031.50→Q:902.044→T:-1→C: -1→M:Prog“H1”：IfEnd (第二缓曲线)If P≥17702.805 And P≤19550:Then 67102.343→A：13799.17→B：17702.805→K:288059031.50→Q: Prog“ZZ”：IfEnd （直线）。

卡西欧5800p双坐标系转换计算程序程序名：ZBZH（坐标转换）“X1=”?A:“Y1=”?B:“X11=”?C:“Y11=”?D:“X2=”?R:“Y2=”?S:“X22=”?U:“Y22=”?V:Pol((R-A),(S-B)):If J≥0:Then J→P:Else J+360°→P:IfEnd↙Pol((U-C),(V-D)):If J≥0:Then J→Q:Else J+360°→Q:IfEnd↙If Q≥P:Then Q-P→M:IfEnd↙If Q＜P:Then 360°-P+Q→M:IfEnd↙Lb1 0:“X00=”?E:“Y00=”?F:Pol((E-C),(F-D)):I→I:J→J:J-M→K:I×cos（K）+A →H:I×sin（K）+B→L↙“JJ=”:M?DMS◢“X0=”:H◢“Y0=”:L◢Goto 0↙说明：以上程序按5800计算器编程规则输入即可,“↙”表示程序内的回车符。

要将坐标在两套坐标系（A整体坐标系和B局部坐标系）内进行转化，必须已知两套坐标系的一条公共直线上的两点分别在两套坐标系内的坐标值，且两坐标系为同手坐标系（均为左手或右手系）。

运行程序后按提示输入数据：X1,Y1为第一个点在A坐标系内的坐标；X11,Y11为第一个点在B坐标系内的坐标；X2,Y2为第二个点在A坐标系内的坐标；X22,Y22为第二个点在B坐标系内的坐标；输完上述步骤后提示输入X00和Y00，此处（X00，Y00）坐标值为B坐标系内的点的坐标值，输入完毕进行计算后可将B坐标系内的（X00,Y00）转换成A坐标系内的（X0,Y0），“JJ”表示两套坐标系的夹角。

eg：在整体坐标系和局部坐标系内有公共线段AB，A在整体坐标系和局部坐标系内的坐标分别为[(1062.4905,367.7365),(0,0)];B在整体坐标系和局部坐标系内的坐标分别为[(1276.1609,890.6324),(461.8813,325.1781)]；求B坐标系内的（-453.2052,487.1620）转化成A坐标系内的坐标。

线路测量坐标正、反算计算原理及卡西欧fx-5800P程序说明一、计算原理在各测量书中对于坐标正算的相关计算式均有说明，故在此不做详解，仅对正算过程中需要用到的原理及公式做一汇总。

对于坐标反算，虽然都采用无限趋近原理进行计算，但计算方式各有差别，本文仅对其中一种自认为相对简单易懂并便于操作的原理进行详解。

1.1 坐标转换[1]如图1，设X P、Y P为P点在国家控制网坐标系中的坐标；x P、y P为P点在工程独立控制网坐标系中的坐标。

X O、Y O为工程独立坐标系原点o在国家坐标系中的坐标，Δα为两坐标系纵坐标轴的交角。

如果一条边在国家坐标系中的坐标方位角为A，而在工程独立坐标系中的坐标方位角为α，则：∆α=A−α(1-1)当由工程独立坐标系坐标换算至国家坐标系坐标时，换算公式为：X=x cos∆α−y sin∆α+X O(1-2)Y=x sin∆α+y cos∆α+Y O}当由国家坐标系坐标换算至工程独立坐标系坐标时，也可使用式(1-2)，此时应将X、Y与x、y互换，且∆α=α−A。

1.2 坐标方位角关系计算1.2.1 正、反坐标方位角[2]一条直线的坐标方位角与直线的前进方向有关，沿直线前进方向的坐标方位角称为正坐标方位角，与其相反方向的坐标方位角称为反坐标方位角。

如图2，由于轴子午线之间是互相平行的，因此同一直线的正、反坐标方位角相差180°，即：α正=α反±180°(2-1)当α反<180°时，取“+”号；当α反>180°时，取“-”号。

1.2.2 坐标方位角的推算[3] 1.2.2.1 转折角为右角如图3(a)，α12为已知边坐标方位角，α23为推算边的坐标方位角，β右为该两边所夹的右角，则：α23=α12±180°−β右=α21−β右 (2-2)1.2.2.2 转折角为左角如图3(b)，α12为已知边坐标方位角，α23为推算边的坐标方位角，β左为该两边所夹的右角，则：α23=α12+β左±180°=α21+β左(2-3) 无论用右角还是左角推算，如遇出现负数的情形，应加上360°。

卡西欧5800计算器坐标正反算程序主程序坐标正反算程序名称：ZBZFSLb1 A↙30→Dim Z :"0=ZS,1=FS"?Z ↙If Z=0:Then Goto B:IfEnd↙ (Z=0进入里程点坐标正算）If Z=1:Then Prog"FSLCZ "↙ (Z=1进入反算里程边桩）Lb1 B ↙Prog”ZBZS” ↙子程序反算里程桩名称：FSLCZLb1 1↙"XK="?X:"YK="?Y↙（输入任意测点的X、Y坐标）“K=”?K ↙（试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）Lb1 2↙Prog"SJK ":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos(W)-(X-Z[18])*Sin(W))→S↙If S<0.0001:Then Goto 4:Else Goto 3:Ifend↙Lb1 3↙K+S→K: Prog"SJK":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z[18])*Sin (W)→Q↙If Q<0.0001 :Then Q→S: Goto 4:Else if Q<S: Then K+Q→K:Goto 2 :Else if Q>S :Then K-Q→K:Goto 2:Ifend:Ifend:Ifend↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙子程序坐标正算名称：ZBZSLb1 0 ↙“XHS="?G(后视点X):"YHS="?L(后视点Y):"XZJ="?M(置镜点X):"YZJ="?N(置镜点Y):Pol(G-M,L-N):"DH=":I(后视距)◢J<0=>J+360→J:"FH=":JDMS◢(后视方位角)Lb1 1↙“K=”?K ：（输入所需计算里程）Prog"SJK "↙XI : Z[18]↙YI: Z[19] ↙Pol(Z[18]-M,Z[19]-N): J<0=>J+360→J↙“PJ=”?P↙ (输入桩与线路夹角)“PD=”?D↙ （输入桩距中线的距离）Z[18]+D*Cos(T+P) →Z[20] ↙Z[19]+D*Sin(T+P) →Z[21] ↙“X=”: Z[20] ◢ （放样坐标X）“Y=”: Z[21] ◢ （放样坐标Y）Pol(Z[20]-M,Z[21]-N):"D=":I◢（放样距）J<0=>J+360→J:"F=":JDMS◢（放样方位角）Goto 1↙子程序数据库名称：SJKif K<本段曲线终点里程And K≥上段曲线终点里程：Then 本段曲线终点里程→Z[1] : 上段曲线终点里程（第一段曲线输起点的里程）→Z[2] ：1→O （注：左偏曲线输入-1→O,右偏曲线输入1→O）: 半径→R :曲线偏角→A:第一缓和曲线→Z[6] : 第二缓和曲线→Z[7] : 交点X→B :交点Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : Pr og”JSPB”:Return: Ifend↙if…………Prog”JSPB”:Retur n:Ifend（曲线段分段输入）↙补充直线段输入如下 (只需输线路的最后一段直线数据）if K<本段直线终点里程And K≥本段直线起点里程:Then 本段直线终点里程→Z[3]:终点X→Z[16]:终点Y→Z[17]:方位角→E: Z[16]+ (K- Z [3])*Cos(E)→Z[18]:Z[17]+ (K- Z[3])*Sin(E)→Z[19] : Return:Ifend ↙子程序计算判别名称：JSPBLb1 2 ↙（曲线要素计算）Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ↙ （M1）Z[7]/2- Z[7]^3/(240*R^2)+ Z[7]^5/(34560*R^4) →Z[9] ↙ （M2）Z[6]^2/(24*R)- Z[6]^4/(2688*R^3) →Z[10] ↙ （P1）Z[7]^2/(24*R)- Z[7]^4/(2688*R^3) →Z[11] ↙ （P2）π*A*R/180+0.5*( Z[6]+ Z[7])→Z[25] ↙ （曲线总长）90* Z[6]/(R*π) →Z[14] ↙ （第一缓和曲线总偏角）90* Z[7]/(R*π) →Z[15] ↙（第二缓和曲线总偏角,可以省略）Z[8]＋（R+Z[10])Tan(A/2)-(Z[10]-Z[11] )/Sin( A)→Z[12]↙ (切线T1)Z[9]＋（R+Z[11])Tan(A/2)+(Z[10]-Z[11] )/Sin (A）→Z[13]↙ (切线T2)B+ Z[12]*Cos (E+180)→ Z[16] ↙ （ZH点X）C+ Z[12]*Sin(E+180)→ Z[17] ↙ （ZH点Y）Z[1]- Z[25]→Z[3] ↙ (ZH点里程)Z[3]+ Z[6]→Z[4] ↙ (HY点里程)Z[1]- Z[7]→Z[5] ↙ (YH点里程)Goto 3 ↙LB1 3 ↙(判断里程点与曲线关系)if K≤Z[3] And K> Z[2] : Then Goto 4 : Ifend ↙if K≤Z[4] And K> Z[3] : Then Goto 5 : Ifend ↙if K≤Z[5] And K> Z[4] : Then Goto 6 : Ifend ↙if K≤Z[1] And K> Z[5] : Then Goto 7 : Ifend ↙Lb1 4 ↙（里程小于直缓点直线独立坐标）K- Z[3] →Z[23] : 0→Z[24] : E→T : Goto 8↙Lb1 5 ↙（第一缓和曲线独立坐标）K- Z[3] →H ↙H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →Z[23] ↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Z[24] ↙90*H^2/( R*π* Z[6]) →T ↙if O>0 ：Then T +E→T : Eles E-T →T : T<0=>360+T→T : Ifend ↙Got o 8 ↙Lb1 6 ↙（圆曲线独立坐标）K- Z[4] →H ↙H*180/( R*π)+ Z[14]→T ↙R*Sin(T)+ Z[8]→Z[23] ↙R*(1-Cos(T))+ Z[10]→Z[24] ↙if O>0 ：Then T +E→T : Eles E-T →T : T<0=>360+T→T : Ifend ↙Goto 8 ↙Lb1 7 ↙（第二缓和曲线独立坐标）Z[1] -K →H ↙H-H^5/(40*R^2* Z[7]^2)+H^9/(3456*R^4* Z[7]^4) →U↙H^3/(6*R* Z[7])-H^7/(336*R^3* Z[7]^3) →V ↙90*H^2/( R*π* Z[7]) →T ↙Z[13]Cos(A)+ Z[12]-U*Cos(A)-V*Sin(A)→Z[23] ↙Z[13]*Sin(A)-U*Sin(A)+V*Cos(A)→Z[24] ↙if O>0 ：Then F-T→T : T<0=>360+T→T : Else F+T →T : Ifend ↙Goto 8 ↙Lb1 8if O<0 : Then - Z[24]→Z[24] : Ifend ↙Z[16]+Z[23]*Cos(E) -Z[24]*Sin（E）→Z[18] ↙Z[17]+Z[23]*Sin(E)+Z[24]*Cos(E)→Z[19] ↙Return↙此程序大约占计算器2900字节！此程序本人已亲自验算，无误！！子程序反算里程桩名称：FSLCZLb1 1↙"XK="?X:"YK="?Y↙（输入任意测点的X、Y坐标）“K=”?K ↙（试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）Lb1 2↙Prog"SJK ":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos(W)-(X-Z[18])* Sin(W))→S↙If S<0.0001:Then Goto 4:Else Goto 3:Ifend↙Lb1 3↙K+S→K: Prog"SJK":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z [18])*Si n (W)→Q↙If Q<0.0001 :Then Q→S: Goto 4:Else if Q<S: Then K+Q →K: Goto 2 :Else if Q>S :Then K-Q→K:Goto 2:Ifend:Ifen d:Ifend↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:------------差个撒样"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙。