5800道路中边桩坐标计算程序

0年11月14日目录一、坐标正算基本公式...............................................................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、SJK1（平面数据库）......................................................14 11、SJK2（纵面数据库）......................................................14 12、SJK3（左路宽度数据库）......................................................15 13、SJK4（右路宽度数据库）......................................................15 14、SJK5（横坡数据库） (16)15、SJK6（下边坡数据库）......................................................16 16、SJK7（左上边坡数据库）......................................................17 17、SJK8（右上边坡数据库）......................................................18 五、后记 (19)0年11月14日CASIO 5800计算器公路工程测量程序一、正算所涉及的计算公式XR id XαβBd Yd ld βI图表 1在图1中，A 点为回旋曲线起点，B 点为回旋曲线止点，I 点为所求坐标点。

附件：道路中桩、边桩计算公式(fx-5800P程序)主程序：“QXZBJS”“QZH”?B: （线路起点桩号，前一个曲线的HZ或YZ,或是）“LS1”?C:“LS2”?I: （第一缓和曲线与第二缓和曲线，可以不等）”JDZH”?N:”JDX”?G:”JDY”?H: (交点桩号与坐标)”XZH”?M: （下一个曲线的ZH或ZY桩号）“T1”?S:”T2”?J: （第一、二切线长）“QXC”?Z: （曲线全长，含缓和曲线长）“JDPJ”?A: （本交点的转向角）“QDFWJ”?O: （起始点到交点的方位角）“R”?R: （本曲线的圆曲线半径）“W”?W: （曲线偏转信息，左为－１，右为＋１）Lbl 0:“K”?K: （所求断面的桩号）If K<B: （比较所求桩号是否小于起点桩号）Then Goto6: （条件为真，转到Lbl 6）IfEnd:If K>N-S:（判定所求桩号是否大于ZH或ZY点，即所求点是否在曲线段）Then Goto 1:（条件为真，转到Lbl 1）IfEnd:(条件为假时运行下例程序，即所求桩号在第一段直线上)G+(N-K)COS(O+180)→X ▲（以上条件都为假是运行该段，所求桩号在直线段，H+(N-K)Sin(O+180)→Y ▲求其中桩坐标）O→T: （起始方位角赋值与T）Prog “ZI-1”Goto 0Lbl 1:If K>N-S+C: （判定所求桩号是否大于HY点）Then Goto 2: (条件为真，转到Lbl 2，即为所求桩在圆曲线或第二缓和曲线上) IfEnd:(条件为假时运行下例程序，即所求桩号在第一缓和曲线上)((K-N+S)2/(6RC))*(180/π) →Q ▲（所求桩中心所占缓和曲线长度的角度）O+3WQ→T:（ZH到所求点的方位角）√(((K-N+S)-(K-N+S)5/(40R2C2))2+(K-N+S)3/(6RC))2)→D:（所求点与ZH点旋长）G+Scos(180+O)+Dcos(O+WQ)→X ▲（根据ZH点坐标和旋长计算中心点坐标）H+Ssin(180+O)+Dsin(O+WQ)→Y ▲Prog “ZI-1”Goto 0Lbl 2:If K>N-S+Z-I:（判定所求桩号是否大于YH点）Then Goto 3: (条件为真，转到Lbl 3，即为所求桩在第二缓和曲线上)IfEnd:(条件为假时运行下例程序，即所求桩号在圆曲线段上)(90(K-N+S-C))/(πR) →Q:（所求占的圆心角的一半，旋的外角）2Rsin(Q) →D:（所求点与HY点之间的旋长）O+W((C/(2R))*(180/π))+2WQ →T:（HY到所求点的方位角）G+Scos(180+O)+ √((C-C5/(40R2C2))2+(C3/(6RC))2)cos(O+W((C2/(6RC))*(180/π)))+Dcos(O+WQ+W((C/(2R))*(180/π)))→X ▲（根据HY点坐标和旋长计算中心点坐标）H+Ssin(180+O)+ √((C-C5/(40R2C2))2+(C3/(6RC))2)sin(O+W((C2/(6RC))*(180/π)))+Dsin(O+WQ+W((C/(2R))*(180/π)))→Y ▲Prog “ZI-1”Goto 0:Lbl 3:If K>N-S+Z:（判定所求桩号是否大于HZ或YZ点）Then Goto 4: (条件为真，转到Lbl 4，即为所求桩在第二段直线上)TfEnd:(条件为假时运行下例程序，即所求桩号在第二缓和曲线段上)((N-S+Z-K)2/(6RI))*(180/π) →Q（所求桩中心所占缓和曲线长度的角度）O+WA-180-3WQ→T（HZ到所求点的方位角）√((N-S+Z-K-( N-S+Z-K)5/(40R2I2))2+( N-S+Z-K)3/(6RI))2) →D（所求点与HZ 点旋长）G+Jcos(O+WA)+Dcos(OWA-180-WQ) →X ▲（根据HZ点坐标和旋长计算中心点坐标）H+Jsin(O+WA)+Dsin(OWA-180-WQ) →Y ▲Prog “ZI-2”:Goto 0:Lbl 4:If K>M:（判定所求桩号是否大于本次计算的桩号范围）Then Goto 6: (条件为真，转到Lbl 6，即为所求桩超出范围)IfEnd:(条件为假时运行下例程序，即所求桩号在第二段直线上)G+Jcos(O+WA)+(K-(N-S+Z))cos(O+WA) →X ▲（中心坐标）H+Jsin(O+WA)+(K-(N-S+Z))sin(O+WA) →Y ▲O+WA→T:Prog “ZI-1”Goto 0:Lbl 6:“END”：子程序“ZI-1”(求边桩坐标)“L”?L:（输入中心至左右边桩的距离L，左为负，右为正）X+Lcos(T+90) →E ▲Y+Lsin(T+90) →F ▲“ZI-2”(求边桩坐标)“L”?L:（输入中心至左右边桩的距离L，左为负，右为正）X+Lcos(T-90) →E ▲Y+Lsin(T-90) →F ▲。

道路桩号算中边桩坐标高程计算程序道路桩号是指道路上的标志桩，用于表示道路上的位置和距离。

在道路规划、设计和施工中，需要根据桩号来确定道路的线形和纵断面，并计算出桩号对应的坐标和高程。

道路桩号的计算程序可以分为以下几个步骤：1.确定基准点：选择一个具备准确坐标和高程的点作为道路的起点，确定其坐标和高程。

2.确定桩号起点：确定一个参考点作为桩号的起点，通常选择道路的起点或其他规定的地点。

为了方便计算，可以选择一个整数作为起点桩号，如0、100等。

3.桩号计算：根据道路设计和实际情况，确定桩号的计数方式和间隔。

通常情况下，桩号以米为单位，从起点开始递增或递减。

4.桩号与坐标的关系：桩号与坐标之间存在一定的数学关系，可以根据道路的几何特征和设计参数进行计算。

例如，对于一条平直无坡道路，可以使用线性插值法计算桩号对应的坐标。

5.桩号与高程的关系：桩号与高程之间也存在一定的数学关系，可以根据道路的纵断面和地形特征进行计算。

例如，对于一条按规定坡度设计的道路，可以使用坡比法计算桩号对应的高程。

6.精度控制：在桩号计算过程中，需要考虑测量误差和计算方法的精度。

为了提高计算结果的准确性，可以采用较精确的测量方法和计算算法，并进行误差修正。

7.应用场景：道路桩号的计算程序可以应用于道路工程中的位置控制、导线布设、测量定位、横断面绘制等方面，为道路规划、建设和维护提供准确的空间位置和高程信息。

总结起来，道路桩号的计算程序是根据道路的设计和实际情况，通过选择基准点和起点桩号，确定桩号计算方式和间隔，以及桩号与坐标、高程之间的关系，计算出桩号对应的坐标和高程。

这个程序可以应用于道路工程中的各个环节，为道路的设计、施工和维护提供准确的空间位置和高程信息，提高工程质量和效率。

本程序非本人原创，结合上传者长处，进行改编，愿收获和大家共享，不足之处多多指教。

希望过路神仙把竖曲线加进去带算隧道超欠挖，谢谢曲線任意里程中邊樁坐標正反算5800p本程序由内蒙@文哥和师父阿玖哥共同改编完成主程序 ZBJSLbl 6："1.SZ→ XY"："2.XY→ SZ"：？S：？N：Prog“SUB0”1÷P→C：(P-R)÷(2HPR) →D：180÷π→E:N=1 =>Goto 1: Goto 2Lbl 1:?Z: Abs(S-O)→W: Prog"SUB1":"XS=":X◢"YS=":Y◢"FS=":F-90 ►DMS◢“XC” ?C：“YC” ? E：0→I：0→J：Pol（X-C，Y-E）:“I=”：I◢If J∠0: Then “J=”: J+360►DMS◢IfEndIf J >0：Then “J=” ：J►DMS◢IfEnd：Goto 6Lbl 2?X:?Y:X→I:Y→J:Prog"SUB2":O+W→S:Cls:"S=":S◢"Z=":Z◢Goto 6正算子程序SUB10.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F：1-K→M：U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bc os(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD)))→X：V+W(Asin(G+QEKW(C+KWD))+Bsi n(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QE MW(C+MWD))) →Y：G+QEW(C+WD)+90→F：X+Zcos（F）→X：Y+Zsin（F）→Y反算子程序 SUB2G-90→T(Y-V)cosT-(X-U)sin(T) →WAbs(W)→W:0→ZLbl6：Prog "SUB1"T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→ZIF Abs（Z）<1E-6：Then0→Z:Prog "SUB1":(J-Y)÷sin（F）→Z:Else W+Z→W：G oto6:IfEnd数据库子程序SUB0Goto 1(线元可输入多条，分离式可在前多加一位，匝道一样。

5800道路坐标（正反算）程序（终极版）5800道路坐标（正反算）程序1. ZBJSCX (主程序)Lbl 4："ZHENG SUAN=1 "："FAN SUAN=2"：？N：“ZH”？S：Prog“SJK0”↙1÷P→C: (P-R)÷(2HPR) →D:180÷∏→E:N=1 => Goto1： Goto2：↙Lbl 1:“BD”? Z： Abs(S-O) →W：Prog "ZS"： "X="：X◢"Y="：Y◢F-90→F：“FWJ=”：F▲DMS◢Goto4↙Lbl 2：？X：？Y：X→I：Y→J：Prog“FS”：O+W→S：“ZH=”：S◢“BD=”：Z◢Goto4↙2. ZS (正算子程序)0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F： 1-K→M：U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))) →X：V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin( G+QEMW(C+MWD))) →Y：G+QEW(C+WD)+90→F： X+Zcos（F）→X： Y+Zsin（F）→Y3. FS (反算子程序)G-90→T： Abs((Y-V)cosT-(X-U)sinT) →W：0→Z：Lbl 0：Prog "ZS"：T+QEW(C+WD) →L： (J-Y)cos（L）-(I-X)sin（L）→Z：IF Abs（Z）<0.000001：ThenGoto1:ElseW+Z→W:IfEnd:Goto0↙Lbl1:0→Z:Prog "ZS": (J-Y)÷sin（F）→Z4,SJK0( 数据库子程序)Goto1↙ (同时保存多个曲线时的指针)Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd↙Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1（直线0、左偏-1、右偏1）→Q：Return:IfEnd……………..为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。

Casio5800交点法程序（输入版）（歪哥哥2009版）本程序由一个主程序JD和三个子程序（JDA、JDB、JDC）构成，运行时只需运行主程序即可！本程序适用于单交点对称型、不对称型、无缓和曲线单圆曲线型一个交点范围内（含交点前后有直线段时）的曲线要素核对和坐标计算，手工输入要素，对设计图纸的“直线、曲线转角表”中交点数据进行复核验证，并为线元法程序提供起点坐标起点切线方位角等数据！当然本程序也可单独逐交点输入进行放样计算用！鉴于5800计算器的空间和以上所述本程序的主要目的，故此程序不修改为数据库版本！需要的自行修改结合XY框架自己修改为数据库反算程序等！主程序名：5.JD24→Dimz↙Cls :＂XC＂?U :＂YC＂?V :＂K（JD）＂?K :＂X（JD）＂?X :＂Y（JD）＂?Y :＂LS1＂?B :＂LS2＂?C : ?R :＂（ZH）FWJ°＂?M : ＂α（Z-,Y+）°＂?O : M+O→N :Prog ＂JDA＂↙Cls :＂T1=＂:＂T2=＂:＂L=＂:＂L Y=＂: Locate 4,1,S : Locate 4,2,T : Locate 4,3,L : Locate 4,4,Q◢Cls :＂E=＂:＂K（ZH）=＂: Locate 7,1,E : Locate 7,2,Z[1]◢Cls : ＂K（HY）=＂:＂K（QZ）=＂:＂K（YH）=＂:＂K（HZ）=＂: Locate 7,1, Z[2] : Locate 7,2, Z[3] : Locate 7,3, Z[4] : Locate 7,4, Z[5]◢LbI 0 : ＂K×+×××＂?P : ＂Z＂?D : If D≠0 :Then ＂RJ＂?H : IfEnd : Prog ＂JDB＂↙If D＜0 :Then Cls : ＂X(L)=＂:＂Y(L)=＂: Locate 6,1,F : Locate 6,2,G◢Pol(F-U,G-V : Cls : ＂S(L)=＂: Locate 6,1,I : ＂F(L)=＂:360Frac((J+360)÷360▼DMS◢Goto 0 : IfEnd↙If D=0 :Then Cls : ＂X(Z)=＂:＂Y(Z)=＂: Locate 6,1,F : Locate 6,2,G : ＂QXFWJ(Z)=＂: Z▼DMS◢Pol(F-U,G-V : Cls : ＂S(Z)=＂: Locate 6,1,I : ＂F(Z)=＂:360Frac((J+360)÷360▼DMS◢Goto 0 : IfEnd↙If D＞0 :Then Cls : ＂X(R)=＂:＂Y(R)=＂: Locate 6,1,F : Locate 6,2,G◢Pol(F-U,G-V : Cls : ＂S(R)=＂: Locate 6,1,I : ＂F(R)=＂:360Frac((J+360)÷360▼DMS◢Goto 0 : IfEnd↙子程序1名：JDAIf O＜0 :Then -1→W : Else 1→W : IfEnd : WO→A ↙B2 ÷24÷R-B^(4)÷2688÷R ^(3) →Z[6] ↙C2 ÷24÷R-C^(4)÷2688÷R ^(3) →Z[7] ↙B÷2-B^(3)÷240÷R2 →Z[8] ↙C÷2-C^(3)÷240÷R2 →Z[9] ↙Z[8]+((R+Z[7]-(R+Z[6])cos(A))÷sin(A))→S↙Z[9]+((R+Z[6]-(R+Z[7])cos(A))÷sin(A))→T↙RAπ÷180+(B+C) ÷2→L↙RAπ÷180-(B+C) ÷2→Q↙(R+(Z[6]+Z[7])÷2)÷cos(A÷2)-R→E↙K-S→Z[1] ↙↙Z[1]+B→Z[2] ↙↙Z[2]+Q÷2→Z[3]↙Z[1]+L-C→Z[4]↙Z[4]+C→Z[5]↙子程序2名：JDBX-Scos(M)→Z[19]:Y-Ssin(M)→Z[20]↙X+Tcos(N)→Z[21]:Y+Tsin(N)→Z[22]↙If P＞Z[1]:Then Goto 1 :IfEnd↙Z[1]-P→L↙X-(S+L)cos(M)+Dcos(Z+H)→F↙Y-(S+L)sin(M)+Dsin(Z+H)→G↙M→Z : Goto 5↙LbI 1 : If P＞Z[2]:Then Goto 2 :IfEnd↙P-Z[1]→L:L→Z[12]:B→Z[13]:Prog＂JDC＂↙Z[19]+Z[14]cos(M)-WZ[15]sin(M)+Dcos(Z+H)→F↙Z[20]+Z[14]sin(M)+WZ[15]cos(M)+Dsin(Z+H)→G↙M+90WL2 ÷(BRπ)→Z↙Goto 5↙LbI 2 : If P＞Z[4]:Then Goto 3 :IfEnd↙P-Z[1]→L:90(2L-B)÷R÷π→Z[11]↙Rsin(Z[11])+Z[8]→Z[14]:R(1-cos(Z[11]))+Z[6]→Z[15]↙Z[19]+Z[14]cos(M)-WZ[15]sin(M)+Dcos(Z+H)→F↙Z[20]+Z[14]sin(M)+WZ[15]cos(M)+Dsin(Z+H)→G↙M+WZ[11]→Z↙Goto 5↙LbI 3 : If P＞Z[5]:Then Goto 4 :IfEnd↙Z[5]-P→L:L→Z[12]:C→Z[13]:Prog＂JDC＂↙Z[21]-Z[14]cos(N)-WZ[15]sin(N)+Dcos(Z+H)→F↙Z[22]-Z[14]sin(N)+WZ[15]cos(N)+Dsin(Z+H)→G↙N-90WL2 ÷(CRπ)→Z↙Goto 5↙LbI 4 : P-Z[5]→L↙X+(T+L)cos(N)+Dcos(Z+H)→F↙Y+(T+L)sin(N)+Dsin(Z+H)→G↙N→Z↙Goto 5↙LbI 5 : 360Frac((Z+360)÷360→Z↙子程序3名：JDCIf Z[12]=0 :Then 0→Z[14]: 0→Z[15]:Else↙Z[12]- Z[12]^(5)÷40÷(RZ[13])2+ Z[12]^(9)÷3456÷(RZ[13])^(4) →Z[14]↙Z[12]^(3)÷6÷(RZ[13])-Z[12]^(7)÷336÷(RZ[13])^(3)+ Z[12]^(11) ÷42240÷(RZ[13])^(5)→Z[15] ↙IfEnd↙程序说明：已知数据输入：XC ? 测站X坐标YC ? 测站Y坐标K（JD）?交点桩号X（JD）?交点X坐标Y（JD）?交点Y坐标LS1 ?第一缓和曲线长度LS2 ?第二缓和曲线长度R ? 圆曲线半径（ZH）FWJ°?交点前（即前交点至本交点也即ZH点）的正切线方位角α（Z-,Y+）?本交点处线路转角（左转为负，右转为正，度分秒输入）K×+×××? 待求桩号Z ?待求桩号距中距离（左负值，右正值，中为0）RJ ?斜交右角（线路切线前进方向与边桩右侧夹角）计算结果显示：T1=第一切线长T2=第二切线长L=曲线总长L Y=圆曲线长E=曲线外距K(ZH)=直缓点桩号K(HY)=缓圆点桩号K(QZ)=曲中点桩号K(YH)=圆缓点桩号K(HZ)=缓直点桩号X= Y=待求点的坐标（其中：L-左Z-中R-右）QXFWJ(Z)=待求点的中桩切线方位角（当求中桩坐标时显示）S= F=测站至待求点的水平距离、方位角（其中L-左Z-中R-右）（歪哥哥郑重声明：代码无误，输入注意！使用者使用不当所造成的一切损失概不负责！！）。

交点法解算线路中边桩坐标5800程序摘要：采用本程序，中边桩计算只需三个交点坐标及JD里程，边桩计算只需输入里程及离对应中桩的距离即可，程序内容本身使用语言种类少，逻辑运行简单。

关键词：交点法；坐标；程序；卡西欧5800在线路施工中，全站仪配合卡西欧5800计算机放样，仍然是目前采用的主要方法之一。

现以单一交点为例，介绍线路中边桩坐标计算的方法。

线路平曲线设计数据：交点坐标：JD1（74586.235，45658.223），JD2（75482.321，45586.235），JD3（75658.243，45422.321），JD2里程K1+156.365,曲线半径1000m，缓和曲线长度60m。

主程序命名为“ZBJS”,输入设计参数后程序内容为：100→DimZ：74586.235→Z[1]:45658.223→Z[2]:75482.321→Z[3]:45586.235→Z[4]:75658.243→Z[5]: 45422.321→Z[6]:1156.365→Z[70]:1000→R:60→L:Prog“JDF”：Prog“K1”Prog“K1”代表线路为左偏曲线，若线路为右偏曲线，上方最后语句修改为Prog“K2”。

Prog“JDF”程序内容为：Pol(Z[3]-Z[1]，Z[4]-Z[2]）:I→Z[7]:J→Z[8]:Pol(Z[5]-Z[3]，Z[6]-Z[4]）:I→Z[9]:J→Z[10]:Abs（Z[8]-Z[10]）→Z[11]:L2÷24÷R÷2688÷R＾(3)→P:L÷L＾(3)÷240÷R2→M:90L÷π÷R→Z[12]:（R+P）tan（Z[11]÷2）+M→T：πR(Z[11]-2Z[12]）÷180→Z[13]:Z[70]-T→Z[14]:Z[13]+Z[14]+L→Z[15]:Z[3]+Tcos(Z[8]+180°）→Z[21]:Z[4]+Tsin(Z[8]+180°）→Z[22]:Z[3]+Tcos(Z[10]）→Z[23]:Z[4]+Tsin(Z[10]）→Z[24]:Cls:“A1=”:Locate4,1,Z[8]:“A2=”:Locate4,2,Z[10]:“ZH=”:Locate4,3,Z[14]:“YH=”:L ocate4,4,Z[15]◢Cls:“XZH=”:Locate5,1,Z[21]:“YZH=”:Locate5,2,Z[22]:“XHZ=”:Locate5,3,Z[23]:“YHZ=”:Locate5,4,Z[24]◢Prog“K1”程序内容为：“K?”K:Z[8]→A:IfK≤Z[14]:ThenZ[1]→N:Z[2]→E:K-Z[70]+Z[7]→D:Prog“ZX”:ElseZ[21]→N:Z[22]→E:IfK≤Z[14]+L:ThenK-Z[14]→D:Prog“H1Z”:ElseIfK≤Z[15]:ThenK-Z[14]→D:Prog“YZ”:ElseZ[23]→N:Z[24]→E:IfK≤Z[15]+L:ThenZ[10]+180°→A:Z[15]+L-K→D:Prog“H2Z”:ElseZ[10]→A:K-Z[15]-L→D:Prog“ZX”Prog“K2”程序内容为：“K?”K:Z[8]→A:IfK≤Z[14]:ThenZ[1]→N:Z[2]→E:K-Z[70]+Z[7]→D:Prog“ZX”:ElseZ[21]→N:Z[22]→E:IfK≤Z[14]+L:T henK-Z[14]→D:Prog“H1Y”:ElseIfK≤Z[15]:ThenK-Z[14]→D:Prog“YY”:ElseZ[23]→N:Z[24]→E:IfK≤Z[15]+L:ThenZ[10]+180°→A:Z[15]+L-K→D:Prog“H2Y”:ElseZ[10]→A:K-Z[15]-L→D:Prog“ZX”Prog“K1”及Prog“K2”中的子程序在作者所著写的论文《积木法编制卡西欧5800线路坐标计算程序》中皆有详细编写。

5800计算器程序输入格式YOUXIAN―K=‖ ? L: If L <= 108451.609 : Then Prog ―GT –BQU‖ : IfEnd 回车换行ZUOXIAN―K=‖ ? L: If L <= 391200 : Then Prog ―GT –AQU‖ : IfEnd 回车换行一、主程序（ALLINE）―K=‖ ? L: If L <= 240987.601 : Then Prog ―GT – AQU‖ : IfEnd 回车换行If L <= 243536.746 : Then Prog ―GT – BQU‖ : IfEnd二、1--1级子程序(GT – AQU)001（此为编号，不用输入，下同）If L<= 228634.196 : Then 227400 → A : 228634.196 → B : 0→R : 0→V :2978464.863 →C :487166.400 →D :129-01-13.17 →E（方位角采用度分秒的格式输入）:Prog‖H-SIMP‖ : Return :IfEnd回车换行002 If L<=228964.196 : Then 228634.196 → A : 228964.196 → B : 0 →R: -1 /5000 →V :2977687.820 →C :488125.276 →D : 129—01—13.17 →E（方位角采用度分秒的格式输入）:PROG‖H-SIMP‖ : Return :Ifend 换行符003 If L <= 233231.971 : THEN 228964.196 → A : 233231.971 → B : -1/5000 →R : -1 /5000 →V :2977482.896 →C :488383.918→D :127—07—46.59 →E（方位角采用度分秒的格式输入）:PROG‖H-SIMP‖ : Return :Ifend 换行符004If L<= 233561.971 : THEN 233231.971 →A : 233561.971 →B : -1/5000→R:0→V:2976574.472 →C :492422.451 →D :78-13-28.23 →E （方位角采用度分秒的格式输入）:Prog‖H-SIMP‖ : Return :IfEnd 换行符005IF L<= 239146.929 : THEN 233561.971 → A : 239146.929→ B : 0 →R:0 →V :2976648.916 →C :492743.884 →D :76-20-01.33 →E（方位角采用度分秒的格式输入）:PROG‖H-SIMP‖ : Return :Ifend 换行符006IF L<= 239316.929 : THEN 239146.929 → A : 239316.929 → B : 0 →R:-1/8000 →V :2977822.207→C :497569.249 →D :76-20-01.33→E （方位角采用度分秒的格式输入）:PROG‖H-SIMP‖: Return :Ifend 换行符007IF L<= 239700.272 : THEN 239316.929 → A : 239700.272 →B:-1/8000→R:-1/8000→V:2977862.957→C:497734.292→D :72-43-29.7 7→E（方位角采用度分秒的格式输入）:PROG‖H-SIMP‖: Return :Ifend 换行符008IF L<= 239870.272 : THEN 239700.272 →A : 239870.272 →B:-1/8000→R:0→V:2977966.344→C:498103.392→D :72-58-46→E（方位角采用度分秒的格式输入）:PROG‖H-SIMP‖ : Return :Ifend 换行符009 IF L<= 240987.601 : THEN 239870.272 →A : 240987.601 →B:0→R:0→V:2978017.255→C:498265.587→D :72-22-14.44→E（方位角采用度分秒的格式输入）:PROG‖H-SIMP‖ : Return :Ifend 换行符注：L 代表桩号，A 代表这一段的起点桩号，B代表终点桩号，R代表起点曲率，V代表终点曲率（左偏为—，右偏为+），C代表起点坐标x，D代表起点坐标y，E 代表起点方位角。