下载文档原格式
曲线道路坐标计算§1 曲线要素计算缓和曲线是在不改变直线段方向和保持圆曲线半径不变的条件下，插入到直线段和圆曲线之间的。
其曲率半径ρ从直线的曲率半径∞（无穷大）逐渐变化到圆曲线的半径R ，在缓和曲线上任意一点的曲率半径ρ与缓和曲线的长度l 成反比，以公式表示为：l1∝ρ 或 C l =⋅ρ（C 为常数，称曲线半径变更率）。
当o l l =时，R =ρ，应有o l R l C ⋅=⋅=ρ以上几式是缓和曲线必要的前提条件。
在实际应用中，可采取符合这一前提条件的曲线作为缓和曲线。
常用的有辐射螺旋线及三次抛物线，我国采用辐射螺旋线。
为了在圆曲线与直线之间加入一段缓和曲线o l ，原来的圆曲线需要在垂直于其切线的方向移动一段距离p ，因而圆心就由'O 移到O ，而原来的半径R 保持不变，如图。
由图中可看出，缓和曲线约有一半的长度是靠近原来的直线部分，而另一半是靠近原来的圆曲线部分，原来圆曲线的两端其圆心角o β相对应的那部分圆弧，现在由缓和曲线所代替，因而圆曲线只剩下缓圆点(HY )到圆缓点(YH )这段长度即y l 。
o β为缓和曲线的切线角，即缓圆点或圆缓点切线与直缓点或缓直点切线的交角，亦即圆曲线HY→YH 两端各延长2ol 部分所对应的圆心角。
γ为缓和曲线总偏角，即从直缓点(ZH )测设缓圆点(HY )或从缓直点(HZ )测设圆缓点(YH )的偏角。
q 为切线增量(切垂距)，即ZH (或HZ )到从圆心O 向ZH (或HZ )的切线作垂线垂足的距离。
p 为圆曲线内移值，即垂线(从圆心O 向ZH (或HZ )的切线作垂线)长与圆曲线半径R 之差。
§1.1 不等长缓和曲线要素计算：在铁路曲线测设中，线路曲线一般是由相等的两条缓和曲线中间加一个圆曲线构成，有时还会出现由两个不等长的缓和曲线中间加一个圆曲线构成的特殊情况，如图：缓和曲线长分别为1o l 、2o l ， 切线长分别为1T 、2T ，曲线偏角(线路转角)为α，圆曲线半径为R ，圆曲线长为y l ，曲线长为L ，外矢距为E ，切曲差为J ，(缓和曲线后)圆曲线内移值分别为1p 、2p ，(缓和曲线)切线增量分别为1q 、2q ，缓和曲线偏角分别为1o β、2o β ， 回旋线参数分别为121o Rl A =、222o Rl A =各曲线要素计算公式如下：231112402R ll q o o -=232222402Rll q o o -=341211268824R lR l p o o -=342222268824RlR l p o o -=ααsin )(2)(12111p p tgp R q T -+++=ααsin )(2)(21222p p tgp R q T -+++=Rl R l o o o ππβ111901802=⋅=Rl R l o o o ππβ222901802=⋅=︒--++=180)(2121Rl l L o o o o πββα从以上公式可以看出，当21o o l l =时，就是等长(对称)缓和曲线的情况。
纯Excel公式开发的施工测量数据计算表版本号：1.4功能：可以根据输入的起点桩号、各交点坐标、各交点的圆曲线半径、缓和曲线长（支持不对称缓和曲线）直接生成标准的《直线曲线及转角表》和积木法平曲线参数表。
可以计算非对称基本型（直-缓1-圆-缓2）以及下列由其变化而成共6种线型的任意交角中边桩坐标。
1、单一完整缓和曲线；分缓1和缓2两种情况，此曲线为前直线、圆曲线和另一缓和曲线长度均为0的特殊曲线。
2、单一圆曲线，即前直线和缓和曲线长度均为零的曲线。
3、卵型曲线，即圆曲线长为0的曲线（前直线有无不限，此类曲线在程序中要加入防止R被作除数时出错）4、缓1+圆曲线型曲线，即缓2长度为零的曲线（前直线长度不限）5、圆曲线+缓2型曲线，即前直线和缓1长度为零的曲线。
6、单一直线，即所有缓和曲线和圆曲线的长度均为零的曲线（常用于路线终点直线计算，此类曲线需防止已知量为零并被当作被除数时出错）可以根据输入的桩号自动生成用于AutoCAD脚本绘图的代码，拷贝后保存为*.scr文件，在任意版本AutoCAD中加载展绘中线和桩号。
可以用近似公式计算竖曲线高程。
对于以变量赋值形式作为数据库的CASIOfx-5800P计算器程序，可自动生成交点法、积木法、竖曲线数据库。
使用说明：此表用纯Excel公式编成，理论上可计算多达30000个交点，但考虑到实际情况，只拉出了200个交点的表格，所有的已知数据在此表中输入，结果在《直曲表》中查看。
此表格经多条100公里以上的线路计算，误差均小于1mm，可用于设计数据复核。
因为很多线路的桩号，缓和曲线参数，主点桩号等数值均只取2位小数，而交点坐标一般会给出小数点后三位，半径和缓和曲线长一般均为整数，用此表计算可以帮助有程序的朋友在建立数据库时提高计算精度。
又有朋友手上只有积木法的程序，对于交点法不想另编，本表又提供了一种从交点法表格直接转为积木法的功能，已知数据还是没有变，根据交点个数，拉了800行。
用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。
路线坐标计算软件
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快速计算任何难度的道路中线、边线坐标主程序Public i As IntegerPublic pi As DoubleSub TP()Dim ii As IntegerDim k(1000) As DoubleDim xzq, yzq, kq, xzh, yzh, kzh, xjd, yjd, kjd, khy, kyh As Double'直线区域pi = 3.14159265358979xzq = 71862.642yzq = 63474.651kq = 0 '因为直线连接终点为ZH点，与圆曲线起点为同一点，所以在直线区域不定义ZH点参数'直线区域'曲线区域xzh = 71858.3267yzh = 63375.2684kzh = 99.4763xhz = 71909.3687yhz = 63283.8076 '曲线区域定义内容有：ZH（坐标、里程）、HZ（坐标、里程）、JD（坐标、里程）khz = 212.3392 'R（半径）、LS（缓和曲线长度）、HY（里程）、YH（里程）xjd = 71855.658yjd = 63313.806kjd = 160.9966ls = 30r = 75khy = 129.4763kyh = 182.3385'曲线区域i = 2 '从第二格开始读取数据所以定义I=2ii = 1 '桩号从第一个开始启用，所以定义II=2Dok(ii) = Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 1) '定义桩号等于读取数据If Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 1) = "" Then Exit Do '当没有数据读取时退出循环If k(ii) < kq Then '若计算点超过计算起点给予提示并退出程序MsgBox ("猪啊！！你的输入的桩号居然超过计算起点桩号")Exit SubElseIf k(ii) <= kzh Then '若计算点在ZH点前，则进入直线程序Call zx(xzq, yzq, kq, xzh, yzh, kzh, k(ii))ElseIf kzh < k(ii) And k(ii) <= khy Then '若计算点在ZH和HY之间则调入前段缓和曲线程序Call qhhqx(xzh, yzh, kzh, xhz, yhz, khz, xjd, yjd, kjd, ls, r, k(ii))ElseIf khy < k(ii) And k(ii) <= kyh Then '若计算点在HY和YH之间则调入圆曲线程序Call yqx(xzh, yzh, kzh, xhz, yhz, khz, xjd, yjd, kjd, ls, r, k(ii))ElseIf kyh < k(ii) And k(ii) <= khz Then '若计算点在YH和HZ之间则调入后段缓和曲线程序Call hhhqx(xzh, yzh, kzh, xhz, yhz, khz, xjd, yjd, kjd, ls, r, k(ii))ElseMsgBox ("笨啊！！数据已超出计算范围了") '若出现超出范围的桩号则给与提示并退出程序Exit SubEnd Ifi = i + 1ii = ii + 1LoopEnd Sub直线模块Sub zx(ByV al xzq As Double, ByV al yzq As Double, ByV al kq As Double, ByV al xzh As Double, ByV al yzh As Double, ByV al kzh As Double, ParamArray k())fw = fwj(xzh, xzq, yzh, yzq) '首先调入方位角程序计算直线方位角x = xzq + (k(ii) - kq) * Cos(fw) '然后根据桩号和长度计算出坐标值y = yzq + (k(ii) - kq) * Sin(fw)zdfm = dfm(fw) '将弧度形式的前进方位角转换度分秒形式'输出坐标值以弧度和度分秒形式的前进方位角Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 2) = xWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 3) = yWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 4) = fwWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 5) = zdfmEnd Sub圆曲线模块Sub yqx(ByV al xzh As Double, ByV al yzh As Double, ByV al kzh As Double, ByV al xhz As Double, ByV al yhz As Double, ByV al khz As Double, ByV al xjd As Double, ByV al yjd As Double, ByV al kjd As Double, ByV al ls As Double, ByV al r As Double, ParamArray k())l = Abs(k(ii) - kzh) '计算ZH点（因为以直缓点起算）到待求桩号的弧度长度ly = l - ls / 2 '计算圆弧长度p = ls ^ 2 / 24 / r - ls ^ 4 / 2688 / r ^ 3 '曲线内移值m = ls / 2 - ls ^ 3 / 240 / r ^ 2 '曲线切线长增量u = r * Sin(ly / r) + m '偏量坐标计算v = r * (1 - Cos(ly / r)) + p'调入方位角fwq = fwj(xjd, xzh, yjd, yzh) '计算ZH点方位角fwh = fwj(xhz, xjd, yhz, yjd) '计算HZ点方位角（此角作用是用来推算曲线是左偏还是右偏）'调入偏角判定nq = n(fwq, fwh) '计算偏角方向，左偏为-1右偏为1'计算坐标x = u * Cos(fwq) - nq * v * Sin(fwq) + xzhy = u * Sin(fwq) + nq * v * Cos(fwq) + yzhd = (90 * (2 * l - ls) / pi / r) * pi / 180 '计算圆曲线上的偏角（此句要点为角度必须转换为弧度即：pi/180）fw = fwq + d * nq '计算前进方位角zdfm = dfm(fw) '将弧度形式的前进方位角转换度分秒形式'输出坐标值以弧度和度分秒形式的前进方位角Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 2) = xWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 3) = yWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 4) = fwWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 5) = zdfmEnd Sub前缓和段Sub qhhqx(ByV al xzh As Double, ByV al yzh As Double, ByV al kzh As Double, ByV al xhz As Double, ByV al yhz As Double, ByV al khz As Double, ByV al xjd As Double, ByV al yjd As Double, ByV al kjd As Double, ByV al ls As Double, ByV al r As Double, ParamArray k())l = Abs(k(ii) - kzh) '计算测点到ZH点的距离u = l - l ^ 5 / 40 / r ^ 2 / ls ^ 2 + l ^ 9 / r ^ 4 / ls ^ 4 / 3456 '计算偏量v = l ^ 3 / 6 / r / ls - l ^ 7 / 336 / r ^ 3 / ls ^ 3'用公式二则以下两计算可省略(圆曲线和后缓和曲线也同理)Rem t = Atn(v / u)Rem s = Sqr(u ^ 2 + v ^ 2)'调入方位角计算fwq = fwj(xjd, xzh, yjd, yzh) '计算ZH点方位角fwh = fwj(xhz, xjd, yhz, yjd) '计算HZ点方位角（此角作用是用来推算曲线是左偏还是右偏）'调入偏角判定nq = n(fwq, fwh) '计算偏角方向，左偏为-1右偏为1'结果计算Rem x = xzh + s * Cos(fwq + nq * t)Rem y = yzh + s * Sin(fwq + nq * t)x = u * Cos(fwq) - nq * v * Sin(fwq) + xzh '经过测试，计算结果中的两种公式计算结果是一样的y = u * Sin(fwq) + nq * v * Cos(fwq) + yzhd = (90 * l * l / pi / r / ls) * pi / 180 '计算缓和曲线上的偏角（此句要点为角度必须转换为弧度即：pi/180）fw = fwq + d * nq '计算前进方位角zdfm = dfm(fw) '将弧度形式的前进方位角转换度分秒形式'输出坐标值以弧度和度分秒形式的前进方位角Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 2) = xWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 3) = yWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 4) = fwWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 5) = zdfmEnd Sub后缓和段模块Sub hhhqx(ByV al xzh As Double, ByV al yzh As Double, ByV al kzh As Double, ByV al xhz As Double, ByV al yhz As Double, ByV al khz As Double, ByV al xjd As Double, ByV al yjd As Double, ByV al kjd As Double, ByV al ls As Double, ByV al r As Double, ParamArray k())l = Abs(k(ii) - khz) '计算测点到HZ点的距离（后缓和曲线是以HZ点为起点）u = l - l ^ 5 / 40 / r ^ 2 / ls ^ 2 + l ^ 9 / r ^ 4 / ls ^ 4 / 3456 '计算偏量v = l ^ 3 / 6 / r / ls - l ^ 7 / 336 / r ^ 3 / ls ^ 3Rem t = Atn(v / u)Rem s = Sqr(u ^ 2 + v ^ 2)'调入方位角计算fwq = fwj(xjd, xzh, yjd, yzh) '计算ZH点方位角fwh = fwj(xhz, xjd, yhz, yjd) '计算HZ点方位角（此角作用是用来推算曲线是左偏还是右偏）'调入偏角判定nh = n(fwh, fwq) '计算偏角方向，左偏为-1右偏为1（注意：因为是从后HZ点起算，所以必须将HZ点方位角放在前ZH放在后）'结果计算Rem x = xzh + s * Cos(fwq + nq * t)Rem y = yzh + s * Sin(fwq + nq * t)x = xhz - (u * Cos(fwh) - nh * v * Sin(fwh)) '经过测试，计算结果中的两种公式计算结果是一样的y = yhz - (u * Sin(fwh) + nh * v * Cos(fwh))d = (90 * l * l / pi / r / ls) * pi / 180 '计算缓和曲线上的偏角（此句要点为角度必须转换为弧度即：pi/180）fw = fwh + d * nh '计算前进方位角zdfm = dfm(fw) '将弧度形式的前进方位角转换度分秒形式'输出坐标值以弧度和度分秒形式的前进方位角Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 2) = xWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 3) = yWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 4) = fwWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 5) = zdfmEnd Sub偏角模块Function n(ByV al fw1 As Double, ByV al fw2 As Double) As Doublepj = fw1 + pi - fw2 '前进的右角pjIf pj - pi > 0 Then '当右角pj-pi 〉0时为左偏否则为右偏n = -1Else: n = 1End IfEnd Function方位角模块Function fwj(ByV al x1 As Double, ByV al x2 As Double, ByV al y1 As Double, ByV al y2 As Double) As Double'计算增量x0 = x1 - x2y0 = y1 - y2'由增量判断方位角所在象限，不同象限取不同的值If x0 = 0 And y0 > 0 Thenfwj = pi / 2 '当在大地坐标中偏量在X轴上的值时ElseIf x0 = 0 And y0 < 0 Thenfwj = 3 * pi / 2 '当在大地坐标中偏量在负X轴上的值时ElseIf x0 < 0 Thenfwj = Atn(y0 / x0) + pi '当在大地坐标中偏量在第二第三象限上的值时ElseIf x0 > 0 And y0 < 0 Thenfwj = Atn(y0 / x0) + 2 * pi '当在大地坐标中偏量在第四象限上的值时Elsefwj = Atn(y0 / x0) '当在大地坐标中偏量在第一象限上的值时End IfEnd Function度分秒模块Function dfm(ByV al ao As Double) As V ariantao = ao * 180 / pi '将弧度转化为度jd = Int(ao) '计算度jf = Int(ao * 60 - jd * 60) '计算分jmx = (ao - jd - jf / 60) * 3600 '计算秒jm = Left(jmx, 8) '因为拆分出来的秒数经常占到十多位，所以只取秒数的前八位dfm = jd & "°" & jf & "′" & jm & "″" '连接度分秒End Function边桩公式(此公式在电子表格中直接输入)：=B2+J2*COS(D2+RADIANS(L2)+PI()) =C2+J2*SIN(D2+RADIANS(L2)+PI()) =B2+K2*COS(D2+RADIANS(M2)) =C2+K2*SIN(D2+RADIANS(M2))=B3+J3*COS(D3+RADIANS(L3)+PI()) =C3+J3*SIN(D3+RADIANS(L3)+PI()) =B3+K3*COS(D3+RADIANS(M3)) =C3+K3*SIN(D3+RADIANS(M3))=B4+J4*COS(D4+RADIANS(L4)+PI()) =C4+J4*SIN(D4+RADIANS(L4)+PI()) =B4+K4*COS(D4+RADIANS(M4)) =C4+K4*SIN(D4+RADIANS(M4))=B5+J5*COS(D5+RADIANS(L5)+PI()) =C5+J5*SIN(D5+RADIANS(L5)+PI()) =B5+K5*COS(D5+RADIANS(M5)) =C5+K5*SIN(D5+RADIANS(M5))=B6+J6*COS(D6+RADIANS(L6)+PI()) =C6+J6*SIN(D6+RADIANS(L6)+PI()) =B6+K6*COS(D6+RADIANS(M6)) =C6+K6*SIN(D6+RADIANS(M6))=B7+J7*COS(D7+RADIANS(L7)+PI()) =C7+J7*SIN(D7+RADIANS(L7)+PI()) =B7+K7*COS(D7+RADIANS(M7)) =C7+K7*SIN(D7+RADIANS(M7))=B8+J8*COS(D8+RADIANS(L8)+PI()) =C8+J8*SIN(D8+RADIANS(L8)+PI()) =B8+K8*COS(D8+RADIANS(M8)) =C8+K8*SIN(D8+RADIANS(M8))=B9+J9*COS(D9+RADIANS(L9)+PI()) =C9+J9*SIN(D9+RADIANS(L9)+PI()) =B9+K9*COS(D9+RADIANS(M9)) =C9+K9*SIN(D9+RADIANS(M9))=B10+J10*COS(D10+RADIANS(L10)+PI()) =C10+J10*SIN(D10+RADIANS(L10)+PI()) =B10+K10*COS(D10+RADIANS(M10)) =C10+K10*SIN(D10+RADIANS(M10))=B11+J11*COS(D11+RADIANS(L11)+PI()) =C11+J11*SIN(D11+RADIANS(L11)+PI()) =B11+K11*COS(D11+RADIANS(M11)) =C11+K11*SIN(D11+RADIANS(M11))=B12+J12*COS(D12+RADIANS(L12)+PI()) =C12+J12*SIN(D12+RADIANS(L12)+PI()) =B12+K12*COS(D12+RADIANS(M12)) =C12+K12*SIN(D12+RADIANS(M12))=B13+J13*COS(D13+RADIANS(L13)+PI()) =C13+J13*SIN(D13+RADIANS(L13)+PI()) =B13+K13*COS(D13+RADIANS(M13)) =C13+K13*SIN(D13+RADIANS(M13))=B14+J14*COS(D14+RADIANS(L14)+PI()) =C14+J14*SIN(D14+RADIANS(L14)+PI()) =B14+K14*COS(D14+RADIANS(M14)) =C14+K14*SIN(D14+RADIANS(M14))=B15+J15*COS(D15+RADIANS(L15)+PI()) =C15+J15*SIN(D15+RADIANS(L15)+PI()) =B15+K15*COS(D15+RADIANS(M15)) =C15+K15*SIN(D15+RADIANS(M15))=B16+J16*COS(D16+RADIANS(L16)+PI()) =C16+J16*SIN(D16+RADIANS(L16)+PI()) =B16+K16*COS(D16+RADIANS(M16)) =C16+K16*SIN(D16+RADIANS(M16))=B17+J17*COS(D17+RADIANS(L17)+PI()) =C17+J17*SIN(D17+RADIANS(L17)+PI()) =B17+K17*COS(D17+RADIANS(M17)) =C17+K17*SIN(D17+RADIANS(M17))=B18+J18*COS(D18+RADIANS(L18)+PI()) =C18+J18*SIN(D18+RADIANS(L18)+PI()) =B18+K18*COS(D18+RADIANS(M18)) =C18+K18*SIN(D18+RADIANS(M18))=B19+J19*COS(D19+RADIANS(L19)+PI()) =C19+J19*SIN(D19+RADIANS(L19)+PI()) =B19+K19*COS(D19+RADIANS(M19)) =C19+K19*SIN(D19+RADIANS(M19))=B20+J20*COS(D20+RADIANS(L20)+PI()) =C20+J20*SIN(D20+RADIANS(L20)+PI()) =B20+K20*COS(D20+RADIANS(M20)) =C20+K20*SIN(D20+RADIANS(M20))=B21+J21*COS(D21+RADIANS(L21)+PI()) =C21+J21*SIN(D21+RADIANS(L21)+PI()) =B21+K21*COS(D21+RADIANS(M21)) =C21+K21*SIN(D21+RADIANS(M21))=B22+J22*COS(D22+RADIANS(L22)+PI()) =C22+J22*SIN(D22+RADIANS(L22)+PI()) =B22+K22*COS(D22+RADIANS(M22)) =C22+K22*SIN(D22+RADIANS(M22))=B23+J23*COS(D23+RADIANS(L23)+PI()) =C23+J23*SIN(D23+RADIANS(L23)+PI()) =B23+K23*COS(D23+RADIANS(M23)) =C23+K23*SIN(D23+RADIANS(M23))呵呵，这个就是传说已久的EXCEL自动计算曲线坐标．给大家研究研究．欢迎在此基础上开发增加新模块。
EXCEL辅助计算公路中桩与边桩参数的应用- 水利施工简介：在现今的道路桥梁建设中，全站仪以其高精度、高效率的特点，在测量放样中得到广泛的应用，它将传统的二维测量提升到了三维测量的高度。
但大量的坐标计算加重了测量人员的工作强度，也影响了全站仪的使用效率。
本人根据《公路工程施工测量》提供的FX-4500程序，利用EXCEL函数与逻辑循环语来计算公路中桩及边桩参数，大大提高了计算的可靠性和速度。
关键字：EXCEL 计算公路中桩公路边桩一、前言在当前的高速公路施工中，全站仪以其高精度、使用简便的特点在施工测量放样中得到广泛的应用。
它的使用还有一个意义，它将测量从传统的二维测量提升到三维测量的高度，使测量效率大为提高。
全站仪使用时需计算大量的坐标，坐标的正确和精度直接关系到测量结果的准确度；另外，对于现场计算的坐标，能否迅速的得出结果也会影响全站仪的使用效率，所以准确、快速的坐标计算也是全站仪使用中的关键一环。
坐标的计算涉及到大量的数据和公式，计算时极易出错。
测量中使用最多、计算最为复杂的是公路中桩、边桩坐标和方位角，其他的施工测量都是以此作为参照，经过简单运算即可得结果。
现测量员普遍采用可编程计算器帮助计算，使用前需根据不同的计算器特性输入程序，计算时根据不同的路段选择不同的程序输入参数，逐个得到坐标或方位角，还得反复校核。
这对于大批量计算来说是一个不小的工程。
现在，电脑在施工单位已得到使用普及，在工程量及数测量数据处理过程中，Office中EXCEL 经常是测量技术人员利用的工具。
本人根据本人从2002年至2005年的利用电脑处理测量资料的经验结合EXCEL的强大功能，总结出了用电脑计算公路中桩、左侧、右侧坐标的可行方法，并利用EXCEL 的逻辑语句编制了一个小的应用程序。
希望能把测量人员从重复性的劳动中解脱出来。
二、工作原理1.应用程序的目标上文已讲到施工测量中使用最多、计算最为复杂的是公路中线逐点坐标和切线方位角，既中桩的X值、Y值。
利用EXCE L 进行道路中边线逐桩坐标计算张宝利1,薛志勤2(1.西北农林科技大学资源与环境学院,陕西杨凌 712100;2.华南农业大学信息与软件学院,广州 510642)摘要:利用坐标增量的基本思想,将道路曲线坐标的计算归结为3种最基本的形式———直线段、圆曲线段、缓和曲线段,直接计算出中、边线的测量坐标,省去了切线支距法坐标向测量坐标的转换过程,简化了计算思路,节约了计算量。
并利用办公软件Excel 的强大计算功能进行实例计算。
关键词:Excel;道路边线;道路中线;逐桩坐标;圆曲线;缓和曲线中图分类号:U452.1文献标志码:B文章编号:1672-741X (2008)03-0302-03Calcul ati on of Coordi n ates of St akes on Centerli n e and Si deli n e Usi n g ExcelZHANG Baoli 1,XUE Zhiqin2(1.College of R esou rces and Environm en t,N orthw est A &F U niversity,Yangling 712100,Shanxi,Ch ina;2.College of Infor m ation and S oft w are,S outh China A gricultural U n iversity,Guangzhou 510642,China )Abstract:The coordinate calculati on of the high way curves is divided int o three basic types,i .e .,straight secti on,cir 2cular curve and transiti on curve,on basis of the coordinate incre ment p rinci p le .W ithout coordinate transf or mati on,the survey coordinates of the stakes on the alignment of the centerline and sideline are directly worked out,which si m p lifies the calculati on p rinci p le and reduces the calculati on quantity .Further more,coordinate calculati on is made f or p r oject cases by using the str ong calculati on functi ons of Excel .Key words:Excel;high way sideline;highway centerline;coordinates of stakes;circular curve;transiti on curve不论是在铁路还是公路的设计和建设中都不可避免地遇到大量的曲线计算。
用Excel程序快速计算公路工程逐桩坐标H现代公路IGHWAY用Excel程序快速计算公路工程逐桩坐标文/左卿引言在公路工程设计及施工过程中，逐桩坐标的计算是一项必不可少的工作，尤其近些年随着全站仪的普及使用，坐标放样法因速度快、精度高而被广泛应用，逐桩坐标计算的工程量也因此而大大增加。
例如：某高速公路某一合同段有20公里的路基，按照常规的设计及施工方法，不考虑曲线加桩，按每20m一个断面，每个断面左右幅2个点，那么需要计算1000个点2000个坐标。
如此大的计算量假如采用传统的计算器计算，不仅费时，还难免出错。
笔者从事多年施工工作，针对这一情况，利用Microsoft Office中的Excel电子表格的编程计算功能，设计了一套包括直线、缓和曲线和圆曲线的坐标计算程序，该计算程序具有快速、精确、程序简单、使用方便等优点，较其他编程软件编制的计算程序更灵活、更适用。
坐标计算直线段逐桩坐标首先介绍函数Value（MIDB（“ZH”，2，n）MIDB(“ZH”，n+3，7)），n表示公里数的位数，当里程公里数为1位数，那么n=1，当公里数为2位数，n=2，当里程公里数为3位数时，n=3……。
如ZH=K136+700，运行该函数，n=3，结果=136700，其目的为将K136+700文本数据（字符串）去掉“K”和“+”，并转换成数字格式。
那么里程≤ZH的前直线上任意桩号坐标计算：在Excel某一单元格中，坐标计算编辑公式为：X= Xj+（T+Value（MIDB （“ZH”，2，n）MIDB(“ZH”，n+3，7)）－Value（MIDB（“K”，2，n）MIDB(“K”，n+3，7)））*cos((A1+180)*PI()/180)Y= Yj+（T+Value（MIDB（“ZH”，2，n）MIDB(“ZH”， n+3，7)）－Value（MIDB（“K”，2，程序设计程序设计初步打开Microsoft Office Excel2022，新建工作表book1，并将其命名为“公路工程逐桩坐标计算表”。
