经过对多位前辈测量程序的研究及修改，现向广大测友推出鄙人抄袭后的计算程序。本程序经验证准确无误可放心输入，另欢迎留言探讨，如有不足之处还请指教！本程序经过综合考虑５８００的设计缺陷，计算速度较慢故只使用高斯四节点法为计算内核．

计算器主程序：ZHU-CHENG-XU

Lbi0:“1，ZS=FS，2GC，3ZS，4FS，5。。。。”？U:

U=1＝＞Prog”ZS-FS”:  进入公路三维程序

U=2＝＞Porg”GC”:    进入坐标正算程序

U=3＝＞Porg”ZS”:    进入坐标反算程序

U=4＝＞Porg”FS”:   进入大地坐标转施工坐标程序

…………Goto0:

说明：计算器总的主程序，进入选择各种分支计算程序。1为公路三维计算，2为普通正算，3为普通反算，4为大地坐标转施工坐标。。。。选择错误重新选择。此程序可以不用输入，只为给大家一个思路,可以把计算器所有程序集中到一个主程序内管理.

正反算选择程序：ZS-FS

Deg:              //设置角度模式    

“1LZ＝＞XY，2XY＝＞LZ”？U： //正反算选择，正算选1，反算选2,坡口坡脚选3

If U=1: ThenProg”ZS-XH”:IfEnd:             //进入正算循环主程序

IfU=2:Then Prog”FS-XH”:IfEnd:              //进入反算循环主程序 

正算循环主体程序ZS-XH

“XL-XZ“？U：                           //选择线路   1～N

Lbi0:”L=”?L:”Z=”?Z:                  //输入桩号和边距

Porg”ZS-XH-1“：Goto0:   // 运行正算循环子程序

正算循环子程序ZS-XH-1

Prog”PM-XL-XZ”:                      //进入平面线路选择程序对线路线元进行选择

Prog”CHAOXIAN-CL”:                    //里程超限处理

Prog”GL-ZS”:                          //进入坐标正算程序

Prog”GL-BZ-ZS”:                        //进入公路边桩正算程序

Prog”ZS-XS”:                          //进入正算显示程序，提取正算三维结果

反算循环主体程序FS-XH

0→Z:“XL-XZ”？U：                         //选择线路    1～N

Lbi0：

“L”？L：“X=”？T：“Y=”？S：            //输入起算假定桩号，实测X Y

Prog”GL-FS”:                                /进入公路反算程序

Prog”FS-XS”:                                //进入反算显示程序

Goto0:

坐标正算程序GL-ZS

.1739274226→I:.3260725774→J:

.0694318442→K:.3300094782→E :                 //四节点法计算系数录入

1-E→F：1-K→M：

Q+KrO(C+KD) →K: //数据库要用X、Y、Q作为起点坐标方位角

Q+ErO(C+ED) →E: //r是弧度,输入方法FUNCTION 选5

Q+FrO(C+FD) →F:

Q+MrO(C+MD) →M:

X+O(ICos(K)+JCos(E)+JCos(F)+ICos(M)) →X:   //计算出X坐标

Y+O(ISin(K)+JSin(E)+JSin(F)+ISin(M)) →Y:      //计算出Y坐标

Q+Or(C+D) →Q：                            //计算出算点方位角,O是字母,r是弧度:

边桩坐标计算程序GL-BZ-ZS:

X+Zcos(Q+90) →G:                              //计算出边桩X坐标

Y+ZSin(Q+90) →M:                              //计算出边桩Y坐标

公路反算程序GL-FS

Do:Porg”PM-XL-XZ”:                         //进入平面线路选择程序

Porg”CHAOXIAN-CL”:                       //里程超限处理程序

Porg”GL-ZS”:                            //进入公路正算4节点程序

T-X→I:S-Y→J:                          //试算点与测点X，Y距离差

Pol(I,J):J-Q→J:                               //方位角转换为接近试算点到测点的值：

Rec(I,J):                                     //计算出X差和Y差（X里程，J边）

I/(1-CJ-2DJ) →V:                          //计算出里程大致差距（此公式为加速反算）

L+V→L：                                 //修改后的里程

Abs(V) →V:                               //修改值的绝对值

LpWhile V＞0.001:                            //小于一毫米通过，大于一毫米循环计算

正算显示ZS-XS

“FWJ”：Q◢                                    //显示计算里程方位角

“ZHONG-X”:X◢                                //显示中线X坐标

“ZHONG-Y”: Y◢                               //显示中线Y坐标

“B-X”:G◢                                     //显示边桩X坐标

“B-Y”: M◢                                  //显示边桩Y坐标

反算显示FS-XS

“L”：L◢                                      //显示计算里程        

“Z”：J◢                                       //显示计算边距

里程超限CHAOXIAN-CL

If L＜A:Then “L-＜＜＜”；Stop:IfEnd:        //小于本路线最小里程提示，退出程序

If L＞B: Then”L+＞＞＞”:Stop:IfEnd:            //大于本路线最大里程提示，退出程序

数据库部分

1平曲线数据库

平面线路选择程序PM-XL-XZ

U=1＝＞Prog”X1”:                         //进入线路1

U=2＝＞Prog”X2”:                             //进入路线2

U=3＝＞prog”X3”:                              //进入路线3

（P-R）/（2HPR）→D：

L-O→O：                                    //正算程序需要补充部分

OD→D:                                      

1/P→C：                                     

（数据库）线路X1              

线路起点里程→A：线路终点里程→B：

If L≤第一线元止点里程：Then 起点X坐标→X：起点Y坐标→Y：

起点半径（左-，右+）→P：止点半径（左-，右+）→R：本线元长度→H：

线元起点桩号→O：起点方位角→Q：Return:IfEnd:

If L≤第二线元止点里程：。。 。。Return: IfEnd:

  。。  。。

当曲率半径为0时输入10的45次方

二号线路X2

三号线路X3

一，坐标正算方法

1，进入主程序ZHU-CHENG-XU

2，选择使用程序，根据大家录入的代码，本程序，选择1

3，选择坐标1，正算，2，反算，选择1

4，线路选择 1～N，根据计算路线编号进行选者

5，输入计算里程

6，输入计算边距，左-，右正

7，显示计算点方位角，显示计算点中桩XY，显示边桩XY。

二，坐标反算方法

1，进入主程序ZHU-CHENG-XU

2，选择使用程序，根据大家录入的代码，本程序，选择1

3，选择坐标1，正算，2，反算，选择2

4，线路选择 1～N，根据计算路线编号进行选者

5，输入试算点里程L  可输入任意里程，数据库有数据就好，离真实坐标越近越快

6，输入测点X，Y

7，显示测点里程和边距