道路中边桩坐标计算EXCEL程序王中伟

利用EXCEL表格计算线路中边桩坐标高速公路施工中，桥梁、隧道施工以及路沿石施工时对测量要求相当严格，要求总体宽度误差在1cm以内。

这就要求测量工程师必须把线路边桩都计算出来，每隔10米（曲线段）或者20米（直线段）放样线路边桩，以便指导施工。

如果采用常规的计算器计算，不仅繁琐、费力而且容易出现差错。

本人在计算中边桩坐标时，试着利用EXCEL表格功能，编辑函数大批量计算线路中边桩坐标，达到了高效、准确的目的。

由于直线部分相对简单，复曲线又特别繁琐。

现就圆曲线举例计算如下：例：在浙江省龙丽一级公路施工中，左线圆曲线起点里程为K86+966.6，曲线半径为1500m，曲线左偏，起始方位角为147°43 ′58.2″，线路中心和隧道中心偏差15cm。

计算隧道中心线和衬砌边线（半径5.1m）。

1.在B4方格中输入：=$B$3+3000*SIN((A4-$A$3)/3000)*COS(2.5784244-(A4-$A$3)/3000)+ 0.15*COS(4.149220679-2*(A4-$A$3)/3000)按回车键确认即可以计算出该里程隧道中心线X轴坐标。

2.在C4方格中输入：=$C$3+3000*SIN((A4-$A$3)/3000)*SIN(2.5784244-(A4-$A$3)/3000)+0.15*SIN(4.149220679-2*(A4-$A$3)/3000)按回车键确认即可以计算出该里程隧道中心线Y轴坐标。

同理在3.D4方格中输入：=$B$3+3000*SIN((A4-$A$3)/3000)*COS(2.5784244-(A4-$A$3)/3000)+ COS(4.149220679-2*(A4-$A$3)/3000) *（0.15+5.1）4.在E4方格中输入：=$C$3+3000*SIN((A4-$A$3)/3000)*SIN(2.5784244-(A4-$A$3)/3000)+SIN(4.149220679-2*(A4-$A$3)/3000)*（0.15+5.1）其中 2.5784244为方位角147°43 ′58.2″的弧度表达方式，4.149220679为起始方位角147°43 ′58.2″加90°后的弧度值。

高速公路中边桩坐标程序单曲线主程序<R>1、LbI 1：{S}：S“KC=”：prog“ys”：L=S-Z：prog“HQ”回车D=0零：{D}：D“LR=0零=＞推出Z”：D＞0零=＞{H}：H“JIAJIAO”：x=x+Dcos（Q+H）：y=y+Dsin（Q+H）◣空心回车2、D＜0零=＞{H}：H：“Jiajiao小写”：D=AbsD：X=X+cos（Q+H）：y=y+Dsin（Q+H）：PoI（X-E“X0”，Y-O“Y0”）：J＜0零=＞J=J+360◣空心回车J“N”=⊿I“L”=⊿X⊿Y⊿Goto1回车S：里程Z：ZH点里程A：转角R：半径G交点X坐标V交点Y坐标C缓和曲线长B：起点方位角坐标计算子程序（HQ）K=ARπ÷180+C：p=C2÷24R-C∧4÷2688R∧3：T=（R+P）tan(A÷2)+ C÷2- C∧3÷240 R∧2:N=-1:L＞K-C=>推出L=K-L：N=1：◣空心回车L≤0零=>X=L:y=0零:Q=0零: Goto1:◣空心回车L≤C=>X=L-L∧5÷40 R∧2C∧2:y=L∧3÷6Rc-L∧7÷336 R∧3C∧3:Q=90L∧2÷πRC:Goto1◣空心回车L≤K-C=>Q=180（L-C）÷πR+90C÷πR：X=RsinQ+C÷2-C∧3÷240 R∧2：y=R(1-cosQ)+P◣空心回车LbI1:N=1=>Q=A-Q：I=X：J=Y：X=T+（T-I）cosA-JsinA：Y=（T-I）sinA+J cosA：◣空心回车I=X：J=Fy：X=G+Tcos（180+B）+IcosB-JsinB：y=V+Tsin（180+B）+IsinB+JcosB：Q=B+FQ回车K 曲线总长T切线长C 缓和曲线长N第①或第②缓和段P 内移值L所求点到ZH距离A偏角X、y参数方程R 半径Q切线方位角I△x J△y B起点方位角曲线要素子程序（YS）S≤24800.423＝>推出A=75度33分23秒：R=125：C=55：Z=24580.585：G=4666.705：V=218.997：B=49度45分16秒：F=1：Goto2⊿回车S≤24980.783推出A=40度12分03秒：R=200：C=40：Z=24800.423：G=4540.417：V=397.236：B=125度19分14秒：F=-1：Goto2⊿回车LbI2回车注：S≤HZ点里程＝＞A“转角”：R=半径：C=缓和曲线长：Z=ZH 点里程：G=交点X坐标：V=交点Y坐标：B=起点方位角；F为转向左-1、右+1使用说明：KC=：桩号LR=0=＞Z:默认为中桩。

附件：道路中桩、边桩计算公式(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 ▲。

小白自学ExcelVBA道路坐标高程计算程序流程1声明1.1 本人是个工作时间不长测量施工员。

1.2听说想轻松愉快的干好测量得学门编程，又听说ExcelVBA是个简单的可视化的常用的……1.3 本人很懒，面对厚厚的全是“昏天暗地”代码的ExcelVBA教科书籍，实在是没心没力。

其实俺还是懂点儿VB编程的（学校学过“=”是赋值，if、while、end的意思，嘿嘿，还有画窗体图框，其他的就都还给老师了）。

1.4 谷歌、度娘，是咱的好帮手。

1.5 仰望的存在，神一般的存在，无私的化身！道路中边桩坐标计算程序120424.xls（王中伟”教授”QQ：595077）曲线坐标计算程序VBA 4.6.xls（陈超”中铁”QQ：295188316）陈师傅的VBA代码乍一看感觉比王老师的让人头痛，怎么说咱们也是个有理想有目标的人，要啃就啃硬骨头！就选~~~嘿嘿~~~王老师的程序吧！1.6 道路中边桩坐标计算程序120424.xls代码获取。

王老师、陈师傅担心不懂编程的人，不小心修改代码导致程序错误，设置了vba密码（个人认为），其次我跟他们不熟，没有正大光明的要的密码，于是求助谷歌度娘，这个关于VBA 密码破解，于是很不厚道的得到了密码，为了方便本次学习，又很不厚道的把道路中边桩坐标计算程序120424.xls的密码贴出来（df750726），神人勿怪呀！2 代码理解ing…我们通过密码打开vba代码窗口，大家可以先把每个表格打开浏览一下代码，是不是感觉很短呀，感觉应该是超级链接的意思，嘿嘿，主要的代码在模块“lx”里，打开下，代码多吧，这就是咱们要肯的骨头。

现在咱们冲呀！2.1 第一段代码理解貌似第一段代码是定义的意思“定义π值”、'定义桩号点”、“'定义坡段”、“'定义基本曲线”①Public Const PI As Double = 3.14159265358979谷歌翻译：Public Const PI As Double公共常量PI 作为双精度浮点型百度百科：双精度浮点数双精度浮点数（Double）用来表示带有小数部分的实数，一般用于科学计算，用8个字节（64位）存储空间，其数值范围为1.7E-308～1.7E+308，双精度浮点数最多有15或16位十进制有效数字，双精度浮点数的指数用“D”或“d”表示。

路线坐标计算软件
ROAD COORDINATE CALCULATING SOFTWARE
Copyright © 2001
版权所有 All Rights Reserved
1．本计算软件基于EXCEL开发，使用VBA编程，为保证顺利使用，请将宏安全性设为“中”；2．本软件适用于基本型曲线的道路的计算，但不适用于卵形曲线、立交匝道的计算；
3．作者开发有适用于卵形曲线、立交匝道计算的计算软件，需要者请与作者联系；
4．本软件供工程一线技术人员免费使用，但版权和改编权属于作者；
5．本软件使用简单，使用者有任何问题，可与作者联系，作者尽可能提供帮助；
6．用户若有支助意向，请与作者联系。

：Excel 公路中线坐标计算现阶段我国公路工程中已普遍使用大地坐标进行线型的控制及测设，在施工中经常要对中线坐标进行复核、加密，才能满足公路工程施工的需要。

本文是结合公路工程的实际需要，利用Exce电子表格制作的用于由直线、圆曲线、缓和曲线组成的一般公路线型中桩、边桩等计算的通用模版，用于减轻计算工作的劳动强度和提高计算结果的准确度。

一、采用公式1 直线段1.1 中桩坐标计算公式1.2 边桩坐标计算公式2 缓和曲线段2.1 中桩坐标计算公式当P点位于顺时针方向时,其方位角为αE→p=αA→B+90°；当P点位于逆时针方向时,其方位角为αE→p=αA→B-90°。

2.2 边桩坐标计算公式3 圆曲线段3.1 中桩坐标计算公式当E点位于顺时针方向时取“+”,当E点位于逆时针方向时取“-”。

3.2 边桩坐标计算公式XP、YP——未知点P的坐标X1、Y1——各线型起点的坐标（第二曲线段为终点）XA、YA、XB、YB——P点边桩A点、B点的坐标（A为左侧、B为右侧）α1→2——直线段起点的方位角αA→B——各线形起点的切线方位角（第二曲线段为终点）L——P点距各线形起点的长度LS——缓和曲线段缓和曲线长R——各曲线段的半径β——P点的切线角（曲线左转时取“-”、曲线右转时取“+”）T1、T2——P点至边桩A、B的距离（A为T1、B为T2）二、计算模板的建立步骤表1－1表1－21．新建一个工作薄，在其中输入如表1所示的内容；2. 选中工作表A列，打开格式菜单，选中“单元格”，在单元格菜单中选中“数字”栏，自定义单元格格式为“K000＋000.000”。

按此方法分别将其他列设置为如表1所示单元格格式；3. 将“4”行做为路线起点数据行，在“5”行中进行公式编辑；4．在“J5”单元格中输入“=IF(C5=4,RADIANS(IF((G5+H5/60+I5/60/60)<180,(G5+H5/60+I5/60/60)+180,(G5+H5/60+I5/60/60)-18 0)),IF(C5=5,RADIANS(IF(B5=0,G5+H5/60+I5/60/60-E5/2/D5*180/PI(),G5+H5/60+I5/60/60+E5/2/D5*180 /PI())),RADIANS(G5+H5/60+I5/60/60)))”；5．在“K5”单元格中输入“＝A5－$A$4”；6．在“L5”单元格中输入“=$L$4+IF(C5=1,K5*COS(J5),IF(C5=2,2*D5*SIN(K5/2/D5)*IF(B5=0,COS(J5-K5/2/D5),COS(J5+K5/2/D5) ),IF(C5=3,(K5-K5^5/40/D5^2/E5^2)*COS(J5)+(K5^3/6/D5/E5-K5^7/336/D5^3/E5^3)*IF(B5=0,COS(J5-90 *PI()/180),COS(J5+90*PI()/180)),IF(C5=4,(K5-K5^5/40/D5^2/E5^2)*COS(J5)+(K5^3/6/D5/E5-K5^7/33 6/D5^3/E5^3)*IF(B5=0,COS(J5-90*PI()/180),COS(J5+90*PI()/180)),2*D5*SIN(K5/2/D5)*IF(B5=0,COS( J5-K5/2/D5),COS(J5+K5/2/D5))))))”；7．在“M5”单元格中输入“=$M$4+IF(C5=1,K5*SIN(J5),IF(C5=2,2*D5*SIN(K5/2/D5)*IF(B5=0,SIN(J5-K5/2/D5),SIN(J5+K5/2/D5) ),IF(C5=3,(K5-K5^5/40/D5^2/E5^2)*SIN(J5)+(K5^3/6/D5/E5-K5^7/336/D5^3/E5^3)*IF(B5=0,SIN(J5-90 *PI()/180),SIN(J5+90*PI()/180)),IF(C5=4,(K5-K5^5/40/D5^2/E5^2)*SIN(J5)+(K5^3/6/D5/E5-K5^7/33 6/D5^3/E5^3)*IF(B5=0,SIN(J5-90*PI()/180),SIN(J5+90*PI()/180)),2*D5*SIN(K5/2/D5)*IF(B5=0,SIN( J5-K5/2/D5),SIN(J5+K5/2/D5))))))”；8．在“N5”单元格中输入“=IF(B5=0,J5-RADIANS(IF(C5=2,K5/D5*180/PI(),IF(C5=3,K5^2/2/D5/E5*180/PI(),IF(C5=4,K5^2/2/D5 /E5*180/PI(),IF(C5=5,K5/D5*180/PI()))))),J5+RADIANS(IF(C5=2,K5/D5*180/PI(),IF(C5=3,K5^2/2/D5 /E5*180/PI(),IF(C5=4,K5^2/2/D5/E5*180/PI(),IF(C5=5,K5/D5*180/PI()))))))”；9．在“P5”单元格中输入“=IF(C5=4,L5+O5*COS(N5+90*PI()/180),L5+O5*COS(N5-90*PI()/180))”；10．在“Q5”单元格中输入“=IF(C5=4,M5+O5*SIN(N5+90*PI()/180),M5+O5*SIN(N5-90*PI()/180))”；11．在“S5”单元格中输入“=IF(C5=4,L5+O5*COS(N5-90*PI()/180),L5+O5*COS(N5+90*PI()/180))”；12．在“T5”单元格中输入“=IF(C5=4,M5+O5*SIN(N5-90*PI()/180),M5+O5*SIN(N5+90*PI()/180))”；13．选定“5行”向下复制到表格最后一行；14．点击“文件”菜单中的“另存为”选项，输入文件名为“坐标计算”在文件类型中选择为“模版”并点击“保存”。

用EXCEL进行中桩坐标、高程计算（铜汤高速第一合同段）测量计算、平差软件种类繁多，但无论什么软件都不能包罗万象，无法满足众多数据处理的需要，即使某种计算软件的功能非常齐全，但此类软件一般都是针对性开发的，而像小件坐标计算，例如承台坐标，则显得非常烦琐。

本文结合本项目介绍了应用EXCEL进行测量计算的方法和思路，并提供了一些常用的计算公式。

1、 EXCEL基本计算方法打开EXCEL后显示的表格，将相关数据依次输入到相应的表格上中。

当原始数据输入完毕后，再在需要显示结果的相应单元格内输入计算公式，EXCEL则会按照公式自动计算出结果，下面以两点坐标计算距离为例，简单说明下：第一步：在相应位置输入表头、序号、坐标后，图示1所示。

第二步在显示距离格内输入计算公式=SQRT((B4-D4)^2+(C4-E4)^2)并按回车，距离的结算结果就显示出来。

2、在测量计算中，最常见的是三角函数的处理。

EXCEL中备用的三角函数很齐全，但计算中是以弧度为单位来处理的，在计算中角度（度、分、秒）需要转换成弧度后，即可进行三角函数的计算了。

3、结合本项目分别介绍了直线、圆曲线、竖曲线的计算事例铜汤高速公路第一合同段，起点里程K78+060.618,终点里程K82+825,全长4764.4米，其中包括朱家畈大桥，水桥湖一桥、二桥三座桥梁。

全线平面首尾为直线段，中间由三段圆曲线连接，（后附全线平面图），计算原始数据取自设计文件第一册《直线、曲线及转角表》(图号为S3-4);全线包括6段竖曲线，其中凹、凸曲线各3段，计算原始数据采用设计文件第一册《路线纵断面图》（图号为S3-3）。

在表格中，桩号栏中的ZJF8、SQA7、SQB6分别代表朱家畈大桥8号墩中心，水桥湖一桥7号墩中心和水桥湖二桥6号墩中心。

坐标和高程的单位均为m，角度单位为弧度。

一、坐标计算方法：a、直线段：坐标增量法例如：第一段直线范围K78+060.618~K78+773.787坐标计算已知：BP点里程为K78+060.618，坐标（3412768.816，484851.623）起始方位角a＝140°45′46.1″＝2.456774421（弧度）计算直线范围内的K78+200坐标如下：在EXCEL表格中输入相应的表头和需要计算的里程后在C36格内输入公式=3412768.816+（B21-78060.618）* COS（2.456774421）在D36格内输入公式=484851.623+（B21-78060.618）* SIN（2.456774421）X=3412768.816+(78200-78060.618)×COS2.456774421＝3412660.860Y=484851.623+(78200-78060.618)×SIN2.456774421＝484939.787 图示：b、曲线段：偏角法第一段圆曲线K78+773.787～K79+524.582上坐标计算：起点（即ZY点）里程为K78+773.787，同时位于第一段直线上，可通过坐标增量法算出ZY点坐标为（3412216.442，485302.725）；始方位角a＝140°45′46.1″＝2.456774421（弧度）；曲线圆心在线路右侧，曲线半径R＝4500m。

如何计算边桩坐标
编写：路桥原创软件QQ群 74313213
一、 需要用到的功能
1、Excel程序
2、【纬地V5.88版本】—【工具】—【逐桩用地表】
二、 实现步骤
1、首先熟悉 TF 文件的格式，第 1 列为桩号，第 23、24 列分别为左侧和右
侧的用地宽
2、打开excel程序，把需要的桩号放到第1列，如下图，可添加所需桩号
3、在excel的 B 列，添加21个0，注意，不能多不能少，一定是21个
4、在 C 列填写左侧距离，不管正负，此距离为法向；D 列填写右侧距离，
不管正负，此距离为法向
5、在E列填写11个0
6、另存为txt 文件或者直接复制excel里面的内容
7、打开生成的txt 文档或者打开记事本把复制内容粘贴进去
8、在文本最前面加一行数据：HINTCAD5.83_TF_SHUJU，注意是独立行
9、把txt 文档后缀改成TF
10、 【纬地 V5.88 版本】—【项目】—【项目管理器】 ，找到最后的 TF
文件，重新指定为111.TF，保存退出
11、 【纬地 V5.88 版本】—【工具】—【逐桩用地表】即可出表，批量
的出
三、 备注
1、如果有问题可以到QQ群：74313213探讨
2、写得太快，上面弄错了，图例中桩号后面都是 22个0，有误，应该是21
个0。