坐标反算程序
线元法线路坐标正反算程序
经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高, 功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。
功能简介及特点：1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。
2、线元数据可自动从数据库调用，也可手工输入。
3、可管理多条线路，如里程不在线路或线元范围，将警告里程偏大、偏小。
4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。
5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。
第三次及以后试算才要求输入近似里程。
6、程序代码规范简洁，便于阅读、理解。
完整程序清单：ZFS %正反算主程序B=.1739274226:C=.5-B:Lbl 1:U"0 ZS 1 FS"=0=>Prog "ZS":工>U=1=>Prog"FS":工>Goto 1ZS %正算子程序{K}:Prog"ZZ":l=0:{l}:l"L" 丰 0=>"Prog"WY":丰 >Prog"ZB"FS %反算子程序{KVW}:V"XC"W" YC":Lbl 2:Prog "ZZ":I=V-S:J=W-T:Pol(I,J: J=J~F:K二K+Rec(l,J:Absl<1m=>Prog"WZ": 半 >Gota\2M=O:{M}:M"O NEXT"=O二>U=U+1:Goto 2: 半 >U=1ZZ %高斯法中桩子程序(4节点)Prog"XL":M=K-L: O=(P-R) 2PQR:D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D:l=5:Lbl 1:C[l]=A+MrC[l](1 P+OMC[l]:Dsz l:Goto 1:S=X+M(BcosD+CcosE+CcosG+BcosH:T=Y+M(Bsi nD+Csi nE+Csi nG+Bsi nHWY %外移点计算子程序Lbl 1:J=90:{J}:J=F+J"<":F=J:S=S+Rec(l,J:T=T+J: Prog"ZB":l=0:{l}:l"L" 工0=>Goto 1WZ %位置显示子程序"KJ":K:Pause 1:J丄ZB %坐标显示子程序"XY ":S:Pause 1:T AYC %异常处理子程序U=1=>K=L:U=2 △U=3=>K=M:U=4AU=5=>{K}:U=4 △K<L=>"V<!":Z=1△K>M=>">>!":Z=1DL %断链处理子程序"DL":K=L:I>0=>K=L+Q------------------ 以上为程序运算部分,以下为数据库部分XL %线路数据库选择子程序Lbl 1:Z=0:N"0 SD"=0=>Prog"0" △N=1=>Prog"1" △N=2=>Prog"2" △有几条线路仿上行格式输几行Z=1=>{NLX YOPQRK}:Goto 10 %手工输入子程序L"K0"X YAQ"LS"P"R0"R"RN":M二L+Q:Prog" YC"1 %线路一数据库子程序①Lbl B:L=线路起点里程:M=线路终点里程:Prog"YC":Z=1=>Goto E △②Q=线元长:P=起点半径:R=终点半径:K< L+Q二>X=起点X坐标:Y二起点Y坐标:A=起点方位角:Goto E △ L=L+Q:③…… @Q=短链长:K<L+Q=>Prog "DL":Goto B △ L=L+Q:⑤Q=线元长:P=起点半径:R=终点半径:K＜L+Q二＞X=起点X坐标:Y二起点Y坐标:A=起点方位角:Goto E △ L=L+Q:⑥……炉0=线元长:P=起点半径:R=终点半径:X=起点X坐标:丫二起点Y坐标:A=起点方位角丄bl E2 %线路二数据库子程序输入要求和线路一相同。
坐标反算程序
坐标正反算程序(5800)反算主程序:FS1.Prog"SJ"2.90U÷π÷R→H 第一回旋线所对应圆心角β013.90V÷π÷R→T 第二回旋线所对应圆心角β024.Rec((Z[3]+Z[4]÷cosK)÷tanK,F+90G):Z[5]+I→Z[11]:Z[6]+J→Z[12]ZH 、HZ点垂线交点坐标5.F-90G→Z[13]ZH点与ZH、HZ点垂线交点连线方位角6.Z[13]+GH→Z[14]HY点与圆心连线方位角7.Z[13]+G(K-T)→Z[15]YH点与圆心连线方位角8.Z[15]+GT→Z[16]HZ点与ZH、HZ点垂线交点连线方位角9.Lbl 0:?X:?Y10.Pol(X-Z[11],Y-Z[12])11.If GIsin(J-Z[13])<0:Then Pol(N-X,E-Y):Rec(I,F-J)12.Q-I→Z: J→S:Cls:“Z=”:Locate3,1,Z:”S=”:Locate 3,2,S◢13.Goto 0:IfEnd 第一直线14.If GIsin(J-Z[16])>0:Then Pol(X-N,Y-E):Rec(I,J-F-GK)15.I-Z[4]+D→Z: J→S:Cls:“Z=”:Locate3,1,Z:”S=”:Locate 3,2,S◢16.Goto 0:IfEnd 第二直线17.If GIsin(J-Z[13])>0:Then Pol(X-Z[9],Y-Z[10])18.If GIsin(J-Z[14])<0:ThenZ[13]→P:U→H:A→M:1→T:Z[5]→Z[17]:Z[6]→Z[18]:F→Z[19]:Goto 1:IfEnd 第一缓和曲线19.If GIsin(J-Z[15])<0:Then Abs(J-Z[14])→Z[20]:Z[20]>360?Z[20]-360→Z[20]20.B+πRZ[20]÷180→Z: G(R-I)→S:Cls:“Z=”:Locate3,1,Z:”S=”:Locate 3,2,S◢21.Goto 0:IfEnd 圆曲线22.If GIsin(J-Z[15])>0:Then Z[16]→P:V→H:D→M:-1→T:Z[7]→Z[17]:Z[8]→Z[18]:F+GK+180→Z[19]:Goto 1:IfEnd 第二缓和曲线23.Lbl 1:Pol(X-Z[11],Y-Z[12]): Abs(J-P)→J:J>180?J-360→J:JπR÷180→L L为以R为半径,以所求点-(Z[11],Z[12])连线与P的角度差计算得的弧长24.Lbl 2:L-L5÷(40R2H2)+L9÷(3456R4H4)-L13÷(599040R6H6)→O25.L3÷(6RH)-L7÷(336R3H3)+L11÷(42240R5H5)-L15÷(9676800R7H7)→P26.Rec(sqrt(O2+P2),Z[19]+TGtan-1(P÷O)):Z[17]+I→O:Z[18]+J→P27.Pol(X-O+0.00001,Y-P+0.00001):Rec(I,J-(Z[19]+90TGL2÷(πRH)))28.If AbsI>0.001:Then L+I→L:Goto 2:Else M+TL→Z:TJ→S:Cls:“Z=”:Locate3,1,Z:”S=”:Locate 3,2,S◢29.Goto 0:IfEnd子程序(曲线要素数据):SJ1.20→DimZ:?M2.If M=1:Then *→F:*→K:*→G:*→U:*→V:*→R:*→Q:*→N:*→E:Goto 1:IfEnd3.If M=2:Then *→F:*→K:*→G:*→U:*→V:*→R:*→Q:*→N:*→E:Goto 1:IfEnd4.………………5.Lbl 1:U2÷24÷R-U4÷2688÷R3+U6÷506880÷R5→Z[1]Z[1]第一回旋线内移值P16.V2÷24÷R-V4÷2688÷R3+V6÷506880÷R5→Z[2]Z[2]第二回旋线内移值P27.(R+Z[2]-(R+Z[1])cosK)÷sinK→M8.(R+Z[1]-(R+Z[2])cosK)÷sinK→P9.M+U÷2-U3÷240÷R2→Z[3]Z[3]第一切线长T110.P+V÷2-V3÷240÷R2→Z[4]Z[4]第二切线长T211.πRK÷180+(U+V)÷2→L 曲线长(Ls1+圆+Ls2)12.Q-Z[3]→A:A+U→B:A+L→D:D-V→C ZH,HY,HZ,YH里程桩号13.0→I:0→J:Rec(Z[3],F+180):N+I→Z[5]:E+J→Z[6]ZH点坐标14.Rec(Z[4],F+GK):N+I→Z[7]:E+J→Z[8]HZ点坐标15.Rec(sqrt(M2+(R+Z[1])2),F+G(K+tan-1((R+Z[2])÷P))):N+I→Z[9]:E+J→Z[10]Z[9],Z[10] 圆心坐标注:F-起始边方位角K-转角G-线路左转输-1,右转输1U-第一回旋线长V-第二回旋线长R-圆曲线半径Q-交点桩号N,E-交点坐标Z-待求点桩号S-偏距(左偏为-,右偏为+)sqrt为根号。
坐标正算与反算、施工与大地坐标换算程序(4800P)
Goto 1
算例：1点，X=10.000，Y=20.000
2点，X=500.000，Y=600.000
3点，X=700.000，Y=400.000
专业的3S站
DB子程序：V=V+1:B=Z[V]+T:Prog“A”:D=Z[V+N]:Prog“XY”
DEG子程序：P=Frac B*100:B=INTB+INTP/60+Frac P/36
三、 支导线的计算
支导线就是从一个已知点出发，既不符合到另一个已知点也不回到原来的点上。支导线没有检核条件，不易发现错误，过不宜采用。
X？ 输入5点X坐标
Y？ 输入5点Y坐标
S=218500 （1.2.3.4.5点围成的面积）
X？ 输入6点X坐标
Y？ 输入6点Y坐标
S=230500 （1.2.3.4.5.6点围成的面积）
.........
.........
Lbl 3：｛XY｝：X“X=”：Y“Y=”：A=（X－C）cos R+(Y－D) sin R:“A=” ：Pause 0:B=（Y－D）cos R－(X－C) sin R:B“B=” ▲ Goto 0 △
Lbl 4：｛AB｝：A“A=”：B“B=”：X=C+A cos R－B sin R：“X=”：Pause 0:Y=D+A sin R+B cos R：Y“Y=” ▲ Goto 0 △
坐标正算与反算、施工与大地坐标换算程序（4800P）
Deg：Fix4：J=0：C“X=”：D“Y=”：R“R=“：
公路坐标正反算程序
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名： JDB
X-Scos(M)→Z[19]:
(Abs(K-O)) →J：Prog"SUB1"：
(F-M) →F
Return
2. 反算主程序 GSFS
Lbi 0：?X：?Y：X→Z[2]：Y→Z[3]：
“QDXO”？I："QDY0"?S:"QDLC"?O:"QDFWJ "?G："ZDLC"?H："QDR"?P："ZDR"?R：”Q(Z=-1 ZX=0 Y=1)” ?Q:
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]↙
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↙
卡西欧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: Goto4:Else if Q<S: Then K+Q→K: Goto 2 :Else if Q>S :Then K-Q→K:Goto2: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 : Pro g”JSPB”:Return: Ifend↙if…………Prog”JSPB”:Return: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）。
CASIOfx-4500P,4800P程序
说明：I(V)、J(W)中，不带括号的适用于CASIO fx-4500P型计算器，带括号的适用于CASIO fx-4800P型计算器。
作者：徐万鹏
坐 标 反 算 子 程 序
行
MODE EXP 程 序
程 序 说 明
F5
Z
1
程序区及程序名。
:
Y
=
Z
(
B
+
C
Sin
D
+
E
Sin
P
-
F
Cos
P
说明：
1、本子程序5行共160步。
2、计算结果存入X、Y寄存器。
3、G2和H2各为两步。
说明：1、V(I)、W(J)中，不带括号的适用于CASIO fx-4500P型计算器，带括号的适用于CASIO fx-4800P型计算器。
2、为提高计算精度，曲线计算数学模型中的级数展开式取前三项，教课书中一般仅取前一至二项。
-
Y
Cos
P
:
″
计算第一切线长L1、第二切线长L2。
T
1
=
″
◢
Y
:
″
T
2
=
″
说明：1、本程序4行共143步。
2、输入：显示ZX:Y=1,Z=-1输入曲线转向；显示BJ 输入曲线半径；显示PJ输入曲线偏角；显示L1输入第一缓和曲线长；显示L2输入第二缓和曲线长。
4、输出：QXC= 曲线长；T1=第一切线长；T2=第二切线长。
/
2
R
:
D
=
N
求圆曲线坐标增量计算元素C和D。
第六章 坐标正算和反算
第六章→第三节→导线测量内业计算导线计算的目的是要计算出导线点的坐标，计算导线测量的精度是否满足要求。
首先要查实起算点的坐标、起始边的方位角，校核外业观测资料，确保外业资料的计算正确、合格无误。
一、坐标正算与坐标反算1、坐标正算已知点的坐标、边的方位角、两点间的水平距离，计算待定点的坐标，称为坐标正算。
如图6-6 所示，点的坐标可由下式计算：式中、为两导线点坐标之差，称为坐标增量，即：【例题6-1】已知点A坐标，=1000、=1000、方位角=35°17＇36.5"，两点水平距离=200.416，计算点的坐标？35o17＇36.5"=1163.58035o17＇36.5"=1115.7932、坐标反算已知两点的坐标，计算两点的水平距离与坐标方位角，称为坐标反算。
如图6-6可知，由下式计算水平距离与坐标方位角。
（6-3）（6-4）式中反正切函数的值域是-90°～+90°，而坐标方位角为0°～360°，因此坐标方位角的值，可根据、的正负号所在象限，将反正切角值换算为坐标方位角。
【例题6-2】=3712232.528、=523620.436、=3712227.860、=523611.598，计算坐标方位角计算坐标方位角、水平距离。
=62°09＇29.4"+180°=242°09＇29.4"注意：一直线有两个方向，存在两个方位角，式中：、的计算是过A点坐标纵轴至直线的坐标方位角，若所求坐标方位角为，则应是A点坐标减点坐标。
坐标正算与反算，可以利用普通科学电子计算器的极坐标和直角坐标相互转换功能计算，普通科学电子计算器的类型比较多，操作方法不相同，下面介绍一种方法。
【例题6-3】坐标反算，已知=2365.16、=1181.77、=1771.03、=1719.24，试计算坐标方位角、水平距离。
任意曲线坐标放样正反算程序(4800,4850,5800均可以用)
任意曲线坐标放样正反算程序（fx－4800p）正算ZSLbl 0↙｛K,L,Q,S｝:K:L:Q:Prog“FXL”:Prog“SUB”↙“X=”:X=X:Pause0:Y=Y◢K=S+K↙Goto0↙反算FSLbl 0:｛PT｝:P“X0”:T“Y0”↙Lbl 1:｛K｝:L=0:Q=0:Prog“FXL”: Prog“SUB”↙U=I-90↙V=（T-Y）CosU-（P-X）SinU↙AbsV≥0.001=>K=K+V:Goto1:≠>“K0=”:K=K:Pause0:W“L0”=（T-Y）CosI-（P-X）SinI◢Goto0↙子程序SUBNorm↙Deg↙U=（E-D）÷Abs（G-F）:V=Abs（K-F）:W=UV↙J=C+90V（W+2D）÷π:J<0=>J=J+360⊿J≥360=>J=J-360⊿↙M=C+45V（W÷4+2D）÷（2π）:N=C+135V（3W÷4+2D）÷（2π）:O=C+45V（W÷2+2D）÷π↙I=J+Q↙Fix3↙X=A+V÷12×（CosC+4（CosM+CosN）+2CosO+CosJ）+L CosI↙Y=B+V÷12×（SinC+4（SinM+SinN）+2SinO+SinJ）+L SinI↙数据库FXLK≥0=>K≤1211.753=>A=9764.547:B=3344.148:C=142，6，26，:D=0:E=0:F=0:G=1211.753⊿⊿↙K≥1211.753=>K≤1406.455=>A=8808.279:B=4088.39:C=142，6，26，:D=1÷700:E=1÷700:F=1211.753:G=1406.455⊿⊿↙使用说明：K－里程桩号; A,B－每段线元的起点X,Y坐标; C－每段线元的起点的切线方位角; F－每段线元的起点桩号;G－每段线元的终点桩号;D－每段线元的起点半径曲率;E－每段线元的终点半径曲率; （左转取-1÷半径曲率,右转取1÷半径曲率,直线为下一个桩号（例: 赋值20确认,即上个桩号加20）D,E取值说明：当线元为直线段,不论其起,止与何线元相接, D,E都取0.当线元为圆曲线,不论其起,止与何线元相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径.当线元为缓和段, 起,止与直线相接, D,E都取0. 与圆曲线相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径. 与缓和段相接, D,E都取0.任意点的切线方位角C说明：编数据库时,若两曲线相接,例：先编ZH至HY段,运行算出HY点的切线方位角,(即求HY点的坐标,然后终止程序调出C值),然后再以计算出的C编HY至YH段.依此类推.反算桩号时需多次确认所输入的X,Y坐标,程序运算至桩号误差小于允许值时,会同时显示结果:K0-桩号,L0-距离.坐标反算程序:ZBFSC“X1”：D“Y1”：E“X2”：F“Y2”：Fixm：Pol(E－C,F－D：I“S=”◢J≤0=>J=J+360⊿J“J=”↙坐标正算程序ZBZSC“X1”：D“Y1”：L：O“J”：Fixm：X=C+Rec(L,O◢Y=D+J↙竖曲线SQXJ“BPD”：B“BPH”：C“I1”：D“I2”：RW=D－C◢A=W÷AbsW：T=Abs（0.5RW）◢L=2T↙E=T2÷2R◢X“QD”=J－T◢Y“ZD”=J+T◢Lbl1：P“K”：H=B+C（P－J）+A（P－J+T）2÷2R◢W“CONT”：W=1=>Goto1：≠>Goto2：⊿↙Lbl2↙说明:BPD－变坡点桩号BPH－变坡点高程I1－前一个坡度（上“+”,下“-”）I2－后一个坡度（上“+”,下“-”）R－半径E－外距QD－该竖曲线起点桩号ZD－该竖曲线终点桩号K－要求的桩号H－要求桩号的高程任意曲线坐标放样正反算程序（fx－4850p）正算ZSLbl 0↙｛K,L,Q,S｝:K:L:Q:Prog“FXL”:Prog“SUB”↙“X=”:X=X:Pause0:“Y=”:Y=Y◢K=S+K↙Goto0↙反算FSLbl 0:｛PT｝:P“X0”:T“Y0”↙Lbl 1:｛K｝:L=0:Q=0:Prog“FXL”: Prog“SUB”↙U=I-90↙V=（T-Y）CosU-（P-X）SinU↙AbsV≥0.001=>“K0=”:K=K+V:Goto1:≠>“K0=”:K=K:Pause0:W=（T-Y）CosI-（P-X）SinI ↙“L0=”:W=W◢Goto0↙子程序SUBNorm↙Deg↙U=（E-D）÷Abs（G-F）:V=Abs（K-F）:W=UV↙J=C+90V（W+2D）÷π:J<0=>J=J+360⊿J≥360=>J=J-360⊿↙M=C+45V（W÷4+2D）÷（2π）:N=C+135V（3W÷4+2D）÷（2π）:O=C+45V（W÷2+2D）÷π↙I=J+Q↙Fix3↙X=A+V÷12×（CosC+4（CosM+CosN）+2CosO+CosJ）+L CosI↙Y=B+V÷12×（SinC+4（SinM+SinN）+2SinO+SinJ）+L SinI↙数据库FXLK≥0=>K≤1211.753=>A=9764.547:B=3344.148:C=142，6，26，:D=0:E=0:F=0:G=1211.753⊿⊿↙K≥1211.753=>K≤1406.455=>A=8808.279:B=4088.39:C=142，6，26，:D=1÷700:E=1÷700:F=1211.753:G=1406.455⊿⊿↙使用说明：K－里程桩号; A,B－每段线元的起点X,Y坐标; C－每段线元的起点的切线方位角;D－每段线元的起点半径曲率;E－每段线元的终点半径曲率; （左转取-1÷半径曲率,右转取1÷半径曲率,直线为0. ）;L－距中桩距离（左负,右正,中桩为0）; Q－与线路前进方向的右交角（正交为90）; S－赋值加上一个桩号放下一个桩号（例: 赋值20确认,即上个桩号加20）D,E取值说明：当线元为直线段,不论其起,止与何线元相接, D,E都取0.当线元为圆曲线,不论其起,止与何线元相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径.当线元为缓和段, 起,止与直线相接, D,E都取0. 与圆曲线相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径. 与缓和段相接, D,E都取0.任意点的切线方位角C说明：编数据库时,若两曲线相接,例：先编ZH至HY段,运行算出HY点的切线方位角,(即求HY点的坐标,然后终止程序调出C值),然后再以计算出的C编HY至YH段.依此类推.反算桩号时需多次确认所输入的X,Y坐标,程序运算至桩号误差小于允许值时,会同时显示结果:K0-桩号,L0-距离.任意曲线坐标放样正反算程序（fx－4800p）任意曲线坐标放样正反算程序（fx－4850p）。