2024年3月29日发(作者：建瑶岑)

5800公路测量程序使用说明

一、程序使用流程

本程序数据和主程序是分开的，编程时将不同的工程数据存放到不同的数据文件里，如A匝道，文件

名为A，将匝道A所有的曲线线元参数输入A文件里。运行时只要运行文件名A的程序就可以了，具体运

行流程见下图：

二、数据文件的编写

（一）交点法数据文件编辑

交点法编写数据文件必须是对称型的，即直线段→缓和曲线段→圆曲线段→缓和曲线段→直线段，（如

果任意一端没有直线段，则把直线段长度看做是0），另外圆曲线两侧缓和曲线的旋转常数必须相等，并且

和直线段连接处的半径必须是无穷大。

交点法数据文件编写一般是根据设计图纸提供的平面曲线参数一览表提供的参数来编写，每个弯道包

括：弯道起点方位角（C），交点X坐标（D），交点Y坐标（E），缓和曲线长度（F，当没有设缓和曲线

时，F=0），交点转交（G，向左转弯，G为负值，向右转弯，G取正值），交点桩号（H），弯道圆曲线半

径（R）。

下图是一段市政道路设计参数数据。

根据上图提供的数据，可以编辑成如下的数据文件：

文件名:CHLNR

3→DimZ “X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Els

e “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L>0 :Then 98°39°35.12°→C：4474.384→D：2415

.861→E：140→F：31°17°23°→G：410.007→H:600→R:IfEnd: If L>1060 Then 129°56°58.19°→C:

4206.421→D：3093.946→E：70→F：-33°50°48°→G：1285.437→H:600→R:IfEnd:Prog”XLJS”:If Z[3]≠0.

0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（

O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L

◢

“K=”:K

◢

IfEnd:Goto 0: IfEnd

在面程式中，有两个条件转移语句即If L>0：Then 98°39°35.12°→C:4774.384→D: 2415.861→E:140→

F:31°17°23°→G:410.007→H:600→R:IfEnd

If L>1060：Then 129°56°58.19°→C:4206.421→D: 3093.946→E:70→F:-33°50°48°→G:1285.437→H:600

→R:IfEnd

……

如果还有其他弯道，可以继续完后加。在这些存放设计参数的语句前后的程序表达式是固定的。

说明：

编辑曲线参数时，每个曲线参数放在一个If L>***.***（两个弯道中间直线段上的任意桩号） :Then ***°**°

**.**°→C（弯道起点方位角）：****.***→D（交点X坐标）：****.****→E（交点Y坐标）：***.***→F（缓

和曲线长度：**°**°**°→G（转角，向左转为负值，向右转为正值）：***.***→H（交点里程桩号）:***

→R（圆曲线半径）:IfEnd 条件式语句里，如果有多个弯道，一直按上述形式编下去，

变量说明：

1、 L>***.***

弯道参数的起点范围，可以是弯道起点，也可以是弯道前直线段里的任意点桩号

2、 ***°**°**.**°→C 弯道起点方位角

3、****.***→D、****.****→E

交点X坐标、Y坐标

4、***.***→F

缓和曲线长度，如果没设置缓和曲线，则缓和曲线长度看做0

5、**°**°**°→G 转角，向左转为负值，向右转为正值

6、***.***→H、交点里程桩号

7、***→R 圆曲线半径

注意，运行程序时，第一部是输入测站坐标X0，Y0，然后输入近似桩号，主要是为了加快收敛，输入

时注意当根据桩号和距中距离计算所求点坐标时，请不要把小数点第四位输入0.0001，如K21+369.0921，请

把桩号输入21369.092即可，相反，当根据坐标反算桩号和距中距离时，请在输入时把第四位小数数位0.0001，

如K21+200，输入时请输入21200.0001。

（二）线元法数据文件编辑

线元法数据文件编辑适合于各种各样的线型组合，特别是在大型立交区或者是不规则的线型组合地方，

如小区道路、重丘区山路等有特别明显的有点，本程序为了减少使用者设计参数的输入，特别把直线与缓和

曲线以及圆曲线连接的地方只要输入缓和曲线参数就可以了，这样可以减少输入很多设计参数，同时节约大

量的内存。看下图，

该图是一个半喇叭型的互通立交，一共是5条匝道，我们以D匝道为例，来看看线元法数据文件的建立。

下表是D匝道的线元表

：

线元名称 起点桩号 偏向

半径或

缓和曲

线旋转

常数A

55

40

65

线元起点坐标、方位角

线元长度

(m)

39.289

39.216

105.625

X(m)

494341.902

494347.926

494380.677

Y(m)

478027.825

477989. 873

477971.271

方位角

260°38′55.7″

302°18′26.0″

358°28′45.9″

缓和曲线

圆曲线

缓和曲线

DK0+000.00

DK0+039.289

DK0+078.504

右偏

右偏

右偏

D匝道线元数据分析：从上表中可以看出，D匝道由三个线元组成，第一段是缓和曲线段，根据上表判

断，该缓和曲线属于第一段缓和曲线，其最小半径处桩号为J=K0+039.289，该处坐标X=494347.926,Y=477989.

873，半径为R=40米，方位角W=302°18′26.0″，缓和曲线旋转常数P=55，在缓和曲线终点K0+039.289处，

有一圆曲线与其相连接，长度I=39.216。第二段是圆曲线，第三段又是一段缓和曲线，根据上表判断，该曲线

属于第二缓和曲线，曲线最小半径位于与园曲线连接处，桩号J=K0+78.504，坐标为X=494380.677，Y=477971.271,

方位角W=358°28′45.9″，最小半径R=40，缓和曲线旋转常数P=-65（为了判断是第一段还是第二缓和曲线，

我们通过P的正负来判断，第一缓和段P统一用正数表示，第二缓和段统一用负号表示。起点处相连的圆曲线

长度I可以看做0，也可以看做是39.216.由此我们整理下。从整个线型全部是向右转向，所以R全部取正数。

具体如下：

第一段缓和曲线线元参数：

最小半径处坐标X=494347.926 Y=477989. 873

最小半径R=40

最小半径处方位角：W=302°18′26.0″

缓和曲线旋转常数P=55

与之相连的圆曲线长度I=39.216

最小半径处桩号J=39.289

此段数据计算范围可以从K0+000~ K0+78.504

第二段缓和曲线线元参数

最小半径处坐标X=494380.677 Y=477971.271

最小半径R=40

最小半径处方位角：W=358°28′45.9″

缓和曲线旋转常数P=65

与之相连的圆曲线长度I=39.216

最小半径处桩号J= W=358°28′45.9″

此段数据计算范围可以从K0+78.504到匝道结束。

当桩号L位于K0+000~ K0+78.504时，曲线参数按第一段编写，当L>78.504时，按第二段缓和曲线编

写，当然第二段缓和曲线计算范围也可以从圆曲线起点K0+039.289到K0+78.504中间任意桩号开始。下面

是D匝道的数据程序文件（文件名“DZD”）：

3→DimZ：X0”?A:”Y0”?B:“Ln”?L: Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L:

Else “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494347.926→X：477989. 873→Y：40

→R：302°18°26°→W：55→P:39.216→I: 39.289→J：IfEnd: If L≥78.504：Then 494380.677→X：4779

71.271→Y：40→R：358°28°45.9°→W：-65→P:39.216→I: 78.504→J：IfEnd:Prog”FDCH”: Prog”XLJS”:

If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]

－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L

◢

“K=”:K

◢

IfEnd:Goto 0: IfEnd

在上面程式中，每个缓和曲线的线元参数放在一个If … Then …… IfEnd循环语句中，If 后面紧跟计算的

起点范围，终点范围从下段的起点范围结束。

变量定义：

1、X、Y 输入缓和曲线半径最小处的X、Y坐标

2、R 缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。

3、W 输入缓和曲线半径最小处的方位角

4、P 缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不设缓和曲线时，

P=0

5、I 与缓和曲线最小半径端连接的圆曲线长度。

6、J缓和曲线小半径段的桩号。

再看看E匝道数据文件编辑

E匝

元设

据表

线元名

称

直线

缓和曲线

圆曲线

缓和曲线

缓和曲线

圆曲线

起点桩号 偏向

半径或

缓和曲

线旋转

A 常数

120

200

95

75

86.75

线元起点坐标、方位角

线元长

度(m)

71.296

72

60.070

45.125

64.841

17.800

道线

计数

X(m)

494360.619

494373.205

494381.631

494372.527

494355.878

494337.373

Y(m)

477706.293

477776.470

477847.870

477907.018

477948.932

478010.658

方位角

79°49′57.0″

79°49′57.0″

90°08′44.6″

107°21′16.3″

113°49′05.6″

92°24′19.0″

EK0+000.00

EK0+071.296

EK0+143.296

EK0+203.366

EK0+248.491

EK0+313.333

右偏

右偏

右偏

左偏

左偏

在E匝道线元表中，我们可以看出，E匝道式由一段直线段→缓和曲线段→圆曲线段→缓和曲线段→缓

和曲线段→圆曲线段连接起来的。值得注意的是，在位于K0+248.491处，是两段反向的缓和曲线连接点，

在编写数据程序文件时，分三段编写，第一段从K0+000到K0+203.366，第二段从K0+203.366到K0+248.919，

第三段从K0+248.919到线路终点K0+331.13。具体文件如下：（文件名“EZD”）

3→DimZ：X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else

“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494381.6316→X：477847.870→Y：200→R：

90°08°44.6°→W：120→P:60.07→I: 143.296→J：IfEnd: If L≥203.366：Then 494372.527→X：477948.932

→Y：200→R：107°21°16.3°→W：-95→P:60.07→I: 203.366→J：IfEnd: If L≥248.491：Then 494337.373

→X：478010.658→Y：86.75→R：92°24°19°→W：75→P:17.8→I: 313.333→J：IfEnd: Prog”FDCH”: Prog

“XLJS”：If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]¬－Y)sin(Ｏ)→N:L＋N→L:Int(1000L)÷

1000+0.0001→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=

“:L

◢

“K=”:K

◢

IfEnd:Goto 0: IfEnd

为了让使用者熟悉掌握数据程序文件的编写线元法的数据程序文件，现在将C、F、G等匝道的线元数据

列出来，供使用者练习。

C匝道数据

半径或

缓和曲

线旋转

常数A

32

65

140

65

85

85

85

线元起点坐标、方位角

线元长

度(m)

56.876

21.277

30.179

30.431

30.179

422.023

85

72.739

线元名

称

圆曲线

直线

缓和曲线

圆曲线

缓和曲线

直线

缓和曲线

圆曲线

起点桩号 偏向

X(m)

495005.991

494991. 262

494971. 511

494943.127

494913.167

494883.008

494461.552

494379.485

Y(m)

478056.677

478104. 123

478112.036

478122.242

478127.219

478126.740

478104.874

478086.681

方位角

56°19′42.3″

158°09′54.4″

158°09′54.4″

164°20′25.7″

176°47′40.6″

182.°58′11.8″

182.°58′11.8″

211°37′04.2″

CK0+000.00

CK0+056.876

CK0+078.153

CK0+108.332

CK0+138.763

CK0+168.941

CK0+590.964

CK0+675.964

右偏

右偏

右偏

右偏

右偏

右偏

右偏

F匝道数据

线元名

称

直线

缓和曲线

圆曲线

缓和曲线

起点桩号 偏向

半径或

缓和曲

线旋转

常数A

100

110

120

线元起点坐标、方位角

线元长

度(m)

111.525

90.909

101.114

130.909

X(m)

494591.311

494557.835

494542.812

494595.310

Y(m)

478420.033

478313.650

478224.689

478142.420

方位角

252°32′00.2″

252°32′00.2″

276°12′33.6″

328°52′35.7″

FK0+000.00

FK0+111.505

FK0+202.434

FK0+303.548

右偏

右偏

右偏

G匝道数据

缓和曲线 GK0+000.000

GK0+133.152

GK0+282.286

右偏

右偏

右偏

线元名

称

起点桩号 偏向

半径或

缓和曲

线旋转

常数A

175

230

190

线元起点坐标、方位角

线元长

度(m)

133.152

149.124

149.72

X(m)

494718.802

494587.600

494472.335

Y(m)

478112.963

478093.368

478002.890

方位角

182°58′11.8″

199°33′17.5″

236°42′21.1″

圆曲线

缓和曲线

三、道路计算程序

道路计算程序是核心程序，但是道路计算程序是固定的，使用者不需要做任何修改，也不直接运行道路

计算程序分为两部分，第一部分是常规的对成型线型组合计算的程序（文件名：XLJS），利用交点法编辑的

数据文件可以直接调用；第二部分是适用于非对称的线型组合计算文件（文件名：FDCH），在数据文件调

用完“FDCH”文件后，还要调用“XLJS”文件。

线路计算程序（XLJS）

Deg:Fix 4:G÷Abs(G)→V:F2÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:F÷2-F

^(3)÷（240R2）+F^（5）÷(34560R^（4）)-F^（7）÷(8386560R^（6）)+F^（9）÷(3158507520R^（8）)→Q:

Q+(R+P)Vtan(G÷2)→T:πRGV÷180+F→I:IF L≤H-T:Then L-H→J:Goto 1:IfEnd:If L

→Z:1→S:Goto 2:IfEnd:If L≤H-T+I-F:Then L-H+T→J:C→Z:1→S: Goto 3:IfEnd:If L

J:C+G→C:-1→S:Goto 2:If L≥H-T+I:Then L-H-I+2T→J:C+G→C:Goto 1:IfEnd:Lbl 1:D+Jcos C→X:E+Jsin C

→Y:C→O: Goto C:Lbl 2:(J-J^(5)÷40(RF)2+J^(9)÷(3456(RF)^(4))-J^(13)÷(599040(RF)^(6))+J^(17)÷(175472640

(RF)^(8))-T)S→M:(J^(3)÷(6RF)-J^(7)÷(336(RF)^(3))+J^(11)÷(42240(RF)^(5))-J^(15)÷(9676800(RF)^(7))+J^(19)÷

(3530097000(RF)^(9)))V→N:C+90J2SV÷πRF→O:Goto B

Lbl 3:180(J-0.5F) ÷πR→O:Q+Rsin (O)-T→M:(P+R(1-cos (O)))V→N:C+VO→O：Goto B

Lbl B:D+Mcos(C)-Nsin(C)→X:E+Msin(C)+Ncos(C)→Y :Goto C:Lbl C:If Z[3]≠0.0001:Then ?K: "M"?U:K<0

=>180-U →U:O+U×Abs(K)÷(K+10^(-10)) →Z:X+Abs(K)×cos(Z)→X

◢

Y+Abs(K)×sin(Z)→Y

◢

tan-1((Y-B)÷(X-A+

10

^(-10))→W

X-A<0=>W+180→W:W<0=>W+360→W:√((X-A)

2

+(Y-B)

2

)→S

"DIS": S

◢

"FW": W>DMS

◢

IfEnd

非对称线型组合计算程序（FDCH）

Abs(P+

10

^(-10)）÷(P+

10

^(-10))→S：Abs (R)÷R→V:Abs(R)→R:P

2

÷R→F:180(I+F)V÷(πR)→G:90FSV÷(π

R)→O:If S=1:Then W-O→C:Else If S=-1:Then W-O-G→C :IfEnd:IfEnd:If C<0:Then C=C+360→C:IfEndIf

C>360:Then C-360→C:IfEnd:F

2

÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→

P:F÷2-F^(3)÷(240R

2

)+F^5÷(34560R^(4))-F^(7)÷(8386560R^6)+F^9÷(3158507520R^(8))→

Q:Q+(R+P)Vtan(G÷2)→T:If S=1:Then J-F+T→H:C→Z:If S=-1:Then J-I-F+T→H:C+G→Z:IfEnd:IfEnd:(F-F^(3)

÷(40R

2

)+F^(5)÷(3456R^(4))-F^(7)÷(599040R^(6))+F^(9)÷(175472640R^(8))-T)S→M:(F

2

÷(6R)-F^(4)÷

(336R^(3))+F^(6)÷(42240R^(5))-F^(8)÷(9676800R^(7))+F^(10)÷(3530097000R^(9)))V→N:X-Mcos(Z)+Nsin(Z)

→D:Y-Msin(Z)-Ncos(Z)→E

四、程序执行

程序运行分为两种，第一种是根据根据使用者输入的桩号和距中距离计算出所求点的坐标和测站点到该

点的方位角和水平距离，可以直接利用计算器提供的距离和方位角进行极坐标放样。另一种是根据使用者输

入的实测坐标计算出该点的桩号和距中距离，可以根据桩号和距中距离进行高边坡、隧道断面点的设计高程、

以及进行路面高程控制计算方面很方便。

（一）、根据桩号和距中距离计算所求点的坐标、及方位角和水平距离。

要进行程序的线路计算，直接运行数程序文件就可以了，不需要理会计算程序。如果计算某桩号L距离

中线K处的坐标X，Y和该点到测站点的距离和方位角，运行提示如下：

X0? 提示输入测站点X坐标

Y0？ 提示输入测站点Y坐标

（测站点的位数和数据程序文件里坐标取位要一致）

Ln? 计算点近似桩号

近似桩号输入第四位小数后不能等于0.0001，

L？ 输入计算点准确桩号

近似桩号输入第四位小数后不能等于0.0001，

K？ 提示输入距离中桩的桩号

如果计算点位于线路的左幅，则在数据前加入“-”号，如果是计算中桩坐标，则输入0

M？ 输入该点与线路的前进方向的按顺时针转过的角度，常于斜交的涵洞，桥梁等构造物。

X 显示所求点的X坐标

Y 显示所求点的Y坐标

FW 显示所求点至测站点的方位角

DIS 显示所求点至测站点的距离

要继续计算则从L?开始，计算新点的坐标。

具体流程见程序流程图的左边部分。

例：要计算D匝道的DK0+020中桩坐标、DK0+080左边5.3米的坐标、Dk0+240斜交75°58′11.8″、距

中右侧5米处的坐标，并计算出该点至坐标为X=495005.991，Y=478056.677控制点的水平距离和方位角，

运行文件“DZD”提示如 下：

X0? 495005.991 测站X坐标

Y0？ 478056.677 测站Y坐标

Ln？ 0 输入小输点第四位起不等于0.0001的桩号

L？ 20 输入所求点桩号

K? 0 输入距中距离

M? 90 输入斜交角度

X= 494341.482 显示所求点X坐标

Y= 478007.903 显示所求点Y坐标

DIS 666.2959 显示所求点至测站平距

FW 184°11′52.42″显示所求点至测站方位角

L？ 80 输入下点的桩号

K？ -5.3

M？ 90

X= 494382.2290

Y= 477965.9594

DIS 630.324

FW 188°16′29.5″

L？ 240

K？ 5.0

M？ 75°58°11.8°

X= 494455.3110

Y= 478101.5491

DIS 552.505

FW 175°20′29.6″

（二）、根据坐标反算桩号和距中

如果是根据某点的坐标，计算出该点在线路中的桩号和距离，程序运行时提示如下：

X0? 提示输入测站点X坐标

Y0？ 提示输入测站点Y坐标

（测站点的位数和数据程序文件里坐标取位要一致）

Ln? 计算点近似桩号

近似桩号输入第四位小数后不能等于0.0001，

Xp？ 所求点P的X坐标

Yp？ 所求点P的Y坐标

L 显示所求点的桩号

K 显示所求点的距中距离，如果位于线路的左边，则显示结果为负值，如果K为0时，则

该点在线路中线上。

计算新点的桩号和距中距离，则从?Xp 开始。如此循环计算。

具体流程见流程图的右边部分

例：求点X=494382.2290，Y=477965.9594在D匝道中的位置，运行文件“DZD”如下：

X0? 495005.991 测站X坐标(可输入0)

Y0？ 478056.677 测站Y坐标(可输入0)

Ln？ 60.0001 输入小输点第四位起等于0.0001的桩号

Xp? 494382.2290 输入所求点X坐标

Yp？ 477965.9594 输入所求点Y坐标

L= 79.999999 显示所求点桩号

K= -5.2989 显示所求点距中距离, “-“ 表示该点位于中桩左

侧

五、程序应用

（一）本程序在高边坡放样中的应用

高边坡放样是目前公路道路测量工作普遍要做的一项工作，边坡坡顶线或者填方边坡坡脚线是边坡放样

最繁琐的工作，利用本程序进行边坡放样，工作步骤如下：

1、

2、

司镜员走到目的地后，先测出一个点的坐标和高程；

利用本程序求出该点在线路中的桩号和距中距离，然后根据计算出来的桩号和设计图纸，内插

出道路边桩高程，根据道路边桩高程和道路设计宽度、坡度、设计碎落平台的宽度和间隔高度、

实测点的高程计算出按实测高程，边坡距中距离。

3、 根据高程算出的距中距离和程序计算出立镜点至距中距离差值dK，指挥司镜员靠近道路中线或

者往外偏移多少米。

4、 重新做1~3工作，直到dK 小于边坡宽度控制精度要求为止。

dK

K1

h

h

m

m

K2

如上图，司镜员在P点立镜，测出P点坐标和高程后，用本程序计算出P点到的中桩桩号和距离

中桩距离K，根据中桩桩号和有关设计数据内插出边坡坡脚点的设计高程Hs和距离中桩的距离K1，

根据Hc、和Hs的高差计算出K2，按照P点的高程，边坡顶距中桩的理论宽度K

理论

，具体计算如下：

K

理论

=K1+h*n1+m+h*n2+m+（Hc-Hs-2h）*n3

则dK=K- K

理论

dK 求出来后，可以指挥司镜员靠近道路dK米，然后重复以上步骤，直到dK值能满足边坡控制精度。

（二）本程序在隧道断面欠挖、超挖放样中的应用。

利用本程序进行隧道断面开挖的欠挖和超挖计算，看下图：

R2'

R2

用全站仪测出P点的坐标和标高后，根据坐标算出P点的中桩桩号和距中距离K，根据中桩桩号

查看P点的路面高程，算出P点至设计路面的高差h，根据P点所在圆弧的圆心O2在断面中的位置（K2，

h2）算出P点到O2点距离R2’，用R2’-R2的差距就是P点的超挖或者欠挖的数据。

（三）本程序在路面施工控制放样中的应用

路面施工时，路面的轴线偏位、路面宽度、路面高程是测量的重要控制指标，也是道路施工控制最严格

的一道测量控制。利用本程序输入桩号和距中距离，可以直接计算出所求点至测站点的水平距离和方位角，

从而直接通过极坐标法直接放出该点在实地的位置。

（四）处理线路设计中的“断链”问题

断链是线路设计中经常出现的问题，断链分为长链与短链，如果是短链，那么只要弯道的计算范围从

短链处分界就可以了，如果是长链，则需要把线路分成两个文件，从开始出现断链处分，这样可以避免桩号

相同，但位置不同的情况。

（五）本程序在道路竣工验收过程中的应用

道路竣工验收，其主要方面就是轴线偏位和宽度，可以通过现场测量道路中线的坐标，通过坐标反算

出该点和设计中线的偏移值K（轴线偏位值）。

附：

附表一：CHLNR计算出来的逐桩坐标表

C匝道坐标一览表

桩号

CK0+000.00

CK0+020.00

CK0+035.00

CK0+060.00

CK0+145.00

CK0+170.00

CK0+190.00

CK0+210.00

CK0+230.00

CK0+250.00

CK0+270.00

CK0+290.00

X(m) Y(m) 桩号 X(m) Y(m) 桩号 X(m) Y(m)

495005.991 478056.677 CK0+310.00

495011.337 478075.612 CK0+330.00

495007.349 478089.930 CK0+350.00

494988.362 478105.285 CK0+370.00

494906.934 478127.439 CK0+390.00

494881.950 478126.685 CK0+410.00

494742.138 478119.431 CK0+550.00 494502.461 478106.996

494722.165 478118.395 CK0+570.00 494482.488 478105.960

494702.192 478117.358 CK0+590.00 494462.514 478104.924

494682.219 478116.322 CK0+610.00 494442.551 478103.728

494662.246 478115.286 CK0+630.00 494422.682 478101.484

494642.273 478114.250 CK0+650.00 494403.182 478097.112

494861.977 478125.648 CK0+430.00 494622.299 478113.214 CK0+670.00 494384.666 478089.632

494842.004 478124.612 CK0+450.00 494602.326 478112.177 CK0+690.00 494368.192 478078.372

494822.031 478123.576 CK0+470.00 494582.353 478111.141 CK0+710.00 494354.794 478063.585

494802.058 478122.540 CK0+490.00 494557.387 478109.846 CK0+730.00 494345.212 478046.082

494782.085 478121.503 CK0+510.00 494542.407 478109.069 CK0+748.70 494340.176 478028.110

494762.111 478120.467 CK0+530.00 494522.434 478108.032

D匝道坐标一览表

桩号

DK0+000.00

DK0+020.00

DK0+040.00

DK0+060.00

DK0+080.00

DK0+100.00

DK0+120.00

DK0+140.00

DK0+160.00

DK0+184.13

坐标

X

494341.902

494341.483

494348.316

494362.944

494382.173

494401.433

494417.398

494429.112

494437.236

494444.365

Y

478027.825

478007.908

477989.281

477975.948

477971.259

477975.999

477987.874

478004.018

478022.269

478045.316

桩号

坐标

XY

桩号

坐标

XY

E匝道坐标一览表

桩号

坐标

X

Y

桩号

坐标

X

Y

桩号

坐标

X

Y

EK0+000.00

494360.619

477706.293

EK0+200.00

494373.504

477903.797

EK0+020.00

494364.149

477725.979

EK0+220.00

494366.991

477922.7

EK0+040.00

494367.68

477745.665

EK0+240.00

494359.296

477941.159

EK0+060.00

494371.21

477765.351

EK0+260.00

494351.272

477959.479

EK0+080.00

494374.734

477785.038

EK0+280.00

494344.011

477978.109

EK0+100.00

494378.002

477804.768

EK0+305.00

494338.104

478002.36

EK0+120.00

494380.481

477824.612

EK0+331.13

494338.449

478028.394

EK0+140.00

494381.613

477844.574

EK0+160.00

494380.892

477864.553

EK0+180.00

494378.18

477884.36

F匝道坐标一览表

桩号

FK0+000.00

FK0+020.00

FK0+040.00

FK0+060.00

FK0+080.00

FK0+100.00

FK0+120.00

FK0+140.00

FK0+160.00

FK0+180.00

FK0+200.00

FK0+220.00

坐标

X

494591.311

494585.308

494579.305

494573.302

494567.299

494561.296

494555.302

494549.658

494545.116

494542.481

494542.575

494546.095

Y

478420.033

478400.955

478381.878

478362.8

478343.722

478324.644

478305.563

478286.377

478266.905

478247.091

478227.112

478207.452

桩号

FK0+240.00

FK0+260.00

FK0+280.00

FK0+300.00

FK0+320.00

FK0+340.00

FK0+360.00

FK0+380.00

FK0+400.00

FK0+420.00

FK0+434.46

坐标

X

494553.112

494563.394

494576.602

494592.302

494609.954

494628.87

494650.453

494668.378

494688.371

494708.361

494722.801

Y

478188.752

478171.63

478156.649

478144.303

478134.953

478128.512

478124.33

478122.732

478122.373

478122.971

478123.685

桩号

坐标

XY

G匝道坐标一览表

桩号

坐标

XY

494718.802

478112.963

494698.831

478111.883

494678.876

478110.543

494656.977

478108.457

494639.141

478106.042

494619.483

478102.371

494600.109

478097.425

494576.543

478089.117

494558.422

478080.668

494541.103

478070.677

494524.718

478059.22

桩号

GK0+230.00

GK0+255.00

GK0+270.00

GK0+285.00

GK0+305.00

GK0+330.00

GK0+350.00

GK0+370.00

GK0+390.00

GK0+410.00

GK0+432.01

坐标

XY

494505.736

478042.969

494488.629

478024.756

494479.35

478012.974

494470.858

478000.612

494460.776

477983.346

494449.997

477960.797

494442.64

477942.202

494436.211

477923.265

494430.512

477904.096

494425.338

477884.776

494420

477863.424

桩号

坐标

XY

GK0+000.00

GK0+020.00

GK0+040.00

GK0+062.00

GK0+080.00

GK0+100.00

GK0+120.00

GK0+145.00

GK0+165.00

GK0+185.00

GK0+205.00

六、公路纵断面高程计算

公路纵断面计算也是采用计算程序和数据程序分开的原则，数据文件是专门存放道路纵断面设计参数，

对于多个工程，可以建立多个数据文件，这样每个数据文件通过调用计算程序就可以计算出该点的高程。

（一） 数据文件

每个纵断面由纵坡和竖曲线组成，本程序每一个变坡点计算范围是从竖曲线前的纵坡任意桩号开始到

竖曲线后纵坡上的任一桩号结束。

下表是某线路路基纵断面设计的参数一览表：

根据上表编程数据文件如下：文件名BBR

从上面程序看， 这个语句是固定的，“Dh”提

醒使用者输入所求点距离路面的厚度，比如路面底基层距离路面厚度为37cm，则输入0.37，表示所求点是路

面底基层顶面高程，“L=”，输入所求点桩号,”I”=表示输入路面横坡值，“K=”表示输入路面某点至中桩距

离，均为正值。以上数据均在程序运行时输入。

在每一个变坡点处，要求输入的参数有：

C；输入变坡点高程

D：输入边坡点桩号

M：输入变坡点小里程段纵坡坡度

N：输入变坡点大里程侧纵坡坡度

R：竖曲线半径

L>表示竖曲线前纵坡任意里程的桩号，本程序每段竖曲线参数计算范围是从竖曲线前的纵坡上任意点

到竖曲线后纵坡上任意点的桩号。

（二） 高程计算程序（文件名：GCJS）

高程计算文件代码如下：

（三） 运行

运行BBR程序显示如下

Dh？ 提示输入至路面厚度

L=? 提示输入桩号

i？ 提示输入横坡度

K？ 提示输入距中距离，均为正数

H= 提示输出所求点高程

例：根据上表计算K38+660中桩高程、K38+900右侧4.5米高程、K39+600底基层面左侧4.5

米高程（假设路面沥青厚度为6cm，封釉层1cm，水泥稳定层基层厚度为30cm）

运行BBR文件

Dh？ 0

L=？ 38660

i？ 0

K？ 0

H= 76.9412

L=？ 38900

i？ 0.02

K？ 4.50

H= 85.2200

退出程序，重新运行

Dh？ 0.370

L=？ 39600

i？ 0.02

K？ 4.50

H= 93.0509

2024年3月29日发(作者：建瑶岑)

5800公路测量程序使用说明

一、程序使用流程

本程序数据和主程序是分开的，编程时将不同的工程数据存放到不同的数据文件里，如A匝道，文件

名为A，将匝道A所有的曲线线元参数输入A文件里。运行时只要运行文件名A的程序就可以了，具体运

行流程见下图：

二、数据文件的编写

（一）交点法数据文件编辑

交点法编写数据文件必须是对称型的，即直线段→缓和曲线段→圆曲线段→缓和曲线段→直线段，（如

果任意一端没有直线段，则把直线段长度看做是0），另外圆曲线两侧缓和曲线的旋转常数必须相等，并且

和直线段连接处的半径必须是无穷大。

交点法数据文件编写一般是根据设计图纸提供的平面曲线参数一览表提供的参数来编写，每个弯道包

括：弯道起点方位角（C），交点X坐标（D），交点Y坐标（E），缓和曲线长度（F，当没有设缓和曲线

时，F=0），交点转交（G，向左转弯，G为负值，向右转弯，G取正值），交点桩号（H），弯道圆曲线半

径（R）。

下图是一段市政道路设计参数数据。

根据上图提供的数据，可以编辑成如下的数据文件：

文件名:CHLNR

3→DimZ “X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Els

e “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L>0 :Then 98°39°35.12°→C：4474.384→D：2415

.861→E：140→F：31°17°23°→G：410.007→H:600→R:IfEnd: If L>1060 Then 129°56°58.19°→C:

4206.421→D：3093.946→E：70→F：-33°50°48°→G：1285.437→H:600→R:IfEnd:Prog”XLJS”:If Z[3]≠0.

0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（

O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L

◢

“K=”:K

◢

IfEnd:Goto 0: IfEnd

在面程式中，有两个条件转移语句即If L>0：Then 98°39°35.12°→C:4774.384→D: 2415.861→E:140→

F:31°17°23°→G:410.007→H:600→R:IfEnd

If L>1060：Then 129°56°58.19°→C:4206.421→D: 3093.946→E:70→F:-33°50°48°→G:1285.437→H:600

→R:IfEnd

……

如果还有其他弯道，可以继续完后加。在这些存放设计参数的语句前后的程序表达式是固定的。

说明：

编辑曲线参数时，每个曲线参数放在一个If L>***.***（两个弯道中间直线段上的任意桩号） :Then ***°**°

**.**°→C（弯道起点方位角）：****.***→D（交点X坐标）：****.****→E（交点Y坐标）：***.***→F（缓

和曲线长度：**°**°**°→G（转角，向左转为负值，向右转为正值）：***.***→H（交点里程桩号）:***

→R（圆曲线半径）:IfEnd 条件式语句里，如果有多个弯道，一直按上述形式编下去，

变量说明：

1、 L>***.***

弯道参数的起点范围，可以是弯道起点，也可以是弯道前直线段里的任意点桩号

2、 ***°**°**.**°→C 弯道起点方位角

3、****.***→D、****.****→E

交点X坐标、Y坐标

4、***.***→F

缓和曲线长度，如果没设置缓和曲线，则缓和曲线长度看做0

5、**°**°**°→G 转角，向左转为负值，向右转为正值

6、***.***→H、交点里程桩号

7、***→R 圆曲线半径

注意，运行程序时，第一部是输入测站坐标X0，Y0，然后输入近似桩号，主要是为了加快收敛，输入

时注意当根据桩号和距中距离计算所求点坐标时，请不要把小数点第四位输入0.0001，如K21+369.0921，请

把桩号输入21369.092即可，相反，当根据坐标反算桩号和距中距离时，请在输入时把第四位小数数位0.0001，

如K21+200，输入时请输入21200.0001。

（二）线元法数据文件编辑

线元法数据文件编辑适合于各种各样的线型组合，特别是在大型立交区或者是不规则的线型组合地方，

如小区道路、重丘区山路等有特别明显的有点，本程序为了减少使用者设计参数的输入，特别把直线与缓和

曲线以及圆曲线连接的地方只要输入缓和曲线参数就可以了，这样可以减少输入很多设计参数，同时节约大

量的内存。看下图，

该图是一个半喇叭型的互通立交，一共是5条匝道，我们以D匝道为例，来看看线元法数据文件的建立。

下表是D匝道的线元表

：

线元名称 起点桩号 偏向

半径或

缓和曲

线旋转

常数A

55

40

65

线元起点坐标、方位角

线元长度

(m)

39.289

39.216

105.625

X(m)

494341.902

494347.926

494380.677

Y(m)

478027.825

477989. 873

477971.271

方位角

260°38′55.7″

302°18′26.0″

358°28′45.9″

缓和曲线

圆曲线

缓和曲线

DK0+000.00

DK0+039.289

DK0+078.504

右偏

右偏

右偏

D匝道线元数据分析：从上表中可以看出，D匝道由三个线元组成，第一段是缓和曲线段，根据上表判

断，该缓和曲线属于第一段缓和曲线，其最小半径处桩号为J=K0+039.289，该处坐标X=494347.926,Y=477989.

873，半径为R=40米，方位角W=302°18′26.0″，缓和曲线旋转常数P=55，在缓和曲线终点K0+039.289处，

有一圆曲线与其相连接，长度I=39.216。第二段是圆曲线，第三段又是一段缓和曲线，根据上表判断，该曲线

属于第二缓和曲线，曲线最小半径位于与园曲线连接处，桩号J=K0+78.504，坐标为X=494380.677，Y=477971.271,

方位角W=358°28′45.9″，最小半径R=40，缓和曲线旋转常数P=-65（为了判断是第一段还是第二缓和曲线，

我们通过P的正负来判断，第一缓和段P统一用正数表示，第二缓和段统一用负号表示。起点处相连的圆曲线

长度I可以看做0，也可以看做是39.216.由此我们整理下。从整个线型全部是向右转向，所以R全部取正数。

具体如下：

第一段缓和曲线线元参数：

最小半径处坐标X=494347.926 Y=477989. 873

最小半径R=40

最小半径处方位角：W=302°18′26.0″

缓和曲线旋转常数P=55

与之相连的圆曲线长度I=39.216

最小半径处桩号J=39.289

此段数据计算范围可以从K0+000~ K0+78.504

第二段缓和曲线线元参数

最小半径处坐标X=494380.677 Y=477971.271

最小半径R=40

最小半径处方位角：W=358°28′45.9″

缓和曲线旋转常数P=65

与之相连的圆曲线长度I=39.216

最小半径处桩号J= W=358°28′45.9″

此段数据计算范围可以从K0+78.504到匝道结束。

当桩号L位于K0+000~ K0+78.504时，曲线参数按第一段编写，当L>78.504时，按第二段缓和曲线编

写，当然第二段缓和曲线计算范围也可以从圆曲线起点K0+039.289到K0+78.504中间任意桩号开始。下面

是D匝道的数据程序文件（文件名“DZD”）：

3→DimZ：X0”?A:”Y0”?B:“Ln”?L: Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L:

Else “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494347.926→X：477989. 873→Y：40

→R：302°18°26°→W：55→P:39.216→I: 39.289→J：IfEnd: If L≥78.504：Then 494380.677→X：4779

71.271→Y：40→R：358°28°45.9°→W：-65→P:39.216→I: 78.504→J：IfEnd:Prog”FDCH”: Prog”XLJS”:

If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]

－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L

◢

“K=”:K

◢

IfEnd:Goto 0: IfEnd

在上面程式中，每个缓和曲线的线元参数放在一个If … Then …… IfEnd循环语句中，If 后面紧跟计算的

起点范围，终点范围从下段的起点范围结束。

变量定义：

1、X、Y 输入缓和曲线半径最小处的X、Y坐标

2、R 缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。

3、W 输入缓和曲线半径最小处的方位角

4、P 缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不设缓和曲线时，

P=0

5、I 与缓和曲线最小半径端连接的圆曲线长度。

6、J缓和曲线小半径段的桩号。

再看看E匝道数据文件编辑

E匝

元设

据表

线元名

称

直线

缓和曲线

圆曲线

缓和曲线

缓和曲线

圆曲线

起点桩号 偏向

半径或

缓和曲

线旋转

A 常数

120

200

95

75

86.75

线元起点坐标、方位角

线元长

度(m)

71.296

72

60.070

45.125

64.841

17.800

道线

计数

X(m)

494360.619

494373.205

494381.631

494372.527

494355.878

494337.373

Y(m)

477706.293

477776.470

477847.870

477907.018

477948.932

478010.658

方位角

79°49′57.0″

79°49′57.0″

90°08′44.6″

107°21′16.3″

113°49′05.6″

92°24′19.0″

EK0+000.00

EK0+071.296

EK0+143.296

EK0+203.366

EK0+248.491

EK0+313.333

右偏

右偏

右偏

左偏

左偏

在E匝道线元表中，我们可以看出，E匝道式由一段直线段→缓和曲线段→圆曲线段→缓和曲线段→缓

和曲线段→圆曲线段连接起来的。值得注意的是，在位于K0+248.491处，是两段反向的缓和曲线连接点，

在编写数据程序文件时，分三段编写，第一段从K0+000到K0+203.366，第二段从K0+203.366到K0+248.919，

第三段从K0+248.919到线路终点K0+331.13。具体文件如下：（文件名“EZD”）

3→DimZ：X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else

“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494381.6316→X：477847.870→Y：200→R：

90°08°44.6°→W：120→P:60.07→I: 143.296→J：IfEnd: If L≥203.366：Then 494372.527→X：477948.932

→Y：200→R：107°21°16.3°→W：-95→P:60.07→I: 203.366→J：IfEnd: If L≥248.491：Then 494337.373

→X：478010.658→Y：86.75→R：92°24°19°→W：75→P:17.8→I: 313.333→J：IfEnd: Prog”FDCH”: Prog

“XLJS”：If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]¬－Y)sin(Ｏ)→N:L＋N→L:Int(1000L)÷

1000+0.0001→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=

“:L

◢

“K=”:K

◢

IfEnd:Goto 0: IfEnd

为了让使用者熟悉掌握数据程序文件的编写线元法的数据程序文件，现在将C、F、G等匝道的线元数据

列出来，供使用者练习。

C匝道数据

半径或

缓和曲

线旋转

常数A

32

65

140

65

85

85

85

线元起点坐标、方位角

线元长

度(m)

56.876

21.277

30.179

30.431

30.179

422.023

85

72.739

线元名

称

圆曲线

直线

缓和曲线

圆曲线

缓和曲线

直线

缓和曲线

圆曲线

起点桩号 偏向

X(m)

495005.991

494991. 262

494971. 511

494943.127

494913.167

494883.008

494461.552

494379.485

Y(m)

478056.677

478104. 123

478112.036

478122.242

478127.219

478126.740

478104.874

478086.681

方位角

56°19′42.3″

158°09′54.4″

158°09′54.4″

164°20′25.7″

176°47′40.6″

182.°58′11.8″

182.°58′11.8″

211°37′04.2″

CK0+000.00

CK0+056.876

CK0+078.153

CK0+108.332

CK0+138.763

CK0+168.941

CK0+590.964

CK0+675.964

右偏

右偏

右偏

右偏

右偏

右偏

右偏

F匝道数据

线元名

称

直线

缓和曲线

圆曲线

缓和曲线

起点桩号 偏向

半径或

缓和曲

线旋转

常数A

100

110

120

线元起点坐标、方位角

线元长

度(m)

111.525

90.909

101.114

130.909

X(m)

494591.311

494557.835

494542.812

494595.310

Y(m)

478420.033

478313.650

478224.689

478142.420

方位角

252°32′00.2″

252°32′00.2″

276°12′33.6″

328°52′35.7″

FK0+000.00

FK0+111.505

FK0+202.434

FK0+303.548

右偏

右偏

右偏

G匝道数据

缓和曲线 GK0+000.000

GK0+133.152

GK0+282.286

右偏

右偏

右偏

线元名

称

起点桩号 偏向

半径或

缓和曲

线旋转

常数A

175

230

190

线元起点坐标、方位角

线元长

度(m)

133.152

149.124

149.72

X(m)

494718.802

494587.600

494472.335

Y(m)

478112.963

478093.368

478002.890

方位角

182°58′11.8″

199°33′17.5″

236°42′21.1″

圆曲线

缓和曲线

三、道路计算程序

道路计算程序是核心程序，但是道路计算程序是固定的，使用者不需要做任何修改，也不直接运行道路

计算程序分为两部分，第一部分是常规的对成型线型组合计算的程序（文件名：XLJS），利用交点法编辑的

数据文件可以直接调用；第二部分是适用于非对称的线型组合计算文件（文件名：FDCH），在数据文件调

用完“FDCH”文件后，还要调用“XLJS”文件。

线路计算程序（XLJS）

Deg:Fix 4:G÷Abs(G)→V:F2÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:F÷2-F

^(3)÷（240R2）+F^（5）÷(34560R^（4）)-F^（7）÷(8386560R^（6）)+F^（9）÷(3158507520R^（8）)→Q:

Q+(R+P)Vtan(G÷2)→T:πRGV÷180+F→I:IF L≤H-T:Then L-H→J:Goto 1:IfEnd:If L

→Z:1→S:Goto 2:IfEnd:If L≤H-T+I-F:Then L-H+T→J:C→Z:1→S: Goto 3:IfEnd:If L

J:C+G→C:-1→S:Goto 2:If L≥H-T+I:Then L-H-I+2T→J:C+G→C:Goto 1:IfEnd:Lbl 1:D+Jcos C→X:E+Jsin C

→Y:C→O: Goto C:Lbl 2:(J-J^(5)÷40(RF)2+J^(9)÷(3456(RF)^(4))-J^(13)÷(599040(RF)^(6))+J^(17)÷(175472640

(RF)^(8))-T)S→M:(J^(3)÷(6RF)-J^(7)÷(336(RF)^(3))+J^(11)÷(42240(RF)^(5))-J^(15)÷(9676800(RF)^(7))+J^(19)÷

(3530097000(RF)^(9)))V→N:C+90J2SV÷πRF→O:Goto B

Lbl 3:180(J-0.5F) ÷πR→O:Q+Rsin (O)-T→M:(P+R(1-cos (O)))V→N:C+VO→O：Goto B

Lbl B:D+Mcos(C)-Nsin(C)→X:E+Msin(C)+Ncos(C)→Y :Goto C:Lbl C:If Z[3]≠0.0001:Then ?K: "M"?U:K<0

=>180-U →U:O+U×Abs(K)÷(K+10^(-10)) →Z:X+Abs(K)×cos(Z)→X

◢

Y+Abs(K)×sin(Z)→Y

◢

tan-1((Y-B)÷(X-A+

10

^(-10))→W

X-A<0=>W+180→W:W<0=>W+360→W:√((X-A)

2

+(Y-B)

2

)→S

"DIS": S

◢

"FW": W>DMS

◢

IfEnd

非对称线型组合计算程序（FDCH）

Abs(P+

10

^(-10)）÷(P+

10

^(-10))→S：Abs (R)÷R→V:Abs(R)→R:P

2

÷R→F:180(I+F)V÷(πR)→G:90FSV÷(π

R)→O:If S=1:Then W-O→C:Else If S=-1:Then W-O-G→C :IfEnd:IfEnd:If C<0:Then C=C+360→C:IfEndIf

C>360:Then C-360→C:IfEnd:F

2

÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→

P:F÷2-F^(3)÷(240R

2

)+F^5÷(34560R^(4))-F^(7)÷(8386560R^6)+F^9÷(3158507520R^(8))→

Q:Q+(R+P)Vtan(G÷2)→T:If S=1:Then J-F+T→H:C→Z:If S=-1:Then J-I-F+T→H:C+G→Z:IfEnd:IfEnd:(F-F^(3)

÷(40R

2

)+F^(5)÷(3456R^(4))-F^(7)÷(599040R^(6))+F^(9)÷(175472640R^(8))-T)S→M:(F

2

÷(6R)-F^(4)÷

(336R^(3))+F^(6)÷(42240R^(5))-F^(8)÷(9676800R^(7))+F^(10)÷(3530097000R^(9)))V→N:X-Mcos(Z)+Nsin(Z)

→D:Y-Msin(Z)-Ncos(Z)→E

四、程序执行

程序运行分为两种，第一种是根据根据使用者输入的桩号和距中距离计算出所求点的坐标和测站点到该

点的方位角和水平距离，可以直接利用计算器提供的距离和方位角进行极坐标放样。另一种是根据使用者输

入的实测坐标计算出该点的桩号和距中距离，可以根据桩号和距中距离进行高边坡、隧道断面点的设计高程、

以及进行路面高程控制计算方面很方便。

（一）、根据桩号和距中距离计算所求点的坐标、及方位角和水平距离。

要进行程序的线路计算，直接运行数程序文件就可以了，不需要理会计算程序。如果计算某桩号L距离

中线K处的坐标X，Y和该点到测站点的距离和方位角，运行提示如下：

X0? 提示输入测站点X坐标

Y0？ 提示输入测站点Y坐标

（测站点的位数和数据程序文件里坐标取位要一致）

Ln? 计算点近似桩号

近似桩号输入第四位小数后不能等于0.0001，

L？ 输入计算点准确桩号

近似桩号输入第四位小数后不能等于0.0001，

K？ 提示输入距离中桩的桩号

如果计算点位于线路的左幅，则在数据前加入“-”号，如果是计算中桩坐标，则输入0

M？ 输入该点与线路的前进方向的按顺时针转过的角度，常于斜交的涵洞，桥梁等构造物。

X 显示所求点的X坐标

Y 显示所求点的Y坐标

FW 显示所求点至测站点的方位角

DIS 显示所求点至测站点的距离

要继续计算则从L?开始，计算新点的坐标。

具体流程见程序流程图的左边部分。

例：要计算D匝道的DK0+020中桩坐标、DK0+080左边5.3米的坐标、Dk0+240斜交75°58′11.8″、距

中右侧5米处的坐标，并计算出该点至坐标为X=495005.991，Y=478056.677控制点的水平距离和方位角，

运行文件“DZD”提示如 下：

X0? 495005.991 测站X坐标

Y0？ 478056.677 测站Y坐标

Ln？ 0 输入小输点第四位起不等于0.0001的桩号

L？ 20 输入所求点桩号

K? 0 输入距中距离

M? 90 输入斜交角度

X= 494341.482 显示所求点X坐标

Y= 478007.903 显示所求点Y坐标

DIS 666.2959 显示所求点至测站平距

FW 184°11′52.42″显示所求点至测站方位角

L？ 80 输入下点的桩号

K？ -5.3

M？ 90

X= 494382.2290

Y= 477965.9594

DIS 630.324

FW 188°16′29.5″

L？ 240

K？ 5.0

M？ 75°58°11.8°

X= 494455.3110

Y= 478101.5491

DIS 552.505

FW 175°20′29.6″

（二）、根据坐标反算桩号和距中

如果是根据某点的坐标，计算出该点在线路中的桩号和距离，程序运行时提示如下：

X0? 提示输入测站点X坐标

Y0？ 提示输入测站点Y坐标

（测站点的位数和数据程序文件里坐标取位要一致）

Ln? 计算点近似桩号

近似桩号输入第四位小数后不能等于0.0001，

Xp？ 所求点P的X坐标

Yp？ 所求点P的Y坐标

L 显示所求点的桩号

K 显示所求点的距中距离，如果位于线路的左边，则显示结果为负值，如果K为0时，则

该点在线路中线上。

计算新点的桩号和距中距离，则从?Xp 开始。如此循环计算。

具体流程见流程图的右边部分

例：求点X=494382.2290，Y=477965.9594在D匝道中的位置，运行文件“DZD”如下：

X0? 495005.991 测站X坐标(可输入0)

Y0？ 478056.677 测站Y坐标(可输入0)

Ln？ 60.0001 输入小输点第四位起等于0.0001的桩号

Xp? 494382.2290 输入所求点X坐标

Yp？ 477965.9594 输入所求点Y坐标

L= 79.999999 显示所求点桩号

K= -5.2989 显示所求点距中距离, “-“ 表示该点位于中桩左

侧

五、程序应用

（一）本程序在高边坡放样中的应用

高边坡放样是目前公路道路测量工作普遍要做的一项工作，边坡坡顶线或者填方边坡坡脚线是边坡放样

最繁琐的工作，利用本程序进行边坡放样，工作步骤如下：

1、

2、

司镜员走到目的地后，先测出一个点的坐标和高程；

利用本程序求出该点在线路中的桩号和距中距离，然后根据计算出来的桩号和设计图纸，内插

出道路边桩高程，根据道路边桩高程和道路设计宽度、坡度、设计碎落平台的宽度和间隔高度、

实测点的高程计算出按实测高程，边坡距中距离。

3、 根据高程算出的距中距离和程序计算出立镜点至距中距离差值dK，指挥司镜员靠近道路中线或

者往外偏移多少米。

4、 重新做1~3工作，直到dK 小于边坡宽度控制精度要求为止。

dK

K1

h

h

m

m

K2

如上图，司镜员在P点立镜，测出P点坐标和高程后，用本程序计算出P点到的中桩桩号和距离

中桩距离K，根据中桩桩号和有关设计数据内插出边坡坡脚点的设计高程Hs和距离中桩的距离K1，

根据Hc、和Hs的高差计算出K2，按照P点的高程，边坡顶距中桩的理论宽度K

理论

，具体计算如下：

K

理论

=K1+h*n1+m+h*n2+m+（Hc-Hs-2h）*n3

则dK=K- K

理论

dK 求出来后，可以指挥司镜员靠近道路dK米，然后重复以上步骤，直到dK值能满足边坡控制精度。

（二）本程序在隧道断面欠挖、超挖放样中的应用。

利用本程序进行隧道断面开挖的欠挖和超挖计算，看下图：

R2'

R2

用全站仪测出P点的坐标和标高后，根据坐标算出P点的中桩桩号和距中距离K，根据中桩桩号

查看P点的路面高程，算出P点至设计路面的高差h，根据P点所在圆弧的圆心O2在断面中的位置（K2，

h2）算出P点到O2点距离R2’，用R2’-R2的差距就是P点的超挖或者欠挖的数据。

（三）本程序在路面施工控制放样中的应用

路面施工时，路面的轴线偏位、路面宽度、路面高程是测量的重要控制指标，也是道路施工控制最严格

的一道测量控制。利用本程序输入桩号和距中距离，可以直接计算出所求点至测站点的水平距离和方位角，

从而直接通过极坐标法直接放出该点在实地的位置。

（四）处理线路设计中的“断链”问题

断链是线路设计中经常出现的问题，断链分为长链与短链，如果是短链，那么只要弯道的计算范围从

短链处分界就可以了，如果是长链，则需要把线路分成两个文件，从开始出现断链处分，这样可以避免桩号

相同，但位置不同的情况。

（五）本程序在道路竣工验收过程中的应用

道路竣工验收，其主要方面就是轴线偏位和宽度，可以通过现场测量道路中线的坐标，通过坐标反算

出该点和设计中线的偏移值K（轴线偏位值）。

附：

附表一：CHLNR计算出来的逐桩坐标表

C匝道坐标一览表

桩号

CK0+000.00

CK0+020.00

CK0+035.00

CK0+060.00

CK0+145.00

CK0+170.00

CK0+190.00

CK0+210.00

CK0+230.00

CK0+250.00

CK0+270.00

CK0+290.00

X(m) Y(m) 桩号 X(m) Y(m) 桩号 X(m) Y(m)

495005.991 478056.677 CK0+310.00

495011.337 478075.612 CK0+330.00

495007.349 478089.930 CK0+350.00

494988.362 478105.285 CK0+370.00

494906.934 478127.439 CK0+390.00

494881.950 478126.685 CK0+410.00

494742.138 478119.431 CK0+550.00 494502.461 478106.996

494722.165 478118.395 CK0+570.00 494482.488 478105.960

494702.192 478117.358 CK0+590.00 494462.514 478104.924

494682.219 478116.322 CK0+610.00 494442.551 478103.728

494662.246 478115.286 CK0+630.00 494422.682 478101.484

494642.273 478114.250 CK0+650.00 494403.182 478097.112

494861.977 478125.648 CK0+430.00 494622.299 478113.214 CK0+670.00 494384.666 478089.632

494842.004 478124.612 CK0+450.00 494602.326 478112.177 CK0+690.00 494368.192 478078.372

494822.031 478123.576 CK0+470.00 494582.353 478111.141 CK0+710.00 494354.794 478063.585

494802.058 478122.540 CK0+490.00 494557.387 478109.846 CK0+730.00 494345.212 478046.082

494782.085 478121.503 CK0+510.00 494542.407 478109.069 CK0+748.70 494340.176 478028.110

494762.111 478120.467 CK0+530.00 494522.434 478108.032

D匝道坐标一览表

桩号

DK0+000.00

DK0+020.00

DK0+040.00

DK0+060.00

DK0+080.00

DK0+100.00

DK0+120.00

DK0+140.00

DK0+160.00

DK0+184.13

坐标

X

494341.902

494341.483

494348.316

494362.944

494382.173

494401.433

494417.398

494429.112

494437.236

494444.365

Y

478027.825

478007.908

477989.281

477975.948

477971.259

477975.999

477987.874

478004.018

478022.269

478045.316

桩号

坐标

XY

桩号

坐标

XY

E匝道坐标一览表

桩号

坐标

X

Y

桩号

坐标

X

Y

桩号

坐标

X

Y

EK0+000.00

494360.619

477706.293

EK0+200.00

494373.504

477903.797

EK0+020.00

494364.149

477725.979

EK0+220.00

494366.991

477922.7

EK0+040.00

494367.68

477745.665

EK0+240.00

494359.296

477941.159

EK0+060.00

494371.21

477765.351

EK0+260.00

494351.272

477959.479

EK0+080.00

494374.734

477785.038

EK0+280.00

494344.011

477978.109

EK0+100.00

494378.002

477804.768

EK0+305.00

494338.104

478002.36

EK0+120.00

494380.481

477824.612

EK0+331.13

494338.449

478028.394

EK0+140.00

494381.613

477844.574

EK0+160.00

494380.892

477864.553

EK0+180.00

494378.18

477884.36

F匝道坐标一览表

桩号

FK0+000.00

FK0+020.00

FK0+040.00

FK0+060.00

FK0+080.00

FK0+100.00

FK0+120.00

FK0+140.00

FK0+160.00

FK0+180.00

FK0+200.00

FK0+220.00

坐标

X

494591.311

494585.308

494579.305

494573.302

494567.299

494561.296

494555.302

494549.658

494545.116

494542.481

494542.575

494546.095

Y

478420.033

478400.955

478381.878

478362.8

478343.722

478324.644

478305.563

478286.377

478266.905

478247.091

478227.112

478207.452

桩号

FK0+240.00

FK0+260.00

FK0+280.00

FK0+300.00

FK0+320.00

FK0+340.00

FK0+360.00

FK0+380.00

FK0+400.00

FK0+420.00

FK0+434.46

坐标

X

494553.112

494563.394

494576.602

494592.302

494609.954

494628.87

494650.453

494668.378

494688.371

494708.361

494722.801

Y

478188.752

478171.63

478156.649

478144.303

478134.953

478128.512

478124.33

478122.732

478122.373

478122.971

478123.685

桩号

坐标

XY

G匝道坐标一览表

桩号

坐标

XY

494718.802

478112.963

494698.831

478111.883

494678.876

478110.543

494656.977

478108.457

494639.141

478106.042

494619.483

478102.371

494600.109

478097.425

494576.543

478089.117

494558.422

478080.668

494541.103

478070.677

494524.718

478059.22

桩号

GK0+230.00

GK0+255.00

GK0+270.00

GK0+285.00

GK0+305.00

GK0+330.00

GK0+350.00

GK0+370.00

GK0+390.00

GK0+410.00

GK0+432.01

坐标

XY

494505.736

478042.969

494488.629

478024.756

494479.35

478012.974

494470.858

478000.612

494460.776

477983.346

494449.997

477960.797

494442.64

477942.202

494436.211

477923.265

494430.512

477904.096

494425.338

477884.776

494420

477863.424

桩号

坐标

XY

GK0+000.00

GK0+020.00

GK0+040.00

GK0+062.00

GK0+080.00

GK0+100.00

GK0+120.00

GK0+145.00

GK0+165.00

GK0+185.00

GK0+205.00

六、公路纵断面高程计算

公路纵断面计算也是采用计算程序和数据程序分开的原则，数据文件是专门存放道路纵断面设计参数，

对于多个工程，可以建立多个数据文件，这样每个数据文件通过调用计算程序就可以计算出该点的高程。

（一） 数据文件

每个纵断面由纵坡和竖曲线组成，本程序每一个变坡点计算范围是从竖曲线前的纵坡任意桩号开始到

竖曲线后纵坡上的任一桩号结束。

下表是某线路路基纵断面设计的参数一览表：

根据上表编程数据文件如下：文件名BBR

从上面程序看， 这个语句是固定的，“Dh”提

醒使用者输入所求点距离路面的厚度，比如路面底基层距离路面厚度为37cm，则输入0.37，表示所求点是路

面底基层顶面高程，“L=”，输入所求点桩号,”I”=表示输入路面横坡值，“K=”表示输入路面某点至中桩距

离，均为正值。以上数据均在程序运行时输入。

在每一个变坡点处，要求输入的参数有：

C；输入变坡点高程

D：输入边坡点桩号

M：输入变坡点小里程段纵坡坡度

N：输入变坡点大里程侧纵坡坡度

R：竖曲线半径

L>表示竖曲线前纵坡任意里程的桩号，本程序每段竖曲线参数计算范围是从竖曲线前的纵坡上任意点

到竖曲线后纵坡上任意点的桩号。

（二） 高程计算程序（文件名：GCJS）

高程计算文件代码如下：

（三） 运行

运行BBR程序显示如下

Dh？ 提示输入至路面厚度

L=? 提示输入桩号

i？ 提示输入横坡度

K？ 提示输入距中距离，均为正数

H= 提示输出所求点高程

例：根据上表计算K38+660中桩高程、K38+900右侧4.5米高程、K39+600底基层面左侧4.5

米高程（假设路面沥青厚度为6cm，封釉层1cm，水泥稳定层基层厚度为30cm）

运行BBR文件

Dh？ 0

L=？ 38660

i？ 0

K？ 0

H= 76.9412

L=？ 38900

i？ 0.02

K？ 4.50

H= 85.2200

退出程序，重新运行

Dh？ 0.370

L=？ 39600

i？ 0.02

K？ 4.50

H= 93.0509