线路测量坐标正、反算计算原理及卡西欧5800程序说明
CASIO5800坐标正反算计算程序
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列表数组中，要想多次调用最好随公式取出结果，并赋给变量。
5800计算器全线坐标计算放样正反算程序
5800计算器全线坐标计算放样正反算程序此处无法提供完整的程序，但是可以提供一些关于全线坐标计算放样正反算的基本知识。
全线坐标计算放样是指将一条曲线按照一定比例进行分割，然后计算每个分段的坐标值，最终得到完整曲线上每个点的坐标。
这个过程可以分为正算和反算两部分。
正算是根据曲线的起点、终点以及每个分段的长度，通过数学公式计算出每个分段末点的坐标，最终得到完整曲线上每个点的坐标。
反算则是根据给定的曲线上某个点的坐标，以及曲线起点到该点的距离，计算出该点在曲线上的位置。
在实际工作中，使用计算器进行全线坐标计算放样非常常见。
通常使用的计算器有5800和2000等型号，这些计算器都有相应的程序可以使用。
以下是大致的步骤：
正算步骤：
1. 选择“坐标计算放样”功能；
2. 输入曲线的起点坐标和终点坐标；
3. 输入曲线分段长度以及分段数量；
4. 计算每个分段的末点坐标；
5. 输出所有点的坐标。
反算步骤：
1. 选择“坐标计算放样”功能；
2. 输入曲线的起点坐标和终点坐标；
3. 输入参考点的坐标；
4. 输入参考点到曲线起点的距离；
5. 计算参考点在曲线上的位置；
6. 输出参考点在曲线上的坐标。
全线坐标计算放样是一项非常重要的工程测量工作，也是计算器的常见应用之一。
需要结合具体需求和计算器机型进行操作。
卡西欧5800测量坐标计算程序
一．使用说明：1.规定：（1）以道路中线的前进方向区分线路的左右偏转方向，线路左偏Q=-1，线路右偏Q=1，直线Q=0。
（2）当所求点位于线路中线时，Z=0；当位于左侧时，Z取负值；当位于右侧时，Z取正值。
（3）当线元为直线时，其起终点的曲率半径为无穷大，取10的45次方代替。
（4）当线元为圆曲线时，起终点曲率半径为圆曲线半径。
（5）当线元为缓和曲线时，起终点曲率半径与直线相接取10的45次方代替，与圆曲线相接取圆曲线半径。
2.输入与显示说明：（1）输入部分：N？：选择计算方式，输入1表示由里程边距计算坐标，输入2表示由坐标计算里程边距。
X0？线元起点X坐标。
Y0？线元起点Y坐标。
S0？线元起点里程。
F0？线元起点切线方位角。
LS？线元长度。
R0？线元起点曲率半径。
RN?线元终点曲率半径。
Q？线元左右偏标志（左偏=-1，右偏=1，直线=0）S?正算时所求点里程。
Z？正算时所求点距中线边距。
位于线路左侧输入负值，在右侧输入正值。
AR?斜交角度。
（为线路右角，在0-180度之间）X?反算时所求点X坐标。
Y？反算时所求点Y坐标。
（2）显示部分：XS=计算结果，所求点X坐标。
YS=计算结果，所求点Y坐标。
FS=计算结果，所求点切线方位角。
S=计算结果，反算时所求点里程Z=计算结果，反算时所求点边距。
二．程序清单1（用于任意曲线元的坐标计算）1.主程序(TYQXJS)″1.SZ=>XY″:″2.XY=>SZ″:?N:″X0″?U:″Y0″?V:″S0″?O: ″F0″?G:″LS″?H:″R0″?P:″RN″?R:?Q:1÷P→C:(P-R) ÷(2HPR) →D:180÷π→E:If N=1: THEN Goto 1:Else Goto 2:IfEnd:Lbl 1:?S:”D(Z-，Y+)=”?Z:90→T:”O=”?T:Abs(S-O) →W:Prog ″SUB1″:″XS″:X◢″YS″:Y◢″FS″:F◢Goto 1:Lbl 2:?X:?Y:I→X:J→Y:Prog ″SUB2″:″S″:O+W→S:◢″Z″:Z◢Goto 22.正算子程序(SUB1)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) →F:X+Zcos(F+T) →X:Y+Zsin(F+T) →Y3.反算子程序(SUB2)G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z:Lbl 0:Prog ″SUB1″:T+QEW(C+WD) →L:(J-Y)cos(L)-(I-X)sin(L)→Z:IF Abs(Z)<1E-6 :Then Goto 1:Else W+Z→W:Goto 0 : IfEnd:Lbl 1:0→Z:Prog ″SUB1″:(J-Y)÷sin(F)→Z 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。
卡西欧5800坐标正反算万能程序
CASIO5800计算器程序曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序一、程序功能本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。
另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。
本程序也可以CASIO fx-5800P计算器计算器上运行。
二、源程序1.主程序(TYQXJS)"1.SZ => XY"："2.XY => SZ"：?N："X0"? U："Y0"? V："S0"? O："F0"? G："LS"? H："R0"? P："RN"? R：?Q←┘1÷P-> C：(P-R)÷(2HPR) -> D：180÷π->EIf N=1：Then Goto1：Else Goto 2：IfEndLbl 1：？S：？Z：Abs(S-O) -> W：Prog "SUB1"："XS"：X◢"YS"：Y◢Goto 1←┘Lbl 2：？X：？Y：X-> I：Y->J：Prog "SUB2"："S"：O+W-> S◢"Z"：Z◢Goto 22. 正算子程序(SUB1)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+K WD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD)))-> Y：G+QEW(C+WD)+90-> F：X+ZcosF-> X：Y+ZsinF-> Y3. 反算子程序(SUB2)G-90-> T：Abs((Y-V)cosT-(X-U)sinT) -> W ：0-> Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD) ->L：(J-Y)cosL-(I-X)sinL-> Z：If Abs(Z) <0.0000001：Then Goto1：Else W+Z ->W：Goto 0：IfEndLbl 1：0 ->Z：Prog "SUB1"：(J-Y)÷sinF-> Z三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
卡西欧fx-5800p计算器测量程序
卡西欧fx-5800p计算器测量程序（作者：水江华）1.坐标反算程序"ZBFS QH1-4"◢Deg：Fix 3←┘"X+YiS(m）="?A←┘Lb1 0:"X+YiE(m)，0=>END="?B←┘Abs（B）=0=>Goto E←┘Arg（B-A)→JJ<0=>J+360→J"DIST(m）="：Abs（B-A)◢"α（DMS)=":J▶DMS◢Goto 0←┘Lb1 E:"ZBFS QH1-4=>END"2.坐标正算程序"ZBZS QH1-3"◢Deg：Fix 3←┘"XS(m)= "？A:"YS(m)= "？B←┘Lb1 1:"DIST(m)= "？L←┘If L≤0:Then Goto 2:IfEnd←┘"α（Deg)="？R←┘A+Lcos（R）→C←┘B+Lsin（R）→D←┘Fix 3："XE(m)= ":C◢"YE(m)= ":D◢Goto 1←┘Lb1 2:"ZBZS QH1-3=>END"3.单一闭、附合导线近似平差程序"BFHDXPC QH1-8"◢Deg:Norm 1:FreqOn←┘n→N←┘If List Y[N]=0:Then N-1→DimZ:Goto 0←┘Else N→DimZ:Goto 1: IfEnd←┘Lb1 0: "CLOSE Or CONECT TRA VE"◢"ANGLE NUM=":N◢"SIDE NUM=":N-1◢"UNKNOWN P NUM=":N-2◢"CLOSE(0)，CONNECT(≠0)="?→Z←┘Fix 3："XA(m)，≤0=>αA→B(Deg)="？A←┘If A>0：Then "YA(m)= "？B：Else "αA→B(Deg)="？R：IfEnd←┘"XB(m)= "？C:"YB(m)= "？D←┘C+Di→U：U→V←┘If A>0：Then Pol（C-A，D-B）：Cls←┘J<0=>J+360→J：J→R←┘"DIST A→B(m)="：I◢"αA→B(DMS)="：R▶DMS◢IfEnd←┘If Z≠0：Then "XC(m)= "？E:"YC(m)= "？F←┘E+Fi→V←┘"XD(m)，≤0=>αC→D(Deg)="？G←┘If G>0：Then "YD(m)= "？H←┘Pol（G-E，H-F）：Cls←┘J<0=>J+360→J：J→S←┘"DIST C→D(m)="：I◢"αC→D(DMS)="：S▶DMS◢Else "αC→D(Deg)="？S：IfEnd←┘Else If R>180：Then R-180→S：Else R+180→S：IfEnd←┘IfEnd←┘∑y→M←┘For 1→I To N←┘If I=1：Then R+List X[I]→L：Else L+List X[I]→L：IfEnd←┘If L>180：Then L-180→L：Else L+180→L：IfEnd←┘Nent←┘3600(L-S)→T←┘"α CLOSE ERROR(S)="：T◢40→W←┘Abs(T)>W=>"α CLOSE ERROR OVRE!"←┘-T÷N÷3600→P：0→O←┘For 1→I To N←┘If I=1：Then R+List X[I]+P→L←┘Else L+List X[I] +P→L：IfEnd←┘If L>180：Then L-180→L：Else L+180→L：IfEnd←┘If I<N：Then List Y[I]<L→Z[I] ：Z[I]+O→O：IfEnd←┘Nent←┘3600(L-S)→T←┘"CHECK α CLOSE ERROR(S)="：T◢U+O-V→O←┘Int(M÷Abs(O))→K←┘"DELTA X(m)= "：ReP(O)◢"DELTA Y(m)= "：ImP(O)◢"RELAT CLOSE ERROR=1÷"：K◢K<4000=>"RELAT CLOSE ERROR OVER! "←┘For 1→I To N-1←┘Z[I]-OList Y[I]÷M→X←┘If I=1：Then U+X→Y：Else Y+X→Y：IfEnd←┘Norm 1："POINT n="：I◢Fix 3："X(m)= "：ReP(Y)◢"Y(m)= "：ImP(Y)◢Nent←┘Y-V→Q←┘"CHECK X (m)= "：ReP(Q)◢"CHECK Y (m)= "：ImP(Q)◢Goto E←┘Lb1 1："NO DIRECT TRA VE"◢"ANGLE NUM="：N-1◢"SIDE NUM="：N◢"UNKNOWN P NUM="：N-1◢"XA(m)= "？A:"YA(m)= "？B←┘"XB(m)= "？C:"YB(m)= "？D←┘A+Bi→U:C+Di→V←┘0→H:U+List Y[1]<H→Z[1]←┘For 1→I To N-1←┘H+List X[1]→H←┘If H>180：Then H-180→H：Else H+180→L：IfEnd←┘以下未完程序由读者完成。
卡西欧5800公路坐标正反算程序
卡西欧5800公路坐标正反算程序卡西欧5800公路坐标正反算程序是一种用于测量和计算公路坐标的工具。
它基于卡西欧5800全站仪的技术和功能，能够准确地测量和计算公路的坐标信息。
下面我将详细介绍卡西欧5800公路坐标正反算程序的工作原理和使用方法。
一、工作原理：卡西欧5800全站仪是一种高精度的测量仪器，它能够通过测量角度和距离来确定地点的坐标。
在公路测量中，我们通常会在公路上设置一系列的控制点，然后使用全站仪测量这些控制点的坐标。
通过这些测量数据，我们可以计算出公路上其他点的坐标。
卡西欧5800公路坐标正反算程序基于这些测量数据，通过一系列的计算和算法，能够准确地计算出公路上其他点的坐标。
具体的计算方法包括正算和反算。
正算是指已知控制点的坐标，通过测量数据和算法计算出其他点的坐标。
反算是指已知某个点的坐标，通过测量数据和算法计算出该点在公路上的位置。
二、使用方法：1. 设置控制点：首先，在公路上选择一些具有代表性的控制点，并使用全站仪进行测量。
在测量过程中，需要确保全站仪的水平和垂直仪器准确，以获得准确的测量数据。
2. 导入测量数据：将测量数据导入卡西欧5800公路坐标正反算程序。
该程序通常会提供一个界面，用于导入和管理测量数据。
3. 正算：如果需要计算公路上其他点的坐标，可以选择正算功能。
在正算功能中，需要输入已知控制点的坐标和测量数据。
程序将根据这些数据进行计算，并给出其他点的坐标。
4. 反算：如果需要计算某个点在公路上的位置，可以选择反算功能。
在反算功能中，需要输入已知点的坐标和测量数据。
程序将根据这些数据进行计算，并给出该点在公路上的位置。
5. 结果输出：计算完成后，卡西欧5800公路坐标正反算程序将给出计算结果。
结果通常以表格或图形的形式呈现，以便用户查看和使用。
三、注意事项：1. 测量准确性：为了获得准确的测量结果，使用卡西欧5800全站仪时需要注意测量准确性。
在测量过程中，需要保持仪器的稳定和准确，并避免测量误差。
卡西欧fx-5800 坐标正算程序
卡西欧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 （直线）。
CASIO5800ZXZFS直线正反算坐标计算程序
CASIO5800ZXZFS直线正反算坐标计算程序CASIO fx-5800P ZXZFS 直线正反算坐标计算程序原理：采用正余弦定理求待求点垂距进行确定里程和偏距。
公式：tan-1((yb-ya)/(xb-xa)) 、√(xB-xA)2+(yB-yA)2、sinα×D、Q÷tan(α)说明：该程序适用于计算器 CASIO fx-5800P，针对直线正反算设计，正算可算直线中边桩坐标及其斜交、反算可推算直线距线路中心的里程偏距。
1、DK(QD)? 输入起点里程2、X(QD)? 输入起点坐标X3、Y(QD)? 输入起点坐标Y4、FWJ? 输入计算方位角5、J? 输入1为正算、输入2为反算6、JSDK? 输入计算里程主程序名：ZXZFSDeg：Fix 4↙＂DK(QD)＂?A：＂X(QD)＂？B：＂Y(QD)＂？C：＂FWJ＂？F↙Lbl 0＂1=>XY.2=>DK＂：＂J＂？W↙If W=1：Then Goto 1：Else If W=2：Then Goto 2：IfEnd↙Lb1 1↙＂JSDK＂？E：＂PJ＂?G：＂PA＂？H↙＂X=＂：B+(E-A) ×cos(F)+G×cos(F+H)→X◢＂Y=＂：C+(E-A) ×sin(F)+G×sin(F+H)→Y◢＂ZJD(X)＂？S：＂ZJD(Y)＂？T↙＂PJ=＂：Pol(X-S，Y-T)：I→M◢＂FWJ=＂：If J<0：Then J+360→N：Else J→N：N▲DMS◢Goto 0↙Lb1 2↙＂X=＂？K：＂Y=＂？L↙Pol(K-B，L-C)：I→D：If J<0：Then J+360→P：Else J→P：IfEnd↙If F<="" p="" p-f→ｏ：else="">Sin(Ｏ)×D→Q↙＂DK=＂：A+Q÷tan(Ｏ)→R◢＂PJ=＂：sin(Ｏ)×D→Q：If F<="" 1×q→q：else="" p=""> Goto 0注：1、◢ 为显示指令。
线路测量坐标正、反算计算原理及卡西欧5800程序说明
线路测量坐标正、反算计算原理及卡西欧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线路坐标正反算程序
卡西欧5800线路坐标正反算程序（测试篇）（ZBZFS）主程序：Lbl A↙Cls ：Fix 4 ：30→Dim Z：“0=ZS,1=FS”?Z↙If Z=0：Then Goto B：IfEnd↙If Z=1：Then Prog“FSLCBZ”：Goto A：IfEnd↙Lbl B↙“K=”?K：Prog“SJL1”：Prog“ZBFS”：GotoA↙子程序（SJK1）If K＜HZ桩号And K≥上段曲线HZ桩号：Then HZ→Z[1]：上段曲线HZ桩号→Z[2]：“1=Y，-1=Z”→O：半径→R：转向角→A：L1→Z[6]：L2→Z[7]：JDx→B：JDy→C：前切线方位角→E：后切线方位角→F：Prog“XLZBJSCX”：Return：IfEnd↙If K＜Z[1]：then Prog”XLZBJSCX”：Return：ifend↙If K≥Z[2]：then Prog”XLZBJSCX”Return：ifend↙If K＜（本段直线终点桩号）And K≥（本段直线起点桩号）：Then (本段直线终点桩号)→Z[3]：本段直线终点X→Z[16]：本段直线终点Y→Z[17]：方位角→E：Z[16]+(K-Z[3])×cosE→Z[18]：Z[17]+(K-Z[3])×sinE→Z[19]：Return：ifend↙子程序（XLZBJSCX）Lbl 2↙Z[6]÷2-Z[6]3÷（240×R2）+Z[6]5÷（34560×R4）→Z[8]↙Z[7]÷2-Z[7]3÷（240×R2）+Z[7]5÷（34560×R4）→Z[9]↙Z[6]2÷（24R）-Z[6]4÷(2688R3)→Z[10]↙Z[7]2÷（24R）-Z[7]4÷(2688R3)→Z[11]↙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]）÷sinA→Z[12]↙Z[9]+(R+Z[11])tan(A÷2)-（Z[10]-Z[11]）÷sinA→Z[13]↙B+Z[12]×cos(E+180)→Z[16]↙C+Z[12]×sin(E+180)→Z[17]↙Z[1]-Z[25]→Z[3]↙Z[3]+Z[6]→Z[4]↙Z[1]-Z[7]→Z[5]↙Goto 3↙Lbl 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↙Lbl 4↙K-Z[3]→Z[23]：0→Z[24]：E→T：Goto 8↙Lbl 5↙K-Z[3]→H↙H-H5÷（40R2×Z[6]2）+H9÷（3456R4×Z[6]4）→Z[23]↙H3÷（6R×Z[6]）-H7÷（336R3×Z[6]3）→Z[24]↙90×H2÷（Rπ×Z[6]）→T↙If O＞0：Then T+E→T：Else E-T→T：T＜0⇒360+T→T：IfEnd↙Goto 8↙Lbl 6↙K-Z[4]→H↙H×180÷（Rπ）+Z14]→T↙R×sinT+Z[8]→Z[23]↙R×（1-cosT）+Z[10]→Z[24]↙If O＞0：Then T+E→T：Else E-T→T：T＜0⇒360+T→T：IfEnd↙Goto 8↙Lbl 7↙Z[1]-K→H↙H-H5÷（40R2×Z[7]2）+H9÷（3456R4×Z[7]4）→U↙H3÷（6R×Z[7]）-H7÷（336R3×Z[7]3）→V↙90×H2÷（Rπ×Z[7]）→T↙Z[13]cosA+Z[12]-U×cosA-V×sinA→Z[23]↙Z[13]sinA-U×sinA+V×cosA→Z[24]↙If O＞0：ThenF-T→T：T＜0⇒360+T→T：Else F+T→T：IfEnd↙Goto 8↙Lbl 8↙If O＜0：Then E-Z[24]→Z[24]：If End↙Z[16]+Z[23]×cosE-Z[24]×sinE→Z[18]↙Z[17]+Z[23]×sinE+Z[24]×cosE→Z[19]↙Return↙子程序反算里程边桩名称：（FSLCBZ）Lb1 1↙"XK="?X:"YK="?Y↙（输入任意测点的XY坐标）:“K=”?K （试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）↙Lb1 2↙Prog"SJK1 ":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"SJK1":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◢（里程）Return↙子程序坐标放样：ZBFYLB1 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=":J→DMS◢(后视方位角) “QXJ=”:T◢（计算里程点切线方位角，可以不显示）“XI=”: Z[18] ◢（中线X）“YI=”: Z[19] ◢（中线Y）Pol(Z[18]-M,Z[19]-N):"DI=":I（中桩放样距）◢J<0⇒J+360→J:"FI=":J→DMS◢（中桩放样方位角）“PJ=”?P◢(输入边桩与线路夹角，左-右+)“PD=”?D◢（输入边桩距）Z[18]+D*COS(T+P) →Z[20] ↙Z[19]+D*SIN(T+P) →Z[21] ↙“XP=”: Z[20] ◢（边桩X）“YP=”: Z[21] ◢（边桩Y）Pol(Z[20]-M,Z[21]-N):"DP=":I◢（边桩放样距）J<0⇒J+360→J:"FP=":J→DMS◢（边桩放样方位角）Return↙。
卡西欧计算器5800程序
If Z≥26157.8 And Z≤27421.915:Then 300000→R:171.915→T:27250→C:32.65→G:-0.00356→I:Goto 0:IfEnd:
If Z≥26915.555 And Z≤27316.952:Then 26915.555→O:11382.792→U:2842.531→V:105。55’51.83”→G:401.396→H:1800→P: 1800→R:1→Q:Goto 0:IfEnd:
……
LbI 0:O:U:V:G:H:P:R:Q:
程序名:1ZD-XY
Lb1 0:Norm 2
F=1：(正反算判别，F=1正算，F=2反算,也可以改F前加？，改F为变量)
Z[1]=90（与路线右边夹角）
Prog＂THB＂：F=1=>Goto 1:F=2=>Goto 2
Lb1 1:Ｆix 3:＂ Ｘ=＂：Locate 6,4,Ｘ◢
＂ Ｙ＝＂：Locate 6,4,Ｙ◢
V+W(Asin(G+QEKW(C+KWＳ))+Bsin(G+QELW(C+LWＳ))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WＳ))) →Y：
G+QEW(C+WＳ)+Ｚ［１］→Ｚ［２］：X+Ｄcos（Ｚ［２］）→Ｘ： Y+Ｄsin（Ｚ［２］）→Ｙ
卡西欧5800公路坐标正反算程序
目录一、坐标正算基本公式………………………………………………………０２二、坐标反算原理 (04)三、高程数据库录入变换 (05)四、计算器程序………………………………………………………０7 ０1、ZBZS(坐标正算)………………………………………………………０702、ZBFS(坐标反算)………………………………………………………０803、ＧＣJF(高程积分） (09)04、PＪFＹ(坡脚放样)………………………………………………………1０05、ＪFCＸ（积分程序) (11)0６、ZＢFY(坐标放样)………………………………………………………１１０７、ＤT(递推） (12)08、H P(横坡） (13)０9、LK(路宽)………………………………………………………１410、SJK1(平面数据库) (14)11、SＪK２(纵面数据库)………………………………………………1４12、SJK3(左路宽度数据库） (15)１3、SＪＫ4(右路宽度数据库) (15)1４、SJＫ5(横坡数据库)......................................................1６１5、SＪK６(下边坡数据库) (16)16、SJK7(左上边坡数据库)………………………………………………1７17、ＳＪK８（右上边坡数据库）………………………………………………1８五、后记 (19)CASIO 58０0计算器公路工程测量程序一、正算所涉及得计算公式 X R i d XαβBd Y d l d βI图表 错误!不能识别的开关参数。
在图1中,A 点为回旋曲线起点,B 点为回旋曲线止点,I 点为所求坐标点。
设: Ａ点得X 坐标为X A ,Y 坐标为Y Ａ ,Ａ点得切线方位角为α,A 点得曲率为ρＡ,A 点得里程为ＬＡ，B 点得曲率为ρＢ,Ｂ点得里程为L Ｂ,I 点得曲率为ρI ,Ｉ点得里程为ＬI 。
卡西欧5800测量坐标计算程序
一．使用说明：1.规定：（1）以道路中线的前进方向区分线路的左右偏转方向，线路左偏Q=-1，线路右偏Q=1，直线Q=0。
（2）当所求点位于线路中线时，Z=0；当位于左侧时，Z取负值；当位于右侧时，Z取正值。
（3）当线元为直线时，其起终点的曲率半径为无穷大，取10的45次方代替。
（4）当线元为圆曲线时，起终点曲率半径为圆曲线半径。
（5）当线元为缓和曲线时，起终点曲率半径与直线相接取10的45次方代替，与圆曲线相接取圆曲线半径。
2.输入与显示说明：（1）输入部分：N？：选择计算方式，输入1表示由里程边距计算坐标，输入2表示由坐标计算里程边距。
X0？线元起点X坐标。
Y0？线元起点Y坐标。
S0？线元起点里程。
F0？线元起点切线方位角。
LS？线元长度。
R0？线元起点曲率半径。
RN?线元终点曲率半径。
Q？线元左右偏标志（左偏=-1，右偏=1，直线=0）S?正算时所求点里程。
Z？正算时所求点距中线边距。
位于线路左侧输入负值，在右侧输入正值。
AR?斜交角度。
（为线路右角，在0-180度之间）X?反算时所求点X坐标。
Y？反算时所求点Y坐标。
（2）显示部分：XS=计算结果，所求点X坐标。
YS=计算结果，所求点Y坐标。
FS=计算结果，所求点切线方位角。
S=计算结果，反算时所求点里程Z=计算结果，反算时所求点边距。
二．程序清单1（用于任意曲线元的坐标计算）1.主程序(TYQXJS)″1.SZ=>XY″:″2.XY=>SZ″:?N:″X0″?U:″Y0″?V:″S0″?O: ″F0″?G:″LS″?H:″R0″?P:″RN″?R:?Q:1÷P→C:(P-R) ÷(2HPR) →D:180÷π→E:If N=1: THEN Goto 1:Else Goto 2:IfEnd:Lbl 1:?S:”D(Z-，Y+)=”?Z:90→T:”O=”?T:Abs(S-O) →W:Prog ″SUB1″:″XS″:X◢″YS″:Y◢″FS″:F◢Goto 1:Lbl 2:?X:?Y:I→X:J→Y:Prog ″SUB2″:″S″:O+W→S:◢″Z″:Z◢Goto 22.正算子程序(SUB1)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) →F:X+Zcos(F+T) →X:Y+Zsin(F+T) →Y3.反算子程序(SUB2)G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z:Lbl 0:Prog ″SUB1″:T+QEW(C+WD) →L:(J-Y)cos(L)-(I-X)sin(L)→Z:IF Abs(Z)<1E-6 :Then Goto 1:Else W+Z→W:Goto 0 : IfEnd:Lbl 1:0→Z:Prog ″SUB1″:(J-Y)÷sin(F)→Z 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。
CASIO 5800道路坐标计算器程序
CASIO 5800计算器程序一、程序功能：本程序由一个主程序（AK）及两个子程序：坐标正算程序（SUB1）、坐标反算程序（SUB2）平曲线数据库（PQX）构成，可以根据曲线段的组成部分：直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、止点坐标、起点切线方位角、线元长度、起点曲率半径、止点曲率半径、里程、边距或坐标），对该曲线段范围内任意里程中边桩坐标进行正反算。
二、源程序：1、主程序（AK）A排在前面，方便运行程序，不用往下翻页了。
Deg：“0→ZX,1→A,2→B”?Z ：“1→JDF,2→XYF”?E：“X(CZ)”?A：“Y”?B 说明：设置角度模式，选择线路0为主线，1为A匝道。
可以自己再增加。
选择计算方法，交点法还是线元法计算坐标，1为交点法，2为线元法。
输入测站坐标，如不需要计算极坐标，即放样的时候需要输入坐标到仪器里面去，则X(CZ)输0，Y不用管可以为任何值，建议把测站坐标输进去，省了放一个点输一次坐标数据到仪器里面去。
只有一条线路可以不需要黄色部分。
“1.Ki→XY,2.XY→Ki”?N：Lbl 0：“Ki”?K：Prog“PQX”IfE=1：Then Prog“JD”：IfEndN=1＝＞Goto 1：N=2＝＞Goto 2Lbl 1?LK-O→W：Prog“SUB1”说明：1为正算，即知道桩号和偏距（注意只能计算边桩，不能计算斜交的点，构造物可以通过后面的ZHIXIAN程序计算），2为反算，即知道坐标，求桩号和偏距。
Ki为桩号，正算直接输，反算需要输入附近桩号，第一次计算最好输好桩号之后，把反算出来的桩号再输进去，重新反算下，如果和第一次反算的桩号一样，那就说明反算的桩号正确。
调用平曲线数据库，输入偏距L（偏距左输负数，右输正数），如果选择的是交点法计算，那需要输入是数据库是交点数据，调用正算子程序SUB1。
If A≠0：Then Pol（X-A,Y-B）：“JULI=”：I◢J<0=>J+360→J：“FWJ=”：J ▶DMS◢Else “Xi=”：X◢“Yi=”：Y◢“FWJ=”：F ▶DMS◢If End ：Goto0说明：正算，如果前面X(CZ)输0，则直接显示X、Y坐标和方位角，否则显示仪器到放样点的距离和方位角。
卡西欧5800计算器的详细编程
CASIO fx—5800计算器工程测量与公路测量程序1、坐标正算”X0=〃?X:”Y0=〃?Y:”I=〃?I:”J=〃?JX+ICos(J)→U:Y+ISin(J)→V”X=〃:U （待求点的X坐标）”Y=〃:V （待求点的Y坐标）说明：X0 Y0：已知点坐标I：两点的距离J：方位角2、坐标反算Lbl 0”X1=〃?X:”Y1=〃?Y:”X2=〃?U:”Y2=〃?VPol(U-X,V-Y):J<0 J+360→J”I=〃:I”J=〃:J DMSGOTO 0说明：X1Y1：第一点的坐标，X2Y2第二点的坐标，I：两点的距离，J：方位角3、圆曲线”X0=〃?X: ”Y0=〃?Y:”FWJ=〃?D:”QDHAO=〃?G: ”ZDHAO=〃?C: ”R=〃?R Lbl 1”DQHAO=〃?LL<G Or L＞C GOTO 2”PY=〃?K180(L-G)÷(πR) →E2RSin(0.5E) →FReC(Abs(F),D+0.5E)：Cls”X=〃:X+I+KCos(D+E+90) →A”Y=〃:Y+J+KSin(D+E+90) →BD+E→o:o<0 360+o→0o”FWJ=〃:o DMSGoto 1Lbl 2”END〃说明：X0Y0：起始点的坐标，FWJ：起始点的方位角，QDHAO：起点里程，ZDHAO终点里程，R：半径，DQHAO：待求点里程，PY：偏移量4、竖曲线”QZHAO =〃?J: ”H+B〃=?B:”I1=〃?C:”I2=〃?D:”R=〃?R:”T=〃?T:((D-C)÷100) ÷Abs((D-C) ÷100)→A:Abs(0.5R(D-C) ÷100)→S”T〃:SJ-T→XJ＋T→YLbl 0”DQ=〃?LIf L<X:Then Goto 1:Else If L＞Y:Then Goto 1:If End:If End”GC=〃:B+C(L-J) ÷100+A(L-J+T) ÷(2R) →HGoto 0Lbl 1”END〃说明：QZHAO：曲中点里程，H+B：曲中点高程，I1：第一坡度，I2：第二坡度，R：半径，T：切线长，DQ：待求点里程。
卡西欧5800计算器的详细编程
CASIO fx—5800计算器工程测量与公路测量程序1、坐标正算”X0=〃X:”Y0=〃Y:”I=〃I:”J=〃JX+ICos(J)→U:Y+ISin(J)→V”X=〃:U （待求点的X坐标）”Y=〃:V （待求点的Y坐标）说明：X0 Y0：已知点坐标I：两点的距离J：方位角2、坐标反算Lbl 0”X1=〃X:”Y1=〃Y:”X2=〃U:”Y2=〃VPol(U-X,V-Y):J<0 J+360→J”I=〃:I”J=〃:J DMSGOTO 0说明：X1Y1：第一点的坐标，X2Y2第二点的坐标，I：两点的距离，J：方位角3、圆曲线”X0=〃X: ”Y0=〃Y:”FWJ=〃D:”QDHAO=〃G: ”ZDHAO=〃C: ”R=〃RLbl 1”DQHAO=〃LL<G Or L＞C GOTO 2”PY=〃K180(L-G)÷(πR) →E2RSin →FReC(Abs(F),D+：Cls”X=〃:X+I+KCos(D+E+90) →A”Y=〃:Y+J+KSin(D+E+90) →BD+E→o:o<0 360+o→0o”FWJ=〃:o DMSGoto 1Lbl 2”END〃说明：X0Y0：起始点的坐标，FWJ：起始点的方位角，QDHAO：起点里程，ZDHAO 终点里程，R：半径，DQHAO：待求点里程，PY：偏移量4、竖曲线”QZHAO =〃J: ”H+B〃=B:”I1=〃C:”I2=〃D:”R=〃R:”T=〃T:((D-C)÷100) ÷Abs((D-C) ÷100)→A:Abs(D-C) ÷100)→S”T〃:SJ-T→XJ＋T→YLbl 0”DQ=〃LIf L<X:Then Goto 1:Else If L＞Y:Then Goto 1:If End:If End”GC=〃:B+C(L-J) ÷100+A(L-J+T) ÷(2R) →HLbl 1”END〃说明：QZHAO：曲中点里程，H+B：曲中点高程，I1：第一坡度，I2：第二坡度，R：半径，T：切线长，DQ：待求点里程。
卡西欧5800程序
卡西欧FX5800计算器公路测量全线中边桩正反算交点法程序修改版1.主程序：一般坐标计算及放样程序-----交点法程序名:1-JDLal 0:10→DimZ:Norm 2:1→A:＂A: XY=1,KMD=2＂? A:A=1=>Goto 1:A=2=>Goto 2: （选择正算还是反算）Lal 1:Lal A:Prog＂QX＂: 90→B ＂PJ＂?B:B→Z[1]: （默认与道路切线方向90度）Lal B:＂KM＂?Z:Z=-1=>Stop:Z=-2=>Break: Prog＂K＂:?D:（输入正算计算里程，输入-1则停止，输入- 2转反算），偏距默认中桩0，左为-，右为+，）Porg＂Z＂:I+Dcos(Z[2]+Z[1]→X:J+Dsin(Z[2]+Z[1]→Y:Fix 3:＂X=＂:Locate 6,4,X: ＂Y=＂:Locate 6,4,Y: ＂FWJ＂=:Locate 6,4,Z[2]◢（显示坐标，切线方位角）Prog “S-Z”: Prog “3-GC”: “H=”：H◢（高程计算并显示，如不需要高程此部分可以省略） Prog＂5-SD＂:＂-CW +QW=＂:M◢（隧道超欠挖，如不需要此部分可以省略）LpWhile Z≠-1(Goto B此处也可以用这个命令)Lal 2 Prog＂QX＂: （选线，默认（线路1） X1=1，（线路2）X2=2.........）Lal C:Do:＂KM＂?Z:Z=-1=>Stop:Z=-2=>Break: （反算时输入大概桩号，输入-1则停止，输入-2正算）Prog＂K＂: ＂X0＂?X: ＂Y0＂?Y: Porg＂B＂Fix 3: ＂KM＂:Lcoate6,4,Z: ＂D=＂:Locate 6,4,D◢Prog “S-Z”: Prog “3-GC”: “H=”：H◢Prog＂5-SD＂:＂-CW +QW=＂:M◢LpWhile Z≠-1(Goto C此处也可以用这个命令)2．线路中线元段判别子程序（K）程序名：KIf Z[3]=0:Then Return: IfEnd:If Z[3]=1:Then Prog “X1”:Return:IfEnd：If Z[3]=2:Then Prog “X2”:Return:IfEnd：If Z[3]=n:Then Prog “Xn”:Return:IfEnd.....（后面可以增加N条线路！！但需要在QX程序里面＂QX:?=0,X1=1,X2=2，Xn=n＂做更改）。
卡西欧5800正反算程序
卡西欧5800P坐标正算、反算程序
坐标正算ZB ZS
已知一个起算的坐标点，根据测得的方位角和平距推算另一点的坐标值。
程序如下:
“X0=”?B: 仪器点X坐标
“Y0=”?C: 仪器点Y坐标
Lbl 0:
“FWJ=”?F: 观测方位角
“PJ=”?D: 观测平距
Rec(D,F):
B+I→X:
C+J→Y:
Cls:
Fix 3:
“Xp=”: 待求点X坐标
“Yp=”: 待求点Y坐标
Locate 7,1,X: 显示在第一行
Locate 7,2,Y◢显示在第二行
Goto 0
注解：
（B，C）是仪器点的坐标，F为方位角，D为平距。
在本程序中，采用了循环（多次）运算多点与仪器点的关系计算法。
如果不必要，只算一个，那么把Lbl 0和Goto 0删除就行了。
坐标反算ZBFS
根据两已知坐标点求算其方位角及距离。
程序如下:
“X0=”?B: 仪器点X坐标
“Y0=”?C: 仪器点Y坐标
“Xp=”?X: 待求点X坐标
“Yp=”?Y: 待求点Y坐标
ol(X-B,Y-C):
J<0 => J+360→J:
I→D:
J→F:
Cls:
Fix 3:
“FWJ=”: 计算的方位角
F►DMS◢按度分秒显示
“PJ=”: 计算的平距
Locate 10,3,D◢平距显示在第三行
注解:
执行程序计算时,首先询问四变量两个点的坐标（B，C）、（X，Y）。
为了明了显示结果具体是何意思，特地使用字符串“PJ=”表示平距，字符串“FWJ=”表示方位角，但要注意是字符串而非变量名。
casio5800积分法线路坐标正反算
CASIO-5800积分法线路坐标正反算—雨中听风施工测量工作中较高档的计算器已很普及，施工放样的数据处理也应该采用程序化计算，以避免出错和提高测量工作的效率。
下面介绍本人炎汝(炎陵至汝城湘粤界)高速公路施工测量工作中应用的CASIO-5800计算器线元积分法坐标正反算程序（程序应用了测量空间里李艳阳老师提供的通用积分公式，在此表示感谢）。
该程序适用于任意线元的坐标、里程及偏距计算，只需更新数据库或者手工输入线元要素便可计之所需。
1.积分坐标正反算主程序(JF.ZBZFS)Lb1 0: Cls↓〝1.SZ => XY〞↓〝2.XY =>SZ〞:?Q↓〝灵活LH.K=0〞↓〝左线ZX.K=1〞↓〝右线YX.K=2〞:?K↓IF K=1:Then Prog〝YS1〞：Else IF K=2:Then Prog〝YS2〞:Else IF K=0:Then 〝QS.FW〞?A: 〝QS.LC〞?O: 〝QS.X〞?U:〝QS.Y〞?V: 〝R0〞?P: 〝RN〞?R: 〝XY.CD〞?L↓Ifend：Ifend：Ifend:？S：S -（字母）O→X：IF X <（数字）0 Or X > L:Then 〝=CUOWU="◢Goto 0 ：Ifend :Q=2=>Goto 2↓Lbl 1：？Z：Prog 〝ZS〞↓Cls：Fix 3：〝N=〞：Locate 3,1,N:〝E=〞：Locate 3,2,E:〝FW=〞：F▲DMS◢Norm 2 :Goto 0↓Lbl 2：〝CX〞？B：〝CY〞？C：B→N：C→E：Prog〝FS〞↓Cls: Fix 3：〝S=〞：Locate 4,1,S:〝Z=〞：Locate 4,2,Z◢Norm 2 :Goto 02. 正算子程序(ZS)0.5(1÷R-1÷P)÷L→D:S-(字母)O→X↓U+∫（cos(A+(X÷P+DX²)×180÷（圆周率）∏)，（数字）0，X）→N↓V+∫（sim(A+(X÷P+DX²)×180÷（圆周率）∏)，（数字）0，X）→E↓A+(X÷P+DX²)×180÷（圆周率）∏→FN+Zcos(F+90)→N:E+Zsim(F+90)→E↓3. 反算子程序(FS)Lb1 1：(数字)0→Z: 1→Q: Prog〝ZS〞↓Pol(N-B+10^(-46) ,E-C+10^(-46)) ↓Isim(F-90-J)→W:S+W→S↓Abs(W)>0.0001=> Goto 1↓Lb1 2: (数字)0→Z : Prog〝ZS〞↓（C-E）÷sin(F+90)→Z数据库子程序(炎汝11标左线)（YS1）Goto 1 ↓Lb1 1 ↓IF S≥45871.346 And S<47089.759:Then194°17°39.6°→A:45871.346→O: 3966.911→U: 733.518→V:10＾45→P: 10＾45→R: 1218.413→L：Return:IfEnd↓IF S≥47089.759 And S<47254.759:Then194°17°39.6°→A:47089.759→O: 2786.220→U: 432.688→V:10＾45→P: -1480→R: 165→L：Return:IfEnd↓IF S≥47254.759 And S<47905.918:Then191°06°01.73°→A:47254.759→O: 2625.621→U: 394.932→V:-1480→P: -1480→R: 651.159→L：Return:IfEnd↓IF S≥47905.918 And S<48070.918:Then165°53°30.93°→A:47905.918→O:1979.924→U:411.883→V:-1480→P: 10＾45→R: 165→L：Return:IfEnd↓IF S≥48070.918 And S<48332.968:Then162°41°53°→A:48070.918→O: 1821.528→U: 458.013→V: 10＾45→P: 10＾45→R: 262.050→L:Return:IfEnd↓IF S≥48332.968 And S<48477.968:Then162°41°53°→A:48332.968→O: 1571.335→U: 535.949→V: 10＾45→P: 800→R: 145→L: Return:IfEnd↓IF S≥48477.968 And S<48915.797:Then167°53°25.75°→A: 48477.968→O: 1431.708→U: 574.858→V: 800→P: 800→R: 437.829→L：Return:IfEnd↓IF S≥48915.797 And S<49060.797:Then199°14°51.64°→A: 48915.797→O:1000.161→U:547.941→V: 800→P: 10＾45→R: 145→L：Return:IfEnd↓IF S≥49060.797 And S<49556.014:Then204°26°24.40°→A: 49060.797→O:866.451→U:491.983→V: 10＾45→P: 10＾45→R: 495.217→L：Return:IfEnd↓IF S≥49556.014 And S≤50236.938:Then204°26°24.40°→A: 49556.014→O:415.608→U:287.091→V: -2500→P: -2500→R:680.924→L：Return:IfEnd↓Lb1 1: A：O：U：V：P：R:L数据库子程序(炎汝11标右线)（YS2）Goto 1 ↓Lb1 1 ↓IF S≥45700.089 And S<47055.633:Then194°17°39.6°→A:45700.089→O:4117.556→U:500744.038→V:10＾45→P: 10＾45→R:1355.544→L:Return:IfEnd↓IF S≥47055.633 And S<47220.633:Then194°17°39.6°→A:47055.633→O: 2803.979→U: 409.351→V:10＾45→P: -1480→R: 165→L：Return:IfEnd↓IF S≥47220.633 And S<47871.792:Then191°06°01.73°→A:47220.633→O: 2643.381→U: 371.595→V:-1480→P: -1480→R: 651.159→L：Return:IfEnd↓IF S≥47871.792 And S<48036.792:Then165°53°30.93°→A: 47871.792→O:1997.684→U:388.546→V:-1480→P: 10＾45→R: 165→L：Return:IfEnd↓IF S≥48036.792 And S<48291.984:Then162°41°53°→A: 48036.792→O: 1839.287→U: 434.676→V: 10＾45→P: 10＾45→R: 255.192→L：Return:IfEnd↓IF S≥48291.984 And S<48436.984:Then162°41°53°→A: 48291.984→O: 1595.643→U: 510.571→V: 10＾45→P: 800→R: 145→L：Return:IfEnd↓IF S≥48436.984 And S<48957.197:Then167°53°25.75°→A: 48436.984→O:1456.015→U:549.480→V: 800→P: 800→R:520.213→L：Re turn:IfEnd↓IF S≥48957.197 And S<49102.197:Then205°08°52.79°→A: 48957.197→O:948.224→U:491.452→V: 800→P: 10＾45→R: 145→L：Return:IfEnd↓IF S≥49102.197 And S<49510.315:Then210°20°25.50°→A: 49102.197→O:820.975→U:422.046→V: 10＾45→P: 10＾45→R: 408.118→L：Return:I fEnd↓IF S≥49510.315 And S<49650.315:Then210°20°25.50°→A: 49510.315→O:468.753→U:215.891→V: 10＾45→P: -1000→R:140→L：Return:IfEnd↓IF S≥49650.315 And S<49909.924:Then206°19°46.96°→A: 49650.315→O:346.337→U:148.025→V: -1000→P: -1000→R:259.609→L：Return:IfEnd↓IF S≥49909.924 And S≤50049.924:Then191°27°18.76°→A: 49909.924→O:101.402→U:64.201→V: -1000→P: 10＾45→R:140→L：Return:IfEnd↓Lb1 1: A：O：U：V：P：R：L数据库变量说明：S≥（线元的起点里程）［可要亦可不要］。
卡西欧5800计算器坐标正反算程序(带详细说明)
卡西欧5800程序完整版＂ZS=1,FS=2＂？N: 1÷P →C:(P-R)÷(2HPH)→D:180÷π→E: If N=1: Then Goto 1: Else Goto 2: IfE nd←┚LbI 1:1 →I:0 →Z：Abs(K-Ｏ) →W: Prog＂ZS＂←┚F-90→F←┚＂JJ＂?K:＂PJ＂?L: ＂ZX=＂：X+Lcos(F+K) →Z[3]◢＂ZY=＂: Y+Lsin(F+K) →Z[4]◢＂QXFWJ=＂:F DMS◢Abs(Z) →Z:Goto 8←┚LbI 2:＂X＂?I:＂Y＂?J: Prog＂FS＂:＂FSZH=＂:Ｏ+W→K◢＂FSJL=＂:Z◢Lbl 8正算子程序（ZS）（不运行）1÷P→C:(P-R)÷(2HPR) →D: 180÷π→E←┚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+ZcosF→X：Y+ZsinF→Y3反算子程序（FS）（不运行）G-90→T：Abs（（Y-V））cos（T）-（X-U）sin（T）→W：0→Z←┚Lbl 0←┚prog＂ZS＂：T+QEW（C+WD）→L：（J-Y）cosL-（I-X）sinL→Z：If Abs(Z)＜0.001: Then Goto1：E lse W+Z→W：IfEnd :Goto 0←┚Lbl 1：0→Z：Prog＂ZS＂：（J-Y）÷sinF→ZZX-WN（运行主程序，该程序适用于匝道，对全段完整非对称曲线计算太过麻烦）＂CZX＂?S(仪站X坐标):＂CZY＂?W（仪站Y坐标）LbI 0: ＂K＂?K: If K≤126.9（第1条单一曲线的终点下面同上）: Then Goto A：Else If K≤176.9 And K＞126.9: Then Goto B:Else If K≤259.562 And K＞176.9: Then Goto C:Else If K≤309.562 And K ＞259.562: Then Goto D：IfEnd:IfEnd:IfEnd:IfEnd←┚LbI A:92094.336（起点X坐标）→U:8287.962（起点Y坐标）→V:0起点桩号→O:145°11＇54＂（方位角）→G:126.9（曲线长度）→H:10^(45) （起点半径）→P:10^(45) （止点半径）→R:+1（路线左右偏，左负右正，直线输0）→Q:Prog ＂WN＂:Goto WLbI B:91990.134→U:8360.388→V:126.9→O:145°11＇54＂→G:50→H:10^(45) →P:150 →R:+1→Q:Pr og ＂WN＂:Goto WLbI C:91947.609→U:8386.570→V:176.9→O:154°44＇51＂→G:82.662→H:150 →P:150 →R:+1→Q:Pr og ＂WN＂:Goto WLbI D:91867.101→U:8399.992→V:259.562→O:186°19＇19＂→G:50→H:150 →P: 10^(45) →R:+1→Q: Prog ＂WN＂:Goto WLbI W:＂PJ1＂?T←┚（该项是只左右边桩再左右偏距，输0为不再左右偏，可以直接出方位角、距离，适用于桥梁桥台、盖梁、涵洞八字墙等放样，以路线前进方向，前-，后+）If T=0：Then Z[3]→Z: Z[4]→M: Goto P: Else Goto V: IfEnd←┚LbI V:＂JJ1＂? O←┚（左右边桩再偏夹角，用180°减去路线夹角）＂X1=＂: Z[3]+Tcos（F+N+ O）→U◢＂Y1=＂:Z[4]+Tsin（F+N+ O）→V◢U→Z: V→M: GotoP←┚LbI P:Pol(Z-S,M-W): If J≤0：Then J+360°→J: IfEnd←┚＂JL=＂:I◢（仪站至测点的距离）＂FWJ=＂:J DMS◢（仪站至测点的方位角）Goto0备注：1．规定以道路中线前进方向（即里程增大的方向）区分左右；当线元往左偏时Q=-1；当元线往右偏时，Q=1：当元线为直线时，Q=0。